《弹性力学与有限元》第5章平面问题有限元分析Word格式.docx
《《弹性力学与有限元》第5章平面问题有限元分析Word格式.docx》由会员分享,可在线阅读,更多相关《《弹性力学与有限元》第5章平面问题有限元分析Word格式.docx(50页珍藏版)》请在冰豆网上搜索。
在平面应力问题中,最简单、最常用的单元类型是三角形三节点单元。
如果我们选择这种类型的单元,就可以用一组网格线把求解区域图5-2实线部分划分为若干个三角形单元,如图5-3所示。
网格的交点或单元的角点就是节点。
每个单元有三个节点,把每个节点沿坐标轴的位移取为基本未知点数。
由于是平面问题,所以每个节点有两个自由度,
所有节点都
图5-2结构的简化
图5-3有限元剖分图
可以取为铰接。
对于非节点上的外载荷按静力等效的原则移置到节点上去,成为节点载荷。
这样就把平面连续体(或结构离散成为由有限个单元和节点所组成的等效集合体。
有了计算模型,就可以用处理杆系问题的基本思想来分析单元和结构特性。
3.结构离散化之后,对单元和节点分别进行编号
将结构离散化之后,把所有的单元及节点按一定顺序分别进行编号,以便进行计算。
编号时应注意:
单元号及节点号均不能有错漏或重复。
需要指出,编号的顺序不影响计算结果,原则上可以任意编排。
但是为了节省计算机内容,减少计算时间,单元中每个节点的编号与周围节点的编号应尽可能接近,即同一个单元内的节点编号差值的最大值尽量的小,这样就可以使在求解线性方程组时,使系数矩阵中非零元素的宽度(即带宽为最小。
4.定义单元
单元和节点分别编号后,还要定义单元(即定义每个单元的i,j,m,且单元一经定义后,整个计算中不允许再改变了。
最好把求解区域分成内部单元和边界单元。
对于内部单元,三个节点i,j,m的位置可以是任意的。
对于边界单元,为了使计算公式简化,格式统一起见,只准有一条单元边界处于求解区域边界上,并且规定节点j,m处于边界上。
5.1.3离散化应注意的问题
划分有限单元网格时,必须注意以下几点:
1.单元数目的确定应兼顾精度、经济性和计算机容量。
从有限元法本身来讲,单元划分得越细,节点布置得越多,越接近真实求解区域,因而计算精度越高。
但是随之而来的就是计算时间、计算费用和计算机内存的增加。
所以要考虑各因素来恰当地确定单元的数目。
其原则是,在满足工程精度要求的前提下,单元数目应尽量划分得少一些。
2.在初步分析的基础上,再离散化。
首先根据求解区域的形状、载荷分布情况和边界条件大致分析,在此基础上进行有限元分割(或剖分,划分。
3.如果求解区域具有不同的厚度或者由两种以上的材料组成时,不要把厚度不同或材料不同的区域划分在一个单元中。
4.单元的边长尽量接近,以便提高计算精度。
对于三角形单元,不应出现过大的钝角或过小的锐角。
对于四边形单元,长度和宽度也不宜相差太大。
5.任意一个单元的角点(顶点或节点必须同时也是相邻单元的角点,而不能是相邻单元边上的内点。
5.2单元位移模式的插值函数
常应变三角形单元是最早提出、最简单的一种单元,对于从杆系结构转入连续体问题,理解有限单元法的原理、方法和步骤,从而更好地掌握有限单元法是十分必要的。
这种单元虽然简单,但它也可用来拟合复杂边界体形。
虽然边界为曲线时,存在以三角形直边来代替而带来离散误差,但比用矩形单元进行离散,离散误差要小得多。
5.2.1单元结点位移和结点力
图5-4为任一典型单元。
单元局部结点编号记作1、2、3,逆时针进行标记,其对应的整体结点编号记作i、j、m。
由于平面问题局部坐标和整体坐标是一致的,因此没有坐标转换等问题,故也可只标
记整体编号以便形成定位向量。
如图所示,单元每个结点有两个位移分量,称为结点位移,记为:
iiiuvδ⎛⎞
=⎜⎟⎝⎠
将三个结点位移按结点编号排在一起称为单元结点位移,记为:
(T
eTTT
ijmδδδδ=
与结点位移相对应,每个结点上受有2个其它单元对它作用的力,称为结点力,记为
ixiiyFFF⎛⎞=⎜⎟
⎝⎠
将三个结点力按结点编号排在一起称为单元结点力,记为
(
Te
T
TTij
m
FFFF
=
单元上作用的体积力记为
bxbye
b
FFF⎛⎞=⎜⎟⎝⎠
若单元的边界是物体边界,并且该边界有表面力作用的话,该表面力记为
jx
图5-4典型单元
sxsyesFFF⎛⎞=⎜⎟⎝⎠
体积力和表面力表达式中的矩阵元素均是沿坐标方向的分布荷载集度。
5.2.2插值函数的建立
如图5-5所示是一个典型三节点三角形平面单元,三个节点为mji、、,按逆时针方向排列。
节点位移分量用一个向量表示为:
{}⎪
⎪⎪
⎪⎭⎪⎪⎪
⎪⎬⎫⎪⎪⎪⎪⎩⎪
⎪⎪⎪⎨⎧=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧=nmjji
imjie
vuvuvuδδδδ只知道节点位移,并不能直接求出单元内的应变和应力。
因此,必须假定单元内任一点的位移为坐标的某种函数。
有限元中,单元的位移模式一般采用多项式,即
y
xyxvyxyxu654321,(,(αααααα++=++=
将三个节点的坐标及位移代入上式得到六个多项式:
mmjjjiiimmmjjjiiiyxvyxvyxvyxuyxuyxu654654654321321321αααααααααααααααααα++=++=++=++=++=++=,,,,
按照线性代数中的克莱姆法则,从上面式中解出系数iα,得到单元位移函数如下:
X
Y
图5-5三节点三角形平面单元
((({}((({}mmmmjjjjiiiimmmmjjjjiiiivycxbavycxbavycxbaA
vuycxbauycxbauycxbaAu++++++++=++++++++=
21
其中
i
jmjimijjimmijimjmiimjjmimjijmmjixxcyybyxyxaxxcyybyxyxaxxcyybyxyxa−=−=−=−=−=−=−=−=−=,,,,,,m
jjiiy
x
yxyxA1
1
121=
按照解析几何的理论,A等于三角形ijm的面积,为了保证面积A为正,所以节点
mji、、的次序必须是逆时针方向。
为了简化位移函数的表达式,记
A
cxbaNAycxbaNAycxbaNmmmmjjjjiiii222++=++=++=
,
代入上面式子,得到位移函数的简易表达式如下:
jjmmiijjmmiivNvNvNvuNuNuNu++=++=,(5-2
写成矩阵形式如下:
{}[
]
{}[]{}eemj
immjjiimj
mji
NINININvuvuvuNNNNNNvuuδδ==⎪
⎪⎪⎪⎭⎪⎪
⎪⎪⎬⎫⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧⎥⎦⎤⎢⎣⎡=⎭⎬⎫⎩⎨
⎧=0
000
(5-3
[]N是坐标的函数,它们反映了单元的位移形态,称为单元的形函数或插值函数。
上式按虚线所示分块,其各子块分别记作
22222200010[];
[];
[]00001j
iijmjimNNNNINININNNN×
×
⎡⎤⎡⎤⎡⎤
⎡⎤====⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦
⎣⎦式中[I]为二阶单位矩阵。
由式(5-2可知,当1iiuv==,而其它节点位移为零时,单元内任何一点处的位移为
ii
uNvN==
就是说节点I发生位移时,函数iN就表示单元内部的位移分布形状,故称iN,jN,mN为位移的形函数(或形状函数。
形函数是定义于单元内部的坐标的连续函数。
对于三角形三节点单元,其形函数是坐标的线性函数,因此位移函数式(5-2也是坐标的线性函数。
线性位移函数的一个特点是:
单元上的任意一条直线在变形后,仍然是一条直线。
所以相邻单元的公共边线,在变形后仍是一条直线,且只要相邻两个单元在公共节点处保持位移相等,则公共边线在变形后仍然保持密合。
因此,对于三角形单元来说,选择线性位移函数是能够满足协调性要求的。
5.2.3插值函数的性质
(1在单元节点上,插值函数的值为
⎩
⎨
⎧≠==jij
iyxNjji当当01,(;
(i,j,m即有(
1iiiNxy=,(
(,,0ijjimmNxyNxy==。
(2在单元中任意一点各插值函数之和应等于1,即
ijmNNN++=(3.2.12
(3插值函数是坐标x、y的线性函数,因此单元内及边界上的位移也是线性变化的。
5.3单元应变和应力
平面问题单元内具有三个应变分量xyyxγεε、、,可用矩阵表示如下:
{}[]{}[]{}mmjjiimj
mjievuvuvuNNNNNNxyyxNxyyxuLδε⎪
⎪⎪⎪⎭⎪⎪⎪⎪⎬⎫⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧⎥⎦⎤⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡∂∂∂∂∂∂∂∂=⎥⎥⎥
⎥
⎥⎥⎥⎦⎤⎢⎢
⎢⎢⎢⎢
⎢⎣
⎡∂∂∂∂∂∂∂∂
==0
0000
[]{}e
mmjjiimmjjiimji
mjiBvuvuvuxNy
Nx
Ny
NyNyNyNxNxNxNδ=⎪⎪⎪
⎪⎭
⎪⎪⎪⎪
⎬
⎫⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂=000000
(5-3其中
[][]
⎥⎥⎥⎦
⎤⎢⎢⎢
⎣
⎡=
=mm
j
imj
ji
bcbcbccccbbbABBBB000000
21[][][](mjibccbA
NLBii
ii,,0021,⎥⎥⎥⎦
⎣⎡==(5-4
其余两式脚码用j,m代换。
由于单元面积A以及系数iicb、都是常量,所以矩阵[B]的元素都是常量,因此应变{ε}的元素也是常量,也就是说,在每个单元中,三个应变分量都是常量,这种单元称为常应变单元。
在求出单元应变之后,根据前面对平面问题的研究不难求出单元应力如下:
{}[]{}[][]{}[]{}eeSBDDδδεσ===(5-5
其中[]D称为弹性矩阵,[]S称为应力矩阵,由弹性矩阵和应变矩阵相乘而得。
对于平面应力问题
[]⎥⎥⎥⎥
⎦
⎡−−=210
0010112
µ
ED(5-6[]210000100002(1100
2i
jmi
mijmi
ijjm
mbbbE
SssscccAcbcbcbµ
⎡
⎤⎢⎥⎡⎤⎢⎥⎢⎥⎡⎤==⎢⎥⎢⎥⎣⎦−⎢⎥⎢⎥−⎣⎦
⎢⎥⎣
(5-7
[]2,(,,2(1112
2
iiiiiiibcE
SbcijmAcbµ
⎡⎤⎢⎥
⎢⎥=
⎢⎥−⎢⎥
−−⎢⎥⎣⎦(5-8
平面应变问题把1µ
−,2
1E
Eµ
−代入上式(5-6,可的得到平面应变公式中的[]D:
[]D=
⎥⎥
⎥⎥⎥⎥⎦⎤
⎢⎢⎢⎢⎢⎢
⎡−−−−−+−1(2210
00110
1121(1(1(µ
E
同理,可得到平面应变公式中的[]iS:
[]⎥⎥
⎥⎦
⎤
⎢⎢
⎢⎣⎡
−−−−−−−+−=iii
iii
ibccbcbES1(2211(2211121(1(21(µ
因为以上式子中iicbE、、、µ
等都是常量,所以矩阵[S]的元素都是常量,也就是说,单元内的应力都是常量。
同样,节点力[]e
F可由虚功方程导出:
{}
[]{}T
ee
e
mFFFFkδ⎡⎤==⎣⎦
其中:
[]
ii
ijime
jijjjmmimj
mmkkkkkkkkkk⎡⎤⎢⎥=⎢⎥⎢⎥⎣
[]⎥
⎥⎥⎦
⎢⎢⎢⎣
⎡−+−+−+
−+−=srsrs
rsrsrsrsrsrrsbbccc
bb
cbccbccbbAEtk21212
1211(42
(r,s=i,j,m。
对于平面应变问题,µ
E做与上式相应的代换:
1µ
−。
5.4单元等效节点载荷列阵
5.4.1平面问题节点力的一些约定:
1、分别用U和V加节点号下角字表示某节点的水平与垂直节点力分量;
2、节点力的方向,假设节点对单元的节点力为正,方向取为坐标轴正方向;
单元对节点的节点力为负,方向取为坐标轴负方向。
5.4.2关于载荷移置
载荷移置就是把弹性体所受的各种形式的载荷按等效原则移置到节点上成为节点载荷。
如果弹性体所受的载荷全都是集中力,那么把所有集中力作用点都取为节点,就不存在移置问题,集中力就是节点载荷,把它们分解到x轴和y轴即可。
但实际工程问题往往受有分布的面力和体力作用,在这种情况下就必须把面力和体力向节点移置。
就是对集中力,若其作用点未被取为节点,也要向节点移置。
单元载荷的移置是载荷移置的第一步也是最主要的一步,只要有了单元载荷的移置结果,总体节点载荷就容易计算了。
要把载荷移置到节点上,首先要遵循静力等效原则,即原载荷与节点载荷在任意轴上的投影之和以及对任一轴的力矩之和都相等。
对连续弹性体来讲,静力等效原则实质上就是虚功等效原则,即单元的原载荷与移置后的节点载荷在任何虚位移上的虚功都相等。
必须注意的是,这里的单元虚位移与所设选的单元位移函数要采用同一位移模式。
我们前面所选的单元位移函数是线性模式,所以单元虚位移也都采取线性模式。
设单元所受载荷有:
体力{}f,面力{}T;
移置后的节点载荷为{}e
P;
在外载荷作用下
节点产生的任意虚位移为{}e
*δ,则单元内各点的任意虚位移为{}[]{}e
Nu**
δ=。
根据静力
等效原则有:
{}{}{}{}{}{}∫∫+=S
VTe
dSTudVfuP***δ
即
{}{}{}{}{}{}{}{}∫∫+=S
V
dSTNdVfNP***δδδ
由于虚位移
{}e
*δ是任意的,所以有:
{}{}{}{}∫∫+=⎪
⎪⎬⎫⎪⎪⎪⎪⎩⎪⎪⎪
⎪⎨⎧=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧=ST
VTmymxjyjxiyixmjie
dSTNdVfNpppppppppP(5-95.4.3载荷移置举例
1、单元某点处的集中载荷
设在三角形单元的(x,y点受到集中载荷{}[
yx
PPP=,则:
⎪⎭⎪⎪⎪⎪⎬⎫
⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧=⎪⎪⎪⎪⎭⎪⎪
⎪⎪⎬⎫⎪⎪
⎪⎪⎩⎪⎪⎪⎪⎨⎧=ymxmyjxjyiximymxjyjxiyixe
PNPNPNPNPNPNppppppP2、均质等厚单元的自重
设单元厚度为t,材料的比重为γ,则单元自重的合力为StGγ=,并作用在单元形心(即三角形的重心c上,如图5-6所示,我们先求节点载荷iyp。
为了简化计算,我们假设单元发生这样的虚位移:
节点i只沿y方向移动一个单位,即
01**==iiuvδδ,,而其余两个节点都不移动,
0*
***====mmjjvuvuδδδδ。
由于采取的是线性
位移模式,因此单元形心将产生1/3的虚位移,所以由静力等效原则可知:
3
1⋅−=⋅Gpiy
故得
mx
P
3S
tGpiyγ−
=−
=按照上述分析方法同样可得
。
0,3
===−
==mxjxixmyjypppS
tppγ综上,写成矩阵形式为
⎪⎪⎪⎪⎭
⎪
⎪⎪⎪⎬⎫
⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧−
=⎪⎪⎪⎪
⎭⎪⎪⎪⎪⎬⎫⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧=1010103StppppppPmymxjyjxiyixe
γ(5-10由此可见,像重力这种在单元内均匀分布的体力,把它平均地移置到单元的三个节点上即可。
3、单元一边上受有均布侧压
设单元在jm边上受有均匀分布的垂直压力,jm边与x轴的夹角为θ,如图5-7所示,求其节点载荷。
我们先来求节点载荷分量mxp。
为了简化计算同样假设01*
*==mmvuδδ,,其余两个节点都不移动,即0*
*
====jjiivuvuδδδδ。
由图可知,原载荷的两个分量为
sincosjmxjmm
jyjmyyqqq
lxxqqqlθθ−⎫
=−=⎪⎪
⎬−⎪
==⎪⎭
(5-11
其中jml为单元该边的长度。
设s是自节点m到jm边上任一点的距离,由于单元虚位移也采取线性模式,所以当节点m发生水平单位虚位移时,s点的虚位移是
图5-7
根据静力等效原则有
∫+=⋅jm
yxxmtdssvqsuqp((1δδ
imjljmjm
mjljm
jmx
mxtbqyytqdssllyyqt
tdslslqpjm
jm
(2(0
=−=
−−=−=∫
∫按照上述类似的分析方法,可得其它节点力如下:
综上,写成矩阵形式为
⎪⎪⎪⎪⎭⎪⎪⎪⎪⎬⎫⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧=⎪⎪⎪⎪⎭
⎪⎪⎪⎪⎬⎫⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧=iiiimymxjyjxiyixe
cbcbqtppppppP002(5-12
由此可见,对于面力的移置,在采取线性位移模式的情况下,只要把原载荷静力等效地移置到受力边的两个节点上即可,均布压力如此,对于非均布压力也是这样。
3、单元一边上受有x方向的均布载荷
设单元的ij边上受有沿x方向的集中力P,其作用点离i及j点的距离分别为il及jl,则:
((0jmjmlsuslvsδδ−⎧=⎪
⎪=⎩
(22(22(2200
mymjijxjmi
jymji
ixiyqq
pxxtcqqpyytbqqpxxtcpp=−==−==−===
{}⎪⎪⎪⎪⎭
⎬⎫
⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧=⎪⎪⎪⎪⎭
⎪⎪⎪⎪⎬⎫⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧=0000ll
llPppppppPijmymxjyjxiyixe
其中,l为ij边长度。
由此推导:
设单元在ij边上受有均匀分布的x向压力,如图5-8所示,求其节点载荷。
边界上的作用力写成:
Q=⎭
⎬⎫⎩⎨⎧0q
qlttdssuqpijxix2
1(1==⋅∫δ
0====jmimjyiypppp
由以上同样的分析方法可以得出等效节点载荷为:
{}⎪⎪⎪⎪⎭⎪⎪⎪⎪⎬
⎫⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧=⎪⎪
⎪⎪⎪⎪⎬⎫⎪⎪⎪⎪
⎩⎪⎪⎪⎪⎨⎧=0001012qtlppppppPmymxjyjxiyixe
(5-13
4、单元一边上受有x