Damage Analysis and Evaluation of High Strength Concrete Frame Based on Deformation-Energy Damage Model

A new method of characterizing the damage of high strength concrete structures is presented, which is based on the deformation energy double parameters damage model and incorporates both of the main forms of damage by earthquakes: first time damage beyond destruction and energy consumption. Firstly, test data of high strength reinforced concrete (RC) columns were evaluated. Then, the relationship between stiffness degradation, strength degradation, and ductility performance was obtained. And an expression for damage in terms ofmodel parameterswas determined, aswell as the critical input data for the restoring forcemodel to be used in analytical damage evaluation. Experimentally, the unloading stiffness was found to be related to the cycle number.Then, a correction for this changing was applied to better describe the unloading phenomenon and compensate for the shortcomings of structure elastic-plastic time history analysis. The above algorithm was embedded into an IDARC program. Finally, a case study of high strength RCmultistory frames was presented. Under various seismic wave inputs, the structural damages were predicted.The damage model and correction algorithm of stiffness unloading were proved to be suitable and applicable in engineering design and damage evaluation of a high strength concrete structure.


Introduction
Structural concrete is a material with quite outstanding cumulative damage characteristics; however, damage development and damage accumulation under different conditions are fundamental for understanding structural failure.Seismic damage greatly influences the bearing capacity of a building or structure during the follow-up service period and remaining life.Structural damage and its accumulation occur under dynamic loading, and a reasonable damage model must be constructed so that the damage can be described.
The structure damage caused by earthquakes is closely related to the maximum structural deformation and the low cycle fatigue effect caused by the accumulation of damage.In this paper, we have adopted a model based on two major forms of seismic damage [1]: (1) damage which occurs the first time the structure is taken beyond initial cracking and (2) damage due to subsequent energy dissipation.These forms of damage are incorporated into a double parameter damage model based on deformation and energy dissipation.
In this paper, we use the stiffness degradation trilinear model (DT3 model, Figure 1).The DT3 model uses three lines to describe the loading and restoring force skeleton curve and considers the stiffness degradation properties of the reinforced concrete structure or structural component.The model can describe in more detail the real restoring force curve of a reinforced concrete structure.Depending on whether a particular structure or structural member is expected to harden after yielding, the curve is divided into two parts-one part using the maximum top hardening degradation from the trilinear model, while the other part without this consideration.The degradation trilinear model, as used in this paper, is widely applied in engineering projects.
For these structures, the elastic-plastic time history analyses of the structure are carried out in order to determine the effect of accumulated damage on performance.Since the relationship between the unloading stiffness and ductility coefficient is also unknown, IDARC [2] analysis is introduced in this paper (Scheme 1).

Model for Deformation-Energy Consumption Damage.
We chose the Park and Ang's double parameters damage model for this paper [1].This model is based on the test results of 261 groups of reinforced concrete.The seismic damage is modeled via using a double parameter combination of the maximum deformation and cumulative hysteresis energy dissipation, and expressed as is ultimate deformation under monotonic loading;   is maximum deformation under earthquake;   is calculated yield strength (if   is smaller than   ,   =   );  is incremental absorbed hysteretic energy, and ,  are nonnegative parameters.Under elastic response, the value of  should theoretically be zero;  ≥ 1.0 signifies complete collapse or total damage.
Park and Ang represented the effect of cyclic loading on structural damage by the parameter .The absorbed hysteretic energy (excluding potential energy) was integrated up to the failure point for 261 groups of cyclic test data from columns and beams.Only those in which a sudden failure was clearly observed or gradual failure could be identified from the envelope curve were included.
When calculating the damage index, the load-deformation curve for each test is traced up to the failure point; at the point of failure ( = 1), the corresponding value of  is evaluated.In fact, these  values also depend on the value of local plastic rotation supply and the concrete strengths beyond.Based on the calculated  values [1], a negative correlation was observed between  and the confinement ratio,   , and weak positive correlations were observed between  and the shear span ratio, /, longitudinal steel ratio,   , and axial stress,  0 .Consider  = (−0.447+ 0.73   + 0.24 0 + 0.314  ) × 0.7  1 shows the correlation between  values from the model and the extent of damage.

Serviceable Range of Double Parameters Damage Model.
Park and Ang proposed their model based on a large number of test data, covering the following range of parameters: 1.0 < / < 6.6; 0.2 <   < 2.0; 0 <  0 < 0.52; 15.84 N/mm 2 <    < 41.34 N/mm 2 ; and 0.04 <   < 0.45.They used standard compressive tests with cylindrical samples (150 mm in diameter, 300 mm high) to determine the compressive strength and the failure strength.Beolchini et al. [3], taking test data for concrete columns under repeated loading cycles, believe that the value of  depends on   ,   ,  = / (shear span), and /    =  0 (axial compression ratio) as given in the following expression:  = (−0.28+ 0.06 + 0.47 0 + 0.19  ) × 0.66   .(5) More recently, high strength concrete (   > 41.34 N/ mm 2 ) has been widely used in engineering; such concrete is beyond the strength range of Park and Ang's model.So, the parameters for the damage model expression of high strength concrete must be adapted.These parameters can then be assessed to determine the damage mechanism, and an appropriate extension created for the damage calculation program.
Relevant tests of high strength concrete columns have been performed by authors, including component design, seismic tests, design of loading devices, and evaluation of the test data to determine the loading mechanism, failure mode,  − Δ curve, ductility coefficient, and yield displacement.
The tests data are presented in reference [4].Ductility is an important indicator of high strength, constrained, concrete.The comparison between component stiffness degradation, strength degradation, and ductility is shown in Figures 2 and  3.

Calculating Coefficient of Cyclic
Loading Effect .Defining  1 as that calculated using formula (3) and  2 as that calculated using formula (5), failure points can be obtained from  − Δ curves of each test extrapolating to the damage point,  = 1.Then, values of the damage parameter  3 can then be calculated; see    is quite significant.Component tests are therefore necessary for determining the appropriate parameters for high strength concrete.Table 2 shows the parameters obtained by fitting the expression (, ,   ,   ) for high strength concrete columns in the following procedure: (1) calculate formula (1) using  1 evaluated from formula (3); (2) calculate formula (2) using   2 evaluated from formula (5); (3) extrapolate each test curve to the failure point, and finally (4) determine the value at the damage point.The intermediate calculations are given in Table 2. Consider  = (−0.447+ 0.73 Formula ( 6) was used to calculate modified parameters for the deformation-energy damage model.These parameters were then used in the evaluation of testing member, and the new damage index was obtained, as shown in Table 2.The experimentally observed parameters-structural member ruptures, yield, and the final state-correlated well with the damage index.The damage index is close to 1 when the structure fails.The observed phenomena during the tests corresponded well with the predictions from the damage index.While, Park and Ang's model damage index at failure had a value more than 1.0.

Skeleton Curve of High Strength Concrete.
Through research on the restoring force model of high strength reinforced concrete columns, we have established a restoring force model which can adequately explain the nonlinear characteristics, using the principle of simplified calculation, in addition to considering the hysteresis characteristics of component, using the degrading trilinear model (Figure 1).
The route and hysteresis movement of the component can be described from the figure by tracing around the figure in the following sequence:  →  →   →  →   → ⋅ ⋅ ⋅ →  →  →   →   .The  section represents the elastic behavior, the  section is the strengthening,  is degradation,       is the reverse strengthening period, and     is the reverse degradation period.
The data at key points on the skeleton curve were determined experimentally and then used to initialize an analysis program that is implemented in IDARC.

High Strength Reinforced Concrete
Unloading Rigidity and Ductility Coefficient

Elastic Principle of Dynamic Unloading Stiffness.
The IDARC program for structural damage analysis had been used to simulate the earthquake response of a frame structure collocated with a high strength concrete column.To better simulate the dynamic performance of a high strength concrete structure, the dependency of the unloading stiffness and ductility coefficient is analyzed after being introduced into the IDARC program.Loading and unloading standards can be expressed as follows: loading (/  )  > 0; neutral loading (/  )  = 0; unloading (/  )  < 0 ( = 0).Unloading during structural vibration is complicated, as structure unloading stiffness matrix changes with time, leading to a time varying parameter system.At present, structure elastic-plastic time history analysis generally does not consider the influences that the unloading stiffness and motion state have on each other.These influences mean that the dynamic equations represented by tangent stiffness and secant stiffness are not equivalent, as the matrix is constantly changing during this process.
From the deformation properties and the restoring force model tests of high strength concrete columns, Figure 5 shows the relationship between the structural member stiffness and number of unloading cycles; it indicates that the variation of   () shows changes consistent with linear drag.The unloading stiffness   of for all samples changes with cycling, implying that the system is not fully elastic.The most important feature is that   () is approximately constant, although it varies from sample to sample.

Unloading Stiffness Correction Coefficient.
The unloading stiffness degradation used to determine the damage index conforms to observed evolution of structure damage.As the cyclic stress presenting in structural member, the strain  can be divided into reversible strain   and irreversible strain   and   includes the plastic strain   and the residual strain due to unrecoverable damage,   .From the principle of damage, In formula (7), destruction of the member is defined as  = 1.For the damage model, , we have used the deformation-energy double parameter combination; ] is only taken into account by the combination of factors of the maximum displacement damage coefficient (0.5 ≤ ] < 1.0).According to independent state variables principle in thermomechanics, the  − Δ state during damage evolution can be expressed as Formula (8) implies that using unloading stiffness   () to describe the structure damage value  is feasible.For ductile members (2.0 <  < 4.0), assuming that the damage function  =   ( 1 ,  2 , . . .,   ) is a monotonically increasing function, we can obtain the following from formula (8): where for the th cycle, unloading stiffness is    and elastic stiffness ratio is   .At cycle ( + 1), the unloading stiffness is   +1 and elastic stiffness ratio is  +1 . is defined as the ratio for the th and ( + 1)th cycles: Figure 6 charts the relationship between the measured value  and cycle number, .The value of  remained steady between 0.85 ≤  ≤ 0.95, implying that this parameter can be fixed during modeling; we suggest  = 0.9 for HSC.The test specimens reflect the phenomenon of degeneration of unloading stiffness with increasing displacement amplitude and can explain the fact that for an actual component the response excursion is always larger in one direction (i.e., offset), and eventually the accumulated errors as a result of this can lead to computational instability.Dynamic adjustment of unloading stiffness is simple but practically important.

Damage Evaluation of Highly Confined Concrete Structure Damage Evaluation Based on IDARC
4.1.Design Information of High Column Structure.In this paper, we model a structure comprising a six-floor frame of cast-in-situ high strength reinforced concrete.The strength grades of the concrete beams are C60-C80 and the main reinforcement and lateral reinforcement are HRB400 (  ), where the effects on structure ductility while utilizing HRB400 as reinforcement in high strength concrete have been studied to be limited under certain conditions [5,6].The columns use symmetric steel reinforcement, the cover thickness is 30 mm, and the gravitational load of each floor is:  1 = 6800 kN,  2 =  3 =  4 =  5 =  6 = 6600 kN.The earthquake intensity is at level VIII (Liedu scale, CSIS).According to the structure facade height, earthquake intensity, and design ground acceleration, the seismic design parameters are design grade 2, site category II, and design group 1 (feature period   = 0.35 s).Therefore, the average material strength of the high strength concrete is necessary to consider the dynamic response characteristics (including the restoring force) of the structural members; as the natural period of structure is  1 = 0.49 s, the design ground acceleration is 0.15 g.
The axial forces applied to each column in the structure are computed by elastic analysis.This paper adopts two indices to characterize the wave, the design characteristic period of site structure   , and structure basic natural period  1 .To choose the earthquake duration and natural period for the model, we have selected two strong earthquake records P0151 and P0994 from the PEER Strong Motion Database [7], the EL Centro wave (Figure 7) which is widely used in engineering, and also the RH4TG040 artificial seismic wave, ThiTG040 natural seismic wave, and TH2TG040 natural seismic wave from the EPDA module [8].

Damage Programming. IDARC is from the Earthquake
Engineering Research Centre of the State University of New York at Buffalo and is widely used for various types of nonlinear dynamic response time history analysis and damage analysis of structures.Details of our IDARC analysis will be given below.
(1) Hysteretic Rules Model.The program provides a threeparameter Park's model, bilinear hysteretic model, Kelvin's model, Maxwell's model, and smooth hysteretic model.For this paper, we have chosen a three-line attenuation model created by authors; all control information required by the preprocessor is obtained from the experimental skeleton curve.
(2) Nonlinear Dynamic Analysis.Nonlinear dynamic analysis is carried out using a combination of the Newmark-beta integration method, and the pseudoforce method.The equation is solved numerically, according to where [] is the lumped mass matrix of the structure; [] is the viscous matrix of the structure; [  ] is the tangent stiffness matrix; {Δ}, {Δ u }, and {Δ ü } are the incremental vectors for displacement, velocity, and acceleration in the structure; { ℎ } and { V } are the allocation vectors for the horizontal and vertical ground accelerations; Δ ẍ ℎ and Δ ẍ V are the increments in the horizontal and vertical ground accelerations;  corr is a correction coefficient; and {Δ err } is the vector containing the unbalanced forces in the structure.At the end of step  + Δ, the difference between the restoring force calculated using the hysteretic mode {} and the restoring force considering no change in stiffness during the step {  } yields the residual force {Δ err } = {} − {  }.This corrective force is then applied for the next time step of the analysis.The unbalanced forces are computed when moments, shears, and stiffness are being updated in the hysteretic model.We selected a step length of 0.02 seconds, to produce smaller changes at each step and thus avoid larger residual forces, which could lead to unstable calculations.according to the structural basic cycle for the overall structure damage index.This paper used the modified double parameters damage index and calculated the indices using the embedded computing methods (a), (b), and (c) for comparison.

Plastic Hinge
The damage to columns is light, the damage level of most columns is "medium destruction, " and the performance of the high-performance concrete is satisfactory.While the damage in the beams is much serious, the distribution of the damage index reflects the seismic design principle of "strong column, weak beam".Apart from the top floor, all beam ends act as plastic hinges, whilst only the column bottoms and a few centers of columns show hinge behavior.The beams act as hinges earlier than columns, nearly forming an overall symmetrical beam hinge system.The structure performance measured by the quantitative damage index is more reliable than displacement-based design.The shear-displacement hysteresis curves for the bottom of the structure are shown in Figure 11.Qualitative analysis of the hysteric performance of each floor can be performed based on the lateral displacement-shear hysteresis relationships under the action of the shock wave.Preliminary analysis of the graphs shows that the bottom layer has yielded, and the extent of rigidity degradation shows that the first shock transcends beyond destruction and the maximum displacement appears only for a few seconds.Most of the time, the displacement is relatively small, and the plastic energy dissipation and cumulative plastic deformation are small.

Seismic Damage Analysis of Whole Structure.
The modified parameter damage model was used to calculate damage parameters for all columns in each floor of the frame structure under the earthquake waves.The results are shown in Table 3.Both the results and value range of the modified damage model calculation are close to those observed in real damaged structures.Getting instantaneous snapshots from the program during time history analysis, we can calculate the instantaneous damage index for the whole structure as a function of the time period.We find that for this example wave used the damage index is where ( 0 ) initial is the initial fundamental period and ( 0 ) equivalent is the equivalent fundamental period.The overall damage can reflect more accurately the reaction of the whole structure.This is probably because that, using the vibration frequency attenuation method (i.e., the change in stiffness with time), it accurately simulates the real behavior of structures.The results show that the high strength concrete frame had been damaged, but not up to the extent of complete collapse; it reaches the seismic design goal of not collapsing under a major earthquake.The presented algorithm results in a damage index always less than 1.0.This modified double parameters damage model is correctly describing the performance of damaged structures, because it predicts the results that are quite close to those observed in real damaged structures.

Conclusions
In this paper, based on the deformation-energy damage model and using experimental test data, we determined the expressions for the coefficient of the cyclic load effect and values required for the double parameter damage model of members of a high strength reinforced concrete column member.Ductility is the key index of high-performance concrete, as found by the analysis of the relationships between stiffness degradation, strength degradation, and ductility.A model of the restoring force was used; this model allowed extrapolation of the data to the point of failure to obtain the restoring force.A trilinear model was employed to describe damage.Generally, the elastic-plastic time history analysis does not consider the relationship between unloading stiffness and motion state which influence each other.Thus, we evaluate the structure elastic-plastic time history analysis by introducing the unloading correction coefficient, which is a parameter that characterizes the change with time of the unloading stiffness matrix structure, and solve the equivalent problem between the structural dynamic equation indicated by tangential stiffness and secant stiffness.This gives a good description of the elastic-plastic structure vibration dynamic process, which can be used to describe the concrete structure unloading in the dynamic elastic-plastic state.We then implemented these modifications into the elastic-plastic time history analysis in IDARC program.
Presenting the case study of a high strength RC frame, the member dynamic response characteristics of the restoring force was calculated from the material average strengths, where the longitudinal bar and stirrup bar were high strength bar and the concrete strength of the beam-column was high.The structural damages were observed under various seismic waves, including two strong artificial seismic wave records from the PEER library, the El Centro wave, two natural seismic waves, and one artificial seismic wave from EPDA.Dynamic analysis was used to analyze the seismic performance and evaluate the degree of damage.Under a strong earthquake, the maximum displacement and angular displacement between floors are consistent.The maximum angular displacement is more than limits by code, but there is no obvious weak floor; the damage to the base of the structure accumulates up until the maximum impact of reaction.The distribution of the damage index reflects the seismic design criterion of "strong column, weak beam, " nearly forming a symmetrical beam hinge system.The lateralshear hysteresis curve is showing the yielding and the stiffness degradation degree.The displacements only appear to be exceeding initially and to be later small generally.

Figure 3 :
Figure 3: Relationship between strength degradation and ductility coefficient.

Figure 4 :
Figure 4: Comparison of values of  (parameter of damage model) obtained with different models.

( 3 )
Damage Index.IDARC incorporates three models for damage index, (a) classical Park-Ang's model; (b) a model based on low cycle fatigue damage characteristics; and (c)

Figure 8 :O
Figure 8: Maximum floor displacement and story drift angle of building in seismicity basic intensity 8 area.

Figure 9 :
Figure 9: Plastic hinge distribution of building.

Figure 11 :
Figure 11: Base shear-drift hysteretic curve of building first floor.

Table 2 .
Then, the comparison of parameters  1 ,  2 , and  3 for a variety of values of beam

Table 1 :
Damage evaluation rules from Park and Ang model based on value of .

Table 2 :
Calculation of  (parameter of damage model for cyclic loading).