Development of Damage Type Viscoelastic Ontological Model for Soft and Hard Materials under High-Strain-Rate Conditions

: By improving the ZWT model, a principal structure model applicable to both soft and hard materials under dynamic loading conditions was obtained. Dynamic mechanical experiments were conducted using SHPB to obtain stress–strain curves for coal rock and foam concrete. The ZWT intrinsic model was simpliﬁed according to the dynamic impact characteristics of concrete, and the intrinsic model was established by introducing macroscopic damage quantity D and correction factor δ . The stress–strain curves of coal rock, foamed concrete, steel ﬁber concrete, granite, lightweight foamed concrete, and EPS concrete at different strain rates were used to validate the present constitutive model and prove the correctness of the model.


Introduction
Disasters caused by dynamic loads have occurred [1][2][3][4]. Soft and hard materials are commonly used materials in protective works [5][6][7][8][9]. Because dynamic load experiments are often difficult, numerical simulations are the most common means [10][11][12][13][14]. Many scholars have made outstanding contributions to the study of the present constitutive model, which directly affects the accuracy of numerical simulations [15,16]. The selection of the principal structure model directly affects the accuracy of the calculation, and most of the current principal structure models are only applicable to soft or hard materials, lacking a more general dynamic principal structure model. Therefore, it is important to establish a more general intrinsic structure model for the rapid selection of the intrinsic structure model for numerical simulation.
The ZWT instantonal model can describe static and impact loads well and the strain rates described by the ZWT instantonal model range from 10 −4 to 10 3 s −1 . Initially, Tang first proposed a nonlinear viscoelastic instantonal equation that can describe the kinetic behavior of epoxy resin, and the instantonal model can describe eight magnitudes of strain rates, which is called the ZWT instantonal model. Many scholars have improved the ZWT intrinsic model and established the instrinsic model for certain materials. The ZWT intrinsic model has been widely used in the fields of polymer materials, concrete materials, and soils. Luo et al. [17] studied the mechanical properties of TDE86 epoxy resin at different temperatures, established the ZWT intrinsic structure model of epoxy resin containing nine parameters, and verified the applicability of the ZWT intrinsic structure model at different temperatures by developing a subroutine. Liu et al. [18], according to the principle of building the constitutive model of polymer material, added the one-dimensional structure of the ZWT material constitutive model to the dashpot element in parallel, in which strain rate and coefficient of viscosity were introduced and the nonlinear viscoelastic constitutive model of Polycarbonate material was achieved; this intrinsic structure model contains 14 parameters. Xu et al. [19] developed a ZWT nonlinear viscoelastic model, which contains 10 parameters. The standard ZWT nonlinear viscoelastic model was selected to predict the elastic behavior of LNBR/epoxy composites under a wide range of strain rates in this study.
However, the damage of materials such as concrete and soil is guided by the generation of damage such as microcracks and microporosity, and these microscopic defects are the main cause of damage softening of materials. In order to reflect the main cause of microscopic defects on the softening of material damage and to reflect the effect of microscopic defects on the softening of material damage, numerous scholars have considered damage in improving the ZWT intrinsic structure model. Zhang et al. [20] established the dynamic principal structure relationship of BFRC by improving the ZWT model and fitted the stress-strain curve using the new equation, which has 10 parameters. Zhang et al. [21] developed an ontogenetic model of PFRC based on the ZWT ontogeny, which has nine parameters. Li et al. [22] combined the discrete ZWT model with the D-P yielding criterion to develop a permafrost intrinsic model with 10 parameters. Ma et al. [23] conducted dynamic mechanical experiments on frozen sandy soils by SHPB, and developed a dynamic native model of frozen sandy soils containing 12 parameters on the basis of the ZWT native model. Fu et al. [24] developed a viscoelastic intrinsic model with temperature sensitivity based on the ZWT intrinsic model, which has 15 parameters. Zhang et al. [25] tested the dynamic mechanical properties of permafrost at different temperatures and high strain rates using SHPB and established a dynamic intrinsic model of permafrost, which predicts results in good agreement with experimental results, and the model has 11 parameters. It can be seen that the ZWT intrinsic structure model is used in many materials, which proves that the ZWT intrinsic structure model has great applicability in research and industry all over the world, and therefore deserves further study. However, the existing model has many parameters, and the calibration parameters are tedious for researchers, which is not conducive to the promotion of the ZWT intrinsic structure model, so it is necessary to establish a dynamic intrinsic structure model with fewer parameters and which is more general.
In this paper, the ZWT principal structure model was simplified to contain only four parameters according to the kinetic properties, and the applicable range was extended by introducing damage and correction coefficients to establish a principal structure model containing only seven parameters. The fitting effect was verified by the SHPB test data of coal rock and foam concrete, and the obtained goodness of fit was high, which verifies the universal applicability of the intrinsic constitutive model established in this paper.

Establishment of Constitutive Model
The Zhu-Wang-Tang (ZWT) intrinsic model was proposed by Zhu Zhaoxiang, Wang Lili, and Tang Zhiping, the famous Chinese mechanics scholars, which can simulate the nonlinear viscoelastic intrinsic relationship of engineering plastics in the range of quasistatic load to impact load well. The main application objects of the original ZWT principal structure model are engineering plastics, including epoxy resin, plexiglass PMMA, polycarbonate PC, nylon, ABS, and PBT, etc. Concrete and rock is a typical non-homogeneous material, and there are a large number of defects inside the material, such as particle boundaries, pores, fissures and soft media, etc. The influence of damage on dynamic strength needs to be considered under high-strain-rate conditions. For the original ZWT instantonal model, it contains eight parameters, and there will be more parameters after considering the damage, and some equations of the original ZWT instantonal models do not work under the high-strain-rate conditions. Therefore, it is necessary to improve and simplify the ZWT instantonal model so that it can be applied to non-homogeneous materials. The development of a damage type viscoelastic instantonal model for soft and hard materials under high-strain-rate conditions is of practical significance for the study of the mechanical response of non-homogeneous materials under high-strain-rate conditions. Instantonal models are models that describe the short-term mechanical changes of materials and are mainly used in the field of dynamics.
Considering the strain hardening and strain rate hardening of the material in dynamic experiments, the existing Zhu-Wang-Tang constitutive model can be well expressed. The ZWT constitutive model is the equation expression used to express the impact resistance of materials, and the formula is as follows: α, β, E 0 represent the elastic constant of nonlinear spring; two integral terms represent the Maxwell body with different relaxation degrees, where θ 1 represents relaxation time of the Maxwell viscoelastic response at low strain rate, θ 2 represents the relaxation time of the Maxwell viscoelastic response at high strain rate, E 1 , E 2 represent the elastic constants of the Maxwell body. The constitutive model can be expressed by the mechanical model composed of a spring sticky pot, as shown in Figure 1.
icance for the study of the mechanical response of non-homogeneous materials under high-strain-rate conditions. Instantonal models are models that describe the short-term mechanical changes of materials and are mainly used in the field of dynamics.
Considering the strain hardening and strain rate hardening of the material in dynamic experiments, the existing Zhu-Wang-Tang constitutive model can be well expressed. The ZWT constitutive model is the equation expression used to express the impact resistance of materials, and the formula is as follows: α, β, 0 represent the elastic constant of nonlinear spring; two integral terms represent the Maxwell body with different relaxation degrees, where 1 represents relaxation time of the Maxwell viscoelastic response at low strain rate, 2 represents the relaxation time of the Maxwell viscoelastic response at high strain rate, 1 , 2 represent the elastic constants of the Maxwell body. The constitutive model can be expressed by the mechanical model composed of a spring sticky pot, as shown in Figure 1. In view of the actual mechanical properties under dynamic loading summarized in the experiment and the physical significance of the "Zhu-Wang-Tang" constitutive model, the following simplifications or improvements are made as follows: (a) Simplification I. For the first term of the equation of = 0 + 2 + 3 , the meaning is strain rate-independent equilibrium stress, which describes the nonlinear elasticity of the material. Often, the elastic part in dynamic experiments can be considered as linear. The trinomial of the first term can be simplified into monomial, that is, = 0 ; the nonlinear spring in the corresponding ZWT constitutive model becomes a linear spring. With this method, the constitutive model is simplified, and the simplified formula of the ZWT constitutive model is: (b) Simplification II. In the dynamic loading experiment of SHPB in this paper, the impact loading time is very short. The loading time scale of the low-frequency Maxwell body was between 10 0 and 10 2 s. Therefore, the loading time in this experiment is not one order of magnitude as that of the low-frequency Maxwell body, and the former is much smaller than the latter. Therefore, when the dynamic load is completed, the lowfrequency Maxwell body is not able to relax. Therefore, the second medium-lowfrequency Maxwell body in the ZWT constitutive model can be replaced by a simple In view of the actual mechanical properties under dynamic loading summarized in the experiment and the physical significance of the "Zhu-Wang-Tang" constitutive model, the following simplifications or improvements are made as follows: (a) Simplification I. For the first term of the equation of σ e = E 0 ε + αε 2 + βε 3 , the meaning is strain rate-independent equilibrium stress, which describes the nonlinear elasticity of the material. Often, the elastic part in dynamic experiments can be considered as linear. The trinomial of the first term can be simplified into monomial, that is, σ e = E 0 ε; the nonlinear spring in the corresponding ZWT constitutive model becomes a linear spring. With this method, the constitutive model is simplified, and the simplified formula of the ZWT constitutive model is: (b) Simplification II. In the dynamic loading experiment of SHPB in this paper, the impact loading time is very short. The loading time scale of the low-frequency Maxwell body was between 10 0 and 10 2 s. Therefore, the loading time in this experiment is not one order of magnitude as that of the low-frequency Maxwell body, and the former is much smaller than the latter. Therefore, when the dynamic load is completed, the low-frequency Maxwell body is not able to relax. Therefore, the second medium-low-frequency Maxwell body in the ZWT constitutive model can be replaced by a simple spring. The simplified physical model is shown in Figure 2, and the formula of the ZWT constitutive model can be further simplified as: 4 of 15 spring. The simplified physical model is shown in Figure 2, and the formula of the ZWT constitutive model can be further simplified as: The "Zhu-Wang-Tang" model simplified.
As shown in Figure 3, springs E0 and E1 are parallel, and a simple spring can replace the two parallel springs. The equivalent physical model is shown in Figure 3. The ZWT constitutive model is further simplified to: (c) Improvement I. The material cannot be completely evenly distributed in the production process of composite concrete. From a microscopic point of view, there are disordered micropores and micro-cracks inside, so the damage factor must be considered when establishing the constitutive model. The distribution of micro-pores and micro-cracks in concrete is very complex and has great discreteness. In order to simplify the model, the continuum damage mechanics method is adopted, that is, the composite concrete is assumed to be a continuous medium. The macroscopic damage D is introduced to evaluate the damage degree of composite concrete. According to the Lemaitre strain equivalence principle, the basic form of the constitutive model is as follows: In this formula, σa represents effective stress; σr represents the original stress; and D represents the damage variable. As shown in Figure 3, springs E 0 and E 1 are parallel, and a simple spring E a can replace the two parallel springs. The equivalent physical model is shown in Figure 3. The ZWT constitutive model is further simplified to: spring. The simplified physical model is shown in Figure 2, and the formula of the ZWT constitutive model can be further simplified as: The "Zhu-Wang-Tang" model simplified.
As shown in Figure 3, springs E0 and E1 are parallel, and a simple spring can replace the two parallel springs. The equivalent physical model is shown in Figure 3. The ZWT constitutive model is further simplified to: (c) Improvement I. The material cannot be completely evenly distributed in the production process of composite concrete. From a microscopic point of view, there are disordered micropores and micro-cracks inside, so the damage factor must be considered when establishing the constitutive model. The distribution of micro-pores and micro-cracks in concrete is very complex and has great discreteness. In order to simplify the model, the continuum damage mechanics method is adopted, that is, the composite concrete is assumed to be a continuous medium. The macroscopic damage D is introduced to evaluate the damage degree of composite concrete. According to the Lemaitre strain equivalence principle, the basic form of the constitutive model is as follows: In this formula, σa represents effective stress; σr represents the original stress; and D represents the damage variable. (c) Improvement I. The material cannot be completely evenly distributed in the production process of composite concrete. From a microscopic point of view, there are disordered micropores and micro-cracks inside, so the damage factor must be considered when establishing the constitutive model. The distribution of micro-pores and micro-cracks in concrete is very complex and has great discreteness. In order to simplify the model, the continuum damage mechanics method is adopted, that is, the composite concrete is assumed to be a continuous medium. The macroscopic damage D is introduced to evaluate the damage degree of composite concrete. According to the Lemaitre strain equivalence principle, the basic form of the constitutive model is as follows: In this formula, σ a represents effective stress; σ r represents the original stress; and D represents the damage variable.
Strain and strain rate will affect the damage evolution process, and the damage evolution equation is as follows: In this formula, • ε 0 = 10 −4 s −1 represents the reference value of the strain rate; • ε represents the strain rate; and A, α, and β represent undetermined coefficients respectively.
In summary, the Equations (4)-(6) can be obtained as: Considering that SHPB is approximately assumed to be constant strain rate, the constitutive model can be rewritten as follows: Although the improved ZWT constitutive model in this paper considers the damage factor, it still assumes that the material is continuous. In reality, there are a series of uncertainties in the material such as voids, cracks, and in homogeneous weave that may lead to early or delayed arrival of the stress maximum in the stress-strain curve. In order to characterize this uncertainty factor, correction factor δ was introduced to correct the fitting curve. Finally, the dynamic damage constitutive model with the correction factor is obtained as follows:

Experimental Research
In order to verify the general applicability of the viscoelastic constitutive model with damage established in this paper to hard and soft materials, dynamic mechanical experiments of hard and soft materials were carried out in this paper. The hard material was coal rock, and the soft material was selected as foam concrete. The stress-strain curve of the material was obtained by SHPB experiment. In this paper, the stress-strain curve obtained from the experiment was used to verify the constitutive model.

Preparation of Coal Rock
The coal sample for this test was taken from the 13-1 coal seam of the Zhangji Coal Mine of Huainan Mining Group, China. We took out the large coal body just exposed on the coal mining face, wrapped it in a sealed bag, and vacuumized it. After the raw coal arrived in the laboratory, it was cored with a rock drilling rig, and then cut and polished to prepare a coal-rock specimen with a diameter of 75 mm and a height of 35 mm. The drilling of the raw coal specimen is shown in Figure 4, and the dynamic mechanical specimen is shown in Figure 5. Zhang et al. [26] conducted a static compressive test on the same coal rock as that used in this paper, and the test result was 17.67 MPa.

Preparation of Foam Concrete
The composition of foam concrete includes cement, fly ash, water, and foam mixed uniformly. More detailed information on the raw materials is listed in Table 1. The amount of each raw material added is shown in Table 2. After sufficient mixing, concrete specimens with a diameter of 75 mm and a height of 35 mm were supported, as shown in Figure 6. In order to ensure the rationality of preparing foam concrete, this paper refers to the preparation parameters in the paper of He et al. [7], which recorded a static compressive strength of about 6 Mpa.

Preparation of Foam Concrete
The composition of foam concrete includes cement, fly ash, water, and foam mixed uniformly. More detailed information on the raw materials is listed in Table 1. Th amount of each raw material added is shown in Table 2. After sufficient mixing, concret specimens with a diameter of 75 mm and a height of 35 mm were supported, as show in Figure 6. In order to ensure the rationality of preparing foam concrete, this paper re fers to the preparation parameters in the paper of He et al. [7], which recorded a stati compressive strength of about 6 Mpa.

Preparation of Foam Concrete
The composition of foam concrete includes cement, fly ash, water, and uniformly. More detailed information on the raw materials is listed in T amount of each raw material added is shown in Table 2. After sufficient mix specimens with a diameter of 75 mm and a height of 35 mm were supporte in Figure 6. In order to ensure the rationality of preparing foam concrete, th fers to the preparation parameters in the paper of He et al. [7], which reco compressive strength of about 6 Mpa.

Experimental Method
Using the split Hopkinson bar experiment technology, the impact bar, incident bar, and transmission bar were high-strength bars, and the lengths were 400 mm, 4000 mm, and 2500 mm, respectively. The schematic diagram of the SHPB experimental device is shown in Figure 7. In this paper, the power system was driven by high-pressure gas, and data acquisition was conducted by USING a TP-5-12 semiconductor strain gauge and a DHHP-20 ultra-dynamic strain gauge, and DHHP dynamic signal analysis software was used for da-

Experimental Method
Using the split Hopkinson bar experiment technology, the impact bar, incident bar, and transmission bar were high-strength bars, and the lengths were 400 mm, 4000 mm, and 2500 mm, respectively. The schematic diagram of the SHPB experimental device is shown in Figure 7.

Experimental Method
Using the split Hopkinson bar experiment technology, the impact bar, incident bar, and transmission bar were high-strength bars, and the lengths were 400 mm, 4000 mm, and 2500 mm, respectively. The schematic diagram of the SHPB experimental device is shown in Figure 7. In this paper, the power system was driven by high-pressure gas, and data acquisition was conducted by USING a TP-5-12 semiconductor strain gauge and a DHHP-20 ultra-dynamic strain gauge, and DHHP dynamic signal analysis software was used for da- In this paper, the power system was driven by high-pressure gas, and data acquisition was conducted by USING a TP-5-12 semiconductor strain gauge and a DHHP-20 ultradynamic strain gauge, and DHHP dynamic signal analysis software was used for data analysis. The power system device is shown in Figure 8, the TP-5-12 semiconductor strain gauge is shown in Figure 9, the DHHP-20 ultra-dynamic strain gauge is shown in Figure 10, and the DYNAMIC signal analysis of DHHP is shown in Figure 11. Appl. Sci. 2022, 12, x FOR PEER REVIEW 8 of 16 ta analysis. The power system device is shown in Figure 8, the TP-5-12 semiconductor strain gauge is shown in Figure 9, the DHHP-20 ultra-dynamic strain gauge is shown in Figure 10, and the DYNAMIC signal analysis of DHHP is shown in Figure 11.     ta analysis. The power system device is shown in Figure 8, the TP-5-12 semiconductor strain gauge is shown in Figure 9, the DHHP-20 ultra-dynamic strain gauge is shown in Figure 10, and the DYNAMIC signal analysis of DHHP is shown in Figure 11.    Appl. Sci. 2022, 12, x FOR PEER REVIEW 8 of 16 ta analysis. The power system device is shown in Figure 8, the TP-5-12 semiconductor strain gauge is shown in Figure 9, the DHHP-20 ultra-dynamic strain gauge is shown in Figure 10, and the DYNAMIC signal analysis of DHHP is shown in Figure 11.

Verification of Stress Balance
In the SHPB test, the reliability of test results depends on the stress balance before and after the specimen failure, so it is necessary to verify the stress balance. Figure 12 shows the stress balance at both ends of the specimen. The sum of incident wave stress and reflected wave stress was basically equal to transmitted wave stress, so the experimental results are reliable.

Verification of Stress Balance
In the SHPB test, the reliability of test results depends on the stress balance before and after the specimen failure, so it is necessary to verify the stress balance. Figure 12 shows the stress balance at both ends of the specimen. The sum of incident wave stress and reflected wave stress was basically equal to transmitted wave stress, so the experimental results are reliable. Appl. Sci. 2022, 12, x FOR PEER REVIEW 9 of 16 Figure 11. Stress wave data analysis and processing software.

Verification of Stress Balance
In the SHPB test, the reliability of test results depends on the stress balance before and after the specimen failure, so it is necessary to verify the stress balance. Figure 12 shows the stress balance at both ends of the specimen. The sum of incident wave stress and reflected wave stress was basically equal to transmitted wave stress, so the experimental results are reliable.

Data Analysis
The SHPB experiment refers to the elastic stress wave propagation theory and adopts two basic assumptions [27][28][29]: (1) one-dimensional stress wave assumes that the cross section of the rod is always planar; (2) stress uniformity is assumed, that is, the stress wave propagates in the specimen, and stress is equal everywhere in practice. The experimental principle is to measure the incident wave, reflected wave, and pulse in the transmission rod by the strain gauge, and then deduce the stress-strain relationship of the specimen according to one-dimensional stress wave theory.
In this paper, the three-wave method was used to calculate the stress σ(t), strain rate . ε(t), and strain ε(t) of the specimen from the strain waveform in the elastic rod, namely: .
In the formula, E, C 0 , and A are the elastic modulus, elastic wave velocity, and sectional area of the rod, respectively. ε i , ε r , and ε t are incident wave strain, reflected wave strain, and transmitted wave strain, respectively. A 0 and L 0 are the cross-sectional area and length of the specimen, and t is the stress wave continuous specimen.

Experimental Results
Self-made foam concrete was selected as the soft material. The SHPB experiment was carried out by using a Φ50 mm aluminum rod. Hard materials were selected from the coal rock. The impact rod, incidence rod, and transmission rod were selected as the 75 mm steel rod, and the stress-strain curve under different strain rates was obtained. The coal-rock stress-strain curve is shown in Figure 13. The stress-strain curve of foam concrete is shown in Figure 14.
adopts two basic assumptions [27][28][29]: (1) one-dimensional stress wave assumes cross section of the rod is always planar; (2) stress uniformity is assumed, tha stress wave propagates in the specimen, and stress is equal everywhere in pract experimental principle is to measure the incident wave, reflected wave, and puls transmission rod by the strain gauge, and then deduce the stress-strain relation the specimen according to one-dimensional stress wave theory.
In this paper, the three-wave method was used to calculate the stress ̅( rate ̅ ( ), and strain ̅ ( ) of the specimen from the strain waveform in the elas namely: In the formula, E, C0, and A are the elastic modulus, elastic wave velocity, a tional area of the rod, respectively. , , and are incident wave strain, reflecte strain, and transmitted wave strain, respectively. A0 and L0 are the cross-section and length of the specimen, and t is the stress wave continuous specimen.

Experimental Results
Self-made foam concrete was selected as the soft material. The SHPB exp was carried out by using a Φ50 mm aluminum rod. Hard materials were selecte the coal rock. The impact rod, incidence rod, and transmission rod were selecte 75 mm steel rod, and the stress-strain curve under different strain rates was ob The coal-rock stress-strain curve is shown in Figure 13. The stress-strain curve concrete is shown in Figure 14.

Validation of the Viscoelastic Constitutive Model with Damage
In this paper, Formula (9) was used to fit the stress-strain curves of coal ro foam concrete under different strain rates, and the 1stOpt program was used to the calculation. The coal-rock fitting curve is shown in Figure 15, the foam concr ting curve is shown in Figure 16, and the specific parameters obtained are shown ble 3.

Validation of the Viscoelastic Constitutive Model with Damage
In this paper, Formula (9) was used to fit the stress-strain curves of coal rock and foam concrete under different strain rates, and the 1stOpt program was used to realize the calculation. The coal-rock fitting curve is shown in Figure 15, the foam concrete fitting curve is shown in Figure 16, and the specific parameters obtained are shown in Table 3.

Validation of the Viscoelastic Constitutive Model with Damage
In this paper, Formula (9) was used to fit the stress-strain curves of coal rock and foam concrete under different strain rates, and the 1stOpt program was used to realize the calculation. The coal-rock fitting curve is shown in Figure 15, the foam concrete fitting curve is shown in Figure 16, and the specific parameters obtained are shown in Table 3.  Ye et al. [30] obtained the stress-strain curves of steel fiber-reinforced concrete at different strain rates by SHPB. Cai et al. [31] studied the dynamic mechanical properties of granite. Feng et al. [32] studied the dynamic mechanical properties of lightweight foamed concrete. Li et al. [33] studied the impact properties of EPS concrete. In order to further verify the correctness of the constitutive model, this paper cites the experimental data of the above scholars. In this paper, Formula (9) was used to fit the stress-strain curves of steel fiber-reinforced concrete, granite, lightweight foam concrete, and EPS concrete at different strain rates, the 1stOpt program was used to realize the calculation, and the fitting curve is shown in Figure 17. The specific parameters obtained are shown in Table 2. It can be seen from Table 2 that the fitting degree R 2 was greater than 0.9, indicating that the constitutive equation established in this paper has wide applicability. 1stOpt is a mathematical optimization analysis synthesis tool developed by 7D-Soft High Technology Inc. in China, which is widely used in the fields of nonlinear regression, curve fitting, and nonlinear complex engineering mode parameter evaluation solutions, etc. The data in the references in this paper were extracted by GetData.  Ye et al. [30] obtained the stress-strain curves of steel fiber-reinforced concrete at different strain rates by SHPB. Cai et al. [31] studied the dynamic mechanical properties of granite. Feng et al. [32] studied the dynamic mechanical properties of lightweight foamed concrete. Li et al. [33] studied the impact properties of EPS concrete. In order to further verify the correctness of the constitutive model, this paper cites the experimental data of the above scholars. In this paper, Formula (9) was used to fit the stress-strain curves of steel fiber-reinforced concrete, granite, lightweight foam concrete, and EPS concrete at different strain rates, the 1stOpt program was used to realize the calculation, and the fitting curve is shown in Figure 17. The specific parameters obtained are shown in Table 2. It can be seen from Table 2 that the fitting degree R 2 was greater than 0.9, indicating that the constitutive equation established in this paper has wide applicability. 1stOpt is a mathematical optimization analysis synthesis tool developed by 7D-Soft High Technology Inc. in China, which is widely used in the fields of nonlinear regression, curve fitting, and nonlinear complex engineering mode parameter evaluation solutions, etc. The data in the references in this paper were extracted by GetData.   The parameter A is a constant in the principal structure model established in this paper, and determines the magnitudes of the parameters E a and E 2 . The parameter alpha and the parameter beta jointly determine the trend of the whole curve. The parameter alpha has a greater influence on the peak strain and peak stress of the curve. The parameter beta has a greater influence on the direction of the opening of the curve. The parameter E a is determined by the rising section of the curve, and E 2 is determined by the falling section of the curve. The intrinsic constitutive model established in this paper also has certain defects, for example, the model was only able to describe the large trend of stress-strain, and it could not accurately describe the curve wavy fluctuation. The description of material strength and weakness by the intrinsic constitutive model is described by E a and E 2 . Due to the limited experimental data, the model could only prove to be valid for a certain range of material strengths at present.

Conclusions
(a) In this paper, the ZWT constitutive model was simplified to only four parameters according to the dynamic characteristics. In this paper, the scope of application of the constitutive model was expanded by introducing damage factors and correction coefficients, and a constitutive model with only seven parameters was established.
(b) In this paper, the effect of damage on the dynamic load strength of materials under the action of all was considered, the damage theory was introduced into the ZWT intrinsic structure model, the stress-strain curves of a variety of materials were fitted using the improved damage-type viscoelastic intrinsic structure model equations, and the fitted curves and the measured curves had a good fit line. In this paper, the established intrinsic structure model can reflect the complex dynamic mechanical properties of materials such as the strain hardening effect and strain rate strengthening effect.
(c) In this paper, kinetic experiments were conducted on coal rock and foam concrete by SHPB to obtain stress-strain curves, and stress-strain curves of steel fiber concrete, granite, lightweight foam concrete, and EPS concrete by other authors were cited. In this paper, the viscoelastic damage-type intrinsic model equations were used to fit the curves for the abovementioned materials at high strain rates, and the fitted curves were in good agreement with the measured curves. Therefore, the intrinsic constitutive model established in this paper has good applicability, thus providing a certain reference value for the study of the intrinsic constitutive relationship of materials under high-strain-rate and engineering applications.