符号	名称	单位
De	外径	mm
Di	内径	mm
h0	锥高	mm
t	厚度	mm

DIN 2093 显式动力学分析：冲击与瞬态强非线性问题

1. 显式动力学的基本思想

显式动力学采用中心差分法直接对时间积分，无需组装和求逆切线刚度矩阵，也无迭代收敛问题。因此，它特别适合碟形弹簧在瞬态冲击、高速碰撞、成形过程、以及高度不连续的非线性问题（如复雑接触、材料断裂）。

2. 运动方程的显式积分

在时间步 $$n$$ ，离散的动力学方程为：

\mathbf{M} \mathbf{a}_n = \mathbf{F}_{ext,n} - \mathbf{F}_{int,n}

其中 $\mathbf{M}$ 为对角质量阵， $\mathbf{a}_n$ 为节点加速度向量。由此可直接解耦求解加速度。速度和位移通过中心差分格式递推：

\mathbf{v}_{n+1/2} = \mathbf{v}_{n-1/2} + \mathbf{a}_n \Delta t

\boxed{\mathbf{u}_{n+1} = \mathbf{u}_n + \mathbf{v}_{n+1/2} \Delta t}

整个过程无需形成全局刚度阵，亦无迭代，每一步计算成本极低。

3. 临界时间步（Courant‑Friedrichs‑Lewy 条件）

显式方法是条件稳定的，时间步长 $\Delta t$ 必须小于临界值 $\Delta t_{crit}$ 以确保数值稳定。物理上，这意味着信息在一个时间步内不能穿越单个单元。该条件表达为：

\boxed{\Delta t_{crit} \approx \frac{L_{min}}{c_d}}

其中： - $L_{min}$ — 模型中最小单元的特征长度（mm）。对六面体单元，通常取最小边长；对四面体，取最小高。 - $$c_d$$ — 材料的一维纵波波速（mm/s），按下式计算：

c_d = \sqrt{\frac{E \cdot (1-\nu)}{(1+\nu)(1-2\nu)\rho}}

$$E$$ — 弹性模量（MPa）
$\nu$ — 泊松比
$\rho$ — 材料密度（t/mm³ 或 kg/mm³），注意单位一致性。

对于钢制碟簧： $E \approx 206\,000\ \text{MPa}$ ， $\nu=0.3$ ， $\rho \approx 7.85\times10^{-9}\ \text{t/mm}^3$ ，计算得 $c_d \approx 5.2\times10^6\ \text{mm/s}$ 。

若碟簧有限元网格中最小单元尺寸 $L_{min} = 0.2\ \text{mm}$ ，则：

\Delta t_{crit} \approx \frac{0.2}{5.2\times10^6} \approx 3.85\times10^{-8}\ \text{s}

这正是显式动力学中时间步长通常落在 $10^{-8}\sim10^{-7}$ 秒量级的原因。单元越小，材料声速越高，临界步长越短，计算量越大。

4. 碟形弹簧冲击分析的典型设置

项目	建议设置
分析类型	Explicit Dynamics (显式)
时间步控制	自动（基于最小单元），可引入少量质量缩放以增大步长
单元类型	一阶缩减积分六面体（耐大变形，效率高）
网格密度	厚度方向 ≥ 4 层，内外缘加密，防止最小单元过小
材料模型	弹塑性，可含应变率效应（如 Cowper‑Symonds）
接触	自动面-面接触，摩擦系数 $\mu_f$ （0.03~0.10）
加载	定义冲击初速度 $$v_0$$ （或落锤质量），施加在碟簧支承结构上
输出	总反力‑时间，应力‑时间，碟簧压缩量‑时间，动能‑内能转换

5. 质量缩放与计算效率

为提高计算速度，可在保证精度的前提下通过质量缩放人为增大时间步长。基本关系为：由于 $\Delta t_{crit} \propto L / \sqrt{E/\rho}$ ，若将最小单元密度增加 $\alpha$ 倍，则临界步长增大 $\sqrt{\alpha}$ 倍。质量增加应控制在总质量的 < 2%，以免显著改变动态响应。

6. 显式方法与隐式方法的对比

项目	显式动力学	隐式静力/动力学 (Newton‑Raphson)
适用问题	短时瞬态、冲击、强非线性、波动传播	静态、准静态、低频振动、中等非线性
时间步长	极小（ $10^{-8}\sim10^{-7}$ s），单元尺寸决定	较大，载荷步可任意（迭代收敛决定）
每步计算成本	低（无矩阵求逆）	高（组装并分解切线刚度阵）
总时间	适合毫秒级事件；秒级过程极耗时	适合秒级及以上静力或缓慢事件
稳定性	条件稳定	无条件稳定（标准方法）
碟簧典型应用	落锤冲击、爆炸、碰撞缓冲	准静态压缩、预紧力加载、模态分析

7. 计算示例（概念）

碟簧组承受质量块 $m = 50\ \text{kg}$ ，初速 $v_0 = 5\ \text{m/s}$ 的冲击。
最小单元尺寸 $L_{min} = 0.15\ \text{mm}$ ，钢 $c_d \approx 5.2\times10^6\ \text{mm/s}$ → $\Delta t_{crit} \approx 2.88\times10^{-8}\ \text{s}$ 。
总冲击持续时间约 2 ms，需要的增量步数 ≈ $2\times10^{-3} / (2.88\times10^{-8}) \approx 69\,400$ 步。此规模可在几小时内完成。输出最大压缩量、力峰值及能量吸收比，与能量法对比验证。

8. 注意事项

网格不宜过细：极小单元会严重减小时间步，应合理控制最小单元尺寸，平衡精度与效率。
避免沙漏模式：使用缩减积分单元时，需检查沙漏能（< 内能的 5%），必要时增加沙漏控制。
阻尼与材料率效应：冲击问题建议引入材料应变率敏感性，否则可能低估峰值应力。
后处理重点：关注力‑位移滞回环、最大压缩量、应力波传播及可能的局部塑性变形。

总结：显式动力学分析以中心差分法和临界时间步控制为核心，能够精确模拟碟形弹簧在毫秒级冲击下的瞬态响应。正确设定单元尺寸、材料参数和接触，可得到可靠的力、应力与能量数据，为缓冲与防护设计提供依据。

显式动力学

公式表达式

参数列表

需要计算该公式？

详细计算指南