§3-3 虚功原理推导梁单元的（单元）刚度矩阵
设在力P的作用下，梁单元i-j的两端点分别发生了线位移和角位移，用{?e}来表示梁单元的端点位移（又称结点位移）： {?}?vie??ivj?j?
T
使梁单元发生结点位移{?e}的单元结点力（杆端力）为： {F}?Fie?MiFjMj?
T根据材料力学，如果已知梁的两端点位移，则可求出等截面梁上任意一点的位移（挠度）。即梁上任意一点的位移v(x)可以用{?e}表示出来，设二者的关系为：
?vi?????i?v(x)??N1(x)N2(x)N3(x)N4(x)????{N}T{?e}
?vj????j??e两端点的又设由于某种其他原因，该梁发生了变形，引起梁单元○
位移为（用向量形式表示）：
{?}?vi?ie*?vj?j
?*梁中任意一点的位移为：
?vi?????i?v*(x)??N1(x)N2(x)N3(x)N4(x)????{N}T{?e}*
?vj????j??相对于力P引起的位移v(x),称v*(x)为虚位移
*e的外力虚功和内力虚功 计算梁单元○对梁单元来说，两端点的力即是外力，则外力虚功为：
Wex?{Fe}T?{?e}*?({?e}*)T{Fe}
内力虚功 = 虚应变能
2*22*??dvdvdvdv*W?M?d??EI?d?EI?dx ??l in?l22?ldx2?dxdx?dx??d2N1dvTe ∵ 2?{N''}{?}??2dx?dx2dN2dx2d2N2dx22dN3dx2d2N3dx22dN4?{?e}?[B]{?e} 2?dx?d2N4?e*e*{?}?[B]{?} ?2dx?T2T?d2N1d2v*Te* 2?{N''}{?}??2dx?dxllWin??EI[B]{?e}*?[B]{?e}dx??{?e}*T[B]TD[B]{?e}dx∴
?{?e}*T??[B]TD[B]dx?{?e}
l?{?e}*T?[ke]?{?e}eT式中： [k]?l[B]D[B]dx
?虚功原理：系统保持平衡状态的充要条件是外力虚功=内力虚功 即： Wex?Win {?e*T}{F}?{?}[k]{?}
ee*Tee而虚位移为任意、不为零，所以上式等价于：
eee{F}?[k]{?}
§3-4 单元位移函数的基本概念
对于梁和二力杆，已知单元两端点位移（两端点的力），即可求得
单元内任意一点的位移。对于其他类型的单元即使知道单元结点位移（或单元结点力）也难以求得单元内任意一点的位移。
有限元法的解决方法：假设一个单元结点位移与单元内任意一点位移的关系—多项式形式，称这个关系为单元位移函数，又称单元位移场、单元位移模式。多项式的获得采用插值的方法—对单元结点位移插值。
例：用插值方法获得梁单元的单元位移函数v(x),已知梁两端点i,j的位移分别为：vi?ivj?j
解：设v(x)为一个多项式，其阶数根据已知条件的个数确定。
dv?? ∴ 对于v(x)有4个已知条件：
∵
dxdvdv????j i ③ v(l)= vj ④ ① v(0)=vi ②
dxx?0dxx?l∴ 可设：v(x)?a1?a2x?a3x?a4x 其中：a123a2a3a4为待定系数，现将其用向量{a}表示：
a2a3a4?
?a1??a???x3??2? ?a3???a4??T {a}??a1Tv(x)?{x}{a}??1x则v(x)又可以写成：
x2将4个已知条件代入，可得到以4个待定系数为未知量的方程组： ① v(0)=vi a1?vi
dv??i a2??i ②
dxx?023a?al?al?al?vj ③ v(l)= vj 1234dv④
dxx?l??j a2?2a3l?3a4l2??j
将上面的方程组写成矩阵形式：
?1?0? ?1??00??a1??vi??a????100??i?2???????? 23?lll?a3??vj??12l3l2???a4?????j??00012?l1l2003l2?l→ [A]{a}?{?}e
解此方程组，即可求得4个待定系数a1a2a3a4：
?1?0?3{a}?[A]?1{?e}????l2
??2??l30??vi?0????1???i????v?l??j? 1????j??l2??TT?1eTev(x)?{x}{a}?{x}[A]{?}?{N}{?} ∴
因此我们用插值的方法获得了单元内任意一点位移与单元结点（端点）位移之间的关系---单元位移函数。
有了单元位移函数就可以用能量法推得的单元刚度矩阵的公式来计算单元刚度矩阵了。对于直梁，其虚应变能：
d2v?dv*?dv2d2v*d????lEI2?2dx Win??lM?d???lEI2?dxdxdx?dx?*d2vTee?{N''}{?}?[B]{?}2dx d2v*
Te*e*?{N''}{?}?[B]{?}2dxeT[k]?[B] ?D[B]dx
l