# Application of Variance Analyses Comparison in Seismic Damage Assessment of Masonry Buildings Using Three Simplified Indexes

Authors: Qiwang Su; Gaochuang Cai; Hervé Degée
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

The reasonable assessment of potential damage type of masonry structures in seismic-prone zone is very significant to strengthen existing masonry structures and guide the construction of the new building. The primary objective of the study is to propose and determine a reasonable assessment index to predict the damage type of masonry structures in different seismic intensity zones using the survey results of the 2008 Wenchuan earthquake and variance analyses comparison methods. Three potential theory assessment indexes are considered in the evaluation of damage of masonry structures, including wall density indexIw, strength index Isq, and combined index Isd. In order to compare the feasibility of the three indexes, One-way analysis of variance and Scheffé’s method were used for in-depth discussion. Based on the proposed assessment indexes, further analyses and recommendations were provided. Results show the combined index Isd has a high potential to predict the damage levels of masonry structures. Based on the study, several recommendations were provided for the masonry structures in seismic-prone zones.


---

## Body

## 1. Introduction

Earthquake is considered as one of the most serious natural disasters that cause immeasurable damage to the building structures and our living environment. China is one of the countries that has a number of seismic active zones in the world and frequently suffers from catastrophic earthquake such as the 2008 Wenchuan earthquake (8.0 Mw) which has caused 69197 dead and 374176 injured people with 18222 listed as missing as of July 2008 [1]. National Police Agency of Japan [2] has also confirmed that there were 15894 dead, 6152 injured, and 2562 homeless people during the 2011 Great East Japan Earthquake (9.0–9.1 Mw). The World Bank has estimated that the economic cost caused by this earthquake was US$235 billion, making it the costliest natural disaster in the world history [3]. In the seismic-prone regions, it had been a consensus that the main seismic behaviour of building structures should be predicted precisely as much as possible. This is not only to understand the potential damage of new structures to guide our design works, but also in order to provide a series of technological supports to strengthen existing structures to resist a future coming earthquake. On the other hand, as one of common structural types, masonry buildings are widely used in a number of countries such as Europe, Chile, and China [4, 5], for its advantages including energy saving, easy availability of materials, and low construction cost. Generally, masonry structure usually is made of brick and mortar and can be divided mainly into unreinforced masonry (URM), confined masonry (CM) and reinforced masonry (RM) structures [6]. However, there are also a number of masonry structures constructed without referring to any design codes or standards, such as in the rural regions of developing countries. Consequently, a reliable assessment method is very important to enhance the seismic safety of the structures and predict their potential damage. It is also very helpful and meaningful to retrofit and strengthen the existing masonry structures [7–9].In the past three decades, many efforts have been performed to understand the structural characteristic of masonry structures and to establish reasonable masonry material models and effective finite analysis methods for their simulation and design [7, 10–13]. These studies usually involve some small-scale experimental studies and numerical simulation based on various finite analysis models [14–16]. Even though quite a few numerical simulation analyses have presented good evaluation results, earthquake filed investigations are limited and expected. Among them, several studies have attempted to develop simple seismic safety assessment indexes for URM structures [17–20]. For instance, Lourenço et al. [17] proposed three simplified indexes to assess the safety of the historical URM buildings under seismic loads. On the basis of their study, a few studies have been reported to make some modifications to these three indexes for URM or to present some combined analysis methods using the three indexes and finite element method to characterize CM structures [21, 22]. Besides, some vulnerability assessment methods based on uncertain structural systems analysis have been developed recently [18–20]. A few methods have been proposed which are suitable for both URM and CM structures. Moreover, due to the lack of seismic field survey data, there are a few studies (e.g., [17, 21]) that involved the relationship between assessment indexes and actual damage of URM/CM structures. In addition, several previous studies [23, 24] have confirmed that the confinement members of the structural masonry walls in CM structures can prevent effectively the collapse of masonry wall and then improve the seismic resistance of the whole structure system. The enhancement effect of the confinement elements on the deformation and resistance of masonry wall has been illustrated which then improves the seismic behaviour of the whole structure, as in the terms of structural ultimate failure modes. One more important thing is that up to now the primary assessment procedures of masonry structures in most current seismic design codes such as the ones in China and Latin America are established only on the basis of the structural strength design theory which was used to design them. It is suspectable and needs to be improved for enhancing the vulnerability of masonry structures.Therefore, the primary objective of the study is to propose reasonable assessment indexes to predict the potential damage type of masonry structures under seismic effects using the survey results of a number of masonry buildings in the 2008 Wenchuan earthquake. For this, some potential theory assessment indexes are considered including wall-related index, strength-related index, and a combined index considering the effect of confinement elements on masonry walls. Detailed mathematical statistic comparison methods are used to determine which index is more appropriate to assess the damage of masonry structures. Based on the proposed assessment index, some further analyses and recommendations are provided.


## 2. Damage Categories of Masonry Structures during an Earthquake

Generally, due to the fact that the thickness and openings ratio of masonry wall at each floor are identical, the stiffness of masonry structural members at each floor is initially the same. When the structure is subjected to lateral shearing force caused by earthquake effects, the main damage of masonry buildings usually takes place at the first storey. This behaviour has been confirmed by the previous field surveyed results and other researcher [25]. Therefore, the reasonable evaluation of the damage level of bottom walls is considered as a very important process to assess the potential damage degree of whole masonry building under earthquake loads. The damage degrees are determined and divided into four levels according to the damage situation of the wall pieces at the first floor, including fine cracks (width ≦ 1.5 mm), large cracks (width ≧ 1.5 mm), severe cracks (width ≧ 3.0 mm), and collapse [26]. The details of these damage categories are listed in Table 1 while some representative damage degrees are presented in Figure 1.Table 1
Damage categories of masonry wall.
Damage categories Damage description: cracking and collapse Treatments and measures Slight or no damage No obvious damage occurs in any wall pieces, or the number of the wall pieces with small cracks is less than 50% of all wall pieces in the seismic direction Small-scale repair such as surface repair  Moderate damage The number of the wall pieces with small cracks is more than 50% of the total walls in the seismic direction or the number of the wall pieces with large cracks is less than 50% of the total wall, or the number of wall pieces with severe cracks is under 10% of the total one Large-scale repair including partial change  Heavy damage The number of wall pieces with large cracks is more than 50% of all wall pieces in the seismic direction, or the number of wall pieces having severe cracks ranges from 10% to 50% of the total walls in the seismic direction Total/partial reconstruction  Collapse damage The number of wall pieces having the severe cracks, broken or collapse is more than 50% of the total walls in the seismic direction; or total collapse of building structure Total demolition and reconstructionFigure 1
Damage of bottom masonry walls during the 2008 Wenchuan earthquake.


## 3. Potential Damage Assessment Indexes for Masonry Structures

The earthquake-resistance performance of masonry buildings significantly depends on the lateral resistance capacity of masonry walls at the earthquake direction. This means that adequate walls should be provided in each potential earthquake direction of building structures. On the other hand, however, the seismic behaviour of CM structures consisting of confined members such as RC tie columns and beams is significantly affected by their seismic behaviour and strength properties of the masonry walls in the earthquake direction. Therefore, the reasonable assessment indexes to predict the potential damage of masonry structures should take into account the structural strength or the effect of confinement elements or both of them. In this study, three potential assessment indexes are considered to evaluate the damage degree of masonry structures, that is, wall density index in seismic directionIw, strength index Isq, and combined index Isd considering both strength and the effect of confinement elements. The section introduces a summary of their definitions and relevant concepts.


### 3.1. Wall Density per Unit FloorIw

According to the field survey obtained from one of the authors of the paper, the more masonry walls at each floor there are, the higher seismic resistance capacity the masonry structure can resist, especially for the URM structures. In this study, as reported by other researchers [16, 17], wall density per unit floor has been used to quantify the number of walls at each floor of the structure. Due to the fact that the masonry walls and tie columns in masonry structures are the main shear resistance members of masonry structures, the wall density per unit floor index (Iw) is a ratio of the total cross-section areas of all walls at ith direction to the total area of all the floors in structures, which is given by (1)Iwi=Awi+AcinAft,where Awi is the total cross-section areas of all walls in the ith direction; Aci is the total horizontal cross-section areas of all RC tie columns at the ith direction, which will be taken as zero for URM structures; Aft is the total area of each floor of masonry structures; and n is the number of storeys.During the calculation of this nondimensional indexIw, all masonry walls and RC tie columns are considered as the main members to resist the seismic shear forces of structures. Meanwhile, RC tie columns are roughly taken as a part of masonry walls in this process. It is very suitable and helpful to some masonry buildings which lack the detailed information of used construction materials’ strength. In order to predict the potential seismic damage levels of masonry walls, it is necessary to analyse and calibrate the relationship between Iw and the actual damage levels obtained from field survey data of masonry buildings.


### 3.2. Strength IndexIsq

In theory,Iw might lowly evaluate the seismic safety and damage of the masonry structures using high strength masonry materials. Therefore, the strength index of all vertical members at each floor could be an alternative assessment method that considers both wall density and the mechanical properties of wall materials. For example, the lateral resistance capacity of the vertical members of masonry structure at the base floor affects significantly the seismic resistance of the whole masonry building. Therefore, in low-rise wall bearing structure system, referring to the study result of Japan Building Disaster Prevention Association [27], for masonry structure, the strength index Cj of the vertical members (wall and tie column) in certain direction is expressed as(2)Cj=Qj∑w.The equation shows that the factor Cj is determined mainly by the used material strength and the horizontal cross-section of the calculated buildings. Therefore, for the wall and tie columns that can be considered as a group in masonry structures, the structural strength factor Isq in the ith direction is recommended in the study and given by (3)Isqi=Cj=Awiτw+Aciτcw∑Af=Awiτw+AciτcwnAft,τw=ξNfvk,sw.Isqi is strength factor in the ith direction; Qj is lateral resistance of vertical members in group j and in the ith direction, calculated as τjAj; τj is standard shear strength of the vertical members of group j in the ith direction; Aj is the effective area of vertical members of group j at the calculation floor; τc is shear strength of tie columns in masonry structure; τw is shear strength of masonry wall; ξN is effect factor of normal shear stress of masonry structures; fvk,sw is standard shear strength of masonry wall; w is gravity load per unit area of masonry structure, which was taken as 12000 Pa; ∑Af is the total area of the calculation floor and the ones over this floor.


### 3.3. Combined Assessment IndexIsd

In theory, wall density indexIw and strength index Isq can evaluate the damage of URM structures well; however, they are not good for CM structures. Previous research (e.g., [23]) reported that the enhancement effect of confinement elements of CM buildings should be taken into account when predicting the seismic damage of this confined masonry structures. The two indexes Iw or Isq may be improved further into a more effective index to assess the potential damage of masonry structures including URMs and CMs when considering the enhancement effect of confinement elements. Therefore, the study proposes a combined index considering masonry structural strength and the enhancement effects of confinement elements together. Using the above proposed strength index and referring to the seismic evaluation theory of existing RC structures in Japan Building Disaster Prevention Association [28], a combined index Isd in the ith direction considering the effects of strength and confinement’s effects can be defined as(4)Isdi=RiIsqi.Comparing with URM buildings, CM structures usually include RC tie columns and beams which both have a significant effect on the seismic behaviour of the total structures, in particular in structural deformation and lateral resistance performance. As a consequence, considering the effect of RC tie columns, tie beams, and rigid slab on total structural seismic behaviour, an assessment factor is proposed to evaluate the enhancement effectiveness of the confinement elements in confined masonry buildings, which is called the enhancement factor of confinement elements in theith direction Ri and expressed as follows:(5)Ri=γciγbiγsi,where Ri is the enhancement factor of confinement elements in confined masonry structures in the ith direction; γci, γbi, and γsi are three coefficients considering, respectively, the impacts of tie columns, tie beams, and rigid floor slabs on the structural performance of the whole structure in the ith direction. To understand the action of the three kinds of confinement on the masonry wall, a confined masonry structure was supposed to be put upside down. If no confinement values are set, the whole masonry system could collapse or just have limited seismic capacity. Therefore, in this study, they are defined simply as the corresponding volume ratios to masonry wall and expressed by (6)–(8). For URM buildings, γci and γbi both are equal to 1.0. Meanwhile, the value of γsi is also equal to 1.0 when precast RC slabs without strong connection are used in floor because of the insufficient of structural integrity. However, if the precast RC slabs have a good connection such as using concrete reinforced with hooked/welded steel connection, the value of γs should be calculated as per (8).(1)
γci is the enhancement ratio of RC tie columns to masonry wall:(6)γci=1+n1HciAcin1HciAci+HwiAwi.(2)
γbi is the enhancement ratio of RC tie beam to masonry wall:(7)γbi=1+n1∑j=1mLqi,jAqi,jn1∑j=1mLqi,jAqi,j+∑j=1mLwi,jAwyi,j.(3)
γsi is the enhancement ratio of RC floor slab to the masonry structures:(8)γsi=1+n1∑j-1mAb,jhb,j∑j=1mVj,whereHci and Hwi are the total vertical height of tie columns and of masonry walls in the ith direction, respectively; Aqi,j is the total vertical cross-section areas of tie beams at the jth floor in the ith direction; Lqi,j is the average horizontal length of tie beams at the jth floor in the ith direction; Awyi,j is the total vertical cross-section areas of masonry walls at the jth floor in the ith direction; Lwyi,j is the average horizontal length of masonry walls at the jth floor at the ith direction; Ab,j is the in-plane area of floor slab (which has been reliably connected) at the jth floor; hb,j is the thickness of floor slab at the jth floor; m is the total of floors; Vj is the total volume of the jth floor which is calculated using the height (hj) and in-plane area (Af,j) of the jth floor, respectively, as follows:(9)Vj=n1Ab,jhb,j+hj-hb,jAf,j.Meanwhile,n1 is the shear modules ratio of concrete element to masonry wall in masonry structure. This factor can be used to explain the different structural behaviour and damage degree of the buildings with the same structural geometry but with sensible difference regarding mechanical ratio of rebar on nodes. In order to calculate the factor n1, the enhancement effect of confined elements on seismic behaviour of CM buildings should be understood clearly. For example, to understand the effect of the tie column elements to the seismic performance of the whole buildings, masonry building is reversed vertically. If there are no vertical confinement elements such as URM, the vertical connection strength of the building will be determined by the axial ultimate tensile strength of wall. In CM building, however, tie columns significantly enhance the vertical integrity of the structures. As a result, γci, γbi, and γsi also represent an increase ratio of the vertical, horizontal, and spatial connection strength. The factor n1 can be expressed as(10)n1=ftkc-ftkmftkm,where ftkm and ftkc are the characteristic axial tensile strengths of masonry units and tie columns, respectively. According to the axial tensile calculation theory of reinforced concrete and the definition of ftkc, a simplified mechanical equilibrium in CM building is obtained as follows:(11)ftkcAc=fykAs,where fyk is the tensile strength of steel rebar in tie column; As is the sectional areas of steel rebar in per tie column (almost the same in the same building); Ac is the sectional area of per tie column; According to (11), ftkc is obtained which is then used to calculate the factor n1. For ftkm that is determined by the tensile strength of mortar which can be gained by the related codes such as in [4], the ratio of ftkc to ftkm is close to the ratio of ultimate shear strength of concrete to masonry.


## 3.1. Wall Density per Unit FloorIw

According to the field survey obtained from one of the authors of the paper, the more masonry walls at each floor there are, the higher seismic resistance capacity the masonry structure can resist, especially for the URM structures. In this study, as reported by other researchers [16, 17], wall density per unit floor has been used to quantify the number of walls at each floor of the structure. Due to the fact that the masonry walls and tie columns in masonry structures are the main shear resistance members of masonry structures, the wall density per unit floor index (Iw) is a ratio of the total cross-section areas of all walls at ith direction to the total area of all the floors in structures, which is given by (1)Iwi=Awi+AcinAft,where Awi is the total cross-section areas of all walls in the ith direction; Aci is the total horizontal cross-section areas of all RC tie columns at the ith direction, which will be taken as zero for URM structures; Aft is the total area of each floor of masonry structures; and n is the number of storeys.During the calculation of this nondimensional indexIw, all masonry walls and RC tie columns are considered as the main members to resist the seismic shear forces of structures. Meanwhile, RC tie columns are roughly taken as a part of masonry walls in this process. It is very suitable and helpful to some masonry buildings which lack the detailed information of used construction materials’ strength. In order to predict the potential seismic damage levels of masonry walls, it is necessary to analyse and calibrate the relationship between Iw and the actual damage levels obtained from field survey data of masonry buildings.


## 3.2. Strength IndexIsq

In theory,Iw might lowly evaluate the seismic safety and damage of the masonry structures using high strength masonry materials. Therefore, the strength index of all vertical members at each floor could be an alternative assessment method that considers both wall density and the mechanical properties of wall materials. For example, the lateral resistance capacity of the vertical members of masonry structure at the base floor affects significantly the seismic resistance of the whole masonry building. Therefore, in low-rise wall bearing structure system, referring to the study result of Japan Building Disaster Prevention Association [27], for masonry structure, the strength index Cj of the vertical members (wall and tie column) in certain direction is expressed as(2)Cj=Qj∑w.The equation shows that the factor Cj is determined mainly by the used material strength and the horizontal cross-section of the calculated buildings. Therefore, for the wall and tie columns that can be considered as a group in masonry structures, the structural strength factor Isq in the ith direction is recommended in the study and given by (3)Isqi=Cj=Awiτw+Aciτcw∑Af=Awiτw+AciτcwnAft,τw=ξNfvk,sw.Isqi is strength factor in the ith direction; Qj is lateral resistance of vertical members in group j and in the ith direction, calculated as τjAj; τj is standard shear strength of the vertical members of group j in the ith direction; Aj is the effective area of vertical members of group j at the calculation floor; τc is shear strength of tie columns in masonry structure; τw is shear strength of masonry wall; ξN is effect factor of normal shear stress of masonry structures; fvk,sw is standard shear strength of masonry wall; w is gravity load per unit area of masonry structure, which was taken as 12000 Pa; ∑Af is the total area of the calculation floor and the ones over this floor.


## 3.3. Combined Assessment IndexIsd

In theory, wall density indexIw and strength index Isq can evaluate the damage of URM structures well; however, they are not good for CM structures. Previous research (e.g., [23]) reported that the enhancement effect of confinement elements of CM buildings should be taken into account when predicting the seismic damage of this confined masonry structures. The two indexes Iw or Isq may be improved further into a more effective index to assess the potential damage of masonry structures including URMs and CMs when considering the enhancement effect of confinement elements. Therefore, the study proposes a combined index considering masonry structural strength and the enhancement effects of confinement elements together. Using the above proposed strength index and referring to the seismic evaluation theory of existing RC structures in Japan Building Disaster Prevention Association [28], a combined index Isd in the ith direction considering the effects of strength and confinement’s effects can be defined as(4)Isdi=RiIsqi.Comparing with URM buildings, CM structures usually include RC tie columns and beams which both have a significant effect on the seismic behaviour of the total structures, in particular in structural deformation and lateral resistance performance. As a consequence, considering the effect of RC tie columns, tie beams, and rigid slab on total structural seismic behaviour, an assessment factor is proposed to evaluate the enhancement effectiveness of the confinement elements in confined masonry buildings, which is called the enhancement factor of confinement elements in theith direction Ri and expressed as follows:(5)Ri=γciγbiγsi,where Ri is the enhancement factor of confinement elements in confined masonry structures in the ith direction; γci, γbi, and γsi are three coefficients considering, respectively, the impacts of tie columns, tie beams, and rigid floor slabs on the structural performance of the whole structure in the ith direction. To understand the action of the three kinds of confinement on the masonry wall, a confined masonry structure was supposed to be put upside down. If no confinement values are set, the whole masonry system could collapse or just have limited seismic capacity. Therefore, in this study, they are defined simply as the corresponding volume ratios to masonry wall and expressed by (6)–(8). For URM buildings, γci and γbi both are equal to 1.0. Meanwhile, the value of γsi is also equal to 1.0 when precast RC slabs without strong connection are used in floor because of the insufficient of structural integrity. However, if the precast RC slabs have a good connection such as using concrete reinforced with hooked/welded steel connection, the value of γs should be calculated as per (8).(1)
γci is the enhancement ratio of RC tie columns to masonry wall:(6)γci=1+n1HciAcin1HciAci+HwiAwi.(2)
γbi is the enhancement ratio of RC tie beam to masonry wall:(7)γbi=1+n1∑j=1mLqi,jAqi,jn1∑j=1mLqi,jAqi,j+∑j=1mLwi,jAwyi,j.(3)
γsi is the enhancement ratio of RC floor slab to the masonry structures:(8)γsi=1+n1∑j-1mAb,jhb,j∑j=1mVj,whereHci and Hwi are the total vertical height of tie columns and of masonry walls in the ith direction, respectively; Aqi,j is the total vertical cross-section areas of tie beams at the jth floor in the ith direction; Lqi,j is the average horizontal length of tie beams at the jth floor in the ith direction; Awyi,j is the total vertical cross-section areas of masonry walls at the jth floor in the ith direction; Lwyi,j is the average horizontal length of masonry walls at the jth floor at the ith direction; Ab,j is the in-plane area of floor slab (which has been reliably connected) at the jth floor; hb,j is the thickness of floor slab at the jth floor; m is the total of floors; Vj is the total volume of the jth floor which is calculated using the height (hj) and in-plane area (Af,j) of the jth floor, respectively, as follows:(9)Vj=n1Ab,jhb,j+hj-hb,jAf,j.Meanwhile,n1 is the shear modules ratio of concrete element to masonry wall in masonry structure. This factor can be used to explain the different structural behaviour and damage degree of the buildings with the same structural geometry but with sensible difference regarding mechanical ratio of rebar on nodes. In order to calculate the factor n1, the enhancement effect of confined elements on seismic behaviour of CM buildings should be understood clearly. For example, to understand the effect of the tie column elements to the seismic performance of the whole buildings, masonry building is reversed vertically. If there are no vertical confinement elements such as URM, the vertical connection strength of the building will be determined by the axial ultimate tensile strength of wall. In CM building, however, tie columns significantly enhance the vertical integrity of the structures. As a result, γci, γbi, and γsi also represent an increase ratio of the vertical, horizontal, and spatial connection strength. The factor n1 can be expressed as(10)n1=ftkc-ftkmftkm,where ftkm and ftkc are the characteristic axial tensile strengths of masonry units and tie columns, respectively. According to the axial tensile calculation theory of reinforced concrete and the definition of ftkc, a simplified mechanical equilibrium in CM building is obtained as follows:(11)ftkcAc=fykAs,where fyk is the tensile strength of steel rebar in tie column; As is the sectional areas of steel rebar in per tie column (almost the same in the same building); Ac is the sectional area of per tie column; According to (11), ftkc is obtained which is then used to calculate the factor n1. For ftkm that is determined by the tensile strength of mortar which can be gained by the related codes such as in [4], the ratio of ftkc to ftkm is close to the ratio of ultimate shear strength of concrete to masonry.


## 4. Three Proposed Indexes versus Actual Damage Degree of Masonry Structures

Through investigating 130 masonry buildings (mainly including CM) in the seismic intensity IX zone (modified Mercalli scale) during the 2008 Wenchuan earthquake, the relationships between the proposed indexes and the damage degree surveyed building are studied, respectively. In the inspected buildings, the tie columns have a width of 240 mm and a height of 240 mm; both are equal to the thickness of masonry wall. All tie columns usually were arranged wall joints, corners, and margin of openings. They generally were reinforced by 10 mm deformed steel rebar (about 235 MPa yielding strength) and confined by diameter 6 mm stirrup with a spacing of 200 mm per code at that time. These reinforcements have been connected well in the nodes of these inspected buildings, which guarantee the deformation and mechanical performance of RC confined members.(1) Wall Density per Unit Floor I w . Figure 2 shows the relationships between Iw and the actual damage categories of masonry structures at seismic intensity IX zone, respectively. Results present that this index distinguishes these different damage degrees of masonry structures well. When CM buildings are at the high seismic intensity zone, heavy damage and moderate damage occurred in the building with the indexes Iw of less than 2.0% and 2.5%, respectively. On the other hand, there are no masonry structures which had failed as collapse damage when the Iw of buildings is larger than 1.25%, both URM and CM buildings. In addition, as shown in the figure, the increase of Iw has resulted in an obvious decrease in the damage degree of these masonry structures and presents a clear difference between different damage. This indicates that index Iw could be applied to distinguish the potential damage type of masonry structures. Taking the above presented Iw levels as the critical values corresponding to different damage degrees, the simple wall density index Iw is very helpful to assess potential damage categories of CM structures, especially for the ones with limited recorded information.Figure 2
I w and damage categories of masonry structures at seismic intensity IX zone.(2) Strength Index I s q . In order to clearly present the relationship between damage degree and Isq of masonry buildings at same seismic intensity zone, two simplified factor are introduced, that is, x and y. Here, x is a ratio of total wall area in seismic direction to the total plane area of the calculated floor and its upper floors, which is defined as(12)x=Aw∑Af.On the other hand, y is a ratio of the total weight of the calculated floor and its upper floors to the total plane area of tie columns and walls in seismic direction, which represents the average axial compressive strength of these masonry walls in this direction, which is defined as (11). In this equation, the self-weight of per unit volume masonry wall is taken as a constant value (0.012 MPa) according to Chinese code [4].(13)y=0.012∑AfAw+AcinMPa.Therefore, using (12) and (13), (2) can be expressed in the form of(14)y=0.00840.012Isq+0.35xinMPa.According to the (12), as shown in Figure 3, the relationship between x and y can be attained to express the relationship damage degree and Isq of masonry buildings, which comes from the seismic survey data of the 2008 Wenchuan earthquake (in intensity IX zone). The results show that the critical values of the index Isq provide a good distinguishing method for assessing the potential damage degree of the masonry buildings. Based on this figure, the masonry structures for most CM buildings with collapse damage were observed when the values of Isq range from 0.2 to 0.4. Similarly, the heavy damage and moderate damage are, respectively, observed when the index Isq ranges from 0.4 to 0.6 and from 0.4 to 0.8, respectively. When the values of Isq range from 0.6 to 0.8, heavy damage was confirmed in some URM structures, even though the number of the kinds of URM is small.Figure 3
I s q and damage levels of the masonry structures in seismic intensity IX zone.(3) Combined Index I s d . To study the relationship between index Isd and damage level of masonry structures, the same procedures and field investigation data as the ones of indexes Iw and Isq are used. Comparing with the process to calculate these two indexes, more detailed information of masonry buildings are required for index Isd, including mainly the detailed information of tie columns, tie beams, and floor slab. The relationships between Isd and the damage categories of these investigated masonry structures are presented in Figure 4. There is a decline in the damage degree of the masonry structures as index Isd increases. In addition, the results plotted as Figure 4 shows that the damage categories of these masonry building are differentiated by index Isd more easily and clearly. For the collapse buildings, all the values of Isd are not greater than 0.6 and the values of the strength index Isq for these building also are less than 0.4. Similarly, for the CM buildings with heavy damage and moderate damage, the corresponding values of their index Isd are not greater than 0.8 and 1.02, respectively. The results imply that the combined index also presents a high potential to be applied to assess the damage levels of masonry structures in seismic intensity IX zone.Figure 4
I s d and damage categories of the masonry building in seismic intensity IX zone.


## 5. Detailed Mathematical Statistic Comparison

The above comparisons indicate that indexIsd presents a higher rationality and feasibility to predict the potential damage of masonry structures. However, more detailed comparison between these three indexes is expected, for it is good to support the relevant provisions and codes formation based on a more quantitative comparison. Meanwhile, for the view of large sample database, the appearance of the four types of seismic damage of masonry structures is approximately normal distribution. If considering all single samples with the same damage type as a sample group, the differences between these four groups should be reflected on each proposed assessment indexes. Therefore, the obvious difference or distinguishable degree of the assessment indexes can be used to evaluate the superiority of these assessment indexes. It should be noted that the results regarding URM structures in the paper may represent a limited behaviour for the number of URM buildings which is quite small in the seismic intensity IX zone of China. Form the view of mathematical statistics, in order to study whether the control variables (here, four seismic damage levels) had a significant impact on observed variables (three proposed indexes in this study), there are mainly two methods, that is, a single factor analysis and multivariate analysis.


### 5.1. One-Way Analysis of Variance (ANOVA) Method

In order to analyse the differences of the effects of the above three assessment indexes on the potential damage categories of masonry structure, a comprehensive statistic comparison method is applied in the study, which is called one-way Analysis of Variance (ANOVA) [29]. This method is a kind of univariate analysis, which majorly studies if the difference of the one of the control variables has a significant influence on the total observed variables. Based on the theory studies of Lawal and Famoye [30], the mathematical model of one-way ANOVA is explained simply as follows:(1)
Assuming the control variableA has k levels and each level has r samples, the jth sample xij with Ai control variable level is(15)xij=μi+εiji=1,2,3,…,k;j=1,2,…,r,where μi is the theory value of the Ai control variable and εij is random sampling error.(2)
Assumeu is the total theory value of control variable  A, which is expressed as follows:(16)μ=1k∑i=1kμi.(3)
So the additional effect of control variablesAi on the experimental results αi is obtained and the total value of αi is zero, given by (17)αi=μi-μ,i=1,2,…,k,(18)∑i=1kαi=0.(4)
Based on (16) and (17), xij is transferred as (19) and called mathematic model of the one-way ANOVA:(19)xij=μ+αi+εij,i=1,2,…,k;j=1,2,…,r.(5)
In the equation, the unbiased estimators of the total theory value of control variables (μ) and the additional effect of control variables Ai on the experimental results (αi) are(20)∧μ=x¯,∧αi=xi¯-x¯.If control variableA has an influence on observed variables, at least one of the values of αi is not zero. Otherwise, all of the values of αi are zero. The detailed analysis is called F-test analysis, whose completed steps are as follows:(1)
Firstly, two null hypotheses are proposed, that is,Haccept and Hobject.(2)
Then, calculate the selection test statisticF value, which is defined as(21)F=SSA/k-1SSE/m-k=MSAMSE.In this equation, m is total number of samples; SSA is the sum of squares between groups; SSE is the sum of squares in the groups; k-1 and m-k are the freedom degree of SSA and SSE; MSA is the sum of squares of between the average values groups; MSA is the sum of squares of in the average values groups. The standard value of F is from F-distribution having (k-1,m-k) of freedom degree.(3)
CompareF value and probability P value; when the observed F value is greater than 1, that means the control variables have a significant impact on the observed variables; otherwise the impact is not significant. Additionally, the larger F value is, the more significant the effect is.Based on the above, one-way ANOVA is performed to study the effects of the damage categories on the three proposed indexesIw, Isq, and Isd at different seismic intensity zone. Here, as an example, the one in seismic intensity IX zone is presented in the study. A statistical package for the social sciences statistical program (SPSS) was used and four damage groups were created in the one-way ANOVA study (group 1 collapse damage, group 2 heavy damage, group 3 moderate damage, and group 4 slight or no damage). Firstly, it was assumed a single experiment variable (main control variable) has k levels in the ANOVA, that is, the above damage degree. Based on this, the value of k equals 4. Therefore, the test level of overall significance is obtained using the following hypotheses:Hobject: α1=α2=α3=α4=0.Haccept: at least one of the values of αk is not equal to zero.Subsequently, the above hypotheses were performed and determined by comparing computedF value and the critical one in one-way ANOVA [31]. The critical values of F were obtained according to some statistic variables of specimen database. On the other hand, in some cases, the computed F values are close which means the difference is not obvious. For these reasons, in this study, the proportion of variance (η2) was also used in the comparison which was provided by McClave and Sincich [31], given by(22)η2=SSBetweenGroupsSSTotal,where SSBetweenGroups is the sum of squares of between groups and SSTotal is the sum of squares in total.


### 5.2. Multiple Comparisons: Scheffé’s Method

If the hypothesis is rejected as a result of one-way analysis of variance, it namely means that the various levels of assessment indexes are significantly different. However, this only shows that at least a difference between the two levels of assessment indexes is significant and it cannot say that there are significant differences between any two levels of the assessment indexes with the damage degree. For more than three assessment indexes, it often needs to determine where the specific significant differences are, that is, between which levels are the assessment indexes found. Therefore, if a refusal result was obtained from the analysis of one-way ANOVA, the further comparison should be finished for any difference between any two average groups from a number of average groups. This method is generally named multiple-comparison. The common method of multiple comparisons has Tukey–Kramer’s and Scheffé’s method. Tukey–Kramer method is better when only a fixed number of pairwise comparisons are to be made which will result in a more precise confidence interval. Scheffé’s method is more appropriate to the cases having a large number of comparisons, on the other hand, when there are many or all contrasts in samples. Therefore, in this study, Scheffé’s method was used in the following comparison.Firstly, the mathematical model of Scheffé’s method [32, 33] is illustrated as follows:(1)
When the evaluation result of one-way ANOVA is to reject the hypotheses, that is, the various levels of assessment indexes are significantly different, the difference between theith and jth groups can be compared via using a calculated S value which is defined as(23)S=max1≤i,j≤r⁡X¯i-X¯jMSE1/ni+1/nj.(2)
The fractile quantile ofS value (Sa) can be obtained from S-table or from the following equation:(24)Sa=r-1Fαr-1,n-r.(3)
Similar to the steps in one-way ANOVA, the steps of Scheffé’s method are as follows:
(a)
Firstly, a null hypothesis is proposed:
H0: μi=μj, i≠j, i,j=1,2,3,4,…,r.H1: H0 is not correct.(b)
CalculateSa based on (24) or from S-table and then calculate Tij which is defined as(25)Tij=Sani+njninj.(c)
CompareTij and the absolute difference value of dij between sample and average values. dij is defined as(26)dij=xi-xj.If the absolute difference value ofdij is greater than Tij, that means there is an obvious difference between the ith and jth groups; otherwise, their difference is not significant.


### 5.3. Comparison and Discussion

Using the above two comparison methods, the three indexes are compared in detail. In this study, only the results in seismic intensity IX zone are presented, as a representative result. The other comparison can be reached at the detailed research report from the authors [34].


#### 5.3.1. One-Way ANOVA

Table2 shows the one-way ANOVA comparison results of the assessment results of these indexes. There are four groups which were used in the analysis, that is, collapse, heavy damage, moderate damage, and slight or no damage. Results show that the critical values of F were both taken as 2.60 for the two indexes Iw and Isd according to the values of df [31]. The results also indicate that the calculated F value in the ANOVA is 29.74, which is larger than its critical value. Therefore, the hypotheses is rejected strongly which means index Iw has a significant effect on the damage categories of masonry structures. In other words, index Iw can be used to clearly distinguish or predict the potential damage categories of the masonry buildings in seismic intensity IX zone. Comparing with the results of index Iw, however, the computed F values of indexes Isq and Isd are better and reach 40.28 and 176.008, respectively. This implies that index Isd has a more significant influence on the damage categories of surveyed masonry structures. On the other hand, the significance degrees (Sig.) of these indexes all are equal to zero, meaning that they all have significant influences on the damage degree of masonry structures. Besides, the result shows that the value of η2 for indexes Iw, Isq, and Isd are 0.39, 0.47, and 0.79, respectively. This illustrates that the combined assessment index Isd is the best index to predict the potential damage categories of the masonry buildings in earthquake-prone zone. In summary, comparing other indexes, Isd has the largest computed F value, which means the index has the most significant influence on the damage levels of these masonry buildings.Table 2
Comparison of the results of one-way ANOVA (seismic intensity IX zone).
Indexes Sum of squares df Mean square CalculatedF value CriticalF value Proportion of variance (η2) Sig. Iw Between groups 3.25 E + 09 3 1.08 E + 09 29.735 2.6 3.98 E - 01 0 Within groups 4.92 E + 09 135 3.64 E + 07 Total 8.17 E + 09 138 Isq Between groups 3.629 3 1.21 40.282 2.6 4.72 E - 01 0 Within groups 4.054 135 0.03 Total 7.683 138 Isd Between groups 14.772 3 4.924 176.008 2.6 7.96 E - 01 0 Within groups 3.777 135 0.028 Total 18.548 138


#### 5.3.2. Scheffé’s Comparison

Table3 shows the results of Scheffé’s analysis comparison of the three indexes in seismic intensity IX zone. Firstly, regarding index Iw, the results present that there are significant difference between the first group (collapse) and other groups, because their probability P values are less than 0.05, a supposed significant level. Similar results are observed for the fourth group (slight or no damage). For the moderate and heavy damage group, however, they also present significant difference levels with other kinds of damage, except for each other, that is, the difference between them. In other words, they just can differentiate between the buildings with the collapse and slight/no damage groups, respectively. It indicates that Iw is not able to differentiate the buildings with moderate and heavy damage. For the assessment index Isq, similar results are confirmed to the ones of index Iw. The analysis results of the buildings with collapse and slight/no damage all show a significant difference with the corresponding three types of damage. However, comparing with the ones of index Iw, the P values of the moderate damage group were improved which has no obvious difference with the heavy damage. Regarding the combined index Isd, all damage groups present well a significant difference level with the other groups. The appearance probability of no difference between any two groups is zero which means this index effectively divides all kind of damage degrees including the moderate and heavy damage types for masonry structures in seismic intensity IX zone.Table 3
Comparison of the results of Scheffé’s analysis (seismic intensity IX zone).
Indexes ( I ) group ( J ) group Mean difference(I-J) Std. error Sig. 95% Confidence interval Lower bound Upper bound I w 1 2 - 3966.734 1186.070 0.013 - 7324.597 - 608.872 3 - 5809.695 2035.044 0.047 - 11571.072 - 48.318 4 - 15010.000 1601.164 0.000 - 19551.362 - 10485.308 2 1 3966.734 1186.070 0.013 608.872 7324.597 3 - 1843.000 1972.715 0.832 - 7427.880 3741.958 4 - 11050.000 1521.160 0.000 - 15358.130 - 6745.072 3 1 5809.695 2035.044 0.047 48.318 11571.072 2 1842.961 1972.715 0.832 - 3741.958 7427.880 4 - 9208.640 2246.902 0.001 - 15569.805 - 2847.476 4 1 15018.335 1601.164 0.000 10485.308 19551.362 2 11051.601 1521.160 0.000 6745.072 15358.130 3 9208.640 2246.902 0.001 2847.476 15569.805  I s q 1 2 - 0.120 0.034 0.007 - 0.217 - 0.024 3 - 0.209 0.058 0.006 - 0.375 - 0.044 4 - 0.498 0.046 0.000 - 0.629 - 0.368 2 1 0.120 0.034 0.007 0.024 0.217 3 - 0.088 0.057 0.490 - 0.249 0.072 4 - 0.377 0.044 0.000 - 0.501 - 0.254 3 1 0.209 0.058 0.006 0.044 0.375 2 0.088 0.057 0.490 - 0.072 0.249 4 - 0.289 0.065 0.000 - 0.472 - 0.107 4 1 0.498 0.046 0.000 0.368 0.629 2 0.377 0.044 0.000 0.254 0.501 3 0.289 0.065 0.000 0.107 0.472  I s d 1 2 - 0.181 0.033 0.000 - 0.274 - 0.088 3 - 0.376 0.056 0.000 - 0.536 - 0.217 4 - 0.992 0.044 0.000 - 1.118 - 0.867 2 1 0.181 0.033 0.000 0.088 0.274 3 - 0.194 0.055 0.007 - 0.350 - 0.040 4 - 0.811 0.042 0.000 - 0.930 - 0.692 3 1 0.376 0.056 0.000 0.217 0.536 2 0.194 0.055 0.007 0.040 0.350 4 - 0.616 0.062 0.000 - 0.792 - 0.440 4 1 0.992 0.044 0.000 0.867 1.118 2 0.811 0.042 0.000 0.692 0.930 3 0.616 0.062 0.000 0.440 0.792


## 5.1. One-Way Analysis of Variance (ANOVA) Method

In order to analyse the differences of the effects of the above three assessment indexes on the potential damage categories of masonry structure, a comprehensive statistic comparison method is applied in the study, which is called one-way Analysis of Variance (ANOVA) [29]. This method is a kind of univariate analysis, which majorly studies if the difference of the one of the control variables has a significant influence on the total observed variables. Based on the theory studies of Lawal and Famoye [30], the mathematical model of one-way ANOVA is explained simply as follows:(1)
Assuming the control variableA has k levels and each level has r samples, the jth sample xij with Ai control variable level is(15)xij=μi+εiji=1,2,3,…,k;j=1,2,…,r,where μi is the theory value of the Ai control variable and εij is random sampling error.(2)
Assumeu is the total theory value of control variable  A, which is expressed as follows:(16)μ=1k∑i=1kμi.(3)
So the additional effect of control variablesAi on the experimental results αi is obtained and the total value of αi is zero, given by (17)αi=μi-μ,i=1,2,…,k,(18)∑i=1kαi=0.(4)
Based on (16) and (17), xij is transferred as (19) and called mathematic model of the one-way ANOVA:(19)xij=μ+αi+εij,i=1,2,…,k;j=1,2,…,r.(5)
In the equation, the unbiased estimators of the total theory value of control variables (μ) and the additional effect of control variables Ai on the experimental results (αi) are(20)∧μ=x¯,∧αi=xi¯-x¯.If control variableA has an influence on observed variables, at least one of the values of αi is not zero. Otherwise, all of the values of αi are zero. The detailed analysis is called F-test analysis, whose completed steps are as follows:(1)
Firstly, two null hypotheses are proposed, that is,Haccept and Hobject.(2)
Then, calculate the selection test statisticF value, which is defined as(21)F=SSA/k-1SSE/m-k=MSAMSE.In this equation, m is total number of samples; SSA is the sum of squares between groups; SSE is the sum of squares in the groups; k-1 and m-k are the freedom degree of SSA and SSE; MSA is the sum of squares of between the average values groups; MSA is the sum of squares of in the average values groups. The standard value of F is from F-distribution having (k-1,m-k) of freedom degree.(3)
CompareF value and probability P value; when the observed F value is greater than 1, that means the control variables have a significant impact on the observed variables; otherwise the impact is not significant. Additionally, the larger F value is, the more significant the effect is.Based on the above, one-way ANOVA is performed to study the effects of the damage categories on the three proposed indexesIw, Isq, and Isd at different seismic intensity zone. Here, as an example, the one in seismic intensity IX zone is presented in the study. A statistical package for the social sciences statistical program (SPSS) was used and four damage groups were created in the one-way ANOVA study (group 1 collapse damage, group 2 heavy damage, group 3 moderate damage, and group 4 slight or no damage). Firstly, it was assumed a single experiment variable (main control variable) has k levels in the ANOVA, that is, the above damage degree. Based on this, the value of k equals 4. Therefore, the test level of overall significance is obtained using the following hypotheses:Hobject: α1=α2=α3=α4=0.Haccept: at least one of the values of αk is not equal to zero.Subsequently, the above hypotheses were performed and determined by comparing computedF value and the critical one in one-way ANOVA [31]. The critical values of F were obtained according to some statistic variables of specimen database. On the other hand, in some cases, the computed F values are close which means the difference is not obvious. For these reasons, in this study, the proportion of variance (η2) was also used in the comparison which was provided by McClave and Sincich [31], given by(22)η2=SSBetweenGroupsSSTotal,where SSBetweenGroups is the sum of squares of between groups and SSTotal is the sum of squares in total.


## 5.2. Multiple Comparisons: Scheffé’s Method

If the hypothesis is rejected as a result of one-way analysis of variance, it namely means that the various levels of assessment indexes are significantly different. However, this only shows that at least a difference between the two levels of assessment indexes is significant and it cannot say that there are significant differences between any two levels of the assessment indexes with the damage degree. For more than three assessment indexes, it often needs to determine where the specific significant differences are, that is, between which levels are the assessment indexes found. Therefore, if a refusal result was obtained from the analysis of one-way ANOVA, the further comparison should be finished for any difference between any two average groups from a number of average groups. This method is generally named multiple-comparison. The common method of multiple comparisons has Tukey–Kramer’s and Scheffé’s method. Tukey–Kramer method is better when only a fixed number of pairwise comparisons are to be made which will result in a more precise confidence interval. Scheffé’s method is more appropriate to the cases having a large number of comparisons, on the other hand, when there are many or all contrasts in samples. Therefore, in this study, Scheffé’s method was used in the following comparison.Firstly, the mathematical model of Scheffé’s method [32, 33] is illustrated as follows:(1)
When the evaluation result of one-way ANOVA is to reject the hypotheses, that is, the various levels of assessment indexes are significantly different, the difference between theith and jth groups can be compared via using a calculated S value which is defined as(23)S=max1≤i,j≤r⁡X¯i-X¯jMSE1/ni+1/nj.(2)
The fractile quantile ofS value (Sa) can be obtained from S-table or from the following equation:(24)Sa=r-1Fαr-1,n-r.(3)
Similar to the steps in one-way ANOVA, the steps of Scheffé’s method are as follows:
(a)
Firstly, a null hypothesis is proposed:
H0: μi=μj, i≠j, i,j=1,2,3,4,…,r.H1: H0 is not correct.(b)
CalculateSa based on (24) or from S-table and then calculate Tij which is defined as(25)Tij=Sani+njninj.(c)
CompareTij and the absolute difference value of dij between sample and average values. dij is defined as(26)dij=xi-xj.If the absolute difference value ofdij is greater than Tij, that means there is an obvious difference between the ith and jth groups; otherwise, their difference is not significant.


## 5.3. Comparison and Discussion

Using the above two comparison methods, the three indexes are compared in detail. In this study, only the results in seismic intensity IX zone are presented, as a representative result. The other comparison can be reached at the detailed research report from the authors [34].


### 5.3.1. One-Way ANOVA

Table2 shows the one-way ANOVA comparison results of the assessment results of these indexes. There are four groups which were used in the analysis, that is, collapse, heavy damage, moderate damage, and slight or no damage. Results show that the critical values of F were both taken as 2.60 for the two indexes Iw and Isd according to the values of df [31]. The results also indicate that the calculated F value in the ANOVA is 29.74, which is larger than its critical value. Therefore, the hypotheses is rejected strongly which means index Iw has a significant effect on the damage categories of masonry structures. In other words, index Iw can be used to clearly distinguish or predict the potential damage categories of the masonry buildings in seismic intensity IX zone. Comparing with the results of index Iw, however, the computed F values of indexes Isq and Isd are better and reach 40.28 and 176.008, respectively. This implies that index Isd has a more significant influence on the damage categories of surveyed masonry structures. On the other hand, the significance degrees (Sig.) of these indexes all are equal to zero, meaning that they all have significant influences on the damage degree of masonry structures. Besides, the result shows that the value of η2 for indexes Iw, Isq, and Isd are 0.39, 0.47, and 0.79, respectively. This illustrates that the combined assessment index Isd is the best index to predict the potential damage categories of the masonry buildings in earthquake-prone zone. In summary, comparing other indexes, Isd has the largest computed F value, which means the index has the most significant influence on the damage levels of these masonry buildings.Table 2
Comparison of the results of one-way ANOVA (seismic intensity IX zone).
Indexes Sum of squares df Mean square CalculatedF value CriticalF value Proportion of variance (η2) Sig. Iw Between groups 3.25 E + 09 3 1.08 E + 09 29.735 2.6 3.98 E - 01 0 Within groups 4.92 E + 09 135 3.64 E + 07 Total 8.17 E + 09 138 Isq Between groups 3.629 3 1.21 40.282 2.6 4.72 E - 01 0 Within groups 4.054 135 0.03 Total 7.683 138 Isd Between groups 14.772 3 4.924 176.008 2.6 7.96 E - 01 0 Within groups 3.777 135 0.028 Total 18.548 138


### 5.3.2. Scheffé’s Comparison

Table3 shows the results of Scheffé’s analysis comparison of the three indexes in seismic intensity IX zone. Firstly, regarding index Iw, the results present that there are significant difference between the first group (collapse) and other groups, because their probability P values are less than 0.05, a supposed significant level. Similar results are observed for the fourth group (slight or no damage). For the moderate and heavy damage group, however, they also present significant difference levels with other kinds of damage, except for each other, that is, the difference between them. In other words, they just can differentiate between the buildings with the collapse and slight/no damage groups, respectively. It indicates that Iw is not able to differentiate the buildings with moderate and heavy damage. For the assessment index Isq, similar results are confirmed to the ones of index Iw. The analysis results of the buildings with collapse and slight/no damage all show a significant difference with the corresponding three types of damage. However, comparing with the ones of index Iw, the P values of the moderate damage group were improved which has no obvious difference with the heavy damage. Regarding the combined index Isd, all damage groups present well a significant difference level with the other groups. The appearance probability of no difference between any two groups is zero which means this index effectively divides all kind of damage degrees including the moderate and heavy damage types for masonry structures in seismic intensity IX zone.Table 3
Comparison of the results of Scheffé’s analysis (seismic intensity IX zone).
Indexes ( I ) group ( J ) group Mean difference(I-J) Std. error Sig. 95% Confidence interval Lower bound Upper bound I w 1 2 - 3966.734 1186.070 0.013 - 7324.597 - 608.872 3 - 5809.695 2035.044 0.047 - 11571.072 - 48.318 4 - 15010.000 1601.164 0.000 - 19551.362 - 10485.308 2 1 3966.734 1186.070 0.013 608.872 7324.597 3 - 1843.000 1972.715 0.832 - 7427.880 3741.958 4 - 11050.000 1521.160 0.000 - 15358.130 - 6745.072 3 1 5809.695 2035.044 0.047 48.318 11571.072 2 1842.961 1972.715 0.832 - 3741.958 7427.880 4 - 9208.640 2246.902 0.001 - 15569.805 - 2847.476 4 1 15018.335 1601.164 0.000 10485.308 19551.362 2 11051.601 1521.160 0.000 6745.072 15358.130 3 9208.640 2246.902 0.001 2847.476 15569.805  I s q 1 2 - 0.120 0.034 0.007 - 0.217 - 0.024 3 - 0.209 0.058 0.006 - 0.375 - 0.044 4 - 0.498 0.046 0.000 - 0.629 - 0.368 2 1 0.120 0.034 0.007 0.024 0.217 3 - 0.088 0.057 0.490 - 0.249 0.072 4 - 0.377 0.044 0.000 - 0.501 - 0.254 3 1 0.209 0.058 0.006 0.044 0.375 2 0.088 0.057 0.490 - 0.072 0.249 4 - 0.289 0.065 0.000 - 0.472 - 0.107 4 1 0.498 0.046 0.000 0.368 0.629 2 0.377 0.044 0.000 0.254 0.501 3 0.289 0.065 0.000 0.107 0.472  I s d 1 2 - 0.181 0.033 0.000 - 0.274 - 0.088 3 - 0.376 0.056 0.000 - 0.536 - 0.217 4 - 0.992 0.044 0.000 - 1.118 - 0.867 2 1 0.181 0.033 0.000 0.088 0.274 3 - 0.194 0.055 0.007 - 0.350 - 0.040 4 - 0.811 0.042 0.000 - 0.930 - 0.692 3 1 0.376 0.056 0.000 0.217 0.536 2 0.194 0.055 0.007 0.040 0.350 4 - 0.616 0.062 0.000 - 0.792 - 0.440 4 1 0.992 0.044 0.000 0.867 1.118 2 0.811 0.042 0.000 0.692 0.930 3 0.616 0.062 0.000 0.440 0.792


## 5.3.1. One-Way ANOVA

Table2 shows the one-way ANOVA comparison results of the assessment results of these indexes. There are four groups which were used in the analysis, that is, collapse, heavy damage, moderate damage, and slight or no damage. Results show that the critical values of F were both taken as 2.60 for the two indexes Iw and Isd according to the values of df [31]. The results also indicate that the calculated F value in the ANOVA is 29.74, which is larger than its critical value. Therefore, the hypotheses is rejected strongly which means index Iw has a significant effect on the damage categories of masonry structures. In other words, index Iw can be used to clearly distinguish or predict the potential damage categories of the masonry buildings in seismic intensity IX zone. Comparing with the results of index Iw, however, the computed F values of indexes Isq and Isd are better and reach 40.28 and 176.008, respectively. This implies that index Isd has a more significant influence on the damage categories of surveyed masonry structures. On the other hand, the significance degrees (Sig.) of these indexes all are equal to zero, meaning that they all have significant influences on the damage degree of masonry structures. Besides, the result shows that the value of η2 for indexes Iw, Isq, and Isd are 0.39, 0.47, and 0.79, respectively. This illustrates that the combined assessment index Isd is the best index to predict the potential damage categories of the masonry buildings in earthquake-prone zone. In summary, comparing other indexes, Isd has the largest computed F value, which means the index has the most significant influence on the damage levels of these masonry buildings.Table 2
Comparison of the results of one-way ANOVA (seismic intensity IX zone).
Indexes Sum of squares df Mean square CalculatedF value CriticalF value Proportion of variance (η2) Sig. Iw Between groups 3.25 E + 09 3 1.08 E + 09 29.735 2.6 3.98 E - 01 0 Within groups 4.92 E + 09 135 3.64 E + 07 Total 8.17 E + 09 138 Isq Between groups 3.629 3 1.21 40.282 2.6 4.72 E - 01 0 Within groups 4.054 135 0.03 Total 7.683 138 Isd Between groups 14.772 3 4.924 176.008 2.6 7.96 E - 01 0 Within groups 3.777 135 0.028 Total 18.548 138


## 5.3.2. Scheffé’s Comparison

Table3 shows the results of Scheffé’s analysis comparison of the three indexes in seismic intensity IX zone. Firstly, regarding index Iw, the results present that there are significant difference between the first group (collapse) and other groups, because their probability P values are less than 0.05, a supposed significant level. Similar results are observed for the fourth group (slight or no damage). For the moderate and heavy damage group, however, they also present significant difference levels with other kinds of damage, except for each other, that is, the difference between them. In other words, they just can differentiate between the buildings with the collapse and slight/no damage groups, respectively. It indicates that Iw is not able to differentiate the buildings with moderate and heavy damage. For the assessment index Isq, similar results are confirmed to the ones of index Iw. The analysis results of the buildings with collapse and slight/no damage all show a significant difference with the corresponding three types of damage. However, comparing with the ones of index Iw, the P values of the moderate damage group were improved which has no obvious difference with the heavy damage. Regarding the combined index Isd, all damage groups present well a significant difference level with the other groups. The appearance probability of no difference between any two groups is zero which means this index effectively divides all kind of damage degrees including the moderate and heavy damage types for masonry structures in seismic intensity IX zone.Table 3
Comparison of the results of Scheffé’s analysis (seismic intensity IX zone).
Indexes ( I ) group ( J ) group Mean difference(I-J) Std. error Sig. 95% Confidence interval Lower bound Upper bound I w 1 2 - 3966.734 1186.070 0.013 - 7324.597 - 608.872 3 - 5809.695 2035.044 0.047 - 11571.072 - 48.318 4 - 15010.000 1601.164 0.000 - 19551.362 - 10485.308 2 1 3966.734 1186.070 0.013 608.872 7324.597 3 - 1843.000 1972.715 0.832 - 7427.880 3741.958 4 - 11050.000 1521.160 0.000 - 15358.130 - 6745.072 3 1 5809.695 2035.044 0.047 48.318 11571.072 2 1842.961 1972.715 0.832 - 3741.958 7427.880 4 - 9208.640 2246.902 0.001 - 15569.805 - 2847.476 4 1 15018.335 1601.164 0.000 10485.308 19551.362 2 11051.601 1521.160 0.000 6745.072 15358.130 3 9208.640 2246.902 0.001 2847.476 15569.805  I s q 1 2 - 0.120 0.034 0.007 - 0.217 - 0.024 3 - 0.209 0.058 0.006 - 0.375 - 0.044 4 - 0.498 0.046 0.000 - 0.629 - 0.368 2 1 0.120 0.034 0.007 0.024 0.217 3 - 0.088 0.057 0.490 - 0.249 0.072 4 - 0.377 0.044 0.000 - 0.501 - 0.254 3 1 0.209 0.058 0.006 0.044 0.375 2 0.088 0.057 0.490 - 0.072 0.249 4 - 0.289 0.065 0.000 - 0.472 - 0.107 4 1 0.498 0.046 0.000 0.368 0.629 2 0.377 0.044 0.000 0.254 0.501 3 0.289 0.065 0.000 0.107 0.472  I s d 1 2 - 0.181 0.033 0.000 - 0.274 - 0.088 3 - 0.376 0.056 0.000 - 0.536 - 0.217 4 - 0.992 0.044 0.000 - 1.118 - 0.867 2 1 0.181 0.033 0.000 0.088 0.274 3 - 0.194 0.055 0.007 - 0.350 - 0.040 4 - 0.811 0.042 0.000 - 0.930 - 0.692 3 1 0.376 0.056 0.000 0.217 0.536 2 0.194 0.055 0.007 0.040 0.350 4 - 0.616 0.062 0.000 - 0.792 - 0.440 4 1 0.992 0.044 0.000 0.867 1.118 2 0.811 0.042 0.000 0.692 0.930 3 0.616 0.062 0.000 0.440 0.792


## 6. Discussion and Recommendation

Based on the above analyses, the combined indexIsd is suggested to predict the potential damage of masonry structures when the detailed information of the used materials in the inspected structures is available. However, for the existing buildings without the detailed information of materials, indexes Iw and Isq still can be used to predict the potential damage of the structures. Using the same methods reported above, the study also analyses the buildings located at the seismic intensity VIII and X zones of the 2008 Wenchuan earthquake. Through analysing the relationship between three proposed indexes and damage degree, some relevant suggestions of masonry structures built in seismic-prone zone are provided as follows.


### 6.1. Seismic Intensity and Damage

Through studying the relative increase indexes of the proposed three assessment indexes, the total damage characteristics of masonry structures in seismic intensity VIII, IX, and X zones are summarized. The analysis method is as follows:(1)
Suppose the index values of no collapse damage in seismic intensity VIII–X zones asA1–A3 and the ones of no heavy and moderate damage as A4–A6 and A7–A9, respectively. The detailed relationships between three assessment indexes and damage types of masonry buildings are presented in Figures 5–7. The detailed values of the indexes are reported elsewhere [34].(2)
According to these values ofA1–A9 of the above three indexes, the relationships between them and the actual damage types of the masonry structures are obtained.(3)
Calculate and compare the relative increase indexesm1, m2, m3, which are defined, respectively, as follows:(27)m1=A3-A2A2-A1,m2=A6-A5A5-A4,m3=A9-A8A8-A7.Figure 5
I w and the damage type of masonry structures.Figure 6
I s q and the damage type of masonry structures.Figure 7
I s d and the damage type of masonry structures.For indexesIw and Isq, according to the fact that indexes m2 and m3 are both greater than zero, the required values of Iw and Isq of masonry structures should be increased to avoid their similar level of damage when their seismic intensity changes from VIII to IX or IX to X. Meanwhile, m3>m2, m3>1.0, and it indicates that the increase in amplitude of the required values of Iw and Isq of the masonry buildings to avoid their moderate damage is higher than the one to avoid the heavy damage. For the combined index Isd, the values of m1–m3 are 1.5, 2.5, and 3.5, respectively. It means the required levels of Isd of masonry structures to avoid the corresponding damage type increase as the seismic intensity. Meanwhile, the case m3>m2>m1>1.0 indicates that the increase in amplitudes of the required Isd of masonry structures to avoid their moderate damage (heavy damage) is higher than the one of the buildings to avoid their heavy damage (collapse damage). On the other hand, the above studies on the strength and combined indexes also indicate that to reduce the number of building story of masonry structures is a good way to increase the seismic resistance capacity of total masonry building, and to use more confinement element such as tie columns is a good way to improve the capacity and deformation of total structures as well. However, the implementation of the two ways is limited for the wall density which is majorly dependent on the plan distribution and use function of masonry structures; in addition, the number of building storeys is preset by users, and it cannot be changed by designers or engineers. It is difficult to avoid the appearance of damage in the masonry structures in high seismic intensity zone, especially the slight or moderate damage, because the increase of confinement elements and the increased requirement of wall density both are limited.


### 6.2. Application of Tie Columns in Masonry Structures

Comparing with the result in Figure5, the result plotted in Figure 8 clearly shows the positive influence of tie columns in masonry structures. These columns can effectively improve the total properties of masonry structures and enhance the confinement of masonry wall, which both are good for the resistant of masonry walls during an earthquake, in particular for the masonry buildings with potential moderate and heavy damage.Figure 8
I s q and the damage type of masonry structures with/without tie columns.


## 6.1. Seismic Intensity and Damage

Through studying the relative increase indexes of the proposed three assessment indexes, the total damage characteristics of masonry structures in seismic intensity VIII, IX, and X zones are summarized. The analysis method is as follows:(1)
Suppose the index values of no collapse damage in seismic intensity VIII–X zones asA1–A3 and the ones of no heavy and moderate damage as A4–A6 and A7–A9, respectively. The detailed relationships between three assessment indexes and damage types of masonry buildings are presented in Figures 5–7. The detailed values of the indexes are reported elsewhere [34].(2)
According to these values ofA1–A9 of the above three indexes, the relationships between them and the actual damage types of the masonry structures are obtained.(3)
Calculate and compare the relative increase indexesm1, m2, m3, which are defined, respectively, as follows:(27)m1=A3-A2A2-A1,m2=A6-A5A5-A4,m3=A9-A8A8-A7.Figure 5
I w and the damage type of masonry structures.Figure 6
I s q and the damage type of masonry structures.Figure 7
I s d and the damage type of masonry structures.For indexesIw and Isq, according to the fact that indexes m2 and m3 are both greater than zero, the required values of Iw and Isq of masonry structures should be increased to avoid their similar level of damage when their seismic intensity changes from VIII to IX or IX to X. Meanwhile, m3>m2, m3>1.0, and it indicates that the increase in amplitude of the required values of Iw and Isq of the masonry buildings to avoid their moderate damage is higher than the one to avoid the heavy damage. For the combined index Isd, the values of m1–m3 are 1.5, 2.5, and 3.5, respectively. It means the required levels of Isd of masonry structures to avoid the corresponding damage type increase as the seismic intensity. Meanwhile, the case m3>m2>m1>1.0 indicates that the increase in amplitudes of the required Isd of masonry structures to avoid their moderate damage (heavy damage) is higher than the one of the buildings to avoid their heavy damage (collapse damage). On the other hand, the above studies on the strength and combined indexes also indicate that to reduce the number of building story of masonry structures is a good way to increase the seismic resistance capacity of total masonry building, and to use more confinement element such as tie columns is a good way to improve the capacity and deformation of total structures as well. However, the implementation of the two ways is limited for the wall density which is majorly dependent on the plan distribution and use function of masonry structures; in addition, the number of building storeys is preset by users, and it cannot be changed by designers or engineers. It is difficult to avoid the appearance of damage in the masonry structures in high seismic intensity zone, especially the slight or moderate damage, because the increase of confinement elements and the increased requirement of wall density both are limited.


## 6.2. Application of Tie Columns in Masonry Structures

Comparing with the result in Figure5, the result plotted in Figure 8 clearly shows the positive influence of tie columns in masonry structures. These columns can effectively improve the total properties of masonry structures and enhance the confinement of masonry wall, which both are good for the resistant of masonry walls during an earthquake, in particular for the masonry buildings with potential moderate and heavy damage.Figure 8
I s q and the damage type of masonry structures with/without tie columns.


## 7. Conclusions

This paper studied the relationship between three potential assessment indexes with the actual damage degrees of 130 surveyed masonry buildings in the 2008 Wenchuan earthquake and presented the application of variance analyses comparison in the determination of the optimum seismic damage assessment method of masonry buildings using three simplified indexes. For this end, one-way analysis of variance and Scheffé’s method were used to compare comprehensively the feasibility of the three indexes. The main conclusion can be drawn as follows:(1)
From the viewpoint of mathematical statistics, the three assessment indexes proposed in the study all present discrimination in certain degree for the masonry structures having the same damage level. This means the indexes can be used to predict the potential damage types of masonry structures, especially for combined index considering strength of used materials and enhancement influence of confinement elements of masonry structures.(2)
TheF value of the strength index Isq is higher comparing with that of wall density index Iw. This means the index presents a more obvious difference between four types of damage; however, the difference between moderate and heavy damage is small.(3)
After combining with the enhancement factor of confinement elements in masonry structures, indexIsd has the highest F  value. Comparing the other two indexes, when evaluating the potential damage type of masonry structures, this index shows a more obvious distinction degree and difference between each damage level, including all four damage types.(4)
Based on the comparison and analyses, the study indicated that to reduce the number of building storeys of masonry structures is a good way to increase seismic resistance capacity of total building and to use more confinement elements such as more tie columns/beams is a good way to improve the seismic behaviour of total structures. It is difficult to avoid the appearance of damage in the masonry structures in high seismic intensity zone, especially the slight or moderate damage because the increase of confinement elements and the increased requirement of wall density both are limited.(5)
Univariate comparison analysis can be used to determine simply and fast the superiority of the three proposed assessment indexes in distinguishing capacity of the damage degree of masonry buildings such as the one-way analysis of variance. Multiple-comparison analysis is more complex but can provide a more direct and detailed comparison of the difference between any two kinds of damage levels.(6)
To resist of the damage of masonry walls during an earthquake, in particular for the buildings with potential moderate and heavy damage, the addition of tie columns is good way for it can effectively improve the total seismic properties of masonry structures and enhance the confinement of masonry wall.(7)
Since the number of URM buildings is small in the inspected database of the paper, it should be noted that the results regarding URM structures may be verified further, meaning more field investigations are expected.



---

*Source: 3741941-2017-02-20.xml*