《结构分析中的有限元法》2015-有限元习题-参考答案.doc 52页

本科有限元习题参考答案 2015年3月10日作业 1、简述力学课程中介绍的各种力学模型的简化条件、基本假设和适用范围（包括有拉压杆模型、弯曲梁模型、平面应力和平面应变模型、轴对称模型、板模型、壳模型等） 力学模型 简化条件 基本假设 适用范围 拉压杆 线弹性，无塑性变形 拉压平面假设 拉伸或压缩的向力 非线性可以是由材料性质、变形状态和边界接触条件引起的，分别称为材料、几何、边界非线性。材料非线性就是材料的本构关系不是线性的 2015年3月17日作业单元类型 单元图形 结点数 结点自由度 杆单元 2 2 梁单元 2 3 平面单元 3 2 平面四边形 4 2 轴对称问题 3 2 板壳单元 4 6 四面体单元 4 3 基函数：，， 根据拉格朗日插值多项式：。 将带入： 2015年3月24日作业 1、为了保证有限单元法解答的收敛性，位移函数应满足哪些条件？完备协调元、非协调元和完备元分别是什么意思？ 为了保证有限单元法解答的收敛性，位移函数应满足： 1) 位移函数必须包括单元的刚性位移（即常量项）； 2) 位移函数必须包括常量应变（即线性项）； 3) 位移函数在单元内部必须连续（连续性条件）； 4) 位移函数应使得相邻单元间的位移协调（协调性条件） 注：上述四个条件称为有限元解收敛于真实解的充分条件；前三个条件称为必要条件。满足四个条件的位移函数构成的单元称为完备协调元；满足前三个条件的单元称为非协调元；满足前两个条件的单元称为完备元。 2、如下图所示4结点平面应力单元，结点1~结点4对应的结点坐标分别为(0,0), (0,1), (2,0),(2,1)，结点1~结点4对应的结点位移分别为(u1,v1), (u2,v2), (u3,v3), (u4,v4)，试基于拉格朗日插值基函数构造如下单元的位移函数。 解：根据拉格朗日插值基函数： 所以： 同求：。 3、说明有限元方法解误差的主要来源？ 答：影响有限元解的误差：1) 离散误差。边界上以直线代曲线导致离散化模型与实际物体的差异；2) 位移函数误差。一般情况下单元位移函数不可能与实际单元的位移场一致；3) 计算机计算误差。计算机字长的限制、相差悬殊的数值加减运算。 4、说明用有限单元法解题的主要步骤。 答：研究问题的力学建模；结构离散；单元分析；整体分析与求解；结果分析及后处理。 5、推导基于变分原理的总势能泛函极值条件。 解：有积分形式确立的标量泛函有 其中是未知函数，和是特定的算子，是求解域，是的边界。称为未知函数的泛函，随函数的变化而变化。连续介质问题的解使泛函对于微小的变化取驻值，即泛函的“变分”等于零，此为变分法。 将虚功原理用于弹性变形时，总功W要包括外力功(T)和内力功(U)两部分，即：W = T - U；内力功(-U)前面有一负号，是由于弹性体在变形过程中，内力是克服变形而产生的，所有内力的方向总是与变形的方向相反，所以内力功取负值。 根据虚功原理，总功等于零得：T- U = 0，即外力虚功T = 内力虚功U 弹性力学中的虚功原理可表达为：在外力作用下处于平衡状态的弹性体，如果发生了虚位移，那么所有的外力在虚位移上的虚功(外力功)等于整个弹性体内应力在虚应变上的虚功(内力功)。 根据虚功原理得到 其中的即为总势能泛函。由上面变分为零式表明：在所有区域内满足几何关系，在边界上满足给定位移条件的可能位移中，真实位移使系统的总势能取驻值(可证明此驻值为最小值)。此即总势能泛函的极值条件。 2015年3月31日作业 1、给出利用最小势能原理建立单元有限元静力平衡方程的一般推导过程。 解：由单元位移函数： 式中：为插值函数（或称形函数），得到单元内的应变和应力分别为 其中为弹性矩阵，它完全取决于弹性常数和。将位移、应力和应变代入势能泛函有 根据最小势能原理，势能泛函取驻值的必要条件： 式中称为单元刚度矩阵，称为体积力等效结点力，称为面力等效结点力，则上式可写成单元平衡方程： 此式即单元有限元静力平衡方程。 2、利用最小势能原理推导结构（假设包含N个单元）有限元静力平衡方程。 解：利用弹性体所有单元的总势能最小，将所有单元进行组装后即可形成整个结构的有限元静力平衡方程。为此先需要建立单元结点位移向量和总体结点位移向量之间的对应关系，即 式中：表示第个单元的结点位移向量；是总体结点位移向量；矩阵是第个单元的结点位移提取矩阵。对于杆单元可以表示为 利用提取矩阵对单元矩阵进行合同变换与累加后，可得到总体的刚度矩阵、质量矩阵和载荷向量，即 根据最小势能变分原理即可得到，此式即为总