基于逐点时间谐波仿真的薄壁管道面向模态有限元模型修正方法

《Results in Engineering》：Mode-Oriented Finite-Element Model Updating for Thin-Walled Ducts via Point-Wise Time-Harmonic Simulations

　　

未明确材料属性和边界条件是导致薄壁管道振动力学有限元预测不准确的主要原因，而含非可忽略阻尼时的模型修正尤为困难。本文提出一种面向模态(model-oriented)的模型修正流程，仅需在参考测量估计的目标模态频率处进行时间谐波(time-harmonic)仿真评估，从而可直接考虑线性阻尼。通过频率响应函数(Frequency Response Function, FRF)向量的空间相干性最大值——即频域保证准则(Frequency-Domain Assurance Criterion, FDAC)——将仿真频点与实测频谱 bin 配对以构建目标函数。研究人员考察了三种标量化目标函数变体，可选择性地考虑配对后的频率偏移、空间相干性及振动幅值信息。参数辨识采用全局优化结合基于随机森林(Random Forest, RF)回归器的代理模型完成。研究人员在三个含阻尼、复杂度递增的薄壁算例中对方法适用性进行了验证，分别在多个训练集规模下评估：S1在简化组合结构中恢复已知法兰弹性模量与频率常数阻尼；S2为去特征化法兰及安装部位识别5个等效参数以模拟更详细的参考模型；S3以激光多普勒测振仪(Laser Doppler Vibrometry, LDV)锤击激励实测数据为参考，修正真实通风管道模型的14个参数。验证通过对频带内频率平均空间相干性和相对均方误差进行比较来完成。结果表明，所有场景中随训练集增大代理模型预测能力稳步提升，S1成功恢复目标参数，S2较基准模型显著改善，S3取得中等程度改善。最稳健的性能由仅联合频率偏移与FDAC空间相干性惩罚项的目标函数取得。

论文解读：基于逐点时间谐波仿真的薄壁管道面向模态有限元模型修正方法

研究背景与意义

有限元(Finite Element, FE)方法广泛用于薄壁复合结构的振声场仿真，但理想化误差（几何简化、边界/界面条件假设偏差、材料属性不确定）常导致预测失准，尤其当结构中阻尼不可忽略时，传统基于无阻尼特征值问题的模态(model-based)模型修正难以同时处理阻尼参数，且需依赖实验模态分析(Experimental Modal Analysis, EMA)提取真实模态振型，易受空间不完整性和噪声影响。现有频响函数(Frequency Response Function, FRF)法虽可纳入阻尼，但频点选择敏感，模态配对困难。为此，Heidegger P, Czwielong F, Wurzinger A等人于《Results in Engineering》发表本研究，提出一种面向模态(mode-oriented)的FE模型修正工作流：利用仅在参考测量预估目标模态频率处求解的逐点时间谐波(time-harmonic)仿真获取FRF向量，借助频域保证准则(Frequency-Domain Assurance Criterion, FDAC)自动配对仿真与实测频点，构建类模态修正目标函数，并结合随机森林(Random Forest, RF)代理模型辅助差分进化(Differential Evolution, DE)全局寻优，可在不直接提取实测模态的前提下同步修正含非比例阻尼的等效材料参数与边界条件。

主要关键技术方法

研究人员采用三类测试场景验证工作流：(S1/S2)以详细FE仿真结果为参考，(S3)以离心风机测试台通风管道的激光多普勒测振仪(Laser Doppler Vibrometry, LDV)单点锤击激励实测FRF为参考。时间谐波FE方程求解含Rayleigh阻尼或可单元级定义的非比例阻尼；参考FRF峰值频率由单通道SL₂(平方L₂范数)曲线经find_peaks检测获得；仿真与实测FRF向量按最大FDAC值配对；目标函数J⁽¹⁾仅含归一化频率偏移平方项，J⁽²⁾加入FDAC残差(1?FDAC)，J⁽³⁾再加入对数SL₂幅值差项。参数空间用随机置乱Sobol序列采样，openCFS求解FE，pyCFS做后处理插值(k近邻Shepard法)，scikit-learn训练RF回归器(n_estimators=150, min_samples_leaf=2, 按目标值加权样本权重, OOB袋外评估)，SciPy差分进化算法(rand/1/bin策略, maxiter=5000)在RF代理曲面上寻优，最终以频带内平均FDAC系数κ^(FDAC)和相对均方误差κ^(MSE)评价全频段更新质量。

研究结果

8.1 S1) 工作流收敛性验证

研究人员以已知参数的简化组合薄壁结构为参考，仅更新法兰等效弹性模量E和频率常数阻尼比ξ。结果显示RF代理模型OOB R²随训练样本Ns增大单调升高，Ns≥2048时三目标函数均>0.99；更新参数随Ns增大趋近真值，E在Ns=32时均值误差<20%，Ns=4096时J⁽²⁾最佳相对误差0.14%；ξ收敛较慢(RF特征重要性仅0.01~0.036)，Ns≥512时落于10%误差带内，归因于其对目标函数灵敏度较低，SHAP分析佐证此结论。表明工作流在可控场景下可正确收敛并识别含阻尼参数。

8.2 S2) 等效参数识别（去特征化几何）

参考为含详细螺栓法兰与安装块的复杂FE模型，待修正去特征化模型更新5参数：法兰等效E与ξ，安装部位三方向等效法向刚度K_x、K_y、K_z。κ^(FDAC)与κ^(MSE)显示更新模型全面优于零位移BC基准(基线κ^(FDAC)≈0.92)，J⁽²⁾表现最稳(Ns≥32时κ^(FDAC)>0.96，κ^(MSE)≈0.1)。平行坐标图揭示解不唯一(等效参数存在补偿效应，K_x较可辨识)，符合去特征化几何无法被有限等效参数唯一表征的预期；RF特征重要性与SHAP值表明不同目标函数对各参数敏感度权重不同。

8.3 S3) 实际应用——真实通风管道

FE模型含14更新参数(7个部件等效弹性模量E^{(HW/DW/R/MH/FH/FD/C)}，法兰与帆布双频点Rayleigh阻尼比ξ₁^(F/C)、ξ₂^(F/C)，壳体安装法向刚度K_x/y/z)，参考为前表面1453点单次锤击LDV测得的H₁频响函数估计(频带150~800 Hz，因相干性下降截断)。Investigation 1(固定5目标峰，Ns=32~4096)：RF OOB R²最高约0.88(J⁽²⁾)，J⁽³⁾最低(0.66)；κ指标较基线中等改善，J⁽²⁾最稳，但未随Ns持续显著收敛至唯一解。平行坐标证实参数非唯一性(多参数互补偿)，SHAP示E^(HW)、E^(DW)影响最大。Investigation 2(固定Ns=256，变目标峰数2/5/11)：意外发现2目标峰总体评分最优，增多目标峰因参考振型间较高离对角FDAC值(达0.475)致配对歧义增大反使性能下降；配对歧义映射图显示J⁽²⁾在低目标值区配对一致性最好。说明实际复杂结构更新受限于模型简化程度、测量不完整性与目标峰选择，适度少而分离的目标峰可减配对噪声。

讨论与结论翻译

本工作提出了一种用于含阻尼薄壁复合结构有限元模型修正的工作流并详述其开源实现。该工作流由逐点时间谐波仿真构造面向模态的目标函数，并利用随机森林代理模型辅助差分进化全局优化。在三个复杂度递增的薄壁含阻尼算例中验证了参数辨识：S1中成功收敛至已知法兰弹性模量与频率常数阻尼比；S2中更新的等效参数与边界条件使模型大幅优于基准；S3中14参数更新较基线获中等改善，但因真实结构复杂性不能完全消解差异，解具非唯一性。三种目标函数中，仅联合归一化频率偏移与FRF向量空间相干性(FDAC)惩罚项的J⁽²⁾表现最稳健，尤在少量训练样本时；加入振动幅值项(J⁽³⁾)未进一步改善反致代理训练效率降低与配对歧义敏感；仅频率偏移项(J⁽¹⁾)在S2中表现较差。未来可将本工作流与其他修正方法对比，或引入自适应采样策略以细化局部搜索。