一、填空题（5×2分）
1.1单元刚度矩阵中，矩阵为__________，矩阵为___________。
1.2边界条件通常有两类。通常发生在位置完全固定不能转动的情况为_______边界，具体指定有限的非零值位移的情况，如支撑的下沉，称为_______边界。
1.3内部微元体上外力总虚功：
+的表达式中，第一项为____________________的虚功，第二项为____________________的虚功。
1.4弹簧单元的位移函数+=_________。
1.5数学表达式：令=_____， =_____，，则力。
二、判断题（5×2分）
2.1位移函数的假设合理与否将直接影响到有限元分析的计算精度、效率和可靠性。（ ）
2.2变形体虚功原理适用于一切结构（一维杆系、二维板、三位块体）、适用于任何力学行为的材料（线性和非线性），是变形体力学的普遍原理。 （ ）
2.3变形体虚功原理要求力系平衡，要求虚位移协调，是在“平衡、协调”前提下功的恒等关系。 （ ）
2.4常应变三角单元中变形矩阵是或的函数。 （ ）
2.5 对称单元中变形矩阵是或的函数。 （ ）
三、简答题（26分）
3.1列举有限元法的优点。（8分）
3.2写出有限单元法的分析过程。（8分）
3.3列出3种普通的有限元单元类型。（6分）
3.4简要阐述变形体虚位移原理。（4分）
四、计算题（54分）
4.1对于下图所示的弹簧组合，单元①的弹簧常数为10000N/m，单元②的弹簧常数为20000N/m，单元③的弹簧常数为10000N/m，确定各节点位移、反力以及单元②的单元力。（10分）
4.2对于如图所示的杆组装，弹性模量E为10GPa，杆单元长L均为2m，横截面面积A均为2×10-4m2，弹簧常数为2000kN/m，所受荷载如图。采用直接刚度法确定节点位移、作用力和单元②的应力。（10分）
4.3对称桁架如图（a）所示，杆单元弹性模量均为E，横截面面积均为A，单元长度如图，根据对称性，求图（b）的整体刚度矩阵。（12分）
（a） （b）
4.4如图所示的平面桁架，确定转换矩阵，并写出（10分）
4.5确定下图所示梁的各节点位移。梁已按节点编号离散化。梁在节点1固支，节点2有滚柱支撑，节点3作用有垂直向下的力P=50kN。令沿梁弹性模量E=210GPa，I=12×10-4m4，梁单元长L=3m。弹簧常数k=200kN/m。（12分）
参（A）：
一、填空题（5×2分）
1.1变形矩阵或应变矩阵 弹性矩阵或本构关系矩阵 1.2 齐次边界 非齐次边界
1.3 微元体上外力在随基点刚体平移所做虚功 外力在微元体变形虚位移上所做虚功 1.4 1 1 1.5 1 0
二、判断题（5×2分）
2.1 √ 2.2 √ 2.3 √ 2.4 × 2.5 √
三、简答题（26分）
3.1答：优点有：①很容易地模拟不规则形状结构；②可以很方便地处理一般荷载条件；③由于单元方程是单个建立的，因此可以模拟由几种不同材料构成的物件；④可以处理数量不受和各类边界条件；⑤单元尺寸大小可以变化；⑥改变模型比较容易⑦可以包括动态作用⑧可以处理大变形和非线性材料带来的非线性问题。（8分）
3.2答：有限元方法的一般步骤有：①离散和选择单元类型；②选择位移函数；③定义应变位移和应力应变关系；④推导单元刚度矩阵和方程；⑤组装单元方程得出总体方程并引入边界条件；⑥求解未知自由度；⑦求解单元应变和位移；⑧解释结果。（8分）
3.3答：弹簧单元，杆单元，梁单元，轴对称单元，常应变三角单元，线应变三角形单元，四面体单元等。（任意上述三种均可）（6分）
3.4答：变形体虚位移原理：受给定外力的变形体处于平衡状态的充分、必要条件是，对一切虚位移，外力所作总虚功恒等于变形体所接受的总虚变形功。（4分）
四、计算题（54分）
4.1解：沿弹簧建立X坐标：
（A）每个弹簧单元刚度矩阵如下：
总体刚度矩阵：
（B）总体刚度矩阵方程：
边界条件：，，，
解得：，，，
（C）求单元2节点力
解得：，
4.2解：沿杆单元建立X坐标：
（A）每个单元刚度矩阵如下：
N/m
N/m
总体刚度矩阵：
（B）总体刚度矩阵方程：
边界条件：，，，
解得：，，，
（C）单元②的应力
解得：，
= 杆单元②受压
有限元法基础试题（B）
一、填空题（5×2分）
1.1整体刚度矩阵方程中节点荷载由两部分组成，一是__________，二是___________。
1.2常应变三角形单元的位移函数++=_________。
1.3最小势能原理与虚位移原理等价，一个是以_____的形式描述，另一个用____的形式表达。
1.4计算轴对称单元刚度矩阵有三种方法，一是采用数值积分，二是__________________，三是__________________。
1.5基本的三维单元是_____________。
二、判断题（5×2分）
2.1边界条件通常有两类。通常发生在位置完全固定不能转动的情况为非其次边界。（ ）
2.2线应变三角形单元中变形矩阵是或的函数。 （ ）
2.3杆单元的位移函数+=1。 （ ）
2.4单元刚度矩阵中，矩阵为弹性矩阵，矩阵为变形矩阵。 （ ）
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，杆单元长为L，横截面面积为A，试分别计算（a）、（b）总体x-y坐标下的刚度矩阵。（10分）
（a） （b）
4.3对称桁架如图（a）所示，杆单元弹性模量均为E，横截面面积均为A，单元长度如图，根据对称性，求图（b）的整体刚度矩阵。（12分）
（a） （b）
4.4确定下图所示梁的各节点位移。梁已按节点编号离散化。梁在节点2作用有垂直向下的力P=12kN。令沿梁弹性模量E=70GPa，I=2×10-4m4，梁单元长L=4m。弹簧常数k=200kN/m。（10分）
4.5如图所示梁，确定节点位移，以及每一单元的力和反作用力。梁弹性模量E=70GPa，I=3×10-4m4，梁单元长L=4m。作用在梁单元②的均布荷载P为8 kN/m。（12分）
参（B）：
一、填空题（5×2分）
1.1直接节点荷载 等效节点荷载 1
1.3 能 功 1.4 直接积分 对单元中心点计算 1.5 四面体单元
二、判断题（5×2分）
2.1 × 2.2 √ 2.3 √ 2.4 × 2.5 ×
三、简答题（26分）
3.1答：刚度矩阵的特性有：①对称的；②奇异的；③主对角项总是正的。（6分）
3.2答：位移函数的含义：将单元中任意一点的位移近似地表示成该单元节点的函数。当第个单元自由度为1，而所有其他自由度值为0，代表在整个单元域中假定的位移函数形状。（4分）
3.3答：推导弹簧单元刚度矩阵的分析过程为①选择单元类型；②选择位移函数；③定义应变位移和应力应变关系；④推导单元刚度矩阵和方程；⑤组装单元方程得出总体方程并引入边界条件；⑥求解节点位移；⑦求解单元力。（7分）
3.4答：①ABAQUS是一套先进的通用有限元系统，属于高端CAE软件。它长于非线性有限元分析,可以分析复杂的固体力学和结构力学系统，特别是能够驾驭非常庞大的复杂问题和模拟高度非线性问题。ABAQUS不但可以做单一零件的力学和多物理场的分析，同时还可以做系统级的分析和研究，其系统级分析的特点相对于其他分析软件来说是独一无二的。②ANSYS软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件，发展了很多版本，但是它们核心的计算部分变化不大，只是模块越来越多。ANSYS系统擅长于多物理场和非线性问题的有限元分析,在铁道，建筑和压力容器方面应用较多。③LS-DYNA是一个通用显式非线性动力分析有限元程序，最初是1976年在美国劳伦斯利弗莫尔国家实验室由J.O.Hallquist主持开发完成的，主要目的是为核武器的弹头设计提供分析工具，后经多次扩充和改进，计算功能更为强大。LSDYNA长于冲击、接触等非线性动力分析。（9分）
四、计算题（54分）
4.5常应变三角单元（12分）
4.4如图所示梁，确定节点位移，以及每一单元的力和反作用力。梁弹性模量E=210GPa，I=2×10-4m4，作用在梁单元上的均布荷载P为5 kN/m。图中长度L=1m。（10分）
4.5如图所示梁，确定节点位移，以及每一单元的力和反作用力。梁弹性模量E=10GPa，I=2×10-4m4，作用在梁单元上的均布荷载P为1 kN/m。图中长度L=1m。（10分）下载本文
