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 ductility-based 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.
Natural Science Foundation of China50578066Project of National Public Welfare Industry Research201111020-21. 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. The three control parameters, stiffness degradation-based 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.
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 vertex-oriented 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 βe and the ductility-based strength degradation parameter βd. A schematic diagram of the methodology of this paper is shown in Scheme 1.
Programming flowchart.
2. 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 βd, the energy-based strength degradation parameter βe, 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 (y-axis) αMy:(1)RK+=Mcur+αMyK0∅cur+αMy,(2)Kcur=RKK0,(3)Kcur=Mcur+αMyK0∅cur+αMyK0.Mcur is current moment; ∅cur is current curvature; RK is stiffness degradation factor; K0 is initial elastic stiffness; Kcur is current stiffness (wherein, when Mcur,∅cur is on the right side of the elastic segment, My=My+; when Mcur,∅cur is on the left side of the elastic segment, My=My-); α 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 dδm/δu and by the normalized hysteretic energy dissipation increment dE/δuQy. 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]:(4)β=dδm/δudE/δuQy=dδmdE/Qy,β=-0.447+0.037ld+0.24n0+0.314ρt×0.7ρw.
dδm is increment of the maximum displacement; δu is ultimate displacement; dE is hysteretic energy dissipation increment; Qy is yield strength; l/d is shear-span ratio; n0 is axial compression ratio; ρt is longitudinal reinforcement ratio; ρw is stirrup ratio.
The modified Park restoring force model, as introduced, adopts the mathematical forms of the ductility-based strength degradation parameter βd and the energy-based strength degradation parameter βe. 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 βd and βe, based on two types of restoring force hysteretic curve rules:(5)My+/-=My0+/-1-∅max+/-∅u+/-1/βd1-βe1-βeHHult,ΔH=M+M+ΔM2Δ∅-ΔMRKK0.
βd is ductility-based strength degradation parameter; βe is energy-based strength degradation parameter (wherein, the presence of the term βd represents the strength degradation of deformation and ductility increment, and the presence of the term βe represents the strength degradation of hysteretic energy dissipation, and ΔH represents the hysteretic energy dissipation increment); My+/- is positive or negative yield moment; My0+/- is initial positive or negative yield moment; ∅max+/- is positive or negative maximum curvature; ∅u+/- is positive or negative ultimate curvature; M is moment; ΔM is moment increment; H is hysteretic energy dissipation; Hult is hysteretic energy dissipation of monotonic load under nondegraded ultimate curvature, with value ranges of βe∈(0.0,0.6) and βd∈(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):(6)Mγ=γMy,∅γ=γ∅γy+1-γ∅γu.
My is yield moment; ∅u is ultimate curvature; Mγ is moment caused by slip and pinching; ∅γ is curvature caused by slip and pinching; ∅γy is yield curvature caused by slip and pinching; ∅γu is ultimate curvature caused by slip and pinching.
3. 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 δu 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 βd and the energy-based strength degradation parameter βe, 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 βd and βe on internal forces of structural members, structural responses, and structural damage.
4. Simulation Analysis of Test and Numerical Model4.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.
Main information of test components in structural test.
Designed concrete strength grade
Axial compression ration=N/fckbh
Longitudinal reinforcement
Stirrup
Shear-span ratioλ=H/h0
Longitudinal reinforcement ratioρsv
Stirrup ratioρsv
C60
0.20
425
⍎10@150
1.9
3.70%
0.269%
Measured material performance data of test piece.
Elastic modulusEs (GPa)
Yield strengthfy (MPa)
Yield strainεy (με)
Ultimate tensile strengthfb (MPa)
fcu (MPa)
fc (MPa)
ft (MPa)
202
450
2233
585
67.1
53
3.47
Note: Es: elastic modulus of reinforcement; fy: yield strength of reinforcement; εy: yield strain of reinforcement; fb: ultimate tensile strength of reinforcement. fcu: cubic compressive strength of concrete. fc: axial compressive strength of concrete. ft: compressive strength of concrete. με: 10-6 strain.
4.2. 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:(7)σ-=Kfc2ε-εp--ε-εp-2,ε-≤εp-,σ-=Kfc1-Zε--εp-≥0.2Kfc,εp-≤ε-≤εu-,εp-=0.002K,K=1+ρsfyhfc,Z=0.53+0.29fc/145fc-1000+0.75ρsh′/sh-εp-.
εp- represents the strain corresponding to the stress peak; K represents the strength enhancement coefficient of concrete under the confinement of stirrups; Z represents the strain softening angle of concrete; fc represents the compressive strength of concrete; fyh represents the stirrup yield strength; ρs represents the ratio between the stirrup volume and the core concrete volume; sh represents the stirrup spacing. The numerical model values in the program are εp- = 0.203%, fc=53 MPa, fyh=390 MPa, and sh=150 mm; for other values refer to the test data in Tables 1 and 2, while the Default value of Z 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:(8)σs=Es·εs;0≤εs≤εsy,σs=fy;εsy<εs≤εs0,σs=fy+Esp·ε-ε0;Esp=0.025Es,εs0<εs≤εsu,σs=0;εsu<εs.
εs represents the reinforcement strain; σs represents the reinforcement stress; Es represents the elastic modulus of reinforcement; fy represents the yield stress of reinforcement; Esp represents the equivalent elastic modulus of strengthening stage; εs0 represents the strain of strengthening stage; εsu represents the ultimate strain of reinforcement. Numerical model values are εsy = 3% and hardening segment of reinforcement Esp=3416MPa; for other values refer to the test data in Tables 1 and 2.
4.3. 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 βd and the energy-based strength degradation parameter βe 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).
(4) Ductility-Based Strength Degradation Parameter βd. βd 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; βe=0.01, α=200, and γ=1.0 (numbered as βd1~βd11, a total of 11 models). βd was provided with normalization processing, so that βd∈(0.0,1.0).
(5) Energy-Based Strength Degradation Parameter βe. βe 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; βd=0.01, α=200, and γ=1.0 (numbered as βe1~βe11, a total of 11 models). βe was provided with normalization processing, so that βe∈(0.0,1.0).
5. 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.
5.1. Analysis of Influence of Control Parameters <inline-formula>
<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M219">
<mml:mrow>
<mml:mi>α</mml:mi></mml:mrow>
</mml:math></inline-formula>, <inline-formula>
<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M220">
<mml:mrow>
<mml:mi>β</mml:mi></mml:mrow>
</mml:math></inline-formula>, and <inline-formula>
<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M221">
<mml:mrow>
<mml:mi>γ</mml:mi></mml:mrow>
</mml:math></inline-formula> 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 degradation-based 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.
Analytical table of the influence of varying stiffness degradation parameter α on structural responses and damage.
Range of control parameter α
Shear force (kN)
Moment ratio
Displacement (mm)
Acceleration (mm/s^{2})
Total energy dissipation (kN⋅m)
Total damage value
α∈(200~62)
→
→
→
→
~
~
α∈(62~28)
↓
↑
↓
↑
↓
↓
α∈(28~19)
↓
↑↑
↓
↑
→
→
α∈(19~9)
↓
↑↑
↑↑
↓
→
→
α∈(9~4)
↑
↑↑
↑↑↑
↓↓↓
↓↓↓
↑↑↑
Analytical table of the influence of varying closure slip degradation parameter γ on structural responses and damage.
Range of control parameter γ
Shear force (kN)
Moment ratio
Displacement (mm)
Acceleration (mm/s^{2})
Total energy dissipation (kN⋅m)
Total damage value
γ∈(1.0~0.62)
→
→
→
→
~
↑
γ∈(0.62~0.34)
↓
→
↑
↑
↓
↑↑↑
γ∈(0.34~0.15)
↓↓
↑↑
~
↑↑
↓
↑↑
γ∈(0.15~0.05)
↓↓↓
↑↑↑
↓↓↓
↑↑↑
↑↑↑
↑
Analytical table of the influence of varying energy degradation-based parameter β on structural responses and damage.
Range of control parameter β
Shear force (kN)
Moment ratio
Displacement (mm)
Acceleration (mm/s^{2})
Total energy dissipation (kN⋅m)
Total damage value
β∈(0.01~0.13)
↓
↓
~
~
~
↑
β∈(0.13~0.36)
↓
↓
~
↓
↓
↑↑
β∈(0.36~0.48)
↓↓
↓↓
↑↑
↓↓
↑↑
↑↑
β∈(0.48~0.60)
↓↓↓
↓↓↓
↓↓↓
↓↓↓
↑
↑↑
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).
↑↑↑ means significantly increased; ↑↑ means obviously increased; ↑ means slightly increased; → means slightly unchanged.
↓ means slightly decreased; ↓↓ means obviously decreased; ↓↓↓ means significantly decreased; ~ means slightly fluctuated.
5.1.3. Influence Analysis of Energy Degradation Parameter <inline-formula>
<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M267">
<mml:mrow>
<mml:mi>β</mml:mi></mml:mrow>
</mml:math></inline-formula>
See Table 5 and Figure 1.
5.2. Analysis of Influence of Control Parameters <inline-formula>
<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M268">
<mml:mrow>
<mml:mi>β</mml:mi></mml:mrow>
</mml:math></inline-formula>, <inline-formula>
<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M269">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>β</mml:mi></mml:mrow>
<mml:mrow>
<mml:mi>e</mml:mi></mml:mrow>
</mml:msub></mml:mrow>
</mml:math></inline-formula>, and <inline-formula>
<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M270">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>β</mml:mi></mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi></mml:mrow>
</mml:msub></mml:mrow>
</mml:math></inline-formula> 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 ductility-based strength degradation parameter βd, and of the energy-based strength degradation parameter βe. The influences are also listed out for different value ranges of parameters βe and βd in Tables 6 and 7, respectively.
Analytical table of the influence of varying energy-based strength degradation parameter βe on structural responses and damage.
Range of control parameter βe
Shear force (kN)
Moment ratio
Displacement (mm)
Acceleration (mm/s^{2})
Total energy dissipation (kN⋅m)
Total damage value
βe∈(0.01~0.13)
↓
↓
~
~
↓
↑
βe∈(0.13~0.36)
↓
↓
↑
↓
~
↑
βe∈(0.36~0.48)
↓
↓↓
↑
↓↓
~
~
βe∈(0.48~0.60)
↓↓
↓↓
↑
↓↓↓
↑↑↑
~
Analytical table of the influence of varying ductility-based strength degradation parameter βd on structural responses and damage.
Range of control parameter βd
Shear force (kN)
Moment ratio
Displacement (mm)
Acceleration (mm/s^{2})
Total energy dissipation (kN⋅m)
Total damage value
βd∈(0.01~0.13)
~
~
~
~
~
~
βd∈(0.13~0.36)
↓
↓
↓
↓
↓↓
↓
βd∈(0.36~0.48)
↓
↓
↓↓
↓↓
↑↑↑
↓
βd∈(0.48~0.60)
↓↓
↓↓
↓
↓↓
↓↓
~
Influence on main structural responses of the numerical model by varying values of energy degradation parameter β, ductility-based strength degradation parameter βd, and energy-based strength degradation parameter βe. Note: for the convenience of comparative expression of various parameters in the same range, normalization processing was provided here, so that energy degradation parameters ββdβe∈(0.0,1.0).
5.2.1. Influence Analysis of Energy Degradation Parameter <inline-formula>
<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M300">
<mml:mrow>
<mml:mi>β</mml:mi></mml:mrow>
</mml:math></inline-formula>
The analysis results were the same with those provided in Section 5.1.3.
5.3. 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 βd, and the energy-based strength degradation parameter βe, 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.
Influence of the five control parameters of the modified Park restoring force model on structural responses.
Control parameter
Shear value (kN)
Moment ratio
Displacement (mm)
Acceleration (mm/s^{2})
Total energy dissipation (kN⋅m)
Total damage value
Stiffness degradation parameter α
★★
★★
★
★
★★
★
Slip or fracture closure parameter γ
★
★
★★
★★
★
★
Energy degradation parameter β
★★★
★★★
★★★
★★★
★★★
★★★
Energy-based strength degradation parameter βe
★★★
★★★
★★★
★★★
★★★
★★★
Ductility-based strength degradation parameter βd
★★
★
★
★★
★★★
☆
Note: “★★★”: relatively significant influence; “★★”: moderate influence; “★”: relatively insignificant influence; “☆”: no or little influence.
Value ranges proposed for the five control parameters of the modified Park restoring force model.
Control parameter
Default
Mild
Moderate
Severe
Stiffness degradation parameter α
α∈(200~62)
α∈(62~28)
α∈(28~9)
α∈(9~4)
Slip or fracture closure parameter γ
γ∈(1.0~0.62)
γ∈(0.62~0.34)
γ∈(0.34~0.15)
γ∈(0.15~0.05)
Energy degradation parameter β
β∈(0.01~0.13)
β∈(0.13~0.36)
β∈(0.36~0.48)
β∈(0.48~0.60)
Energy-based strength degradation parameter βe
βe∈(0.01~0.13)
βe∈(0.13~0.36)
βe∈(0.36~0.48)
βe∈(0.48~0.60)
Ductility-based strength degradation parameter βd
βd∈(0.01~0.13)
βd∈(0.13~0.36)
βd∈(0.36~0.48)
βd∈(0.48~0.60)
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 βe=0.54, and ductility-based strength degradation parameter βd=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.
6. 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 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 energy-based strength degradation parameter βe 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 ductility-based strength degradation parameter βd 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≤βe≤0.60, and ductility-based strength degradation parameter 0.48≤βd≤0.54.
Competing Interests
The authors of paper declare that there are no competing interests regarding the publican of this paper. The authors do not have a direct financial relation with the commercial identity that might lead to competing interests for any of the authors.
Acknowledgments
This work is supported by Natural Science Foundation of China (Grant no. 50578066) and Project of National Public Welfare Industry Research (201111020-2). The authors sincerely thank them for their support and funding.
KrätzigW. B.PetrynaY. S.Quasistatic seismic damage indicators for RC structures from dissipating energies in tangential subspacesSivaselvanM. V.Hysteretic models with stiffness and strength degradation in a mathematical programming formatRayT.ReinhornA. M.Enhanced smooth hysteretic model with degrading propertiesParkY. J.ReinhornA. M.KunnathS. K.ReinhornA. M.RohH.SivaselvanM. V.SivaselvanM. V.ReinhornA. M.Hysteretic models for deteriorating inelastic structuresCloughR. W.JohnstonS. B.TakedaT.SozenM. A.NielsenN. N.Reinforced concrete response to simulated earthquakesKabeyasawaT.ShioaraH.OtaniS.US-Japan cooperative research on RC full scale building test-Part5:discussion on dynamic response systemProceedings of the 8th World Conference on Earthquake Engineering1984San Francisco, Calif, USA627634ParkY.AngA.Mechanistic seismic damage model for reinforced concreteSivaselvanM. V.ReinhornA. M.Hysteretic models for deteriorating inelastic structuresParkY. J.AngA. H.WenY. K.Damage-limiting aseismic design of buildingsRajabiR.BarghiM.RajabiR.Investigation of Park-Ang damage index model for flexural behavior of reinforced concrete columnsWangC.XiaoJ.SunZ.Seismic damage assessment on cast-in-situ recycled aggregate concrete frame model structureChenL.JiangH.LuX.Modified Park-Ang damage model for reinforced concrete structuresLinH.-B.TangS.-G.LanC.Damage analysis and evaluation of high strength concrete frame based on deformation-energy damage modelLinH.WangQ.Seismic behavior on HRB400 high-strength reinforced concrete columnsLinH.WangQ.Damage evaluation of high-strength renforced concrete Structure based on IDRACKentD. C.ParkR.Flexural members with confined concreteScottB. D.ParkR.PriestleyM. J. N.Stress-strain behavior of concrete confined by overlapping hoops at low and high strain rates