(完整word版)深入浅出的讲清楚有限元法.docVIP

“有限元法基础及应用”补充讲义（一） 顾克秋 （2005 年 3 月） ，、引子一一弹簧单元与弹簧系统 目标：掌握离散结构直接刚度法分析的原理和形式。了解有限元位移法列式的形 式和基本概念。 1、典型弹簧单元分析 已知弹簧力一一位移关系： F二k k- 弹簧刚度 二Uj _ ui —弹簧伸长量 F - 弹簧力，拉伸为正 考虑弹簧变形平衡时的条件和弹簧物理特性，得到下列方程: fi _ _ F - - k（uj _ 5） = kui _ ku j (1-1) fj = F 二 k(uj - uj = -kui kuJ 写成矩阵形式: 写成矩阵符号形式: f 二 kd 式（1-2）、（1-3）为弹簧单元的刚度方程，反映了单元特性：节点力与节点位移 之间的关系。式中： k ――弹簧单元的刚度矩阵 d ――单元节点位移列阵 f ――单元节点力列阵 （注意：单元节点力是节点对单元的作用力） 弹簧单元刚度方程讨论： 1） k有何特点？ 对称、奇异、主对角元素恒正 2） k中元素代表什么含义？ 刚度系数大小等于弹簧刚度；每列元素代表一端固定、另一端产生单位 位移时加在弹簧单元上的节点力。 3）上面单元刚度方程可以求解吗？为什么？ 不可以。刚度方程仅仅表征一个典型单元的弹性特性，单元水平上无法 确定单元节点位移。只有把系统中所有单元特性集成后，在系统水平上才可 能求出所有未知位移和反力。单元水平上，若已知单元的节点位移，可由刚 度方程求出所有单元节点力分量。若节点力已知，单元节点位移不能确定， 单元可作刚体运动（小位移）。这也是单元刚度矩阵奇异性的物理解释。  以右图的一个弹簧系统为例，研究如何由单元特性集成系统特性并建立对系 统进行求解的控制方程。 由前面得到的弹簧单元的刚度方程公式 （1-2），分别写出2个弹簧单元的特 性方程如下:  （注：右端节点力分量的下标1, 2为单元节点的局部编号，上标是单元号） 下面按两个方法完成系统特性的装配和控制方程的建立。并在特定条件下求 1）由节点平衡方程导出： 系统处于平衡时，考虑各节点（1, 2, 3节点）的平衡条件： 节点受到的外载荷与节点受到与其连接的所有单元对其作用力 （单元节点力 的反作用力）之和等于零。因此有下列（节点）平衡方程（组）  把单元特性（1-4），（ 1-5）代入（1-6）得到: k1u2 F2 二 _k1u1 (k「k2)u2 - k2u3 F3 二-k2± k2U3 写成矩阵形式: 或矩阵符号形式: （1-9） 式（1-8 ）, （1-9）就是系统平衡方程，该方程建立了离散系统的外载荷与节 点位移之间的关系，是求解节点位移的控制方程。 弹簧系统的结构总刚度矩阵 D ――系统节点位移列阵 F —系统节点载荷列阵 讨论：（1）K有那些特点和性质? （2）上述方程能求解吗？ 2）由单元刚度方程叠加导出 将单元1, 2的刚度方程（1-4），（ 1-5）进行增广（扩大到系统规模）： 处 一 k] 0 IL /；q 一 0 * i * = J 1 fi （1-10） _ 0 0 oj v 0 「0 0 0 _ 比 0 1 0 k2 -k2 »= 卜 （1-11） 0 -k. k、 上述两个矩阵方程叠加，得: k\ 0 ■ 广 一何 A] +心 -k2 * [0 （1-12） -k2 k2 £ ■ V J j [fl \  上式中代入节点力平衡关系（1-6）,就得到与（1-8 ）相同的节点平衡方程。 上述两种方法都必须考虑1）单元特性集成；2）离散结构的节点上外载荷（系 统外力）与节点力（系统内力）的平衡。因此方程（1-8 ）的本质是节点的力平 衡关系，左边是由节点位移表示的（总）节点力，右边是节点所受外载荷。 3）给定载荷和约束条件下的求解 设边界条件为： 比=0 F F P 2 3 则节点平衡方程（1-8 ）变化为:  0 0 代十处 卜=€ _ 0 P k J (1-14) 该方程组展开后分为2个部分: 第1个方程变化为: 第2, 3个方程变化为： （1-15）   (1-16) 先后解方程（1-15 ）、（ 1-16）得到: 卜be』 卜:=-2 P从而解出了系统的未知位移和未知反力，并可以进一步求弹簧力 3、例题 图1-4所示一个3个弹簧的系统。 =100N/mm,k2 =200N/mm, k3 =100N / mm, P =500N,Ui 二 U4 =0 :（a） 系统总刚度矩阵 (b) 节点2, 3的位移 (c) 节点1、4的反力 (d) 弹簧2中的力 解： （a）: 写出各单元刚度矩阵: 100 -100 -100 100 ■ 200 -200 -200 200 _ 100 -I0CT -100 100 (N/mm) (N/nnn)