非线性动力学方程的求解方法.docVIP

非线性动力学方程的求解方法 1、概述 在工程实际问题中，我们常常面临这样的选择：我们所遇的问题究竟是静力的还是动力的。静力问题与动力问题，从力学的角度看就是是否考虑与加速度有关的力，而从数学求解方法看则是一个三维边值问题还是一个四维边值-初值问题。在这个问题的选择上没有固定的原则，一般取决于我们研究者、分析者对工程问题的判断。一般认为，实际工程大都是处于动力环境之中，因而属于动力问题。但是，由于时间、经费等方面的原因的限制，我们不可能把所有的问题都按照动力问题的方法来分析。对于许多具体的问题，与速度和加速度有关的力足够小，但是又影响结构分析结果的，将采用静力假定来模拟这些力。 线性的动力有限元控制方程如式(1-1)所示。 (1-1) 式中[M][D][K]分别为结构的质量、阻尼和刚度矩阵，{R}为荷载列矢量，、和分别是加速度、速度和位移列矢量。 式(1-1)的解法大体上可以分为两类：直接积分法和模态叠加法。直接积分法在对控制方程进行数值积分之前不对方程做任何形式的变换，直接用数值积分的方法在时域上一步一步地对方程进行积分。模态叠加法是在求解之前对方程进行某种数学变换，使基底降低，或使矩阵的带宽减小，再进行求解。这两种方法在形式上不同，但是密切相关。 上述每一类求解方法中又有许多具体的解法，每一种解法又有各自的特点。因此我们在选择一种方法求解一个问题时，要对该方法的收敛性、稳定性、效率、精度和费用等进行一些分析，讨论它对所求问题的有效性，从而使我们能够针对某一特定的问题，选择合适的方法。 直接积分法基于以下两条： (1)不是在求解时间区间内任意时刻t都满足式(1-1)，而是在相隔△t上的一些离散时刻满足式(1-1)。 (2)对位移、速度和加速度在每一时间区间△t内变化的形式进行假设，事实上若把式(1-1)看成一个常系数微分方程组，便可以用任何一种有限差分格式通过位移来近似表示速度和加速度，因此不同的差分格式就得到不同的方法。 从差分格式上看，分显式与隐式两大类方法，所谓显式差分法就是不必对方程进行求解，而是由前一时刻t的平衡条件直接可求解△t增量后t+△t时刻的各参数解。而隐式差分法则必须对方程进行求解。显式差分法有中心差分法等。隐式差分法有Wilson-法和Newmark法。不同方法解的精度、稳定性、收敛性、效率不同。 2、中心差分法 2.1中心差分法简述 中心差分法的差分格式为： (2-1) (2-2) 认为和在任意时刻t上满足平衡方程(1-1)，即： (2-3) 将式(2-1)和式(2-2)代入到式(2-3)中，可得： (2-4) 2.2中心差分法的特点 (1)的解是基于利用在时刻t的平衡条件，即是在假定t时刻式(2-3)成立的条件下来计算的，该积分过程称为显式积分法，因此中心差分法为显式差分法。但是，用式(2-4)计算时，不仅要用到而且还要用到。所以，计算在时刻t+△t的解，必须存在一个具体的起始过程，即必须利用初始条件。这样，已知和，可由式(2-4)求得，由式(2-1)和式(2-3)可求得，如式(2-5)所示。 (2-5) 我们称式(2-5)为差分格式(2-4)的初始条件。 (2)在求解时不需对刚度矩阵[K]进行三角分解。在应用此法时一般采用集束质量矩阵或称对角质量阵，而阻尼矩阵也通常为对角形式，这样，在运用式(2-4)时，就不需要对等号右端的系数矩阵)进行三角分解，从而可以节省计算时间。 如果不考虑系统的阻尼，则，这时式(2-4)可简化为： (2-6) 其中： (2-7) 从式(2-6)和式(2-7)可以看出，若质量阵[M]为对角矩阵，即在系统中只考虑集束质量的情况，则方程(2-6)实际上是解耦的，即方程组中是相互独立的，因此只需进行矩阵相乘就可得到。 把式(2-6)改写为： 写成分量形式，即： (2-8) 其中，和分别为和的第i个分量，而为质量矩阵的第i个对角线元素。对于，我们可以用静凝聚的方法进行处理。 若对总体刚度和质量矩阵不必进行三角分角，就意味着不必形成总体的[K]和[M]，因此，式(2-7)的计算只需在单元一级进行即可，即： (2-9) 式中[Ke]和[Me]分别为与待求节点自由度有关的单元刚度矩阵与质量矩阵。所以，使用式(2-8)和式(2-9)构成的中心差分具有明显的优点。它不需要计算总体刚度矩阵[K]和质量矩