《有限元法基础试题》由会员分享，可在线阅读，更多相关《有限元法基础试题（8页珍藏版）》请在金锄头文库上搜索。
1、有限元法基础试题（A）、填空题（5X2分）1.1单元刚度矩阵keBrDBd0中，矩阵B为，矩阵D为。Q边界条件通常有两类。通常发生在位置完全固定不能转动的情况为边界，具体指定有限的非零值位移的情况，如支撑的下沉，称为边界。1.2 内部微元体上外力总虚功：d8We+Tx,xxy,y+Fbx）8xy,x+Gy,y+Fby)5vdxdyx,xy,yxy,y+5uJdxdy的表达式中，第一项为的虚功，第二项为的虚功。1.4弹簧单元的位移函数N+N=121.5k数学表达式：令d=，d=，k工j，则力F二k。ijjkiij二、判断题（5X2分）位移函数的假设合理与否将直接影响到有限元分析的计算精度、效率和
2、可靠性。（）2.2变形体虚功原理适用于一切结构（一维杆系、二维板、三位块体）、适用于任何力学行为的材料（线性和非线性），是变形体力学的普遍原理。（）变形体虚功原理要求力系平衡，要求虚位移协调，是在“平衡、协调”前提下功的恒等关系。（）2.4常应变三角单元中变形矩阵是x或y的函数。（）2.5对称单元中变形矩阵是x或y的函数。（）三、简答题（26分）列举有限元法的优点。（8分）3.2写出有限单元法的分析过程。（8分）列出3种普通的有限元单元类型。（6分）3.4简要阐述变形体虚位移原理。（4分）四、计算题（54分）4.1对于下图所示的弹簧组合，单元的弹簧常数为10000N/m，单元的弹簧常数为200
3、00N/m,单元的弹簧常数为10000N/m,确定各节点位移、反力以及单元的单元力。（10分）4.2对于如图所示的杆组装，弹性模量E为10GPa，杆单元长L均为2m,横截面面积A均为2X10-4m2，弹簧常数为2000kN/m，所受荷载如图。采用直接刚度法确定节点位移、作用力和单元的应力。（10分）A.4123LLk5kN24.3对称桁架如图（a）所示，杆单元弹性模量均为E,横截面面积均为A,单元长度如图,根据对称性，求图（b）的整体刚度矩阵。（12分）1/22”Pb)4.4如图所示的平面桁架，确定转换矩阵印，并写出叩丘比卜（10分）4.5确定下图所示梁的各节点位移。梁已按节点编号离散化。梁在
4、节点1固支，节点2有滚柱支撑，节点3作用有垂直向下的力P=50kN。令沿梁弹性模量E=210GPa,I=12X10-4m4,梁单元长L=3m。弹簧常数k=200kN/m。（12分）参考答案（A）：一、填空题（5X2分）1.1变形矩阵或应变矩阵弹性矩阵或本构关系矩阵1.2齐次边界非齐次边界1.3微元体上外力在随基点刚体平移所做虚功外力在微元体变形虚位移上所做虚功1.4丄丄1.5丄0二、判断题（5X2分）V2.2V2.3V2.4X2.5V三、简答题（26分）31答：优点有：很容易地模拟不规则形状结构；可以很方便地处理一般荷载条件；由于单元方程是单个建立的，因此可以模拟由几种不同材料构成的物件；可以
5、处理数量不受限制和各类边界条件；单元尺寸大小可以变化；改变模型比较容易可以包括动态作用可以处理大变形和非线性材料带来的非线性问题。（8分）3.2答：有限元方法的一般步骤有：离散和选择单元类型；选择位移函数；定义应变位移和应力应变关系；推导单元刚度矩阵和方程；组装单元方程得出总体方程并引入边界条件；求解未知自由度；求解单元应变和位移；解释结果。（8分）答：弹簧单元，杆单元，梁单元，轴对称单元，常应变三角单元，线应变三角形单元，四面体单元等。（任意上述三种均可）（6分）答：变形体虚位移原理：受给定外力的变形体处于平衡状态的充分、必要条件是，对一切虚位移，外力所作总虚功恒等于变形体所接受的总虚变形功
6、。（4分）四、计算题（54分）4.1解：沿弹簧建立X坐标：（A）每个弹簧单元刚度矩阵如下：k（1）=k（3）=10000100001000010000k（2）=20000-20002000020000总体刚度矩阵：K=k（1）+kC）+k（3）00100001000010000100000100003000020000020000300000010000（B）总体刚度矩阵方程:rf-100001000000-d、1x1xF1000030000200000dV2x=V2x,F0200003000010000d3x3xF4x_001000010000_d4x边界条件：F=450N,F=0,d=0,
7、d=02x3x1x4x解得：d二0.027m,2xd二0.018m,F=270N,F二180N3x1x4x（C）求单元2节点力f-20000-200001/2x=/J-2000020000d2xd3x解得：/=180N2x，f=180N3x4.2解：沿杆单元建立X坐标:（A）每个单元刚度矩阵如下:AEj-11二1x106-1-11L-11-11k（i）二k（2）二N/m1N/m-11k（3）=2x106-1总体刚度矩阵：-1K=k（i）+k（2）+k（3）K二lx10-1-1-1-2rf-1-1001d、1x1xF-12-10dV2x=1x106V2x,F0-13-2d3x3xF00-22d4
8、x114x（B）总体刚度矩阵方程:二5000N边界条件:F二0，3xF2xd1xd二04x二0.003md3x二0.001m,F=-3000N，F1x4x-2000N厂八Jrff-1-V2x=1x106、f一-11（C）单元的应力d2xd3x解得：d2x解得:f=2000N2x，f=-2000N3x八丄-2000xA2x10-4一10丹杆单元受压有限元法基础试题（B）1.1 一、填空题（5X2分）整体刚度矩阵方程中节点荷载由两部分组成，一是，二是。1.2常应变三角形单元的位移函数N+N+N=。ijm最小势能原理与虚位移原理等价，一个是以的形式描述，另一个用的形式表达计算轴对称单元刚度矩阵有三种
9、方法，一是采用数值积分，二是，三是。1.3 基本的三维单元是。二、判断题（5X2分）边界条件通常有两类。通常发生在位置完全固定不能转动的情况为非其次边界。（）2.2线应变三角形单元中变形矩阵是x或y的函数。（）2.3杆单元的位移函数N+N=1。（）122.4单元刚度矩阵keBiDBd0中，矩阵B为弹性矩阵，矩阵D为变形矩阵。（）Q在梁单元中节点力与位移的方向规定应该是与材料力学中规定是一致的。（）三、简答题（26分）简述刚度矩阵的特性。（6分）3.2写出位移函数的含义。（4分）写出推导弹簧单元刚度矩阵的分析过程。（7分）3.4试列举三种有限元商用软件，并说明各自优点。（9分）四、计算题（54分
10、）4.1对于下图所示的弹簧组合,单元的弹簧常数为2000N/m,单元的弹簧常数为2000N/m,节点3处位移5为0.01m,确定各节点位移、单元力和反力。（10分）4.2如图所示的杆单元，杆单元弹性模量为E,杆单元长为L,横截面面积为A,试分别计算（a）、（b）总体x-y坐标下的刚度矩阵。（10分）12a）4.3对称桁架如图（a）所示，杆单元弹性模量均为E,横截面面积均为A,单元长度如图,根据对称性，求图（b）的整体刚度矩阵。（12分）L3一L一y45P3X11a)b)4.4确定下图所示梁的各节点位移。梁已按节点编号离散化。梁在节点2作用有垂直向下的力P=12kN。令沿梁弹性模量E=70GPa
11、,I=2X10-4m4,梁单元长L=4m。弹簧常数k=200kN/m。（10分）1O2r.O311P4.5如图所示梁，确定节点位移，以及每一单元的力和反作用力。梁弹性模量E=70GPa,I=3X10-4m4，梁单元长L=4m。作用在梁单元的均布荷载P为8kN/m。（12分）P3参考答案（B）：一、填空题（5X2分）1.1直接节点荷载等效节点荷载1.21育艮功1.4直接积分对单元中心点计算1.5四面体单元二、判断题（5X2分）X2.2V2.3V2.4X2.5X三、简答题（26分）3.1答：刚度矩阵的特性有：对称的；奇异的；主对角项总是正的。（6分）3.2答：位移函数的含义：将单元中任意一点的位移
12、近似地表示成该单元节点的函数。当第i个单元自由度为1，而所有其他自由度值为0，N代表在整个单元域中假定的位移函数形状。i（4分）3.3答：推导弹簧单元刚度矩阵的分析过程为选择单元类型；选择位移函数；定义应变位移和应力应变关系；推导单元刚度矩阵和方程；组装单元方程得出总体方程并引入边界条件；求解节点位移；求解单元力。（7分）3.4答：ABAQUS是一套先进的通用有限元系统，属于高端CAE软件。它长于非线性有限元分析,可以分析复杂的固体力学和结构力学系统，特别是能够驾驭非常庞大的复杂问题和模拟高度非线性问题。ABAQUS不但可以做单一零件的力学和多物理场的分析，同时还可以做系统级的分析和研究，其系
13、统级分析的特点相对于其他分析软件来说是独一无二的。ANSYS软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件，发展了很多版本，但是它们核心的计算部分变化不大，只是模块越来越多ANSYS系统擅长于多物理场和非线性问题的有限元分析,在铁道，建筑和压力容器方面应用较多。LS-DYNA是一个通用显式非线性动力分析有限元程序，最初是1976年在美国劳伦斯利弗莫尔国家实验室由主持开发完成的，主要目的是为核武器的弹头设计提供分析工具，后经多次扩充和改进，计算功能更为强大。LSDYNA长于冲击、接触等非线性动力分析。（9分）四、计算题（54分）4.5常应变三角单元（12分）321y41xPP4.4如图所示梁，确定节点位移，以及每一单元的力和反作用力。梁弹性模量E=210GPa,I=2X10-4m4，作用在梁单元上的均布荷载P为5kN/m。图中长度L=1m。（10分）P5L|4L4.5如图所示梁，确定节点位移，以及每一单元的力和反作用力。梁弹性模量E=10GPa,I=2X10-4m4，作用在梁单元上的均布荷载P为1kN/m。图中长度L=1m。（10分）P
