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

非线性动力学方程的求解方法 1、概述 在工程实际问题中，我们常常面临这样的选择：我们所遇的问题究竟是静力的还是动力的。静力问题与动力问题，从力学的角度看就是是否考虑与加速度有关的力，而从数学求解方法看则是一个三维边值问题还是一个四维边值-初值问题。在这个问题的选择上没有固定的原则，一般取决于我们研究者、分析者对工程问题的判断。一般认为，实际工程大都是处于动力环境之中，因而属于动力问题。但是，由于时间、经费等方面的原因的限制，我们不可能把所有的问题都按照动力问题的方法来分析。对于许多具体的问题，与速度和加速度有关的力足够小，但是又影响结构分析结果的，将采用静力假定来模拟这些力。 线性的动力有限元控制方程如式(1-1)所示。 (1-1) 式中[M][D][K]分别为结构的质量、阻尼和刚度矩阵，{R}为荷载列矢量，、和分别是加速度、速度和位移列矢量。 式(1-1)的解法大体上可以分为两类：直接积分法和模态叠加法。直接积分法在对控制方程进行数值积分之前不对方程做任何形式的变换，直接用数值积分的方法在时域上一步一步地对方程进行积分。模态叠加法是在求解之前对方程进行某种数学变换，使基底降低，或使矩阵的带宽减小，再进行求解。这两种方法在形式上不同，但是密切相关。 上述每一类求解方法中又有许多具体的解法，每一种解法又有各自的特点。因此我们在选择一种方法求解