结构力学数值方法：积分法：结构动力学积分法教程.pdf

结构力学数值方法：积分法：结构动力学积分法教程

1绪论

1.1结构动力学的基本概念

结构动力学是力学的一个分支，主要研究结构在动态载荷作用下的响应。

动态载荷可以是周期性的（如风、海浪、机器振动），也可以是非周期性的（如

地震、爆炸）。结构动力学分析通常涉及解决二阶微分方程，这些方程描述了结

构的质量、刚度和阻尼如何影响其动态行为。

1.1.1质量矩阵

质量矩阵反映了结构各部分的质量分布，通常是一个对角矩阵，其中对角

线元素表示各节点的质量。

1.1.2刚度矩阵

刚度矩阵描述了结构在变形时的恢复力，是结构分析中的关键矩阵。在有

限元分析中，刚度矩阵由各个单元的刚度矩阵组合而成。

1.1.3阻尼矩阵

阻尼矩阵表示结构的能量耗散特性，可以是粘性阻尼、库仑阻尼或结构阻

尼。阻尼矩阵的计算较为复杂，有时采用简化模型。

1.2数值积分法在结构动力学中的应用

数值积分法是解决结构动力学问题的重要工具，尤其是对于复杂结构和非

线性问题。常见的数值积分法包括显式积分法和隐式积分法。

1.2.1显式积分法

显式积分法如欧拉法、中央差分法和Newmark-β法中的β=0.5,γ=0.5（也

称为Wilson-θ法，θ=1.0）等，计算速度快，但稳定性条件严格，时间步长受

限。

1.2.2隐式积分法

隐式积分法如Newmark-β法（β0.5,γ0.5）和HHT-α法等，虽然计算

量大，但具有更好的稳定性，可以使用较大的时间步长。

1

1.3结构动力学积分法的发展历程

结构动力学积分法的发展经历了从简单的欧拉法到更复杂的Newmark-β法、

HHT-α法和广义α法等。这些方法的改进主要集中在提高计算的精度、稳定性

和效率。

1.3.1欧拉法

欧拉法是最简单的数值积分法，分为前向欧拉法和后向欧拉法。前向欧拉

法不稳定，而后向欧拉法虽然稳定，但精度较低。

1.3.2Newmark-β法

Newmark-β法由Newmark在1959年提出，通过调整β和γ的值，可以在

精度和稳定性之间找到平衡。当β=0.25,γ=0.5时，Newmark-β法成为线性加

速度法，具有无条件稳定性。

1.3.3HHT-α法

HHT-α法由Hilber、Hughes和Taylor在1977年提出，通过引入一个额外的

参数α，可以进一步提高方法的稳定性，同时保持较高的精度。

1.3.4广义α法

广义α法由Chung和Hulbert在1993年提出，是一种高阶的隐式积分法，

通过调整α的值，可以在精度、稳定性和效率之间找到最佳平衡。

1.3.5代码示例：Newmark-β法

importnumpyasnp

defnewmark_beta(M,C,K,F,u0,v0,dt,beta=0.25,gamma=0.5,t_end=10):

Newmark-β法求解结构动力学问题

:paramM:质量矩阵

:paramC:阻尼矩阵

:paramK:刚度矩阵

:paramF:力向量

:paramu0:初始位移

:paramv0:初始速度

:paramdt:时间步长

:parambeta:Newmark-β法的β参数

:paramgamma:Newmark-β法的γ参数

:paramt_end:模拟结束时间

2

:return:位移、速度和加速度的时间历程

n=len(M)

u=np.zeros((n,int(t_end/dt)+1))

v=np.zeros((n,int(t_end/dt)+1))

a=np.zeros((n,int(t_end/dt)+1))

u[:,0]=u0

v[:,0]=v0

#计算初始加速度

a[:,0]=np.linalg.solve(M,F[:,0]-C@v0-K@u0)

foriinrange(1,int(t_end/dt)+1):

#计算预测位移

u_pred=u[:,i-1]+dt*v[:,i-1]+dt