【JY】浅析各动力求解算法及其算法数值阻尼(人工阻尼)

        hello~

即刻关注

“通俗讲动力数值算法”

    在谈数值阻尼之前，咱们先聊聊直接积分法的稳定性问题。

    直接积分法是不对运动方程进行任何变换，直接运动方程进行积分求解。本质上讲，直接积分法是基于两个基本概念：

    1.是将在求解域0≤t≤T内任何时刻t都应满足运动方程的要求，代之以仅在一定条件下近似地满足运动方程，例如仅在相隔△t的离散时间点0，△t，2△t，……上满足运动方程。

    2.是在一定数目的△t区域内，假设位移D、速度V和加速度A的近似函数形式。

即在假定已知时刻t=0的位移D、速度V，并将时间求解域 0 ~ T等分成n个时间段0，△t，2△t，…，n△t（△t=T/n），且已求得时刻0，△t，2△t，…和t的解，现在要求t+△t时刻的解。建立了求解t+△t时刻解的算法后，即可类似地求得所有离散时刻的解，因此该方法也称逐步积分法。

 

    详情可以看下历史文章【JY】结构动力学初步-单质点结构的瞬态动力学分析。

 

    逐步积分法是一种近似的求解方法，每一步积分计算都会带来误差。解的误差的来源主要有两大类：

    第一类是在计算速度和加速度的近似表达式中略去了高阶小量，称为截断误差，这类误差是由方法本身决定的，一般可以作出估计，它随步长△t增大而增大。

    第二类是由于计算机的字长总是有限位，超过其位数的数必须四舍五入，这类误差称为舍入误差。尽管在每一步中舍入误差并不大，但这类误差在求解过程中可能会被不断放大，以致使计算结果完全失真，这就是直接积分方法的稳定性问题。

    如果无论△t取多大，给定任意初始条件，积分结果都不会无界增大，则称此方法是无条件稳定的。

    反之，如果△t必须小于某个临界值时，积分结果才不会无界增大；则称此方法是条件稳定的。

 

    在建立直接积分法格式时，我们假定已知时刻0，△t，2△t，…，t-△t和t的解，要求时刻t+△t的解，因此直接积分法利用时刻0，△t，2△t，…，t-△t和t的解递推求解t+△t时刻的解，即直接积分法的格式可以写成递推的形式

    其中 Xt+Δt 和 Xt 为含有位移、速度等结果的向量，矩阵D称为递推矩阵或放大矩阵。

    上式右端的第二项与载荷有关，在进行稳定性分析时可令L=0，此时给定初始条件X0后，时刻 t+△t = (n+1)△t 的解 Xt+Δt 可以写为

    对矩阵D进行谱分解可得

    假定矩阵未有重根(重根情况自行讨论)，则可得到

    那么下式即为该递推矩阵的谱半径：

    因此可以看出