單元的剛度矩陣由下式積分得到:(四節(jié)點(diǎn)矩形單元應(yīng)該是8×8)該式中的omiga表示單元的空間域,B是形函數(shù)對(duì)空間坐標(biāo)的偏導(dǎo),D是本構(gòu)矩陣,這些矩陣中都不含節(jié)點(diǎn)位移矢量,各種矩陣相乘后得到的8×8矩陣中每一個(gè)元素都是一個(gè)三元函數(shù)。然而我們?cè)诔绦蛑袥](méi)法對(duì)BT*D*B矩陣每一個(gè)元素進(jìn)行解析積分,只能依靠數(shù)值積分手段。在ABAQUS這個(gè)軟件中,所采取的是高斯積分公式。