3c 虚功原理推导单元刚度矩阵

1、3-3虚功原理推导梁单元的（单元）刚度矩阵设在力P的作用下,梁单元i-j的两端点分别发生了线位移和角位移,用5e来表示梁单元的端点位移（又称结点位移）:5e=eveiijj使梁单元发生结点位移5e的单元结点力（杆端力）为：Fe=fmfmbiijj根据材料力学,如果已知梁的两端点位移,则可求出等截面梁上任意一点的位移（挠度）。即梁上任意一点的位移V（x）可以用6e表示出vieiVejij丿来,设二者的关系为：=NT5ev(x)=N(x)N(x)N(x)N(x)1234又设由于某种其他原因，该梁发生了变形，引起梁单元两端点的位移为（用向量形式表示）：5e*=veveiijj梁中任意一点的位移为：v

2、iv*(x)=n(x)N(x)N(x)N(x)=NT5e*0iV00j相对于力P引起的位移V(x),称v*(x)为虚位移计算梁单元9的外力虚功和内力虚功对梁单元来说，两端点的力即是外力，则外力虚功为：W=FeTq5e*=(5e*)TFeex内力虚功=虚应变能WinfMQdol=fEIld2vddx2厂dv*、dx=fEIldv2dmdx2dx2dxd2vdx2=NT5ed2N1dx2d2N2dx2牛牛5e=B5edx2dx2Jd2v*=NT5e*=dx2W=fEIB5e*口B5edxinl=f5e*TBTDB5edx=5e*TqBtDBdxE(5el=5e*TkeE5ed2Nd2N12dx2d

3、x2d2N3dx2d2N45e*=B5e*dx2J式中：ke=fBTDBdxl虚功原理：系统保持平衡状态的充要条件是外力虚功=内力虚功即：W=Wexin5e*TFe=5e*Tke5e而虚位移为任意、不为零，所以上式等价于：Fe=kee3-4单元位移函数的基本概念对于梁和二力杆，已知单元两端点位移(两端点的力)，即可求得单元内任意一点的位移。对于其他类型的单元即使知道单元结点位移(或单元结点力)也难以求得单元内任意一点的位移。有限元法的解决方法：假设一个单元结点位移与单元内任意一点位移的关系多项式形式，称这个关系为单元位移函数，又称单元位移场、单元位移模式。多项式的获得采用插值的方法对单元结点位

4、移插值。例：用插值方法获得梁单元的单元位移函数v(x),已知梁两端点i,j的位移分别为：V-ivjjiijj解：设v(x)为一个多项式，其阶数根据已知条件的个数确定。dv=dx对于V(x)有4个已知条件:v(0)=vidvdxx=0V(l)=vjdvdxx=l=j可设：v(x)=a+ax+ax2+ax3其中：aaaa为待定系数，现将其用向量a表示:a=a1aa23aT4a1则V(x)又可以写成：v(x)=xTa=1xx21ax3/彳2a3a4a=v1iv(0)=v将4个已知条件代入，可得到以4个待定系数为未知量的方程组：101001l100l22l00l33l2a1avi0av3ja04jAa=e解此方程组，即可求得4个待定系数aa=A-15e=10312213012I112003121100-1112vi0iv00jaa34dvdx=0ia=02ix=0v(l)=va+ai+ai2+ai3=j1234dvdx=0a+2ai+3ai2=0234jx=l将上面的方程组写成矩阵形式：vjiv(x)=xTa=xTA-15e=Nt5e因此我们用插值的方法获得了单元内任意一点位移与单元结点（端点）位移之间的关系-单元位移函数。有了单元位移函数就可以用能量法推得的单元刚度矩阵的公式来计算单元刚度矩阵了。对于直梁，其虚应变能：W=fMttfe*=f