计算固体力学(有限元以及无网格方法)全套教学【121P】PPT课件
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
vi
i(xi , yi ) u i
um vj
uj j(xj , yj )
x O
三角形单元
将位移试函数代入上式，并求偏导数，得
xxyy222111 (((bcciiiuuviii
bjuj cjvj cjuj
bmum) cmvm) cmum)(bivi
bjvj
bmvm)
第二章 平面弹性力学的有限元法
反映了单元的位移形态，称为形函数
vm
m (xm, ym)
vi
i(xi , yi ) u i
um vj
uj j(xj , yj )
x
三角形单元
同理有 vN iv i N jvj N m v m N kv k
则位移向量可表示为
i,j,m
{ } e 单元节点位移向量
ui
vi
{f
}
u v
Ni
0
0 Ni
求
L(u)0
解 域
u aiui
离 散
i
L'(ui) 0
AXB
各种数值方法
ui u(xi)离散节点的变量值
第一章 科学和工程中的数值方法
1.3 几个简单示例
(a) 开孔板力学模型
(b) 力学模型离散化
平面问题有限元法
第一章 科学和工程中的数值方法
BEM的变形
起重机吊钩
FEM的变形
第一章 科学和工程中的数值方法
2.2 三角形常应变单元
y
3 单元中的应变和应力
{}[B]{}e
由于[B]是常量，单元内各点应变分
量也都是常量，这是由于采用了线性位移 O 函数的缘故，这种单元称为三角形常应变 单元。
vm
m (xm, ym)
vi
i(xi , yi ) u i
E
1
1
0 0
1 2 0
0
1
2
第二章 平面弹性力学的有限元法
2.1 平面弹性力学的基本方程
2 平面应变问题
x 0 , x z0 , y z0
但一般情况下 z 0
平面应变问题的弹性矩阵只需将上页中的
E
换成 1
E
2
换成
1
即可。
1
D
E(1) (1)(12)
1
1
1
0
0
0
0
12
2(1)
第二章 平面弹性力学的有限元法
2.1 平面弹性力学的基本方程
无论是平面应力问题还是平面应变问题的应力
与 应变 之间的关系均为：
D 0 ，其中：
x
y
T xyxy源自 xy式中 0 为初应变。
第二章 平面弹性力学的有限元法
2.2 三角形常应变单元 1 结构的离散化
um vj
uj j(xj , yj )
x O
三角形单元
u (x,y)1 2x 3y v(x,y) 4 5x 6y
u iui
i
第二章 平面弹性力学的有限元法
2.2 三角形常应变单元
2 位移试函数
由于位移函数适用于单元中的任意 一点，所以代入3个节点的坐标后，得 出节点处位移函数为
ui12xi3yi uj 12xj3yj um12xm3ym
联立求解，得
121(aiui ajuj amum) 221(biui bjuj bmum) 321(ciui cjuj cmum)
y
vm
m (xm, ym)
vi
i(xi , yi ) u i
um vj
uj j(xj , yj )
x O
三角形单元
1 11
2
xi xj
1 xm
yi yj
12(bicj bjci)
岩体裂隙扩展
第二章 平面弹性力学的有限元法
2.1 平面弹性力学的基本方程 2.2 三角形常应变单元 2.3 单元的总势能 2.4 物体的总势能 2.5 有限元方程组 2.6 面积坐标 2.7 节点荷载向量的计算 2.8 结构总刚度矩阵的性质 2.9 有限元程序设计 2.10 有限元解的收敛性 2.11 矩形单元 2.12 6节点三角形单元 2.13 等参元
计算固体力学
整体 概述
一 请在这里输入您的主要叙述内容
二
请在这里输入您的主要 叙述内容
三 请在这里输入您的主要叙述内容
说
目录
第一章 科学和工程中的数值方法 第二章 平面弹性力学的有限元法 第三章 三维弹性力学的有限元法 第四章 无网格方法
第一章 科学和工程中的数值方法
1.1 科学和工程问题的研究方法 1.2 科学和工程中的数值方法 1.3 几个简单示例
加权余量法（加权残数法）
基础方法
最小二乘法、Galerkin方法
有限差分法
辅助方法
有限元法 边界元法
目前常用的计算固体力学方法
(在科学和工程各个领域得到广泛应用)
无网格方法
目前研究的热点之一
第一章 科学和工程中的数值方法
1.2 科学和工程中的数值方法
科学和工程中的数值方法的一般思路： 连续－离散
求解方程离散
Nj 0
0 Nj
Nm 0
0 Nm
uvjj
[IN i IN j IN m]{}e
um
vm
（位移试函数）
第二章 平面弹性力学的有限元法
2.2 三角形常应变单元 3 单元中的应变和应力
由几何方程知
{ }
x y
xy
x
0
y
0
y
u v
x
y
vm
m (xm, ym)
第一章 科学和工程中的数值方法
1.1 科学和工程问题的研究方法
科学研究的三大方法： 实验研究 (真实性)
实验条件，人力、物力和财力
理论分析 (严密性)
基本基于简单问题，难以适应复杂问题
数值方法 (近似性、虚拟性)
应用广泛，与实验和理论分析结合具有广阔前景
第一章 科学和工程中的数值方法
1.2 科学和工程中的数值方法
第二章 平面弹性力学的有限元法
2.1 平面弹性力学的基本方程
O
平面问题应力状态
平衡方程 几何方程 物理方程 边界条件
最小势能原理
第二章 平面弹性力学的有限元法
2.1 平面弹性力学的基本方程 1 平面应力问题
z 0 , x zz x0 , y zz y0 ,
z 0
平面应力问题的弹性矩阵为：
D
(a) 均匀受力板力学模型
(b) 力学模型离散化
平面问题有限单元法的计算力学模型
第二章 平面弹性力学的有限元法
2.2 三角形常应变单元
y
vm
2 位移试函数
在以位移为基本未知量的有 限元法中，需构造出单元的位移 插值函数，即位移试函数。
采用线性插值，位移试函数 可表示为
m (xm, ym)
vi
i(xi , yi ) u i
ym
（ ijm的面积，ijm逆时针0 ）
ai xjymxmyj bi yj ym ci xj xm
第二章 平面弹性力学的有限元法
2.2 三角形常应变单元
y
2 位移试函数
u N iu i N ju j N m u m N ku k
i,j,m
N i2 1 (aibixciy), (i,j,m )
O
i(xi , yi ) u i
um vj
uj j(xj , yj )
x O
三角形单元
将位移试函数代入上式，并求偏导数，得
xxyy222111 (((bcciiiuuviii
bjuj cjvj cjuj
bmum) cmvm) cmum)(bivi
bjvj
bmvm)
第二章 平面弹性力学的有限元法
反映了单元的位移形态，称为形函数
vm
m (xm, ym)
vi
i(xi , yi ) u i
um vj
uj j(xj , yj )
x
三角形单元
同理有 vN iv i N jvj N m v m N kv k
则位移向量可表示为
i,j,m
{ } e 单元节点位移向量
ui
vi
{f
}
u v
Ni
0
0 Ni
求
L(u)0
解 域
u aiui
离 散
i
L'(ui) 0
AXB
各种数值方法
ui u(xi)离散节点的变量值
第一章 科学和工程中的数值方法
1.3 几个简单示例
(a) 开孔板力学模型
(b) 力学模型离散化
平面问题有限元法
第一章 科学和工程中的数值方法
BEM的变形
起重机吊钩
FEM的变形
第一章 科学和工程中的数值方法
2.2 三角形常应变单元
y
3 单元中的应变和应力
{}[B]{}e
由于[B]是常量，单元内各点应变分
量也都是常量，这是由于采用了线性位移 O 函数的缘故，这种单元称为三角形常应变 单元。
vm
m (xm, ym)
vi
i(xi , yi ) u i
E
1
1
0 0
1 2 0
0
1
2
第二章 平面弹性力学的有限元法
2.1 平面弹性力学的基本方程
2 平面应变问题
x 0 , x z0 , y z0
但一般情况下 z 0
平面应变问题的弹性矩阵只需将上页中的
E
换成 1
E
2
换成
1
即可。
1
D
E(1) (1)(12)
1
1
1
0
0
0
0
12
2(1)
第二章 平面弹性力学的有限元法
2.1 平面弹性力学的基本方程
无论是平面应力问题还是平面应变问题的应力
与 应变 之间的关系均为：
D 0 ，其中：
x
y
T xyxy源自 xy式中 0 为初应变。
第二章 平面弹性力学的有限元法
2.2 三角形常应变单元 1 结构的离散化
um vj
uj j(xj , yj )
x O
三角形单元
u (x,y)1 2x 3y v(x,y) 4 5x 6y
u iui
i
第二章 平面弹性力学的有限元法
2.2 三角形常应变单元
2 位移试函数
由于位移函数适用于单元中的任意 一点，所以代入3个节点的坐标后，得 出节点处位移函数为
ui12xi3yi uj 12xj3yj um12xm3ym
联立求解，得
121(aiui ajuj amum) 221(biui bjuj bmum) 321(ciui cjuj cmum)
y
vm
m (xm, ym)
vi
i(xi , yi ) u i
um vj
uj j(xj , yj )
x O
三角形单元
1 11
2
xi xj
1 xm
yi yj
12(bicj bjci)
岩体裂隙扩展
第二章 平面弹性力学的有限元法
2.1 平面弹性力学的基本方程 2.2 三角形常应变单元 2.3 单元的总势能 2.4 物体的总势能 2.5 有限元方程组 2.6 面积坐标 2.7 节点荷载向量的计算 2.8 结构总刚度矩阵的性质 2.9 有限元程序设计 2.10 有限元解的收敛性 2.11 矩形单元 2.12 6节点三角形单元 2.13 等参元
计算固体力学
整体 概述
一 请在这里输入您的主要叙述内容
二
请在这里输入您的主要 叙述内容
三 请在这里输入您的主要叙述内容
说
目录
第一章 科学和工程中的数值方法 第二章 平面弹性力学的有限元法 第三章 三维弹性力学的有限元法 第四章 无网格方法
第一章 科学和工程中的数值方法
1.1 科学和工程问题的研究方法 1.2 科学和工程中的数值方法 1.3 几个简单示例
加权余量法（加权残数法）
基础方法
最小二乘法、Galerkin方法
有限差分法
辅助方法
有限元法 边界元法
目前常用的计算固体力学方法
(在科学和工程各个领域得到广泛应用)
无网格方法
目前研究的热点之一
第一章 科学和工程中的数值方法
1.2 科学和工程中的数值方法
科学和工程中的数值方法的一般思路： 连续－离散
求解方程离散
Nj 0
0 Nj
Nm 0
0 Nm
uvjj
[IN i IN j IN m]{}e
um
vm
（位移试函数）
第二章 平面弹性力学的有限元法
2.2 三角形常应变单元 3 单元中的应变和应力
由几何方程知
{ }
x y
xy
x
0
y
0
y
u v
x
y
vm
m (xm, ym)
第一章 科学和工程中的数值方法
1.1 科学和工程问题的研究方法
科学研究的三大方法： 实验研究 (真实性)
实验条件，人力、物力和财力
理论分析 (严密性)
基本基于简单问题，难以适应复杂问题
数值方法 (近似性、虚拟性)
应用广泛，与实验和理论分析结合具有广阔前景
第一章 科学和工程中的数值方法
1.2 科学和工程中的数值方法
第二章 平面弹性力学的有限元法
2.1 平面弹性力学的基本方程
O
平面问题应力状态
平衡方程 几何方程 物理方程 边界条件
最小势能原理
第二章 平面弹性力学的有限元法
2.1 平面弹性力学的基本方程 1 平面应力问题
z 0 , x zz x0 , y zz y0 ,
z 0
平面应力问题的弹性矩阵为：
D
(a) 均匀受力板力学模型
(b) 力学模型离散化
平面问题有限单元法的计算力学模型
第二章 平面弹性力学的有限元法
2.2 三角形常应变单元
y
vm
2 位移试函数
在以位移为基本未知量的有 限元法中，需构造出单元的位移 插值函数，即位移试函数。
采用线性插值，位移试函数 可表示为
m (xm, ym)
vi
i(xi , yi ) u i
ym
（ ijm的面积，ijm逆时针0 ）
ai xjymxmyj bi yj ym ci xj xm
第二章 平面弹性力学的有限元法
2.2 三角形常应变单元
y
2 位移试函数
u N iu i N ju j N m u m N ku k
i,j,m
N i2 1 (aibixciy), (i,j,m )
O