非线性分析培训课件.ppt

4、材料非线性-非线性计算方法 4、材料非线性-非线性计算方法 ANSYS 使用Newton-Raphson平衡迭代法 克服了增量求解的问题。 在每个载荷增量步结束时，平衡迭代驱使解回到平衡状态。 F u 位移 载荷 1 2 3 4 一个载荷增量中全 Newton-Raphson 迭代求解。（四个迭代步如图所示） 4、材料非线性-非线性计算方法 Newton-Raphson 法迭代求解使用下列方程： [KT]{?u} = {Fa} - {Fnr} 这里: [KT] = 切向刚度矩阵 {?u} = 位移增量 {Fa} = 施加的载荷矢量 {Fnr} = 内力矢量 目标是迭代至收敛 (后面定义)。 Fa u 1 2 3 4 [KT] 4、材料非线性-非线性计算方法 Newton-Raphson法是ANSYS用于求解非线性方程组的一种数值方法 。 Newton-Raphson法基于增量加载与迭代，使每个载荷增量步达到平衡。 Newton-Raphson 法的优点是对于一致的切向刚度矩阵有二次收敛速度。 也就是每一迭代步的求解误差与前一步误差的平方成正比。 修正和初始刚度N-R求解过程收敛速度比完全N-R法慢，但是这两种方法使用较少进行重新求解切线刚度矩阵及其逆矩阵。 4、材料非线性-非线性计算方法 由方程（11-129）和（11-130）可知，在每次迭代过程中刚度矩阵都被更新。这种求解方法就是完全Newton-Raphson法。 此外，用户可以使用修正的Newton-Raphson方法（NROPT,MODI）来减少刚度矩阵的更新次数。他的特别之处在于，对于静态或瞬态分析，修正的Newton-Raphson方法仅仅在每一个子步的第一次或第二次迭代更新刚度矩阵。如图11-33所示，初始刚度Newton-Raphson方法对切线刚度矩阵不做任何更新。 4、材料非线性-非线性计算方法 Newton-Raphson 法需要一个收敛的度量以决定何时结束迭代。给定外部载荷（Fa），内部载荷（ Fnr ）（由单元应力产生并作用于节点），在一个体中，外部载荷必须与内力相平衡。 Fa - Fnr = 0 收敛是平衡的度量。 4、材料非线性-非线性计算方法 Newton-Raphson 迭代过程如下所示。基于 u0 时的结构构形，计算出切向刚度KT，基于?F 计算出的位移增量是?u ，结构构形更新为 u1。 Fa u1 Fnr ?u ?F u0 位移 载荷 R KT 在更新的构形中计算出内力（单元力） 。 迭代中的Newton-Raphson 不平衡量是: R = Fa - Fnr 4、材料非线性-非线性计算方法 Newton-Raphson不平衡量 (Fa - Fnr) 实际上从未真正等于零。当不平衡量小到误差允许范围内时，可中止Newton-Raphson 迭代，得到平衡解。 在数学上，当不平衡量的范数||{Fa} - {Fnr}||小于指定容限乘以参考力的值时就认为得到收敛。 4、材料非线性-非线性计算方法 4、材料非线性-非线性计算方法 4、材料非线性-非线性计算方法 ANSYS 缺省的收敛判据是力 / 力矩和位移 / 旋转增量。 对于力 / 力矩缺省的容限是0.5%，对于位移 / 旋转增量的容限是 5% 。 经验表明这些容限对于大多数问题具有足够的精确度。缺省的设置对于广泛的工程问题既不“太紧”也不“太松”。 4、材料非线性-非线性计算方法 力收敛判据提供了一个收敛的绝对度量，因为它可直接度量内部力与外部力间的平衡。 基于检查的位移判据只应作为力收敛判据的辅助手段使用。 只依据位移判断收敛在一些情况下将导致错误的结果。 4、材料非线性-非线性计算方法 4、材料非线性-非线性计算方法 非线性求解可按下列三个层次组织： 载荷步 载荷步是顶层，求解选项，载荷与边界条件都施加于某个载荷步内。 子步 子步是载荷步中的载荷增量。子步用于逐步施加载荷。 平衡迭代步 平衡迭代步是ANSYS为得到给定子步（载荷增量）的收敛解而采用的方法。 4、材料非线性-非线性计算方法 在每一增量载荷步中完成平衡迭代步。 载荷步一中有两个子步，载荷步二中有三个子步。 每个载荷步及子步都与 “ 时间 ”相关联。 “时间” 载荷 载荷步 2 载荷 1 子步 两个载荷步的求解 4、材料非线性-非线性计算方法 每个载荷步与子步都与 “ 时间 ”相关联。 子步 也叫时间步。 在率相关分析（蠕变，粘塑性）与瞬态分析中，“ 时间 ”代表真实的时间。 对于率无关的静态分析，“ 时间