基于LS-DYNA的钢筋混凝土梁近场爆炸仿真
关键词:LS-DYNA;钢筋混凝土梁;近场爆炸;ALE流固耦合;破坏模式
一、文章简要介绍
梁在爆炸下是弯坏还是剪坏,工程上一直缺判据,爆点偏离跨中就更难定性。这篇2026年发表在Buildings的研究,做了三组1/2缩尺RC梁的近场空气爆炸试验,又用LS-DYNA的ALE流固耦合做数值模拟,提出用剪切变形比δ划分破坏模式(DOI:10.3390/buildings16112054)。
图1 RC梁近场空气爆炸试验流程图
试验全流程:从浇筑、材性测试到起爆与测量。
图2 试验现场布置图
梁两端各150毫米夹紧,TNT按预设高度悬挂。
二、仿真步骤
梁为1/2缩尺,截面125×250毫米、净跨2700毫米,C40混凝土实测抗压40.5兆帕,上部2C10、下部4C10,箍筋B6@100。三个工况:2千克装药跨中起爆,4千克装药分别偏心0.3米、0.6米,爆高都是0.4米。
图3 RC梁配筋布置图
配筋:上部2C10、下部4C10,箍筋B6@100。
材料本构是爆炸仿真最容易做错的一环。混凝土用*MAT_RHT,三个极限面加P-α状态方程;钢筋用*MAT_PLASTIC_KINEMATIC,应变率按Cowper-Symonds取C=40、P=5;炸药是*MAT_HIGH_EXPLOSIVE_BURN配*EOS_JWL,空气用*MAT_NULL配线性多项式状态方程。
图4 ALE流固耦合有限元模型
欧拉域装空气和炸药,拉格朗日域是梁与支座。
耦合走ALE流固耦合,混凝土、空气和支座用SOLID164,钢筋用BEAM161。网格按2、5、10毫米三档做了收敛对比,误差都在15%以内,最终取5毫米。
判据是全文最有价值的部分:沿跨分段反推各段的弯曲变形,总变形减掉它即为剪切变形,比值定义为δ——δ≤5%判弯曲破坏,20%到35%为弯剪,超过50%即剪切破坏。对标结果:三个工况的损伤区与裂缝长度误差3.8%到9.76%,L1残余位移试验0.74厘米、仿真0.79厘米。
图5 L3工况试验破坏形态
L3工况试验破坏形态,迎爆面被打穿成空腔。
图6 L3工况仿真损伤云图
同工况仿真损伤云图,损伤集中在爆点正对区域。
图7 L1工况跨中残余位移对比
L1残余位移对比,试验0.74厘米、仿真0.79厘米。
三、关键结果与发现
偏心距的影响不单调:2千克装药下,位移从跨中的3.50厘米降到偏心0.9米处的1.37厘米,拉到1.05米又回到1.43厘米。偏心把一部分冲量换成了扭矩,爆点附近材料提前进塑性形成局部塑性铰,刚度软化后位移回升。位移小不等于安全。
图8 2kg药量下不同爆点位置的损伤与δ临界阈值
2千克装药下不同爆点的损伤分布与δ变化。
δ随偏心距呈V形:2千克时从0.3117降到0.0670再回升到0.2635,装药加大,拐点往跨中挪。药量影响更单调,跨中起爆时装药从2千克加到2.5千克,位移由3.50厘米涨到8.22厘米,δ升到0.4428,直接跨进类剪切。作者建议跨中两侧0.45到0.90米加密箍筋。
四、我们提供的服务
我们承接爆炸与冲击方向的仿真外包:ALE流固耦合与CONWEP方案选型,RHT、JWL参数标定,构件抗爆评估与加固比选,交付模型文件与解读报告。
如需相关服务,欢迎通过公众号“320科技工作室”联系我们。










评论0