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进行谱分解可得
假定矩阵未有重根(重根情况自行讨论),则可得到
那么下式即为该递推矩阵的谱半径:
因此可以看出