摘要
金属结构在冲击载荷作用下,最终的变形方向与加载的方向相反,这一现象被称为反直观行为,试验发现背爆面柔性聚脲防护混凝土靶板在75 g TNT爆炸荷载下也发生了反直观行为。为了研究背爆面柔性聚脲防护混凝土靶板在爆炸荷载下的反直观行为,采用ANSYS/LSDYNA有限元软件,建立了背爆面柔性聚脲防护混凝土靶板在爆炸荷载作用下的有限元模型。利用有限元模型分析背爆面柔性聚脲防护混凝土靶板的动态响应规律,从能量的角度研究靶板反直观行为的发生机理。以靶板的中心点位移和挠度为指标,参数化分析了炸药药量和聚脲涂层厚度对靶板反直观行为的影响规律。结果表明:靶板的反直观行为是由聚脲涂层的应变能释放、混凝土损伤破坏的能量耗散及两种材料间能量的相互转化这三种因素共同作用导致的;炸药药量是反直观行为能否发生的关键因素,当药量较低或较高时,靶板都无法发生反直观行为;涂层厚度在2~8 mm范围内,靶板均出现了反直观行为,其弯曲程度随着涂层厚度的增加先增大后减小。
图文摘要
爆炸荷载对土木工程结构的破坏远大于一般静态荷载和低速冲击荷载,土木工程结构中混凝土占比较大,其在爆炸荷载作用下将产生大量的高速破片,造成更为严重的损
国内外针对聚脲防护结构在爆炸荷载下的响应开展了大量的研究。孙鹏飞
由于聚脲分子结构的差异性,不同分子结构的聚脲所表现出的力学性能有明显区别,因而在爆炸荷载等高速荷载作用下的响应特性不同,导致最终的防护效果也不同。Ji
反直观行为是结构动态响应的一个特殊的现象,研究结构的反直观行为可为结构的动力学响应机理带来新的认识。1985年,Symonds
近期课题组在试验中发现,背爆面F‑T126柔性聚脲(以下简称F‑T126)防护混凝土靶板在爆炸荷载下也出现了反直观行为。为研究爆炸荷载作用下背爆面柔性聚脲涂层防护混凝土靶板(以下简称涂层混凝土靶板)的反直观行为发生机理,本研究建立了与试验工况相同的有限元模型。通过数值模拟,分析炸药药量和涂层厚度对涂层混凝土靶板反直观行为的影响规律。研究爆炸荷载作用下涂层混凝土靶板反直观行为的发生机理,可为抗爆复合结构的设计带来新的思路。
试验共制作了两块混凝土靶板,其长宽为300 mm×300 mm,厚度为30 mm,混凝土强度等级为C40。一块混凝土靶板的一侧喷涂聚脲,聚脲涂层尺寸为300 mm×300 mm×4 mm。另一块作为对照组,无聚脲涂层防护。
爆炸试验采用非接触爆炸的形式,试验详细情况如

a. test setup image

b. cross‑section schematic of test set
图1 试验现场布置
Fig.1 Test site layout
无涂层混凝土靶板在爆炸荷载作用下的试验结果如

a. concrete fragmentation after blast load

b. retrieved concrete fragments
图2 无涂层防护靶板试验照片
Fig.2 Test photograph of uncoated concrete target plate
涂层混凝土靶板的爆炸损伤情况如

a. blast surface

b. back‑blast surface

c. side surface

d. cross‑section bending deformation diagram
图3 涂层混凝土靶板的试验照片
Fig.3 Test photograph of coated concrete target plate
背爆面F‑T126的损伤如
涂层混凝土靶板侧面有明显开裂,裂缝多贯穿靶板侧面(
本研究使用有限元软件ANSYS/LS‑DYNA建立与爆炸试验尺寸相同的有限元模型,整体模型如

图4 有限元模型示意图
Fig.4 Diagram of finite element model
涂层和基材的粘结强度对结构的防护性能有重要影
(1) |
式中,σn、σs分别为接触面上的法向应力和剪切应力,MPa;NFLS和SFLS表示界面的法向强度和剪切强度。NFLS和SFLS的数值由厂家提供,NFLS和SFLS的数值分别为10 MPa和4 MPa。
混凝土的材料模型采用*MAT_RHT模型,该模型用于模拟岩石、混凝土等脆性材料在爆炸冲击荷载作用下的力学行为,可模拟混凝土靶板的裂纹损伤分布和靶后崩落等现象,材料参数见文献[
F‑T126的材料模型选用*MAT_PLECEWISE_LINEAR_PLASTICITY,该模型能够有效模拟聚脲的动态力学性能。采用MZ‑4000D1万能试验机,按照GB/T 528-2009《橡胶、硫化或热塑性橡胶拉伸应力应变性能的测定》对F‑T126进行准静态力学性能测试。在加载速度为500 mm·mi

图5 不同速度下F‑T126的拉伸应力应变曲线
Fig.5 Tensile stress‑strain curves of F‑T126 flexible polyurea with different loading speeds
基于
ρ / g·c | E / GPa | ν | σ0 / GPa | Et / GPa | FAIL |
---|---|---|---|---|---|
0.96 | 0.22 | 0.4 | 0.00778 | 0.00231 | 1.3 |
Note: ρ is density. E is elasticity modulus. ν is poisson ratio. σ0 is yield strength. Et is tangent modulus. FAIL is failure strain.
使用关键字*EOS_LINEAR_POLYNOMIAL定义空气的状态方程通过
ρ / g·c | C0‑C3 | C4 | C5 | e / J·c |
---|---|---|---|---|
0.00129 | 0 | 0.4 | 0.4 |
2.5×1 |
Note: ρ is density. C0‑C5 is polynomial equation coefficient. e is internal energy per unit volume.
(2) |
式中,e为单位体积能量,J·c
炸药为TNT装药,其本构模型采用*MAT_HIGH_EXPLOSIVE_BURN表述,使用*EOS_JWL来描述其状态方程,该方程可描述炸药爆炸产物压力,其表达式如
(3) |
式中,p是爆轰产物的压力,GPa;V为相对体积;A、B为常数,GPa;R1、R2、ω为无量纲常数。其中方程右端第一项在高压段起主要作用、第二项在中压段起主要作用,第三项代表低压段。
ρ / g·c | D / m· | PCJ / GPa | A / GPa | B / GPa |
---|---|---|---|---|
1.58 | 6880 | 19.4 | 307 | 3.898 |
Note: D is detonation velocity of TNT. PCJ is Chapman‑Jouguet pressure. A,B are empirical constant.
根据现场爆炸试验的观察,钢架无明显的位移和变形。因此,在模拟中将钢架简化为完全约束的刚体,以提高计算效率。

图6 无涂层靶板的模拟结果
Fig.6 Numerical result of uncoated concrete target plate

a. blast surface

b. back‑blast surface

c. side surface
图7 涂层混凝土靶板的试验和数值模拟破坏模式对比
Fig.7 Comparison of experimental and simulated failure modes of coated concrete target plate
基于建立的有限元模型,分析了涂层混凝土靶板的动态响应过程,并从能量的角度研究了涂层混凝土靶板的反直观行为发生机理。

a. blast surface

b. back‑blast surface

c. coating
图8 靶板损伤过程
Fig.8 Damage process of target plate
爆炸发生后,冲击波向四周传播,在0.07 ms时,冲击波到达靶板,此时涂层混凝土靶板背爆面混凝土开始产生损伤,迎爆面无明显变化(

图9 靶板迎爆面中心点的压力时程曲线
Fig.9 Pressure time history curve of the center point of the concrete target plate facing blast

图10 涂层中心点的位移曲线
Fig.10 Displacement curve of the coating center point
为进一步研究涂层混凝土靶板的反直观行为发生机理,

a. total energy

b. internal energy

c. kinetic energy
图11 材料的能量吸收曲线
Fig.11 Energy absorption curves of materials
对于涂层混凝土靶板,根据能量的变化可以将荷载作用后靶板的响应过程分为三个阶段,第一个阶段为涂层变形吸能阶段,第二个阶段为涂层弹性势能释放阶段。第三个阶段为混凝土损伤耗能阶段。第一个阶段为靶板的正向变形过程,第二个阶段和第三个阶段为靶板的反向回弹过程。
在第一个阶段(0.25~3.19 ms),混凝土总能量达到最大值后开始下降,涂层总能量随着涂层的变形开始增加(
在第二阶段(3.19~7.00 ms),混凝土的动能增加(
在第三个阶段(7.00~15 ms),混凝土的动能减小(
上述结果表明,涂层在爆炸荷载下能够通过高速变形吸收较多的能量,从而快速降低靶板的动能,减少混凝土的损伤程度。涂层在变形时,大部分能量用于粘性耗散和塑性耗散。在达到最大变形后,涂层中储存的弹性势能会传递给整个靶板,使得靶板发生回弹,涂层释放的弹性势能可用涂层最大内能和最终内能的差值表示,即
为了研究不同荷载对涂层混凝土靶板反直观行为的影响规律,开展了不同药量爆炸荷载作用于涂层混凝土靶板的数值模拟。TNT药量分别为25,50,75,100 g和150 g,不同药量的药柱长度一致,均为68 mm,药柱的半径分别为8.6,12.2,14.9,17.2 mm和21.1 mm。

图12 不同药量下涂层中心点的位移曲线
Fig.12 Displacement curves of coating center point for different charge masses
由

a. internal energy change

b. released elastic potential energy
图13 不同药量下涂层的内能曲线
Fig.13 Internal energy curves of coating with different charge masses

图14 粘结层失效
Fig.14 Coating adhesion failure
为了研究涂层厚度对涂层混凝土靶板反直观行为的影响规律,开展了不同涂层厚度的涂层混凝土靶板在75 g TNT爆炸下的数值模拟,涂层的厚度从2 mm到8 mm不等。

a. 2 mm

b. 4 mm

c. 6 mm

d. 8 mm
图15 不同涂层厚度下涂层的损伤
Fig.15 Damage of coatings with different coating thicknesses
由
coating thickness / mm | Smax / mm | Sf / mm | Wf / mm |
---|---|---|---|
2 | 39.26 | -10.50 | -10.50 |
4 | 29.81 | -16.05 | -16.05 |
6 | 29.87 | -29.87 | -21.63 |
8 | 28.56 | -26.02 | -15.39 |
Note: Smax is the maximum displacement of coating center point. Sf is the final displacement of coating center point. Wf is the final deflection of target plate.

a. internal energy change

b. released elastic potential energy
图16 不同涂层厚度下涂层的内能
Fig.16 Internal energy curves of coatings with different coating thicknesses
在75 g TNT爆炸荷载作用下,不同涂层厚度下涂层的挠度均与荷载作用方向相反(

a. 2 mm

b. 4 mm

c. 6 mm

d. 8 mm
图17 不同涂层厚度下模型的最终形态
Fig.17 The final morphology of the model with different coating thicknesses
本研究通过数值模型研究了背爆面喷涂柔性聚脲的混凝土靶板在爆炸荷载作用下的动态响应规律,从能量变化的角度重点分析了反直观行为的发生机理。在此基础上,研究了炸药药量和聚脲涂层厚度对背爆面柔性聚脲涂层防护混凝土靶板反直观行为的影响规律,具体结论如下:
(1)背爆面柔性聚脲涂层防护靶板的反直观行为是由柔性聚脲弹塑性变形引起的应变能释放、混凝土靶板脆性破坏引起的能量耗散及两种材料间能量的相互转化共同作用导致的;
(2)炸药药量是决定靶板反直观行为能否发生的关键因素,在中等药量(50~100 g)下,会发生反直观行为,当药量过低时(25 g)时,由于涂层释放的应变能较少而无法发生反直观行为,当药量过高(150 g)时,由于涂层脱粘影响两种材料的能量转化,也无法发生反直观行为;
(3)在本研究的工况下,不同涂层厚度下靶板均发生了反直观行为。靶板反直观行为的弯曲程度与涂层厚度并非正相关,当涂层厚度较低(低于6 mm)时,随着涂层厚度的增加,涂层释放的弹性势能增加,靶板的弯曲程度增加;当涂层厚度较高(高于6 mm)时,随着涂层厚度的增加,在回弹过程中,涂层对靶板弯曲变形的限制作用增强,导致靶板的弯曲程度降低。
参考文献
冯加和, 董奇, 张刘成, 等. 聚脲弹性体在爆炸防护中的研究进展 [J]. 含能材料, 2020, 28(4): 277-290. [百度学术]
FENG Jia‑he, DONG Qi, ZHANG Liu‑cheng, et al. Review on using polyurea elastomer for enhanced blast‑mitigation[J]. Chinese Journal of Energetic Materials(Hanneng Cailiao), 2020,28(4):277-290. [百度学术]
ZHANG R, HUANG W B, LYU P, et al. Polyurea for blast and impact protection: A review[J]. Polymers, 2022, 14(13):2670. [百度学术]
WANG J G, REN H Q, WU X Y, et al. Blast response of polymer‑retrofitted masonry unit walls[J]. Composites Part B: Engineering, 2017,128:174-181. [百度学术]
FALLON C, MCSHANE G J. Impact mitigating capabilities of a spray‑on elastomer coating applied to concrete[J]. International Journal of Impact Engineering, 2019,128:72-85. [百度学术]
LI Z, CHEN L, FANG Q, et al. Experimental and numerical study on CFRP strip strengthened clay brick masonry walls subjected to vented gas explosions[J]. International Journal of Impact Engineering, 2019,129:66-79. [百度学术]
WU J, LIU Z C, YU J, et al. Experimental and numerical investigation of normal reinforced concrete panel strengthened with polyurea under near‑field explosion[J]. Journal of Building Engineering, 2022,46:103763. [百度学术]
黄微波. 喷涂聚脲弹性体技术[M].北京:化学工业出版社, 2005. [百度学术]
HUANG Wei‑bo. Spray Polyurea Elastomer Technology[M]. Beijing: Chemical Industry Press, 2005. [百度学术]
孙鹏飞,吕平,黄微波,等. 喷涂抗爆型聚脲钢筋混凝土板抗爆性能研究 [J]. 材料导报, 2021,35(S2): 642-648. [百度学术]
SUN Peng‑fei, LYU Ping, HUANG Wei‑bo, et al. Research of explosion resistance of sprayed anti‑blast polyurea reinforced concrete slab[J]. Materials Reports, 2021,35(S2): 642-648. [百度学术]
WANG W, WEI G S, WANG X, et al. Structural damage assessment of RC slab strengthened with POZD coated steel plate under contact explosion[J]. Structures, 2023,48:31-39. [百度学术]
WU G, WANG X, JI C, et al. Experimental and numerical simulation study on polyurea‑coated fuel tank subjected to combined action of blast shock waves and fragments[J]. Thin‑Walled Structures, 2021,169:108436. [百度学术]
SUN Y X, WANG X, JI C, et al. Experimental investigation on anti‑penetration performance of polyurea‑coated ASTM1045 steel plate subjected to projectile impact[J]. Defence Technology, 2021,17(4):1496-1513. [百度学术]
ZHANG L, WANG Y T, WANG X, et al. Investigation on the influence mechanism of polyurea material property on the Blast Mitigation of polyurea‑steel composite plate[J]. Structures, 2022,44:1910-1927. [百度学术]
刘保华,徐文龙,王成,等.接触爆炸下聚脲涂层增强钢板的抗爆性能 [J]. 兵工学报, 2024, 45(5):1637-1647. [百度学术]
LIU Bao‑hua, XU Wen‑long, WANG Cheng, et al. Contact explosion resistance of steel plates reinforced with polyurea coating[J]. Acta Armamentarii, 2024, 45(5):1637-1647. [百度学术]
SYMONDS P S, YU T X. Counterintuitive behavior in a problem of elastic‑plastic beam dynamics[J]. Journal of Applied Mechanics, 1985,52 :517-522. [百度学术]
GALIEV S U. Distinctive features of counter‑intuitive behavior of plates and shells after removal of impulse load[J]. International Journal of Impact Engineering, 1997,19(2):175-187. [百度学术]
DONG Q, LI Q M, ZHENG J Y. Further study on counter‑intuitive response of single‑degree‑of‑freedom structures[J]. International Journal of Impact Engineering, 2011,38(5):305-308. [百度学术]
李海旺,树学峰,秦冬祺.球面薄壳在径向朝外脉冲作用下反直观动力响应研究[J].北京工业大学学报,2011,37(1):98-103. [百度学术]
LI Hai‑wang,SHU Xue‑feng,QIN Dong‑qi. A study of counter‑intuitive behaviors of pin‑ended thin spherical shell subjected to rectangle impulse loadings[J]. Journal of Beijing University of Technology, 2011,37(1):98-103. [百度学术]
康婷,许金余,白应生,等.爆炸冲击荷载作用下拱结构的弹塑性动力响应研究[J].兵工学报,2013,34(9):1097-1102. [百度学术]
KANG Ting, XU Jin‑yu, BAI Ying‑sheng, et al. Elastic‑plastic analysis on dynamic response of arch subjected to explosive impact[J]. Acta Armamentarii, 2013,34(9):1097-1102. [百度学术]
ACKLAND K, ANDERSON C, NGO T D. Deformation of polyurea‑coated steel plates under localised blast loading[J]. International Journal of Impact Engineering, 2013,51:13-22. [百度学术]
BORRVALL T, RIEDEL W. The RHT concrete model in LS‑DYNA[C]//Proceedings of The 8th European LS‑DYNA user conference. 2011: 23-24. [百度学术]