Abstract

To accurately and efficiently carry out numerical calculations for blasting demolition projects, based on the basic principles of engineering mechanics, a simplified concrete model is established. This model represents the postpeak behavior of concrete in compression and the peak behavior in traction. Based on the mechanical properties of steel and concrete materials, a new reinforced concrete model is proposed. The proposed reinforced concrete model has been implemented in the computer program LS-DYNA and used to simulate a real blasting demolition. The result shows that this model can accurately reflect the collapse and movement process of the chimney.

1. Introduction

In current engineering applications, concrete is considered to be in a safe state only when its peak strength is not exceeded. In the field of blasting demolition, the process of blasting demolition is accompanied by the failure of materials. The failure process of materials is related to their postpeak mechanical properties, and it is necessary to study the mechanical characteristics of concrete materials on the postpeak stage. Especially for demolition blasting, this postpeak characteristic determines whether the lower seat of demolition blasting occurs.

Furthermore, for blasting simulation modeling, a simple reinforced concrete material which can represent the overall behavior of concrete and rebar is often preferred due to its computational efficiency.

Wei and Yong [1] proposed a tensile constitutive model of reinforced concrete based on the existing continuous damage theory. Bifeng and Junjie [2] introduced the influence of steel bars based on the concrete microplane model and proposed a reinforced concrete dynamic constitutive model. Lin et al. [3] constructed a dynamic damage constitutive model of reinforced concrete based on its viscoelastic properties. Xinjin et al. [4] adopted the assumption of strain coordination and strength equivalence [5] and gave a constitutive model describing the damage of reinforced concrete materials, and they observed that after the reinforcement yielded, the material strength immediately entered the decline stage but did not consider the tension stiffening effect and compression stirrup restraint effect caused by the interaction of the concrete and reinforcement. Chao et al. [6] used the principle of stiffness equivalence to treat reinforced concrete as a continuous homogeneous material and simplified the constitutive relationship of reinforced concrete materials by using the piecewise linear method. He further modified the tensile stiffness effect caused by the interaction of reinforcement and concrete but did not consider the restraint effect of stirrups and the influence of reinforcement ratio on the tensile stiffness effect. At present, in the field of demolition blasting, the research on the mechanical properties of concrete mainly focuses on the prepeak stage of the material, and there is little research on the mechanical properties of the postpeak stage of the material. To perform more accurate numerical calculations for demolition blasting projects, these postpeak characteristics must be studied.

Commonly, the representation of a reinforced concrete element requires separate material models for concrete and steel and a proper mesh discretization. The need for a dense discretization grid could have a detrimental effect on calculation efficiency; moreover, the traditional finite element algorithm often forces the calculation to be terminated. In the finite element software LS-DYNA, this happens because existing materials are not able to reflect the large deformation of the postpeak response of concrete. In the numerical calculation of demolition blasting, the MAT_JOHNSON_HOLMQUIST_CONCRETE model and MAT_RHT model [7] are often used to calculate the mechanical properties of concrete materials. However, the HJC model uses a unified strain rate effect enhancement factor, which cannot accurately reflect the tensile and compressive strain rate effects of materials; although the RHT model can reflect material strain rate effects and less obvious postpeak softening effect, it cannot accurately reflect the softening characteristics of the material after the peak in the demolition blasting.

In this paper, the authors focus on the tensile and compressive behavior of concrete and steel with special attention to the postpeak behavior of concrete. Based on this knowledge, the envelope curve of a reinforced concrete material is proposed and used to set a predefined material in the computer program LS-DYNA.

2. Study on the Mechanical Properties of Concrete

2.1. Simplified Mechanical Properties of Concrete Materials Subjected to Axial Tensile Stress

The study of the mechanical properties of concrete materials in the prepeak stage needs to meet three conditions:

①The material is uniform. ②In the experimental, the material is in a state of uniform stress and uniform strain at all times. ③The size of the specimen does not change much in the experimental.

After applying axial tension to the specimen, corresponding stress and deformation are generated in the concrete material. When the deformation exceeds a certain degree, microfractures begin to occur in the specimen and the occurrence of these microfractures is random. The presence of these microfractures is represented by the damage variable D.

Assume that the original cross-sectional area of the specimen is A, and the effective cross-sectional area after the damage is . Then, the damage on the section of the material can be expressed as follows:

In the equation, is the damage variable, and its value range is [0, 1]. When D = 0, the material is not damaged, when D = 1, the material is destroyed; is the effective area of the material.

The nominal stress on the cross section of the concrete material based on the above three assumptions can be expressed as the following equation:

The true stress on the cross section of the material can be expressed as the following equation:

Then, the relationship between nominal stress and real stress can be obtained by equations (1)–(3) and it can be expressed as the following equation:

In the equation, is the nominal stress, is the real stress, and is the damage variable, and it can be calculated according to the expression in the Chinese Code for the design of concrete structures (GB 50010-2010) [8].

The damage variable has a clear physical meaning, but the true effective area of material under uniaxial action is difficult to accurately determine. The equivalent strain hypothesis holds that the strain produced by the nominal stress acting on the damaged material is equivalent to the strain produced by the effective stress acting on the nondamaged specimen. Therefore, either by replacing the nominal stress in the nondamaged specimen with the effective stress after damage or replacing the elastic modulus of the nondamaged specimen with the effective elastic modulus of the damaged material, the constitutive equation of the damaged material can be obtained and expressed as the following equation:

In the equation, E is the elastic modulus of the nondamaged specimen and is the elastic modulus of the damaged material.

The constitutive relationship when the size of the concrete material does not change under uniaxial action can be expressed as the following equation:

Equation (6) is derived based on the above three assumptions. It has the same form as the concrete constitutive expression in the Chinese Code for the design of concrete structures (GB 50010-2010).

Under the action of axial tensile stress, concrete materials can withstand little deformation and exhibit brittle characteristics. Before reaching the peak tensile strength of the concrete material, its stress-strain relationship can be expressed by equation (6). It is assumed that under the action of axial tensile tension, when the strain at peak stress is exceeded, the postpeak stress directly becomes 0.

The mechanical properties of the concrete material under axial tension can be expressed as the following equation:

Under axial tension, the simplified mechanical properties of concrete materials are shown in Figure 1.

2.2. Optimization of Mechanical Properties of Concrete Materials under Axial Compression

Compared with the axial tensile mechanical properties, concrete materials can withstand greater loads and deformations under axial compression. With the continuous increase of the axial pressure of the concrete material, the internal damage of the material continues to increase simultaneously. When the deformation exceeds the strain at peak stress, the microfractures begin to expand and the material enters the postpeak state, which is characterized by a significant deformation of the material until it fails macroscopically. In the finite element calculation, it often happens that if the deformation of the material is too large, the calculation is forced to terminate. In this case, under the conditions of uniform material and uniform stress and strain, it is necessary to further study the mechanical properties of concrete material for large postpeak deformations in compression.

In general, the initial state of strain-softening reflects the transformation of the material from a continuum to a structure and the change in the geometric size of the material. The study of the mechanical properties of concrete materials under axial compression needs to meet the following three assumptions: ①the material is uniform; ②in the experimental, the material is in a state of uniform stress and uniform strain at all times; ③the size of the specimen occurs a macroscopic change in the test.

Under the action of axial compression, the size of the material changes greatly; especially in the postpeak stage, the size of the concrete material changes significantly and is accompanied by the destruction of the material. Therefore, under axial compression, the longitudinal microstrain of the concrete material can be expressed as the following equation:

In the equation, is the current length of the concrete material. Integrating the , the true axial strain of the concrete material can be obtained as the following equation:where is the initial length of the concrete material.

Equation (9) reflects the superposition of concrete material deformation at any time, which can accurately represent the real deformation of concrete materials under axial compression and reflect the softening behavior characteristics of concrete materials.

To accurately study the softening behavior characteristics of concrete materials under compression, the influence of the change of the material section due to the Poisson effect cannot be ignored. According to the conservation of volume, it can be written as follows:where is the initial cross-sectional area of the concrete material, is the initial length of the concrete material, is the uniaxial compression cross-sectional area of the concrete material at any time, and is the length of the concrete material at any time.

Then, when the concrete material undergoes softening behavior under axial compression, its stress-strain relationship can be expressed as the following equation:

From equations (9) and (11), the true stress of concrete material under axial compression can be expressed as the following equation:

The nominal strain of concrete material under axial compression can be expressed as the following equation:

From equations (9) and (13), the true strain of concrete material under axial compression can be obtained as the following equation:

Combining equations (6)–(12) and (14), the envelope curve of concrete under axial compression can be expressed with the following equation:

The optimized envelope curve of concrete materials under compression is shown in Figure 2.

Equation (15) considers the superposition of concrete material strain at any time and considers the change of material shape due to the Poisson effect under compression. It can be applied to the study of the mechanical characteristics of concrete material in the postpeak deformation.

3. Constitutive Model of Concrete Confined by Stirrups

The stirrups have a strengthening effect on the volume of concrete enclosed by the “cage” constituted by the longitudinal reinforcements and the stirrups themselves. Based on the Mander confinement model, the strengthening effect of the stirrups on the concrete is accounted for in the proposed material model as an increase in strength and ductility of the material.

The Mander confined concrete constitutive model [9] was proposed by Mander in 1988, and in his paper, the enveloper curve of concrete proposed by Popovics is adopted, which is suitable for the representation of concrete strength under low strain rate and cyclic loading. The shape of the envelope curve is shown in Figure 3 for two different confinement conditions.

The full stress-strain curve of the concrete in the confinement zone can be expressed as the following equation:where is the peak stress of confined concrete and is the stress of confined concrete.

Define x as the ratio of confined concrete longitudinal compressive strain to peak compressive strain, and it can be expressed as the following equation:

In the equation, is the longitudinal compressive strain of confined concrete and is the peak compressive strain of confined concrete.

The correlation coefficient of concrete brittleness r is defined, and it is given by the following equation:where is the initial elastic modulus of concrete and is the secant modulus of constrained concrete.

The longitudinal peak compressive strain of confined concrete can be expressed as the following equation:

In the equation, is the peak strength of confined concrete, is the peak strength of unconfined concrete, and is the strain at peak stress of unconfined concrete.

The strain of confined concrete when the stirrup breaks can be expressed as the following equation:where is the confined concrete strain when the stirrup breaks, is the volume reinforcement ratio of the stirrup, is the yield strength of the stirrup, and is the ultimate strain of the constrained stirrup.

4. A New Simplified Model of Reinforced Concrete

4.1. Steel Modeling

In this study, the constitutive behavior of reinforcing steel is represented by the elastic-perfectly plastic law represented in Figure 4:

4.2. Calculation of Equivalent Elastic Modulus and Equivalent Stress

In the axial compression analysis, the stirrups have a strengthening effect on the confined concrete. To simplify the modeling workload, the strengthening effect of stirrups on the concrete can be represented by the increase in concrete strength. Assume that the section size of the reinforced concrete element is , and the area of the longitudinal reinforcement is .

The longitudinal reinforcement ratio can be calculated with the following equation:

In the equation, is the longitudinal reinforcement ratio, is the number of longitudinal bars, is the sectional area of a single longitudinal bar, b is the section width, and h is the section length.

The concrete area can be expressed as the following equation:

Under the action of the axial force, provided that the strain in the reinforcements and the concrete is the same, according to the force balance condition, it is

From equations (23) and (24), it can be found that when the strains of reinforcement and concrete are equal, the equivalent elastic modulus of reinforced concrete can be expressed as the following equation:where is the equivalent elastic modulus of reinforced concrete, is the elastic modulus of reinforcement, is the elastic modulus of concrete, and is the longitudinal reinforcement ratio.

The equivalent stress of the reinforced concrete material can be expressed with the following equation:where is the equivalent stress of reinforced concrete, is the stress of concrete, and is the stress of reinforcement.

4.3. Mechanical Properties of Reinforced Concrete Subjected to Axial Tensile Stress

The tensile mechanical properties of concrete materials are shown in Figure 1, and it is shown in equation (7). Before the strain at peak stress of concrete, the deformation of both reinforcement and concrete is small. It is considered that both reinforcement and concrete are in the elastic stage, and the strains of reinforcement and concrete are equal. Then, the tensile mechanical properties of reinforced concrete at the strain at peak stress of concrete can be expressed as the following equation:where is the peak tensile stress of reinforced concrete materials under axial tension and is the peak tensile strain of reinforced concrete materials (it equals the peak tensile strain of concrete). is the equivalent elastic modulus of reinforced concrete under tension.

When the concrete strain exceeds its strain at peak stress, it begins to enter the softening stage, which can be regarded as the concrete material failure under tensile tension. After that, the axial tensile tension is sustained entirely by the longitudinal reinforcement, and then the reinforcement enters the yielding state. At this time, the mechanical properties of reinforced concrete can be expressed as the following equations:where is the brittle rupture stress of reinforced concrete materials under axial tension and is the ultimate tensile strain of reinforced concrete under axial tension. is the yield stress of the reinforcement under axial tension; is the yield strain of the reinforcement under axial tension.

After the reinforcement enters the yield stage, necking begins to occur. To simplify the calculation, it is considered that tensile fracture occurs when the reinforcement reaches the yield strength. At this time, the reinforced concrete material can be truncated by tensile stress. The mechanical model of reinforced concrete material under axial tension can be expressed in Figure 5.

4.4. Simplified Mechanical Properties of Reinforced Concrete under Axial Compression

When the reinforced concrete material is subjected to axial pressure, the stirrups have a restraint strengthening effect on concrete. This effect can be expressed by the confinement model illustrated in the second section. The material is compressed and deformed under the action of axial pressure, from the beginning until it is broken. It is assumed that plane sections remain plane, the strain at each position on the section is equal, a perfect bond exists between concrete and reinforcement, and no external bulging phenomenon occurs. The strain of reinforced concrete material under axial compression can be expressed as the following equation:where is the compressive strain of reinforcement under axial compression and is the compressive strain of concrete under axial compression.

When both reinforcement and concrete are in the elastic deformation stage, the mechanical properties of reinforced concrete material under axial compression can be expressed as the following equation:where is the elastic limit stress of reinforced concrete material under axial compression, is the equivalent elastic modulus of reinforced concrete material under axial compression, and is the elastic compressive strain of concrete confined by stirrups.

When the confined concrete material begins to undergo plastic deformation under axial compression, although the stress continues to increase, its rate of increase gradually decreases, but the reinforcement is still in the elastic stage, and the proportion of pressure carried by the reinforcement gradually increases. When the confined concrete reaches its strain at peak stress, to simplify the calculation, it is considered that the reinforcement also enters the yield state; at this time, the peak stress of reinforced concrete material can be expressed as the following equation:where is the peak stress of reinforced concrete material under axial compression, is the peak stress of confined concrete, which can be calculated according to equation (16), is the yield stress of longitudinal reinforcement under axial compression, and is the longitudinal reinforcement ratio.

At this time, the strain at peak stress of reinforced concrete material under axial compression can be calculated by equation (19).

After the longitudinal reinforcement yields, its stress remains unchanged. Due to the restraint of the stirrups, the postpeak ductility characteristics of concrete are remarkable. When the stirrup fracture occurs, it is regarded as the reinforced concrete material compression failure. The mechanical properties of reinforced concrete materials under axial compression can be expressed as the following equations:where is the failure stress of reinforced concrete materials under axial compression, is the failure strain of reinforced concrete materials under axial compression, and and are the yield strength and reinforcement ratio of longitudinal reinforcement, respectively.

The failure strain of reinforced concrete materials under axial compression can be calculated by equation (20).

The stress-strain law of the proposed uniaxial reinforced concrete material under compression is shown in Figure 6.

5. Model Parameter Determination and Verification

The analytical application of this study is based on the computer software LS-DYNA, which is an explicit dynamic analysis program, based on the Lagrange algorithm, combined with ALE and Euler algorithm. Its explicit algorithm can simulate the large displacement deformation problem. The proposed envelope curve for reinforced concrete has been implemented in LS-DYNA. The implementation consisted in editing the material template MAT_PLASTICITY_COMPRESSION_TENSION [7]. This material template allows the user to assign the envelope curve for a custom material for both tensile and compressive stress. The other parameters required to define the material are density (RO), elastic modulus (E), Poisson’s ratio (PR), material failure parameter (FAIL), custom curve LCSRT and LCSRC, custom curve LCIDT, and LCIDC. LCSRT and LCSRC define a unique strain rate custom curve for tension and compression, and LCIDT and LCIDC define unique yield stress versus effective plastic strain for tension and compression.

The simulation of a tensile and compression test is performed in LS-DYNA on a simple reinforced concrete prismatic element, similar in shape to a C30 standard specimen, to obtain the stress-strain response of the proposed material. The parameter values are shown in Table 1, and the complete stress-strain curve is shown in Figure 7. It accurately reproduces the postcompressive mechanical properties of concrete reflected in equation (15) and the tensile stress truncation characteristics reflected in equation (7).

5.1. Axial Compression Experimental Verification of the Reinforced Concrete Simplified Model

In this section, the proposed material model is validated by comparing the force deformation response of actual reinforced concrete specimens with the ones obtained numerically. Details about the tests on the real specimens could be found in [10, 11]. Geometrical dimensions of reinforced concrete axially compressed short columns: length is 150 mm; width is 150 mm; height is 450 mm; concrete compressive strength fc = 58.6 MPa; longitudinal reinforcement yield strength is 530 MPa, and longitudinal reinforcement ratio is 0.84%; the protective layer thickness is 25 mm; stirrup strength grade is HPB235, with the stirrup ratio is 0.82%, and the hoop characteristic value is 0.03. The axial loading rate is 10 mm/s. In the finite element calculation of the axial compression short column, the concrete adopts the simplified model, and the Solid164 element is used to build the model. Both the upper and lower plates adopt rigid models, the lower plate is fully constrained, and the upper plate can only move in the longitudinal direction. By loading the upper plate, the centroids of the plate block and the column are kept on the same axis. Automatic single-sided contact between the plate and the column is adopted. The stress-strain point of the envelope curve of the simplified model is shown in Table 2, and the finite element model is shown in Figure 8.

Incorporating the above parameters into the numerical model, the numerical calculation result and experimental result of the load-strain curve of the simplified model can be obtained in Figure 9.

It can be seen from Figure 9 that there is a certain difference in the axial stiffness and maximum load capacity between the numerical calculation result and the experimental result. This is because the experimental material and the equivalent material have a certain difference in material uniformity. However, the overall trend of the numerical calculation result and the experimental result are consistent.

5.2. Compression and Bending Experimental Verification of the Reinforced Concrete Simply Supported Beam Simplified Model

The geometric dimensions of the reinforced concrete simply supported beam are as follows: length is 1400 mm, width is 150 mm, height is 200 mm, and support spacing is 1200 mm. The concrete strength grade is C25, and the protective layer thickness is 25 mm. The longitudinal reinforcement grade is HRB335, and the reinforcement ratio is 1.05%. The stirrup grade is HPB235, and the stirrup ratio is 0.628%. A two-point monotonic loading method is adopted, and the loading rate is 10 mm/s. In the finite element calculation of the simply supported beam compression and bending test, the simply supported beam adopts the simplified model, and the Solid164 element is used to build the model. Both the loading plate and the roller adopt the rigid body model. The loading plate can move in the vertical direction and rotate in the plane. Automatic single-sided contact is adopted in the test. The mechanical parameters of the simplified model are shown in Table 3, and the finite element model is shown in Figure 10.

Incorporating the above parameters into the numerical model, the numerical calculation result and experimental result of the load-displacement curve of the reinforced concrete simply supported beam simplified model are shown in Figure 11.

The numerically calculated load-displacement curve has certain volatility after the elastic phase, but after averaging the vibration data, it can be seen from Figure 11 that the numerical calculation result can reflect the changing trend of the experimental result.

6. Application

In this section, the proposed reinforced concrete model is applied to simulate a real blasting demolition project. The analytical simulation is then compared to the real demolition event using high-speed photographic observation. The experimental result of chimney blasting demolition comes from the literature [12].

The height of a reinforced concrete chimney is 210 m, its bottom outer radius is 9.97 m, wall thickness is 600 mm, top outer radius is 4.23 m, wall thickness is 200 mm, and the concrete strength grade is C30. The average reinforcement rate of longitudinal reinforcements is 0.50%, and the average reinforcement rate of stirrups is 0.45%. The upper incision is a regular trapezoidal and was opened at 110 m, the incision central angle is 210°, the incision height is 2.5 m, and the angle of the directional window is 45°. The lower incision is a regular trapezoid and was opened at the bottom of the chimney, the incision central angle is 216°, the incision height is 3.5 m, and the angle of the directional window is 45°. The upper incision detonates 3 s before the lower incision.

The reinforced concrete chimney finite element equivalent model is shown in Figure 12. The model of the chimney is constituted of three parts: the upper blasting incision, the lower blasting incision, and the equivalent reinforced concrete cylinder. The equivalent reinforced concrete cylinder adopts the proposed material, and the Solid164 element is used to build the model. The ground adopts the SOLID164 element and rigid body model, the rigid body density is 2000 kg/m3, the elastic modulus is 30 GPa, and Poisson’s ratio is 0.2. Gravity is achieved through LOAD_BODY_Y. The contact between the cylinder and ground is defined through CONTACT_ERODING_SINGLE_SURFACE. The blasting incision is realized by the time control parameter through ∗MAT_ADD_EROSION. The elements in the blasting incision are deleted by one-time control, without considering the impact of the explosive explosion on the formation of the incision. The parameters of the equivalent reinforced concrete cylinder are shown in Table 4.

Bringing the above parameters into the numerical model, the numerical calculation result and the experimental result of chimney blasting demolition can be obtained in Figure 13.

By controlling the detonation time of the upper incision and the lower incision of the model, making the initiation time of the model incision is the same as that of the actual incision. It can be seen from Figure 13 that the measured collapse process of the chimney is the same as the numerical calculation result before the 4th second. After the 4th second, the upper structure of the numerical model shifted slightly to the right compared with the actual collapse of the chimney, which was caused by the unevenness of the chimney material and weathering. However, on the whole, the measured collapse process of the chimney has a high degree of consistency with the numerical calculation result. The user-defined function of the simplified model allows users to perform equivalent processing of complex materials. On the premise of ensuring accuracy, the computational efficiency has increased by 30% and 35%, respectively, compared with the HJC model and the RHT model. At the same time, the correctness and rationality of the selection of the parameters of the simplified model are also verified.

7. Conclusion

(1)The simplified reinforced concrete model illustrated in this paper can not only reflect the softening behavior characteristics of concrete but also reflect the characteristics of concrete material tensile stress truncation. On the premise of ensuring the calculation accuracy, the modeling workload is greatly simplified.(2)The simplified model has the function of self-defining compression and tension parameters, which can accurately reflect the constitutive relationship of the reinforced concrete model, and greatly improves the efficiency of finite element calculation.(3)Demolition blasting is nonrepeatable and highly dangerous, and the effect of numerical calculation as an effective analysis method is becoming more and more obvious. The numerical calculation result of the simplified reinforced concrete model used in this paper is highly consistent with the actual high-speed photographic observation result. Therefore, the numerical calculation technique can provide a reference for similar blasting demolition projects.

Data Availability

The numerical data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

Acknowledgments

This work was supported by the National Science Foundation of China (grant nos. 50474040 and 50674078) and the National High Technology Research and Development Program of China (grant no. 2012AA061201).