有限元学习心得.docx
- 文档编号:15229105
- 上传时间:2023-07-02
- 格式:DOCX
- 页数:27
- 大小:263.18KB
有限元学习心得.docx
《有限元学习心得.docx》由会员分享,可在线阅读,更多相关《有限元学习心得.docx(27页珍藏版)》请在冰点文库上搜索。
有限元学习心得
有限元学习心得
吴清鸽车辆工程50110802411
短短八周的有限元课已经结束。
关于有限元,我一直停留在一个很模糊的概念。
我知道这是一个各个领域都必须涉及的点,只要有关于CAE分析的,几乎都要涉
及有限元。
总体来说,这是一门非常重要又有点难度的课程。
有限元方法(finiteelementmethod)或有限元分析(finiteelementanalysis),是求取复杂微分方程近似解的一种非常有效的工具,是现代数字化科技的一种重要基础性原理。
将它用于在科学研究中,可成为探究物质客观规律的先进手段。
将它应用于工程技术中,可成为工程设计和分析的可靠工具。
本课程教学基本内容有固体力学和结构力学简介;有限元法基础;桁架、梁、刚架、二维固体、板和壳、三维固体的有限元法;建模技术;热传导问题的有限元分析;PATRAN软件的使用.
通过有限元分析课程学习使我了解和掌握了一些有限元知识:
1.简要了解二维和三维固体以及桁架、梁和板结构的三组基本力学方程,即表示位移-应变关系的几何方程,表示应力-应变关系的本构方程和表示内力-外力关系的平衡方程。
2.了解利用能量法形成有限元离散系统方程的基本原理,即哈密尔顿原理。
掌握有限元分
析的基本方法及步骤,包括域的离散、位移插值、构造形函数、单元有限元方程的建立、坐标变换、整体有限元方程的组装、整体有限元方程的求解技术。
3.具体深入的了解并掌握桁架结构、梁结构、刚架结构、二维固体、板和壳结构、三维固体的有限元法分析技术,包括他们具体的形函数构造,应变矩阵,局部坐标系和整体坐标系中的单元矩阵。
各种结构的实例研究。
4.了解并掌握建立高质量建模所涉及的各种关键技术。
包括单元类型的选择,单元畸形的限制,不同阶数单元混用时网格的协调性问题,对称性的应用(平面对称、轴对称、旋转对称、重复对称),由多点约束方程形成刚域及应用(模拟偏移、不同自由度单元的连接、网格协调性的施加)等,以及多点约束方程的求解。
以PATRAN有限元通用软件为例了解一般商业有限元软件的组成及结构。
掌握PATRAN软件的基本使用。
利用PATRAN软件上机实践完成两个上机练习:
刚架结构有限元分析和三维固体有限元分析。
课程的具体学习内容:
内容:
1.三节点三角形单元:
单元分析、总刚度矩阵组装、引入约束条件修正总刚度矩阵、载荷移置、方程求解;
2.四边形单元分析、四节点四面体单元分析、八节点六面体单元分析;
3.其他常用单元形函数、自由度
1、三节点三角形单元
1.1.单元分析
1.1.1分析步骤
单元分析的任务是建立单元平衡方程,形成单元刚度矩阵。
不失一般性,从图1-1三角形离散结构中任取一个单元,设单元编号为e,单元节点按右手法则
顺序编号为i,j,m,在定义的坐标系xOy中,节点坐标分别为(xi+yi),(xj+yj),
其中,转换矩阵称为单元刚度矩阵。
单元分析的主要目的就是要求出单
元刚度矩阵
单元分析的步骤可表示如下:
结点位移
(1)
a二二三■
内部各点位移
(2)
应变
⑶
———■
应力
⑷
as■
结点力
单元分析
1
1.1.2位移模式和形函数
对于平面问题,单元任意一点的位移可用位移分量u,v描述,他们是坐标
x,y的函数。
假定三节点单元的位移函数为x,y的线性函数,六个节点位移只能确定六个多项式的系数,所以平面问题的3结点三角形单元的位移函数如下:
ua1a2xa3y
va4a5xa6y
所选用的这个位移函数,将单元内部任一点的位移定为座标的线性函数,位移模式很简单。
位移函数写成矩阵形式为:
Uf
V
a2
1xy000a3
0001xya4
3s
将水平位移分量和结点坐写成矩阵:
代入位移函数第一式:
Ui
a1
a2Xi
asYi
Ui
1
Xi
Yi
4
Uj
q
a2Xj
a3Yj
Uj
1
Xj
Vj
a2
Um
令则有
%
a2Xm
a3Ym
Um
1
Xm
Ym
a3
1
Xi
yi
a1
Ui
1
Xj
yj
T
a2
T
1
Uj
1
Xm
ym
a3
Um
A为三角形单元
令则有[T]
ai
b
Ci
aiaj
am
a1
1
ai
aj
am
Ui
aj
bj
Cj
bibj
bm
a2
b
bj
bm
Uj
2A
am
Cm
CiCj
Cm
a3
G
Cj
Cm
Um
T2A
XjYm
*
XmYj
YjYm
XmXj
[T]的伴随矩阵为
TXmYi
XiYm
YmYi
XiXm
XiYj
XjYi
YiYj
XjX
同样,将垂直位移分量与结点坐标代入位移插值公式:
最终确定六个待定系数
ai
a2
a3
1
2A
ai
b
c
aj
bj
Cj
—[(ai
2A
2A
biX
bx
a4
ai
aj
am
V
1
a5
bi
bj
bm
Vj
2A
a6
C
C
Cm
Vm
amUi
a4
A
ai
aj
amVi
bmUj
as
1
b
bj
bmVj
2A
CmUm
a6
Cj
CmVm
CiY)Ui(ajbjXCj
Y)Uj
(am
bmXCmY)Um]
CY)V(ajbjX
CY)Vj
(am
bmX
CmY)
Ui
V
Ni0Nj0Nm
0
Uj
0Ni0Nj0
Nm
Vi
Vm]
Um
Vm
换)
Ui
Vi
i
Uj
j
Vj
m
Um
[N]称为形态矩阵,Ni称为位移的形态函数
1.1.3位移函数的收敛性
选择单元位移函数时,应当保证有限元法解答的收敛性,即当网格逐渐加密时,有限元法的解答应当收敛于问题的正确解答。
因此,选用的位移模式应当满足下列两方面的条件:
(1)必须能反映单元的刚体位移和常量应变。
6个参数a1到a6反映了三个刚体位移和三个常量应变。
(2)必须保证相邻单元在公共边界处的位移连续性。
(线性函数的特性)
1.1.4应变矩阵和应力矩阵
利用几何方程、物理方程,实现用结点位移表示单元的应变和单元的应力
用结点位移表示单元的应变的表达式为:
b
0
b
j
0
bm
0
1
m
Uj
2A
0
Ci
0
q
0
Cm
Vj
c
bi
Cj
bj
Cm
bm
Um
Vm
1
b
0
Bi
0
Ci
2A
Ci
bi
U
x
V
yuV
yx
{}[B]{}e
Ui
Vi
e
SSiSjSm
由物理方程,可以得到单元的应力表达式:
BBiBjBm
[B]矩阵称为几何矩阵
为应力矩阵
DBi
E
2A(12)
1.1.5单元刚度矩阵
Ci
Ci
1
eT
Fe
讨论单元内部的应力与单元的结点力的关系,导出用结点位移表示结点力的表
U
i
达式。
由应力推算结点力,需要利用平衡方租。
0用虚dyd程表示出平衡方程。
考虑上图三角形单元的实际受力,任意虚设位移,节点位移结点力和内部应力为:
与内部应变为:
x
*
y
gxy
Um
代入虚功方程
eTF
Fe
接上式,将应力用结点位移表示出
[B]t
(T
eTt
[B]t
(Ttdxdy
DBS
tdxdy
有令则
Fe
Ke
Fe
微小矩形的内力虚功为
[B]T[D][B]tdxdyS
[B]T[D][B]tdxdy
K:
e
建立了单元的结点力与结点位移之间的关系,称为单元刚度矩阵。
它
是6*6矩阵,其元素表示该单元的各结点沿坐标方向发生单位位移时引起的结点力,它决定于该单元的形状、大小、方位和弹性常数,而与单元的位置无关,即不随单元或坐标轴的平行移动而改变。
1.2总刚度矩阵组装
整体刚度矩阵K是单元刚度矩阵ke的集成。
1、刚度集成法的物理概念:
冈寸度矩阵中的元素是刚度系数,即由单位结点位移引起的结点力。
2、刚度矩阵的集成规则:
e
先对每个单元求出单元刚度矩阵k,然后将其中的每个子块kj送到结构刚度矩阵中的对应位置上去,进行迭加之后即得出结构刚度矩阵[K]的子块,从而得出结构刚度矩阵e[K]。
关键是如何找出k中的子块在[K]中的对应位置。
这需要了解单元中的结点编码与结构中的结点编码之间的对应关系。
结构中的结点编码称为结点的总码,各个单元的三个结点又按逆时针方向编为i,j,m,称为结点的局部码。
单元刚度矩阵中的子块是按结点的局部码排列的,而结构刚度矩阵中的子块
[Kjj]
(2)
(:
[Kjm]
)[Kji](
2)
[Kmj](:
)
【Kmm】(
wr
2)
[Kij]
(2)
【Kim】"
)[kb](
2)
T
a
j二
a
j2
2
3
4
5
m2
1
2
j2
是按结点的总码排列的。
因此,在单元刚度矩阵中,把结点的局部码换成总码,并把其中的子块按照总码次序重新排列。
以单元②为例,局部码i,j,m对应于总码5,2,4,因此k⑵中的子块按照总码重新
排列后,得出扩大矩阵为上图所示:
⑴⑶⑷
、K、K
将各单元的扩大矩阵迭加,即得出结构刚度矩阵[K]:
KK⑴K
(2)K⑶K⑷K(e)
集成规则包含搬家和迭加两个环节:
.e
1、将单元刚度矩阵k中的子块搬家,得出单元的扩大刚度矩阵
2、将各单元的扩大刚度矩阵Ke迭加,得出结构刚度矩阵[K]。
Ke
j总码
ji1
mi
j22
13
ii
m33
j4
m24
i2
j35
m4
l46
[Kjj](i)
[Kim](i)
[KJ(i)
[Kmm]U)[kJ2)[Kii](3)
[Kml]⑴
[Km](3)
[Kjm]
(2)
[Kj|]
(2)
[Kij](3)
[Kii](l)
[Kmm](3)
[心]⑷
[Kmj]®
[Kim](4)
[Qi]⑷
[Kmm]
(2)
[Km|]
(2)
[Kii]
(2)[K»](3)[Kmm](4)
[Km|](4)
[Ku](4)
23456
i.3引入约束条件修正总刚度矩阵
整体刚度矩阵[K]求出后,结构的结点力{F}可表示为:
F[K]3
在无支杆的结点处,结点力就等于已知的结点载荷。
在有支杆的结点处,则求结点力时,还应把未知的支杆反力考虑在内。
如果用{P}表示结点载荷和支
杆反力组成的向量,则结点的平衡方程为
[K]3P
根据支承条件对平衡方程加以处理。
先考虑结点n有水平支杆的情况。
与
结点n水平方向对应的平衡方程是第2n-i个方程,
K2ni,iuiK2ni,2Vi■■-K2ni,2niUnK2ni,2nVn■■-En
根据支承情况,上式应换成,即在旳中,第2n-i行的对角线元素应
改为iK2该行全部非对角线元素应改为0。
在{P}中,第2n-i个元素应改为0。
此外,为了保持矩阵[K]的对称性,则第2n-i列全部非对角线元素也改为0。
同理,如果结点n有竖向支杆,则平衡方程的第2n个方程应改为,
为此,在矩阵[K]中,第2n行的对角线元素改为i,该行全部非对角线元素改为
00U2p(2
0
0
0
0
0
0
i
0
0
0
0
0Un
0
0
0
0
0
0
0
0
i
0
0
0
0vn
0
0
0
0
0
00
00
v2Py2
00
00
00
i.4载荷移置
将载荷移置到结点上,必须遵循静力等效的原则。
静力等效是指原载荷与结点载荷在任意虚位移上做的虚功相等。
在一定的位移模式下,移置结果是唯一的,且
{R}e[N]T{p}tdxdy
*eTe*Tt
由虚功相等可得:
{}{R}{}[N]{p}tdxdy
2、四边形、四节点四面体、八节点六面体单元分析
2.1四边形单元分析
四边单元可取矩形作为研究对象,矩形单元也是一种常用的单元,它采用了比常应变三角形单元次数更高的位移模式,因而可以更好地反映弹性体中的位移状态和应力状态。
]V4(V4)
V3(V3)
U4(U4)
41
3
U3(U3)
2b
o~•
■
1
2
U2(U2)
U1(U1)'
2aJ
在局部坐标系中,节点i的坐标是(i,i),
其值分别为土
1。
取位移模式:
V1(V1)
V2(V2)
x
矩形单元1234如图所示,其边长分别为2a和2b,两边分别平行于x、y轴。
若取该矩形的四个角点为节点,因每个节点位移有两个分量,所以矩形单元共有8个自由度。
采用三角形单元中的方法,同样可以完成对这种单元的力学特性分析。
然而,如果我们引入一个局部坐标系、,那么就可以推出比较简洁的
结果
用与三角形单元相似的方法建立形函数,则位移模式可写成:
式中:
NiNi1,110,i
ui
(i
1,2,3,4)
01
vi
由几何方程可以求得单元的应变
U
1
u
b」
xx
a
V
y
1
V
丄
Va-
gy
b
ab
gxy_u_v
1u
1
V
a_
u.b
y
x
b
a
可推出:
式中亠
B1B2B3B4
Bi
1
ab
aNi
1
4ab
由虎克定律我们可以得出用节点位移表示的单元应力,即
式中:
D
对于平面应力问题
有:
S1
S2S3
S4
Si
Bi
4ab1
bi1
bi1
ai1
若将单元刚度矩阵写成分块形式:
ai1
ai1
bi1
2
是,
素,
2.2
kn
k12
k13
k14
k21
k22
k23
k24
k31
k32
k33
k34
k41
k42
k43
k44
1
k
2
则其中的子矩阵可按下式进行计算:
BTDBjtdxdy
四边形单元的节点位移与单元节点力之间的关系仍为:
其中载荷列阵{R}e与上节中的(c)式相同,仍可按上式计算等效节点力
要注意的是,矩形单元有四个节点(1,2,3,4),所以{R}e具有8个元即:
ReU1V1U2V2U3V3U4V4T
rfn
keRe
。
但
四节点四面体单元分析
2.2.1单元划分及位移模式
如图1所示的四面体单元,单元结点的编码为i,j,m,n。
每个结点的位移具有三个分量u,v,w。
这样单元结点的位移列阵可表示成:
Tuvwuvwuvwuvw
iiiiiiiiiiii
单元的位移模式采用线性多项式:
3y4Z
7y8Z
11y12z
u12X
v56X
w910X
式中,为待定系数,由单元结点的位移和坐标决定。
将四个结点的坐标(Xi,
yi,z)、(为,yi,Zj)、(Xm,ym,zm)、(xn,yn,zn)和结点位移(Ui,Vi,Wi)、(uj,Vi,Wi)、
(Um,vm,Wm)、(Un,Vn,Wn)代入
(2)式可得12个联立方程,解方程组便可求出。
将这十二个系数回代到式,则得到由结点位移和形函数表示的单元内任一点的位
移表达式:
u
NiUi
NM
NmUm
NnUn
V
NiVi
Njvj
Nmvm
NnVn
w
Niwi
NjWj
Nmwm
NnWn
式中
1
Ni
6V
ai
biX
CiydiZ
1
N
j
6V
-aj
bjx
CjydjZ
1
N
m
6V
am
bmX
CmydmZ
1
N
n
6V
an
bnX
CnydnZ
Ni,Nj,Nm,Nn为四面体单元的形函数位移模式可以用矩阵形式表示:
式中,[I]为三阶单位阵,[N]为形函数矩阵。
上式即为单元结点位移和单元任意点位移之间的关系。
2.2.2单元应变和应力
知道单元内任意一点位移后,可利用几何方程确定单元内该点的应变。
将位移矩阵式代入空间问题几何方程得:
其中
上式表明几何矩阵[B]中的元素都是常量,因此单元中的应变也是常量。
也就是说,采用线性位移模式的四面体单元是常应变单元。
将上式代入物理方程,就得到单元的应力列阵:
Si
6A3
Arbi
Ac
di
V
A2Ci
Azbi
0
0
A2di
A2&
A2di
0
A2bi
DBi
bA1ciA1di
AbgAidi
式中:
[S]为四面体单元的应力矩阵,其分块形式为:
2.2.3单元刚度矩阵
对于四面体单元,利用虚功原理,采用类似平面问题的处理方法可以得到其单刚矩阵。
A3
V
brbsA2(CrCsdrds)
AcbsAzbrCs
AdbsAzbrds
A1brCsA2Crbs
CrCsA2(drdsbrbs)
A1drCsA2Crds
A1brbsAzdrds
AcdsAzdrCs
drdsA2(brbsCrCs)
BTDBdxdydz
其中:
[K]e为单元刚度矩阵
eT
KBD
Bdxdydz
T
B
DBV
写成分块形式为:
kii
kij
kim
kin
ekji
kjj
kjm
kjn
K
kmi
kmj
k
mm
k
mn
kni
knj
k
nm
knn
式中子矩阵[Krs]由下式计算:
T
krsBrDBsV
可以看出,单元刚度矩阵是由单元结点的坐标和单元材料的弹性常数所决定的,是一个常数矩阵。
2.3八节点六面体单元分析
2.3.1形函数与坐标变换形函数
简记为:
2.3.3单元刚度矩阵与等效节点载荷向量
坐标变换
2.3.2位移插值函数与几何矩阵
x
Ni
r,s,t
Xi
e
1
i1
8
u
x,y,x
N1
0
0
N2
0
0
n8
0
0
e
2
y
Ni
r,s,t
yi
V
x,y,x
0
N1
0
0
N2
0
0
N8
0
e
3
i1
8
w
x,y,x
0
0
N1
0
0
N2
0
0
N8
z
Ni
r,s,t
乙
e
24
8
i1
e
将上式中的
Ke
ve
x,y,z替换为
dv
Bdxdydz
r,S有t
Jdrdsdt
单兀刚度矩阵可以表示为:
N(ri,Sj,tk)Fb(ri,Sj,tk)J(ri,Sj,tk)hihjhk
Ke
写成高斯积分形式为:
n
Pbe
3、其他常用单元形函数、自由度
3.1轴对称单元
轴对称结构体可以看成由任意一个纵向剖面绕着纵轴旋转一周而形成。
此旋转轴即为对称轴,纵向剖面称为子午面,如图表示一圆柱体的子午面abed被分
割为若干个三角形单元,再经过绕对称轴旋转,圆柱体被离散成若干个三棱圆环单元,各单元之间用圆环形的铰链相连接。
对于轴对称问题,采用圆柱坐标较为方便。
以弹性体的对称轴为z轴,其约束及外载荷也都对称于z轴,因此弹性体内各点的各项应力分量、应变分量和位移分量都与环向坐标B无关,只是径向坐标r和轴向坐标z的函数。
也就是说,在任何一个过z轴的子午面上的位移、应变和应力的分布规律都相同。
因此轴对称问题可把三维问题简化为以(z,r)为
自变量的二维问题。
轴对称问题的物理方程可写为:
常采用的单元是三角形截面的整圆环。
在单元类型确定之后,单元剖分可以在子午面内进行,如图5-1表示的abede子午面被分割为若干个三角形,绕对称轴z旋转后即形成若干个三棱圆环单元。
相邻的单元由圆环形的铰链相连接。
单元的棱边都是圆,故称为结圆。
每个结圆与rz平面的交点称为结点。
这样,各单元在子午面rz平面上形成三角形网格,就如同平面问题中在xy平面上的网格一样。
采用位移法有限元分析,其基本未知量为结点位移。
单元的结点位移列阵如下:
eTTTTT
ijmUiWiUjWjUmWm
对于每一个环形单兀,性位移模式:
需要假定其位移模式。
仿照平面三角形单元,取线
uu(r,z)13Z
ww(r,z)46Z
类似于平面三角形单元的推导,即将单元的结点坐标及结点位移代入式中,可以解出六个待定系数。
再将这些待定系数回代到式中,就可以得到由结点位移和形函数所表示的单元内任一点的位移表达式:
有上面分析可知,轴对称单元自由度有六个
以上就是关于课程总结的全部内容。
通过这八周
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 有限元 学习心得