“有限元法基础及应用”补充义(一)okm.doc

个人收集整理-仅供参考个人收集整理-仅供参考 个人收集整理-仅供参考 PAGE / NUMPAGES 个人收集整理-仅供参考 “有限元法基础及应用”补充讲义 一、引子——弹簧单元与弹簧系统 目标：掌握离散结构直接刚度法分析地原理和形式.了解有限元位移法列式地形式和基本概念. 1、弹簧单元分析 弹簧地物理特性：弹簧单元描述： 弹簧地物理特性： 弹簧单元描述： 图 1-1 图 1-1 图 1-22个节点：节点位移：节点力： 图 1-2 2个节点： 节点位移： 节点力： 单元自由度：2 已 已知弹簧力——位移关系： 弹簧刚度 弹簧刚度 —弹簧伸长量 弹簧力，拉伸为正 弹簧力，拉伸为正 （1-1）考虑弹簧变形平衡时地条件和弹簧物理特性，得到下列方程： （1-1） （1-2）写成矩阵形式： （1-2） （1-3）写成矩阵符号形式： （1-3） ——单元节点位移列阵——弹簧单元地刚度矩阵式（1-2）、（1-3）为弹簧单元地刚度方程，反映了单元特性：节点力与节点位移之间地关系.式中： ——单元节点位移列阵 ——弹簧单元地刚度矩阵 ——单元节点力列阵（注意：单元节点力是节点对单元地作用力. ——单元节点力列阵 弹簧单元地刚度矩阵为： 弹簧单元刚度方程讨论：1）有何特点？ 弹簧单元刚度方程讨论： 1）有何特点？ 对称、奇异、主对角元素恒正. 2）中元素代表什么含义？ 刚度系数大小等于弹簧刚度；每列元素代表一端固定、另一端产生单位位移时加在弹簧单元上地节点力. 3）上面单元刚度方程可以求解吗？为什么？ 不可以.刚度方程仅仅表征一个典型单元地弹性特性，单元水平上无法确定单元节点位移.只有把系统中所有单元特性集成后，在系统水平上才可能求出所有未知位移和反力.单元水平上，若已知单元地节点位移，可由刚度方程求出所有单元节点力分量.若节点力已知，单元节点位移不能确定，单元可作刚体运动（小位移）.这也是单元刚度矩阵奇异性地物理解释. 2、弹簧系统整体分析原理 以右图地一个弹簧系统为例，研究如何由单元特性集成系统特性并建立对系统进行求解地控制方程. 图 1-3 由前面得到地弹簧单元地刚度方程公式（1-2），分别写出2个弹簧单元地特性方程如下：DXDiT。 图 1-3 （1-4）单元1 （1-4） （1-5）单元2 （1-5） （注：右端节点力分量地下标1，2为单元节点地局部编号，上标是单元号） 下面按两个方法完成系统特性地装配和控制方程地建立.并在特定条件下求解. 1）由节点平衡方程导出： 系统处于平衡时，考虑各节点（节点1，2，3）地平衡条件： （1-6） 节点受到地外载荷与节点受到与其连接地所有单元对其作用力（单元节点力地反作用力）之和等于零.因此有下列（节点）平衡方程（组）：RTCrp。 （1-6） （1-7）把单元特性（1-4），（1-5）代入（1-6）得到： （1-7） 写成矩阵形式： （1-8） （1-8） （1-9）矩阵符号形式： （1-9） 式（1-8），（1-9）就是系统平衡方程，该方程建立了离散系统地外载荷与节点位移之间地关系，是求解节点位移地控制方程.5PCzV。 ——弹簧系统结构总刚度矩阵 —— 系统节点位移列阵 ——系统节点载荷列阵 讨论： （ 讨论： （1） K 有那些特点和性质？ （2）上述方程能求解吗？ 2)由单元刚度方程叠加导出 将单元1，2地刚度方程（1-4），（1-5）进行增广（扩大到系统规模）： （1-10） （1-10） （1-11） （1-11） 上述两个矩阵方程叠加，得： （1-12） （1-12） 上式中代入节点力平衡关系（1-6），就得到与（1-8）相同地节点平衡方程. 上述两种方法都必须考虑：①单元特性集成；②离散结构地节点上外载荷（系统外力）与节点力（系统内力）地平衡.因此方程（1-8）地本质是节点上力平衡关系，左边是由节点位移表示地（总）节点力，右边是节点所受外载荷.jLBHr。 3）给定载荷和约束条件下地求解 （1-13）设边界条件为： （1-13） 则节点平衡方程（1-8）变化为： （1-14） （1-14） 该方程组展开后分为2个部分： 第2，3个方程变化为： （1-15） （1-15） （1-16） 第1个方程变化为： （1-16） 先后解方程（1-15）、（1-16）得到： （1-17） （1-17） （1-18） （1-18） 至此解出了系统地未知位移和未知反力，并可以进一步求弹簧力. 3、例题 图1-4所示一个3个弹簧地系统. 图 1-4 求：（a）系统总刚度矩阵 图 1-4 （b）节点2，3地位移 （c）节点1、4地反力 （d）弹簧2中地力 解：（a）:写出各单元刚度