首页/文章/ 详情

有限元基础知识:线性化

2小时前浏览41

接着上一次的偏导与方向导数,微分与变分,有限元基础知识:偏导、方向导数、微分与变分

这次讲一下线性化,线性化顾名思义其实其实就是一种”以直代曲“的方法,英文叫Linearization,对于非线性问题的求解极为重要。

大家上大学学高等数学的时候都学过泰勒展开,当时可能就在想“这玩意有啥用?”,那其实线性化就是一个对于复杂非线性问题的泰勒级数展开,只是这里只取其一次项进行线性近似,忽略其他高阶项。

首先我们来说一下牛顿法的基本流程,牛顿法的基本思想就是根据当前的状态求出切线矩阵,并求出残余向量,进行位移增量的求解,即:

 

这本质上就是一种“以直代曲”的思想去逐步迭代接近真实值,如下图所示:

那么何为切线矩阵呢?其实这里我们已经隐含了线性化的概念,现在我们看如下的一个非线性方程

 

这里在非线性分析中    是    的函数,    也是     的函数(如几何非线性的情况下),这个式子非常难以求解,那么就需要我们对其进行线性化,具体做法为:

 

那这就是一个典型对于复杂     在     附近的泰勒展开,但我们忽略高阶项,以上公式则变为

 

这里为了公式简便,先暂时不考虑    随着     变化的情况,那么进一步展开上边的式子我们可以得到

 

而我们可以认为    , 可以认为是在初始位移向量对应的内力,那么通过移项,我们就得到了最初的Newton法的基本公式,所以这也就是经常所说的牛顿法和泰勒级数展开与线性化都有很强的关系,也就是牛顿法所谓2阶收敛性的原因,因为他近似只取到泰勒级数的一阶,误差是    

这样对于材料非线性、接触非线性、几何非线性的问题我们就通过以下公式求解其切线矩阵

 

现在可以回想下我所常用的材料:Mises塑性2 中进行的那些应力增量的求解大家一定特别理解,因为本质上内力和应力有如下的关系

 

那么对于求解切线刚度矩阵的时候对于内力求导,在    不变的情况下则可以转化为

 

最后大家需要注意的是,这里的线性化是非线性方程求解的线性化,与应力线性化是完全不同的概念。



来源:大狗子说数值模拟
非线性材料
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-01-22
最近编辑:2小时前
大狗子说数值模拟
博士 传播国际一流的数值模拟算法
获赞 4粉丝 5文章 41课程 0
点赞
收藏
作者推荐

我所常用的材料:Mises塑性2

接上回:我所常用的材料:Mises塑性上次已经说到,应变变化率与应力变化率之间的关系由于这里与当前的应力状态相关,硬化系数也与当前的应力状态相关,总的来说我们需要在单元积分点处根据其应变增量得到其真实的应力状态,通常采用迭代的方法确定真实的应变增量,具体可以分为两个类型的方法:向后Euler法(需要积分点级别的Newton迭代)向前Euler法一般来说对于显式分析我们可能采用向前欧拉法,根据当前的应力状态,应变增量逐步计算出来应力增量,而不进行残差方面的校核(或只在最后进行一个简单的修正)。而一般对于隐式分析,我们则一般采用向后Euler方法,这里也重点介绍这个。现在我们认为有如下的公式:其中代表着第步的应变、应力状态,,为第步骤的试应力,也就是说未必是真实的,需要进行修正,后边的则就是对其的修正,具体可以结合上图理解为,我先试着用当前的应变与上一步的塑性应变进行一个试应力的计算,不行我再沿着n方向修正回来,这里的n为一个2阶张量,在上一期中又说,对于Mises塑性就为f的梯度。那么现在结合着最开始的塑性公式,在仅考虑Isotropichardening的情况下(强化因子可以由应力应变状态之后计算出来,而Kinematic则不行,这里先以Isotropichardening为例,比较简单),那么我们有7个未知数,7条公式(应力有6个公式,有1个),且是个非线性方程组(与当前的应力状态相关),那么我们可以将上述方程联立进行Netwon-Raphson求解在牛顿法中经典的思路就是对函数的变量求偏导,乘上变量的增量,去逐步的减少残差,以直代曲,所以我们有:最后得到这样通过逐步的迭代,消除残差,当的时候,最终我们就得到了步的应力、应变、弹性应变、塑性应变状态。而这样计算出来的应力应变才是我们有限元计算中塑性材料应力应变的由来,而并不是给定一条材料曲线直接去根据材料曲线查表。所以大家可以看到对于Mises塑性材料的非线性分析,这里是有两个Newton-Raphson迭代的,外层的迭代进行整体级别的迭代,计算,而内层的迭代发生在每个单元积分点级别,通过迭代确定每个积分点处真实的应变应力状态,进而基于这个真实应力应变状态才能确定全局级别的切线刚度矩阵。而对于后续的切线刚度矩阵有consistent与continumm之分,下次再说说这个。来源:大狗子说数值模拟

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈