东南大学计算力学习题及答案汇总版文档格式.docx

1、kargJk 心？= falpliaLp?则近似解的位移总体上小于精确解的位移解释如下：单元原是连续体的一部分，具有无限多个自由度，在假定了单元的位移函数后，自由度限制为只有以结 点位移表示的有限自由度，弓I入了更多的约束和限制，使得单元刚度较实际连续体加强了，连续体的整体刚度随之 增加，所以有限元解整体上较真实解偏小。3请分别阐述单元刚度矩阵和整体刚度矩阵中任一元素的物理意义。在单刚l-K e中，坏表示单元第j个位移产生一单位位移，其它位移为零时，第 i个位移方向上引起的节点力。在整体刚度中，Kj表示第j个自由度产生一单位位移，其它自由度为零时，第 i个自由度上引起的节点力。4简述虚功原理，

2、且使用虚功原理导出外荷载与节点荷载的等效关系式。虚功原理：变形体中任意满足平衡的力系在任意满足协调条件的变形状态上作的虚功等于零，即体系外力的虚 功与内力的虚功之和等于零。设g f为外荷载（此处为体力），：pf为节点荷载，twf为单元内位移场，匕为结点位移场T T根据虚功原理：e p -w ： q fdVV由于we =N佃e 故 jwefqedV = j3eTNTqedV =eT JNTqedVV V V则門pe=汀 JNt qedVn pe = JNT qe dVV V5试述弹性力学中按位移求解与有限单元法中按位移求解之间的异同点。弹性力学有限单元法物理模型连续体离散化结构基本方程几何方程物

3、理方程平衡微分方程结点平衡方程解法解微分方程解代数方程解答形式用函数表示用数值表示解答精度精确解近似解6如果三节点三角形单元绕其中某一个节点作小的刚体转动，其转角为 二，证明单元内所有的应力均为零。在三角形单元中 士T - id HbiIJ由于三角形单元绕其中某一个节点作小的刚体转动，各节点的位移可表示为:则可知节点位移向量 J = ：0,0, -*j, Xj , -VymJXm由于弹性矩阵D 为常量矩阵，应变向量:为零向量，故CdK；为零向量，即单元内所有的应力为零。7二维单元在X,y坐标内平面平移到不同位置，单元刚度矩阵相同吗？在平面内旋转时又怎样？试证明之。 解：二维单元在x, y坐标内

4、平面移到不同位置时，单元平移或旋转时， b，G不变，故单元刚度矩阵不变。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

5、 1Ux1X1 1“亠42.55.0-2.5Vy1Py110Px24.5Ux2一 2.52.5 一W 一1 1Py2 一Ux1 用坐标变换芳 = T 佃则丿= “vy1y1卜Ux2 cosot |2r i企 si net J101其中T =cosa-si na0_试建立以Um Uyi, Ux2（与图中巳2侗向的位移）及PXl , Py1, Px2来表示的刚度矩阵方程。由 k IKPIK It =10 0 00 10 0 10 -2.5 2.964 041.83 2.5 2.5 00 0工0=5-2.5 4.5 0.5 032.5 -2.5 -0.5 0一0 0 0L5 一10 1.83.03所

6、示。试求结点2的等效-2.5 2.964 Ux1 PX12.5 2.5 皿十120 j|uxj Px2”10某平面结构采用四节点矩形单元和三节点三角形单元建立有限元计算模型，其如图 荷载列阵 *2：?o单元， 荷载作用于1-2边上，故等效节点力只与 1、号节点有关则 Py = jNydsuql j2d匕=l 0 3丨在1-2边上，傀.；工siqrjpl故节点2的等效荷载列阵R2，ql 213J1,且设Ke表示第e个单元的单元11试求如图4所示的有限元网格的整体刚度矩阵，假设每个节点的自由度数为刚度矩阵（注意：结果应该用 k：表示）。单元刚度矩阵K -f8-6-2616-12-4-30、K（1）

7、 =对13.59-7.5s（1） =-1-1.5.2-0.51.59.5称5.5 j12图5中两个三角形单元组成平行四边形，已知单元按局部编码i, j,m的单元刚度矩阵按图5示单元的局部编码写出 K，S。K和应力矩阵S是由图可知 m i（2） ,i（1） j ,j（1） m（2）S（2）=1-0.51.513如图6 所示8结点矩形单元（每边中点为结点）,3点为坐标原点, 元厚为t。1求该单元的位移函数和形函数和并检验其是否满足收敛性条件。2求在2-6-3边作用均布水平荷载 q时的等效结点荷载。（1）位移函数：-L u =5 Lx ： 3y ：jx2 ： 5xy ： 6y2 ： 7x2y ： 8

8、xy2 2 22 2V ToX 11y ： 12X ： 13xy ： 14y ： 15X y 16xy引入无量纲的局部坐标X1 X3 屮 y则 n = 3, X2 ,、2 ：2 2故 1=0, 2 , 3=1, 1 =0, 2 :11 =2（ -1）（3=1-1）,12 一（ -1）3=2（ _2）,p2（ -）（ -1），p-4（ -1），p32 （-）则 n =2时， 0, 2 =1, 4=0, 2 =1 h =1 -2 二，P1 =1 -，P2 虫：则角节点的形函数为11 “11N C -）（ -）,N4 C -）（ -）C -1）22 2 2N3=4（ -2）（ 一1）（ _才）（_1

9、）,N4 =4 （ ）（ -1）（ -1）边中节点的形函数为N5 =一4 （ -1）,N -4 （ 一1）（1一 ）, N7 =-4 （ 一1）（1 一 ）,N8 = -4 （ -1）证明收敛性：位移函数中厂 2 2 2 2+a5xy+G6y +a7xy+a8Xy2 2 2 2丫=。9+。10*+口11丫+0（12乂 +ct13xy+c（14y +a15xy+c（16xy1, -S表示刚体位移，-：2,二3和二9,二10表示常应变，故位移函数具有完备性设相邻单元公共边界上的直线方程是 y = b （或x = a），代入位移函数中r 2 2 2u 予 +a3b +6b +（a2 +5b +8b

10、）x +（o4 +a7b）x2 2 2y =C（9 +c（11b +c（14b +（O（10 +c（13b） +c（16b X + （。伐 + Gb）X为x （或y）的2次函数，而边界上三点确定的位移函数为也为二次曲线，故单元在公共边界连续， 故位移函数收敛（2）荷载作用在2 -3边上，故等效节点力只与 2,6,3号节点有关11 1 1N2 =4 （ -一）（ 一一）（ 一 1）, N6=4（ - 1）（1一），N3=4（ -一）（ 一 1）（ 一一）一1）N 2 ?N6 ：N3亠 4 4（2）y2=0，仝=0匹卫匹ds 二（；X）2 （；y）2d =2d乐=t. Mqxds =2qt .弘在 -qt Bx =t . Wqxds