有限元法基础试题

1、有限元法基础试题一、填空题（5X 2分）单元刚度矩阵keBTDBd 中，矩阵B为，矩阵D为。边界条件通常有两类。 通常发生在位置完全固定不能转动的情况为 边界， 具体指定有限的非零值位移的情况，如支撑的下沉，称为 边界。内部微元体上外力总虚功：d Wex,x xy,yFbxu xy,x y,yFbyv dxdy+ xu,xyv,yxyu,yu,xdxdy 的表达式中， 第一项为 的虚功，第二项为 的虚功。弹簧单元的位移函数 N1+ N2 = 。kj数学表达式：令 dj =, dk=, k j，则力Fi kj。二、判断题（5X 2分）位移函数的假设合理与否将直接影响到有限元分析的计算精度、效率和

2、可靠性。（）变形体虚功原理适用于一切结构（一维杆系、二维板、三位块体）、适用于任何力学行为的材料（线性 和非线性），是变形体力学的普遍原理。（）变形体虚功原理要求力系平衡，要求虚位移协调，是在“平衡、协调”前提下功的恒等关系。（）常应变三角单元中变形矩阵是x或y的函数。（）对称单元中变形矩阵是x或y的函数。（）三、简答题（ 26 分）列举有限元法的优点。（ 8 分） 写出有限单元法的分析过程。（ 8 分） 列出 3 种普通的有限元单元类型。（ 6 分） 简要阐述变形体虚位移原理。（ 4 分）四、计算题（ 54 分）对于下图所示的弹簧组合，单元的弹簧常数为10000N/m,单元的弹簧常数为 20

3、000N/m,单元的弹簧常数为10000N/m,确定各节点位移、反力以及单元的单元力。（10分）对于如图所示的杆组装，弹性模量E为10GPa，杆单元长L均为2m，横截面面积 A均为2 X 10-4m2，弹簧常数为2000kN/m，所受荷载如图。采用直接刚度法确定节点位移、作用力和单元的应力。（10 分）对称桁架如图（a）所示，杆单元弹性模量均为E，横截面面积均为 A，单元长度如图，根据对称性，求图（b）的整体刚度矩阵。（12分）（ a）（ b）如图所示的平面桁架，确定转换矩阵T1 ，并写出 T1 K T1 （10 分）确定下图所示梁的各节点位移。梁已按节点编号离散化。梁在节点1 固支，节点 2

4、有滚柱支撑，节点 3作用有垂直向下的力P=50kN。令沿梁弹性模量 E=210GPa ,1=12 X 10-4m4，梁单元长 L=3m。弹簧常数k=200kN/m 。（12分） 参考答案 （A ）：一、填空题 （ 5X 2 分）变形矩阵或应变矩阵弹性矩阵或本构关系矩阵齐次边界非齐次边界微元体上外力在随基点刚体平移所做虚功外力在微元体变形虚位移上所做虚功11/0二、判断题（5 X 2分）V "" X "三、简答题（26分）答：优点有：很容易地模拟不规则形状结构；可以很方便地处理一般荷载条件；由于单元方程是 单个建立的，因此可以模拟由几种不同材料构成的物件；可以处理数

5、量不受限制和各类边界条件； 单元尺寸大小可以变化；改变模型比较容易可以包括动态作用可以处理大变形和非线性材料带来 的非线性问题。（8分）答：有限元方法的一般步骤有：离散和选择单元类型；选择位移函数；定义应变位移和应力应变关系；推导单元刚度矩阵和方程；组装单元方程得岀总体方程并引入边界条件；求解未知自由度;求解单元应变和位移；解释结果。（8分）答：弹簧单元，杆单元，梁单元，轴对称单元，常应变三角单元，线应变三角形单元，四面体单元等。（任意上述三种均可）（6分）答：变形体虚位移原理：受给定外力的变形体处于平衡状态的充分、必要条件是，对一切虚位移，外力 所作总虚功恒等于变形体所接受的总虚变形功。（4

6、分）四、计算题（54分）解：沿弹簧建立X坐标：（A ）每个弹簧单元刚度矩阵如下：总体刚度矩阵：（B）总体刚度矩阵方程：边界条件：F2x 450N，F3x 0， d1x0， d4x 0解得:d2x0.027m，d3x0.018m，F1x270N，F4x 180N(C)求单元2节点力解得:fL180N ， ?x180N解：沿杆单元建立 X坐标:（A）每个单元刚度矩阵如下:k1k2 AE11061N/mk32 106 11N/m总体刚度矩阵：（B）总体刚度矩阵方程:边界条件：F2x 5000N， F3x 0，d1x 0，d4x 03000N，F4x2000N解得：d2x 0.003m，d3x 0.0

7、01m，F1x（C）单元的应力解得：f2x 2000N , f?x2000N?x=2000A =2 10 410MPa杆单元受压有限元法基础试题（B）一、填空题（5X 2分）整体刚度矩阵方程中节点荷载由两部分组成，一是 ，二是常应变三角形单元的位移函数Ni+ Nj + Nm=。最小势能原理与虚位移原理等价，一个是以 的形式描述，另一个用 的形式表达。计算轴对称单元刚度矩阵有三种方法，一是采用数值积分，二是，三是基本的三维单元是。二、判断题（5X 2分）边界条件通常有两类。通常发生在位置完全固定不能转动的情况为非其次边界。（）线应变三角形单元中变形矩阵是 X或y的函数。（）杆单元的位移函数 2+

8、2=1。（）单元刚度矩阵keBTDBd 中，矩阵B为弹性矩阵，矩阵 D为变形矩阵。（）在梁单元中节点力与位移的方向规定应该是与材料力学中规定是一致的。（）三、简答题（26分）简述刚度矩阵的特性。（6分）写岀位移函数的含义。（4分）写出推导弹簧单元刚度矩阵的分析过程。（7分）试列举三种有限元商用软件，并说明各自优点。（9分）四、计算题（54分）对于下图所示的弹簧组合，单元的弹簧常数为 2000N/m,单元的弹簧常数为 2000N/m，节点3处位移S为0.01m，确定各节点位移、单元力和反力。（10 分） 如图所示的杆单元，杆单元弹性模量为 E,杆单元长为L,横截面面积为 A，试分别计算（a）、（

9、 b） 总体x-y坐标下的刚度矩阵。（10分）E，横截面面积均为（b）A，单元长度如图，根据对称性，求（a）对称桁架如图（a）所示，杆单元弹性模量均为 图（b）的整体刚度矩阵。（12分）(a)确定下图所示梁的各节点位移。梁已按节点编号离散化。梁在节点（b）2作用有垂直向下的力P=12kN。令沿梁弹性模量 E=70GPa，I=2 X 10-4m4，梁单元长 L=4m。弹簧常数k=200kN/m。（ 10分）如图所示梁，确定节点位移，以及每一单元的力和反作用力。梁弹性模量E=70GPa，I=3 X 10-4m4,梁单元长L=4m。作用在梁单元的均布荷载P为8 kN/m。（ 12分）参考答案（B ）

10、:一、填空题（5 X 2分）直接节点荷载等效节点荷载1能 功 直接积分 对单元中心点计算四面体单元二、判断题（5 X 2分）x V V x x三、简答题（26分）答：刚度矩阵的特性有：对称的；奇异的；主对角项总是正的。（6分）答：位移函数的含义：将单元中任意一点的位移近似地表示成该单元节点的函数。当第i个单元自由度为1而所有其他自由度值为0, Ni代表在整个单元域中假定的位移函数形状。（4分）答：推导弹簧单元刚度矩阵的分析过程为选择单元类型；选择位移函数；定义应变位移和应力应变关系；推导单元刚度矩阵和方程；组装单元方程得岀总体方程并引入边界条件；求解节点位移；求解单元力。（7分）答：ABAQU

11、S是一套先进的通用有限元系统，属于高端CAE软件。它长于非线性有限元分析，可以分析复杂的固体力学和结构力学系统，特别是能够驾驭非常庞大的复杂问题和模拟高度非线性问题。ABAQUS不但可以做单一零件的力学和多物理场的分析，同时还可以做系统级的分析和研究，其系统级分析的特点相对于其他分析软件来说是独一无二的。ANSYS软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件，发展了很多版本，但是它们核心的计算部分变化不大，只是模，在铁道，建筑和压力容器方面应E=210GPa，1=2 x 10-4m4，作E=10GPa，I=2 x 10-4m4,作用块越来越多。ANSYS系统擅长于多物理场和非线性问题