课件ansys非线性求解的理论基础.pptVIP

ANSYS非线性求解的理论基础 * 自动步长法－－控制位移法（5） 将解　　分为两部分： （15） （16） 对其中的控制分量有： （17） 由式（6）可得： （14） ANSYS非线性求解的理论基础 * 自动步长法－－控制位移法（6） （18） 从第二步迭代开始，因为　　为定值，故　　　　，荷载增量因子为：　　 （19） 第一次迭代可定出荷载增量因子： ANSYS非线性求解的理论基础 * 自动步长法－－控制位移法（7） 控制位移法的算法步骤： 求解式（15）得到　　，挑出指定位移分量　　。 求解不平衡力　　及式（16）得到　　，挑出指定位移分量　　。 如果是第一次迭代由式（18）求　　，如果是第二次以后的迭代，由式（19）求　　。 由式（14）求　　。 用迭代收敛准则检验是否得到收敛解，如果未得到，i=i+1，重复2到5步骤；如果收敛，j=j+1，转入下一位移增量步。 ANSYS非线性求解的理论基础 * 自动步长法－－弧长法（1） 　　每施加一个荷载步时直接控制走多长一段荷载－位移曲线（多维情况为空间超曲线）。 　　为此，需要对荷载增量比例参数加一个约束方程。对由 j 至 j+1 荷载步的第一次迭代，施加如下约束： （20） 式中　　称为广义“弧长”，是由用户控制的值。 ANSYS非线性求解的理论基础 * 自动步长法－－弧长法（2） 　　为保证迭代最近地走向收敛，迭代路径应沿 j 点切线的垂直面上，也即切线　与该面上矢量　 的标量积为零： （21） 式中 （22） （23） 用矩阵表示式（21）为： （24） ANSYS非线性求解的理论基础 * 自动步长法－－弧长法（3） 　　与控制位移法一样，将未知矢量　　分解为两个部分： （25） 式中　　及　　分别由式（15）和（16）求得。 　　将式（25）代入式（24）可解出未知的荷载增量因子： （26） ANSYS非线性求解的理论基础 * 自动步长法－－弧长法（4） 　　如果用M.N.R，在迭代中只需对式（15）求解一次，也即， （27） 　　而在第一次迭代开始前，可以认为结构处于平衡状态，因此　　　，由式（16）可知　　　　，从而由式（25）可得：　　　 （28） ANSYS非线性求解的理论基础 * 自动步长法－－弧长法（5） 　　将式（28）代入式（26）中，于是得到弧长法的最终形式： （29） 弧长法与控制位移法的区别： 　　求荷载增量因子　　的公式不同！ 　　弧长法具有更好的适应性。 ANSYS非线性求解的理论基础 * 自动步长法－－弧长法（6） 弧长法原理图（一维）： ANSYS非线性求解的理论基础 * 控制位移法－－ANSYS中的应用 Main MenuSolutionLoad Step Opts 　　　　　Time/FrequencTime - Time Step ANSYS非线性求解的理论基础 * 弧长法－－ANSYS中的应用 Main MenuSolutionLoad Step Opts 　　　　　NonlinearArc-Length Opts ANSYS非线性求解的理论基础 * The End 谢 谢 大 家！ * ANSYS非线性求解的理论基础 ANSYS非线性求解的理论基础 * 主要内容 NROPT, , Option, --, Adptky —— 非线性方程组的解法 CNVTOL, Lab, VALUE, TOLER, NORM, MINREF 　—— 收敛准则 AUTOTS, Key 　ARCLEN, Key, MAXARC, MINARC 　—— 过极值点的方法 ANSYS非线性求解的理论基础 * 非线性方程组的解法－－概述（1） 按照几何、材料、状态非线性理论建立的最后方程都是非线性的。 在非线性分析中很难找到一种适合种类型非线性及各种非线性程度的解法。 各种解法均有各自的适用范围，选择不当可能引起收敛困难甚至发散。 因此一般通用的非线性有限元程序都提供几种解法以备用户选用。 ANSYS非线性求解的理论基础 * 非线性方程组的解法－－概述（2） 　　　求解一个结构的平衡问题通常等于求解结构的总位能的驻值问题。结构总位能Π ： 　　　　　　　　　　　　　　　　　　公式（1） 　　　用虚功原理使公式（1）变分为零即得到有限元的平衡方程： 　　　　　　　　　　　　　　　　　　公式（2） 　　　　　　是　　的函数。 ANSYS非线性求解的理论基础 * 非线性方程组的解法－－增量法（1） 　　　顾名思义，增量法就是将荷载分成一系列的荷载增量，即ANSYS中的荷载步或荷载子步。 　　　要点：在每一个荷载增量求解完成后，继续进行下一个荷载增量之前，调整刚度矩阵以反映结构刚度的变化。 　　　　　　　　　　　　　　　　　　公式（3）