1第三章1如图所示一三角形钢板,两个结点固定,对第三个结点施以单位水平位移,测出所施加的力,从而得出相应的刚度系数。其他点依此类推,这样测得的刚度系数所组成的刚度矩阵,是否与按照常规三角形单元刚度矩阵计算公式所得结果一样?用这样实测所得的刚度矩阵能否进行有限元分析?为什么?解:不一样。单元刚度矩阵中每个元素的物理意义:ijk表示单元第j个自由度产生单位位移,其它自由度固定时,第i个自由度产生的节点力。单元刚度矩阵是在单元处于平衡状态的前提下得出的,单元作为分离体看待,作用在它上面的外力(单元力)必是平衡力系,然而研究单元平衡时没有引入约束承受平衡力系作用的无约束单元,其变形是确定,但位移是不能确定的,即单元可发生任意的刚体位移。不能。因为与有限元中单元与单元之间的约束情况不一样,不能进行有限元分析。2以位移为基本未知量的有限元法其解具有下限性质,试证明之。解:系统总位能的离散形式12TTpaKaaP将求解的方程KaP带入可得1122TTTpaKaaKaaKaU在平衡情况下,系统总位能等于负的应变能。在有限元解中,由于假定的近似位移模式一般来说总与精确解有差别的。设近似解为p、U、[]K、a、KaP,真实解为p、U、[]K、a、KaP且根据最小势能原理,得到的系统的总位能总会比真正的总位能要大,故pp则UUTTTTaKaaKaaPaP则近似解的位移总体上小于精确解的位移解释如下:单元原是连续体的一部分,具有无限多个自由度,在假定了单元的位移函数后,自由度限制为只有以结点位移表示的有限自由度,引入了更多的约束和限制,使得单元刚度较实际连续体加强了,连续体的整体刚度随之增加,所以有限元解整体上较真实解偏小。3请分别阐述单元刚度矩阵和整体刚度矩阵中任一元素的物理意义。解:在单刚eK中,eijk表示单元第j个位移产生一单位位移,其它位移为零时,第i个位移方向上引起的节点力。在整体刚度中,ijK表示第j个自由度产生一单位位移,其它自由度为零时,第i个自由度上引起的节点力。24简述虚功原理,且使用虚功原理导出外荷载与节点荷载的等效关系式。解:虚功原理:变形体中任意满足平衡的力系在任意满足协调条件的变形状态上作的虚功等于零,即体系外力的虚功与内力的虚功之和等于零。设eq为外荷载(此处为体力),ep为节点荷载,ew为单元内位移场,e为结点位移场根据虚功原理TTeeeeVpwqdV由于[]eewN故[][]TTTeeeeeeTTVVVwqdVNqdVNqdV则[][]TTeeeeeeTTVVpNqdVpNqdV5试述弹性力学中按位移求解与有限单元法中按位移求解之间的异同点。解:弹性力学有限单元法物理模型连续体离散化结构基本方程几何方程物理方程平衡微分方程几何方程物理方程结点平衡方程解法解微分方程解代数方程解答形式用函数表示用数值表示解答精度精确解近似解6如果三节点三角形单元绕其中某一个节点作小的刚体转动,其转角为,证明单元内所有的应力均为零。解:在三角形单元中DB0000001100000022ijmjmmiijijmjmmiijiijjmmjmjmmimiijijbbbyyyyyyBcccxxxxxxAAcbcbcbxxyyxxyyxxyy由于三角形单元绕其中某一个节点作小的刚体转动,各节点的位移可表示为:00uuyvvx则可知节点位移向量0,0,,,,TTjjmmyxyx3故应变000000110000220immiijjjmmiijjjmimmimiijijmmyyyyyyyBxxxxxxxAAxxyyxxyyxxyyyx由于弹性矩阵D为常量矩阵,应变向量为零向量,故D为零向量,即单元内所有的应力为零。7二维单元在,xy坐标内平面平移到不同位置,单元刚度矩阵相同吗?在平面内旋转时又怎样?试证明之。解:二维单元在,xy坐标内平面移到不同位置时,刚度矩阵相同。在平面内旋转时,刚度矩阵也相同。刚度矩阵21122114(1)22rsrsrsrsTrsrsrsrsrsrsbbccbccbEhAkBDBhAcbbcccbb单元平移或旋转时,,iibc不变,故单元刚度矩阵不变。8判断有限元网格离散合理性a)对图1(a)所示的有限元网格,评论网格的优劣性,指出模型中的错误,并加以改正。b)评论图1(b)的网格划分合理吗?为什么?请加以改正。图1解:(a)网格划分不合理。1)无过渡单元2)无边界条件3)夹角区应力集中,应适当加密风格4)对称结构网格应对称划分(b)不合理。1)左部网格应适当加密2)由于三角形单元会造成局部精度不够,过渡区可采用其它单元划分3)右部单元的长宽比较大,就进行适当调整。9如图2所示,平面三角形构件以x-y坐标系表示的刚度矩阵方程如下:2211221145.25.25.25.25.25.25.45.25.20.55.283.15.283.15.21010yxyxyxyxPPPPvuvu4试建立以1xu,1yu,,2xu(与图中,2xP同向的位移)及1xP,1yP,'2xP来表示的刚度矩阵方程。解:用坐标变换'T则111111''222''222cossinxxxyyyxxxyxxuuuvvvTuuuvuu其中1000010000cos000sin0T,由'KPKTP1000102.51.832.5102.52.964001001.832.55.02.51.832.52.5040002.54.52.52.52.54.50.5052.52.52.52.532.52.50.500005KT11411''22102.52.964101.832.52.5020xxyyxxuPvPuP10某平面结构采用四节点矩形单元和三节点三角形单元建立有限元计算模型,其如图3所示。试求结点2的等效荷载列阵2R。荷载作用于12边上,故等效节点力只与12、号节点有关解:单元①,形函数12,(1)NN,在1边上,121,1NN12120,0NNNNxyllmm则22()()xydsdld线性分布面力0qq则121103yylqlPNqdsqld5图5形函数1231,,0ssNNNll在1-2边上,0sqsql单元③,010000[]01000120000233TTSSllTssllFNqdsdssssqlllql故节点2的等效荷载列阵2023Rql11试求如图4所示的有限元网格的整体刚度矩阵,假设每个节点的自由度数为1,且设eK表示第e个单元的单元刚度矩阵(注意:结果应该用eijk表示)。解:单元刚度矩阵(1)(1)(1)(1)11121415(1)(1)(1)(1)(1)21222425(1)(1)(1)(1)41424445(1)(1)(1)(1)51525455kkkkkkkkKkkkkkkkk,(2)(2)(2)(2)22232625(2)(2)(2)(2)(2)32333635(2)(2)(2)(2)52636665(2)(2)(2)(2)62535655kkkkkkkkKkkkkkkkk(3)(3)(3)555758(3)(3)(3)(3)757778(3)(3)(3)858788kkkKkkkkkk,(4)(4)(4)555658(4)(4)(4)(4)656668(4)(4)(4)858688kkkKkkkkkk图46整体刚度矩阵:(1)(1)(1)(1)11121415(1)(1)(2)(2)(1)(1)(2)(2)2122222324252526(2)(2)(2)(2)32333536(1)(1)(1)(1)41424445(1)(1)(2)(2)(1)(1)(2)(3)(4)(5152525354555555555600000000000000kkkkkkkkkkkkkkkkkkkkKkkkkkkkkkk2)(4)(3)(3)(4)56575858(2)(2)(2)(4)(2)(4)(4)62636565666668(3)(3)(3)757778(3)(4)(4)(3)(3)(4)858586878888000000000000kkkkkkkkkkkkkkkkkkkk12图5中两个三角形单元组成平行四边形,已知单元①按局部编码,,ijm的单元刚度矩阵K①和应力矩阵S①是(1)806626166126413.597.5313.531.59.535.5K对称(1)003030040301201.5-1.5-0.51.5S---按图5示单元②的局部编码写出K②,S②。解:由图可知(1)(2)(1)(2)(1)(2),,miijjm则由iiijimjjjmijjjjmmmmimmiiKKKKKKKKKKKK得到(2)9.53267.535.56431.580661661213.5913.5K对称(2)3000-300-1040-3-0.51.520-1.5-1.5S13如图6所示8结点矩形单元(每边中点为结点),3点为坐标原点,a=b=2,单元厚为t。①求该单元的位移函数和形函数和并检验其是否满足收敛性条件。②求在2-6-3边作用均布水平荷载q时的等效结点荷载。解:(1)位移函数:2222123456782222910111213141516uxyxxyyxyxyvxyxxyyxyxy引入无量纲的局部坐标,xyabbaxy37481526q7则1313223,,22xxyynxy故123123110,,1,0,,12212312311112()(1),4(1),2(),2()(1),4(1),2()2222lllppp则2n时,12120,1,0,112121,,1,l