Abstract

This paper introduces a probability-based damage state evaluation methodology for shear walls deriving from a data-driven calculation. A previously proposed damage quantification index, formulated in the time domain, which is capable of tracking damage progression based on the availability of structural monitoring data, is adopted here for the quantification of structural hysteresis damages. The probability of the structure lying in a specific damage state is determined on the basis of the derived damage index and the limit state definitions, yielding valuable information for postearthquake decisions. In this study, a database of the hysteretic behavior of 1,000 shear walls, considering different structural parameters, is generated utilizing OpenSees. Four limit states of seismic performance are defined based on material properties and in relation to the simulated stress and strain data. Accordingly, the exceeding probabilities of damage can be estimated by fitting statistical models to the damage index values grouped in terms of the axial load ratio. Followingly, an informative mapping, between the monitoring-derived damage index and postearthquake damage levels, is established considering structural uncertainty. Illustrative examples, including two shear walls subjected to cyclic loading and a seven-story shear wall slice subjected to a shaking table test, are investigated to show the capability and feasibility of the proposed performance evaluation method on damage state evaluation.

1. Introduction

Reinforced concrete (RC) shear walls form the most important lateral load resisting components in high-rise buildings, dominating residential and office building construction in urban areas [1, 2]. Strong earthquake events oftentimes induce severe structural damage, including concrete cracking and crushing of structural RC members and yielding or buckling of structural steel members [3]. Such a deterioration of structural condition, which may even develop progressively under a series of seismic events, can negatively affect structural safety and compromise future performance while exposing residents and societies to significant risks of life and economic loss. An illustrative instance of such an occurrence is the collapse of an RC shear-wall building during the 2003 Bingöl earthquake [4], while further, several RC frame-shear wall buildings experienced severe damage in the 2011 Christchurch earthquake prior to collapsing during aftershocks [5]. Under this prism, it becomes evident that a rapid and informative assessment of earthquake-excited structures can offer essential information for decision-making and provide crucial knowledge for both performance-based earthquake engineering (PBEE) and resilience-based seismic design [3, 6].

Currently, engineering demand parameters (EDPs), such as deformation, ductility, and hysteretic energy, derived on the basis of measured structural response, form the most widely adopted indicators in PBEE for characterizing and quantifying structural performance and deterioration [7], and determining damaged states according to guidelines [8]. However, the extremely limited availability of nonlinear structural response records, as earthquakes form rare (extreme) events, constitutes a primary hurdle in the process of deriving such data-driven indicators for the safety evaluation of shear wall structures in practice [9]. In view of this, traditional visual inspection still plays an important role in the postearthquake assessment of buildings. Guidelines relying on visual inspection are developed worldwide, such as in the USA [10], Italy [11], and New Zealand [12], albeit such processes being possibly unsafe, expensive, time-consuming, and subjective even for well-trained inspectors [13, 14].

Recent advances in vibration-based structural health monitoring (SHM) techniques [15] offer a feasible solution to shift the current inspection mode to one that exploits the availability of continuously gained information on in-situ structural performance [1619]. Such methods have been demonstrated to tackle a broad spectrum of infrastructure‐related issues, including structural identification, damage detection, and quantification [20, 21]. With the utilization of the monitoring data, changes in the vibration response and physical characteristics of instrumented structures can be observed at an early damage stage [22, 23]. Correspondingly, different kinds of data-driven damage-sensitive features (DSFs) have been formulated [24], including the popular modal property-based methods [2527] or more indirect signal processing-based methods [2830]. On the other hand, time-domain damage identification schemes exist, on the basis of which acceleration and displacement-based damage indices (DI) [3133] can be derived. Finally, traditional approaches, such as fitting of typical constitutive models, such as the Park–Ang model [34], have also been utilized for decades.

However, distinguished from the operational monitoring of in-service systems, one of the key challenges in SHM is the difficulty of gathering experimental or field test data from a structure in its damaged state, during which strongly nonlinear hysteretic phenomena and severe structural deterioration are experienced [35, 36]. Consequently, limited real-world demonstrations are found for data-driven damage detection in postearthquake damage assessment due to restrictions imposed from incomplete observations, linear system-based analysis algorithms, and limited records of nonlinear structural behavior [13, 37]. Currently, the determination of seismic performance levels, classified as immediate occupancy (IO), life safety (LS), and collapse prevention (CP), is still based on the qualitative criteria of typical EDPs (such as interstory drifts, inelastic element deformations, and element forces) [38], which may not fully consider the peak and cumulative damage effects. It is thus desirable to exploit data-driven monitoring methods and extend these with mappings that link the derived damage indices to damage states and performance levels. In this way, data can be harnessed to directly feed postearthquake assessment and decision support by capitalizing on knowledge from both the PBEE and SHM domains.

On the other hand, uncertainty, as an essential output of damage evaluation results, should be considered for the probabilistic evaluation of building performance under seismic loading. Uncertainty considerations lie at the core of procedures for damage and loss estimation in performance-based earthquake engineering (PBEE) [39]. A quantification of the confidence we can allot to damage evaluation results, in place of single deterministic values, is more informative for decision-making [40]. Aligning with the seismic fragility approach, which dominates the field of PBEE, a postearthquake assessment method adopting SHM data should yield a probabilistic evaluation of structural performance and corresponding damage states. In this respect, the DI values should be complemented with an estimate of the exceedance probabilities of specific limit states to more efficiently support postearthquake decisions.

Toward such a goal, a damage evaluation method for shear wall structures, which takes uncertainty into account, is proposed in this paper. Compared to those purely data-driven damage detection methods [41, 42], which considers uncertainty initially while proposing the method, the probabilistic metric adopted in this paper is employed directly to the calculated DI values of the shear wall components because the presented DI model is a model reference method in a data-driven manner. First, a simulated database of the hysteretic behavior of 1,000 RC shear wall components across different design parameters is generated by means of a finite element analysis. Material and section properties are utilized to define limit states and, correspondingly, determine associated damage states. The damage index (DI) values, which are calculated at each limit state for the whole database, are then used to drive seismic evaluation as an alternative form of EDPs. An informative mapping is established between the computed DI values and the determined seismic performance by clustering the data with structural design parameters and fitting the probabilistic models of the corresponding DI values. Two RC shear walls subjected to cyclic loading and a seven-story shear wall slice subjected to a shaking table test are adopted to demonstrate and validated the suggested data-driven assessment scheme.

2. Methodology

2.1. Design Framework

To evaluate structural damage under availability of SHM data, a data-driven method based on the so-called model-reference concept was previously proposed by the authors. The method relies on the availability of a reference linear system model, subjected to the same dynamic excitation, as the actual system. The associated structural damage is evaluated and further exploited for damage state classification. In this work, the previous design framework of the structural damage evaluation and seismic performance assessment is extended to consider uncertainty [43]. The specific steps to the proposed framework are outlined as follows:(1)The generation of a large synthetic structural response database: Over 1,000 shear wall FE models were built, based on the OpenSees platform, and was subjected to cyclic loading to extract records of structural hysteresis behavior(2)Critical damage state determination: The damage states in terms of seismic behavior are defined on the basis of prescribed structural drift thresholds suggested by earthquake engineering guidelines or empirical investigations. Correspondingly, the simulated material stress and strain data were utilized to determine four key points in the hysteretic behavior, as shown in the four light blue circles in Figure 1(3)Damage quantification for each damage state: Damage indices reflecting the four characteristic damage states, depicted in Figure 1, were calculated on the basis of the previously proposed damage model [44] and the derived hysteresis curves. This yields a database associating damage index values to damage states for more than 1,000 data points(4)Exceedance probability results for each damage state: The DI results were classified into groups on the basis of shear wall design parameters. Consequently, the exceedance probabilities for each damage state, on the basis of the associated damage index values, are obtained by fitting probabilistic distribution models.

The data-driven probabilistic assessment scheme summarized in this study delivers an informative mapping between seismic performance and monitored structural damage, enabling a monitoring-based approach to structural seismic performance assessment.

2.2. Data-Driven Damage Model

The proposed method allows for derivation of structural hysteresis curves, on the basis of available SHM data, which is then exploited to quantify structural dynamic damage [4547]. A brief introduction to and an essential description of the exploited data-driven damage model is presented below.

As inspired by the widely accepted Park–Ang model, the evaluation of the peak and cumulative damage effects is here exploited for seismic response assessment. This should be accomplished under appropriate normalization to allow for comparison. First, the peak deformation effect related to structural ductility is proposed to be computed as the minimum secant stiffness recorded during the earthquake response, which correlates to the maximum deformation point . This can be determined from the monitored shear wall hysteretic behavior and is defined aswhere is the effective stiffness, is the maximum deformation during the loading process and is the hysteresis force corresponding to .

Adopting the model-reference concept described in [44], a linear system is necessarily constructed as a reference model for comparing the difference between healthy and damaged status, subjected to the same dynamic excitation. At this point, it is useful to explain that certain assumptions are here put in place in terms of monitored quantities; namely, that (i) the force at the base of the wall is measured, often approximated through a measurement of the acceleration at the top of the wall/story, (ii) that the top relative displacement is recoded, and (iii) that a measurement of the seismic excitation is available. These quantities can be measured, although often nontrivially; however, these here serve as initiating assumptions, for proving the applicability of such a framework, which we intend to relax in future work. Theoretically, the preyield stiffness for the linear reference model of the shear wall can be determined by conducting linear regression for response data derived from low-amplitude forcing, when the structure behaves almost linearly. For the shear wall hysteretic behavior, the reference behavior of an assumed linear model, under the same excitation , can be calculated aswhere denotes the structural restoring force in the hysteresis curve.

According to a previous study of the authors [47], the tracking error , which can identify the deviation of the actual structural response [] to the reference (linear) response [], has been proved capable of tracking the occurrence of structural nonlinearity and can be adopted to evaluate the cumulative damage effect. The tracking error for shear wall hysteresis behavior can be defined aswhere and are the reference and actual displacement response in real time, respectively.

However, the determination of of equation (1) is only depending on the peak displacement point, which is bound to be influenced by measurement noise and model uncertainty when using real-world data. To reflect uncertainties, an additional variable, , is computed on the basis of a hysteresis loop that contains the peak response point . By performing a linear regression analysis of one such entire hysteresis loop, the value of can be determined, which offers a measure of the manifested nonlinearity and is conceptually equal to 1 for a linear model.

A data-driven index for the shear-wall damage evaluation has been proposed [44] as follows:where is the effective stiffness; is the linear reference stiffness of the healthy structural state; is the additional linear regression statics of the peak loop; is the root-mean-square (RMS) value of displacement tracking error; is the RMS value of the actual displacement response; , and are three model parameters controlling the weighting balance between the peak and cumulative damage effect with a suggested relationship of . The proposed DI model can evaluate nonlinearity for any system once the hysteresis behavior can be determined.

In the authors’ previous study, the proposed DI exhibits an approximately linear relationship with the reliable Park–Ang model regarding the hysteresis behavior of RC elements under the parameter set up of , , and , which makes it a feasible and recommended value for potential damage evaluation for RC elements and substructures. Accordingly, this parameter setup can also be changed corresponding to different materials or failure types to balance the damage effect. Meanwhile, this model is improved in this study by substituting the RMS value of reference response with the real response to assure the damage index value can converge to 1. Such improvement can bring more information for structural damage quantification and the threshold determination for shear wall components.

3. Database Generation

3.1. OpenSees Modeling

Finite element models for a large shear wall database, across numerous design parameters, have been established using the Open System for Earthquake Engineering Simulation (OpenSees) [48] for simulating the nonlinear shear wall behavior under cyclic loads. Each wall model is comprised of generic displacement beam-column elements consisting of five integration points and the widely-used fiber section model considering the shear effect, which has been shown to improve computational stability and efficiency [49]. The dispersed fiber section comprises four phases, namely unconfined concrete, confined concrete on both sides, and two longitudinal reinforcement steel bars. In terms of material constitute models, the concrete of the fiber section is endowed with the Concrete02 OpenSees material based on the Kent and Park constitutive model [50], while the Steel 02 OpenSees material was adopted to model steel reinforcement based on the Menegotto and Pino model [51]. The shear walls were considered fully constrained at the bottom with a non-zero-length rotating spring and free at the top. The detailed modeling information is illustrated in Figure 2.

Prior to generating this extensive synthetic database, two experimental studies of shear walls with public data were adopted to verify the feasibility and validity of the proposed modeling scheme and eliminate the model error. More specifically, the Dazio et al. [52] and Greifenhagen and Lestuzzi [53] specimens are selected from the shear wall database of Seismic Engineering Research Infrastructures for European Synergies (SERIES) and labeled as Test 1 and 2 in this paper. Each experimental study is complete with detailed information on design, construction, testing, measurement, and availability of related reports. The fundamental information about the two adopted shear walls is summarized in Table 1. As illustrated in Figure 3, the base shear and lateral displacement curves from the OpenSees simulations sufficiently approximate the data collected from the experimental study. For Test 1, Figure 3(a) indicates that the modeling method can simulate the structural envelope curve, with exception of the approximation of the softening behavior and concrete pinching phenomena. Figure 3(b) reflects a better approximation of the structural degradation and softening behavior for Test 2. Even though no modeling approach is expected to perfectly reconstruct diverse experimental results and different approaches lead to nonidentical local damages, the suggested modeling method does prove capable in sufficiently approximating typical shear wall behavior, which will here be exploited to assess the global damage of such walls.

To generate a large-scale database of shear wall seismic behavior, 1,000 shear wall models with diverse degrading and hysteretic characteristics were numerically analyzed. Six critical design parameters of the walls were selected according to authoritative design specification [1], including the axial load ratio, span-to-depth ratio, confinement zone ratio, longitudinal reinforcement strength, transverse reinforcement strength, and transverse reinforcement interval, vary within the suggested and empirical design values, such as the axial load value ranging from 0.1 to 0.35 and the span-to-depth ratio varies from 1 to 3 [54]. A reasonably designed shear wall model can be set up with the combination of these parameters, which were classified into different groups represented by different letters (A1, A2 to E1, E2). Table 2 shows the fundamental values of each group.

3.2. Damage State Determination

The definition and classification of structural damage states (DS) according to the simulated nonlinear structural behavior is accomplished by relating the computed DIs to individual structural performance level. Four DS ranges are defined on the basis of four critical limit states (LS), motivated by the commonly adopted performance level assessment in widely used seismic guidelines [8], as illustrated in Table 3. The shear wall is assumed to exhibit no damage (DS0) prior to the appearance of cracks (DS1). When the steel bar reaches yielding (LS2), while concrete spalling does not appear (LS3), the shear wall is classified into DS2, which corresponds to the performance level of immediate occupancy (IO). Once concrete crushing appears, we assume that DS4 is reached, which requires activation of collapse prevention (CP) mechanisms.

Since the 1,000 shear wall models are synthetic, the LS determination, which is needed for labelling and classification, is straightforward according to the criteria defined in the last column of Table 3. A typical structural hysteresis curve along with the four critical LS points for each shear wall, are shown in the top-right subplot of Figure 1. Figure 4 illustrates the stiffness and strength degradation factors in terms of the reached displacement and force for each LS, which is used to determine the distributions of structural performance situations for the four defined LSs. The stiffness degradation is defined by the instantaneous stiffness at LS4 divided by the linear stiffness at LS1, while the strength degradation is defined by the ratio of the force at LS4 divided by the maximum force during the loading. Figure 4 illustrates the degradation parameters corresponding to displacement and force for four limit states. It can be found that the scatters from LS1 to LS4 can be clearly distinguished and accordingly clustered. Meanwhile, the displacement increases from LS1 to LS4, as shown in Figure 4(a), and the force increases from LS1 to LS3 and decreases at LS4, as illustrated in Figure 4(b).

3.3. Damage Quantification

The DIs computed on the basis of the simulated hysteretic curves will be utilized to subsequently evaluate the structural damage at any point during the loading process. The DI values computed by equation (4) at the interfaces between LSs can be used as thresholds to define corresponding DSs. As evidence by equation (4), the DI computation contains a peak term related to stiffness degradation and a cumulative term computed from the hysteresis time history. Figure 5 shows that the values for these two individual contributing terms and the combined DI value increase from DS1 to DS4, as expected for a measure that is to be used for quantifying seismic performance. However, as illustrated in Figures 5(a) and 5(c), outliers are noted outside the 25th and 75th percentiles (marked in the box plots), which indicates that it would be doubtful to do data processing for the whole database directly.

The previous observation can be attributed to the fact that DI results may considerably vary across the considered design parameters, which are listed in Table 2. As the axial load ratio (variable A) plays a primary role in structural damage development [55], the original dataset is grouped into six categories per axial load ratio, which are labeled as A1 to A6 in the following study; each group comprises over 150 DI values corresponding to the four defined DSs. Due to the near-zero values of DIs in DS1, which is regarded as an operational status with effectively negligible damage occurring, Figure 6 only illustrates the DI values for DS2 to DS4 for these six groups.

A significant improvement in terms of classification is illustrated in the box plot in Figure 6, fewer outliers occur, apart from some irregular DI values in DS2 under high axial ratios, as illustrated in Figure 6(a), in DS3 under moderate axial ratios, as shown in Figure 6(b), and in DS4 for low to moderate axial ratios. In addition, Figure 6(c) indicates that the DI values corresponding to DS4 are more dispersed as the ultimate near-collapse state caused by extensive deformation or multiple damage modes leads to a variability of the recorded DI values.

4. Relationship of Damage and Performance

To generate an informative mapping between performance and structural damage, the corresponding thresholds of the data-driven DIs should be determined. On the other hand, a probabilistic damage evaluation result, which quantifies confidence in the estimation, is more informative for building owners or engineers than a deterministic damage index value. We here exploit the grouped DI values at four discrete DSs to explore their statistical features and develop probabilistic models for structural performance evaluation.

4.1. Probabilistic Model Selection

First, among candidate probabilistic models, the lognormal distribution is preliminary considered because of its broad use within the earthquake engineering community for the purpose of vulnerability analysis. Figure 7(a) indicates the resulting empirical and cumulative distribution functions (CDFs) of DS2 to DS4, and the probability density functions (PDFs) along with the data distributions are shown in Figure 7(b). Generally, the lognormal distribution model can fit the CDF and PDF effectively for DS2 and DS4, including an accurate peak estimation of the PDF, and provide an acceptable evaluation for DS3.

Subsequently, the classical Kolmogorov–Smirnov test (K-S test) was utilized by setting the mean and variance of the lognormal distribution equal to the grouped DI samples to implement the goodness of fit test. Figure 8 illustrates the K-S results on group A6; the empirical and cumulative distribution functions of DS2 to DS4 are shown, accompanied by the bounds with a 95% confidence interval for the K-S test. For DS2 and DS3, all data points meet the requirements and lie within the upper and lower bound, as shown in Figures 8(a) and 8(b). However, the DI values in DS4 are unable to meet these criteria, with violation of the lower bound of the K-S test at the mid-range of the DI, as illustrated in Figure 8(c).

We wish to therefore test the efficacy of the lognormal model versus alternate statistical models. The Akaike information criterion (AIC) and the Bayesian information criterion (BIC) form popular methods for comparing the adequacy of multiple, possibly nonnested models. Figure 9 illustrates the utilization results of AIC and BIC for the six groups considering DS2 to DS4 and different commonly adopted modeling options; namely, noraml [56], lognormal [57], beta [58], gamma [59], extreme [60], generalized extreme [61], logistic [62], and Weibull [63]. It is evident that the lognormal model delivers the lowest score, which is indicated to be the optimal one, in most cases and the best overall performance among all tested statistical models. The lognormal model is thus deemed acceptable for the grouped DI database estimation.

4.2. Damage State Evaluation

By adopting the lognormal distribution model, the cumulative distribution probability can be calculated for the six groups (A1 to A6) of DI values. A conversion is subsequently carried out by subtracting the exceeding probability of two adjacent DSs to acquire a normalized probability result of lying in each DS. Ultimately, the relationship between computed DI values and the exceeding probability of DS can be determined, as illustrated in Figures 10(a)–10(f).

Figure 10(b) plots the results for group A2, corresponding to an axial ratio of 0.15. The shear wall has a 35% probability of lying in DS2 (immediately occupied) and a 65% probability of entering DS3 (limit state) when the DI value is equal to 0.4. Once the DI value increase to 0.8, the probability of lying in DS4 (near collapse) reaches 90% with only a 10% probability of remaining in DS3 (limit state), which indicates severe structural damage to the shear wall. By contrast, Figure 10(d) illustrates that the shear wall in Group A4 is characterized by a probability of 50% of lying in DS3 and 50% of lying in DS4 for a DI value of 0.8, corresponding to the range of moderate to serious damages, respectively. Obviously, for the same damage value of 0.8, the probability of lying in DS4 has a big gap between group A2 with an axial ratio of 0.15 and group A4 with the ratio of 0.25.

These results can be used to develop a mapping between structural damage evaluation (through the computed DIs) and seismic performance levels for shear wall structures. For a rapid and practical postearthquake evaluation, the DI model can be computed, offering a quantification of the induced damage, and allowing for the corresponding exceedance probability of each discrete limit state to be determined, per Figure 10, by obtaining the basic design parameter. Then, the result containing DS levels and relevant probabilities can provide a beneficial reference for rapid postearthquake evaluation and decision-making because taking uncertainty into consideration is more practical in real-world applications.

5. Experimental Validation

5.1. Cyclic Loading Test

To investigate the application of the proposed method, the two shear walls introduced before and overviewed in Table 1 are utilized for seismic performance evaluation. Different design parameters, including the amount of boundary reinforcement, the axial force ratio, and the material properties, lead in different failure modes for the two walls. During the loading process corresponding to Test 1, flexure-shear cracks occurred with the contribution of the flexural and shear deformation. Failure is induced by fracture of the web and boundary reinforcing bars, which limit the deformation and energy dissipation performance, as illustrated in Figure 11(a). The wall in Test 2 suffered shear-related damage, with diagonal corner cracks initially observed and eventually suffering significant loss of shear capacity. The deformation capacity appears; however, not restricted, which favored ductility, as shown in Figure 12(a). The sudden concrete crushing led in termination of the test, indicating a flexural failure mode.

Once the hysteretic relationship between base shear and lateral displacement is determined, the proposed DI can be utilized to quantify structural damage. Treatment of windows of the response data over time allows tracking the DI evolution over time. By adopting four DI values at different experimental stages and the information using corresponding to the shear wall axial load ratio, the exceedance probability of DS for each critical point is determined. The results are summarized in Table 4.

For the shear wall Test 1, four critical levels are determined by test observation, as shown in Figure 11(a). Figure 11(b) illustrates the developed crack and damages after failure. The time for damage development between each state from DS1 to DS3 is stable, with the DI value increasing from near zero to 0.45, as shown in Figure 11(c). With the contribution of the confinement zone on both sides of the section, the wall shows good ductility as damage progression from DS3 to DS4 evolves over a substantial time period, reaching a DI of 0.78. After severe damage, such as core concrete crushing occurs leading to DS4, the shear wall of Test 1 fails at DI = 0.8.

Figure 12(a) shows the hysteretic behavior of the shear wall in Test 2 and Figure 12(b) illustrates the detailed failure performance after DS4. Since the shear walls have no confinement section, damage rapidly evolves from a slight to a moderate level, as shown in Figures 12(a) and 12(c). The DI value is retained low during the first two states with concrete cracking and steel yielding occurring. For DS3, the damage is underestimated to some extent with a computed DI value of only 0.2 in Figure 12(c), which is not conservative in terms of damage assessment. However, the wall shows favorable ductility, which provides enough developing time from DS3 to DS4 (near collapse) with DI increasing to 1.0. The evolution of the DI value over time, shown in Figures 11(c) and 12(c), illustrates that the suggested indicator ensures adequate damage tracking performance for shear wall components.

The estimated exceedance probability results corresponding to the four states determined by experimental observation are summarized in Table 4 for two shear walls. Generally, the probability results of the shear wall in Test 1 indicate a better correlation with the four actual structural performance levels related to experimental observations. When the shear wall enters DS2, the 55% probability of entering DS2 and a 20% probability of entering DS3 are computed, indicating a margin of safety. The core concrete crushing causes the wall to enter DS4 (near collapse), corresponding to an estimated 58% probability of entering DS3 and 42% for DS4. The overall evaluation results are conservative and thus secure for this shear wall component, and some actual serious seismic damage is overestimated by the proposed DI in some cases.

For Test 2, the performance level evaluation results for DS1 and DS4 fit perfectly with the actual damage observation, as both hold a 100% probability. While, for the performance level DS2 (immediately occupied) and DS3 (limit state), the proposed evaluation method results in overestimation of damage, with a low assigned probability on DS2 and only a 7% possibility allotted on entering DS3. The proposed approach seems able to warn of and prevent structural damage and provide enough margins for moderate seismic damage.

5.2. Shake Table Test

In this section, the feasibility and practical implementation of the damage state determination method with the damage evaluation model is further analyzed by adopting the shake table test of a 7-story RC shear wall slice available from DesignSafe [64]. The building was 19.2 m in height and consisted of a 3.65m-wide rectangular web wall, which provides lateral resistance along the applied direction of ground motion. The building was installed on the shake table with a dense array of sensors, including 139 accelerometers, 88 displacement transducers, and 314 strain gauges [65]. Figure 13 illustrates the experimental setup and the structural damage observation after the sequential tests of EQ1 to EQ4, with an increasing intensity of PGAs of 0.15 g, 0.27 g, 0.35 g, and 0.91 g, respectively. The collected data was adopted for reconstructing the illustrated hysteretic curve, as shown in Figures 14(a)–14(d), and computing the associated critical EDPs for the structure. The controlling parameters of the damage model were assigned as , , and, . Table 5 shows the damage quantification results when adopting the proposed damage index; these present a similar trend to the observed drift ratio and recorded concrete strain, which increase from EQ1 to EQ4 correlates to a developing structural damage.

The axial ratio should first be determined with geometric and loading information to determine the probability for each damage state. The structure herein can be classified into group A1, with a calculated axial load ratio of 0.1007. Subsequently, the calculated DI values in Table 5 and the relationship from Figure 10(a) can be adopted to determine the structural damage state under four earthquake cases. The heatmap shown in Figure 14(e) shows the relationship between the development of structural damage states and the increasingly excitation from EQ1 to EQ4. As illustrated, for the case of EQ1 corresponding to a low-amplitude excitation, the shear wall structure is attributed a probability of 89.3% to stay within DS1. In terms of EQ2 and EQ3, the major damage state is evaluated as DS2 with a probability of around 0.57; however, the EQ3 presents a larger probability for DS3. For the most severe earthquake case (EQ4), the result assigns a 97.58% possibility for DS3, which matches well with the structural damage development and the experimental observations such as concrete spalling and bar yielding, as well as the structural hysteresis behavior. Promisingly, the damage quantification and state probability determination methods interpret the dynamic monitoring data deeply with a clear correlation to the seismic evaluation, regarding the inherent structural uncertainty under nonstationary excitations.

6. Conclusion

A probability-based damage assessment method for shear walls is presented in this paper to determine the structural damage states by adopting a previously proposed data-driven damage tracking index. Four major steps are described in the proposed the methodological framework for seismic assessment: (1) the database generation of target structural behavior; (2) critical damage state determination; (3) data-driven damage quantification; and (4) database-based statistical mapping establishment between structural damage and seismic performance.

The structural hysteretic behavior database is first generated to derive and express mechanical knowledge of the target structural behavior. Numerical simulations for over 1,000 RC shear walls are performed by means of a FE analysis to extract the resulting hysteresis behavior and support the classification of critical DSs. The suggested monitoring-based DI can be utilized to quantify the structural damage at four critical LSs for each shear wall. By fitting the lognormal model on the grouped data divided by the axial load ratio, the exceedance probabilities of DSs corresponding to DI values can be established. The regressed results can be regarded as thresholds of the DI value corresponding to seismic performances considering structural uncertainty.

Two RC shear walls subjected to cyclic loading and a seven-story shear wall slice, subjected to a shaking table test, are adopted as component-level and system-level experimental studies, respectively, to demonstrate the practical capability of the proposed method. The DI model can accurately track and evaluate the structural damage development over time, with a favorable correlation with the actual hysteretic behavior and experimental observation. The DS and the exceedance probabilities, corresponding to four specified structural damage levels can be further determined by means of the presented structural performance evaluation method.

This study presents a valuable structural seismic performance evaluation method for shear wall structures, resulting in an assessment of the probability that the structure enters a specific DS. Such an informative mapping between structural monitoring data and performance level offers valuable knowledge for rapid postearthquake evaluation and decision-making, which would be beneficial for structural recovery and resilience. Noticeably, the proposed method can be adapted to the structural or substructural level once the hysteresis behavior can be obtained, which demonstrates strong application capabilities.

Despite the good evaluation performance, limitations remain to be addressed and need further investigation in future research, such as the parameter selection method for the DI model. Future work will focus on (1) the utilization of system identification and response estimation methods to perform story-level damage evaluation; (2) the extension of the current method to account for diversified structural failure types, such as shear-flexure coupling and further expansion to different structural types, such as bridges and steel structures; and (3) the clarification of uncertainty factors and the corresponding effect, such as model error.

Data Availability

Some or all of the data, models, or code that support the findings of this study are available from the corresponding author upon reasonable request

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was sponsored by the National Natural Science Foundation of China (Grant nos: 51878483 and 52278312), the Shanghai Rising-Star program (Grant no. 20QC1400700), the Shanghai Qi Zhi Institute (Grant no. SYXF0120020109), the Peak Discipline Construction Project of Shanghai (No. 2021-CE-03), and the Shanghai Municipal Science and Technology Major Project (2021SHZDZX0100).