北京科技大学有限元最近3年有限元真题

1、北京科技大学2009 2010学年硕士研究生“工程中的有限元方法”试题姓名学号班级成绩说明：1-5题为笔试题，每题 10分。上机题结合实验报告共50分。1、简述弹性力学四边形四节点等参元的收敛性质以及由该单元刚度矩阵装配成的总刚度矩阵的性质。3在单元分析已经提出有限单元解的收敛性要求，即，单元必须是完备的和协调的。对于等参单元：1. 完备性：对于C0型单元，由于等参单元的形函数中包含有常数项和线性项，满足完备性的要求。2. 协调性：由于单元之间的公共边上有完全相同的节点，同时每一单元沿这些边的坐标和未知函数均采用相同的插值函数加以确定。因此，只要在划分网格时，遵守单元选择和节点配置的要求，则等

2、参单元满足协调性的要 求。2总刚的性质1）对称性2）奇异性，需引入合适的位移约束。3）稀疏，（存在许多零元素）4）非零元素呈带状分布5）主元恒正根据物理意义可得此性质，正常情况下，主元占优2、分析图示的两个单元在什么条件下其连接关系正确。要求说明所采用单元的类型和连接方法。采用四边形等参元附加多点约束方程过渡。4边形5节点Serendipity 过渡单元约束方程：U6=（U 2+U3）/2卜八LJ1盼 4 510698713、对于右图所示三节点网格，设每个节点具有一个自由度。其： 最大带宽=（9-1）*仁8最大波阵宽=31,2,109,2,109.2.39.8.34.8.34、某非协调板单元，

3、单元长度为2 2，节点基本未知量为： 在图示的坐标系下，其关于 w的插值函数形式为： 其中：1试：（1咕说明此插值函数属于哪一族插值函数? 0（ 插；值函数及其导函数均1有2古1的（性质+1），式中：0i， 5、三维实体元 如图所示。根据已知形函数，写出它们在图示坐标系下是哪个节点 形函数，依此规律写出节点 7、 18的形函数形式。22 23N7=1/64-10+9(xA2+yA2 ) ( 1+x) (1+y)(1-z)N18=9/64(1-xA2)(1+9z)(1+x)(1+1/3y)北京科技大学2010 2011学年硕士研究生“工程中的有限元方法”试题姓名学号班级成绩说明：17题为笔试题，

4、共50分。上机题结合实验报告共50分。3、等参元有何特点？在四边形等参元网格的划分中，为什么要保证其内角小于180。？ （8分）等参数单元（简称等参元）就是对单元几何形状和单元内的参变量函数采用相同数目的节点参数和相同的形函数 进行变换而设计出的一种新型单元。优点：由于等参变换的采用使等参单元的刚度、质量、阻尼、荷载等特性矩阵的计算仍在前面所表示单元的规则 域内进行，因此不管各个积分形式的矩阵表示的被积函数如何复杂，仍然可以方便地采用标准化的数值积分方法 计算。也正因为如此，等参元已成为有限元法中应用最为广泛的单元形式。具有计算精度高和适应性好的特点， 是有限元程序中主要采取的单元形式。为了保

5、证等参变换的一一对应性质，应当避免单元一边上的两点退化成一个节点，还要防止单元的任意两边的夹角接近180度。更不允许夹角等于或大于180度。X；一族插值函数?,；qn宓4分）2、在图1所示的坐标系下，某插值函数形式为: 其中：试：（1）（2）此插值函数具有什么基本性质？（ Hermite插值函数，满足这有些实际的插值问题不但要求在节点上函数值相等，有些实际的插值问题不但要求在节点上函数值相等， 种要求的插值多项式就是 埃尔米特插值多项式3、构造一个四边形 5节点Sere ndipity单元（长度：2 x 2）,其中，5节点为等距离边内点。要求写出单元的1节点和5节点所对应的形函数的具体形式。（

6、8分） 已知：4结点矩形单元的插值函数：4、三维四面体10节点元如图所示。补全已给形函数中的下标，以表明它们在图示坐标系下n哪个节点形函数。 （6分）对于图示的单元，其形函数为：试证明在这样两个单元的公共边上不满足位移协调条件。（1 + n ）6分） 5、对于右图所示I四边形 4节点网格,2设每个节点具有 2个自由度。分）3层为隔热材料，外 热系数不同。假定烟囱2内表面的温度为 100 c,烟囱外表面暴露在温度为 为h。要求利用对1称性画出有限元网；一、简要回答下列各题（每小题 7分）二、1）试叙述有限元分析的主要步骤。三、它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个

7、合适的 个域总的满足条件（如结构的平衡条件），从而得到问题的解。在选取单元位移函数时，应遵循哪些原则？ 位移形状函数有哪些性质？4节点网格，设每个节点具有 . A ,最大波阵宽=_7、正方形截面的烟囱如图所示，烟囱壁壁由两层45、对于右图所示I四边形 网格的最大带宽生_料构成。n6匚层为钢筋混凝土。内10 C的在空气中。 什么边界条件？（9的导外外层烟囱与空气换热系数（较简单的）（近似解，然后推导求解这四、五、六、4）等参元有何特点？七、5 ）在四边形单元的划分中，为什么要保证其内角小于180。？八、6）对单元的节点编号时，为什么要使同一单元的节点号差最小？九、二、如图所示的结构已划分为 3个

8、3节点三角形单元，设各节点的自由度为1，各单元的十、刚度矩阵相同，先定义每个单元的节点号，然后求出结构的刚度矩阵。（15分）十一、单元刚度矩阵为：（e =1，2，3）十二、三、对于图示的两单元，其位移函为，试证明在单元的公共边上不满足位移协调条件。（10分）十三、 四、图为8节点四边形单元（在局部坐标系中进行分析），先选取适当的位移函数，然十四、后构造岀各节点的形函数。（15分）（已知拉格朗日一维插值函数：十五、2点插值十六、3点插值）十七、十八、五、在对结构进行离散化处理时，若单元划分得不合理将会大大降低计算结果的精度，甚十九、至产生错误的结果。试列举岀 3种不合理的单元划分，并指岀不合理的

9、原因北京科技大学20122013学年硕士研究生“工程中的有限元方法”试题一.是非题(认为该题正确，在括号中打 V;该题错误，在括号中打 X。)(每小题2分)(1) 用加权余量法求解微分方程，其权函数V和场函数u的选择没有任何限制。()(2) 四结点四边形等参单元的位移插值函数是坐标x、y的一次函数。()(3 )在三角形单元中，其面积坐标的值与三结点三角形单元的结点形函数值相等。()(4) 二维弹性力学问题的有限元法求解，其收敛准则要求试探位移函数C1连续。()(5) 有限元位移法求得的应力结果通常比应变结果精度低。()(6) 等参单元中Jacobi行列式的值不能等于零。()(7) 在位移型有限

10、元中，单元交界面上的应力是严格满足平衡条件的。()(8) 四边形单元的Jacobi行列式是常数。()(9) 利用高斯点的应力进行应力精度的改善时，可以采用与位移插值函数不同结点的形函数进行应力插值。()(10) 维变带宽存储通常比二维等带宽存储更节省存储量。()二单项选择题(共20分，每小题2分)1在加权余量法中，若简单地利用近似解的试探函数序列作为权函数，这类方法称为。(A)配点法(B)子域法(C)伽辽金法2等参变换是指单元坐标变换和函数插值采用 的结点和的插值函数。(A)不相同，不相同(B)相同，相同(C)相同，不相同(D)不相同，相同3有限元位移模式中，广义坐标的个数应与 相等。(A)单

11、元结点个数(B)单元结点自由度数(C)场变量个数4采用位移元计算得到应力近似解与精确解相比较，一般 。(A)近似解总小于精确解(B)近似解总大于精确解(C)近似解在精确解上下震荡(D)没有规律5如果出现在泛函中场函数的最高阶导数是m阶，单元的完备性是指试探函数必须至少是 完全多项式。(A) m1 次(B) m次 (C) 2m1 次6与高斯消去法相比，高斯约当消去法将系数矩阵化成了 形式，因此，不用进行回代计算。(A)上三角矩阵(B)下三角矩阵(C)对角矩阵7对称荷载在对称面上引起的 分量为零。(A)对称应力(B)反对称应力(C)对称位移(D)反对称位移8对分析物体划分好单元后，会对刚度矩阵的半

12、带宽产生影响。(A)单兀编号(B)单兀组集次序(C)结点编号9 n个积分点的高斯积分的精度可达到 阶。(A) n-1( B) n (C) 2n-1(D) 2n10引入位移边界条件是为了消除有限元整体刚度矩阵K的。(A)对称性(B)稀疏性(C)奇异性三. 简答题(共20分，每题5分)1、简述有限单元法结构刚度矩阵的特点。2、简述有限元法中选取单元位移函数(多项式)的一般原则。3、简述有限单元法的收敛性准则。4、 考虑下列三种改善应力结果的方法(1)总体应力磨平、(2)单元应力磨平和(3)分片应力磨平，请分别将 它们按计算精度(高 低)和计算速度(快 慢)进行排序。四. 计算题(共40分，每题20

13、分)1、如图1所示等腰直角三角形单元，其厚度为t，弹性模量为E，泊松比V-0 ;单元的边长及结点编号见图中所示。求形函数矩阵N（2）应变矩阵B和应力矩阵S单元刚度矩阵Ke2、图2 （a）所示为正方形薄板，其板厚度为t，四边受到均匀荷载的作用，荷载集度为1N/m2，同时在y方向相应的两顶点处分别承受大小为2N/m且沿板厚度方向均匀分布的荷载作用。设薄板材料的弹性模量为E，泊松比 0。试求（1）利用对称性，取图（b）所示1/4结构作为研究对象，并将其划分为4个面积大小相等、形状相同的直角三角形单元。给出可供有限元分析的计算模型（即根据对称性条件，在图（b）中添加适当的约束和荷载，并进行单元编号和结

14、点编号）。（2）设单元结点的局部编号分别为i、j、m，为使每个单元刚度矩阵 Ke相同，试在图（b）中正确标出每个单元的合理局部编号；并求单元刚度矩阵Ke。（3）计算等效结点荷载。（4）应用适当的位移约束之后，给出可供求解的整体平衡方程（不需要求解）。.是非题（认为该题正确,華括号中打 V图该题错误，在括号中打 X。）（每小题2分） x V.单项选择题（共C B吊C B2N/mV xC C、题5分D C3 ）主对角元恒正;V x x V V 小题2 分)1、答：（答（1）2、答：t一般原则有非零元素带状分布(1)广义坐标的个数应该与结点自由度数相等； 选取多项式时，常数项和坐标的一次项必须完备;

15、多项式的选取应由低阶到高阶；尽量选取完全多项式以提高单元的精度。(b)（4）3、答：完备性要求，协调性要求 具体阐述内容4、答：计算精度（1） （3）计算速度（2） （3）四.计算题1、解:设图1所示的各点坐标为点 1 （a,(2分)(3分)于是，可得单元的面积为0），点 2（ a，a），点Ala2，及23（0，0）(1) 形函数矩阵N为a1 2(0 a(7分)N1N1ax - ay)0Lx ay);N =IN1=N1I n2 I n3 1N 2 N 3 1简答题（共或0分，每 对前I外给4分） 寸称性；（2）奇异性;二 2 (a2 -ax 0Ly)a应变矩阵B和应力矩阵S分别为-a , B 2 =a(712a1 a , B ;-aB 2B3 Ja-a-a，S 2 二E一 2aS E，S 3 _ 2aa1 a -=aia!a1-1 a22