ADAMS求解方案.pptxVIP

ADAMS求解方案;;采用修正的Newton－Raphson迭代算法求解 对刚性系统，采用变系数BDF（Backwards Differentiation Formula）刚性积分程序，它是自动变阶、变步长的预估矫正法，在积分的每一步采用了修正的Newton－Raphson迭代算法；对调频系统采用坐标分配法和ABAM（Adams-Bashforth-Adams-Moulton）方法。与之对应，ADAMS/Solver中包含了三个功能强大的求解器。 ;积分器(Integrator).ADAMS中常用的刚性积分器有三种：GSTIFF、WSTIFF和DBF等3种不同的刚性积分器。计算的效率和稳定性也不相同，计算的稳定性为DBFWSTIFF GSTIFF,计算效率为GSTIFFWSTIFFDBF. GSTIFF是最常用的积分器，也是系统默认的积分器。DBF积分器计算时间最长，但有时可求解GSTIF求解失败的问题。 积分格式（Formulation）。ADAMS提供了三种积分格式：I3，SI2和SI1.I3格式求解速度快，但精度差；SI2格式可以避免Jacobian矩阵的病态，而且考虑了约束方程，虽然求解时间变长，但求解精度高；SI1求解精度比SI2还高，但计算量太大，一般不采用。 ;?;ADAMS中用的最多的积分器是GSTIFF。它运用至少两种求解方程，一个是I3方程，另一个是SI2方程（将ADAMS/Solver选为Fortran时还有一个I1）。I3 GSTIFF方程已经有了二十年历史。它已经发展的颇为完善了。但是，也有一些模型用WSTIFF积分器可以求解而不能用GSTFF积分器求解。所以，在用GSTIFF积分器不能求解模型时，可以尝试用WSTIFF积分器，或者用C++求解器中HHT积分器等新型积分器求解。;拿GSTIFF SI2积分器和GSTIFF I3积分器来对比，当提及中断操作时，I3积分器会容易求解一些。与I3积分器相比，SI2积分器能够得到更好的结果，但是它会用到更多的CPU时间，除非放松它的错误设置。所以一般来说SI2积分器的错误设置点最好是I3的10倍。;如果希望让积分器以固定时间步运行，可以适当设置HMAX为较小值而将ERROR值设置得较大。如果不希望积分器跳过一些较短的时间历程（如两个薄板直接的接触），也可以适当这样设置。;将INTERPOLATE属性设置为ON可以阻止一些曲线尖峰以改善仿真结果。典型的，在用GSTIFF积分器中将INTERPOPLATE设置为ON可以减缓仿真速度。在用C++求解器中的HHT积分器时，INTERPOLATE属性的减缓效果会非??小。;在用Newton-raphson迭代方法时，如果设置PATTERN=T，求解器会在每个迭代步计算雅克比矩阵（这会非常费时！！！），但是收敛速度仅是平常的二次方。另一方面，设置PATTERN=F:F，那么仅在收敛失败时重新计算雅克比矩阵。还需要注意的是，FORTRAN和C++求解器在这方面有一个不同的地方。用C++求解器时，设置PATTERN=F表示用积分器来控制计算雅克比矩阵的频率。;要记住，用的次数太多会费时，但太少又会导致收敛失败。当设置PATTERN=F:F，C++求解器会试着去自动寻找适合的计算雅克比矩阵的频率。总的来说，在C++求解器中设置PATTERN=F：F能包含FORTRAN求解器中设置PATTERN＝F的功能。;在模型用Newton－Raphson程序且难收敛时，CORRECTOR属性会比较有用。选择MODIFIED　corrector，积分器可更好的令求解收敛。这种效果在有接触和不连续的模型中更为明显。用这种方法需要很细心，因为它也可能会导致计算结果与真实值不一致。 ;;