深入浅出讲清楚有限元法.docVIP

(完好word版)深入浅出的讲清楚有限元法 (完好word版)深入浅出的讲清楚有限元法 PAGE / NUMPAGES (完好word版)深入浅出的讲清楚有限元法 “有限元法基础及应用”增补讲义（一） 顾克秋 2005 年 3 月） 一、引子——弹簧单元与弹簧系统 目标：掌握失散结构直接刚度法剖析的原理和形式。 认识有限元位移法列式的形式和基本观点。 1、典型弹簧单元剖析 弹簧单元描绘： 弹簧的物理特征： 图 1-1 2个节点： i, j 图 1-2 节点位移： ui ,u j 已知弹簧力——位移关系： F k k 弹簧刚度 节点力： fi , f j u j ui —弹簧伸长量 单元自由度： 2 F 弹簧力，拉伸为正 考虑弹簧变形均衡时的条件和弹簧物理特征，获得以下方程： fi F k(u j ui ) kui ku j f j F k(u j ui ) kui （ 1-1） ku j 写成矩阵形式： fi k k ui f j k k （ 1-2） u j 写成矩阵符号形式： f kd （ 1-3） 1 式（ 1-2）、（1-3）为弹簧单元的刚度方程，反应了单元特征：节点力与节点位移 之间的关系。式中： k ——弹簧单元的刚度矩阵 d ——单元节点位移排阵 f ——单元节点力排阵 （注意：单元节点力是节点对单元的作使劲） 弹簧单元刚度方程议论： 1） k 有何特色？ 对称、奇怪、主对角元素恒正 2） k 中元素代表什么含义？ 刚度系数大小等于弹簧刚度；每列元素代表一端固定、另一端产生单位位移时加在弹簧单元上的节点力。 3）上边单元刚度方程能够求解吗？为何？ 不可以够。刚度方程只是表征一个典型单元的弹性特征，单元水平上没法确立单元节点位移。只有把系统中所有单元特征集成后，在系统水平上才可能求出所有未知位移和反力。单元水平上，若已知单元的节点位移，可由刚度方程求出所有单元节点力重量。若节点力已知，单元节点位移不可以确立，单元可作刚体运动（小位移） 。这也是单元刚度矩阵奇怪性的物理解说。 2、弹簧系统整体剖析原理 以右图的一个弹簧系统为例， 研究怎样由单元特征集成系统特征并成立对系统进行求解的控制方程。 由前面获得的弹簧单元的刚度方程公式 （1-2 ），分别写出 2 个弹簧单元的特征方程以下： 图 1-3 单元 1 （ 1-4） 单元 2 （ 1-5） （注：右端节点力重量的下标 1， 2 为单元节点的局部编号，上标是单元号） 下边按两个方法达成系统特征的装置和控制方程的成立。 并在特定条件下求 解。 1）由节点均衡方程导出： 系统处于均衡时，考虑各节点（ 1，2，3 节点）的均衡条件： 节点遇到的外载荷与节点遇到与其连结的所有单元对其作使劲 （单元节点力 的反作使劲）之和等于零。所以有以下（节点）均衡方程（组） ： F f 1 1 1 F f 1 f 2 2 2 1 （ 1-6） F f 2 3 2 把单元特征（ 1-4 ），（1-5 ）代入（ 1-6 ）获得： F1 k1u1 k1u2 F2 k1u1 （ 1-7） (k1 k2 )u2 k2u3 F3 k2u2 k2u3 写成矩阵形式： （ 1-8） 3 或矩阵符号形式： KD F （ 1-9） 式（1-8 ），（ 1-9 ）就是系统均衡方程，该方程成立了失散系统的外载荷与节点位移之间的关系，是求解节点位移的控制方程。 K D F  —— 弹簧系统的结构总刚度矩阵 —— 系统节点位移排阵 —— 系统节点载荷排阵 议论：（ 1） K 有那些特色和性质？ 2）上述方程能求解吗？ 由单元刚度方程叠加导出 将单元 1，2 的刚度方程（ 1-4 ），（ 1-5 ）进行增广（扩大到系统规模） ： （ 1-10） （ 1-11） 上述两个矩阵方程叠加，得： （ 1-12） 4 上式中代入节点力均衡关系（ 1-6 ），就获得与（ 1-8 ）同样的节点均衡方程。上述两种方法都一定考虑 1）单元特征集成； 2）失散结构的节点上外载荷（系统外力）与节点力（系统内力）的均衡。所以方程（ 1-8 ）的实质是节点的力均衡关系，左侧是由节点位移表示的（总）节点力，右侧是节点所受外载荷。 3）给定载荷和拘束条件下的求解 设界限条件为： u1 0 （ 1-13） F2 F3 P 则节点均衡方程（ 1-8 ）变化为： （ 1-14） 该方程组睁开后分为 2 个部分： 第 2，3 个方程变化为： （ 1-15） 第 1 个方程变化为： （ 1-16） 先后解方程（ 1-15 ）、（1-16 ）获得： （ 1-17） （ 1-18） 5 进而解出了系统的未知位移和未知反力，并能够进一步求弹簧力。 3、例题 图 1-4 所示一个 3 个弹簧的系统。 k1 100N