ANSYS非线性解的理论基础.ppt

ANSYS非线性求解的理论基础 ANSYS非线性求解的理论基础 主要内容 NROPT, , Option, --, Adptky —— 非线性方程组的解法 CNVTOL, Lab, VALUE, TOLER, NORM, MINREF 　—— 收敛准则 AUTOTS, Key 　ARCLEN, Key, MAXARC, MINARC 　—— 过极值点的方法 非线性方程组的解法－－概述（1） 按照几何、材料、状态非线性理论建立的最后方程都是非线性的。 在非线性分析中很难找到一种适合种类型非线性及各种非线性程度的解法。 各种解法均有各自的适用范围，选择不当可能引起收敛困难甚至发散。 因此一般通用的非线性有限元程序都提供几种解法以备用户选用。 非线性方程组的解法－－概述（2） 　　　求解一个结构的平衡问题通常等于求解结构的总位能的驻值问题。结构总位能Π ： 　　　　　　　　　　　　　　　　　　公式（1） 　　　用虚功原理使公式（1）变分为零即得到有限元的平衡方程： 　　　　　　　　　　　　　　　　　　公式（2） 　　　　　　是　　的函数。 非线性方程组的解法－－增量法（1） 　　　顾名思义，增量法就是将荷载分成一系列的荷载增量，即ANSYS中的荷载步或荷载子步。 　　　要点：在每一个荷载增量求解完成后，继续进行下一个荷载增量之前，调整刚度矩阵以反映结构刚度的变化。 　　　　　　　　　　　　　　　　　　公式（3） 　　　　　　　　　　　　　　　　　　公式（4） 非线性方程组的解法－－增量法（2） 增量法的优点：可以追踪结构变形历程，这对于材料或几何非线性（特别是极限值屈曲分析）十分有用。 增量法的缺点：随着荷载步增量的增加而产生积累误差，导致荷载－位移曲线飘移。 对飘移进行平衡修正，可以大大提高增量法的精度。现在应用最广的就是在每一级荷载增量上用Newton-Raphson或其变形的迭代法。 荷载－位移曲线飘移示意图 非线性方程组的解法－－迭代法 　迭代法公式： 　　　　　　　　　　　　　　　　　　公式（5） 　割线刚度法：　　　　　　，收敛性差，因此很少应用。 　切线刚度法：　　　　　　，令　　　则得到Newton-Raphson迭代法。 Full Newton-Raphson（F.N.R） Modified Newton-Raphson（M.N.R） Quasi-Newton-Raphson（Q.N.R） 非线性方程组的解法－－总结 计算效率： M.N.R　高于　Q.N.R　 高于　 F.N.R 收敛性： F.N.R　 高于　Q.N.R　 高于　 M.N.R 因此： 　当非线性程度不高（一般为加载初期）时用M.N.R，当非线性程度较高（一般为加载后期）时用Q.N.R或F.N.R。 非线性方程组的解法－－ANSYS中的应用 非线性方程组的解法－－ANSYS中的应用 收敛准则－－位移准则 位移收敛准则有时不可靠！ 收敛准则－－不平衡力准则 不平衡力收敛准则有时不可靠！ 收敛准则－－能量准则 收敛准则－－ANSYS中的应用（1） 收敛准则－－ANSYS中的应用（2） 收敛准则－－ANSYS中的应用（3） 过极值点的方法－－非线性屈曲 比例加载 自动步长法－－概述（1） 自动步长法－－概述（2） 自动步长法－－控制位移法（1） 自动步长法－－控制位移法（2） 自动步长法－－控制位移法（3） 自动步长法－－控制位移法（4） 自动步长法－－控制位移法（5） 自动步长法－－控制位移法（6） 自动步长法－－控制位移法（7） 自动步长法－－弧长法（1） 自动步长法－－弧长法（2） 自动步长法－－弧长法（3） 自动步长法－－弧长法（4） 自动步长法－－弧长法（5） 自动步长法－－弧长法（6） 控制位移法－－ANSYS中的应用 弧长法－－ANSYS中的应用 The End 将解　　分为两部分： （15） （16） 对其中的控制分量有： （17） 由式（6）可得： （14） （18） 从第二步迭代开始，因为　　为定值，故　　　　，荷载增量因子为：　　 （19） 第一次迭代可定出荷载增量因子： 控制位移法的算法步骤： 求解式（15）得到　　，挑出指定位移分量　　。 求解不平衡力　　及式（16）得到　　，挑出指定位移分量　　。 如果是第一次迭代由式（18）求　　，如果是第二次以后的迭代，由式（19）求　　。 由式（14）求　　。 用迭代收敛准则检验是否得到收敛解，如果未得到，i=i+1，重复2到5步骤；如果收敛，j=j+1，转入下一位移增量步。 　　每施加一个荷载步时直接控制走多长一段荷载－位移曲线（多维情况为空间超曲线）。 　　为此，需要对荷载增量比例参数加一个约束方程。对由 j 至 j+1 荷载步的第一次迭代，施加如下约束： （20） 式中