Earthquake Damage Assessment for RC Structures Based on Fuzzy Sets

A global damage index based on multiple linear force-deformation curves in pushover analysis is presented to evaluate the integrated damage of reinforced concrete structure. The modified coefficient is provided considering the cyclic load and hysteresis energy.The number of inelastic cycles and the coefficient of hysteresis energy concentration are also introduced as damage indices. Hence, multiple damage indices about displacement and energy for performance-based design are considered. The relation of multiple damage indices or factors and the fuzzy damage set is presented by comprehensive fuzzy evaluation; hence, a performancebased multiple fuzzy seismic damage-assessment method for reinforced concrete frame structures is established. The method can be accomplished based on pushover analysis, code spectrum, and capacity spectrummethod.The fuzzy seismic damage-assessment method is verified through nonlinear analysis four different structures and the corresponding results and assessment conclusions are accurate.


Introduction
The reinforced concrete structure will inevitably be damaged and enter elastoplasticity deformation stage when it is subjected to strong ground motion.In order to determine the appropriate and economical selective intervention or rehabilitation strategy for an existing structure, it is necessary to assess accurately the actual resistance and the potential failure modes.It is also essential to carry out the damage assessment for a new designed building in order to verify that the performance and capacity meet the target specification of the seismic design guideline.Based on broad categories for advanced seismic design approaches, performance-based seismic design (PBSD) is at the forefront of earthquake engineering research [1][2][3].One of the primary aspects of PBSD is the realistic characterization of seismic structural damage and its direct incorporation in the design or performance evaluation methodology.The damage assessment based on PBSD is to ensure that the structures achieve the expected targets for different performance requirements.
Structural ductility, storey drift, element and connection rotation, dissipated energy, and fatigue of the structure are the parameters considered for traditional damage assessment [4].Based on the fundamental parameters, various local and global damage indices have been proposed to evaluate reinforced concrete structures.Research on the global damage indices is essential because the integral damage assessment should be accurate and involve the damage properties of local elements.However, the global damage indices presented until now have the disadvantages such as limited convergence and not considering cumulative damage.In addition, many damage indices are obtained from pushover analysis and capacity spectrum method whose purpose is to evaluate structural performance under earthquake by overprinting seismic response spectrum curve and structural capacity curve [5].A significant amount of research has shown that capacity spectrum method can simplify computation, consider the nonlinear performance of the structural elements and finally obtain effective results.Because capacity spectrum method has major simplifications, such as a multi-degree of freedom system being transformed into a single-degree of freedom system, nonlinear time-history analysis being replaced by pushover method, and accumulated dissipative effect not sufficiently considered, thus the evaluation is inexact in some cases, especially for irregular or complex structures.Hence, it is necessary to present a comprehensive Mathematical Problems in Engineering method to combine multiple damage parameters and obtain more accurate results.
Moreover, the uncertainty always affects the modeling and evaluation of real physical phenomena with different intensity, the damage of structures in earthquake is a vague concept because there exists not explicit classification criterion to evaluate different damaged members and the structure with different damaged information.The values from different damage-assessment methods are also different and even distinct.Hence, the seismic damage assessment and the decision problem are very complex for large and complex structure, which needs the integration based on multiclass information, and the combination of different information is a valid method to reduce the deviation and the uncertainty.In the past decades several alternative mathematical approaches to deal with uncertainty have been developed and the main categories methods include fuzzy set theory, artificial neural networks, interval analysis, and convex modeling.
In classical set theory, the membership of elements in a set is assessed in binary terms according to a bivalent condition; that is, an element either belongs or does not belong to the set.By contrast, fuzzy set theory permits the gradual assessment of the membership of elements in a set [6].This is described with the aid of a membership function valued in the real unit interval [0, 1].Fuzzy sets generalize classical sets, since the indicator functions of classical sets are special cases of the membership functions of fuzzy sets, if the latter only takes value 0 or 1.The fuzzy logic and the fuzzy set provide a powerful tool to solve the uncertainty problem and describe the fuzzy characteristics of the target.Especially, fuzzy logic is fit for describing and processing the uncertain information from multiple sources.In recent years, the methods such as the fuzzy logic system, artificial neural networks, and other soft computing methods develop rapidly, which provides a solid theoretical basis for information fusion and system assessment.
In this paper, the multiple damage indices on displacement and energy for reinforced concrete structures are discussed and applied based on performance-based seismic design and the comprehensive evaluation method is proposed by applying fuzzy set theory.The relationship between the multiple damage factors and the fuzzy damage set is established and an effective and accurate seismic damageassessment method is presented.

Damage-Assessment Index Based on
Displacement and Energy The state of damage of a component, a story, or the whole structure may be represented by an index.The damage index is used as an indicator to describe the state of the lateral load-carrying capacity and the reserve capacity of existing structures.Thus, the study on damage index and its availability is necessary.Some damage indices are calculated for each component of the building (local damage index).The component damage indices may be integrated using a weighting procedure to provide the global damage index for the structure.These damage indices have been formulated using response parameters of the structure that are obtained through analytical evaluation of structural response.There are several techniques and approaches for damage analysis of structures, such as pushover analysis, nonlinear time history analysis, and vulnerability analysis.The typical response-based damage indices include ductility ratio, interstory drift, slope ratio, maximum drift, flexural damage ratio, low cycle fatigue, final softening index and Park-Ang index.The damage indices such as interstory drift and maximum drift, are fundamental and essential for representing the displacement or deformation.
Though various local and global damage indices were presented, the most acceptable damage index is the Park-Ang damage index combining both ductility and cumulative hysteretic energy demand [7,8]: where   is the maximum experienced deformation,   is the ultimate deformation of the element,   is the yield strength of the element, ∫  ℎ is the hysteretic energy absorbed by the element during the response history, and  is the model constant, parameter, which was suggested to be 0.1 for nominal strength deterioration.The Park-Ang damage index also has deficiencies even though many modified method were proposed [9][10][11][12].A perfect damage index typically normalizes the damage on a scale of 0 to 1, where zero represents undamaged state while unity represents the collapse state of the building.The maximum of the Park-Ang damage index is over 1 and nearly close to 2 in some cases.In addition, the Park-Ang damage index is not directly suitable for nonlinear static pushover analysis or normal capacity spectrum method because the cumulative damage does not occur in this case.
Based on pushover analysis, a stiffness damage index of the global frame was presented by Ghobarah et al. [4]: where   is the initial slope of the base shear-top deflection relationship resulting from the pushover analysis of the frame before subjecting it to the earthquake ground motion and   is the initial slope of the same relationship but after subjecting the frame to the earthquake time history.The values of these damage indices range from zero to one depending on the amount of damage experienced.
Although the stiffness damage index has some advantages such as being concise and using pushover analysis twice without any other dynamic analysis, the limitation is also obvious.The damage evaluation is not accurate for moderate damage and collapse stage and the value may exceed 1 and the index is not obtained taking into account cumulative effect.In general, a new and comprehensive damage index needs to be proposed by combining the stiffness damage index and Park-Ang damage index.
Referring to (2), a similar stiffness damage index is proposed as where  0 is the initial stiffness of elastic stage and   is the current stiffness.
It is worth noticing that the damage index in (2) and (3) can also be transferred as the index presented by Dipasquale and Cakmak [5].
where  2 0 is the initial fundamental period and  2  is the current fundamental period.The equation above is suitable especially for damage assessment based on time-history analysis because the instantaneous fundamental period can be acquired conveniently.
According to (3), the relationship of current force   , the current displacement   and the damage index are as follows Based on the form of (5) The typical skeleton curve or capacity curve of reinforced concrete structure behaves nonlinear or multiple linear, as shown Figure 1, and the current force between the th point and the  + 1th point can be expressed as where A general global damage index for capacity curve or skeleton curve is obtained by substituting to (6) in ( 5) with the following result Compared to (2), the general global damage index presented in this paper has some advantages as follows: (a) the damage index varies from 0 to 1, rigorously, and has definite mechanical significance.(b) The damage index is calculated based on only once capacity curve that is, only one pushover analysis is carried out (c) the damage index is suitable for the capacity curve and the damage index of any displacement or any force can be obtained.

Damage-Index Adjustment Considering Hysteretic Energy.
During strong earthquakes, the deformation capacity of structures is reduced due to the cumulative damage caused by cyclic load.However, as the deficiency of pushover analysis, the cumulative damage is not embodied sufficiently and the corresponding modified coefficient is essential.A research on equivalent (reduced) ductility factors has been proposed by Fajfar and the cumulative damage effect is taken into account [13].
The total structure can be simplified as an equivalent single-degree-of-freedom ideal elastoplastic system and the following relations are found or defined:   =  2   ,   =   /  , and   =   /  , where   is the yield strength of the system,  is the mass,  is the natural frequency,   is the yield displacement, and   is the normal ductility.  is the monotonic displacement or current displacement and   is the monotonic ductility factor.
The parameter , which is a function of dissipated hysteretic energy, maximum displacement, and the natural frequency of the structural system, is introduced as follows: where  ℎ is the total hysteretic energy.The transformation can be obtained as Substituting to ( 8) in ( 1), the Park-Ang damage index is rewritten as Thus, the relationship between displacement ductility factor on cumulative damage in earthquake and the displacement ductility factor in monotonic load is as follows: For normal structure, the ratio   /  ranges from 0.65 to 0.85.The result of (11) can also be applied to (7) in order to compute the damage index considering dissipated hysteretic energy.For the capacity point which means the crossover point of the capacity curve and the ductility spectrum curve, the corresponding displacement ductility is viewed as   , but the damage status in pushover analysis is related to the monotonic ductility   , so the displacement on the capacity point is advised to be amplified according to (11) in order to represent the actual damage.

Number of Inelastic Cycles.
From another point of view, the cyclic demand spectrum for use in preliminary seismic assessment and design is necessary and Kunnath proposed a method based on estimating the number of equivalent cycles at a specified ductility [14].The methodology is determined using well established relationships between seismic input energy and dissipated hysteretic energy.If symmetric behavior is assumed, a commonly used idealization of the non-linear behavior is the elastic perfectlyplastic behavior shown in Figure 2. Assuming   is the yield force,   is the yield displacements and  max is the peak displacements.The energy dissipated by the yielding element per cycle is where   =  max /  , is the cyclic displacement ductility factor and  ℎ is an energy shape factor that depends on both the ductility and the shape of the hysteresis loops.For a nondegrading system subjected to   full cycles at the same ductility, the total energy dissipated is the product of   and   , and the following expression between the number of inelastic cycles and dissipated hysteretic energy can be derived: Considering the dynamic characteristics and seismic design parameter, the number of inelastic cycles can also be rewritten as where  ℎ = 1.13(  −1) 0.82 /  is the ratio of hysteresis energy to the total input energy   ,   is the force-reduction factor and the following expression is used where  is fundamental period and   is the characteristic period.
is the design spectral acceleration; V  is the equivalent input energy velocity, which is equal to the Fourier amplitude spectrum of the ground acceleration, and the typical expression is where V *  is the peak amplification factor for the input energy spectrum,   is the duration for the strong motion phase of the ground motion, and  is a parameter that characterizes the spectral shape of the input energy spectrum for periods longer than the predominant period of the ground motion.
It is demonstrated that the force-deformation characteristics of the system, the ductility-based force-reduction factor   , and the ground motion characteristics play a significant role in the cyclic demand imposed on a structure during severe earthquakes.

Coefficient of Hysteresis Energy Concentration.
The new and effective global damage indices are described as above.Furthermore, the deformation and hysteresis energy of the special weak stories is important because the structural weakest-part would suffer from severe damages due to its inability to withstand the different types of loads (lateral, vertical and moments) produced by the ground motion [15,16].The failure of the weak story can reduce the structural performance though other floors are merely in moderate damage stage.For irregular structures and complex structures, the damage assessment on weak floors is necessary besides the normal analysis carried out.
The coefficient of story yield shear is defined as   =   /  , where   is the yield shear force of a story in the The input energy of a multidegree-of-freedom system sustained to an earthquake depends on the structural mass and initial stiffness, in spite of the yield strength value.The story yield first almost is the story that the elastic displacement is the largest and the inelastic interstory drift is also greater than others.The hysteresis energy of the weak story is also greater than that of other stories.
Based on a large number of nonlinear time-history analyses considering various ground motions, fundamental periods, coefficients of story yield shear, and different number of weak stories, the fitting equation of the coefficient of story yield shear and the coefficient of hysteresis energy concentration is as follows: where  is the story number of structure and   ,   , and   are the parameters varying with .The statistical values for the structures with different stories are given in Table 1.
The coefficient of story yield shear in each story used to be over 0.3 even if the normal structures are subjected to severe earthquake.Thus the coefficient of hysteresis energy concentration generally decrease when the coefficient of story yield shear increases.Hence, the energy distribution can be obtained through calculating the coefficient of hysteresis energy and measuring the damage-degree by energy is achieved.

Damage Grade and Damage
Target.According to the discussion above and the damage indices presented, the global damage index and corresponding adjustment are viewed as an accurate damage-evaluation approach and the number of inelastic cycles can fully embody the integral hysteresis.These two damage indices represent the integral damage states in a reinforced concrete structure.On the other hand, the maximum interstory drift and the coefficient of hysteresis energy concentration can represent the local damage details.Hence, all the four damage indices are advised to be combined and used in damage assessment generally.
Taking into account the damage degree description in Table 2 and the empirical values from nonlinear analysis, the damage degree and the limit values for different damage indices are listed in Table 3.

Damage Assessment Model
Based on Fuzzy Set.In many practical cases, the structural damage is a fuzzy and vague state because there are no obvious and clear grading marks in the actual status; furthermore, the results of damage assessment based on different damage indices disagree in special cases [17].These problems may occur when data in Table 3 is used without identification and consideration.Thus, this is more scientific and objective for structural damage according to fuzzy set theory by introducing the concept of fuzzy subset and fuzzy membership [18][19][20][21].In classical set theory, the membership of elements in a set is assessed in binary terms according to a bivalent condition; that is,an element either belongs or does not belong to the set.By contrast, fuzzy set theory permits the gradual assessment of the membership of elements in a set; this is described with the aid of a membership function valued in the real unit interval [0, 1].Fuzzy sets generalize classical sets, since the indicator functions of classical sets are special cases of the membership functions of fuzzy sets, if the latter only take the value 0 or 1.The fuzzy decision fusion method realizes the integration fusion by using fuzzy transform principle and maximum membership principle and considering the various factors relevant to the evaluation objects.In the evaluation process, the factors constitute the factor set as  = { 1 ,  2 , . . .,   } according to different extent and the decision is expressed as the decision set as  = { 1 ,  2 , . . .,   }.
First, the single factor decision is made for the single factor   ( = 1, 2, . . ., ) in the factor set , then the membership degree   in the decision degree V  ( = 1, 2, . . ., ) is determined according to   .Thus, the single factor decision set for the factor   is obtained as   = { 1 ,  2 , . . .,   }, which is the fuzzy subset of the decision set.Hence, considering all the  factors, the general fuzzy decision matrix is constructed as where  is the fuzzy relationship between factor domain and decision domain and   is the membership degree.Multi-factor decision should consider the weighted grade for each factor in the fusion process.The fuzzy subset  = { 1 ,  2 , . . .,   } is the weighted set in factor domain , where   ( ≤   ≤ 1) is the weight value of the single factor   in the total decision.When  and  are determined, the fusion is carried out by fuzzy transformation  = I,  is equivalent fuzzy subset in the evaluation set , the element   ( = 1, 2, . . ., ) is the fusion result for grade V  , and I is the fuzzy synthesis operator.The membership function is the normal distribution function () =  −((−)/) 2 , where the parameters  and  could be determined according to the following principles: (1) the membership value is 1 when the factor is equal to the mean value, (2) the membership value of the top and bottom limitations in the universe of a fuzzy set can be adjusted to the membership value on the boundary in the universe of a fuzzy set is 0.5.

Damage Assessment Based on Fuzzy Set.
Considering various damage indices for reinforced concrete structures, the damage assessment should use one or multilevel fuzzy fusion.According to the importance and the accuracy of each index, the factor fuzzy vector can be determined by the expert experience or analytic hierarchy process, and the membership function is Gaussian model.Taking a reinforced concrete structure as an example, the global damage index, the number of inelastic cycles, the maximum interstory drift, and the coefficient of hysteresis energy concentration are the main parameters that accurately reflect the structural damage information and they are easy to obtain.Therefore, these parameters are selected in for fuzzy decision.
The weighted values for the parameters can be determined according to the experts grading method or the analytic hierarchy process (AHP).The AHP is a theory of measurement through pairwise comparisons and relies on the judgments of experts to derive priority scales [22].
Considering the importance and calculation accuracy of the grade index about damage state in the factor domain and the result from AHP, the weight value matrix of the damage index is A = [0.4,0.25, 0.2, 0.15].The structural damage grade is divided into 5 grades: intact, minor damage, moderate damage, severe damage, and collapse.The damage assessment and grade index are shown in Table 1.
After the fuzzy damage subset is determined, the membership value for the specific damage status in each damage grade can be calculated.In order to meet the requirements of performance-based design, the quantitative description of damage degree is needed.To highlight the dominance hierarchy, the general fuzzy damage index for decisionmaking is suggested as where   is the membership value of the -level damage,   is the median value in the -level damage and  is the comprehensive adjustment index, where it is recommended that  is taken as 2.

Damage-Assessment
Procedure.In summary, the procedure of damage assessment based on fuzzy set is as follows: (1) Determine the decision set as  = { 1 ,  2 , . . .,   } in order to describe the different damage grade and the get the equivalent fuzzy subset  in the evaluation set .
(3) The nonlinear response is obtained by pushover analysis and capacity spectrum method and the damage indices are calculated.
(4) The membership function is calculated and the fuzzy relationship matrix is obtained.
(5) The equivalent fuzzy subset in the evaluation set  is calculated by  = I.
(6) Compute the general fuzzy damage index (GFD) and determine the damage degree of the structure.

Case Studies
In order to verify the damage-assessment method based on fuzzy set, four three-dimensional reinforced concrete frames sustained to strong earthquake are selected to be analyzed by pushover analysis and nonlinear dynamic time history analysis using the three-dimensional nonlinear program CANNY [23].These 4 structures are as follows: 5-story regular frame, 15-story regular frame, 9-story eccentric frame and 4-story partially symmetric frame, as shown in Figure 3.The concrete compressive strength of the columns is 4 × 10 7 KN/m 2 , the compressive strength of the beams is 3 × 10 7 KN/m 2 , the yield strength of the main steel bars is 2.35 × 10 8 KN/m 2 , the reinforcement ratio of the columns and the beams are 2.2% and 1.8%, respectively.According to the corresponding design demand, the dead loads and live loads in all the structures are 1.0 KN/m 2 and 2.0 KN/m 2 respectively, and the loads are transformed as node force in the finite element models.The seismic fortification intensity is 8 and the field type is medium soil.The lateral load distribution of pushover analysis is adopted as the forms based on the stiffness distribution of each story.The capacity curves of the frames are obtained based on the capacity spectrum method.For the capacity spectrum method, an equivalent single-degree-of-freedom (ESDOF) system is a simplistic representation of the actual MDOF model, based on properties of the real structure; a bilinear idealization of the pushover curve can be obtained based on the principle of equality of energy.
For the time-history analysis, each group of waves has two orthogonal horizontal components and one vertical component, the three peak accelerations of each group are adjusted into 4.0 m/s 2 , 3.4 m/s 2 and 2.6 m/s 2 , respectively, and the duration is 15 seconds.
Taking Frame 1 and Frame 3 as an example, the capacity curves and demand curves are shown in Figures 4 and 5 and the global damage curves are in Figures 6 and 7 shown respectively.In the figures, the global damage index on the capacity point is adjusted according to (11).The performance indices are calculated based on ( 7), ( 11), (14), and (17).The results are listed in Table 4. Damage assessment based on fuzzy set is carried out and the membership function is established according to the method previously described and specified index in Table 3.For example, the membership function of the global damage index is shown in Figure 8.
For frame 1, the fuzzy relation matrix is calculated according to (17) and, combined with the weight value matrix, the fuzzy damage vector is Frame 1 is identified as the structure with severe damage according to the fuzzy damage vector, then collapse and finally moderate damage.The first story is weak story according to the results on interstory drift.Equation ( 4) is applied as the average damage index for frame 1 in seismic timehistory analysis, the value is 0.776, and the damage degree and locations are similar to the results of capacity spectrum analysis.The general fuzzy damage index (GFD) is 0.755 and the final assessment value quantized the damage status accurately.
Frame 3 is identified as the structure collapse after earthquake.The top story and the story with irregular stiffness are weak stories according to the results on the interstory drift.The average damage index according to (4) for seismic timehistory analysis is 0.856.The general fuzzy damage index (GFD) is 0.858 and the final assessment value is similar to the damage information from the time-history analysis for strong earthquakes.
The important performance results and the fuzzy damage vectors of all the four frames are shown in Table 3 and Figure 9.It is evident that the end deformation of the beam and column in frame 2 and frame 3 is obvious, the weak stories are in failure, and the whole performance is deficient and will collapse in strong earthquakes, so the maintenance and retrofitting should be carried out.The damage of frame 1 and frame 4 are both severe but the structures still have the capacity and the performance requirements are met.It is obvious that the seismic damage-assessment method based on fuzzy sets is accurate and effective.

Conclusion
The damage index is the necessary basis of damage assessment for reinforced concrete structures.Both the global damage index which represents the integrated performance and the local damage index which highlight the damage details are discussed and presented.The response parameters based on displacement and hysteresis energy are all vital especially for the structures subject to earthquake, and the analysis method is suitable for performance-based seismic design and assessment.Though some damage indices such as Park-Ang damage index, are comprehensive and acceptable for normal  damage assessment, it is unfortunate that these indices do not converge to unity and are unsuitable for pushover analysis if used directly.A global damage index based on multiple linear force-deformation curves in pushover analysis is presented to evaluate the integrate damage of reinforced concrete structure.The modified coefficient is provided considering the cyclic load and hysteresis energy.The number of inelastic cycles is introduced for the consideration of cyclic loads and hysteresis energy.The coefficient of hysteresis energy concentration is presented as a damage index which reflects the property and failure of the weak story.Thus, the global damage index, the number of inelastic cycles, the maximum interstory drift, and the coefficient of hysteresis energy concentration are all used as the damage indices of damage assessment.
In order to harmonize the different effects and weights, and reduce the uncertainty in the damage assessment, the fuzzy set theory is introduced, and the factor domain, the decision set, the fuzzy subset and the membership function are presented.The relation between the multiple damage  indices or factors and the fuzzy damage set is established by comprehensive fuzzy evaluation based on fuzzy mathematics.The fuzzy seismic damage-assessment method is verified through nonlinear analysis for four different structures and the corresponding results and assessment conclusions are accurate, especially when compared with the responses from nonlinear time-history analysis.
The performance-based fuzzy seismic damage-assessment method presented in this paper can be improved and developed in some aspects such as determining the weight values by analytic hierarchy process (AHP) method and considering the multiple dimensional ground motion and the torsional effect in the irregular structures.Furthermore, the performance-based assessment including the foundation and the nonstructural components can also be established based on fuzzy evaluation system.

Figure 3 :
Figure 3: Side elevation of the reinforced concrete frames.

Figure 6 :
Figure 6: Curves of global damage index of frame 1.

Figure 7 :
Figure 7: Curves of global damage index of frame 3.

Figure 8 :Figure 9 :
Figure 8: Membership grade of the global damage index.

Table 1 :
Statistical values for story yield shear and hysteresis energy concentration.  is the elastic shear force of the story, that is, the shear force for the same deformation if the structure is elastic.The coefficients of story yield shear in each story are equivalent for a regular structure; otherwise the values diverge.The story belongs to a weak story if the coefficient of story yield shear exceeds the threshold value.

Table 3 :
Damage degree and the limit value for different damage indices.

Table 4 :
The results of performance indexes of example structures.