非线性有限元的启蒙:基础概念与应用初探

在有限元分析中,我们经常会和非线性打交道,如材料非线性、几何非线性、边界非线性。非线性有限元一直是有限元中较为困难的一部分,在非线性有限元中我们经常碰到诸如牛顿-拉夫森迭代法,切线刚度阵等概念,今天贴主就简单的介绍一下非线性吧。

1、简单实例

首先看一个简单的弹簧杆件结构,如图所示,中间节点作用一个F的力,会产生一个位移v,

  非线性有限元初识的图1


由静力平衡关系可得到

非线性有限元初识的图2


该方程为典型的非线性方程,对于这个方程,如果给定一个位移v就能求得F,如下图所示,从图中曲线可以看到非线性的含义了。图中不同k对应的曲线,可以看到k比较小时,杆内力起主要作用,呈现出几何非线性,K较大时,弹簧起主要作用,呈现出弹簧的线弹性。

非线性有限元初识的图3


2、牛顿迭代法

但是在实际中,我们往往是不知道位移v的,而是知道F,那么给定一个F,怎么求v呢?这时候牛顿迭代法就要上场了。牛顿迭代法的思想是将非线性方程线性化,以线性方程的解逼近非线性方程的解,具体操作如下:

非线性有限元初识的图4

非线性有限元初识的图5

牛顿迭代法图形解释

对于非线性方程f(x)=的迭代解法有如下格式


4、非线性有限元迭代格式

虽然上文只是简单的一维问题,但是我们可以把它当做位移法有限元的原型,对于一般有限元,离散平衡方程一般具有如下形式:


 对于试探解、一般有


 该方程的求解有如下形式

(1)直接迭代法


非线性有限元初识的图10

非线性有限元初识的图11


非线性有限元初识的图13

直接迭代法中要求K矩阵为u的显式函数,只适用于和变形历史无关的非线性问题。该迭代法每次迭代都需要对新的求逆,计算量较大,于是有了如下改进的的常系数矩阵方法

非线性有限元初识的图17

(2)牛顿-辛普森迭代法 newton-Paphson method


运用泰勒展开:

(切线刚度阵)


非线性有限元初识的图22

同理,也可以得到修正的Newton-Paphson 方法


牛顿迭代法一般具有较好的收敛性,但是对于一些从小被分在二班的非线性同学,他也有很大的局限性

比如对于这个问题,牛顿只好呵呵了

非线性有限元初识的图24

对于下面问题,牛顿直接哭晕在厕所,当然这种问题只有等我们的arc-length兄来解决了。

非线性有限元初识的图25

再来看看我们上面的问题:

蓝色曲线为精确解,红色点点为固定载荷增量下求得的位移,k=1000时,牛顿迭代法能够很好地跟踪载荷位移路径,得到所有的位移响应。而当k=100时,曲线有下降段,此时牛顿迭代法就没法得到这个区域的位移响应了。  blob.png



免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删

QR Code
微信扫一扫,欢迎咨询~

联系我们
武汉格发信息技术有限公司
湖北省武汉市经开区科技园西路6号103孵化器
电话:155-2731-8020 座机:027-59821821
邮件:tanzw@gofarlic.com
Copyright © 2023 Gofarsoft Co.,Ltd. 保留所有权利
遇到许可问题?该如何解决!?
评估许可证实际采购量? 
不清楚软件许可证使用数据? 
收到软件厂商律师函!?  
想要少购买点许可证,节省费用? 
收到软件厂商侵权通告!?  
有正版license,但许可证不够用,需要新购? 
联系方式 155-2731-8020
预留信息,一起解决您的问题
* 姓名:
* 手机:

* 公司名称:

姓名不为空

手机不正确

公司不为空