有限元理论与方法-第7讲.doc

有限元理论与方法-第7讲

PAGE

PAGE7

讲授内容

备注

第7讲(第7周)

2.应变矩阵

确定了单元位移后，可以很方便地利用几何方程和物理方程求得单元的应变和应力。作为平面问题，单元内具有3个应变分量εx、εy、γxy（各符号的意义见附录1），用矩阵表示为

将（2.1.4）式代入上式中，得到

或

(2-1

式中B称为应变矩阵，写为分块形式，即

B=[BiBjBm](2-

而其子阵为

(2-1

3节点三角形单元的B是常量阵，所以称为常应变单元。在应变梯度较大(也即应力梯度较大)的部位，单元划分应适当密集，否则将不能反映应变的真实变化而导致较大的误差。

上述应变中包括与应力有关的应变和与应力无关的应变两部分，无关的应变ε0又称为初应变

ε0由温度变化、收缩、晶体生长等因素引起，对工程结构一般只考虑温度应变，无论线性和非线性温度，计算时可近似地采用平均温度

式中，Ti、Tj、Tm分别为节点i、j、m的温度，Tref为参考温度。

对于平面应力问题，温度引起的初始应变为

其中，α为线膨胀系数。

由于温度变化在各向同性介质中不引起剪切变形，所以γxy0=0。以后所述问题，除非特别说明，都指各向同性介质。

对平面应力问题，温度引起的初始应变为

5.等效节点载荷

有限单元法分析只采用节点载荷，作用于单元上的非节点载荷都必须移置为等效节点载荷。可依照静力等效原则，即原载荷与等效节点载荷在虚位移上所作的虚功相等，求等效节点载荷。

(1)集中力的移置。设单元ijm内坐标为(x，y)的任意一点M受有集中载荷f=[fxfy]T，移置为等效节点载荷Pe=[XiYiXjYjXmYm]T。假想单元发生了虚位移，其中，M点虚位移为u*=NT(δ*)e，其中(δ*)e为单元节点虚位移。按照静力等效原则有

则

(2-1-17)

(2)体力的移置。设单元承受有分布体力，单位体积的体力记为q＝[qxqy]T，其等效节点荷载为

(2-1-18)

(3)面力的移置。设在单元的某一个边界上作用有分布的面力，单位面积上的面力为p=[pxpy]T，在此边界上取微面积tds，对整个边界面积分，得到

(2-1-19)

例2-1求单元在以下受力情况下的等效节点荷载：y方向的重力为G、图2-2所示ij边受x方向均布力p、图2-3所示jm边受x方向线性分布力。

图2-2均布力图2-3线性分布力

求解说明

利用上述公式求等效节点载荷，当原载荷是分布体力或面力时，进行积分运算是比较繁琐的。但在线性位移模式下，可以按照静力学中力的分解原理直接求出等效节点载荷，上述三种情况等效节点荷载分别为

6.整体分析

结构的整体分析就是将离散后的所有单元通过节点连接成原结构物进行分析，分析过程是将所有单元的单元刚度方程组集成总体刚度方程，引进边界条件后求解整体节点位移向量。

总体刚度方程实际上就是所有节点的平衡方程，由单元刚度方程组集总体刚度方程应满足以下两个原则：

(1)各单元在公共节点上协调地彼此连接，即在公共节点处具有相同的位移。由于基本未知量为整体节点位移向量，这一点已经得到满足。

(2)结构的各节点离散出来后应满足平衡条件，也就是说，环绕某一节点的所有单元作用于该节点的节点力之和应与该节点的节点载荷平衡。

每一节点统一使用整体节点编号（如图2-4所示），第4单元节点编号i、j、m统一依次改为8、7、5。确定各单元的大域变换矩阵，如第4单元为

其中，I为2×2阶单位矩阵。

图2-4节点与单元编号

求出各单元刚度矩阵，利用大域变换法求出结构整体刚度矩阵[K]，引入边界条件，得到结构的节点平衡方程

(2-1-19)

进而，求解节点位移、单元应力和应变。

例2-2如图2-5所示，一悬臂梁，自由端受合力为P的均布力作用，梁厚t=1，μ=1/3，求节点位移。

图2-5结构与离散

求解说明

结构为平面应力问题，划分为2个三角形单元①、②，有四个节点1、2、3、4，坐标分别为（0，0）、（2，0）、（2，1）、（0，1）。单元①、②节点顺序分别取3、1、2和1、3、4，刚度矩阵完全一样。

对单元①：bi=0，bj=－1，bm=1，ci=2，cj=0，cm=－2。

对单元②：bi=0，bj=1，bm=－