Welcome to Acta Armamentarii ! Today is

Acta Armamentarii ›› 2023, Vol. 44 ›› Issue (12): 3897-3908.doi: 10.12382/bgxb.2023.0285

Special Issue: 爆炸冲击与先进防护

Previous Articles     Next Articles

The Loading Discrepancies in CONWEP and Fluid-structure Interaction Methods and the Dynamic Response Characteristics of Masonry Wall

SHANG Yulu1, XU Xuan1,*(), ZHANG Di2, YANG Jun1   

  1. 1 State Key Laboratory of Explosion Science and Technology, Beijing Institute of Technology, Beijing 100081, China
    2 Aerospace Science and Industry Intelligent Operation and Information Security Research Institute Co., Ltd., Beijing 100081, China
  • Received:2023-03-31 Online:2023-12-30
  • Contact: XU Xuan

Abstract:

The dynamic response behavior of masonry walls under large equivalent blast loading is investigated using a full-scale numerical model of windowed and windowless double-sided masonry walls, which is etablished by a cohesive zone method (CZM). A field experiment with a TNT equivalent of 500kg and blast distance of 13m is used as the basis for the study. The numerical simulations are carried out using two blast load application methods, conventional weapons effects program (CONWEP) and coupled-Eulerian-Lagrangian (CEL). The results indicate that the peak incident overpressures calculated by CONWEP and CEL methods are consistent with the experimental result compared to the swept overpressure data measured at the ground surface. However, the arrival time calculated by CEL method is 7.1% earlier than the test, and the rise time is four times that of the test. On the other hand, the overpressure decay calculated by CONWEP method is slower, and the positive pressure time and impulse volume are 50.7% and 42.56% of the test, respectively. It can be seen by comparing the load curve differences of the two blast loading methods that the rise time of overpressure-time curve calculated by CONWEP method is always constant, but the arrival time and rise time of shock wave calculated by CEL method are affected by the distance from explosive source and mesh size, and the larger the distance from explosive source and mesh size are, the longer the rise time is. The differences in the spatial and temporal distributions of shock wave flow field are reflected in the ideal hemispherical shape of wavefront calculated by CONWEP method and the oblate spherical shape of wavefront loaded by CEL method. For the reflected overpressure peak applied to the wall, the peak of CONWEP method loading is larger in each region and the decay rate is the same in all directions, while the decay rate loaded by CEL method in the y (vertical) direction is greater than that in the z (horizontal) direction. In the study of progressive damage law of wall, CEL loading can simulate the local damage characteristics of the wall more accurately. Moreover, the final damage pattern of the wall under CEL loading is closer to that of the test, while CONWEP method has a greater degree of damage.

Key words: CONWEP method, fluid-structure interaction, masonry wall, dynamic response characteristics

CLC Number: