有限元分析——平面问题.pptVIP

1平面问题的有限元法主讲人：杨红林

2目录一、平面问题的定义1、平面应力问题2、平面应变问题二、平面问题有限元法1、结构离散2、三角形单元分析3、整体分析总体刚度矩阵4、非节点载荷移置5、边界条件处理求解三、简单算例

3一、平面问题的定义1、平面应力问题平面应力问题满足以下两个条件。(1)几何条件结构是一很薄的等厚度薄板；(2)载荷条件作用于薄板上的载荷平行于板平面、沿厚度方向均匀分布，而在两板面上无外力作用。YXZOt结论：板面不受力，那么有σZZ=+t/2=0τYZZ=+t/2=0τZXz=+t/2=0因板很薄，载荷又不沿厚度变化，应力沿板的厚度方向是连续分布的，可以认为，在整个板内各点都有σZ=0τYZ=0τZX=0根据剪应力的互等性、物理方程，可得描述平面应力问题的八个独立的根本变量为图1平面应力问题

4σ=[σXσYτXY]Tε=[εXεYγXY]Td=[μν]T它们仅为x、y的函数而与z无关。2、平面应变问题满足以下两个条件的弹性力学问题为平面应变问题。(1)结构是长柱体，横截面沿长度方向不变；〔2〕载荷平行于横截面且沿纵向方向均匀分布、两端不受力。ZYXOt图2平面应变问题结论：结构不能发生沿Z轴方向的位移，那么有ω=0μ=μ(x,y)ν=ν(x,y)

5根据几何方程、物理方程可得，描述平面应变问题的独立变量也是八个，且与平面应力问题的一样。只是弹性矩阵变为D=而平面应力问题的弹性矩阵为D=

6二、平面问题有限元法1、结构离散结构离散化过程：连续体结构有限单元的结合体代替原连续体平面问题用二维区域表示可用不同形状的单元，此处用三角形单元离散

72、三角形单元分析⑴单元位移模式形函数根据位移函数选择方法，三节点三角形单元的位移函数μ=μ(x,y)=α1+α2x+α3yν=ν(x,y)=α4+α5x+α6yOYX节点1(x1,y1)μ1ν1节点2(x2,y2)μ2ν2节点3(x3,y3)μ3ν3图3三节点三角形单元将三个节点的位移代入，整理得α1=α2=α3=

8α4=(ɑ1ν1+ɑ2ν2+ɑ3ν3)α5=(b1v1+b2v2+b3v3)α6=(c1v1+c2v2+c3v3)其中A=ɑ1=b1=-c1=(1,2,3)上式表示下标轮换，即12，23，31同时更换。

9重写位移函数，并以节点位移的形式进行表达，有其中形函数矩阵为N=其中Ni=(ɑi+bix+ciy)，i=1、2、3。⑵单元的应变与应力单元应变ε=Bqe式中应变矩阵B为B=节点位移列阵qeqe=[u1v1u2v2u3v3]T

10单元应力σ=Dε=DBqe⑶单元分析单元刚度矩阵根据虚位移原理，可得单元刚度方程Fe=Keqe其中单元刚度矩阵为Ke=对于三节点等厚三角形单元，B、D均为常数矩阵，那么单元刚度矩阵可表示为Ke=BTDBtA3、非节点载荷移置有限元模型是一组仅在节点连接、仅靠节点传力、仅受节点载荷、仅在节点处受约束的单元组合体。只有节点是可以承受载荷与约束的。⑴集中力的移置单元内任意一点作用集中力P=[PxPy]T

11123PPxPyR1XR1YR2XR2YR3XR3YYXO图4集中力作用的单元根据虚位移原理，可得移置到节点后的载荷Re=NTP此处的N为载荷作用点的形函数值。虚功原理如下：单元原载荷在虚位移上做的虚功=移置后节点载荷在相应虚位移上做的虚功。⑵体力的移置单元所受的均匀分布体力为PV=[XY]T，那么由虚功原理得Re=⑶面力的移置在单元的边上分布有单位面积上的面力PS=[]T，那么由虚功原理得Re=

124、整体分析整体刚度矩阵整体刚度矩阵组装的根本步骤：先求出各个单元的单元刚度矩阵；将