使用集中参数模型进行铁路耐撞性分析的降阶建模耐撞性评估

《Modelling》：Crashworthiness Assessment Using Lumped Parameter Models for Reduced-Order Modelling in Railway Crashworthiness Analysis

　　

铁路客车的设计必须满足严格的认证要求，尤其是在欧洲标准EN 15227下的耐撞性分析中。使用全尺寸有限元模型（FEM）进行此分析极其耗费时间、计算能力和工程资源，即使采用大型服务器集群也是如此。为提高效率，简化结构中受外部载荷影响较小的区域是有益的。在此方法中，研究人员将较不关键的部分替换为柔性一维单元，减少自由度数量，同时保留车辆的主要动态行为。通过聚焦于特定的中跨段，模型变得更稳健且易于管理。引入校准单元以准确再现被移除结构部件的质量和刚度。该方法还集成了质量和刚度单元，以在更宽的频率范围内捕捉结构响应。随后应用一种迭代无梯度校准程序调整等效刚度和质量分布，使简化模型再现全尺寸参考模型的响应。结果表明，该策略有效，仿真时间减少了77.6%，同时保持了可靠的精度。然而，该过程仍为劳动密集型，且在大变形条件下其性能可能下降。

论文解读文章

**研究背景与问题**
轨道车辆耐撞性（crashworthiness）是保障乘客安全的关键被动安全措施。随着列车运行速度提高，碰撞后果愈发严重，因此必须满足严格认证标准（如EN 15227:2020）。全尺寸有限元（FEM）模型虽能精确预测碰撞行为，但计算耗时巨大，即使使用高性能计算集群也效率低下。现有简化方法（如Huňady等采用刚度衬套、Sch?ffer等基于刚度和失效的替代模型）虽能降低计算成本，但缺乏对轴向和弯曲刚度的同时保留，且未完全适配标准碰撞场景。Xie等针对高速列车的简化模型存在对齐和一致性不足的问题。因此，亟需一种系统化降阶建模（reduced-order modelling）方法，在保持关键动态响应的前提下大幅减少仿真时间，并嵌入标准化评估框架。

**研究内容与结论**
研究人员提出并验证了一种专用简化方法论：将客车的非关键中间段替换为等效集中参数单元（lumped-parameter elements），同时保留被移除部分的质量分布、轴向刚度和弯曲刚度。通过迭代无梯度校准程序调整等效刚度与质量，使简化模型再现全尺寸参考模型的响应。验证结果表明：简化模型在18 km/h刚性墙碰撞仿真中，仿真时间降低77.6%，内能、动能、反力及垂直位移的平均误差均低于4%，等效塑性应变分布与全尺寸模型高度一致。将该方法应用于EN 15227场景1的完整列车碰撞模型时，尽管在0.2秒后出现非物理行为（如缓冲器错位、底架脱开），但0.2秒内的结果仍可用于评估耐撞性合规性。合规性评估揭示该客车存在刚度分布不均、负载路径不足等问题，导致生存空间缩小和减速度超标，需引入能量吸收装置（如碰撞盒）。该研究发表在《Modelling》。

**主要关键技术方法**
1. **区域划分策略**：基于全尺寸模型初步仿真，识别高应力、塑性应变和显著变形的区域（保留为详细FE子结构），而弹性主导、仅贡献整体刚度和惯性的区域（如中间段）被替换为等效单元。
2. **等效刚度识别**：通过静态压缩试验（2000 kN加载）获取客车整体轴向和弯曲刚度，利用欧拉-伯努利梁（Euler-Bernoulli beam）等效模型推导轴向刚度（k₁₁ = 958.92 N/mm）和弯曲刚度（k₃₃ = 1,702,020,535 N·mm/rad）。
3. **迭代无梯度校准**：逐步将中间段替换为弹簧单元组，调整各弹簧的刚度和分布质量（按地板、车顶、侧墙几何分配），使简化模型在能量、反力、位移上逼近全尺寸模型。
4. **材料与连接模型**：采用Johnson-Cook塑性模型（不含应变率和温度效应），焊点连接基于最弱材料假设，使用Material Law 59和耦合法向-剪切位移失效准则。

**研究结果**
**4.1 简化模型开发**
通过比较全尺寸模型在基准碰撞工况下的应力、塑性应变和变形分布，将除前端和后端外的中间非关键区域（图10）替换为等效弹簧单元。弹簧的轴向刚度和弯曲刚度来自全局刚度评估，质量按原分布分配到三组弹簧（地板、屋顶、侧壁）。经迭代调整后，简化模型实现77.6%仿真时间缩减。

**4.2 全局刚度评估**
静力分析中施加2000 kN压缩载荷，结合弯曲偏心距计算弯矩，利用2节点Euler-Bernoulli梁有限元平衡方程识别出等效轴向刚度k₁₁和弯曲刚度k₃₃等参数，为弹簧赋值提供依据。

**4.3 简化模型验证**
对比全尺寸模型与简化模型在18 km/h刚性墙碰撞下的响应：内能平均误差1.13%（最大2.25%），动能平均误差2.43%（最大6.75%），反力平均误差2.67%；垂直位移节点1、2误差0%，节点3、4平均误差分别为2.27%和3.26%。等效塑性应变和von Mises应力分布一致，表明简化模型保留了主要碰撞行为。

**5.1 EN 15227合规性评估**
**5.1.1 爬车防止（Overriding Prevention）**
至少两个不同转向架的轮对保持轨道接触，垂直位移未超过轮缘高度的75%（21 mm），满足标准要求。

**5.1.2 生存空间保持（Survival Space Preservation）**
乘客生存空间在前后端5 m段内纵向缩短超过50 mm，等效塑性应变超过10%阈值，不满足要求；车辆前端变形超标，后端合格；临时占用区域变形极小，满足要求。

**5.1.3 减速度限制（Deceleration Limitation）**
30 ms移动平均减速度超过允许值（10 g）两倍以上，120 ms移动平均减速度约为允许值（5 g）三倍，未达标。

**5.1.4 总结**
刚度分布不均衡，缺少能量吸收特征（如受控折叠区），负载路径缺乏冗余，底架连接处应力集中导致过早屈曲，乘客舱区域未参与能量耗散。

**讨论与结论**
讨论指出当前简化方法虽实现77.6%仿真时间缩减且平均误差低于4%，但在大变形场景下表现不稳定，且开发过程仍属劳动密集型（弹簧和连接器手工定义）。建议未来引入应变率效应、真实连接材料参数，并自动化弹簧建模流程。最终结论强调该方法可接入标准化评估流程，但现有客车设计需优化耐撞性，通过集成能量吸收装置和改进负载路径来提升生存空间与减速度合规性。