《最新有限元法基础试题》由会员分享，可在线阅读，更多相关《最新有限元法基础试题（8页珍藏版）》请在金锄头文库上搜索。
1、有限元法基础试题（A）、填空题（5X2分）1.1单元刚度矩阵keBTDBd Q中，矩阵B为，矩阵D为。1.2 边界条件通常有两类。通常发生在位置完全固定不能转动的情况为边界，具体指定有限的非零值位移的情况，如支撑的下沉，称为边界。1.3 内部微元体上外力总虚功：d8We+ Tx,xxy, y+ F)8 u + (t+a +)8bxxy,xy, ybyv dxdy +(5 ux,xy, yxy , y+ 5 u /J dxdy的表达式中，第一项为的虚功，第二项为的虚功。14弹簧单元的位移函数Ni + N2=1.5 k数学表达式：令d =，d =，k工j，则力F = k。ijjki ij二、判断题
2、（5X2分）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简要阐述变形体虚位移
3、原理。（4分）四、计算题（54 分）4.1对于下图所示的弹簧组合，单元的弹簧常数为10000N/m，单元的弹簧常数为20000N/m,单元的弹簧常数为10000N/m,确定各节点位移、反力以及单元的单元力。（10 分）14.2对于如图所示的杆组装，弹性模量E为lOGPa，杆单元长L均为2m,横截面面积A均 为2X10-4M2，弹簧常数为2000kN/m，所受荷载如图。采用直接刚度法确定节点位移、作用 力和单元的应力。 （10分）4.3对称桁架如图（a）所示，杆单元弹性模量均为E,横截面面积均为A,单元长度如图, 根据对称性，求图（b）的整体刚度矩阵。（12分）y3x21Pb)4.4如图所示的平
4、面桁架，确定转换矩阵卜丿，并写出叩丘旳卜 （10 分）x4.5确定下图所示梁的各节点位移。梁已按节点编号离散化。梁在节点1 固支，节点2有滚 柱支撑，节点3作用有垂直向下的力P=50kN。令沿梁弹性模量E=210GPa, I=12X10-4m4, 梁单元长L=3m。弹簧常数k=200kN/m。（12分）k123 1 F2/参考答案（A）：一、填空题（5X2分）1.1变形矩阵或应变矩阵 弹性矩阵或本构关系矩阵1.2齐次边界 非齐次边界1.3微元体上外力在随基点刚体平移所做虚功外力在微元体变形虚位移上所做虚功1.4丄丄1.5 丄 0二、判断题（5X2分）2.1 V 2.2 V 2.3 V 2.4
5、X 2.5 V三、简答题（26 分）31答：优点有：很容易地模拟不规则形状结构；可以很方便地处理一般荷载条件； 由于单元方程是单个建立的，因此可以模拟由几种不同材料构成的物件；可以处理数量不 受限制和各类边界条件；单元尺寸大小可以变化；改变模型比较容易可以包括动态作 用可以处理大变形和非线性材料带来的非线性问题。（8分）3.2答：有限元方法的一般步骤有：离散和选择单元类型；选择位移函数；定义应变 位移和应力应变关系；推导单元刚度矩阵和方程；组装单元方程得出总体方程并引入边 界条件；求解未知自由度；求解单元应变和位移；解释结果。（8分）3.3 答：弹簧单元，杆单元，梁单元，轴对称单元，常应变三角
6、单元，线应变三角形单元， 四面体单元等。 （任意上述三种均可） （6分）3.4 答：变形体虚位移原理：受给定外力的变形体处于平衡状态的充分、必要条件是，对一 切虚位移，外力所作总虚功恒等于变形体所接受的总虚变形功。（4分）四、计算题（54分）4.1解：沿弹簧建立X坐标：（A）每个弹簧单元刚度矩阵如下：k（1）= k（3）=1000010000-k（2）=2000020000-_1000010000 _2000020000 _总体刚度矩阵：K 二 k（1）+ k（刃 + k（3）10000100001000030000200002000030000100001000010000边界条件：F 二
7、450N ,2xF3x1x=0 , d4x解得： d2x=0.027 m,3x二0.018m , F =270N , F 二 180N1x4 x（B）总体刚度矩阵方程:r f 100001000000-d、1x1xF1000030000200000dV2 x=V2 x SF0200003000010000d3x3xF001000010000d4x114x（C）求单元2节点力20000200002000020000解得：f = 180N , f =180N2 x3x4.2解：沿杆单元建立X坐标:（A）每个单元刚度矩阵如下:k（3）= 2 x 106j-1=1x106-1-1111AEk（1）二
8、k（2）二丝LN/mN/m111总体刚度矩阵：K = k（i）+ k（2）+ k（3）1 10 01 2 1 0 K = 1x10601320 0 2 2r f -1100 -d、1x1xF1210dV2 x= 1x 106V2 x F0132d3x3xF0022d4x114x（B）总体刚度矩阵方程:边界条件：F = 5000N， F = 0， d = 0， d = 02x3x1x4 x解得：d = 0.003m2xd = 0.001m3xF =3000N，F =2000N1x4 x（C）单元的应力=1x1061111解得：f = 2000N2x= 2000N3xb（2）= -3x-xA200
9、02 x10-4= 10 MPa杆单元受压有限元法基础试题（B）一、填空题（5X2分）1.1 整体刚度矩阵方程中节点荷载由两部分组成，一是，二是。1.2常应变三角形单元的位移函数N + N . + N =。i j m1.3 最小势能原理与虚位移原理等价，一个是以的形式描述，另一个用的形式表达1.4 计算轴对称单元刚度矩阵有三种方法，一是采用数值积分，二是，三是。1.5 基本的三维单元是。二、判断题（5X2分）2.1 边界条件通常有两类。通常发生在位置完全固定不能转动的情况为非其次边界。（ ）2.2线应变三角形单元中变形矩阵是x或y的函数。（）2.3杆单元的位移函数N + N =1。（）122.
10、4单元刚度矩阵ke =f BTDBdQ中，矩阵B为弹性矩阵，矩阵D为变形矩阵。（）2.5 在梁单元中节点力与位移的方向规定应该是与材料力学中规定是一致的。（）三、简答题（26 分）3.1 简述刚度矩阵的特性。（6 分）3.2写出位移函数的含义。（4分）3.3 写出推导弹簧单元刚度矩阵的分析过程。（7 分）3.4试列举三种有限元商用软件，并说明各自优点。（9分）四、计算题（ 54 分）4.1对于下图所示的弹簧组合,单元的弹簧常数为2000N/m,单元的弹簧常数为2000N/m, 节点3处位移为0.01m，确定各节点位移、单元力和反力。（10分）4.2如图所示的杆单元，杆单元弹性模量为E,杆单元长
11、为L,横截面面积为A,试分别计 算（a）、（b）总体x-y坐标下的刚度矩阵。（10分）12a）x4.3对称桁架如图（a）所示，杆单元弹性模量均为E,横截面面积均为A,单元长度如图, 根据对称性，求图（b）的整体刚度矩阵。（12分）a)4.4确定下图所示梁的各节点位移。梁已按节点编号离散化。梁在节点2作用有垂直向下的 力P=12kN。令沿梁弹性模量E=70GPa, I=2X10-4m4,梁单元长L=4m。弹簧常数k=200kN/m。（10分）P4.5如图所示梁，确定节点位移，以及每一单元的力和反作用力。梁弹性模量E=70GPa, I=3X10-4m4，梁单元长L=4m。作用在梁单元的均布荷载P为
12、8 kN/m。（12分）1 2 |P31参考答案（B）：一、填空题（5X2分）1.1直接节点荷载 等效节点荷载1.211.3 育艮 功 1.4 直接积分 对单元中心点计算 1.5 四面体单元二、判断题（5X2分）2.1 X 2.2 V 2.3 V 2.4 X 2.5 X三、简答题（26分）3.1答：刚度矩阵的特性有：对称的；奇异的；主对角项总是正的。（6分）3.2答：位移函数的含义：将单元中任意一点的位移近似地表示成该单元节点的函数。当第i个单元自由度为1 ，而所有其他自由度值为0， N 代表在整个单元域中假定的位移函数形状。i（4 分）3.3答：推导弹簧单元刚度矩阵的分析过程为选择单元类型；
13、选择位移函数；定义应 变位移和应力应变关系；推导单元刚度矩阵和方程；组装单元方程得出总体方程并引入 边界条件； 求解节点位移； 求解单元力。 （7分）3.4答：ABAQUS是一套先进的通用有限元系统，属于高端CAE软件。它长于非线性有 限元分析,可以分析复杂的固体力学和结构力学系统，特别是能够驾驭非常庞大的复杂问题 和模拟高度非线性问题。ABAQUS不但可以做单一零件的力学和多物理场的分析，同时还 可以做系统级的分析和研究，其系统级分析的特点相对于其他分析软件来说是独一无二的。 ANSYS软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件， 发展了很多版本，但是它们核心的计算部分变化不大，只是模块越来越多ANSYS系统擅 长于多物理场和非线性问题的有限元分析 ,在铁道
