3*3NDxyz=JacobiINV*NDL;%利用雅可比行列式的逆計算形函數(shù)對結(jié)構(gòu)坐標(biāo)的導(dǎo)數(shù)[DN1Dx DN2Dx DN3Dx;DN1Dy DN2Dy DN3Dy;……]end這樣求出形函數(shù)對物理坐標(biāo)的導(dǎo)數(shù)后就可以代入公式(2)幾何方程求出應(yīng)變場矩陣B,經(jīng)過能量原理的推導(dǎo)可以得到單元剛度矩陣的表達(dá)式,注意剛度矩陣的表達(dá)式是一個積分的運算,由于被積函數(shù)較為復(fù)雜,如果代數(shù)積分進(jìn)行計算要消耗大量的計算資源