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

《有限元分析——平面问题.ppt》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、平面问题的有限元法2014.09.161目录一、平面问题的定义1、平面应力问题2、平面应变问题二、平面问题有限元法1、结构离散2、三角形单元分析3、整体分析总体刚度矩阵4、非节点载荷移置5、边界条件处理求解三、简单算例2一、平面问题的定义1、平面应力问题平面应力问题满足以下两个条件。(1)几何条件结构是一很薄的等厚度薄板；(2)载荷条件作用于薄板上的载荷平行于板平面、沿厚度方向均匀分布，而在两板面上无外力作用。YXZOt结论：板面不受力，则有σZZ=+t/2=0τYZZ=+t/2=0τZXz=+t/2=0因板很薄，载荷又不沿厚度变化，应力沿板的厚度方向是连续分布的，可以认为，在整个板内各点都

2、有σZ=0τYZ=0τZX=0根据剪应力的互等性、物理方程，可得描述平面应力问题的八个独立的基本变量为图1平面应力问题3σ=[σXσYτXY]Tε=[εXεYγXY]Td=[μν]T它们仅为x、y的函数而与z无关。2、平面应变问题满足以下两个条件的弹性力学问题为平面应变问题。(1)结构是长柱体，横截面沿长度方向不变；（2）载荷平行于横截面且沿纵向方向均匀分布、两端不受力。ZYXOt图2平面应变问题结论：结构不能发生沿Z轴方向的位移，则有ω=0μ=μ(x,y)ν=ν(x,y)4根据几何方程、物理方程可得，描述平面应变问题的独立变量也是八个，且与平面应力问题的一样。只是弹性矩阵变为D=而平面应力

3、问题的弹性矩阵为D=5二、平面问题有限元法1、结构离散结构离散化过程：连续体结构有限单元的结合体代替原连续体平面问题用二维区域表示可用不同形状的单元，此处用三角形单元离散62、三角形单元分析⑴单元位移模式形函数根据位移函数选择方法，三节点三角形单元的位移函数μ=μ(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=7α4=(ɑ1ν1+ɑ2ν2+ɑ3ν3)α5=(b1v1+b2v2+b3v3)α6=(c1v1+c2v2+c

4、3v3)其中A=ɑ1=b1=-c1=(1,2,3)上式表示下标轮换，即12，23，31同时更换。8重写位移函数，并以节点位移的形式进行表达，有其中形函数矩阵为N=其中Ni=(ɑi+bix+ciy)，i=1、2、3。⑵单元的应变与应力单元应变ε=Bqe式中应变矩阵B为B=节点位移列阵qeqe=[u1v1u2v2u3v3]T9单元应力σ=Dε=DBqe⑶单元分析单元刚度矩阵根据虚位移原理，可得单元刚度方程Fe=Keqe其中单元刚度矩阵为Ke=对于三节点等厚三角形单元，B、D均为常数矩阵，则单元刚度矩阵可表示为Ke=BTDBtA3、非节点载荷移置有限元模型是一组仅在节点连接、仅靠节点传力、仅受节点

5、载荷、仅在节点处受约束的单元组合体。只有节点是可以承受载荷与约束的。⑴集中力的移置单元内任意一点作用集中力P=[PxPy]T10123PPxPyR1XR1YR2XR2YR3XR3YYXO图4集中力作用的单元根据虚位移原理，可得移置到节点后的载荷Re=NTP此处的N为载荷作用点的形函数值。虚功原理如下：单元原载荷在虚位移上做的虚功=移置后节点载荷在相应虚位移上做的虚功。⑵体力的移置单元所受的均匀分布体力为PV=[XY]T，则由虚功原理得Re=⑶面力的移置在单元的边上分布有单位面积上的面力PS=[]T，则由虚功原理得Re=114、整体分析整体刚度矩阵整体刚度矩阵组装的基本步骤：先求出各个单元的单

6、元刚度矩阵；将单元刚度矩阵中的每个子块放在整体刚度矩阵中的对应位置上，得到单元的扩大刚度矩阵；将全部单元的扩大矩阵相加得到整体刚度矩阵。不失一般性，仅考虑模型中有四个单元，如图所示，四个单元的整体节点位移列阵为其中：45123①②③④yox图5四个单元的模型对每个单元写出相应的单元刚度方程，对于①号单元，有12为了便于组装整体刚度矩阵，将上式以整体节点位移表示，即同理，对于单元②，有对于单元③，有13对于单元④，有整体节点载荷列阵为14整体刚度矩阵为5、边界条件处理求解整体刚度矩阵是奇异的，必须在整体刚度方程中引进位移边界条件，消除整体结构的刚体位移，再求解整体刚度方程获得节点位移的位移解。

7、下面介绍边界条件处理与计算结果整理。⑴对角元素置1法下面以一个只有四个方程的简单例子加以说明，方程如下：15假定系统中节点位移u2=0，则引入节点的已知位移后，方程变为然后，用这组维数不变的方程来求解所有的节点位移。显然，其解答仍为原方程的解答。⑵乘大数法大数M一般取108——1010。已知位移分量u2=c，把此方法用于上面的例子，则原方程变成16⑶计算结果整理一般情况，节点由多个单元所共有，而不同单元得到的