第三章
1如图所示一三角形钢板，两个结点固定，对第三个结点施以单位水平位移，测出所施加的力，从而得出相应的刚度系数。其他点依此类推，这样测得的刚度系数所组成的刚度矩阵，是否与按照常规三角形单元刚度矩阵计算公式所得结果一样？用这样实测所得的刚度矩阵能否进行有限元分析？为什么？
解：不一样。单元刚度矩阵中每个元素的物理意义：kij表示单元第j个自由度产生单位位移，其它自由度固定时，第i个自由度产生的节点力。单元刚度矩阵是在单元处于平衡状态的前提下得出的，单元作为分离体看待，作用在它上面的外力（单元力）必是平衡力系，然而研究单元平衡时没有引入约束承受平衡力系作用的无约束单元，其变形是确定，但位移是不能确定的，即单元可发生任意的刚体位移。
不能。因为与有限元中单元与单元之间的约束情况不一样，不能进行有限元分析。 2以位移为基本未知量的有限元法其解具有下限性质，试证明之。 解：系统总位能的离散形式?p?1TT?a?K?a???a??P? 2将求解的方程?K??a???P?带入可得
?p?11TTTaKa?aKa???????????????a??K??a???U 22在平衡情况下，系统总位能等于负的应变能。在有限元解中，由于假定的近似位移模式一般来说总与精确解有差别
的。
设近似解为?p、U、[K]、?a?、??K???a???P?，真实解为?p、U、[K]、?a?、?K??a???P? 且根据最小势能原理，得到的系统的总位能总会比真正的总位能要大，故?p??p则U?U
?a?T??K???a???a??K??a???a?TT?P???a??P?
T则近似解的位移总体上小于精确解的位移
解释如下：单元原是连续体的一部分，具有无限多个自由度，在假定了单元的位移函数后，自由度限制为只有以结点位移表示的有限自由度，引入了更多的约束和限制，使得单元刚度较实际连续体加强了，连续体的整体刚度随之增加，所以有限元解整体上较真实解偏小。
3 请分别阐述单元刚度矩阵和整体刚度矩阵中任一元素的物理意义。
e解：在单刚?K?中，kij表示单元第j个位移产生一单位位移，其它位移为零时，第i个位移方向上引起的节点力。
e在整体刚度中，Kij表示第j个自由度产生一单位位移，其它自由度为零时，第i个自由度上引起的节点力。
1
4 简述虚功原理，且使用虚功原理导出外荷载与节点荷载的等效关系式。
解:虚功原理：变形体中任意满足平衡的力系在任意满足协调条件的变形状态上作的虚功等于零，即体系外力的虚
功与内力的虚功之和等于零。
设?q?为外荷载（此处为体力），?p?为节点荷载，?w?为单元内位移场，???为结点位移场 根据虚功原理
eeee??????p?????w???q?dV
eTeeTeVe由于?w??[N]???故
eVeT???w???q?eTedV?????Ve??eT[N]?q?dV????Tee??eTVT?[N]?q?dV
e则
??????p???????eeTVTT[N]qdV?p?[N]???????q?dV
Ve5 试述弹性力学中按位移求解与有限单元法中按位移求解之间的异同点。 解： 物理模型 基本方程 弹性力学 连续体 几何方程 物理方程 平衡微分方程 解法 解答形式 解答精度 解微分方程 用函数表示 精确解 有限单元法 离散化结构 几何方程 物理方程 结点平衡方程 解代数方程 用数值表示 近似解 6 如果三节点三角形单元绕其中某一个节点作小的刚体转动，其转角为?，证明单元内所有的应力均为零。 解：在三角形单元中?????D??B????
?bi1?B?0???2A?ci?0cibibj0cj0cjbjbm0cm?yj?ym0??1?cm??0?2A??xj?xmbm???0?xj?xmyj?ymym?yi0?xm?xi0?xm?xiym?yiyi?yj0?xi?xj???xi?xj? yi?yj??0?u?u0??y由于三角形单元绕其中某一个节点作小的刚体转动，各节点的位移可表示为：?
v?v??x0?则可知节点位移向量????0,0,??yj,?xj,??ym,?xmT??
T 2
?yi?ym1?0故应变?????B?????2A???xj?xm?0?xj?xmyi?ymym?yi0?xm?xi0?xm?xiym?yiyi?yj0?xi?xj?0??0??0???0?????y1????j??xi?xj?????0? ?x2Aj???0?yi?yj???????y?m???x??m??由于弹性矩阵?D?为常量矩阵，应变向量???为零向量，故?????D????为零向量，即单元内所有的应力为零。 7 二维单元在x,y坐标内平面平移到不同位置，单元刚度矩阵相同吗？在平面内旋转时又怎样？试证明之。
解：二维单元在x,y坐标内平面移到不同位置时，刚度矩阵相同。在平面内旋转时，刚度矩阵也相同。
?1????brcs?kEhA?brbs?1??crcs刚度矩阵T22c?rbs?rs???Br??D??Bs?hA?4(1??2)?????c1??rbs?2brcsc1???
rcs?brbs?2??单元平移或旋转时，bi,ci不变，故单元刚度矩阵不变。
8 判断有限元网格离散合理性
a) 对图1(a)所示的有限元网格，评论网格的优劣性，指出模型中的错误，并加以改正。 b) 评论图1(b)的网格划分合理吗?为什么?请加以改正。
图1 解：（a）网格划分不合理。 1）无过渡单元 2）无边界条件
3）夹角区应力集中，应适当加密风格 4）对称结构网格应对称划分 （b）不合理。
1）左部网格应适当加密
2）由于三角形单元会造成局部精度不够，过渡区可采用其它单元划分 3）右部单元的长宽比较大，就进行适当调整。
9 如图2所示，平面三角形构件以x-y坐标系表示的刚度矩阵方程如下：
?10?2.51.832.5??u?x1??Px1?104?1.832.55.0?2.5?????v??y1???P?y1???2.54.52.5?2.5? ?2.5?2.5?2.52.5??ux2??Px2?????v???y2????Py2??
3
试建立以ux1,uy1,ux2(与图中Px2同向的位移)及Px1,Py1,Px2' 来表示的刚度矩阵方程。
,,?ux1??ux1??ux1??v??????y1??vy1??vy1?'解：用坐标变换?????T????则????'???T??'?
u?x2??ux2cos???ux2?''??????vusin?uy2x2x??2?????1?0其中?T????0??000100cos?0sin?0?0?? ， 0??0?由?K???'???P???K??T??????P?
?1?10?2.51.832.5??0?1.832.55.0?2.5????0?K??T?????2.54.52.5?2.5??????2.5?2.5?2.52.5??0???10?2.52.964??ux1??Px1???v???P?
104?1.832.52.5???y1??y1??u'??P'????20?x2??0?x2?01000045350??10?2.52.9640?0???2.50????1.832.5?
0???2.54.50.50????2.5?2.5?0.50??0???
10 某平面结构采用四节点矩形单元和三节点三角形单元建立有限元计算模型，其如图3所示。试求结点2的等效荷载列阵
?R2?。
2号节点有关 荷载作用于1?2边上，故等效节点力只与1、解：单元①，
形函数N1???,N2??(1??)，在??1边上，
?N1?N?1,2??1 ?????N?N?N?N?x?y?x?y?l1?02?l,?m1?m2?0则ds?()2?()2d??ld? ????????????????线性分布面力?q????0?? ?q??12??d??0则P1y?N1qyds?qll?ql 3 4
单元③，
?0?ss??形函数N1?1?,N2?,N3?0 在1-2边上，?qs???s?
llq???l?sl0?000??s00??l?T?s0?1?lT?FS???[N]?qS?ds???sll?01??l??ql?10?02?323?00??T?0??s?ds?q??l?
?0???R2??ql??2??故节点2的等效荷载列阵?3??
11 试求如图4所示的有限元网格的整体刚度矩阵，假设每个节点的自由度数为1，且设K表示第e个单元的单元
e刚度矩阵（注意：结果应该用kij表示）。
e
图5
(1)?k11?(1)k??21(1)?k41?(1)??k51(3)?k55?(3)??k75(3)?k85?
(1)k12(1)k22(1)k42(1)k52(3)k57(3)k77(3)k87(1)k14(1)k24(1)k44(1)k54解：单元刚度矩阵?K?(1)(1)(2)??k22k15?(2)(1)?k25(2)?，?K???k32(1)?(2)?k52k45?(2)(1)?k55????k62(4)k56(4)k66(4)k86(2)k23(2)k33(2)k63(2)k53(2)k26(2)k36(2)k66(2)k56(2)?k25(2)?k35? (2)?k65(2)?k55???K?(3)(3)(4)??k55k58(4)?(4)(3)?k78K?，????k65(3)?(4)?k85k88??(4)?k58(4)?k68? (4)?k88? 5