Abstract

Based on the concept of environmental protection of solid waste utilization, material testing is conducted to achieve native improvement using coal gangue-based limestone-calcined clay cement (LC3). Finite element (FE) models of rural raw-soil architecture with a colored-steel roof (RACSR) were established. The effect of modified soil type and seismic character on the vulnerability of single-story raw-soil structures was investigated using probabilistic seismic demand (PSD) analysis. The seismic response characteristics of 80 representative sequences were comparatively investigated when subjected to northwest clay (raw soil) of China, fiber and stone-improved clay (modified soil), and coal gangue-based limestone-calcined clay cement (LC3 soil). The maximum interstory drift angle (ISDAmax) was lower in the LC3 soil model and the modified soil model compared to the raw-soil model. The use of LC3 soil improves structural resistance and reduces the damage probability of a structure, and the influence of different ultimate failure states on the vulnerability of the raw-soil structure was studied.

1. Introduction

Earthquakes are natural catastrophes that severely impact society owing to the associated loss of life and damage to building facilities. To date, most studies have focused on performance-based seismic design and risk assessment, and research on seismic vulnerability has primarily considered steel, concrete, and masonry structures. It is necessary to study the seismic vulnerability of rural buildings as part of earthquake disaster loss assessment and structural reinforcement. Raw-soil structures are local material, facilitate simple construction, have low cost, are warm in winter and cool in summer, and promote energy saving and environmental protection. They are still widely used in rural areas of China, especially in economically underdeveloped regions [1]. However, owing to the poor ductility, shear resistance, and bending strength of rural buildings, if coupled with injudicious site selection, poor foundation treatment, poor quality masonry, and planning and layout defects, the degree of seismic damage is much higher compared to that of other structures [2].

To reduce the casualties and damage associated with seismic events, numerous studies on the antiseismic performance of raw-soil structures have been conducted. The emphasis has been on structural construction measures and the improvements of raw-soil material. Miccoli et al. [3] conducted tests of rammed Earth walls strengthened with polyester fabric strips and discussed the effect of cyclic loading. Ma et al. [4] investigated adobe masonry strengthened with ferrocement mortar based on theoretical analysis, experimental research, and finite element simulation and found that the failure mode is a shear failure and shear-compression failure. Klingner et al. [5] conducted a static pushover test of an adobe house strengthened with a wooden frame. The results showed that the integrity and bearing capacity can be effectively improved. Based on the investigation and analysis of the earthquake damage mechanism of raw-soil buildings, Pavan et al. [6] analyzed the interfacial behavior of cement-stabilized rammed Earth. Zhou et al. [7] conducted a shaking table test and analyzed the scale model of two types of adobe house reinforcement modes. The results showed that the proposed plaster-adobe-wall structure and the wooden-frame-adobe-wall structure could resist earthquakes of magnitude 8.0. Bu et al. [8] proposed a mixed raw-soil structure and conducted a static test and finite element simulation on the wall opening. The results showed that the raw-soil-mixed wall can effectively address the brittle failure of traditional raw-soil walls during earthquakes, and the opening rate of the wall had a significant impact on the seismic performance. Bui et al. [9] focused on the antiseismic design parameters of raw-soil structures and proposed a formula for seismic calculations. According to the experimental and numerical assessment, Li et al. [10] investigated effective antiseismic measures that can be applied to adobe walls. Chen et al. [11] analyzed raw-soil building performance and damage associated with the Nepal earthquake. Illampas et al. [12] analyzed the influence of horizontal loading on the performance of adobe masonry building by laboratory testing and FE simulation. Ghasemalizadeh and Toufigh [13] conducted an experimental study on the durability of rammed building walls and explored the modification of ram soil materials. Tripura and Singh [14] discussed the mechanical behavior of columns reinforced with steel and bamboo. Tripura et al. [15] further discussed the flexural strength and failure trend of bamboo and coir reinforced cement-stabilized rammed Earth wallets. An efficacious design and construction system for raw-soil structures is not currently available, and many self-built houses lack corresponding structural measures, which will significantly affect the seismic performance of raw-soil structures. Research on raw-soil materials mainly focuses on standardized production processes, structural construction measures, and methods to improve raw-soil materials. Structural construction measures are implemented to improve the overall stability and rigidity of the raw-soil structure. However, studies on the improvement of raw-soil materials have mainly focused on local materials and solid waste utilization, such as the use of natural fibers and life/industrial waste to improve the raw soil. Therefore, it is necessary to further improve the theoretical basis of seismic design and risk assessment methods for raw-soil structures.

Seismic vulnerability refers to the possible different failure states of engineering structures under different earthquake intensities, which are used for pre- and postearthquake loss assessments. Seismic vulnerability analysis was first used to evaluate the seismic performance of important structures such as nuclear power plants. With the increase of property losses in casualties as earthquake damage to buildings, seismic vulnerability research has been gradually applied to the evaluation of the seismic performance of common engineering structures. From the perspective of area-wide estimation, An et al. [16] redefined the seismic vulnerability analysis based on the decision tree method. Zhao et al. [17] established 18 structural models of village and town buildings according to the three categories of rural, town, and county and obtained the corresponding vulnerability curves to evaluate the rationality of the three levels of fortification standards. Li et al. [18] discussed the seismic vulnerability and spatial distribution of native soil buildings in western China and proposed a multidimensional information acquisition survey method. He et al. [19] conducted a statistical analysis of the damage and economic losses caused by 80 earthquake disasters in Yunnan Province from 1993 to 2015 and obtained vulnerability curves for different seismic intensities and building structures. Considering the damage samples of buildings in the Wenchuan earthquake area, Sun and Zhang [20] analyzed the seismic damage characteristics of typical structures and performed seismic vulnerability analysis and seismic capacity studies. The results showed that the masonry structure was the most seriously damaged among the buildings in the Wenchuan earthquake-stricken areas. Based on the incremental dynamic analysis (IDA) collapse vulnerability method, Ren et al. [21] compared and discussed the influence of the magnitude, epicentral distance, site conditions, and other ground motion parameters associated with earthquakes on the structural vulnerability curve of three different structures. Lourenço and Roque [22] proposed a seismic vulnerability index for traditional masonry structures that included the statistical analysis of the data for 58 historical churches in Portugal. Lagomarsino and Giovinazzi [23] proposed two models and applied them to the analysis of different types of buildings to evaluate the vulnerability and to examine damage assessment, which can be applied to European towns and regions. Based on the equivalent frame method and SAP2000 software, Pasticier et al. [24] obtained seismic vulnerability curves for typical two-story masonry structures by linearly fitting the data generated from incremental dynamic analysis. Karbassi and Nollet [25] selected 14 synthetic and real ground motion records and utilized them in a 3D model of an unreinforced masonry structure based on the applied element method to analyze seismic vulnerability. The estimation model of horizontal aftershock ground motion in the active crustal region was developed by Byungmin and Moochul [26], which is a useful technique for the vulnerability analysis of raw-soil buildings. Giovanni and Claudio [27] investigated several methods and concluded that the influence of masonry structures under earthquake sequences is important as an engineering reference. Zhang et al. [28] established a transient seismic structure-water-sediment-rock (SWSR) interaction model and discussed its application. Wang et al. [29] obtained the dynamic resistance factor of the soil layer and proposed a substructure method for seismic responses that consider nonlinear pile-soil dynamic interactions. Pang et al. [30] proposed a uniform design- (UD-) based Gaussian process regression (GPR) method for rapid damage assessment and vulnerability estimates of bridges. Pang and Wang [31] presented an efficient method for converting vulnerability curves from the cloud to IDA and MSA and from IDA to MSA toward an efficient and high-fidelity resilience assessment.

However, there is a dearth of research on the vulnerability of raw-soil structures, which is critical for improving seismic design and seismic risk assessment methods of raw-soil structures. In particular, raw-soil structures in rural areas are typically not optimally designed, and many self-built houses lack the necessary confining members; therefore, these raw-soil structures often exhibit inferior seismic performance. To address this problem, the seismic response characteristics and seismic vulnerability of raw-soil structures are investigated in rural areas of China, based on dynamic time history analysis and probabilistic seismic demand analysis. Considering a typical single-story raw-soil structure from a rural region in China as an example, finite element models of raw-soil architecture with a colored steel roof (RACSR) were established. Based on the Mw-R strip method, data for 80 real ground motions were selected for analysis. The seismic response characteristics of raw-soil structures under earthquake action were analyzed by considering the interstory drift angle (ISDA) as a seismic performance index. A probabilistic seismic demand model of the rural raw-soil structure was established, and the vulnerability curves were obtained under different working conditions. On this basis, the effects of different raw-soil improvement methods on the seismic response and fragility of raw-soil structures were rigorously investigated, and the influence of different limit failure states on the seismic vulnerability of rural raw-soil structures was also examined.

2. The Modification of Raw-Soil Materials

To improve raw-soil materials, the physical and mechanical defects of soil should be significantly enhanced by the appropriate selection of modification materials. It is difficult to promote high-cost programs in rural areas. Raw soil has been used as a building material because of its incomparable environmental protection. Therefore, the modification materials should be carefully chosen. If the selected materials negatively impact the ecological environment, this is contrary to the intent of the original project. Therefore, in the selection of modification materials, effectiveness, cost-effectiveness, and environmental protection are the basic criteria. Standardized production technology, mechanical construction method of structural components, and chemical and physical improvement methods of raw materials are the main areas of research of raw-soil structure design and construction, especially the use of natural fiber and domestic/industrial waste to improve the raw soil. Wang [32] added recycled starch, quicklime, and hydrated lime to the raw soil, which were used in raw soil in mountainous areas. The compression and shear properties of the modification materials were tested and analyzed, and the results showed that the recycled starch was suitable for addition to the original soil. Liu et al. [33] used red clay in northern Sichuan of China as raw-soil material, polypropylene fiber, and gravel as physical modification materials, and emulsion powder as a chemical dispersion modification material to test the composite improvement effect. Based on comparative studies of raw and modified soil [33], this study focused on limestone-calcined clay-cement-modified raw soil.

Limestone-calcined clay cement (LC3) is a novel ternary blend. It is a new type of cement made from calcined clay (30%), limestone (15%), and gypsum (5%), instead of (50%) clinker in traditional cement. Coal gangue is a type of waste generated in the production of coal mining and tunneling. Owing to its low calorific value and high ash, it is an industrial solid waste that is difficult to use. Coal gangue has a specific calorific value, contains a certain amount of silicon aluminum components, which is similar to clay, and can be used as an alternative raw material to clay. The annual global cement production is 4 billion tons. China accounts for 58.13% of the world's consumption. The production of 1t of cement generates 0.82t of carbon dioxide [34]. Therefore, the replacement of part of cement clinker with coal gangue can significantly reduce the carbon footprint in cement production. Coal gangue can partly or completely replace clay when improving and strengthening soil. If further sintered into bricks or lightweight aggregate, coal can be saved and the weight of buildings can be reduced by 20%.

The cement used in experiments was coal gangue-based LC3, which was mixed in a laboratory with Portland cement clinker, coal gangue powder, limestone powder, and dehydrated gypsum in the proportion of 50 : 30 : 15 : 5. Portland cement clinker was made by grinding an undisturbed clinker ball using a ball mill for 1 h, followed by filtering using a 200-mesh standard sieve, as shown in Figure 1(a). Coal gangue powder was calcined in an industrial furnace at 850°C and then filtered using a 200-mesh standard sieve, as shown in Figure 1(b). Limestone is white crystalline powders with a content (CaCO3) of more than 99.0%, which meets the requirements of the national standard of China. Dehydrated gypsum is white crystalline powders with a content (CaSO4·2H2O) of more than 99.0%, which meets the requirements of the national standard of China. The kaolin soil produced in the factory is shown in Figure 1(c). The chemical composition of each component is listed in Table 1. Kaoling soil has a specific gravity (Gs) of 2.71, a liquid limit of 36.6%, a plastic limit of 21.2%, a water content (ω) of 0.12%, and a pH value of 7.83.

A batch of coal gangue-based LC3 specimens with different cement contents and curing ages were prepared for mechanical and microstructural tests. First, the kaolin soil and the prepared coal gangue-based LC3 were fully mixed in proportion, the water content was maintained at 38%, and a cement mortar mixer was used to perform uniform mixing. Second, the cement soil was added to the mold after mixing and then placed on a vibration table for compaction. The surface was scraped, covered with a layer of preservative film, and allowed to stand at 20 ± 5°C for 48 h. Finally, the samples were demolded, then placed in a self-sealing bag, and cured in water at 20 ± 1°C. The curing periods were 7 d, 14 d, 28 d, and 60 d.

The unconfined compressive strength (UCS) test is an index that summarizes the engineering properties of cement-stabilized soil. This strength characteristic represents the ability to carry vertical loads under no side limit pressure conditions, and the results can be used to calculate mechanical parameters such as elastic modulus. The loading rate was controlled during the test, and the load was applied evenly to the specimens at 0.03–0.15 kN/s until damage occurred, and the failure load was recorded accurately to ±0.01 kN. The UCS test of coal gangue-based LC3 soil blocks (70.7 mm × 70.7 mm × 70.7 mm) with different cement contents (5%, 10%, 15%, and 20%) and curing ages (7 d, 14 d, 28 d, and 60 d) was performed. Figure 2 shows the test process and some LC3 soil sample failure modes. In this study, the test data at 28 days of curing and 5% content, which is similar to that of fiber-reinforced soil, were selected to analyze seismic vulnerability. The Australian Native Building Code HB 195 stipulates that the elastic modulus of native building materials should be higher than 200 MPa. The test results show that raw clay and modified and LC3 soil met this condition. When the LC3 content was 5%, the cement strength satisfied the requirements of the standard project.

The microstructure of the cement soil was studied using scanning electron microscopy. The test block was broken after the UCS test, and several fragments were cut using a knife. The selected observation surface area was approximately 5 mm × 5 mm and 3 mm thick and dried for standby. Before the test, the surface of the soil sample was coated with a thin layer of gold to increase its conductivity. To observe the overall structure and local microstructural changes of the test sample, the magnification of the SEM image was set at 1000 and 5000 times, respectively. Figure 3 shows that the particles of the LC3 soil have small pores that are evenly distributed and form dense lumps, making the sample uniform and dense.

3. The Description of the Structural Models

3.1. Finite Element Models of the Raw-Soil Structures

In the modeling scheme of the numerical simulation, a three-slot symmetrical layout of the raw-soil architecture with a colored-steel roof (RACSR) was utilized. As shown in Figure 4, the investigated raw-soil structures have a dimension of 13.9 m × 6.4 m. The axis of the bay was 4.5 m × 6 m, the height of the soil transverse wall was 4.28 m, the height of the longitudinal wall was 3 m, and the height of the back wall without a window was 3 m. The roof truss was a symmetrical double-slope roof, which was used to bear the load of the ram soil cross-wall. The wall thickness of the model was 0.4 m. Windows were present in the left and right rooms. The windowsill height was 0.9 m and the window size was 1.2 m × 1.2 m. The size of the entrance door was 1.5 m × 2.1 m, and the size of the room door was 1.2 m × 2.1 m. Finite element models of the raw-soil structures were established using the ABAQUS platform. The structures were simulated using solid elements. After many trial calculations, 200 mm grids were finally selected. The model and grid of the RACSR are shown in Figure 5.

3.2. Constitutive Model of Soil Material

The Mohr–Coulomb model is the most widely used constitutive model in geotechnical engineering. Figure 6(a) shows the Mohr–Coulomb yielding model, which expresses the state when the failure line of the Mohr–Coulomb criterion is tangent to the Mohr stress circle. The parameters σ1, σ2, and σ3 are the maximum principal stress, medium principal stress, and minimum principal stress, respectively. The Mohr–Coulomb criterion can be expressed as follows:

The classical expression that is used in soil mechanics is as follows:where τ is the shearing stress, σ is the normal stress, c is the cohesion, φ is the internal friction angle of the soil, and k is a constant.

It is expressed the form of the principal stress as follows:

The yield surface shapes of the three dimensions and the π-plane shape are shown in Figures 6(b) and 6(c), and the Mohr–Coulomb yield condition can be expressed as follows:where I1 is the invariant of the first principal stress, J2 is the second deviatoric stress invariant, φ is the internal friction angle of the soil, and θσ is the stress Rhode angle.

The yield surface of the Mohr–Coulomb yield condition is an irregular pyramid in three-dimensional stress space. The yield surface in the π-plane was a hexagon. The plastic potential function gives the direction of the next step of the plastic strain, whereas the yield function gives the magnitude of the next step of the plastic strain. When the Mohr–Coulomb yield criterion is adopted, the uncorrelated flow law is used. This model is suitable for soil materials with different tensile and compressive strengths.

In this report, the raw-soil structure of typical villages and towns in China was selected as the research object, and the soil wall was plain rammed earth. The main parameters are listed in Table 2 based on the tests performed and on [33]. LC3 soil has a density of 2150 kg/m³, a Poisson’s ratio of 0.30, a cohesion of 550 kPa, an internal friction angle of 50°, and an elastic modulus of 720 MPa. The colored-steel roof has a density of 10 kg/m³, a Poisson’s ratio of 0.25, and an elastic modulus of 2-E8 MPa.

4. Nonlinear Seismic Response Analysis of Raw-Soil Structures

To obtain the dynamic characteristics of the RACSR structures, finite element models were established to calculate the vibration modes and natural frequencies of interest. A modal analysis was conducted using ABAQUS to obtain the fundamental modes of the RACSR models. The first three modal shapes are shown in Figure 7(a), and a summary of the frequency and modes that were utilized are listed in Table 3. Figure 7 indicates that the first two modes are translational motions, and the third mode is a twisting motion.

Given that ground motion is a very complex natural phenomenon, it has a high degree of uncertainty. Therefore, the selection of appropriate ground motion for dynamic time history analysis of a structure is necessary to achieve reliable calculated results. The specific earthquake types and classification criteria are presented in Table 4; the earthquakes were classified into five categories (SMSR, LMSR, SMLR, LMLR, and NEAR) according to the magnitude and epicentral distance. Referring to Zhang et al. [35] and using the magnitude-epicentral distance (Mw-R) band method, the ground motion was selected in a wide Mw-R range considering the influence of near-fault ground motion on structural vulnerability. A total of 80 real ground motion records were selected from the strong ground motion database of the Pacific Earthquake Engineering Research Center (PEER). Figure 8 shows the Mw-R distribution of the selected ground motion records from PEER. These records include all types of earthquakes, including 11 SMSR records, 10 SMLR records, 25 NEAR records, 22 LMSR records, and 12 LMLR records. Therefore, the uncertainty of the earthquake was considered in the simulation.

The 80 real ground motions were input to analyze the nonlinear seismic response of RACSR for different types of soil materials. The maximum interlayer displacement of raw-soil RACSR ranged from 2.3 mm to 72.4 mm. Figure 9 shows a representative seismic wave (LMSR, G03-000, PGA = 0.547 g) and the corresponding roof displacement of the RACSR, where PGA denotes the peak ground acceleration. Figure 9(b) shows that the maximum displacement for the raw soil, modified soil, and LC3 soil and the positive and negative maximum values are 0.057, 0.052, and 0.046 m and −0.0380, −0.0366, and −0.0362 m, respectively. The corresponding displacement control ratios compared to the raw-soil RACSR were 8.77%, 19.29%, 3.68%, and 4.74%, respectively.

Table 5 lists 25 out of 80 maximum values of ISDA and the decreasing ratio (DR). ISDA-1, ISDA-2, and ISDA-3 are the interstory drift angle of raw-soil RACSR, modified soil RACSR, and LC3 soil RACSR, respectively. The DR is a ratio that describes the decreasing amplitude. DR-1 and DR-2 are the decreasing amplitude ratios of ISDA-1 and ISDA-2, and ISDA-1 and ISDA-3, respectively.

As shown in Table 5, with the increase in the magnitude and PGA and the decrease in the epicentral distance, the maximum of the ISDA increases. When the magnitude and epicentral distances are different, the ISDA of the structure under the action of seismic waves with similar PGA is very close. For the same magnitude and different epicentral distances, the PGA increases in a stepwise manner, and the ISDA increases with the increase in the PGA. This shows that PGA is one of the main factors that affect the ISDA of raw-soil structure buildings. According to Figure 9(b), the displacement of the post-peak time shows that the global integrity of the stone-modified soil structure is different from that of the other two soils. For most seismic waves, the ISDA decreases from ISDA-1 to ISDA-3. This proves that the improvement and modification of raw soil can enhance the rigidity of the RACSR.

5. Seismic Vulnerability Analysis of Raw-Soil Structures

5.1. Probabilistic Seismic Demand Model

The probabilistic seismic demand (PSD) model characterizes the relationship between the engineering demand parameter (EDP) and ground motion intensity measure (IM). In the seismic vulnerability analysis of structures, the commonly used engineering demand parameters (EDPs) mainly include the top floor displacement, maximum value of ISDA, and the damage index. The ISDA is simple and intuitive in the analysis process and can better reflect the seismic demand of a structure. In the present study, the maximum ISDR and PGA were selected as the EDP and IM, respectively, for the PSD analysis of the three structures. Taking 80 real ground motions as the seismic excitation, the nonlinear dynamic time history analysis of RACSR was performed, and the corresponding relationship between seismic intensity parameters and the engineering demand parameters (IMi-EDPi) under different conditions was obtained. Cornell et al. [36] proposed that the relationship between EDP and IM should satisfy the following logarithmic linear relation:

The PSD model under each condition was obtained by fitting the results from the nonlinear dynamic time history analysis of the structure. The logarithmic standard deviation of the seismic demand can be expressed as follows:where is the logarithmic standard deviation of earthquake demand, N is the number of sample points in the regression analysis, Di is the peak value of the ith seismic demand, IMi is the PGA of the ith ground motion, and a and b are the regression parameters.

Figures 1012 show the PSD models of the RACSR under the action of a single main earthquake. Table 6 lists the mathematical expressions and related parameters of the PSD models under different conditions.

5.2. Vulnerability Analytical Method

The fragility function can be expressed as follows [37]:where D and C are the seismic demand and structural capacity, respectively, IM is the ground motion intensity measure, μD and μC are the medians of D and C, respectively, and and σC are the standard deviations corresponding to D and C, respectively.

The national standard of the People’s Republic of China (GB/T24335-2009) outlines a description of the degree of structural damage for five grades of soil and stone structure houses. The building damage was classified into five grades, as shown in Table 7.

The seismic vulnerability of a structure refers to the conditional probability that the structure reaches or exceeds a certain limit state under the action of earthquakes of different strengths; thus, it describes the probability distribution of all limit states of the structure. Therefore, the definition of the limit state is an important component of structural vulnerability analysis. The PGA was used as the IM and the ISDA was used as the EDP. According to the categorization of failure grade and performance level in Table 7, and referring to the ISDA limit values of each performance point for unreinforced clay brick walls in reference, four limit states of the raw-soil structure are defined: slight damage (LS-1), moderate damage (LS-2), severe damage (LS-3), and collapse (LS-4). The ISDA corresponding to the four limit states is listed in Table 8. Under the limit state of slight damage (LS-1), it has a grade I intact state, and above the limit state of collapse (LS-4), it has a grade V collapsed state.

5.3. Vulnerability Results and Discussion

Figure 13 shows a comparison of the RACSR vulnerability curves for different limit states and earthquakes. The vulnerability curves of the raw-soil RACSR, modified soil RACSR, and LC3 soil RACSR model are presented for 80 real ground motions. The damage probability values of the raw-soil RACSR that was subjected to an earthquake were significantly higher than those of the modified soil RACSR and LC3 soil RACSR that were directly subjected to an earthquake of the same intensity.

As shown in Figure 13(a), considering the case of PGA = 0.1 g as an example, the damage probabilities of raw-soil RACSR corresponding to the limit states of slight damage (LS-1), moderate damage (LS-2), severe damage (LS-3), and collapse (LS-4) are 99.19%, 93.89%, 77.37%, and 51.94%, respectively. The damage probabilities of LC3 soil RACSR corresponding to the limit states of slight damage (LS-1), moderate damage (LS-2), severe damage (LS-3), and collapse (LS-4) are 95.22%, 79.1%, 50.72%, and 24.76%, respectively. At PGA = 0.2 g, the damage probabilities of raw-soil RACSR corresponding to the limit states of slight damage (LS-1), moderate damage (LS-2), severe damage (LS-3), and collapse (LS-4) were 99.91%, 98.84%, 93.32%, and 78.06%, respectively. The damage probabilities of LC3 soil RACSR corresponding to the limit states of slight damage (LS-1), moderate damage (LS-2), severe damage (LS-3), and collapse (LS-4) are 99.26%, 94.31%, 78.49%, and 53.54%, respectively. Thus, compared to the raw-soil RACSR model, the LC3 soil reduced the values of the damage limit states of LS-1, LS-2, LS-3, and LS-4 by factors of 1.04, 1.19, 1.53, and 2.1, respectively, for the same earthquake intensity. As shown in Figure 13(b), considering the case of PGA = 0.1 g as an example, the damage probabilities of modified soil RACSR corresponding to the limit states of slight damage (LS-1), moderate damage (LS-2), severe damage (LS-3), and collapse (LS-4) are 93.67%, 75.53%, 46.73%, and 22.19%, respectively. At PGA = 0.2 g, the damage probabilities of the modified soil RACSR corresponding to the limit states of slight damage (LS-1), moderate damage (LS-2), severe damage (LS-3), and collapse (LS-4) are 99.58%, 96.42%, 84.80%, and 63.47%, respectively. That is, compared to the raw-soil RACSR model, the modified soil reduced the values of the damage limit states of LS-1, LS-2, LS-3, and LS-4 by factors of 1.01, 1.03, 1.11, and 1.23, respectively, for the same earthquake intensity.

Figure 14 compares the fragility curves of the three models (raw-soil RACSR, modified soil RACSR, and LC3 soil RACSR) for different limit states (LS-1, LS-2, LS-3, and LS-4). As shown in Figure 14(a), taking PGA = 0.05 g as an example, the damage probability values of the slight damage limit state (LS-1) of the three models are 98.9%, 81.8%, and 89.2%, respectively. As shown in Figure 14(b), taking PGA = 0.1 g as an example, the damage probability values of the moderate damage limit state (LS-2) of the three models are 93.9%, 75.5%, and 79.1%, respectively. Figure 14(c) shows that taking PGA = 0.1 g as an example, the damage probability values of the severe damage limit state (LS-3) of the three models are 77.4%, 46.7%, and 50.7%, respectively. When the PGA = 0.2 g, the three models are 93.0%, 84.8%, and 78.5%, respectively. Figure 14(d) shows that taking PGA = 0.1 g, the damage probability values of the collapse limit state (LS-4) of the three models are 51.9%, 22.2%, and 24.7%, respectively. Taking PGA = 0.2 g, the damage probability values of the collapse limit state (LS-4) of the three models are 78.1%, 63.5%, and 53.5%, respectively. When PGA = 0.4 g, the damage probability values of the collapse limit state (LS-4) of the three models are 93.6%, 93.3%, and 83.7%, respectively.

6. Conclusion

In this study, the seismic fragility of raw-soil structures that are representative of rural areas of China was investigated via numerical simulations. The effects of raw and modified soil on the fragility of rural structures were investigated in detail, considering the influence of LC3 soil. The main conclusions are as follows:(1)Coal gangue-based LC3 was mixed with silicate cement clinker, coltan powder, limestone powder, and dehydrated gypsum in proportion 50 : 30 : 15 : 5. The LC3 soil strength with 5% or 10% cement can meet the demands of standard and engineering projects. The replacement of some cement clinker with coltan can significantly reduce the carbon footprint in cement production, and national policies of China encourage the reuse of solid waste. This type of modification material is effective, cheap, and environmentally friendly and can be used as a rural modification material(2)Compared to the raw-soil RACSR model, the roof displacement and ISDA of the modified soil and LC3 soil RACSR models are both reduced to varying degrees, and the damage probability values for different damage limit states of the RACSR model are significantly lower. Therefore, it is necessary to improve the soil mechanics in rural raw-soil structures to enhance their structural integrity, reduce their seismic responses, and effectively mitigate the risk of severe structural damage or collapse(3)For the seismic analysis of the raw-soil RACSR model, PGA is one of the main factors that influence the ISDA of raw-soil structure buildings. With the increase in magnitude and PGA and the decrease in the epicentral distance, the maximum value of the ISDA increased. When the magnitude and epicentral distance are different, the ISDA of the structure under the action of seismic waves with similar PGA is very close. The magnitude is the same, the epicentral distance is different, the PGA increases stepwise, and the ISDA increases with increasing PGA(4)With an increase in the PGA, the damage probability value of each damage limit state of the structures increases. Taking PGA = 0.1 g, the damage probabilities for raw-soil RACSR corresponding to the limit states of LS-1, LS-2, LS-3, and LS-4 are 99.19%, 93.89%, 77.37%, and 51.94%, respectively. At PGA = 0.2 g, the damage probabilities are 99.91%, 98.84%, 93.32%, and 78.06%. Compared to the raw-soil RACSR model, the modified soil and LC3 soil reduced the damage probability. For limit states of LS-1, LS-2, LS-3, and LS-4, the modified soil is reduced by factors of 1.01, 1.03, 1.11, and 1.23, respectively. The LC3 soil is reduced by factors of 1.04, 1.19, 1.53, and 2.1, respectively, under the same earthquake intensity of 0.2 g as that of LS-1 to LS-4(5)Relatively few modification methods were used in the study to improve the raw soil. The displacement of the post-peak time shows that the global integrity of the modified soil structure should be investigated in detail based on the soil failure mechanism. In addition, it is necessary to study the seismic performance and fragility of structures of various forms. This work serves as a theoretical reference for seismic design and performance improvement of raw-soil structures

Data Availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this study.

Acknowledgments

The authors are grateful for the financial support from the National Key R&D Program of China (2018YFD1100402) and the Project of Liaoning Provincial Department of Science & Technology (2020-BS-161).