本帖最后由 seerainbow 于 2011-8-24 16:13 编辑
已经知道ANSYS利用Newton-Raphson 法利用方程:
[KT]{Du} = {F} - {Fnr} 迭代求出收敛解
程序反复求解方程直到残差{R}(不平衡力),即{R}={F} - {Fnr}足够小。允许的最大残差值称为收敛准则(convergence norm)。
收敛准则数学上表达为:若||{R}|| < (εR *Rref),则:求解收敛。
式中
||{R}|| 是残差{R}的矢量范数 (范数是一个将矢量转换为单一标量值的算子)。
残差的 L1 范数:||{R}||1 =∑|Ri|
残差的 L2 (SRSS) 范数:||{R}||2 = (∑R2i)1/2,即F L2。
残差的无穷范数: ||{R}||¥ = max(|Ri|)
(εR *Rref) 是力收敛准则,即F CRIT= εR *Rref。
εR是容差因子。Rref 是参考力值,Rref 可以是所有施加的力和反作用力的范数 (自动地将准则调整为载荷的量级)。
收敛容差准则命令:
CNVTOL, Lab, VALUE, TOLER, NORM, MINREF
Lab是检查收敛的项目如力,位移等。
Value是残差||{R}||的典型参考值Rref。最小参考值为MINREF。MINREF默认为0.01。
TOLER是容差因子εR。默认为0.005。
NORM是残差的矢量范数选择。L1,L2,或者无穷范数。
缺省准则对大部分的工程问题适用: ||{R}||2 < (0.5% * ||{F}||2) 。
特殊情况下可以改变准则。可以加紧或放松收敛准则。加紧准则能提高精度, 但收敛困难。也可选择其它项目检查收敛性。可用力, 力矩, 位移及转动准则。另外, 可以改变用来度量收敛项目的范数(L1, L2或无限范数)。
改变任何收敛准则,程序将删除所有缺省准则! 第一次改变准则时,出现缺省值的菜单列表。注:
动作按钮为Replace。记住,若改变任何收敛准则,则程序删除所有缺省准则! 因此如果定义了位移收敛准则,应再定义力收敛准则,力准则才有效。 为什么必须再建立力收敛准则? 因为基于位移的检查是收敛的相对度量,它应仅作为基于力收敛的补充。基于力收敛提供收敛的绝对度量,因为它是内力与外力之间平衡的度量。仅依据位移收敛有时会导致错误的结果。 Solution>load step opts>Grph solu Track打开收敛历史曲线图 命令:/GST,ON/OFF 迭代历史曲线图中CRIT表示力收敛准则值,L2表示力残差的L2范数。 |