合肥工业大学力学班有限元试题08A答案

合肥工业大学工程力学专业《工程结构数值方法》2008年试卷A答案2008年5月11日上午考试(共8题，满分100分)，考试时间：120分钟适用专业：工程力学05—1(本科)，出题人：牛忠荣1.(8分)线弹性力学静力学问题有限元法计算列式的理论推导是如何采用弹性力学问题基本方程？答：弹性力学有限元的基本过程是：假设单元的位移场模式{f}[N]{e}代入到几何方程得到{}[B]{e}[Fe][k]{e}[Fe][k]{e}代入到物理方程得到{}[D][B]{」代入到虚功方程，得到单元刚度方程叠加到总刚阵，得到结构的平衡方程引入位移边界条件后，解第5步得到的方程组，可以得到结点位移2・(12分)用有限元法计算图示平面刚架时，[Fe][k]{e}[Fe][k]{e}如何进行结点编号，使系统总刚度矩阵K的带宽最小？计算在你的结点编号下的系统总刚阵K的半带宽。在结点编号确定后，按此顺序进行自由度编号，则图中A结点水平位移对K中的列和行位置是多少？那些单元对该元素应的主对角元素在总刚阵的数值有贡献K中的列和行位置是多少？那些单元对该元素争：Iy91317题2图见图；半带宽=3(41)15;(3)A结点水平位移对应的主对角元素是k28,28；A点周围的4个单元对k28,28有贡献。Iy91317题2图（10分）弹性力学有限元中，平面等参数单元中的“等参数”概念是何意思？该单元在跨相邻单元时，位移场连续吗？应力场连续吗？答：在有限单元法中最普遍采用的是等参变换，即单元几何形状的变换和单元内的场函数采用相同数目的节点参数及相同的插值函数进行变换。采用等参变换的单元称之为等参元。所谓“等参元”是指几何形状插值形函数和单元上的位移插值形函数相同，参数个数相等。相邻等参元之间，位移场是连续的，应力场不连续。（10分）回答下列问题：（1）弹性力学平面问题8节点等参元，其单元自由度是多少个？单元刚阵元素是多少个？（2）弹性力学空间问题20节点等参元，其单元自由度是多少个？单元刚阵元素是多少个？（3）一般的平面刚架结构梁单元的自由度是多少个？单元刚阵元素是多少个？一（4）一般的空间刚架结构梁单元的自由度是多少个？单元刚阵元素是多少个？一答：1）自由度16个，刚阵元素16X16=256；2）自由度60个，刚阵元素60X60=3600；3）自由度6个，刚阵元素6X6=36；4）自由度12个，刚阵元素12X12=144.（10分）结构振动问题有限元离散的无阻尼自由振动方程为（K2m冲0式中Knn是刚度矩阵，Mnn是质量矩阵，是结构固有频率，中是振型向量。试问为什么从上式求出的特征对<“Q>（i1,2,,n）中，只有前若干低阶频率和相应振型是可靠的，误差较小。答：在有限单元法中，采用低阶多项式拟合振型。结构的低阶振型曲线与低阶多项式比较通配，结构的高阶振型曲线与低阶多项式曲线有着显著的差异。因而，有限元法中求出的低阶频率和振型是可信的，而所求出的高阶频率和振型误差较大，甚至无效。（7分）用有限元法求结构的动力响应通常有振型叠加法和逐步积分法，试问该两种方法的优缺点。答：振型叠加法适应于求解线性振动问题，需要计算结构的固有频率和振型。适合于时间较长求解。逐步积分法可用于非线性振动问题的求解，不必需要计算结构的固有频率和振型。适合于时间较短的问题求解，当时间长时计算量大。

7・(25分)图示矩形薄板2m1m，厚度h0.1m。载荷及约束信息如图示，泊松比0,EE0N/m2,p10N，板的重量密度为30N/m3,沿y轴负向。试采用三角形常应变元求：(1)建立适宜的有限元计算模型；(2)结构系统的等效结点力列阵；米用划行划列法引入已知结点位移，计算出结点3和4的位移；单元1和2的各应力分量。(建议单元编号：ijm单元1:342单元2:213)1mP53「4164题7图1m1m1m解：计算模型P/2312-21m题7图单元1的等效结点力:结点311F-3hA(010单元2的等效结点力:T1)0.5(0结点21(010F3hA结构整体等效结点力结点1F(0T1)(0单元1:.,biciyjymxm单元2:cixjx1应力矩阵结点34200101E0010100.500.50.5005.,bi0单元S1ym0.5xjyj0,bj1,c0,bj1,c单元2应力矩阵2130010100101000.500.50.500.5结点,S2E0单元1刚度矩阵:kii，bm1,cm,bm1,cm05011.,biyjym,bj,bmcixjxm1,cj1,cm00.50E0h0.50.5Eoh00.5,kij,kim—01201200kE0h1.50.5，kjmEoh10.5jj20.51.5200.5kEoh_10mm200.5单元1:E0h20.500.50.500.53010100Ehk1.50.510.5420.51.500.5010200.53422130.500.50.500.52010100Eh2101.50.510.51kk20.51.500.5010300.5叠加总刚阵：结点号：1234201已知结点位移：U]U2U3U4E20已知结点位移：U]U2U3U4E0h21111111i111U1V1U200.501.511111•-・0.51V2=1•,00.5U30011.501—L—1V36|i11.50.5u40|i10.51.50.5V]V2021.500.50.5301.5011.50.540.51.5米用划行划列，得到1.5v3V420v31.5V4解得结点位移：v3152/E0代入：得单元应力（常应力）：(―05)246N/m，xy54xy8・（18分）图示正方形薄板，厚度h108/E0。540.1m，边长a1m。约束信息如图示，板的质量密度为30kg/m3，泊松比0，弹性模量E2N/m2，取g10m/s2。不计系统阻尼，试采用图示的2个三角形常应变元求:（1）引入已知结点位移后（2）结构自由振动方程；,结构系统的刚度矩阵和一致质量矩阵;y1m（3）计算结构的固有频率1m解：0.50.50.50.5h0.50.5M11200.5100.50.500.50100.500.50100.500.5020100.500.5M2_h0.50100.501；1200.50100.530.500.501000.500.5011MM1M2hI11111I11111111111ii212---20…0.5"=“030200.5TOC\o"1-5"\h\z0.5010400.501aIKE0K1.5KE0K1.5111.512011.51011.520.5120.5M———1200.