Research Article  Open Access
Rushan Liu, Deyuan Tian, "Study on Analogy Calculation Method for Seismic Vulnerability of EarthWood Structure Houses", Mathematical Problems in Engineering, vol. 2020, Article ID 3525896, 10 pages, 2020. https://doi.org/10.1155/2020/3525896
Study on Analogy Calculation Method for Seismic Vulnerability of EarthWood Structure Houses
Abstract
Earthwood structure houses often cause large casualties and economic losses in historical earthquakes. Therefore, it is estimated that the seismic vulnerability of civil structures in areas that have not experienced earthquakes is of great significance for earthquake prevention and disaster reduction. In this paper, an analogy calculation method was proposed for calculating the seismic vulnerability of earthwood structure houses in unknown regions from the seismic vulnerability of earthwood structure houses in known regions. Firstly, the main factors affecting the seismic capacity of earthwood structure houses were determined, and the weights of influence of influencing factors on the overall seismic capacities of buildings under different seismic intensities were determined by using the fuzzy analytical hierarchy approach; secondly, the relative seismic capacities of the main influencing factors are analogically scored by considering the differences in influencing factors in different regions. Finally, based on the vulnerability matrix of earthwood structure houses in existing regions and comprehensive evaluation of relative seismic capacities of earthwood structures, the vulnerability matrix of unknown earthwood structures was calculated by analogy. The results of trial calculation showed that this method has higher reliability and effectiveness. At the end of the article, the vulnerability analysis of earthwood structure houses in Shigatse, Tibet, was carried out using this method, and the vulnerability curve was initially given.
1. Introduction
The seismic vulnerability of buildings is to estimate the degree of damage that may occur when buildings are affected by a certain intensity earthquake [1]. The seismic disaster of a region or city can be effectively reduced through correctly and reasonably analyzing the seismic vulnerability of all kinds of building structures and taking measures such as seismic reinforcement and reconstruction to improve the seismic performance of building structures that do not meet the requirements of seismic fortification [2].
Although with the continuous development of social economy and the gradual reconstruction of houses in rural areas in China, the number of earthwood structure houses is decreasing year by year in rural areas in China, and there are still more earthwood structure houses in rural areas in economically backward regions and ethnic minority autonomous regions in China. Due to the fact that seismic fortification is not considered in these earthwood structure houses and the building age is old, it is extremely easy for them to be damaged under the action of earthquakes, and their seismic resistances are poor. Therefore, it is still very important to study the seismic vulnerability of earthwood structure houses to effectively carry out the work of earthquake prevention and disaster reduction [3].
Western China is an earthquakeprone area in China; multiple earthquakes of magnitude 7 and above have occurred in history, including the Wenchuan earthquake (, May 12, 2008) and the Yushu earthquake (, April 14, 2010). There are a large number of earthwood structure houses in rural areas in western China, which have been severely damaged by earthquakes in the past [4]. The vulnerability matrix of earthwood structure houses has been obtained from the actual earthquake damage investigation in some western regions where the earthquake occurred. However, due to the lack of experiment and numerical simulation study on earthwood structures, the vulnerability matrix of earthwood structure houses cannot be established, and accurate vulnerability evaluation of earthwood structures cannot be carried out in the regions where no earthquake has occurred. Therefore, this paper proposed a calculation method for seismic vulnerability of earthwood structure houses, that is, based on the vulnerability matrix of earthwood structures obtained from seismic damage statistics in some regions of western China, considering the effects of the differences in structural components of earthwood structure houses in different regions on seismic vulnerability, through comprehensive evaluation of seismic capacities of earthwood structures in different regions, inverse distance weighted interpolation and beta distribution were used to calculate and fit the seismic vulnerability matrix of earthwood structure houses in the region where no earthquake has occurred.
After field investigations, it was found that earthwood structure houses were widely distributed in rural areas of Shigatse, Tibet. According to the structural characteristics of earthwood structure houses in Shigatse, Tibet, this paper uses the method to analyze the vulnerability of earthwood structure houses in this area.
2. General Thinking of the Matrix Simulation Method
It was assumed that the seismic vulnerability matrix of earthwood structure houses in regions (benchmark regions) was known; general thinking of deriving the vulnerability matrix of earthwood structure houses in unknown regions (target regions) from the vulnerability matrix of earthwood structure houses in regions was as follows: (1) the main factors affecting the seismic capacity of earthwood structures were determined, and the fuzzy complementary judgment matrix was established by using the fuzzy analytical hierarchy approach; the influence weights of the main influencing factors on the seismic capacities of the earthwood structure houses under different seismic intensities were calculated; (2) the seismic capacities of the main influencing factors were evaluated, and the seismic capacity scores of the earthwood structure houses in different regions were calculated combining with the influence weights of the main influencing factors on the seismic capacities of the earthwood structure houses under different seismic intensities; (3) the expected values and standard deviations of the seismic damage indexes of the target region were calculated by referring to the expected values, standard deviations, and seismic capacity scores of the seismic damage indexes under different seismic intensities on the basis of the vulnerability matrix in the benchmark region through inverse distance weighted interpolation, the expected values and standard deviations of the seismic damage indexes of the target region were taken as parameters, and the beta distribution function was used to fit the vulnerability matrix of the earthwood structure houses; and (4) the rationality and accuracy of this method were verified by example calculation and error analysis.
3. Selection and Weight Definition of the Main Influencing Factors
3.1. Selection of the Main Influencing Factors
As no seismic fortification measures are taken, earthwood structure houses will suffer obvious seismic damages such as cracks in walls and destroyed internal supporting structures when the seismic intensity is above 6° [5]. For the same type of structure, there are many factors that affect its seismic performance. For example, based on the comprehensive consideration of structural characteristics, seismic damage analysis results, and structural mechanical model, the influencing factors of seismic damage of multistory masonry structures can be summarized as 14 articles in the following: fortification standard, construction age, bearing wall thickness, mortar strength grade, number of floors, usage, roof type, roof weight, floor type, masonry method, horizontal and vertical regularization, housing status, building site, and site soil type [6]. In view of the above influencing factors, based on the study and analysis of the actual seismic damage of historical earthwood structure houses [7] and referring to the existing study results, as shown in Table 1, the number of floors, bearing walls, and internal supporting structures were determined as the main factors affecting the seismic capacity of earthwood structure houses in this paper.

3.2. Determination of Weights of Influencing Factors under Different Seismic Intensities
Through the study of the previous actual seismic damage data of earthwood structures, it was found that the influences of influencing factors on the seismic capacities were different under different seismic intensities. For example, when the seismic intensity was 6°, the seismic damage of the earthwood structure houses was mainly the damage of the wall, and the damage grade was mainly affected by the condition of the wall; however, when the seismic intensity was 8° or above, whether overall earthwood structure houses collapsed was mainly determined by the supporting strength of the internal supporting structures, and the damage grade was mainly affected by the internal supporting structures. In order to determine the influences of the main factors on the seismic capacities under different intensities, the fuzzy analytical hierarchy approach was used in this paper to construct the fuzzy judgment matrix and determine the influence weights of the main factors on the seismic capacities under different seismic intensities.
Fuzzy analytical hierarchy approach [8, 9] is used to combine the ideas and methods of fuzzy mathematics with the analytic hierarchy process to obtain the fuzzy judgment matrix according to the importance degree of one factor to another when different factors are compared and judged. In order to quantitatively describe the importance degree of the pairwise comparison between two factors, the scale method in Table 2 is usually used to quantify it [10].

Thereby, the judgment matrix of the evaluation index is constructed as :
The fuzzy judgment matrix of the main influencing factors under different seismic intensities is shown in Tables 3–6, among which the selection of the comparison coefficient was mainly referred to the literature studies of seismic engineering experts such as Yin et al. [5, 11–13]; meanwhile, this was referred to a large number of actual seismic damage examples and the judgment of experts in the same field by means of questionnaire scoring.




According to the obtained fuzzy judgment matrix, the general formula to solve the weight of the fuzzy judgment matrix proposed by Xu [8] was used to calculate the influence weights of the main influencing factors on the seismic capacities of earthwood structure houses under different intensities. The specific solution formula is as follows:
In formula (2), is the influence weight of the th main influencing factor.
The calculation results are shown in Table 7.

Whether the weight value calculated by the above formula is reasonable still needs to be checked for consistency. In this paper, according to the compatibility of the fuzzy judgment matrix [14], the consistency of weight was checked by the judgment matrix and its characteristic matrix. Characteristic matrix B is
Then, the weight compatibility index is
When the compatibility index was less than or equal to the attitude of the decision maker, it was considered that the judgment matrix satisfied the consistency. The smaller was, the higher the requirement of the decision maker for the consistency of the fuzzy judgment matrix was; generally, was taken.
In this paper, formulas (7), (8), and (9) were used to check whether the weight values in Table 7 were reasonable. The calculation results are shown in Table 8.

All compatibility indexes of the weights of the main influencing factors under different intensities were less than 0.1, so it was considered that the fuzzy judgment matrices of the influencing factors under different intensities were consistent, and the distribution of the weights was reasonable.
4. Comprehensive Evaluation of Seismic Capacities of EarthWood Structures
Because there are some differences in the specific conditions of the main factors affecting the seismic capacities in different regions, for example, the thickness of the walls is different among some regions, the seismic capacities of earthwood structure houses in different regions are different. According to the actual situation of the main influencing factors, the relative seismic capacities of the main influencing factors are determined by using the expert scoring method. Experts participating in scoring mainly come from the fields of structural engineering, earthquake engineering, and disaster prevention and mitigation engineering and have extensive experience in earthquake site inspections. In this paper, the relative seismic capacities of the main influencing factors were divided into five grades, and the better the relative seismic capacity was, the higher the score given by relevant experts was. For example, if the damage grade of the earthwood structure with thick walls is lower than that of the earthwood structure with thin walls during the earthquake [15], the expert score is higher; if the damage grade of the doublefloor earthwood structure is higher than that of the singlefloor earthwood structure during the earthquake, the expert score is lower. The scoring standard for the relative seismic capacity of each influencing factor is shown in Table 9.

According to the housing census data, we could understand the specific situations such as the number of floors of earthwood structures, bearing walls, and internal supporting structures in a certain region. Based on the specific conditions of the main influencing factors of earthwood structures, relevant experts were consulted for evaluating the relative seismic capacities of the main influencing factors. According to the influence weights of the main influencing factors on the seismic capacities of the building under different intensities, the comprehensive scores of the relative seismic capacities of the earthwood structures in a certain region under different intensities were obtained as follows:
In the formula, is the comprehensive score of the relative seismic capacity under seismic intensity ; is the score of the relative seismic capacity of the \th main influencing factor; and is the comprehensive score of the relative seismic capacity of the th main influencing factor under seismic intensity .
5. Seismic Damage Matrix Simulation
Seismic engineering experts such as Hu [16] put forward the concept of seismic damage indexes. The values between 0 and 1 are used to express the degree of seismic damage of the structure. Table 10 shows the corresponding seismic damage index range of five damage grades [1].

The expected value of the seismic damage index is the average value of the seismic damage indexes of buildings in a region. The dispersion of the damage grade of buildings in a region can be expressed by the standard deviation of the seismic damage index. Therefore, the expected value and standard deviation of the seismic damage index of a certain type of structure can be obtained according to its vulnerability matrix.
Based on the difference in the comprehensive score of the relative seismic capacity between earthwood structure houses in the benchmark region and the target region, the expected values and standard deviations of the seismic damage indexes in the benchmark region under different intensities were used to calculate inverse distance weighted (IDW) interpolation [17], and thus, the expected values and standard deviations of the seismic damage index in the target region were obtained. The weighted distance isand thus,
In the formula, is the weighted distance between the relative seismic capacities in the target region and the benchmark region under seismic intensity ; is the comprehensive score of the relative seismic capacity of the target region under seismic intensity ; is the comprehensive score of the relative seismic capacity of the th benchmark region under seismic intensity ; is the expected value of the seismic damage index of the vulnerability matrix of the target region under seismic intensity ; is the expected value of the seismic damage index for the vulnerability matrix in the th benchmark region under seismic intensity ; is the standard deviation of the seismic damage index for the vulnerability matrix in the target region under seismic intensity ; and is the standard deviation of the seismic damage index for the vulnerability matrix in the th benchmark region under seismic intensity ; Liu et al. [18] proposed a method for the use of the beta distribution function to fit the vulnerability matrix by taking the seismic damage index as a continuous variable and taking the expected value and the standard deviation of the seismic damage index as parameters. In this paper, this method was used to fit the vulnerability matrix of earthwood structure houses in the benchmark region.
According to the relationship between the expected value and the standard deviation of parameters and of the seismic damage index, the values of parameters and of the beta distribution function under intensity are calculated as follows:
Thereby, the beta probability density distribution of the seismic damage index under seismic intensity could be obtained as follows:
According to formula (13), the damage probabilities of the earthwood structure houses in the target region under different damage degrades can be fitted, and finally, the vulnerability matrix of the earthwood structure houses in the target region can be obtained.
In the formula, is the continuous variable of the seismic damage index; is the probability of grade damage under seismic intensity ; is the lower limiting value of the range of the seismic damage index for grade i damage; and is the upper limit value of the range of the seismic damage index for grade damage.
6. Analysis of Examples
The vulnerability matrix of earthwood structure houses in Gansu province was fitted using the calculation method proposed in this paper.
In this paper, Sichuan province, Xinjiang Uygur Autonomous Region, and Yunnan province in western China are selected as the benchmark regions, and the main factors affecting the seismic capacity of earthwood structure houses in the benchmark regions and Gansu region are shown in Table 11.

According to the specific conditions of the main influencing factors in each region in Table 11, the expert scores of the main influencing factors in each region are shown in Table 12.

According to formula (6), the comprehensive scores of relative seismic capacities of earthwood structure houses in Gansu province and the benchmark regions under different seismic intensities are calculated in Table 13.

Based on the scores and the expected values and standard deviations of the seismic damage indexes under different intensities in each benchmark region in Table 14, the expected values and standard deviations of seismic damage indexes under different intensities in the Gansu region were calculated according to formulas (8) and (9), and the calculation results are shown in Table 15.


Formulas (10) and (11) were used to calculate the beta distribution parameters under different seismic intensities, and the beta distribution function curve of continuous seismic damage index under different intensities is shown in Figure 1.
According to the beta probability density distribution of the obtained seismic damage indexes, the vulnerability matrix was fitted by using formula (13). The fitting results and the vulnerability matrix obtained through the actual seismic damage investigation are shown in Tables 16 and 17.


In order to verify the reliability and effectiveness of this method, the vulnerability results of earthwood structure houses estimated by this method were compared with the actual vulnerability results of seismic damage investigation. As shown in Table 18, the error between the expected value of the seismic damage index obtained by simulation in Table 15 and the expected value of the actual seismic damage index was relatively small, which basically met the requirements of vulnerability evaluation. The error between the simulated vulnerability matrix and the actual seismic damage vulnerability matrix basically met the expected results. Therefore, this method is feasible.

7. Vulnerability Analysis of Civil Structure Houses in Shigatse, Tibet
Through field research, it was found that the earthwood structure houses in Shigatse, Tibet, were mainly twostory traditional Tibetanstyle buildings with thick walls and good internal support structure, as shown in Figure 2.
(a)
(b)
Based on the actual situation of earthwood structure houses in Shigatse, Tibet, the method of this paper is used to analyze the vulnerability of earthwood structure houses in this area. Sichuan province, Xinjiang Uygur Autonomous Region, and Yunnan province in western China are selected as the benchmark regions, and the main factors affecting the seismic capacity of earthwood structure houses in the benchmark regions and Shigatse, Tibet, are shown in Tables 11 and 19.

According to the specific conditions of the main influencing factors in Shigatse region in Table 19, the expert scores of the main influencing factors in Shigatse region are shown in Table 20.

Based on the scores and the expected values and standard deviations of the seismic damage indexes under different intensities in each benchmark region in Table 14, the expected values and standard deviations of seismic damage indexes under different intensities in Shigatse region are calculated according to formulas (6), (8), and (9), and the calculation results are shown in Table 21.

Formulas (10) and (11) were used to calculate the beta distribution parameters under different seismic intensities and obtain the beta distribution function of continuous seismic damage index under different intensities. Finally, formula (13) is used to fit the vulnerability matrix in Table 22.

According to the vulnerability matrix of earthwood structure houses in Shigatse region calculated using the method in this paper, the vulnerability curve of earthwood structure houses in this region can be obtained, as shown in Figure 3.
The calculation results show that the moderate and above damage degree of earthwood structure houses in Xigaze area occupied the main part when the earthquake intensity is 8°; when the earthquake intensity is 10°, the earthwood structure houses in the area are basically destroyed. Comparing the expected values of the earthquake damage index of earthwood structures in Shigatse and other benchmark regions, because the load bearing wall of the earthwood structure is thicker and good overall condition in Xigaze, the expected value of the earthquake damage index of earthwood structure houses in Xigaze is relatively low in lowintensity areas; in highintensity areas, due to the poor internal support structure compared with Yunnan, the expected value of the earthquake damage index for structural houses is higher than Yunnan.
It can be seen that the vulnerability analysis of the earthwood structure houses in Shigatse, Tibet, using this method meets the expected results.
8. Conclusion
The vulnerability matrix of earthwood structure houses in Gansu province was fitted using the calculation method proposed in this paper. In this paper, according to the needs of vulnerability estimation of earthwood structure houses, a quick calculation method for fitting the vulnerability matrix of earthwood structure houses is proposed. Based on the vulnerability matrix of the earthwood structures in the benchmark region, considering the influence weights of the main influencing factors on the seismic capacity of the building under different seismic intensities and the difference in the actual situation of the influencing factors in different regions, the comprehensive scores of relative seismic capacities of the earthwood structure houses in the benchmark region and the target region are calculated, respectively, and then, the expected values and standard deviations of the seismic damage indexes in the target are obtained by using the inverse distance weighted interpolation method. Finally, the beta distribution function is used to fit the seismic damage matrix in the simulated region. The earthwood structures in Gansu province are checked by using the method in this paper, and the error between the calculated results and the actual seismic damage results is within a reasonable range.
The method proposed in this paper can quickly evaluate the vulnerability of earthwood structure houses in a certain region, which considers the influence of different influencing factors on the seismic capacities of the structures under different seismic intensities compared with the previous methods. The accuracy of this method is affected by the selection of the benchmark region and the reasonableness of the relative seismic capacity scores evaluated by relevant experts on the main influencing factors. Generally, the region with earthwood structures similar to those of the target region should be selected as the benchmark region, and the average value of the scores of multiple relevant experts should be used for evaluation. In this paper, only the vulnerability of earthwood structures is evaluated, and how this analogy method can also be applied to houses of other types of structures needs further study.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest regarding the publication of this paper.
Acknowledgments
This study was supported by Special Fund for Basic Scientific Research Expenses of Institute of Engineering Mechanics (2018A02), National Key R&D Program of China (2018YFC1504602), and Business Expenses of China Earthquake Administration (19024205).
References
 Z. Yin, in Method of Prediction for Seismic Disaster and Loss, Seismological Press, Beijing, China, 1996.
 S. Hu, B. Sun, and D. Wang, “Approach to making empirical earthquake damage matrix,” Earthquake Engineering And Engineering Vibration, vol. 27, no. 06, pp. 46–50, 2007, in Chinese. View at: Google Scholar
 B. Chen, “Analysis of earthquake damage of rural buildings in western China and study on seismic resistant measures,” Construction Quality, vol. 32, no. 06, pp. 23–27, 2014, in Chinese. View at: Google Scholar
 Y. Yuan, “Loss assess ment of wenchuan earthquake,” Earthquake Engineering and Engineering Vibration, vol. 28, no. 05, pp. 10–19, 2008, in Chinese. View at: Google Scholar
 Z. Yin and S. Yang, Seismic Loss Analysis and Fortification Criterion, Seismological Press, Beijing, China, 2004, in Chinese.
 G. Zhang and B. Sun, “A method for earthquake damage prediction of building groups based on multiple factors,” World Earthquake Engineering, vol. 26, no. 01, pp. 26–30, 2010, in Chinese. View at: Google Scholar
 B. Sun and G. Zhang, “Study on vulnerability matrices of masonry buildings of mainland China,” Earthquake Engineering and Engineering Vibration, vol. 17, no. 2, pp. 251–259, 2018. View at: Publisher Site  Google Scholar
 Z. Xu, “Algorithm for priority of fuzzy complementary judgement matrix,” Journal of Systems Engineering, vol. 16, no. 04, pp. 311–314, 2001, in Chinese. View at: Google Scholar
 B. Sun and S. Hu, “Seismic damage prediction for a single building based on a fuzzy analytical hierarchy approach,” Engineering Mechanics, vol. 35, no. 12, pp. 185–193, 2018, in Chinese. View at: Google Scholar
 D. Ji, B. Song, and T. Yu, “The method of decisionmaking based on FAHP and its application,” Fire Control and Command Control, vol. 32, no. 11, pp. 38–41, 2007, in Chinese. View at: Google Scholar
 B. Sun and G. Zhang, “A method for earthquake damage prediction of building group considering seismic fortification level,” China Civil Engineering Journal, vol. 46, no. 5, pp. 117–121, 2013, in Chinese. View at: Google Scholar
 B. Sun and G. Zhang, “The Wenchuan earthquake creation of a rich database of building performance,” Science China Technological Sciences, vol. 53, no. 10, pp. 2668–2680, 2010. View at: Publisher Site  Google Scholar
 X. Yao, Analysis of Tianjin Rural Residential Vulnerability and Seismic Capacity distribution, Institute of Engineering Mechanics, China Earthquake Administration, Beijing, China, 2016, in Chinese.
 H. Chen and J. Zhao, “Research on compatibility of fuzzy judgement matrices,” Operations Research and Management Science, vol. 13, no. 01, pp. 44–47, 2004, in Chinese. View at: Google Scholar
 S. Luo, Antiseismic and Practice of Rural Housing in Western China, Xi’an University of Architecture and Technology, Xi'an, China, 2011, in Chinese.
 Y. Hu, Earthquake Engineering, Seismological Press, Beijing, China, 1988, in Chinese.
 Z. Tang, 3D Data Field Visualazation, Tsinghua University Press, Beijing, China, 1999, in Chinese.
 R. Liu, S. Hu, and Y. Wu, “Research on structure vulnerability expression way based on seismic motion parameter,” Earthquake Engineering and Engineering Vibration, vol. 29, no. 6, pp. 102–107, 2009, in Chinese. View at: Google Scholar
Copyright
Copyright © 2020 Rushan Liu and Deyuan Tian. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.