4. 材料模型与阻尼定义
瞬态结构动力学与静力学分析在材料输入上最大的差别,不是弹性模量或泊松比,而是密度和阻尼。静力学里密度几乎不影响结果,但在瞬态分析中,密度直接决定质量矩阵 M,而质量矩阵与刚度矩阵 K 一起决定了结构的固有频率与惯性响应;阻尼则决定振动能量如何随时间衰减。如果这两项设置不当,即使网格、边界条件和载荷都正确,位移峰值、衰减速度和相位也可能与真实情况相差很大。本节围绕线弹性材料参数、密度与弹性模量的输入,以及 Rayleigh 阻尼、恒定阻尼比和材料阻尼三类阻尼的适用场景与设置方法展开。
一、线弹性材料参数:瞬态分析真正需要哪些
Ansys Mechanical 的 Engineering Data 中,线弹性材料至少需要输入三项:密度 Density、杨氏模量 Young's Modulus(或各向同性弹性模量)、泊松比 Poisson's Ratio。对于各向同性线弹性材料,这三项就足以构造质量矩阵和刚度矩阵。若材料是正交各向异性(如复合材料层合板),则需要输入 E1、E2、E3、ν12、ν13、ν23、G12、G13、G23 共九个独立常数,此时密度仍然只有一个标量。
密度输入的单位必须与项目单位制一致。例如采用 mm-t-s 单位制时,钢的密度应输入 7.85e-9 t/mm³,而不是 7850 kg/m³;若误填 7850,质量矩阵会被放大近 10¹² 倍,求解器要么报错,要么给出完全失真的低频响应。杨氏模量同理,mm-t-s 制下钢为 2.1e5 MPa(即 2.1e5 N/mm²)。这是初学者最常见的错误来源之一,建议在 Engineering Data 中先确认单位制,再逐项核对数量级。
此外,瞬态分析中若考虑重力、离心力等惯性载荷,密度还直接参与体力计算;若考虑大变形或接触,材料模型可能需要从线弹性扩展为双线性等向强化等塑性模型,但本节只讨论线弹性范围。
二、三类阻尼的物理含义与适用场景
Ansys 瞬态动力学中可用的阻尼主要包括 Rayleigh 阻尼、恒定阻尼比(Constant Damping Ratio)和材料阻尼(Material Damping)。它们并非互斥,可以叠加使用,但叠加后总阻尼会变大,需要谨慎。
Rayleigh 阻尼是最常用的形式,其阻尼矩阵写作 C = αM + βK,其中 α 为质量比例系数,β 为刚度比例系数。α 主要抑制低频段响应,β 主要抑制高频段响应。若已知结构在两个频率 ω₁、ω₂ 处的目标阻尼比 ζ₁、ζ₂,可由下式反算:
α = 2ω₁ω₂(ζ₁ω₂ − ζ₂ω₁) / (ω₂² − ω₁²)
β = 2(ζ₂ω₂ − ζ₁ω₁) / (ω₂² − ω₁²)
工程上常取两个覆盖主要响应频段的模态频率,并令 ζ₁ = ζ₂ = ζ,此时 α = 2ζω₁ω₂/(ω₁+ω₂),β = 2ζ/(ω₁+ω₂)。Rayleigh 阻尼适合结构阻尼随频率变化不剧烈、且主要关心某一频段的场景,例如一般金属机架、支架的瞬态冲击响应。
恒定阻尼比是直接给定一个与频率无关的 ζ 值,Ansys 在模态叠加瞬态分析中按各阶模态分别施加该阻尼比。它适合已知结构整体阻尼水平、且各阶模态阻尼差异不大的情况,例如焊接钢结构常取 ζ = 0.02~0.05,混凝土结构常取 0.05~0.07。恒定阻尼比在完全法瞬态分析中不能直接使用,需转换为 Rayleigh 阻尼或材料阻尼。
材料阻尼(Material Damping)通过定义材料的损耗因子或阻尼系数,把阻尼直接绑定到材料本构上,适合多材料装配体中不同材料阻尼差异明显的场景,例如橡胶隔振器与钢支架组合。材料阻尼通常以损耗因子 η 或复模量形式给出,在瞬态分析中会转换为等效的刚度比例阻尼。
三、在 Mechanical 中的设置路径
Rayleigh 阻尼在 Analysis Settings 的 Damping Controls 中设置,选择 Rayleigh Damping,输入 Alpha(质量比例)和 Beta(刚度比例)。恒定阻尼比在同一位置选择 Constant Damping Ratio,输入 ζ。材料阻尼则在 Engineering Data 的材料属性中添加 Damping 相关项,或在 Mechanical 中通过 Command 对象以 APDL 命令 MP,DAMP 定义。
需要特别注意的是,Rayleigh 阻尼的 α、β 必须与所选频率单位一致。若模态频率以 Hz 给出,而公式中使用的是 rad/s,需先乘以 2π。例如 ω₁ = 2π×50 rad/s、ω₂ = 2π×200 rad/s、ζ = 0.03,则 α ≈ 2×0.03×2π×50×2π×200 / (2π×50+2π×200) ≈ 15.1 rad/s,β ≈ 2×0.03 / (2π×250) ≈ 3.82e-5 s。这些数值直接填入 Damping Controls 即可。
四、设置后的验证与常见问题
阻尼设置完成后,不能只看求解是否收敛,而应通过结果曲线判断阻尼是否合理。在瞬态分析结果中提取某个关注节点的位移—时间曲线,观察自由振动段的衰减包络。若包络几乎不衰减,说明阻尼过小或未生效;若第一个峰值就被压得很低,说明阻尼过大。可用对数衰减率 δ = ln(A₁/A₂) 反算等效阻尼比 ζ = δ / √(4π² + δ²),与设定值对比,判断设置是否与预期一致。
常见问题包括:Rayleigh 阻尼的 α、β 与频率单位不匹配,导致阻尼量级错误;恒定阻尼比在完全法瞬态分析中被忽略;材料阻尼未在 Engineering Data 中正确关联到对应材料;以及多材料模型中只给一种材料定义了阻尼,其余材料默认无阻尼。排查时建议先检查 Analysis Settings 中 Damping Controls 的实际生效项,再核对材料级别的阻尼定义,最后用位移衰减曲线做定量验证。
总结来说,瞬态结构动力学的材料与阻尼设置可以归纳为一条主线:密度和弹性模量决定“结构本身怎么振”,阻尼决定“振动怎么衰减”。Rayleigh 阻尼适合频段明确的单材料结构,恒定阻尼比适合模态叠加法下的整体阻尼控制,材料阻尼适合多材料装配体。三者可以组合,但必须保证单位一致、频率一致、量级合理,并通过位移衰减曲线做最终验证。