有限元考试试题及答案——第一组.pdf

有限元考试试题及答案1、简答题（5 道，共计 25 分）。1.1.有限单元位移法求解弹性力学问题的基本步骤有哪些？有限单元位移法求解弹性力学问题的基本步骤有哪些？（5 5 分）分）答：（1）选择适当的单元类型将弹性体离散化；（2）建立单元体的位移插值函数；（3）推导单元刚度矩阵；（4）将单元刚度矩阵组装成整体刚度矩阵；（5）代入边界条件和求解。2.在划分网格数相同的情况下，为什么八节点四边形等参数单元精度大于四在划分网格数相同的情况下，为什么八节点四边形等参数单元精度大于四边形矩形单元？边形矩形单元？（5 分）分）答：在对于曲线边界的边界单元，其边界为曲边，八节点四边形等参数单元边上三个节点所确定的抛物线来代替原来的曲线，显然拟合效果比四边形矩形单元的直边好。3.轴对称单元与平面单元有哪些区别？轴对称单元与平面单元有哪些区别？（5 分）分）答：轴对称单元是三角形或四边形截面的空间的环形单元，平面单元是三角形或四边形平面单元；轴对称单元内任意一点有四个应变分量，平面单元内任意一点非零独立应变分量有三个。4.有限元空间问题有哪些特征？有限元空间问题有哪些特征？（5 分）分）答：（1）单元为块体形状。常用单元：四面体单元、长方体单元、直边六面体单元、曲边六面体单元、轴对称单元。（2）结点位移 3 个分量。（3）基本方程比平面问题多。3个平衡方程，6 个几何方程，6 个物理方程。5.简述四节点四边形等参数单元的平面问题分析过程。简述四节点四边形等参数单元的平面问题分析过程。（5）分）分）答：（1）通过整体坐标系和局部坐标系的映射关系得到四节点四边形等参单元的母单元，并选取单元的唯一模式；（2）通过坐标变换和等参元确定平面四节点四边形等参数单元的几何形状和位移模式；（3）将四节点四边形等参数单元的位移模式代入平面问题的几何方程，得到单元应变分量的计算式，再将单元应变代入平面问题的物理方程，得到平面四节点等参 数单元的应力矩阵；（4）用虚功原理求得单元刚度矩阵，最后用高斯积分法计算完成。2、论述题（3 道,共计 30 分）。1.简述四节点四边形等参数单元的平面问题分析过程。简述四节点四边形等参数单元的平面问题分析过程。（10 分）分）答：（1）通过整体坐标系和局部坐标系的映射关系得到四节点四边形等参单元的母单元，并选取单元的唯一模式；（2）通过坐标变换和等参元确定平面四节点四边形等参数单元的几何形状和位移模式；（3）将四节点四边形等参数单元的位移模式代入平面问题的几何方程，得到单元应变分量的计算式，再将单元应变代入平面问题的物理方程，得到平面四节点等参 数单元的应力矩阵；（4）用虚功原理求得单元刚度矩阵，最后用高斯积分法计算完成。2.轴对称问题的简单三角形单元是否是常应力，常应变？为什么？（轴对称问题的简单三角形单元是否是常应力，常应变？为什么？（10 分）分）答：不是常应力和常应变。因为应变与位移分量的关系式为:,这里除含有微分算符外，还包含了 rwu0010rrururzzrrzzrrwzuzw的倒数项 1/r,则即使位移模式为线性的，但由于该项的存在，使得应变与坐标有关，即不会是常应变。应力应变的物理关系为，由于应变不是常应 D变，则所求得的应力也不会是常应力。3.在薄板弯曲理论中做了哪些假设？薄板单元和厚板单元的基本假设有什么在薄板弯曲理论中做了哪些假设？薄板单元和厚板单元的基本假设有什么不同？（不同？（10 分）分）答：四种假设：1）变形前的中面法线在变形后仍为弹性曲面的法线。2）变形前后板的厚度不变。3）板变形时，中面无伸缩。4)板内各水平层间互不挤压。不同点：薄板单元假设横向纤维无挤压，板的中面法线变形后仍保持为直线，该直线垂直于变形后的中面，但是厚板单元的假设考虑横向变形的影响，板的中面法线变形后仍 基本保持为直线，但该直线不再垂直于变形后的中面，法线绕坐标轴的转角不再是挠度的导数，而是独立的变量。3、计算题（3 道，共计 45 分）。1.如图所示等腰直角三角形单元，其厚度为如图所示等腰直角三角形单元，其厚度为，弹性模量为，弹性模量为，泊松比，泊松比；tE0单元的边长及结点编号见图中所示。求单元的边长及结点编号见图中所示。求(1)(1)形函数矩阵形函数矩阵 N(2)(2)应变矩阵应变矩阵和应力矩阵和应力矩阵 BS(3)单元刚度矩阵单元刚度矩阵 （12 分）eK 解：设图 1 所示的各点坐标为点 1（a，0），点 2（a，a），点 3（0，0）于是，可得单元的面积为 12A 2a，及(1)形函数矩阵N为12122121(0aa)a1(00a)a1(aa0)aNxyNxyNxyAA ；123123 NNNNIIINNN(2)应变矩阵B和应力矩阵S分别为12a010-aa-aaB，220010aaa0B，32-a0100a0-aB；123BBBB12a00-aa11-aa22ES，22000aa1a02ES，32-a000a10-a2ES；123123 SD BBBSSS(3)单元刚度矩阵eK111213T2122233132333110211312011110014020200200020111001eEttAKKKKB DBKKKKKK2.如图所示的四结点矩形单元，求出节点如图所示的四结点矩形单元，求出节点 3 的位移。设的位移。设厚度厚度t1m，0，E 为常量。为常量。（13 分）分）注：对于四节点矩形单元有：1141114111411141.14321NNNN )4,3,2,1()1)(1(41iNiii，eTAekkkkkkkkkkkkkkkkyxtBDBk44434241343332312423222114131211dd.2 )4,3,2,1,(3111311a212123111311218dddd21111 jibabbaabEtBDBabtyxtBDBkjijijijijijijijijijijijijTijTAiij解：对于四节点矩形单元有：11411141114111414321NNNN)4,3,2,1()1)(1(41iNiii，eTAekkkkkkkkkkkkkkkkyxtBDBk44434241343332312423222114131211dd )4,3,2,1,(3111311a212123111311218dddd21111 jibabbaabEtBDBabtyxtBDBkjijijijijijijijijijijijijTijTAiij，代入边界条件 1122440，将对应的行和列划掉没剩下的 eeeRk方程为：，PPk-3333又，且，a=1,b=1)4,3,2,1()1)(1(41iNiii1133，所以 411483113112113113112833EEk所以 PPk-3333PPE-4114833解得1158-33EP3.3.有一如图有一如图 3(a)3(a)所示的剪力墙，墙顶作用竖向荷载所示的剪力墙，墙顶作用竖向荷载 P P。将该剪力墙划分为两个。将该剪力墙划分为两个三结点三角形常应力单元，单元和结点编号如图三结点三角形常应力单元，单元和结点编号如图 3(b)3(b)所示，并将荷载所示，并将荷载 P P 分成两分成两个个 P/2P/2 作用在作用在 3 3、4 4 结点。已知单元厚度为结点。已知单元厚度为 t t，弹性模量为，弹性模量为 E E，泊松比，泊松比=1/3=1/3。求结点求结点 3 3 和结点和结点 4 4 的位移，以及单元的位移，以及单元的应变和应力。的应变和应力。（2020 分）分）解：建立直角坐标系(注 Y 轴向下为正)，单元i,j,m 对应的节点编号为 3,1,4，单元对应的节点编号为 2,4,1。对于单元：i(0,0),j(0,4),m(2,0)bi=yi-ym=4；bj=ym-yi=0；bm=yi-yj=-4 ci=xm-xj=2；cj=xi-xm=-2；cm=xj-xi=0 三角形面积 A=1/2*2*4=4几何矩阵B=40000204024202000481弹性矩阵D=2100010112E单元刚度矩阵=tBDBkT1314)1(8001684)1(4)1(8)1(40816048)1(2)1(4)1(2)23(4)1(4)9(2)1(162Et132)1(8001684)1(4)1(8)1(40816048)1(2)1(4)1(2)23(4)1(4)9(212kk然后合成总刚K。整体节点力矢量为02022211PPFFFFFyxyx节点位移矢量为44330000vuvud，采用缩减矩阵法划去位移为零的行与列，得dKF 02/02/)23(4)1(4)1(88)9(2)1(416)23(40)9(2)1(1644332PPvuvuEt解得TEtPvuvu88.199.850.142.84433单元的应变TEtPdB56.047.0285.011单元的应力 TtPdBD84.068.158.011