最新NSYS非线性求解的理论基础(共43张PPT课件).pptx

《最新NSYS非线性求解的理论基础(共43张PPT课件).pptx》由会员分享，可在线阅读，更多相关《最新NSYS非线性求解的理论基础(共43张PPT课件).pptx（43页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、ANSYS非线性求解的理论非线性求解的理论(lln)基础基础第一页，共四十三页。主要主要(zhyo)内容内容NROPT, , Option, -, Adptky 非线性方程组的解法CNVTOL, Lab, VALUE, TOLER, NORM, MINREF 收敛(shulin)准则AUTOTS, KeyARCLEN, Key, MAXARC, MINARC 过极值点的方法2第二页，共四十三页。非线性方程组的解法非线性方程组的解法(ji f)概述（概述（1）1.按照几何、材料、状态非线性理论建立的最后方程都是非线性的。2.在非线性分析中很难找到一种适合种类型非线性及各种非线性程度的解法。3.各

2、种解法均有各自的适用范围，选择不当可能引起(ynq)收敛困难甚至发散。4.因此一般通用的非线性有限元程序都提供几种解法以备用户选用。3第三页，共四十三页。非线性方程组的解法非线性方程组的解法(ji f)概述（概述（2）求解一个结构的平衡问题通常等于求解结构的总位能的驻值问题。结构总位能 ：公式（1）用虚功原理使公式（1）变分为(fn wi)零即得到有限元的平衡方程：公式（2）是的函数。 ( )TU uuF ( )( )0P uFK uuF ( )K u u4第四页，共四十三页。非线性方程组的解法非线性方程组的解法(ji f)增量法（增量法（1）顾名思义，增量法就是将荷载分成一系列的荷载增量，即

3、ANSYS中的荷载步或荷载子步。要点：在每一个荷载增量求解完成后，继续进行下一个荷载增量之前，调整刚度矩阵以反映(fnyng)结构刚度的变化。公式（3）公式（4） 1111iiiTiiiKuFuuu 5第五页，共四十三页。非线性方程组的解法非线性方程组的解法(ji f)增量法（增量法（2）l增量法的优点：可以追踪结构变形历程，这对于材料或几何非线性（特别是极限值屈曲分析）十分有用。l增量法的缺点：随着荷载步增量的增加而产生积累误差，导致荷载位移曲线飘移。l对飘移进行平衡(pnghng)修正，可以大大提高增量法的精度。现在应用最广的就是在每一级荷载增量上用Newton-Raphson或其变形的迭

4、代法。6第六页，共四十三页。荷载位移曲线荷载位移曲线(qxin)飘移示意图飘移示意图Fu真实解.误差.有限元解7第七页，共四十三页。非线性方程组的解法非线性方程组的解法(ji f)迭代法迭代法迭代法公式：公式（5）割线刚度(n d)法：，收敛性差，因此很少应用。切线刚度法：，令则得到Newton-Raphson迭代法。 11()iiiiiiuuKKuF11iiFKK11iiTKK18第八页，共四十三页。Full Newton-Raphson（F.N.R）Fu在一个荷载步或子步内每次迭代后重新(chngxn)形成一次iTK缺点：消耗机时！优点：适用于高度(god)非线性问题。9第九页，共四十三页

5、。Modified Newton-Raphson（M.N.R）Fu缺点：收敛性较差。优点(yudin)：消耗机时少。在一个(y )荷载步或子步内仅形成一次或每隔一定数目的迭代次数重新形成一次iTK10第十页，共四十三页。FuQuasi-Newton-Raphson（Q.N.R）0u1u2u3u0F1F2F3F0G1G2GiiiGuF 优缺点介于(ji y)前两者之间。11第十一页，共四十三页。非线性方程组的解法非线性方程组的解法(ji f)总结总结计算(j sun)效率：M.N.R高于Q.N.R 高于 F.N.R收敛性：F.N.R 高于Q.N.R 高于 M.N.R因此：当非线性程度不高（一般为

6、加载初期）时用M.N.R，当非线性程度较高（一般为加载后期）时用Q.N.R或F.N.R。12第十二页，共四十三页。非线性方程组的解法非线性方程组的解法(ji f)ANSYS中的应用中的应用Newton-Raphson 选项：选项：Main MenuSolutionAnalysis TypeAnalysis Options13第十三页，共四十三页。非线性方程组的解法非线性方程组的解法(ji f)ANSYS中的应用中的应用NROPT，AUTO：程序基于非线性程度选用。NROPT，FULL：如果自适应下降关闭（Adptky为OFF），每次迭代均使用切线刚度；如果自适应下降打开（Adptky为ON），

7、则迭代保持稳定时使用切线刚度，一旦程序在迭代中探测到发散倾向，将使用切线刚度和割线刚度的加权组合(zh)。NROPT，MODI：程序在每一个子步中修正切线刚度。自适应下降不可用。NROPT，INIT：程序在所有迭代中都使用初始刚度。自适应下降不可用。14第十四页，共四十三页。收敛准则收敛准则(zhnz)位移准则位移准则(zhnz) 11iDiiuuu D位移(wiy)收敛容差，一般取0.0010.005之间；iu某种范数，一般(ybn)为二范数，即：12211NiikkuuN位移收敛准则有时不可靠！15第十五页，共四十三页。位移收敛位移收敛(shulin)准则有时不可靠！准则有时不可靠！Fuu

8、F很小的位移变化(binhu)引起很大的力的变化(binhu)！16第十六页，共四十三页。收敛收敛(shulin)准则不平衡力准则准则不平衡力准则 iFFF iiFFP不平衡力 F外荷载(hzi)矢量 iP第 i 次迭代完成(wn chng)时与内力相平衡的节点力矢量，即： iiiPKu不平衡力收敛准则有时不可靠！17第十七页，共四十三页。不平衡力收敛准则有时(yush)不可靠！FuuF很小的力的变化(binhu)引起很大的位移变化(binhu)！18第十八页，共四十三页。收敛准则收敛准则(zhnz)能量准则能量准则(zhnz) 11TTiiEuFPuFP把每次迭代时内能(ni nn)的增量（

9、不平衡力在位移增量上做的功）与初始的内能(ni nn)比较：或者把上述(shngsh)内能增量与当前的总能量比较： TiiEiuFPU 11112TiiiiiUUuPP其中Ui为势能：需要更多的计算！19第十九页，共四十三页。收敛收敛(shulin)准则准则ANSYS中的应用（中的应用（1）Main MenuSolutionLoad Step OptsNonlinearConvergence Crit20第二十页，共四十三页。收敛收敛(shulin)准则准则ANSYS中的应用（中的应用（2）缺省的收敛准则：将不平衡力的范数与VALUETOLER进行比较，来对力或者力矩进行收敛检查。VALUE的

10、缺省值是所加荷载（或位移）的范数和MINREF（缺省为1.0）中的较大值。如果定义(dngy)了任何收敛准则，缺省准则将“失效”！一般总要使用力的收敛准则！21第二十一页，共四十三页。收敛准则收敛准则(zhnz)ANSYS中的应用（中的应用（3）用户收敛准则：CNVTOL, Lab, VALUE, TOLER, NORM, MINREFq使用严格的收敛准则将提高结果的精度，但以更多的平衡迭代为代价。q一般(ybn)情况下VALUE应使用默认值，而通过调整TOLER改变收敛准则。q应当确保MINREF=1.0在分范围内有意义。22第二十二页，共四十三页。过极值点的方法过极值点的方法(fngf)非

11、线性屈曲非线性屈曲非线性屈曲的主要难点：当荷载增至极值（临界）荷载Fcr时，结构的切线刚度趋于奇异(qy)，方程呈病态因而求解困难。因此，如何求得尽可能接近于真实的临界荷载就成为非线性屈曲研究的重点。荷载的两种类别：比例加载和非比例加载。23第二十三页，共四十三页。比例比例(bl)加载加载对大多数静机械荷载或稳态温度情况，在荷载施加的过程中它们的分布函数不变，而各点荷载值是按比例增减。根据这一特点(tdin)，可以把任一时刻的平衡状态选为参考状态，结构外荷载矢量为，其他时刻的荷载可以写为，求解极值荷载转化为求极值时的比例参数。下面仅讨论比例加载情况下过极值点（过屈曲）的方法：当前的主流方法是自

12、动步长法。rFttrFFcrFcr24第二十四页，共四十三页。自动自动(zdng)步长法概述（步长法概述（1）如前所述，非线性方程组的解法主要是增量加迭代法，即外荷载是逐级施加的，每施加一级荷载增量后进行平衡迭代直至解的结果满足允许容差。因此，荷载增量大小的控制对求解有很重要的意义：增量过大，解难于收敛甚至发散，这一点(y din)在接近极值点时特别明显；增量太小，意味着求解次数过多，使计算工作量增大。25第二十五页，共四十三页。自动自动(zdng)步长法概述（步长法概述（2）在一个问题的求解过程中，没有必要采用等荷载增量。对于非线性程度不高的阶段（一般是荷载作用的初始阶段），步长可选得大些；

13、对于非线性程度较高的阶段（一般是加载后期）步长需选得小些。当荷载接近极值(j zh)点时，很小的荷载增量都会引起很大位移，此时可不再直接施加荷载，而是改由控制其他参数使计算可以安全超过极值(j zh)点，人们仅给出一些控制参数，步长由算法控制，这类方法就称为自动步长法。26第二十六页，共四十三页。自动自动(zdng)步长法控制位移法（步长法控制位移法（1）控制(kngzh)位移法是将有限元中的荷载作为变量而位移作为独立量，即给出位移增量，反求荷载增量，就加载意义而言是一种自动步长法。下面是具体的推导过程：( )( )1112111( )2122222jjrjiriiiKKuFRKKuRF 式中

14、右上标(shn bio)为第 j 荷载步，右下标为第 i 次迭代。( )( )( )jjjRPF第 j 荷载步初始不平衡力。（6）27第二十七页，共四十三页。自动自动(zdng)步长法控制位移法（步长法控制位移法（2）假定位移(wiy)分量在屈曲时反应明显，将其选为人为的控制位移(wiy)分量，令，为给定值。由式（6）交换变量得：2u22uu2 u( )( )( )11111122( )( )( )221222jjrjjjrjiiiiuKFRKuRKFK（7）式（7）展开(zhn ki)： ( )( )( )( )11111122jjjjriiiiiKuFRKu ( )( )( )( )211

15、22222jjjjriiiiiKuFRKu （8）（9）28第二十八页，共四十三页。自动步长法控制自动步长法控制(kngzh)位移法（位移法（3）在一个迭代步中，（8）式为线性的，可以将其解分为(fn wi)两个部分：( )III111jiiiiuuu （10）将式（10）代入（8）中，分别(fnbi)得到：( )I1111jriiKuF ( )( )( )II1111122jjjiiiiKuRKu（11）（12）由式（11）、（12）分别可以求得、I1iuII1iu29第二十九页，共四十三页。自动自动(zdng)步长法控制位移法（步长法控制位移法（4）将由、表示(biosh)的代入式（9），

16、可以求得荷载增量因子：I1iuII1iu1iu( )( )( )2211222( )( )2211jjjIIiiijiijrIiiRKuKuFKu（13）( ) ji上述方法要求将刚度矩阵分解为、等矩阵，在接近极值点处不是很合适。因此可仿造M.N.R，在一系列迭代步中，一直采用初始刚度矩阵，将它分解为、，在迭代中不再(b zi)修改。 ( ) jiK( )11jiK( )12jiK( )21jiK( )22jiK ( )1jK( )111jK( )121jK( )211jK( )221jK30第三十页，共四十三页。自动自动(zdng)步长法控制位移法（步长法控制位移法（5）将解分为(fn wi

17、)两部分：iu ( )1jIriKuF ( )( )1jjIIiiKuR（15）（16）对其中(qzhng)的控制分量有：( )III2222jiiiiuuuu （17）由式（6）可得：( ) jIIIiiiiuuu （14）31第三十一页，共四十三页。自动自动(zdng)步长法控制位移法（步长法控制位移法（6）22( )1121IIjIuuu （18）从第二步迭代开始，因为为定值，故，荷载增量(zn lin)因子为： 2u20u2( )22,3,IIjiiIiuiu （19）第一次迭代可定出荷载增量(zn lin)因子：32第三十二页，共四十三页。自动步长法控制自动步长法控制(kngzh)位

18、移法（位移法（7）控制位移法的算法(sun f)步骤：1.求解式（15）得到，挑出(tio ch)指定位移分量。2.求解不平衡力及式（16）得到，挑出指定位移分量。3.如果是第一次迭代由式（18）求，如果是第二次以后的迭代，由式（19）求。4.由式（14）求。5.用迭代收敛准则检验是否得到收敛解，如果未得到，i=i+1，重复2到5步骤；如果收敛，j=j+1，转入下一位移增量步。Iiu2Iiu ( ) jiRIIiu2IIiu( )1j( ) jiiu33第三十三页，共四十三页。自动自动(zdng)步长法弧长法（步长法弧长法（1）每施加一个荷载步时直接控制走多长一段荷载位移(wiy)曲线（多维情

19、况为空间超曲线）。为此，需要对荷载增量比例参数加一个约束方程。对由 j 至 j+1 荷载步的第一次迭代，施加如下约束： 22111Tuuds （20）式中称为广义“弧长”，是由用户(yngh)控制的值。ds34第三十四页，共四十三页。自动自动(zdng)步长法弧长法（步长法弧长法（2）为保证迭代(di di)最近地走向收敛，迭代(di di)路径应沿 j 点切线的垂直面上，也即切线与该面上矢量 的标量积为零：tu0tu （21）式中 11jTietue （22） jTiiieuue （23）用矩阵(j zhn)表示式（21）为： 1102,3,Tiiuui （24）35第三十五页，共四十三页。

20、自动自动(zdng)步长法弧长法（步长法弧长法（3）与控制(kngzh)位移法一样，将未知矢量分解为两个部分：IIIiiiiuuu （25）式中及分别(fnbi)由式（15）和（16）求得。将式（25）代入式（24）可解出未知的荷载增量因子：IiuIIiuiu111TIIiiTIiuuuu （26）36第三十六页，共四十三页。自动自动(zdng)步长法弧长法（步长法弧长法（4）如果用M.N.R，在迭代(di di)中只需对式（15）求解一次，也即，1IIiuu （27）而在第一次迭代开始前，可以认为结构(jigu)处于平衡状态，因此，由式（16）可知，从而由式（25）可得： 0iR10IIu1

21、1Iiuu（28）37第三十七页，共四十三页。自动自动(zdng)步长法弧长法（步长法弧长法（5）将式（28）代入式（26）中，于是得到弧长法的最终(zu zhn)形式： 121111TIIiiTuuuu （29）弧长法与控制位移法的区别：求荷载增量因子的公式不同(b tn)！弧长法具有更好的适应性。i38第三十八页，共四十三页。自动自动(zdng)步长法弧长法（步长法弧长法（6）uds 1 2JJ+1tu弧长法原理图（一维）：39第三十九页，共四十三页。控制位移控制位移(wiy)法法ANSYS中的应用中的应用Main MenuSolutionLoad Step OptsTime/Freque

22、ncTime - Time Step40第四十页，共四十三页。弧长法弧长法ANSYS中的应用中的应用(yngyng)Main MenuSolutionLoad Step OptsNonlinearArc-Length Opts41第四十一页，共四十三页。The End谢 谢 大 家！42第四十二页，共四十三页。内容(nirng)总结ANSYS非线性求解的理论基础。ANSYS非线性求解的理论基础。 非线性方程组的解法。 收敛准则。 过极值点的方法。在一个(y )荷载步或子步内仅形成一次或每隔一定数目的迭代次数重新形成一次。不平衡力收敛准则有时不可靠。NonlinearConvergence Crit。将不平衡力的范数与VALUETOLER进行比较，来对力或者力矩进行收敛检查。如果定义了任何收敛准则，缺省准则将“失效”第四十三页，共四十三页。