Control Parametric Analysis on Improving Park Restoring Force Model and Damage Evaluation of High-Strength Structure

In the dynamic time-history analysis of structural elastoplasticity, it is important to develop a universal mathematical model that can describe the force-displacement characteristics for restoring force. By defining three control parameters (stiffness degradation, slip closure γ, energy degradation β), the Park restoring force mathematical model can simulate various components. In this study, the Park restoring force has been improved by adding two control parameters (energy-based strength degradation β e and ductilitybased strength degradation β d ). Based on the testing data, the constitutive model is input and 55 numerical models are developed to analyze the effects of various parameters on structural behavior.Conclusion. (1) β has determinative effect on structural behavior; the effect of β e is basically consistent with that of β; α has significant effect on shear forces and bending moments; γ has significant effect on displacements and accelerations; β d has significant effect on shearing forces, acceleration, and total energy consumptions. (2) Based on the classification of four types of damage level, the recommended values for α, γ, β, β e , and β d are presented. (3) Based on the testing data of high-strength columns, the recommended values for the five control parameters of the improved Park restoring force model are presented.


Introduction
In structural seismic analysis, the mechanical responses of materials and structural members under the action of elastoplastic dynamic load are important [1].And in these elastoplastic analyses, the restoring force model is playing the key role.Restoring force models are mathematical models established to describe restoring force and deformation according to the simplified relation curve between restoring force and deformation, where the most essential issue is how to accurately define and program a mathematical model for the restoring force in order to describe the structural member levels in the program.Currently, the most commonly adopted restoring force mathematical models are including the Polygonal Hysteretic Model (PHM) and the Smooth Hysteretic Model (SHM) [2,3].For restoring force models targeting structural member levels, mathematical expressions are mainly created through experimental studies on structural members, in which case sufficient preliminary testing data are needed to establish restoring force models through numerical fitting.When developing the program Inelastic Damage Analysis of Reinforced Concrete (IDARC), Park proposed the Triparameter-Park restoring force model [4], a mathematical model for restoring force that takes into account stiffness degradation, strength degradation, pinching slip, and other structural effects.The model is also subjected to the synchronous control by three parameters, that is, the stiffness degradation-based parameter , the energy degradation-based parameter , and the pinching degradation-based parameter .
The program IDARC adopts PHM and SHM [5] too, the former of which is mainly used to describe the restoring force performance of concrete structures, while the latter [6] is applied to the simulation of the restoring force hysteretic performance of steel structures.IDARC is a program jointly developed by University at Buffalo, the State University of New York, and The Earthquake Engineering Research Center, which is applied for nonlinear dynamic response time-history analysis and damage analysis in various types of structures.parameter , energy degradation-based parameter , and pinching degradation-based parameter , can be employed to determine stiffness degradation, energy degradation, and pinching effects.In particular, when  ⇒ ∞,  = 0, and  ⇒ ∞, the Clough model can be obtained [7], while when  = 2.0,  = 0.1, and  ⇒ ∞, the Takeda model can be obtained [8].Through parameters adjustment, the model can be also employed to simulate the T-Beam restoring force model, the shear wall restoring force model, and the masonry structure restoring force model: with  = 2.0,  = 0.1, and  = 0.5, the T-Beam model can be stimulated [5]; with  ⇒ 0,  ⇒ 0, and  ⇒ ∞, the origin-oriented model can be simulated [9]; and with  ⇒ ∞,  ⇒ 0, and  ⇒ 0, the slip restoring force model can be simulated [5].Therefore, this is a universal model that can be applied to the simulation of various types of structural members.

Modified Park Restoring Force Model and Five Control Parameters
The Triparameter-Park restoring force model is a mathematical model controlled by three parameters, that is, the stiffness degradation parameter , the strength degradation parameter , and the pinching degradation parameter .The modified Park restoring force model is controlled by five control parameters, that is, the stiffness degradation parameter , the ductility-based strength degradation parameter   , the energy-based strength degradation parameter   , the pinching degradation parameter , and the strength degradation parameter .
(1) Stiffness Degradation Parameter .It is used to control the stiffness degradation process of the unloading segment of the restoring force model of structural members.All the unloading lines intersect in opposite directions at one point with vertical coordinate (-axis)   : cur is current moment; 0 cur is current curvature;   is stiffness degradation factor;  0 is initial elastic stiffness;  cur is current stiffness (wherein, when ( cur , 0 cur ) is on the right side of the elastic segment,   =  +  ; when ( cur , 0 cur ) is on the left side of the elastic segment,   =  −  );  is stiffness degradation parameter, with a value range of  ∈ (200.0,0.0).
(2) Energy Degradation Parameter .It is used to set the ratio between the damage increment caused by the increment of the maximum displacement response   /  and by the normalized hysteretic energy dissipation increment /    .Park et al. [4] defined an energy degradation parameter when programming the restoring force model of IDARC, and they also adopted this parameter  in the subsequent energy-deformation two-parameter damage model used for the calculation of structural damage [10]: is increment of the maximum displacement;   is ultimate displacement;  is hysteretic energy dissipation increment;   is yield strength; / is shear-span ratio;  0 is axial compression ratio;   is longitudinal reinforcement ratio;   is stirrup ratio.
The modified Park restoring force model, as introduced, adopts the mathematical forms of the ductility-based strength degradation parameter   and the energy-based strength degradation parameter   .The parameters in the literature [5,11], based on the definitions of energy degradation parameters, are used to describe the strength degradation process of structural members caused by different factors, such as energy dissipation and ductility.The parameter  is employed to define   and   , based on two types of restoring force hysteretic curve rules: is ductility-based strength degradation parameter;   is energy-based strength degradation parameter (wherein, the presence of the term   represents the strength degradation of deformation and ductility increment, and the presence of the term   represents the strength degradation of hysteretic energy dissipation, and Δ represents the hysteretic energy dissipation increment);  +/−  is positive or negative yield moment;  +/− 0 is initial positive or negative yield moment; 0 +/− max is positive or negative maximum curvature; 0 +/−  is positive or negative ultimate curvature;  is moment; Δ is moment increment;  is hysteretic energy dissipation;  ult is hysteretic energy dissipation of monotonic load under nondegraded ultimate curvature, with value ranges of   ∈ (0.0, 0.6) and   ∈ (0.0, 0.6).
(3) Gap Closure Slip Degradation Parameter .It is introduced to describe the typical effects of fracture opening and closure, crack closure and bond slip caused by slip and pinching, and so forth.Slip is a phenomenon in which simulation is adopted, as the target point of loading for the crack closure point and the toward point of the loading force is a partial value of the yield moment, with a value range of  ∈ (1.0, 0.0): is yield moment; 0  is ultimate curvature;   is moment caused by slip and pinching; 0  is curvature caused by slip and pinching; 0  is yield curvature caused by slip and pinching; 0  is ultimate curvature caused by slip and pinching.

Control Parameter Analysis of the Modified Park Restoring Force Model
The modified Park restoring force model, based on different valuation methods for , , and , is able to simulate different structural member characteristics, such as the Clough model, the Takeda model, the origin-oriented model, and the slip model.This model has relatively high universality when it comes to mathematical models describing the restoring force, so that it applies to universally describe the different restoring force models and to establish the relation between generalized force and generalized displacement in the elastoplastic dynamic time-history analyses of different structures.
The modified Park model provides four control parameters to describe the constitutive relation, which is incomparable to other typical restoring force models.In the comparison carried out, the Smooth model [3] shows an advantage in terms of reducing calculation loops in the synchronous processing of motion equations, but it only applies to the description of state space.
(1) Stiffness Degradation Parameter .It is used to describe the degree of stiffness degradation, and its value variations describe the degree of stiffness degradation of the restoring force in the hysteretic process in four ranges; in the value range of (0∼200),  can satisfy the simulation of the restoring force performance of conventional structural members.Some studies have adopted the stiffness degradation parameter  to distinguish between bent members and shear members or to reach the conclusion that it exerts no influence on the strength or stiffness of structural members.There has been no systematic study devoted to quantitatively exploring the degree of influence of stiffness degradation parameter .
(2) Fracture Slip Closure Parameter .It is used to describe the degree of fracture slip closure of the restoring force of structural members in the hysteretic process in four ranges; in the value range of (0∼1.0), can satisfy the simulation of the restoring force performance of conventional structural members.There has been no systematic study devoted to quantitatively exploring the degree of influence of fracture slip closure parameter .
(3) Energy Degradation Parameter .A fitted expression of  is given in studies on the Park deformation-energy dissipation two-parameter damage model.The expressions of parameters   and  are obtained through the fitting analysis of the test data of 402 rectangular-section concrete columns and 132 H-section steel columns [12].In a study conducted in 2013, Rajabi et al. [13] proposed that the value of  should be set at 0.27 based on probability, with a standard deviation of 0.6, and they claimed that putting forward the fitted formula in the value range of (0.1-0.6) in the program would make it more applicable to the simulation of the structural members tested.Through some tests, the study [14] has obtained the expression of ; Chinese scholars [15] have also conducted some studies on energy parameter-based valuations.In the technical report of IDARC 7.0 [5], the value range of  is (0∼0.6).With the constant increase of both the experimental reinforcement ratio of concrete columns and the strength of concrete, in the latest structural test of the studying team [16], the deficiencies of energy degradation parameter  have been constantly verified by structural tests.However, to obtain the fitted and universal expression of , the different tests and systematic studies and analyses are necessary.When the parameter  adopts the ductility-based strength degradation parameter   and the energy-based strength degradation parameter   , the restoring force model and damage evaluation are better described.Nevertheless, it is necessary to analyze the influence of the parameters on internal forces of structural members (shear and moment ratio), structural responses (displacement and acceleration values), and structural damage (structural energy dissipation and damage value) and the influence of the values of energy degradation   and   on internal forces of structural members, structural responses, and structural damage.

Simulation Analysis of Test and Numerical Model
4.1.Test Data of the Numerical Model.Based on the testing data obtained in the preliminary stage, a test was conducted on the seismic performance of HRB400 high-strength concrete columns, and a nonlinear analysis was carried out in quasistatic loading of structural members.The test is detailed in the reference papers of Lin and Wang [17].Tables 1 and  2 provide the main test data of structural members and the measured data of concrete and reinforcement materials.IDARC preprocessor was adopted to set model parameters, and El Centro seismic wave was selected from the PEER Strong Motion Database.In the preliminary work done by Lin and Wang [18], IDARC program employed test data to conduct numerical simulation verification, the results of which revealed the feasibility of applying IDARC program to the damage analysis and evaluation of high-performance concrete structures.

Test Data of the Numerical Model
(1) Constitutive Relation of Concrete.The concrete, to adopt an unrestricted constitutive model, exploits the Kent-Park constitutive model [19] modified by Scott et al. [20], consisting of an ascending segment and a descending segment: −  represents the strain corresponding to the stress peak;  represents the strength enhancement coefficient of concrete under the confinement of stirrups;  represents the strain softening angle of concrete;   represents the compressive strength of concrete;  ℎ represents the stirrup yield strength;   represents the ratio between the stirrup volume and the core concrete volume;  ℎ represents the stirrup spacing.The numerical model values in the program are  −  = 0.203%,   = 53 MPa,  ℎ = 390 MPa, and  ℎ = 150 mm; for other values refer to the test data in Tables 1  and 2, while the Default value of  is calculated based on the program.
(2) Constitutive Relation of Reinforcement.The reinforcement adopts the ideal elastoplastic constitutive model, without taking into account the stress growth caused by strain hardening: represents the reinforcement strain;   represents the reinforcement stress;   represents the elastic modulus of reinforcement;   represents the yield stress of reinforcement;    represents the equivalent elastic modulus of strengthening stage;  0 represents the strain of strengthening stage;   represents the ultimate strain of reinforcement.Numerical model values are   = 3% and hardening segment of reinforcement    = 3416 MPa; for other values refer to the test data in Tables 1 and 2.

Numerical Modeling.
The constitutive relation was implemented based on test results in the preliminary stage [16], and the IDARC was adopted to conduct an elastoplastic dynamic time-history analysis on the test pieces.Programming different restoring force model control parameters, shear, moment ratio, displacement response, acceleration response, structural total energy dissipation, and structural cumulative damage value of high-performance columns under different parameter levels were obtained.The influence of different parameter levels on the above six structural responses was also analyzed.After establishing a total of 55 models, the data of the first three groups, that is, (1)-( 3), were adopted to comparatively analyze the relations between the energy degradation parameter  on one hand and the stiffness degradation parameter  and the fracture slip closure parameter  on the other hand, while the data of the last three groups, that is, (3)-( 4), were adopted to comparatively analyze the relations between the energy degradation parameter  on one hand and the ductility-based strength degradation parameter   and the energy-based strength degradation parameter   on the other hand.
(1) Stiffness Degradation Parameter . was given a total of 11 independent parameter controlled variables, that is, 200, 135, 91, 62, 42, 28, 19, 13, 9, 6 and 4; the fracture slip pinching parameter  was set at its Default value, and the energy degradation parameter  was also set at its Default value (numbered as  1 ∼ 11 , a total of 11 models).For the convenience of comparative expression in the same range,  was provided with normalization processing, so that  ∈ (1.0, 0.0).
(2) Slip Pinching Parameter . was given a total of 11 independent parameter variables, that is, 1.0, 0.91, 0.81, 0.72, 0.62, 0.53, 0.43, 0.34, 0.24, 0.15, and 0.05; the stiffness degradation parameter  was set at its Default value, and the energy degradation parameter  was also set at its Default value (numbered as  1 ∼ 11 , a total of 11 models).A value range was given to ensure that  ∈ (1.0, 0.0).
(3) Energy Degradation Parameter . was given a total of 11 independent parameter variables, that is, 0.01, 0.07, 0.13, 0.19, 0.25, 0.31, 0.36, 0.42, 0.48, 0.54, and 0.60; the stiffness degradation parameter  was set at its Default value, and the slip pinching parameter  was also set at its Default value (numbered as  1 ∼ 11 , a total of 11 models).For the convenience of comparative expression in the same range,  was provided with normalization processing, so that  ∈ (0.0, 1.0).

Analysis of Influence of Control Parameters on Structural Responses and Damage
After performing elastoplastic dynamic time-history analyses on the 55 numerical models established, the results are analyzed and summarized in the following paragraphs.↓ means slightly decreased; ↓↓ means obviously decreased; ↓↓↓ means significantly decreased; ∼ means slightly fluctuated.

Influence Analysis of Stiffness Degradation Parameter 𝛼.
See Table 3. Figure 1: Influence on main structural responses of the numerical model by varying values of energy degradation parameter , stiffness degradation parameter , and fracture slip closure parameter .Note: for the convenience of comparative expression of various parameters in the same range, normalization processing was provided here, so that the stiffness degradation parameter  ∈ (1.0, 0.0) and the energy degradation parameter  ∈ (0.0, 1.0).

Influence Analysis of Closure Slip Degradation Parameter
. See Table 4.

Influence Analysis of Energy Degradation Parameter 𝛽.
See Table 5 and Figure 1.  6 and 7, respectively.

Influence Analysis of Energy Degradation Parameter 𝛽.
The analysis results were the same with those provided in Section 5.1.3.

Influence Analysis of Energy Degradation-Based Strength
Degradation Parameter   .See Table 6.

Influence Analysis of Energy Degradation-Based Strength
Degradation Parameter   .See Table 7 and Figure 2.

Analysis of the Influence of Restoring Force Model
Parameters on Structural Performance.Based on both the comparative analysis of the energy degradation parameter , the stiffness degradation parameter  and the slip closure parameter , and the comparative analysis of the energy degradation parameter , the ductility-based strength degradation parameter   , and the energy-based strength degradation parameter   , the degrees of influence of the various parameters on the main structural performance are summarized and the results of the evaluation and analysis are shown in Table 8.Meanwhile, on the basis of a comprehensive analysis of the influence of various parameters on structural responses in their variation process, value ranges are categorized into four grades (i.e., Default, Mild, Moderate, and Severe) for the five control parameters, as shown in Table 9.After all, by the parameters values of stiffness degradation parameter  = 9.0, slip closure parameter  = 0.15, energy degradation parameter  = 0.54, energy-based strength degradation parameter   = 0.54, and ductility-based strength degradation parameter   = 0.48, the obtained hysteretic restoring force curve and the damage evaluation of the structural members have the closest results to that of the high-performance structural member test.

Conclusions
Through studies on the most essential restoring force model within the elastoplastic dynamic time-history analysis, this paper comments on the five control parameters of the modified Park restoring model and adopts the preliminary  Note: "‰‰‰": relatively significant influence; "‰‰": moderate influence; "‰": relatively insignificant influence; "f": no or little influence.quasistatic test data to analyze the influence of these five parameters on the structural performances of high-strength concrete column members, such as shear value, moment value ratio, displacement value, acceleration value, total energy dissipation, and total damage value.
The results of the study indicates that the energy degradation parameter  exerts an essential control effect on the main performances of structural members; the energybased strength degradation parameter   exerts a degree of influence consistent with ; the stiffness degradation parameter  relatively significantly influences the shear value and moment value of structural members; the slip or fracture closure parameter  relatively significantly influences their displacement and acceleration response; and the Advances in Materials Science and Engineering ductility-based strength degradation parameter   relatively obviously influences their shear value, moment value, and total energy dissipation.
At the end, on the basis of a comprehensive analysis of the influence of various parameters on structural responses in their variation process, value ranges are categorized into four grades (i.e., Default, Mild, Moderate, and Severe) for the five control parameters as shown in Table 9. Considering also the quasistatic test of high-performance columns, it proposed the value ranges of the five control parameters, that is, stiffness degradation parameter 6 ≤  ≤ 9, slip closure parameter 0.15 ≤  ≤ 0.05, energy degradation parameter 0.54 ≤  ≤ 0.60, energy-based strength degradation parameter 0.54 ≤   ≤ 0.60, and ductility-based strength degradation parameter 0.48 ≤   ≤ 0.54.
The modified Park model is to be presented as follows.Based on Clough's and Takeda's models, it is a new universal PHM model established by borrowing the strength of the Triparameter-Park restoring force model, featured by vertexoriented and yield-oriented behaviors.The modified Park restoring force model also takes into account the strength degradation due to energy dissipation and ductility, and it leads to a five-parameter control of the restoring force mathematical model by introducing two additional parameters, that is, the energy-based strength degradation parameter   and the ductility-based strength degradation parameter   .A schematic diagram of the methodology of this paper is shown in Scheme 1.
Parameters , , and  on Structural Responses and Damage.As depicted in Figure 1, there are the influences on various main structural responses of the numerical model due to varying values of the energy degradation-based parameter , of the stiffness degradationbased parameter , and of the fracture slip closure parameter .The influences are also listed out for different value ranges of parameters , , and , in Tables 3, 4, and 5, respectively.↑↑↑means significantly increased; ↑↑ means obviously increased; ↑ means slightly increased; → means slightly unchanged.
Stiffness degrading parameter- Slip or crack-closing parameter- Strength degrading factor  Stiffness degrading parameter- Slip or crack-closing parameter- Strength degrading factor  Stiffness degrading parameter- Slip or crack-closing parameter- Strength degrading factor Parameters ,   , and   on Structural Responses and Damage.Depicted in Figure 2 there are the influences on various main structural responses of the numerical model due to varying values of the energy degradation-based parameter , of the ductilitybased strength degradation parameter   , and of the energybased strength degradation parameter   .The influences are also listed out for different value ranges of parameters   and   in Tables

2 )Figure 2 :
Figure 2: Influence on main structural responses of the numerical model by varying values of energy degradation parameter , ductilitybased strength degradation parameter   , and energy-based strength degradation parameter   .Note: for the convenience of comparative expression of various parameters in the same range, normalization processing was provided here, so that energy degradation parameters     ∈ (0.0, 1.0).

Table 1 :
Main information of test components in structural test.

Table 2 :
Measured material performance data of test piece.  : elastic modulus of reinforcement;   : yield strength of reinforcement;   : yield strain of reinforcement;   : ultimate tensile strength of reinforcement. cu : cubic compressive strength of concrete.  : axial compressive strength of concrete.  : compressive strength of concrete.: 10 −6 strain. Note:

Table 3 :
Analytical table of the influence of varying stiffness degradation parameter  on structural responses and damage.

Table 4 :
Analytical table of the influence of varying closure slip degradation parameter  on structural responses and damage.

Table 5 :
Analytical table of the influence of varying energy degradation-based parameter  on structural responses and damage.

Table 6 :
Analytical table of the influence of varying energy-based strength degradation parameter   on structural responses and damage.

Table 7 :
Analytical table of the influence of varying ductility-based strength degradation parameter   on structural responses and damage.

Table 8 :
Influence of the five control parameters of the modified Park restoring force model on structural responses.

Table 9 :
Value ranges proposed for the five control parameters of the modified Park restoring force model.