Underwater anti-explosion mechanism and damage grade prediction of different corrugated steel-concrete slab composite structures
-
摘要: 为探究不同波纹钢-混凝土板复合结构的水下抗爆机理,采用光滑粒子流体动力学与有限单元(FEM-SPH)耦合方法模拟混凝土板在水下接触爆炸下的损伤过程,并与试验结果对比,以验证数值方法的有效性;采用FEM-SPH方法探究不同防护方案下墙面板毁伤过程及失效模式,揭示其水下防爆机理,并构建出墙面板损伤等级预测曲线。研究结果表明:模拟结果与试验结果较为吻合,验证了模拟方法的有效性;在含12 mm厚波纹钢的复合结构(T-12)、含75°夹角的波纹钢复合结构(A-75)和含厚70 mm波纹钢的复合结构(WH-70)三种防护方案下,墙面板的毁伤范围较未加固墙面板最大降幅分别为83.0%、81.6%和82.5%;预测曲线可以直观评估出炸药量和复合结构中波纹钢波高变化对墙面板损伤等级的影响。
-
关键词:
- 波纹钢-混凝土板复合结构 /
- 毁伤演化特征 /
- 水下抗爆机理 /
- 损伤等级预测
Abstract: In order to explore the underwater anti-explosion mechanism of different corrugated steel-concrete slab composite structures, the damage process of concrete slab under underwater contact explosion was simulated by smoothed particle hydrodynamics and finite element method (FEM-SPH), and the validity of the numerical method was verified by comparing with the experimental results. The FEM-SPH method was used to explore the damage process and failure mode of the wall panel under different protection schemes, to reveal the underwater explosion-proof mechanism, and to construct the prediction curve of the damage grade of the wall panel. The results show that the simulation results are in good agreement with the experimental results, which verifies the effectiveness of the simulation method. Under different protection schemes, the damage range of the wall panel with 12 mm thick corrugated steel composite structure (T-12), 75o angle corrugated steel composite structure (A-75) and 70 mm corrugated steel composite structure (WH-70) is 83%, 81.6% and 82.5% lower than that of the unreinforced wall panel, respectively. In the composite structure, the explosion shock wave propagates to the corrugated steel in the form of incident wave and then propagates in the structure in the form of transmitted wave and reflected wave. When the transmitted wave reaches the lower surface of the corrugated steel, part of the shock wave will continue to propagate to the wall panel, while the remaining shock wave is reflected to form reflected longitudinal wave and reflected transverse wave, which further attenuates the transmitted shock wave acting on the wall panel to achieve the effect of wave clipping and energy absorption. The prediction curve can directly evaluate the influence of explosive amount and wave height change of corrugated steel in composite structure on the damage grade of wall panel. -
建筑设施的大多数结构都采用复合材料或高性能材料进行加固。墙面板是目前较为常见的一种承载构件,但其在抵御爆炸荷载或冲击荷载时容易产生损坏而影响整体结构的安全性,因此对墙面板的防爆加固尤为重要。研究人员通过试验和数值模拟的方法探究了钢筋混凝土梁、板、柱的抗爆性能[1-3]。普通混凝土结构不足以承受爆炸荷载的冲击作用,因此学者们逐步开展混凝土结构加固方面的研究,如:杨谨鸿等[4]采用工程水泥基复合材料加固砌体填充墙并探究其抗爆性能;廖维张等[5]通过开展高强钢丝绳网片-聚合物砂浆加固混凝土板的爆炸试验,分析不同因素对混凝土板损伤程度的影响;Zhao等[6]对玄武岩纤维增强聚合物加固隧道进行抗爆试验以探究其承载能力。目前大多数研究主要聚焦于空气爆炸荷载下普通混凝土结构或经加固混凝土结构的损伤特性,而经过前人的研究发现,水下爆炸对结构的损伤较空气中爆炸更为严重[7-8],且水下爆炸下混凝土结构的损伤特性研究颇少。因此,本文将探究水下爆炸下波纹钢-混凝土板复合结构加固墙面板的抗爆性能及损伤机理。
如今,一些数值算法日益成熟,逐渐成为除试验外最有力的方法,如arbitrary Lagrangian-Eulerian (ALE)、光滑粒子流体动力学(smoothed particle hydrodynamics, SPH)和有限元-光滑粒子流体动力学耦合算法(FEM-SPH),这些方法在爆炸方面已得到广泛的应用,如:吕晋贤等[9]采用ALE算法探究高层钢筋混凝土框架结构在爆炸荷载下的抗倒塌性能;Elveli等[10]采用ALE法模拟弹道冲击薄钢板的过程;Lai等[11]采用SPH法探究侵彻和爆炸下超高性能混凝土和功能梯度胶凝复合材料的动态过程;杨建华等[12]采用FEM-SPH算法探究岩石爆破开裂与地应力的关系;Cui等[13]采用FEM-SPH法模拟土中浅埋爆炸的毁伤过程并分析爆坑与炸药量之间的关系。目前,已有研究发现FEM-SPH法能更好地模拟爆炸的损伤过程[14-15],但该方法用于分析水下接触爆炸下结构的损伤研究较少,仍需进一步研究。
本文中采用LS-DYNA中有限元-光滑粒子流体动力学耦合(FEM-SPH)算法模拟混凝土板在水下接触爆炸荷载下的损伤过程,并将数值结果与试验结果进行对比,以验证该算法的有效性;探究水下接触爆炸荷载下不同波纹钢-混凝土板复合结构对墙面板的毁伤演化过程、失效模式等抗爆性能的影响,揭示其水下防爆机理;依据墙面板的损伤模式建立损伤等级预测模型,预测不同防护方案下墙面板的损伤程度。
1. 爆炸模拟有效性验证
1.1 数值算法
典型的数值计算方法的基本原理如图1所示。ALE法在分析结构运动上具备拉格朗日方法的特点,在网格上具备欧拉方法的优点,可以理解为两层网格重叠在一起,空间网格可以任意运动且分析结构也可随意运动,但在解决大变形问题时容易出现畸变现象。SPH是一种无网格算法,其本身无需使用单元,而是使用固定质量的质点,因此在水下爆炸、高速撞击等材料极度变形的问题上不会产生网格畸变,可用来模拟结构破碎、成坑等现象,但太多粒子会加大计算量。FEM-SPH法则是在大变形区域设置为粒子,而在小变形区域设置为有限元单元,该方法综合了两种方法的优点,既保证计算的准确性,也可以提高工作效率。
1.2 试验概况
本文以文献[16]中水下接触爆炸试验为背景,采用FEM-SPH法模拟抗压强度为28.2 MPa的混凝土板在6 g炸药下的爆炸试验过程。试验选取60 mm厚的钢筋混凝土板(尺寸为500 mm × 500 mm),板的配筋采用直径为6 mm且间距为100 mm的钢筋,其保护层厚度为20 mm。钢筋的屈服强度为600 MPa,杨氏模量为200 GPa。混凝土板钢筋布置图如图2(a)所示。在接触爆炸试验中,混凝土板上表面距自由水面0.5 m且放置在钢架上,炸药位于板的中心位置并采用雷管引爆,试验模型如图2(b)所示。
结合现场试验,采用FEM-SPH法建立水下接触爆炸下钢筋混凝土板数值模型如图3所示。模型中,炸药、炸药附近的混凝土和水体均采用SPH粒子,远处混凝土和水体采用有限元单元,钢筋采用梁单元。空气、混凝土板、钢筋和炸药的网格尺寸均为5 mm,共有3 496 156个实体单元和15 444个SPH粒子。混凝土和钢筋之间采用CONSTRAINED_BEAM_IN_SOILD耦合。在模型对称面施加对称边界,在水体外围施加无反射边界条件。
1.3 材料模型及参数
1.3.1 混凝土
混凝土采用*MAT_CONCRETE_DAMAGE_REL3模型[17],并通过*MAT_ADD_EROSRION关键字定义失效准则,最大主应变取为0.01[18]。该模型只需将密度、泊松比、单轴抗压强度和压力-应变率提高系数关系曲线输入,便可自动生成材料模型的其他参数和状态方程的参数。无侧限抗压强度、密度和泊松比分别为30 MPa、2.4 g/cm3和0.19。混凝土材料的应变率效应可用动力增大系数表示。
混凝土抗压和抗拉强度动力增大系数分别为:
ηc=fcfcs={(˙εc/˙εcs)1.026α˙εc≤30s−1γ(˙εc/˙εcs)1/3˙εc>30s−1 (1) ηt=ftfts={(˙εt/˙εts)δ˙εt≤1s−1β(˙εt/˙εts)1/3˙εt>1s−1 (2) 式中:fc和fcs分别为动、静态抗压强度;
˙εc 和˙εcs 分别为动静态压缩应变率,˙εcs=300s−1 ;α和γ为压缩应变率系数,logγ=6.156α−2 ,α=(5+9fcs/10)−1 ;ft和fts分别为动、静态抗拉强度;˙εt 和˙εts 分别为动、静态拉伸应变率;β和δ为拉伸应变率系数,logβ=6δ−2 ,δ=1/(1+8fcs/10) 。1.3.2 炸药TNT
炸药TNT采用*MAT_HIGH_EXPLOSIVE_BURN模型和*EOS_JWL状态方程描述[19],模型参数如表1所示。
表 1 炸药材料参数Table 1. Material parameters of the explosiveρ1/(g·cm−3) A1/MPa B1/MPa R1 R2 ω1 1.63 373.77 3.75 4.15 0.9 0.35 pd=A1(1−ω1R1V)e−R1V+B1(1−ω1R2V)e−R2V+ω1E0V (3) 式中:
pd 为爆轰压力;A1、B1为材料压力参数;R1 、R2 、ω1 为状态方程量纲系数;E0 为体积内能。1.3.3 水体
水体采用*MAT_NULL模型和*EOS_GRUNEISEN状态方程[20],压缩状态和膨胀状态的压力为
p={ρ0c2μ[1+(1−γ0/2)μ−aμ2/2][1−(S1−1)μ−S2μ2/(μ+1)−S3μ2/(μ+1)2]2+(γ0+aμ)Eρ0c2μ+(γ0+aμ)E (4) 式中:p为水中的压力,ρ0为初始密度,E为初始体积的热力学能,μ为相对体积,c为vs-vp曲线的截距,S1、S2、S3为vs-vp曲线斜率系数,γ0为系数,a为一阶体积修正。模型参数见表2所示。
表 2 水体材料参数Table 2. Material parameters of waterρ2/(g·cm−3) c/(m·s−1) S1 S2 S3 γ0 a 1.0 2417 1.41 0 0 1.0 0 1.3.4 波纹钢及钢筋
钢筋采用能描述高应变率和高温下材料大变形及屈服应力变化的Johnson-Cook模型和Grüneisen状态方程[21],模型参数见表3所示。
表 3 钢筋材料参数Table 3. Material parameters corrugated steelρ3/(g·cm−3) A2/GPa B2/GPa N C m Tm 7.85 0.345 0.336 0.42 0.026 1.4 1720 注:A2、B2、N为参考应变率˙ε0和参考温度Troom下的材料初始屈服应力、应变硬化模量和硬化指数,C为材料应变率强化参数,m为材料热软化参数,Troom为室温,Tmelt为熔点。 1.4 耦合方法有效性验证
水下接触爆炸荷载下混凝土板的试验和数值模拟结果对比如图4所示。由图4可知,在爆炸冲击荷载作用下,混凝土板的迎爆面出现一条贯穿整个板的裂缝,且伴随有一些微小裂缝;在板的背爆面出现震塌剥落现象,且出现多条径向裂缝;对比试验结果与数值结果可知,二者的最终毁伤模式基本一致,产生误差的可能原因是:(1) 数值模拟中混凝土被视作理想均质材料,但现场试验中混凝土则为非均质材料;(2) 仿真模拟中炸药布置定位准确,但现场试验中炸药位置可能会出现人为偏差等不可控因素。综上可知,该仿真模拟方法能够有效表征水下接触爆炸荷载下混凝土板的毁伤过程及失效模式。
2. 波纹钢-混凝土板复合结构防爆模型及防护方案
2.1 数值模型建立
采用FEM-SPH法建立水下接触爆炸模型,如图5所示。模型中墙面板的尺寸为1000 mm ×500 mm × 60 mm,其内部单层钢筋直径为6 mm,纵向钢筋间距为77 mm,横向钢筋间距为107 mm,保护层厚度为20 mm,具体如图6所示。在墙面板上部采用波纹钢-混凝土板复合结构进行加固,其中波纹钢位于混凝土中间,复合结构的尺寸为1000 mm ×500 mm × 100 mm。波纹钢板、墙面板和离炸药中心远的混凝土板及水体采用实体单元,炸药、离炸药中心近的混凝土板(150 mm ×150 mm × 50 mm)和炸药周围200 mm的水体采用SPH粒子;混凝土板、波纹钢、钢筋和炸药网格尺寸均为5 mm,水体网格尺寸随距炸药中心的距离呈梯度分布(5、10、20 mm);模型中共有1 048 410个实体单元和142 790个SPH粒子。混凝土SPH粒子和波纹钢之间用关键字*CONTACT_AUTOMATIC_NODES_TO_SURFACE进行耦合,混凝土和钢筋用关键字*CONSTRAINED_BEAM_IN_SOILD耦合。此外,为更好地接近于实际工程,设定TNT炸药当量为120 g,水体外围设置为自由边界,混凝土板两端为固定约束。该模型中混凝土、炸药、水体和波纹钢及钢筋均采用1.3节中的材料模型和参数;而波纹钢采用Johnson-Cook模型和Grüneisen状态方程定义其材料参数,其强度值为600 MPa。
2.2 波纹钢-混凝土板复合结构防护方案设计
为探究不同波纹钢-混凝土板复合结构的水下抗爆机理,选取12种波纹钢形状参数(厚度、夹角和波高),其中T-3、A-45和WH-50为对照组,防护方案中波纹钢截面如图7所示,其中T、A、WH分别代表波纹钢厚度、夹角和波高,例如:T-3代表波纹钢厚度为3 mm,A-30代表波纹夹角为30°,WH-10代表波高为10 mm。
3. 毁伤演化特征及失效机理
3.1 毁伤过程及机理
为更好探究水下接触爆炸荷载下波纹钢-混凝土板复合结构加固墙面板的毁伤过程,分别对整体结构、波纹钢和墙面板进行分析,限于篇幅仅展示T-3防护方案(对照组)的毁伤过程。
如图8所示,炸药在水下爆炸后,结构迎爆面出现损伤,而背爆面因冲击波尚未到达而无损伤发生;在0.08 ms时,压缩波传播至结构背爆面,一部分冲击波穿透结构后仍在水中传播,其余冲击波在结构背爆面发生反射而导致背爆面发生拉伸破坏,经反射的冲击波向结构迎爆面方向传播致使迎爆面的毁伤范围进一步扩大;当冲击波传播至背爆面,即0.12 ms时,出现震塌现象;而0.25 ms时迎爆面爆坑和背爆面震塌范围的最大直径分别为324和221 mm。
如图9所示,在0.03 ms时波纹钢跨中位置开始受力,且波纹钢的受力范围随着冲击波传播时间增加而逐渐增大;在0.25 ms时,波纹钢整体受力且跨中位置变形达到最大。
如图10所示,墙面板迎爆面在0.06 ms时开始出现损伤,而背爆面无变化;直至0.08 ms时,墙面板背爆面因冲击波发生反射形成的拉伸波而开始产生拉伸破坏,迎爆面冲击压力也随反射波作用而逐渐增大;在0.13 ms时刻,由于波纹钢和混凝土的相互挤压,导致混凝土材料达到失效准则而形成裂缝;随着冲击波传播时间持续增加,墙面板的损伤在0.25 ms时达到最大,其中迎爆面呈现出多条裂缝,背爆面爆坑最大直径和最大宽度分别为221和180 mm。
图11显示了波纹钢-混凝土板复合结构加固墙面板的毁伤机理。在波纹钢-混凝土板复合结构中,爆炸冲击波以入射波的形式传播至波纹钢后,由于波纹钢和混凝土二者的波阻抗不同,在分界面出现两种波的形式:一种以透射波的形式继续向下传播,另一种则在分界面处以反射波的形式在结构内传播;之后,当透射波到达波纹钢的下表面时,一部分冲击波继续向墙面板传播,而其余冲击波则再经反射最终形成反射纵波和反射横波,进一步衰减了作用于墙面板的透射冲击波。由此可知,爆炸冲击波因复合结构中波纹钢结构形式的存在而增加了冲击波在结构内的反射次数,进而达到分散和消耗传播至被防护结构的冲击波的作用,最终达到削波效果。炸药在水中起爆后,冲击波直接作用在结构的迎爆面,导致迎爆面出现轻微压缩破坏;随着冲击波传播,结构迎爆面的受压破坏区域进一步扩大,且结构背爆面出现轻微剥落现象;当传入结构的冲击波到达结构下表面时,一部分被反射形成拉伸波,当拉应力超过混凝土的抗拉强度,结构背爆面出现大范围震塌,且爆轰产物与结构上表面的相互作用会在水体和结构的交界面处呈现空化现象。
3.2 毁伤破坏模式
为更好探究不同波纹钢-混凝土板复合结构的防护效果和结构在遭受爆炸冲击荷载后的安全性,本文将爆坑的最大直径作为测量标准。未加固墙面板毁伤模式如图12所示。墙面板迎爆面出现严重损伤,形成一个大爆坑;在背爆面出现严重的震塌现象,毁伤区的最大直径和最大宽度分别为330和320 mm,且墙面板整体发生贯穿破坏。
图13~15给出了不同厚度波纹钢复合结构防护方案下各部分的毁伤模式。可以看出,与未加固墙面板的毁伤模式相比,T-12防护方案下墙面板的毁伤范围最大降幅为83%,说明增加复合结构能够有效减小墙面板的损伤。不同防护方案下结构迎爆面和背爆面的毁伤特征类似,但损伤程度存在较大差异。整体结构、波纹钢和墙面板的毁伤随着波纹钢厚度的增大而逐渐减小。
由图15可知,T-6防护方案下墙面板迎爆面损伤区域由板中心向外开展且有裂缝生成;背爆面出现环向和径向破坏,表现为震塌剥落现象,墙面板剥落区域最大直径和最大宽度分别为210和151 mm。T-9防护方案下墙面板迎爆面裂缝的数量和长度明显减小,背爆面剥落范围也有所减小且剥落区的最大直径和宽度分别为191和130 mm。T-12防护方案下墙面板迎爆面仅有两条裂缝,背爆面剥落区域最大直径和最大宽度分别为180和100 mm。可见,墙面板迎爆面的损伤区和背爆面的剥落区随波纹钢厚度的增大而逐渐减小,且T-12防护方案下墙面板的毁伤范围较其余三种防护方案最小,说明水下爆炸荷载条件下T-12防护方案具有更好的抗爆防护效果。
图16~18给出了不同夹角波纹钢的复合结构防护方案下各部分毁伤模式。A-75防护方案下墙面板的毁伤范围较未加固墙面板最大降幅为81.6%,说明增加波纹钢夹角能够显著减小墙面板的损伤。整体结构、波纹钢和墙面板的毁伤程度均随着波纹钢夹角的增大而逐渐减小。
由图18可知,A-30防护方案下墙面板迎爆面中部位置呈现大范围损伤,且裂缝几乎贯穿墙面板的整个宽度方向;墙面板背爆面呈现震塌剥落现象,整个损伤区最大直径和最大宽度分别为285和154 mm。A-60防护方案下墙面板的迎爆面损伤区域和裂缝的数量和长度明显减小,背爆面混凝土剥落范围也有所减小,剥落区最大直径和最大宽度分别为211和131 mm。A-75防护方案下,由于波纹钢夹角的增大而使波纹钢受到爆炸冲击波作用与墙面板相互挤压,最终在墙面板迎爆面出现大量应力集中现象;背爆面剥落区域明显减小且最大直径和最大宽度分别为185和105 mm。随着波纹钢夹角的增大,墙面板迎爆面受到的损伤和背爆面的剥落区均逐渐减小,且A-75防护方案下墙面板的损伤最小,说明水下爆炸荷载条件下A-75防护方案抗爆效果最好。
图19~21给出了不同波高波纹钢复合结构防护方案下各部分毁伤模式。可以看出:与未加固墙面板的毁伤模式相比,WH-70防护方案下墙面板的毁伤降幅为82.5%,说明增加复合结构能够有效减小墙面板的损伤,且整体结构、波纹钢和墙面板的毁伤随着波纹钢波高的增大而减小。
如图21所示,WH-10防护方案下墙面板的迎爆面出现大面积损伤,而裂缝现象较少,主要原因是波纹钢波高较小且接近于平钢板,爆炸荷载下波纹钢与墙面板的挤压作用较小,导致失效的单元较少;墙面板的背爆面呈现大范围的震塌,损伤区最大直径和最大宽度分别为281和251 mm。WH-30防护方案下墙面板的迎爆面损伤区域有所减小,而裂缝的数量和长度增加,主要是因波纹钢折角处对墙面板的挤压导致接触处的单元达到失效形成裂缝;墙面板背爆面剥落范围也有所减小,剥落区的最大直径和最大宽度分别为241和220 mm。WH-70防护方案下墙面板的迎爆面裂缝的长度和数量均减小;在板的背爆面剥落区域明显减小,最大直径和最大宽度分别为181和105 mm。随着波纹钢波高的增大,墙面板迎爆面受到的损伤和背爆面的剥落区均逐渐减小,且WH-70防护方案下墙面板的损伤最小,说明增加波纹钢波高在一定程度上会减小被防护结构的损伤。
为更直观地对比不同防护方案下墙面板背爆面的毁伤模式,将剥落区的最大直径和最大宽度以表格的形式给出,详见表4。由表4可知,T-12防护方案中各部分的损伤较小,建议选用T-12防护方案进行水下结构防爆加固设计。
表 4 墙面板背爆面剥落区Table 4. Spalling area of back surface of wall panel方案 最大直径/mm 最大宽度/mm 方案 最大直径/mm 最大宽度/mm 方案 最大直径/mm 最大宽度/mm T-3 221 180 A-30 285 154 WH-10 281 251 T-6 210 151 A-45 221 180 WH-30 241 220 T-9 191 130 A-60 211 131 WH-50 221 180 T-12 180 100 A-75 185 105 WH-70 181 105 注:未加固结构剥落区的最大直径和最大宽度分别为330和320 mm。 3.3 墙面板的钢筋轴力
不同防护方案下钢筋轴力如图22~25所示。无加固结构钢筋中部位置产生较大受力,最大轴力为15.5 kN。随着波纹钢厚度、夹角和波高的增加,钢筋轴力均一致下降,说明增加复合结构防护层能够有效减小墙面板中钢筋受力,T-12、A-75和WH-70防护方案下钢筋轴力较无加固结构的最大降幅分别为24.5%、22.6%和14.8%。综上可知,T-12防护方案下钢筋轴力相对较小,建议选取T-12防护方案进行水下结构的防爆加固。
4. 损伤等级预测
4.1 墙面板损伤等级
水下接触爆炸荷载下被防护结构墙面板迎爆面损伤主要为横向裂缝,而背爆面由于拉伸波的影响形成爆坑。根据Zhou等[22]和Li等[23]对结构损伤程度的划分,并结合本研究不同防护方案下墙面板的毁伤特征,将水下接触爆炸荷载下不同波纹钢-混凝土板复合结构加固墙面板的损伤程度大致划分为三类(图26):(1) 轻度破坏,板没有破坏或出现一些可见的裂缝, 出现浅层层裂到破坏深度为1/3板厚;(2) 中度破坏,破坏深度介于1/3板厚到2/3板厚之间;(3) 重度破坏,破坏深度介于1/3板厚到2/3板厚之间。
图27~29给出了不同防护方案下墙面板的破坏深度。如图27所示,墙面板最大破坏深度分别约为43、38、35和29 mm,说明波纹钢的厚度增加会有效降低墙面板的破坏程度。依据图26的划分标准,T-3防护方案下墙面板的最大破坏深度约为板厚的71.7%,表明该防护方案下墙面板的破坏程度为重度破坏,T-6、T-9和T-12防护方案下墙面板的最大破坏深度分别约为板厚的63.3%、58.3%和48.3%,表明这三种防护方案下墙面板的破坏程度均为中度破坏。如图28所示,墙面板的破坏深度随波纹钢夹角的增大而减小,根据图26所划分的损伤等级评判标准,A-30和A-45防护方案下墙面板的破坏深度大于2/3板,表现为重度破坏。A-60和A-75防护方案下墙面板的破坏深度介于1/3~2/3之间,表现为中度破坏。如图29所示,不同防护方案下墙面板的破坏深度随波纹钢波高的增大而减小,根据图26的划分标准,WH-10、WH-30和WH-50防护方案下墙面板的破坏深度均大于2/3板,表现为重度破坏;WH-70防护方案下墙面板破坏深度介于1/3~2/3之间,表现为中度破坏。
4.2 墙面板损伤等级预测
为了得出不同防护方案下墙面板损伤等级与波纹钢厚度的关系,分别模拟出30、60、90、120、150、180和210 g等7种TNT当量的炸药下墙面板的毁伤模式,并将其破坏深度进行分类,最终依据损伤等级拟合出墙面板的损伤等级预测曲线,如图30所示。预测曲线1为轻度破坏和中度破坏的分界线,其拟合公式为:
h=13.103lnw−50.4483mm≤h≤12mm (5) 预测曲线2为中度破坏与重度破坏的分界线,其拟合公式为:
h=15.509lnw−70.2943mm≤h≤12mm (6) 式中:h为波纹钢厚度,mm;w为炸药的TNT当量,g。
为得出墙面板的损伤等级与波纹钢夹角之间的关系,模拟7种炸药当量下墙面板的损伤过程,最终得出不同防护方案下墙面板的破坏深度并将其进行划分。在此基础上,拟合出墙面板的损伤等级预测曲线,如图31所示。预测曲线1为轻度破坏与中度破坏的分界线:
θ=41.403lnw−110.13430∘≤θ≤75∘ (7) 预测曲线2为中度破坏与重度破坏的分界线:
θ=50.07lnw−194.72930∘≤θ≤75∘ (8) 式中:
为波纹钢夹角,(°)。 为得出波纹钢-混凝土板复合结构加固下墙面板损伤等级与波纹钢波高之间的关系,将数值结果按图26中墙面板的损伤等级进行分类,拟合出墙面板的损伤等级预测曲线,如图32所示。轻度破坏与中度破坏的分界线为预测曲线1:
H=54.947lnw−176.371,10mm≤H≤90mm (9) 中度破坏与重度破坏的分界线为预测曲线2:
H=57.158lnw−225.886,10mm≤H≤90mm (10) 式中:H为波纹钢波高,mm。
通过模拟结果拟合出墙面板损伤程度的预测曲线,可以清晰看出波纹钢波高和炸药量对被防护结构墙面板损伤程度的影响。根据预测曲线,当炸药量一定时,为保证被防护墙面板的损伤程度达到最小,可优先选取恰当波高的波纹钢。
5. 结 论
采用LS-DYNA软件,分析波纹钢形状参数变化对波纹钢-混凝土板复合结构水下抗爆机理及毁伤程度的影响,研究结论如下:
(1) FEM-SPH方法能够有效模拟混凝土板水下接触爆炸过程;
(2) T-12、A-75和WH-70防护方案下墙面板的毁伤范围较未加固墙面板最大降幅分别为83%、81.6%和82.5%;结构的毁伤和钢筋轴力随波纹钢厚度(夹角、波高)增大而逐渐减小,且未加固墙面板较T-12、A-75和WH-70防护方案下钢筋轴力最大降幅分别为24.5%、22.6%和14.8%;对比分析不同防护方案可知,建议选用T-12防护方案;
(3) 依据损伤等级评判标准构建出墙面板损伤等级的预测曲线,可以清晰评估出炸药量和波纹钢厚度、夹角和波高的变化对被防护结构墙面板损伤等级的影响。
-
表 1 炸药材料参数
Table 1. Material parameters of the explosive
ρ1/(g·cm−3) A1/MPa B1/MPa R1 R2 ω1 1.63 373.77 3.75 4.15 0.9 0.35 表 2 水体材料参数
Table 2. Material parameters of water
ρ2/(g·cm−3) c/(m·s−1) S1 S2 S3 γ0 a 1.0 2417 1.41 0 0 1.0 0 表 3 钢筋材料参数
Table 3. Material parameters corrugated steel
ρ3/(g·cm−3) A2/GPa B2/GPa N C m Tm 7.85 0.345 0.336 0.42 0.026 1.4 1720 注:A2、B2、N为参考应变率˙ε0和参考温度Troom下的材料初始屈服应力、应变硬化模量和硬化指数,C为材料应变率强化参数,m为材料热软化参数,Troom为室温,Tmelt为熔点。 表 4 墙面板背爆面剥落区
Table 4. Spalling area of back surface of wall panel
方案 最大直径/mm 最大宽度/mm 方案 最大直径/mm 最大宽度/mm 方案 最大直径/mm 最大宽度/mm T-3 221 180 A-30 285 154 WH-10 281 251 T-6 210 151 A-45 221 180 WH-30 241 220 T-9 191 130 A-60 211 131 WH-50 221 180 T-12 180 100 A-75 185 105 WH-70 181 105 注:未加固结构剥落区的最大直径和最大宽度分别为330和320 mm。 -
[1] 许迎亮, 刘彦, 闫俊伯, 等. 双装药同步爆炸钢筋混凝土梁毁伤效应 [J]. 兵工学报, 2023, 44(12): 3719–3732. DOI: 10.12382/bgxb.2023.0286.XU Y L, LIU Y, YAN J B, et al. The damage effect of reinforced concrete beams subjected to synchronous explosion of double charges [J]. Acta Armamentarii, 2023, 44(12): 3719–3732. DOI: 10.12382/bgxb.2023.0286. [2] 赵小华, 刘树参, 方宏远, 等. 水下接触爆炸下高聚物层对钢筋混凝土板的防护效果 [J]. 爆炸与冲击, 2023, 43(12): 125102. DOI: 10.11883/bzycj-2023-0033.ZHAO X H, LIU S S, FANG H Y, et al. Protective effect of polymer layer on reinforced concrete slabs under an underwater contact explosion [J]. Explosion and Shock Waves, 2023, 43(12): 125102. DOI: 10.11883/bzycj-2023-0033. [3] 杨坤, 张玮, 李营, 等. 水下爆炸作用下复合材料圆柱壳结构失效模式分析 [J]. 中国舰船研究, 2023, 18(2): 55–63. DOI: 10.19693/j.issn.1673-3185.02835.YANG K, ZHANG W, LI Y, et al. Failure mode analysis of composite cylindrical shell structure under underwater explosion [J]. Chinese Journal of Ship Research, 2023, 18(2): 55–63. DOI: 10.19693/j.issn.1673-3185.02835. [4] 杨谨鸿, 李秀地, 张波, 等. 近距离爆炸下ECC涂层加固砌体填充墙抗爆性能研究 [J]. 振动与冲击, 2023, 42(16): 10–18. DOI: 10.13465/j.cnki.jvs.2023.16.002.YANG J H, LI X D, ZHANG B, et al. Anti-explosion performance of ECC coating for strengthening masonry infilled wall under close-in explosion [J]. Journal of Vibration and Shock, 2023, 42(16): 10–18. DOI: 10.13465/j.cnki.jvs.2023.16.002. [5] 廖维张, 刘锴鑫, 张春磊, 等. 高强钢丝绳网片-聚合物砂浆加固RC板抗爆性能试验研究 [J]. 振动与冲击, 2021, 40(2): 235–242. DOI: 10.13465/j.cnki.jvs.2021.02.032.LIAO W Z, LIU K X, ZHANG C L, et al. Blast resistant performance of reinforced concrete slabs strengthened with high strength steel wire mesh and polymer mortar [J]. Journal of Vibration and Shock, 2021, 40(2): 235–242. DOI: 10.13465/j.cnki.jvs.2021.02.032. [6] ZHAO C J, TANG Z X, WANG P, et al. Blast responses of shallow-buried prefabricated modular concrete tunnels reinforced by BFRP-steel bars [J]. Underground Space, 2022, 7(2): 184–198. DOI: 10.1016/J.UNDSP.2021.07.004. [7] YANG G D, WANG G H, LU W B, et al. Cross-section shape effects on anti-knock performance of RC columns subjected to air and underwater explosions [J]. Ocean Engineering, 2019, 181: 252–266. DOI: 10.1016/j.oceaneng.2019.04.031. [8] ZHAO X H, WANG G H, LU W B, et al. Damage features of RC slabs subjected to air and underwater contact explosions [J]. Ocean Engineering, 2018, 147: 531–545. DOI: 10.1016/j.oceaneng.2017.11.007. [9] 吕晋贤, 吴昊, 方秦. 爆炸作用下高层框架结构倒塌分析与设计建议 [J]. 建筑结构学报, 2023, 44(11): 114–128. DOI: 10.14006/j.jzjgxb.2022.0454.LÜ J X, WU H, FANG Q. Collapse analysis and design recommendations of high-rise frame structures under blast loadings [J]. Journal of Building Structures, 2023, 44(11): 114–128. DOI: 10.14006/j.jzjgxb.2022.0454. [10] ELVELI B S, VESTRUM O, HAUGE K O, et al. Thin steel plates exposed to combined ballistic impact and partially confined airblast loading [J]. Engineering Failure Analysis, 2023, 144: 106943. DOI: 10.1016/J.ENGFAILANAL.2022.106943. [11] LAI J Z, LI H J, YIN X X, et al. Properties and simulation of UHPC and FGCC subjected to the coupling of penetration and explosion [J]. Journal of Materials in Civil Engineering, 2021, 33(6): 04021109. DOI: 10.1061/(ASCE)MT.1943-5533.0003665. [12] 杨建华, 吴泽南, 姚池, 等. 地应力对岩石爆破开裂及爆炸地震波的影响研究 [J]. 振动与冲击, 2020, 39(13): 64–70,90. DOI: 10.13465/j.cnki.jvs.2020.13.010.YANG J H, WU Z N, YAO C, et al. Influences of in-situ stress on blast-induced rock fracture and seismic waves [J]. Journal of Vibration and Shock, 2020, 39(13): 64–70,90. DOI: 10.13465/j.cnki.jvs.2020.13.010. [13] CUI Y, LI Z J, FANG J, et al. Crater effects of shallow burial explosions in soil based on SPH-FEM analysis [J]. Frontiers in Earth Science, 2023, 10: 1114178. DOI: 10.3389/FEART.2022.1114178. [14] PAN G C, SU H, LI X X, et al. Coupled FEM-SPH simulation of the protective properties for metal/ceramic composite armor [J]. International Journal of Lightweight Materials and Manufacture, 2023, 6(4): 543–551. DOI: 10.1016/J.IJLMM.2023.05.007. [15] WANG Z L, HUANG Y P, LI S Y, et al. SPH-FEM coupling simulation of rock blast damage based on the determination and optimization of the RHT model parameters [J]. IOP Conference Series: Earth and Environmental Science, 2020, 570(4): 042035. DOI: 10.1088/1755-1315/570/4/042035. [16] ZHAO X H, WANG G H, LU W B, et al. Experimental investigation of RC slabs under air and underwater contact explosions [J]. European Journal of Environmental and Civil Engineering, 2021, 25(1): 190–204. DOI: 10.1080/19648189.2018.1528892. [17] ZHAO C F, LU X, WANG Q, et al. Experimental and numerical investigation of steel-concrete (SC) slabs under contact blast loading [J]. Engineering Structures, 2019, 196: 109337. DOI: 10.1016/j.engstruct.2019.109337. [18] ZHAO C F, HE K C, ZHI L H, et al. Blast behavior of steel-concrete-steel sandwich panel: experiment and numerical simulation [J]. Engineering Structures, 2021, 246: 112998. DOI: 10.1016/J.ENGSTRUCT.2021.112998. [19] XU Q, CHEN J Y, LI J, et al. Numerical study on antiknock measures of concrete gravity dam bearing underwater contact blast loading [J]. Journal of Renewable and Sustainable Energy, 2018, 10(1): 014101. DOI: 10.1063/1.4986330. [20] HAI L, REN X D. Computational investigation on damage of reinforced concrete slab subjected to underwater explosion [J]. Ocean Engineering, 2020, 195: 106671. DOI: 10.1016/j.oceaneng.2019.106671. [21] 李超. 柱面网壳结构在内爆炸下的失效机理和防爆方法 [D]. 泉州: 华侨大学, 2016.LI C. Study on failure mechanism and explosion-proof method for cylindrical shell under internal explosion [D]. Quanzhou: Huaqiao University, 2016. [22] ZHOU X Q, HAO H. Numerical prediction of reinforced concrete exterior wall response to blast loading [J]. Advances in Structural Engineering, 2008, 11(4): 355–367. DOI: 10.1260/136943308785836826. [23] LI J, HAO H. Numerical study of concrete spall damage to blast loads [J]. International Journal of Impact Engineering, 2014, 68: 41–55. DOI: 10.1016/j.ijimpeng.2014.02.001. -