力学WORD教案CHAPTER3 弹性力学平面问题的有限元法.docx
《力学WORD教案CHAPTER3 弹性力学平面问题的有限元法.docx》由会员分享,可在线阅读,更多相关《力学WORD教案CHAPTER3 弹性力学平面问题的有限元法.docx(31页珍藏版)》请在冰豆网上搜索。
力学WORD教案CHAPTER3弹性力学平面问题的有限元法
力学WORD教案-CHAPTER3弹性力学平面问题的有限元法
3弹性力学平面问题的有限元法
本章包括以下的内容:
3.1弹性力学平面问题的基本方程
3.2单元位移函数
3.3单元载荷移置
3.4单元刚度矩阵
3.5单元刚度矩阵的性质与物理意义
3.6整体分析
3.7约束条件的处理
3.8整体刚度矩阵的特点与存储方法
3.9方程组解法
3.1弹性力学平面问题的基本方程
弹性力学是研究弹性体在约束和外载荷作用下应力和变形分布规律的一门学科。
在弹性力学中针对微小的单元体建立基本方程,把复杂形状弹性体的受力和变形分析问题归结为偏微分方程组的边值问题。
弹性力学的基本方程包括平衡方程、几何方程、物理方程。
弹性力学的基本假定如下:
1)完全弹性,2)连续,3)均匀,4)各向同性,5)小变形。
3.1.1基本变量
弹性力学中的基本变量为体力、面力、应力、位移、应变,各自的定义如下。
体力
体力是分布在物体体积内的力,例如重力和惯性力。
面力
面力是分布在物体表面上的力,例如接触压力、流体压力。
应力
物体受到约束和外力作用,其内部将产生内力。
物体内某一点的内力就是应力。
图3.1
如图3.1假想用通过物体内任意一点p的一个截面mn将物理分为?
、?
两部分。
将部分?
撇开,根据力的平衡原则,部分?
将在截面mn上作用一定的内力。
在mn截面上取包含
A,Ap点的微小面积,作用于面积上的内力为。
Q
A令无限减小而趋于p点时,的极限S就是物体在p点的应力。
Q
Qlim,S,A,0,A
应力S在其作用截面上的法向分量称为正应力,用σ表示;在作用截面上的切向分量称为剪应力,用τ表示。
显然,点p在不同截面上的应力是不同的。
为分析点p的应力状态,即通过p点的各个截面上的应力的大小和方向,在p点取出的一个平行六面体,六面体的各楞边平行于坐标轴。
图3.2
将每个上的应力分解为一个正应力和两个剪应力,分别与三个坐标轴平行。
用六面体表面的应力分量来表示p点的应力状态。
应力分量的下标约定如下:
第一个下标表示应力的作用面,第二个下标表示应力的作用方向。
,第一个下标x表示剪应力作用在垂直于X轴的面上,第二个下标y表示剪应力指xy
向Y轴方向。
表示正应力作用于垂直于X轴正应力由于作用表面与作用方向垂直,用一个下标。
x
的面上,指向X轴方向。
应力分量的方向定义如下:
如果某截面上的外法线是沿坐标轴的正方向,这个截面上的应力分量以沿坐标轴正方向为正;
如果某截面上的外法线是沿坐标轴的负方向,这个截面上的应力分量以沿坐标轴负方向为正。
剪应力互等:
,,,,,,,,,,xyyxyzzyzxxz
,,,,,物体内任意一点的应力状态可以用六个独立的应力分量、、、、、yxyyzxzxz
来表示。
位移
位移就是位置的移动。
物体内任意一点的位移,用位移在x,y,z坐标轴上的投影u、v、w表示。
应变
物体的形状改变可以归结为长度和角度的改变。
各线段的单位长度的伸缩,称为正应变,用ε表示。
两个垂直线段之间的直角的改变,用弧度表示,称为剪应变,用γ表示。
物体内任意一点的变形,可以用六个应变分量表示。
、,、,、,、,、,xyzxyyzzx
3.1.2平面应力和平面应变问题
弹性体在满足一定条件时,其变形和应力的分布规律可以用在某一平面内的变形和应力的分布规律来代替,这类问题称为平面问题。
平面问题分为平面应力问题和平面应变问题。
1)平面应力问题
设有很薄的等厚薄板,只在板边上受到平行于板面并且不沿厚度变化的面力,体力也平行于板面且不沿厚度变化。
图3.3
设板的厚度为t,在板面上:
tt,,,,,,0,,0,,,,,,0tzzxzyz,,z,,z,,222
由于平板很薄,外力不沿厚度变化,因此在整块板上有,
,0,,0,,,,0zxzzy
、,、,剩下平行于XY平面的三个应力分量未知。
xyxy
2)平面应变问题
设有很长的柱形体,支承情况不沿长度变化,在柱面上受到平行于横截面而且不沿长度变化的面力,体力也如此分布。
图3.4
以柱体的任一横截面为XY平面,任一纵线为Z轴。
假定该柱体为无限长,则任一截面都可以看作对称面。
由对称性,
,0,,,,0w,0zxzy
,0由于没有Z方向的位移,Z方向的应变。
z
、,、,,物体在Z方向处于自平衡状未知量为平行于XY平面的三个应力分量xyxy
态。
3.1.3平衡方程
弹性力学中,在物体中取出一个微小单元体建立平衡方程。
平衡方程代表了力的平衡关系,建立了应力分量和体力分量之间的关系。
对于平面问题,在物体内的任意一点有,
,,,yxx,,X,0,x,y(3-1),,,,yxy,,Y,0,y,x
3.1.4几何方程
由几何方程可以得到位移和变形之间的关系。
对于平面问题,在物体内的任意一点有,
u,,x,x
v,(3-2),y,y
u,v,,,xy,y,x
刚体位移
由位移u=0,v=0可以得到应变分量为零,反过来,应变分量为零则位移分量不为零。
应变分量为零时的位移称为刚体位移。
刚体位移代表了物体在平面内的移动和转动。
u,0,x
v由,0,y
u,v,,0,y,x
可以得到刚体位移为以下形式,
u,u,y0v,v,,x0
u,v,0,,0由可得,,x,y
u,f(y),v,f(x)12
u,v,,0f,f将代入可得,12,y,x
df(y)df(x)12,,,,dydx
积分后得到,
f(y),u,y10f(x),v,,x20
得到位移分量,
u,u,y0v,v,,x0
u,0,v,0,,,0u当时,物体内任意一点都沿x方向移动相同的距离,可见代000
表物体在x方向上的刚体平移。
u,0,v,0,,,0v时,物体内任意一点都沿y方向移动相同的距离,可见代表当000
物体在y方向上的刚体平移。
u,0,v,0,,,0当时,可以假定,此时的物体内任意一点P(x,y)的位,,000
移分量为,
u,,,y,v,,x
P点位移与y轴的夹角为α,
yytg,,,,tg,xx,
P点合成位移为,
222222u,v,(,,y),,(x),,x,y,,r
r为P点到原点的距离,可见ω代表物体绕z轴的刚体转动。
3.1.5物理方程
弹性力学平面问题的物理方程由广义虎克定律得到。
1)平面应力问题的物理方程
1,,,,,,,,xxyE
1,,(3-3),,,,,,yyxE
2(1,),,,xyxyE
平面应力问题有,
,0z
,,,,,,,,zxyE
2)平面应变问题的物理方程
2,,,,1,,,,,,,,xxy,,E1,,,,
2,,1,,,,,(3-4),,,,,yyx,,E1,,,,
2(1,),,,xyxyE
平面应变问题有,
,0z
,,,,,,,,zxy
E,在平面应力问题的物理方程中,将E替换为、替换为,可以得到平面21,,1,,
E(1,2),、替换为应变问题的物理方程;在平面应变问题的物理方程中,将E替换为2(1,),,,可以得到平面应力问题的物理方程。
1,,
图3.5
求解弹性力学平面问题,可以归结为在任意形状的平面区域内已知控制方程、在位
SS移边界上约束已知、在应力边界上受力条件已知的边值问题。
然后以应力分量为基本u,
未知量求解,或以位移作为基本未知量求解。
如果以位移作为未知量求解,求出位移后,由几何方程可以计算出应变分量,得到物体的变形情况;再由物理方程计算出应力分量,得到物体的内力分布,就完成了对弹性力学平面问题的分析。
3.2单元位移函数
根据有限元法的基本思路,将弹性体离散成有限个单元体的组合,以结点的位移作为未知量。
弹性体内实际的位移分布可以用单元内的位移分布函数来分块近似地表示。
在单元内的位移变化可以假定一个函数来表示,这个函数称为单元位移函数、或单元位移模式。
对于弹性力学平面问题,单元位移函数可以用多项式表示,
22u,a,ax,ay,ax,axy,ay,...123456
22(3-5)v,b,bx,by,bx,bxy,by,...123456
多项式中包含的项数越多,就越接近实际的位移分布,越精确。
具体取多项,由单元形式来确定。
即以结点位移来确定位移函数中的待定系数。
图3.6
(x,y)如图3.6所示的3结点三角形单元,结点I、J、M的坐标分别为、、(x,y)iijj
uvu(x,y)uvv,结点位移分别为、、、、、。
六个节点位移只能确定六个多jjmmmiim
项式的系数,所以3结点三角形单元的位移函数如下,
u,a,ax,ay,123(3-6),v,a,ax,ay456,
将3个结点上的坐标和位移分量代入公式(3-6)就可以将六个待定系数用结点坐标和位移分量表示出来。
将水平位移分量和结点坐标代入(3-6)中的第一式,
u,a,ax,ayi12i3i
u,a,ax,ayj12j3j
u,a,ax,aym12m3m
写成矩阵形式,
u1xya,,,,,,iii1,,,,,,u1xya,(3-7),,,,jjj2,,,,,,,,u1xyammm3,,,,,,
xy1,,ii,,,,1xy,T令,jj,,
,xy1mm,,
au,,,,1i,,,,,1a,,Tu,则有(3-8),,,,2j
,,,au3m,,,,
*[T],1[T],T
,A为三角形单元的面积。
T,2A
[T]的伴随矩阵为,
T,,xyxyyyxx,,,jmmjjmmj,,*xyxyyyxx,,T,,,,(3-9)miimmiim,,
,xyxyyyxx,,,ijjiijji,,
T,,abcaaa,,iiiijm,,,,*abcbbb,,[T]令(3-10)jjjijm,,,,
,,,abccccmmmijm,,,,
,,aaaau,,,,1ijmi1,,,,,,abbbu,则(3-11),,,,2ijmj,,2A,,,,,,acccu3ijmm,,,,,,
同样,将垂直位移分量与结点坐标代入公式(3-6)中的第二式,可得,
,,aaaav,,,,4ijmi1,,,,,,abbbv,(3-12),,,,5ijmj,,2A,,,,,,acccv6ijmm,,,,,,
将(3-11)、(3-12)代回(3-6)整理后可得,
1u,[(a,bx,cy)u,(a,bx,cy)u,(a,bx,cy)u]iiiijjjjmmmm2A
1v,[(a,bx,cy)v,(a,bx,cy)v,(a,bx,cy)v]iiiijjjjmmmm2A
1令(下标i,j,m轮换)N,(a,bx,cy)iiii2A
u,,i,,vi,,
N0N0N0,,uu,,,,,,ijmj可得(3-13),,,,,,,0N0N0Nvvijm,,j,,,,
,um,,,,vm,,
u,,单元内的位移记为f,,,,,v,,
u,,i,,vi,,,,,i,,u,,,,ej,,,,单元的结点位移记为,,,,,,jvj,,,,,m,,,,um,,,,vm,,
单元内的位移函数可以简写成,
e,,,,,,f,N,(3-14)
把[N]称为形态矩阵,N称为形态函数。
i
选择单元位移函数应满足以下条件:
1)反映单元的刚体位移与常量应变。
2)相邻单元在公共边界上的位移连续,即单元之间不能重叠,也不能脱离。
由(3-6)可以将单元位移表示成以下的形式,
aaaa,,5353uaaxyy,,,,1222
aaaa,,5353vaayxx,,,,4622
反映了刚体位移和常应变。
单元位移函数是线性插值函数,因此单元边界上各点的位移可以由两个结点的位移完全
确定。
两个单元的边界共用两个结点,所以边界上的位移连续。
形态函数N具有以下性质:
i
1)在单元结点上形态函数的值为1或为0。
2)在单元中的任意一点上,三个形态函数之和等于1。
用来计算三角形面积时,要注意单元结点的排列顺序,当三个结点i,j,m取逆时T
11针顺序时,;当三个结点i,j,m取顺时针顺序时,。
A,T,0A,T,022
例题:
如图3.7所示等腰三角形单元,求其形态矩阵[N]。
a,xy,xy解:
由ijmmj
b,y,yijm
c,x,ximj
在公式中轮换下标可以计算得
,,a,xy,xy,0,0,0,a,0b,y,y,a,0,aijmmjijm
c,x,x,0,0imj
b,y,y,0,0,0,,a,xy,xy,0,0,a,0,0jmijmiim
c,x,x,a,0,ajim
2,a,xy,xy,a,a,0,0,ab,y,y,0,a,,amijjimij
c,x,x,0,a,,amji
2aA,三角形积为2
形态函数为
x11Nabxcyax,(,,),(0,,0),iiii2Aaa2
y11Nabxcyay,(,,),(0,0,),jjjj2Aaa2
xy112Nabxcyaaxay,(,,),(,,),1,,mmmm2Aaaa2
形态矩阵为
xyxy,,001,,0,,aaaa[N],,,xyxy,,0001,,aaaa,,
三角形面积的计算公式可得,
1xy1a0ii1112,,,A1xy10aajj2221xy100mm
如果把三个结点按顺时针方向排列,即i(a,0),j(0,0),m(0,a)
1xy1a0ii1112,,,,A1xy100ajj2221xy10amm
3.3单元载荷移置
有限元法的求解对象是单元的组合体,因此作用在弹性体上的外力,需要移置到相应的
结点上成为结点载荷。
载荷移置要满足静力等效原则。
静力等效是指原载荷与结点载荷在任
意虚位移上做的虚功相等。
e*单元的虚位移可以用结点的虚位移表示为,,,,
e**(3-15),,,,f,[N],
X,,i,,Yi,,
,X,,ej,,令结点载荷为R,,,Yj,,
,Xm,,,,Ym,,
1)集中力的移置
P,,x如图3.7所示,在单元内任意一点作用集中力,,P,,,Py,,
图3.8
由虚功相等可得,
TT*ee*eT,,,,{,}{R},{,}[N]{P}
eT{R},[N]{P}由于虚位移是任意的,则(3-16)例题1:
在均质、等厚的三角形单元ijm的任意一点p(x,y)上作用有集中载荷。
pp
0XN,,,,ii,,,,0YNii,,,,
P,,,,,,0XN,,xjj,,,,,,,P0YNyjj,,,,,,
,,,XN0mm,,,,,,0YN,,mm,,,,(,)xypp
2)体力的移置
x,,令单元所受的均匀分布体力为{p},,,y,,
由虚功相等可得,
TT*ee*T,,,,R,Nptdxdy{,}{}{,}[]{},,
eT{}[]{}R,Nptdxdy(3-17),,
3)分布面力的移置
T,,P,[X,Y]设在单元的边上分布有面力,同样可以得到结点载荷,
eT{R}[N]{P}tds,(3-18),s
例题2:
设有均质、等厚的三角形单元ijm,受到沿y方向的重力载荷q的作用。
求均布体y
力移置到各结点的载荷。
XN0,,,,ii,,,,Y0Nii,,,,
0,,,,0XN,,,,jj,tdxdy,,,,,,,,q0YNyjj,,,,,,
,,,XN0mm,,,,,,YN0,,mm,,,,
X,0,X,0,X,0ijm
Y,Nqtdxdy,qtNdxdyiiyyi,,,,
1,(,,)Ndxdyabxcydxdyiiii,,,,2A111,[,,],(,,),aAbAxcAyAabxcyAiiciciicic223AA
1,YqAtiy3
11,,,同理,YqAtYqAtjymy33
例题3:
在均质、等厚的三角形单元ijm的ij边上作用有沿x方向按三角形分布的载荷,求
移置后的结点载荷。
XN0,,,,ii,,,,Y0Nii,,,,
,,,XN0q,,,,jjxtdxdy,,,,,,,,Y0N0,,jjs,,,,
,,,XN0mm,,,,,,Y0N,,mm,,,,
取局部坐标s,在i点s=0,在j点s=l,L为ij边的长度。
在ij边上,以局部坐标表示
的插值函数为,
ssN,0,,NN,1,,mijLL
s载荷为qq,xL
L231ssssL
(1)()X,,qtds,qt,,qtLi2,026LLL3L0
L31sssLX,qtds,qt,qtLj2,03LL3L0
3.4单元刚度矩阵
根据单元的位移函数,
u,,i,,vi,,
N0N0N0,,uu,,,,,,ijmj,,,,,,,0N0N0Nvvijm,,j,,,,
,um,,,,vm,,
由几何方程可以得到单元的应变表达式,
u,,i,,,u,,,,vi,,,x,,b0b0b0,,ijm,,u,v1,,,,,,j,,,,,0c0c0c(3-19),,,,ijm,,v,y2Aj,,,,,,cbcbcbiijjmm,,,u,v,,,,um,,,,,,y,x,,,,vm,,
e{,},[B]{,}记为,[B]矩阵称为几何矩阵。
,,,B,BBB[B]矩阵可以表示为分块矩阵的形式ijm
b0,,i1,,,,B,0c(3-20)ii,,2A,,cbii,,
由物理方程,可以得到单元的应力表达式,
e,,,,,,,,,,,,,,D,,DB,(3-21)
[D]称为弹性矩阵,对于平面应力问题,
,,,10,,E,,,,10,,D2,,,(1,)1,,,,002,,
,,,,,S,DB定义为应力矩阵。
将应力矩阵分块表示为,
,,,S,SSSijm
,,,bc,,iiE,,,,,,,,,(3-22)S,DB,bciiii2,,,2A(1,)1,,1,,,,cbii,,22,,
应用虚功原理可以建立单元结点位移与结点力的关系矩阵,单元刚度矩阵。
虚功原理:
在外力作用下处于平衡状态的弹性体,如果发生了虚位移,则所有外力在虚
位移上做的虚功等于内应力在虚应变上做的虚功。
eT,,,,单元的结点力记为F,UVUVUViijjmm
e**单元的虚应变为,,,,,,,,B,
单元的外力虚功为,
Tee*,,,,,,,F
单元的内力虚功为,
T*,,,,tdxdy,,,,
由虚功原理可得,
TeTe**,,,,,,,,,,,,,F,tdxdy(3-23),,
TTee**T*T,,,,,,(,,B,,,,,),,[B]
ee,,,,,,,,[S],,[D][B],
TTeeee**T,,,,,,,,,,,F,,,,[B][D][B]tdxdy,,,
eeT(3-24),,,,,F[B][D][B]tdxdy,,,
eT定义为单元刚度矩阵。
K,BDBtdxdy,,[][][],,
在3结点等厚三角形单元中[B]和[D]的分量均为常量,则单元刚度矩阵可以表示为,
eT(3-25),,K,[B][D][B]tA
单元刚度矩阵表示为分块矩阵:
,,[K][K][K]iiijim,,eKKKK,[][][],,jijjjm,,
,[K][K][K]mimjmm,,
T[K],[B][D][B](3-26)rsrs
KK,,rx,sxrx,syK[],(3-27),,rsKKry,sxry,sy,,
3.5单元刚度矩阵的性质与物理意义
(一)单元刚度矩阵的物理意义
eT假设单元的结点位移如下:
,,,[100000]
eee{F},[K]{,}由,得到结点力如下:
KU,,,,ix,ixi,,,,KViy,ixi,,,,
,K,,U,,,,jx,ixj,(3-28),,,,KVjy,ixj,,,,
,,,KUmx,ixm,,,,K,,V,,my,ixm,,,,
表示i结点在水平方向产生单位位移时,在结点i的水平方向上需要施加的结点力。
Kix,ix
K表示i结点在水平方向产生单位位移时,在结点i的垂直方向上需要施加的结点力。
iy,ix
选择不同的单元结点位移,可以得到单元刚度矩阵中每个元素的物理含义:
表示s结点在水平方向产生单位位移时,在结点r的水平方向上需要施加的结点力。
Krx,sx
K表示s结点在水平方向产生单位位移时,在结点r的垂直方向上需要施加的结点力。
ry,sx
K表示s结点在垂直方向产生单位位移时,在结点r的水平方向上需要施加的结点力。
rx,sy
K表示s结点在垂直方向产生单位位移时,在结点r的垂直方向上需要施加的结点力。
ry,sy
因此单元刚度矩阵中每个元素都可以理解为刚度系数,即在结点产生单位位移时需要施
加的力。
(二)单元刚度矩阵的性质
1)对称性
利用分块矩阵的性质证明如下:
T[K],[B][D][B]rsrs
T[K],[B][D][B]srsr
TTTTTT[K],([B][D][B]),[B][D][B],[B][D][B],[K]srsrrsrsrs
eeT[K],([K])即
2)奇异性
eK,0即单元刚度矩阵的行列式为零,。
eT{,},[101010]将定单元产生了x方向的刚体移动,,此时对应的单元结点力为零。
01,,,,
,,,00,,,,
,,,01,,,,eK,[],,,,00,,,,
,,,01,,,,,,,,00,,,,
可以得到,在单元刚度矩阵中1,3,5列中对应行的系数相加为零,由行列式的性质可
eK,0。
知,
eT{,},[010101]同样如果假定单元产生了y方向上的刚体位移,可以得到,在单元刚度矩阵中2,4,6列中对应行的系数相加为零。
3.6整体分析
得到了单元刚度矩阵后,要将单元组成一个整体结构,根据结点载荷平衡的原则进行分析,即整体分析。
整体分析包括以下4个步骤:
1)建立整体刚度矩阵,
2)根据支承条件修改整体刚度矩阵,
3)解方程组,求出结点的位移,
4)根据结点位移,求出单元的应变和应力。
在这里把结点位移作为基本未知量求解。
如何得到整体刚度矩阵,基本方法是刚度集成法,即整体刚度矩阵是单元刚度矩阵的集成。
图3.9
如图3.9所示,一个划分为6个结点、4个单元的结构。
得到了每个单元的单元刚度矩阵后,要集成为整体刚度矩阵。
3.6.1刚度集成法的物理意义
由单元刚度矩阵的物理意义可知,单元刚度矩阵的系数是由单元结点产生单位位移时引起的单元结点力。
在如图3.9所示的结构中,使结点3产生单位位移时,在单元
(1)中的结点2上引起结点力。
由于结点2、3同时属于单元
(1)、(3),在单元
(2)中的结点2上同样也引起结点力,因此,在整体结构中当结点3产生位移时,结点2上的结点力应该是单元
(1)、
(2)在结点2上的结点力的叠加。
刚体集成法即结构中的结点力是相关单元结点力的叠加,整体刚度矩阵的系数是相关单元的单元