有限元法基础试题

1、有限元法基础试题（A）一、填空题（5X2分）1.1 单元刚度矩阵ke=hBTDBdC中，矩阵B为,矩阵D为。1.2 边界条件通常有两类。通常发生在位置完全固定不能转动的情况为边界，具体指定有限的非零值位移的情况，如支撑的下沉，称为边界。1.3 内部微元体上外力总虚功：d激=（a,x+7xy,y+Fbx向+（5xWy,y+Fby歹vIdxdy+ffx$u,x+oy6v,y+&y6u,y+§u,x）|dxdy的表达式中，第一项为的虚功，第二项为的虚功。1.4 弹簧单元的位移函数N1+N2=。1.5 kj数学表达式：令dj=,dk=,k#j，则力Fi：%。二、判断题（5X2分）2.

2、1 位移函数的假设合理与否将直接影响到有限元分析的计算精度、效率和可靠性。（）2.2 变形体虚功原理适用于一切结构（一维杆系、二维板、三位块体）、适用于任何力学行为的材料（线性和非线性），是变形体力学的普遍原理。（）2.3 变形体虚功原理要求力系平衡，要求虚位移协调，是在“平衡、协调”前提下功的恒等关系。（）2.4 常应变三角单元中变形矩阵是x或y的函数。（）2.5 对称单元中变形矩阵是x或y的函数。（）三、简答题（26分）3.1 列举有限元法的优点。（8分）3.2 写出有限单元法的分析过程。（8分）3.3 列出3种普通的有限元单元类型。（6分）3.4 简要阐述变形体虚位移原理。（4分）四、计

3、算题（54分）4.1对于下图所示的弹簧组合，单元的弹簧常数为10000N/m,单元的弹簧常数为20000N/m,单元的弹簧常数为10000N/m,确定各节点位移、反力以及单元的单元力。（10分）123450N12344.2对于如图所示的杆组装，弹性模量E为10GPa,杆单元长L均为2m,横截面面积A均为2X10-4m2,弹簧常数为2000kN/m,所受荷载如图。采用直接刚度法确定节点位移、作用力和单元的应力。（10分）5kN1L2L:124.3 对称桁架如图（根据对称性，求图a）所示，杆单元弹性模量均为b）的整体刚度矩阵。（12分）E,横截面面积均为A,单元长度如图,3522P(a)(b)4.

4、4 如图所示的平面桁架，确定转换矩阵,】，并写出工KJIT/(10分)31<4.5 确定下图所示梁的各节点位移。梁已按节点编号离散化。梁在节点1固支，节点2有滚柱支撑，节点3作用有垂直向下的力P=50kN。令沿梁弹性模量E=210GPa,I=12x10-4m4,梁单元长L=3m。弹簧常数k=200kN/m。（12分）LL41P参考答案（A）:一、填空题（5X2分）1.1变形矩阵或应变矩阵弹性矩阵或本构关系矩阵1.2齐次边界非齐次边界1.3微元体上外力在随基点刚体平移所做虚功外力在微元体变形虚位移上所做虚功1.4JL1.5±0二、判断题（5X2分）2.1V2.2V2.3V2.4X

5、2.5V三、简答题（26分）3.1 答：优点有：很容易地模拟不规则形状结构；可以很方便地处理一般荷载条件；由于单元方程是单个建立的，因此可以模拟由几种不同材料构成的物件；可以处理数量不受限制和各类边界条件；单元尺寸大小可以变化；改变模型比较容易可以包括动态作用可以处理大变形和非线性材料带来的非线性问题。（8分）3.2 答：有限元方法的一般步骤有：离散和选择单元类型；选择位移函数；定义应变位移和应力应变关系；推导单元刚度矩阵和方程；组装单元方程得出总体方程并引入边（8分）界条件；求解未知自由度；求解单元应变和位移；解释结果。3.3 答：弹簧单元，杆单元，梁单元，轴对称单元，常应变三角单元，线应变

6、三角形单元，四面体单元等。（任意上述三种均可）（6分）3.4 答：变形体虚位移原理：受给定外力的变形体处于平衡状态的充分、必要条件是，对一（4分）切虚位移，外力所作总虚功恒等于变形体所接受的总虚变形功。四、计算题（54分）4.1解：沿弹簧建立X坐标：（A）每个弹簧单元刚度矩阵如下：W=kej10000一1000011000010000卜2=20000-IL-200002000020000总体刚度矩阵:K=k1k2k3一10000-1000000-1000030000-2000000-2000030000-10000-00-1000010000（B）总体刚度矩阵方程:F1xF2xF3xF4x10

7、000-10000-1000000130000-200000d2x-2000030000-10000d3x0-1000010000_p4x.边界条件：F2x=450N,F3x=0,dix=0,d4x=0解得：d”=0.027m,d3xXXxOX-0.018m,F1x=270N,F4x=180N（C）求单元2节点力?2xI|20000?=2000()-20000Ild2x20000200001d3x：解得:屐=180N,f?x-180N4.2解:沿杆单元建立X坐标:(A)每个单元刚度矩阵如下:km二11061一一1一1N/mk3=21061总体刚度矩阵：-11N/mK=k1k2k3=11-1-1

8、2-1(B)总体刚度矩阵方程:F1xF2x110F3x-1-1-1-1F4x0.0-1-2边界条件:F2x=5000N,解得：d2x-0.003m,d3x.0d2x3x2Jd4xjF3x=0,d1x0,d4x=00.001m,F1x-3000N,F4x-2000N(C)单元的应力K61-1d2x22x>=1父106H$IfZxjULU解得：Rx=2000N,1=-2000N2x3x?仃)=T=4=-10MPa杆单元受压A210有限元法基础试题（B）一、填空题（5X2分）1.1 整体刚度矩阵方程中节点荷载由两部分组成，一是,二是。1.2 常应变三角形单元的位移函数Ni+Nj+Nm=。1.3

9、 最小势能原理与虚位移原理等价，一个是以的形式描述，另一个用的形式表达。1.4 计算轴对称单元刚度矩阵有三种方法，一是采用数值积分，二是,1.5 基本的三维单元是。二、判断题（5X2分）2.1 边界条件通常有两类。通常发生在位置完全固定不能转动的情况为非其次边界。（）2.2 线应变三角形单元中变形矩阵是x或y的函数。（）2.3 杆单元的位移函数N+N2=1。（）2.4 单元刚度矩阵ke=JBTDBdC中，矩阵B为弹性矩阵，矩阵D为变形矩阵。（）2.5 在梁单元中节点力与位移的方向规定应该是与材料力学中规定是一致的。（）三、简答题（26分）3.1 简述刚度矩阵的特性。（6分）3.2 写出位移函数

10、的含义。（4分）3.3 写出推导弹簧单元刚度矩阵的分析过程。（7分）3.4 试列举三种有限元商用软件，并说明各自优点。（9分）四、计算题（54分）4.1 对于下图所示的弹簧组合，单元的弹簧常数为2000N/m,单元的弹簧常数为2000N/m,节点3处位移8为0.01m,确定各节点位移、单元力和反力。（10分）4.2 如图所示的杆单元，杆单元弹性模量为巳杆单元长为L,横截面面积为A,试分别计算（a）、（b）总体x-y坐标下的刚度矩阵。（10分）4.3对称桁架如图（a）所示，杆单元弹性模量均为E,横截面面积均为根据对称性，求图（b）的整体刚度矩阵。（12分）A,单元长度如图,-x4522P1(a)

11、(b)4.4 确定下图所示梁的各节点位移。梁已按节点编号离散化。梁在节点2作用有垂直向下的力P=12kN。令沿梁弹性模量E=70GPa,I=2x10-4m4,梁单元长L=4m。弹簧常数k=200kN/m。（10分）4.5 如图所示梁，确定节点位移，以及每一单元的力和反作用力。梁弹性模量E=70GPa,I=3X10-4m4,梁单元长L=4m。作用在梁单元的均布荷载P为8kN/m。（12分）p参考答案（B）:一、填空题（5X2分）1.1直接节点荷载等效节点荷载1.2J1.3能_1.4直接积分对单元中心点计算1.5四面体单元二、判断题（5X2分）2.1X2.2V2.3V2.4X2.5X三、简答题（2

12、6分）3.1 答：刚度矩阵的特性有：对称的；奇异的；主对角项总是正的。（6分）3.2 答：位移函数的含义：将单元中任意一点的位移近似地表示成该单元节点的函数。当第i个单元自由度为1,而所有其他自由度值为0,Nj代表在整个单元域中假定的位移函数形状。（4分）3.3 答：推导弹簧单元刚度矩阵的分析过程为选择单元类型；选择位移函数；定义应变位移和应力应变关系；推导单元刚度矩阵和方程；组装单元方程得出总体方程并引入边界条件；求解节点位移；求解单元力。（7分）3.4 答：ABAQUS是一套先进的通用有限元系统，属于高端CAE软件。它长于非线性有限元分析，可以分析复杂的固体力学和结构力学系统，特别是能够驾

13、驭非常庞大的复杂问题和模拟高度非线性问题。ABAQUS不但可以做单一零件的力学和多物理场的分析，同时还可以做系统级的分析和研究，其系统级分析的特点相对于其他分析软件来说是独一无二的。ANSYS软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件，发展了很多版本，但是它们核心的计算部分变化不大，只是模块越来越多。ANSYS系统擅长于多物理场和非线性问题的有限元分析，在铁道，建筑和压力容器方面应用较多。LS-DYNA是一个通用显式非线性动力分析有限元程序，最初是1976年在美国劳伦斯利弗莫尔国家实验室由J.O.Hallquist主持开发完成的，主要目的是为核武器的弹头设计提供分析工具，后经多次扩充和改进，计算功能更为强大。LSDYNA长于冲击、接触等非线性动力分析。（9分）四、计算题（54分）4.5常应变三角单元（12分）23aPy02i74-P4.4 如图所示梁，