内容发布更新时间 : 2024/12/24 9:22:20星期一 下面是文章的全部内容请认真阅读。
Ti?Ti (在S?内) (i?1,2,3) (1.3.4)
其中 Ti??ijnj,nj是外界法线n的三个方向余弦。
5. 位移边界条件
ui?ui (在Su上) (i?1,2,3) (1.3.5)
6. 应变能和余能
单位体积应变能
U(?mn)?1Dijkl?ij?kl (1.3.6) 2单位体积余能
1V(?mn)?Cijkl?ij?kl (1.3.7)
2 1.3.2 平衡方程和几何方程的等效积分“弱”形式——虚功原理
虚功原理是虚位移原理和虚应力原理的总称。
作为弹性力学微分方程的等效积分形式,虚位移原理与虚应力原理分别是平衡方程与力的边界条件和几何方程与位移边界条件的等效积分形式。在导出它们的过程中都未涉及到物理方程,所以它们不仅可以用于线弹性问题,而且可以用于非线性弹性以及弹塑性等非线性问题。
将物理方程引入虚位移原理和虚应力原理可以分别导出最小位能原理和最小余能原理。它们本质上和等效积分的伽辽金“弱”形式相一致。这是建立弹性力学有限元方程的理论基础。弹性力学最小位能原理和最小余能原理都属于自然变分原理。
2 弹性力学问题有限元方法的一般原理和表达格式
通过弹性力学变分原理建立弹性力学问题有限元方法表达格式的基本步骤。最小位能原理的未知场变量是位移,以节点位移为基本未知量,并以最小位能原理为基础建立有限单元为位移元。它是有限元方法中应用最为普遍的单元。
对于一个力学或物理问题,在建立其数学模型以后,用有限元方法对它进行分析的首要步骤是选择单元形式。平面问题3节点三角形单元是有限元方法最早采用,而且至今仍经常采用的单元形式。以它作为典型,讨论如何应用广义坐标建立单元位移模式与位移插值函数,以及如何根据最小位能原理建立有限元求解方程的原理、方法与步骤,并进而引出弹性力学问题有限元方法的一般表达式。 2.1 弹性力学平面问题的有限元格式 2.1.1 单元位移模式及插值函数的构造
y
x v
um
0 vi
j
ui
vm
u
j
图2.1 3节点三角形单元
1. 单元的位移模式和广义坐标
在有限元方法中单元的位移模式或称位移函数一般采用多项式作为近似函数,因为多项式运算简便,并且随着项数的增多,可以逼近任何一段光滑的函数曲线。多项式的选取应有低次到高次。
3节点三角形单元位移模式选取一次多项式
u = ?1 + ?2x + ?3y
v = ?4 + ?5x + ?6y (2.1.1) 其中?1~?6是待定系数,称之为广义坐标。6个广义坐标可由单元的6个节点位移来表示。在(2.1.1)的1式中带入节点i的坐标(xi,yi)可得到节点i在x方向的位移ui,同理可得uj和um。它们表示为
ui??1??2xi??3yi
uj??1??2xj??3yj
um??1??2xm??3ym (2.1.2)
2. 位移插值函数
将求得的广义坐标?1~?6代入(2.1.1),可将位移函数表示成节点位移的函数,即
u?Niui?Njuj?Nmum
v?Nivi?Njvj?Nmvm (2.1.3) 其中
Ni?1(ai?bix?ciy) (i,j,m) (2.1.4) 2ANi,Nj,Nm称为单元的插值函数或形函数,对于当前情况,它是坐标x、y的一次函数,其中的bi,ci,...,cm是常数,取决于单元的3个节点坐标。 2.2.2 利用最小位能原理建立有限元方程
对于离散模型,系统总位能的离散公式
?p????U???dV???e?dS (2.2.1)
eVeeS?将结构总位能的各项矩阵表达成各个单元总位能的各对应项矩阵之和,隐含着要求单元各项矩阵的阶数(即单元的节点自由度数)和结构各项矩阵的阶数(即结构的节点自由度数)相同。为此需要引入单元节点自由度和结构节点自由度的转换矩阵G,从而将单元节点位移列阵ae用结构结点位移列阵a表示,即
ae?Ga (2.2.2)
则离散形式的总位能可表示为
?1?p?aT?GT??BTDBdVGa??BTD?0dVVe?2Vee
~??BT?0dV??NTfdV??eNTTdS??VeVeS??1?aTKa?aTP (2.2.3) 2由于离散形式的总位能?p的未知变量是结构的结点位移a,根据变分原理,泛函?p取驻值的条件是它的一次变分为零,? ? p=0,这样就得到有限元的求解方程
Ka?P (2.2.4) 其中
K??GTKeG P??GTPe (2.2.5)
eE它们都是有单元敢赌K和P分别称之为结构整体刚度矩阵和结构结点载荷列阵。矩阵Ke和单元等效结点载荷列阵Pe集合而成。
需要注意,将单元刚度矩阵和单元等效结点载荷列阵集成为结构刚度矩阵和
结构等效载荷列阵时,实际执行的并不是如(2.2.5)式所示需通过转换矩阵G的运算,而是将单元矩阵或列阵的元素直接“对号入座”,叠加到结构矩阵或列阵而成。
以上表述的是基于弹性力学最小位能原理形成的有限元求解方程的一般原理。
2.2.3 引入位移边界条件
最小位能变分原理是具有附加条件的变分原理,它要求场函数u满足几何方程和位移边界条件。现在离散模型的近似场函数在单元内部满足几何方程,因此由离散模型近似的连续体内几何方程也是满足的。但是在选择场函数的试探函数(多项式)时,却没有提出在边界上满足位移边界条件的要求,因此必须将这个条件引入有限元方程,使之得到满足。
可以引入边界条件的方法有直接代入法、对角元素改1法和对角元素乘大数法。直接代入法要重新组合方程,组成的新方程阶数降低了,但结点位移的顺序性已被破坏,这给编制程序带来了一些麻烦;对角元素改1法引入强制边界条件比较简单,不改变原来方程的阶数和结点未知量的顺序编号。但这种方法只能用于给定零位移;对角元素乘大数法使用简单,对任何给定位移(零值或非零值)都适用。采用这种方法引入强制边界条件时方程阶数不变,结点位移顺序不变,编制程序十分方便,因此在有限元法中经常采用。 2.3 广义坐标有限元法一般格式
2.3.1 广义坐标有限元位移模式的选择和插值函数的构造
1. 选择广义坐标有限元位移模式的一般原则
(1)广义坐标的个数应与单元结点自由度数相等,否则待定广义坐标?无法以单元结点位移来表示。例如,3结点三角形单元有6个自由度,因此其广义坐标个数只能是6,每个方向3个。
(2)多项式中常数项和坐标的一次项必须完备,目的是确保所选位移模式能反映单元的刚体位移和常应变特性。