没有合适的资源?快使用搜索试试~ 我知道了~
非线性有限元方法及实例分析,梁军,,对在地下工程稳定性分析中常用的非线性方程组的求解方法进行研究,讨论了非线性计算的迭代收敛准则,并利用非线性有限元方法分析
资源推荐
资源详情
资源评论
http://www.paper.edu.cn
非线性有限元方法及实例分析
梁军
河海大学水利水电工程学院,南京(210098)
摘 要:对在地下工程稳定性分析中常用的非线性方程组的求解方法进行研究,讨论了非线
性计算的迭代收敛准则,并利用非线性有限元方法分析了一个钢棒单轴拉伸的实例。
关键词:非线性有限元,方程组求解,实例分析
1 引 言
有限单元法已成为一种强有力的数值解法来解决工程中遇到的大量问题,其应用范围从
固体到流体,从静力到动力,从力学问题到非力学问题。有限元的线性分析已经设计工具被
广泛采用。但对于绝大多数水利工程中遇到的实际问题如地下洞室等,将其作为非线性问题
加以考虑更符合实际情况。根据产生非线性的原因,非线性问题主要有 3 种类型
[1]
:
1.材料非线性问题(简称材料非线性或物理非线性)
2.几何非线性问题
3.接触非线性问题(简称接触非线性或边界非线性)
2 非线性方程组的求解
在非线性力学中,无论是哪一类非线性问题,经过有限元离散后,它们都归结为求解一
个非线性代数方程组[2]:
()
()
()
0
0
0
21
212
211
=
……
=
=
nn
n
n
δδδ
ψ
δδδ
ψ
δ
δ
δ
ψ
Λ
Λ
Λ
(1.1)
其中
n
δ
δ
δ
,,,
21
Λ
是未知量,
n
ψ
ψ
ψ
,,,
21
Λ
是
n
δ
δ
δ
,,,
21
Λ
的非线性函数,引用矢量记
号
[]
T
n
δδδδ
Λ
21
=
(1.2)
[]
T
n
ψψψψ
Λ
21
=
(1.3)
上述方程组(1.1)可表示为
()
0=
δ
ψ
(1.4)
可以将它改写为
() () ()
0
=
−≡−≡ RKRF
δ
δ
δ
δ
ψ
(1.5)
其中
()
δ
K
是一个 的矩阵,其元素 是矢量的函数,
nn×
ij
k
R
为已知矢量。在位移有限
元中,
δ
代表未知的结点位移,
(
)
δ
F
是等效结点力,
R
为等效结点荷载,方程
()
0=
δ
ψ
表
示结点平衡方程。
在线弹性有限元中,线性方程组
-1-
http://www.paper.edu.cn
0=-RK
δ
(1.6)
可以毫无困难地求解,但对线性方程组
(
)
0
=
δ
ψ
则不行。一般来说,难以求得其精确
解,通常采用数值解法,把非线性问题转化为一系列线性问题。为了使这一系列线性解收敛
于非线性解,曾经有过许多方法,但这些解法都有一定的局限性。某解法对某一类非线性问
题有效,但对另一类问题可能不合适。因而,根据问题性质正确选用求解方法成为非线性有
限元的一个极重要的问题。常见的求解非线性方程组的数值方法有迭代法、增量法和混合法
[3][4][5]。
2.1. 迭代法
在每次迭代过程中都施加全部荷载,但逐步修改位移和应变,使之满足非线性的应力-
应变关系。迭代法还分为直接迭代法、Newton-Raphson 方法、修正的 Newton-Raphson 方法
和拟 Newton 法。
①直接迭代法
对非线性方程组
()
0=-RK
δ
δ
(1.7)
设其初始的近似解为 ,由此确定近似的
0
δδ
=
K
矩阵
(
)
00
δ
KK =
(1.8)
可得改进的近似解
()
RK
1
01
−
=
δ
(1.9)
重复这一过程,以第
i
次近似解求出第
1
+
i
次近似解的迭代公式为
(
)
ii
KK
δ
=
(1.10)
()
RK
ii
1
1
−
+
=
δ
(1.11)
直到 变得充分小,即近似收敛时,终止迭代。
iii
δδδ
−=Δ
+1
对于单变量问题,这一迭代过程是收敛的,单对于多自由度情况,由于未知量通过矩阵
K
耦合,迭代过程可能不收敛。
在迭代过程中,得到的近似解一般不会满足
(
)
0=-RK
δ
δ
即
(
)
(
)
0≠−≡ RK
iii
δδδ
ψ
(1.12)
()
δ
ψ
作为对平衡偏离的一种度量,称为失衡力。
② Newton-Raphson 方法
Newton-Raphson 方法是求解非线性方程组
-2-
剩余8页未读,继续阅读
资源评论
weixin_38675797
- 粉丝: 3
- 资源: 968
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
- (源码)基于Arduino和蓝牙技术的机器人控制系统.zip
- C#医院管理系统源码 门诊医生管理系统源码数据库 SQL2008源码类型 WinForm
- python项目9.一键生成QQ个人历史报告.zip
- python项目8.一键生成微信个人专属数据报告(了解你的微信社交历史).zip
- (源码)基于Java的账号管理系统.zip
- python项目6.爬取豆瓣排行榜电影数据(含GUI界面版).zip
- 微信小程序下拉选择组件
- pythonx项目5.爬取5K分辨率超清唯美壁纸.zip
- python项目4.每天不同时间段通过微信发消息提醒女友.zip
- python课程设计3.淘宝已买到的宝贝数据爬虫(已模拟登录).zip
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功