“有限元法基础及应用”补充讲义(一)-ok-m.docVIP

“有限元法基础及应用”补充讲义 一、引子——弹簧单元与弹簧系统 目标：掌握离散结构直接刚度法分析的原理和形式。了解有限元位移法列式的形式和基本概念。 1、弹簧单元分析 考虑弹簧变形平衡时的条件和弹簧物理特性，得到下列方程： 写成矩阵形式： 写成矩阵符号形式： 式（1-2）、（1-3）为弹簧单元的刚度方程，反映了单元特性：节点力与节点位移之间的关系。式中： （注意：单元节点力是节点对单元的作用力。） 弹簧单元的刚度矩阵为： 2、弹簧系统整体分析原理 以右图的一个弹簧系统为例，研究如何由单元特性集成系统特性并建立对系统进行求解的控制方程。 由前面得到的弹簧单元的刚度方程公式（1-2），分别写出2个弹簧单元的特性方程如下： 单元1 单元2 （注：右端节点力分量的下标1，2为单元节点的局部编号，上标是单元号） 下面按两个方法完成系统特性的装配和控制方程的建立。并在特定条件下求解。 1）由节点平衡方程导出： 系统处于平衡时，考虑各节点（节点1，2，3）的平衡条件： 节点受到的外载荷与节点受到与其连接的所有单元对其作用力（单元节点力的反作用力）之和等于零。因此有下列（节点）平衡方程（组）： 把单元特性（1-4），（1-5）代入（1-6）得到： 写成矩阵形式： 矩阵符号形式： 式（1-8），（1-9）就是系统平衡方程，该方程建立了离散系统的外载荷与节点位移之间的关系，是求解节点位移的控制方程。 —— 弹簧系统结构总刚度矩阵 —— 系统节点位移列阵 —— 系统节点载荷列阵 2)由单元刚度方程叠加导出 将单元1，2的刚度方程（1-4），（1-5）进行增广（扩大到系统规模）： 上述两个矩阵方程叠加，得： 上式中代入节点力平衡关系（1-6），就得到与（1-8）相同的节点平衡方程。 上述两种方法都必须考虑：①单元特性集成；②离散结构的节点上外载荷（系统外力）与节点力（系统内力）的平衡。因此方程（1-8）的本质是节点上力平衡关系，左边是由节点位移表示的（总）节点力，右边是节点所受外载荷。 3）给定载荷和约束条件下的求解 设边界条件为： 则节点平衡方程（1-8）变化为： 该方程组展开后分为2个部分： 第2，3个方程变化为： 第1个方程变化为： 先后解方程（1-15）、（1-16）得到： 至此解出了系统的未知位移和未知反力，并可以进一步求弹簧力。 3、例题 图1-4所示一个3个弹簧的系统。 求： （a） 系统总刚度矩阵 （b） 节点2，3的位移 （c） 节点1、4的反力 （d） 弹簧2中的力 解：（a）:写出各单元刚度矩阵： 应用叠加法直接得到系统总刚度矩阵： 或： 该矩阵具有如下特点：对称、奇异、稀疏、非零元素沿主对角线呈带状分布。（b）: 参考前面做法（（1-8）式）及已求出的总刚度矩阵，写出系统节点平衡方程： 考虑到位移边界条件： 则平衡方程组（1-19）第2，3方程化为： 求解上式得： （c）： 由（1-19）的方程1，4得： （d）： 弹簧2内力为： 4、练习题 二、杆单元目标：通过杆单元特性方程的建立，初步掌握有限元法单元分析的过程和原理。了解杆系结构分析的原理。 研究2节点等截面杆单元。 L— 杆长 A— 截面积 E— 弹性模量 单元的力学量和基本关系如下： 应变—位移关系： 应力—应变关系： 下面研究杆单元的单元特性。 1）直接法导出杆单元特性 应用材料力学基本知识对单元进行力学分析。 由于轴向变形模式下，杆单元的行为与弹簧单元相同，因此可比照弹簧单元的刚度方程（1-2），考虑到（2-7），直接写出杆单元的刚度方程： 写成符号形式： 杆单元刚度矩阵为： 2）公式法导出杆单元特性 步骤如下： （1）在单元上定义近似位移函数（位移模式） 把一个单元上的位移分布假设为简单多项式函数。有限元法中用插值法通过节点位移分量作为待定参数来构造单元位移函数。 对图2-1的杆单元，方便起见引入局部坐标 由于该杆单元只有2个未知位移分量，因此单元上假设的简单位移函数采用一次多项式。故对单元的节点位移进行线性插值。 容易定义出节点的插值函数如下： 因此单元上近似位移函数的插值形式为： 该位移函数也称为单元的位移模式，这里是线性位移模式。 式（2-11）中的插值函数又称为形状函数，简称形函数。 式（2-12）写成矩阵形式为： 上式中称为单元的形函数矩阵。 式（2-13）是有限元法中最重要的关系式之一，通过该式把单元上的近似位移分布函数用节点位移来表示，为进行单元层次上的分析打下了基础。 （2）单元应变和单元应力 由杆一维变形的应变——位移方