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

“有限元法基础及应用”补充讲义（2005年3月）一、引子——弹簧单元与弹簧系统目标：掌握离散结构直接刚度法分析的原理和形式。

了解有限元位移法列式的形式和基本概念。

1、弹簧单元分析考虑弹簧变形平衡时的条件和弹簧物理特性，得到下列方程：写成矩阵形式：写成矩阵符号形式：2个节点：节点位移：节点力：单元自由度：2j i ,j i u u ,j i f f ,弹簧的物理特性：已知弹簧力——位移关系：∆=k F-F 弹簧力，拉伸为正i j u u -=∆ —弹簧伸长量弹簧单元描述：ji i j i ku ku u u k F f -=--=-=)(ji i j j ku ku u u k F f +-=-==)(⎭⎬⎫⎩⎨⎧⎥⎦⎤⎢⎣⎡--=⎭⎬⎫⎩⎨⎧j i j i u u k kk k f f （1-1）（1-2）（1-3）图 1-1图 1-2-k 弹簧刚度kdf =式（1-2）、（1-3）为弹簧单元的刚度方程，反映了单元特性：节点力与节点位移之间的关系。

式中：（注意：单元节点力是节点对单元的作用力。

） 弹簧单元的刚度矩阵为：fd k ——单元节点力列阵——单元节点位移列阵 ——弹簧单元的刚度矩阵 弹簧单元刚度方程讨论： 1）k 有何特点？ 对称、奇异、主对角元素恒正。

2）k 中元素代表什么含义？刚度系数大小等于弹簧刚度；每列元素代表一端固定、另一端产生单位位移时加在弹簧单元上的节点力。

3）上面单元刚度方程可以求解吗？为什么？ 不可以。

刚度方程仅仅表征一个典型单元的弹性特性，单元水平上无法确定单元节点位移。

只有把系统中所有单元特性集成后，在系统水平上才可能求出所有未知位移和反力。

单元水平上，若已知单元的节点位移，可由刚度方程求出所有单元节点力分量。

若节点力已知，单元节点位移不能确定，单元可作刚体运动（小位移）。

这也是单元刚度矩阵奇异性的物理解释。

⎥⎦⎤⎢⎣⎡--=k k k k k2、弹簧系统整体分析原理以右图的一个弹簧系统为例，研究如何由单元特性集成系统特性并建立对系统进行求解的控制方程。

由前面得到的弹簧单元的刚度方程公式（1-2），分别写出2个弹簧单元的特性方程如下： 单元1 单元2（注：右端节点力分量的下标1，2为单元节点的局部编号，上标是单元号）下面按两个方法完成系统特性的装配和控制方程的建立。

并在特定条件下求解。

1）由节点平衡方程导出：系统处于平衡时，考虑各节点（节点1，2，3）的平衡条件：节点受到的外载荷与节点受到与其连接的所有单元对其作用力（单元节点力的反作用力）之和等于零。

因此有下列（节点）平衡方程（组）：把单元特性（1-4），（1-5）代入（1-6）得到：写成矩阵形式：22321122111f F f f F f F =+==（1-4）（1-5） （1-6）322233222111221111)(u k u k F u k u k k u k F u k u k F +-=-++-=-=（1-7）图 1-3矩阵符号形式：式（1-8），（1-9）就是系统平衡方程，该方程建立了离散系统的外载荷与节点位移之间的关系，是求解节点位移的控制方程。

K —— 弹簧系统结构总刚度矩阵D —— 系统节点位移列阵F —— 系统节点载荷列阵2)由单元刚度方程叠加导出将单元1，2的刚度方程（1-4），（1-5）进行增广（扩大到系统规模）：上述两个矩阵方程叠加，得：上式中代入节点力平衡关系（1-6），就得到与（1-8）相同的节点平衡方程。

上述两种方法都必须考虑：①单元特性集成；②离散结构的节点上外载荷（系统外力）与节点力（系统内力）的平衡。

因此方程（1-8）的本质是节点上力平衡关系，左边是由节点位移表示的（总）节点力，右边是节点所受外载荷。

（1-8）FK D （1-9） 讨论： （1） K 有那些特点和性质？ （2）上述方程能求解吗？（1-10）（1-11）（1-12）3）给定载荷和约束条件下的求解设边界条件为：则节点平衡方程（1-8）变化为：该方程组展开后分为2个部分： 第2，3个方程变化为：第1个方程变化为：先后解方程（1-15）、（1-16）得到：至此解出了系统的未知位移和未知反力，并可以进一步求弹簧力。

3、例题图1-4所示一个3个弹簧的系统。

求： （a ） 系统总刚度矩阵 （b ） 节点2，3的位移 （c ） 节点1、4的反力 （d ） 弹簧2中的力（1-13）（1-14）（1-15）（1-16）图 1-4,500,/100,/200,/10041321======u u N P mm N k mm N k mm N k （1-17）（1-18）PF Fu===3210解： （a ）:写出各单元刚度矩阵：应用叠加法直接得到系统总刚度矩阵： 或：该矩阵具有如下特点：对称、奇异、稀疏、非零元素沿主对角线呈带状分布。

（b ）:参考前面做法（（1-8）式）及已求出的总刚度矩阵，写出系统节点平衡方程：考虑到位移边界条件：041==u u（1-19）则平衡方程组（1-19）第2，3方程化为：求解上式得： （c ）：由（1-19）的方程1，4得： （d ）：弹簧2内力为：4、练习题())(20023200)(232222N u u k k F =-⨯=-=∆=（拉力）图 1-5对图示弹簧系统，试用叠加法求其总刚度矩阵。

并根据节点平衡方程的含义，尝试由各单元刚度矩阵的元素直接写出总刚度矩阵。

二、杆单元目标：通过杆单元特性方程的建立，初步掌握有限元法单元分析的过程和原理。

了解杆系结构分析的原理。

1、等截面杆单元及其刚度矩阵研究2节点等截面杆单元。

L — 杆长 A — 截面积 E — 弹性模量单元的力学量和基本关系如下：应变—应力—下面研究杆单元的单元特性。

1）直接法导出杆单元特性应用材料力学基本知识对单元进行力学分析。

)()()(x x x u u σσεε===——杆单元位移 ——杆单元应变 ——杆单元应力（2-1）（2-2）单元节点位移： 单元节点力： ij u u -=∆杆单元伸长量： （2-3）⎭⎬⎫⎩⎨⎧=j i f f f ⎭⎬⎫⎩⎨⎧=j i u u d图 2-1由于轴向变形模式下，杆单元的行为与弹簧单元相同，因此可比照弹簧单元的刚度方程（1-2），考虑到（2-7），直接写出杆单元的刚度方程：写成符号形式：杆单元刚度矩阵为：2）公式法导出杆单元特性 步骤如下：（1）在单元上定义近似位移函数（位移模式）把一个单元上的位移分布假设为简单多项式函数。

有限元法中用插值法通过节点位移分量作为待定参数来构造单元位移函数。

对图2-1的杆单元，方便起见引入局部坐标由于该杆单元只有2个未知位移分量，因此单元上假设的简单位移函数采用一次多项式。

故对单元的节点位移进行线性插值。

容易定义出节点的插值函数如下：杆应变：杆应力：杆内力：杆的轴向刚度：（2-4） （2-5）（2-6）（2-7）（2-8）kdf=（2-9）（2-10）因此单元上近似位移函数的插值形式为：j j i i u N u N u x u )()()()(ξξξ+==该位移函数也称为单元的位移模式，这里是线性位移模式。

式（2-11）中的插值函数又称为形状函数，简称形函数。

式（2-12）写成矩阵形式为：上式中N 称为单元的形函数矩阵。

式（2-13）是有限元法中最重要的关系式之一，通过该式把单元上的近似位移分布函数用节点位移来表示，为进行单元层次上的分析打下了基础。

（2）单元应变和单元应力由杆一维变形的应变——位移方程（2-1）和单元的位移函数（2-13）求出单元的应变分布和节点位移的关系：式中：B 称为单元的位移——应变转换矩阵，简称应变矩阵。

由一维杆的应力——应变关系（2-2），得单元应力和单元节点位移的关系：（3）用弹性体的虚位移原理导出杆单元刚度方程 变形体的虚位移：假想在弹性体上发生的，满足位移许可条件（内部连续，边界协调）的微小、任意位移场。

可以理解为某个位移场的微小扰动（变分）。

（2-11）（2-12）[]Nd =⎭⎬⎫⎩⎨⎧=j i j iu u N N u （2-13）（2-14）（2-15）BdE E ==εσ（2-16）虚位移的特征：1）假想的，与真实位移无关；2）几何上是许可的（可能的）：连续、协调； 3）微小、任意大小。

虚位移原理：弹性体受力平衡时，若发生虚位移，则外力虚功等于弹性体的虚应变能（应力在虚应变上做的虚功）。

下面把虚位移原理应用在所研究的杆单元上。

定义杆单元的虚位移：节点虚位移→单元虚位移→单元虚应变 对杆单元应用虚功原理，那么上述节点力（外力）虚功等于虚应变能，因此有下列关系：考虑到d δ的任意性，从上式可以得到：上式就是杆单元的刚度方程，杆单元的刚度矩阵为：⎭⎬⎫⎩⎨⎧=j i u u δδδd 节点虚位移： 单元虚位移：dN δδ=u fd TδdB B d Bd B d T⎪⎪⎭⎫ ⎝⎛==⎰⎰⎰dV E dV E dV VVVT T T T δδσδε单元虚应变：那么，节点力虚功为： 单元虚应变能：dB B d f d T ⎪⎪⎭⎫⎝⎛=⎰dV E V TTδδ（2-17）kd d B B f T =⎪⎪⎭⎫⎝⎛=⎰dV E V （2-18）⎰=VdVE B B k T （2-19）这就是单元刚度矩阵的通式，其导出原理和计算方法可推广到其他类型的单元。

具体计算如下：显然，与前面直接法得到的单元刚度矩阵（2-8）式相同。

3）杆单元讨论a. 只有拉伸、压缩变形的杆单元在局部坐标系下是一维问题，2节点单元只有2个节点位移分量——单元有2个自由度，单元刚度方程、刚度矩阵为2阶。

b. 单元刚度矩阵元素的物理意义： 设单元刚度方程为：令：上式代入（2-21）得到：上式表明，单元刚度矩阵第一列元素就是当单元节点位移满足式（2-22）时的单元节点力分量。

如果能设法求出此时的节点力，就得到第一列的刚度系数。

因此，一般地，单元刚度矩阵的第i （这里i=1,2）列元素等于当维持单元的第i 个自由度位移为１，其它自由度位移为0时，施加在单元上的所有节点力分量。

可以用此方法直接导出杆单元的刚度矩阵元素，试练习。

c. 显然，杆单元的刚度矩阵对称、奇异、主对角元恒正。

（2-20）⎭⎬⎫⎩⎨⎧⎥⎦⎤⎢⎣⎡=⎭⎬⎫⎩⎨⎧j i j i u u k k k k f f 22211211⎭⎬⎫⎩⎨⎧=⎭⎬⎫⎩⎨⎧01j i u u （2-21）⎭⎬⎫⎩⎨⎧=⎭⎬⎫⎩⎨⎧2111k k f f j i （2-22）（2-23）4）例题例1求图示杆中的应力。

考虑图2-2中的约束031==u u 和载荷情况后，方程（2-24）变化为：则上式第2个方程为：求解该方程后得到系统的位移解：计算单元应力：图 2-2解：结构划分为２个杆单元，单元之间在节点２铰接。