钢筋混凝土悬臂梁含四根与六根纵向钢筋布置方案的三维有限元模态分析对比

《Eng》：Comparative Finite Element Modal Analysis of Reinforced Concrete Cantilever Beams with Four- and Six-Bar Longitudinal Reinforcement Layouts

　　

研究人员采用三维有限元分析（FEA）与一项探索性机器学习（ML）练习，比较了含四根与六根直径12 mm纵向钢筋的钢筋混凝土悬臂梁的弹性模态特性。梁几何尺寸、材料属性、粘结假定及边界条件均保持不变。网格由名义单元尺寸50 mm细化至25 mm时，基频变化0.18%，满足所采用的1%收敛判据。Fine网格FEA基频36.048 Hz与欧拉–伯努利换算截面估计值35.52 Hz相差1.49%。Block Lanczos提取给出前六阶固有频率及相应归一化振型。纵向钢筋面积由452.39 mm2增至678.58 mm2时，计算频率变化为?0.02%至+0.19%，两种布置在假定完好、线弹性、完全粘结条件下振型拓扑相似。以模态序号与钢筋根数为输入、固有频率为输出，对12条相关模态-布置记录拟合四种回归算法。所报R2、MAE与RMSE描述全数据集拟合或插值，未确立对未见布置的泛化能力。因特征向量缩放任意且未执行带定义激励与阻尼的受迫响应分析，故未报任何物理模态位移、应变或应力幅值。

研究背景方面，近数十年钢筋混凝土结构在地震、风、交通、机器振动与冲击等动载下的性能备受关注，固有频率与振型是抗共振设计与健康监测的基础参数。纵向钢筋同时贡献抗弯刚度与质量，对频率产生竞争影响，需在受控数值条件下系统分离。既有文献多关注损伤识别、连接刚度与三维FEA校验，但少有在几何与材料固定时孤立考察纵向钢筋数量对悬臂梁前六阶模态的定量影响，也缺乏对小样本相关仿真数据上多回归算法拟合行为的透明说明。该研究由此在《Eng》发表，旨在以三维FEA为主、探索性ML为辅，回答“钢筋根数变化是否显著改变模态”这一问题。

关键技术方法上，研究人员在ANSYS Mechanical 2023 R1中建立250×500×5000 mm悬臂梁，混凝土用SOLID186二十节点六面体单元、钢筋用LINK180杆单元并以嵌入约束实现平动协调（完全粘结），材料取线弹性（混凝土E_c=35 GPa、ρ_c=2400 kg/m³；钢筋E_s=200 GPa、ρ_s=7850 kg/m³），一端全固支。做100/50/25 mm三级网格收敛，以基频<1%变差为判据选25 mm Fine网格（约320416单元）。用Block Lanczos解无阻尼自由振动特征值问题提取前六阶。解析端以换算截面法（模量比n=E_s/E_c）算欧拉–伯努利悬臂基频校验FEA。ML端由12条记录（两布置×六阶）以模态号与钢筋数（4或6）为输入、频率为输出，在全样本上拟合高斯过程（GP）、随机森林（RF）、梯度提升（GB）与多层感知机（NN），仅报内样本R²/MAE/RMSE，不做交叉验证。

研究结果各节要点如下。3.1计算框架概述：分模型定义、离散化校验、模态提取、数据集生成、探索回归与双布置直比五阶段。3.2几何与布置：四根（A_s=452.39 mm²，配筋率0.362%）与六根（A_s=678.58 mm²，0.543%）12 mm筋、保护层40 mm，截面常数。3.3材料：均质各向同性线弹性，仅用于特征值提取。3.4有限元建模：SOLID186+LINK180嵌入平动兼容，无粘结滑移模型。3.5边界与加载：固支端UX=UY=UZ=0，无外载，仅归一化振型。3.6模态流程：解K?=ω²M?，Block Lanczos取前六阶。3.7解析校验：换算截面得f_1,theory=35.52 Hz，FEA=36.048 Hz，偏差1.49%。3.8生成式AI使用：仅ChatGPT辅助语言与图件润色，不涉及数据与结论。

4.1数据集与特征：12条相关记录，输入仅用先验可得的模态号与钢筋数，排除任何FEA派生量防泄漏。4.2 ML模型：GP(RBF核)、RF(100树)、GB(η=0.1,深度3)、NN(64-32 ReLU,Adam,L2=0.001)。4.3评估：R²/MAE/RMSE全样本计算，明确定义为内样本拟合非泛化。4.4受控FEA直比：四变六根时各阶频率相对变?0.02%~+0.19%，最大绝对差0.30 Hz（六阶），模态内一→六阶跨度四根增150.08 Hz、六根增150.31 Hz。

5.1网格收敛：50→25 mm基频变0.18%<1%，Fine网格采纳；FEA较欧拉解高1.49%源于三维实体含剪变形与转动惯量。5.2模态结果：前六阶四根36.05/69.39/97.39/171.78/178.47/186.13 Hz，六根略高；振型序列为侧向弯→竖向弯→二阶侧弯→弯扭耦合→二阶竖弯→三阶侧弯，归一化仅看相对形貌。5.3双布置对比：频率差<0.2%、振型拓扑同序，因钢面积占截面0.362%~0.543%，刚度与质量增量近抵消，√（k/m）基本不变。5.4探索回归：GP/GB/NN全样本R²=1.0000，RF=0.5202；高R²属12点小样本插值而非泛化。5.5两离散布置描述性比较：增50%钢面积不改振型序、不提供优化或安全排序。5.6工程释义：仅改筋数非调频有效手段，调频需动截面/材料/边界；本结果限模态对比不作配筋推荐。

讨论部分强调：小样本相关模态记录不可支撑未见布置验证；特征重要性不解读为物理因果；完全粘结线弹性范围外（裂后、滑移、非线性）结论不适用；物理应力幅值须靠受迫响应分析。结论段译文如下：本研究比较了含四根与六根12 mm纵向筋钢筋混凝土悬臂梁的弹性模态特性。名义单元50→25 mm基频变0.18%，满足1%网格收敛；Fine网格基频36.048 Hz与换算截面欧拉–伯努利估计35.52 Hz差1.49%；钢筋面积452.39→678.58 mm²使前六阶频率变?0.02%~+0.19%，最大差0.30 Hz；完好线弹性完全粘结下两模型归一化振型拓扑相似。上述仅限两离散布置，不立最优配筋。探索ML仅以模态号与筋数为入、频率为出，R²/MAE/RMSE基于12条全样本描述拟合或插值而非泛化；未评物理位移/应变/应力幅值，故不就应力放大、五阶临界性、疲劳、利用率、安全系数或配筋效率下结论。