结构动力响应 返回
Structural Dynamics

结构动力响应模拟器(单自由度系统)

对单自由度 (SDOF) 结构施加正弦、冲击、随机或 El Centro 地震激励,通过 RK4 数值积分实时计算响应,可视化位移、加速度时程与频率响应函数 (FRF)。

参数设置
激励类型
质量 m
kg
刚度 k
阻尼比 ζ
激励频率 f_ex
Hz
幅值 F₀
N
计算结果
ωₙ [rad/s]
固有频率 fₙ [Hz]
最大位移 [mm]
动力放大系数 β
ωd [rad/s]
阻尼固有周期 Td [s]
最大加速度 [m/s²]
阻尼比 ζ
时程响应
频率响应函数 |H(ω)| vs ω/ωₙ
理论与主要公式

单自由度系统运动方程(基底激励时加入 $\ddot{u}_g$):

$$m\ddot{x}+ c\dot{x}+ kx = F(t)$$

其中 $c = 2\zeta\sqrt{km}$。固有角频率 $\omega_n = \sqrt{k/m}$,阻尼固有角频率 $\omega_d = \omega_n\sqrt{1-\zeta^2}$

频率响应函数:$|H(\omega)| = \dfrac{1}{\sqrt{(1-r^2)^2+(2\zeta r)^2}}$,$r=\omega/\omega_n$

共振时($r=1$):$|H| \approx \dfrac{1}{2\zeta}$(动力放大系数最大值)

什么是结构动力响应模拟器(单自由度系统)

🙋
“结构动力响应”听起来好复杂,它到底是什么呀?
🎓
简单来说,就是看一个简单的结构(比如一个重物挂在弹簧上)被“晃”或者“敲”的时候,它会怎么动。在实际工程中,这就像是研究一座桥在车开过时,或者一栋楼在地震时是怎么摇晃的。你可以试着在模拟器里把“激励频率”的滑块拖到5Hz附近,看看当外界的晃动频率接近结构本身的“固有频率”时,会发生什么有趣的现象。
🙋
诶,真的吗?那为什么接近那个“固有频率”就会很危险呢?
🎓
这就涉及到“共振”了。想象一下你推秋千,如果每次都正好在秋千荡到最高点时推一下,它就会越荡越高。结构也一样,当外界的晃动节奏和它自己“想”晃的节奏一致时,能量会不断累积,导致晃动幅度变得非常大。工程现场常见的是,一座桥如果设计时没考虑好,风一吹的频率正好对上桥的固有频率,就可能发生剧烈晃动甚至破坏。你可以在模拟器里把“阻尼比”调小到0.01,再让激励频率等于固有频率,看看位移响应曲线会变得多“夸张”。
🙋
原来如此!那旁边那个“El Centro地震波”选项又是干嘛的?和正弦波晃起来不一样吗?
🎓
非常不一样!正弦波是规律、重复的晃动,而地震波是真实记录下来的、杂乱无章的剧烈震动。比如在汽车碰撞试验中,冲击力也是这种突然、不规则的。选择El Centro波,你就能看到结构在真实、复杂的地震力下的表现了。改变参数后你会看到,即使地震波里没有正好等于固有频率的成分,结构也可能因为其他频率的激励而产生不小的响应,这就是为什么抗震设计不能只看共振。

物理模型与关键公式

这是描述单自由度系统运动的牛顿第二定律核心方程,它平衡了惯性力、阻尼力、弹性恢复力和外部激励力。

$$m\ddot{x}+ c\dot{x}+ kx = F(t)$$

$m$:质量 (kg);$c$:阻尼系数 (N·s/m);$k$:刚度 (N/m);$x, \dot{x}, \ddot{x}$:位移、速度、加速度;$F(t)$:随时间变化的外力。

动力放大系数公式,用于量化系统在简谐激励下,动力响应幅值与静力响应幅值的比值,是评估共振风险的关键。

$$\beta = |H(\omega)| = \frac{1}{\sqrt{(1-r^2)^2 + (2\zeta r)^2}},\quad 其中 \ r = \frac{\omega}{\omega_n}$$

$\beta$:动力放大系数;$r$:频率比;$\zeta$:阻尼比;$\omega_n = \sqrt{k/m}$:固有角频率。当 $r=1$ 且阻尼很小时,$\beta \approx 1/(2\zeta)$,响应会急剧放大。

现实世界中的应用

建筑抗震设计:工程师将复杂的建筑简化为单自由度模型,使用类似本模拟器的工具进行初步分析,估算在地震波(如El Centro波)作用下的最大层间位移和加速度,为结构设计提供依据。

机械振动隔离:对于精密仪器(如光刻机)或重型设备,通过调整系统的质量、刚度和阻尼(即模拟器中的 m, k, ζ),使其固有频率远离工作环境的干扰频率,从而有效隔离振动,保护设备。

车辆悬挂系统调校:汽车的悬挂系统可以看作一个单自由度系统。工程师通过模拟不同刚度(弹簧)和阻尼(减震器)组合下的路面冲击响应,来平衡乘坐舒适性与操控稳定性。

CAE软件验证与教学:在使用Ansys、Abaqus等大型CAE软件进行复杂结构动力分析前,可以用此单自由度模型的计算结果来验证软件设置和线性动力分析模块的正确性,同时也是学习动力学的绝佳可视化工具。

常见误解与注意事项

开始使用这个模拟器时,有几个容易误解的地方。首先,人们常认为“增大质量(m)会使结构更不易摇晃,从而更安全”,但这种想法只对了一半。虽然固有频率 $f_n = \frac{1}{2\pi}\sqrt{k/m}$ 确实会随质量增大而降低,但当地震这类“加速度输入”作用时,情况就不同了。观察运动方程 $m\ddot{x}+ c\dot{x}+ kx = -m\ddot{x}_g$ 的右侧,会发现外力与质量成正比增大。也就是说,过重的质量块反而可能承受更大的惯性力。平衡是关键。

其次,认为阻尼比ζ“越大越好”也是一个误区。虽然振动会更快衰减,但由于阻尼力与速度成正比,如果将阻尼系数 $c$ 设置得过大,虽然能抑制位移响应的峰值,但作用在构件上的阻尼器自身力($c\dot{x}$)可能会变得非常大。例如,将ζ设为0.3时,El Centro波作用下的最大位移会减小,但如果不计算该时刻的阻尼力,支撑部位可能无法承受。设计总是需要权衡的。

最后,不要过度自信地认为“因为是用单自由度计算的,所以实际结构就会完全按此运动”。实际结构具有多个模态(振动形态),这里仅仅考察了第一阶模态。例如,将高耸储罐用单自由度模型简化时,即使将质量(m)设为总重量、用基础弹簧刚度(k)来代表整体刚度,也无法捕捉顶部的液体晃动现象——那完全是另一种物理现象。请记住,这个工具旨在“帮助理解现象本质、体会参数敏感性”,而非直接用于实际设计。