Seismic Vulnerability Analysis of Rural Modified Raw-Soil Structures

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. ,e 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. ,e 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). ,emaximum interstory drift angle (ISDAmax) was lower in the LC3 soil model and the modified soil model compared to the rawsoil model. ,e 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.


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. ey 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. e 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. e 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. e results showed that the proposed plaster-adobewall 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. e 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. erefore, 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. e 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. e 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 rawsoil structures. In particular, raw-soil structures in rural areas are typically not optimally designed, and many selfbuilt 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 2 Shock and Vibration 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. e 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.

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.
erefore, 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. erefore, 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. e 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. e annual global cement production is 4 billion tons. China accounts for 58.13% of the world's consumption. e production of 1t of cement generates 0.82t of carbon dioxide [34]. erefore, 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%. e 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 (CaCO 3 ) of more than 99.0%, which meets the requirements of the national standard of China. Dehydrated gypsum is white crystalline powders with a content (CaSO 4 ·2H 2 O) of more than 99.0%, which meets the requirements of the national standard of China. e kaolin soil produced in the factory is shown in Figure 1(c). e chemical composition of each component is listed in Table 1. Kaoling soil has a specific gravity (Gs) of 2.71, a liquid limit (W L ) of 36.6%, a plastic limit (W P ) 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. e 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. e curing periods were 7 d, 14 d, 28 d, and 60 d. e unconfined compressive strength (UCS) test is an index that summarizes the engineering properties of cement-stabilized soil. is 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. e 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. e 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 Shock and Vibration 3 to analyze seismic vulnerability.
e Australian Native Building Code HB 195 stipulates that the elastic modulus of native building materials should be higher than 200 MPa. e 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.
e microstructure of the cement soil was studied using scanning electron microscopy. e test block was broken after the UCS test, and several fragments were cut using a knife. e 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.

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. e 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. e roof truss was a symmetrical double-slope roof, which was used to bear the load of the ram soil cross-wall. e wall thickness of the model was 0.4 m. Windows were present in the left and right rooms. e windowsill height was 0.9 m and the window size was 1.2 m × 1.2 m. e 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. e structures were simulated using solid elements. After many trial calculations, 200 mm grids were finally selected. e model and grid of the RACSR are shown in Figure 5.

Constitutive Model of Soil Material.
e 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. e parameters σ 1 , σ 2 , and σ 3 are the maximum principal stress, medium principal stress, and minimum principal stress, respectively. e Mohr-Coulomb criterion can be expressed as follows: e 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: (3) e 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:

Shock and Vibration
where I 1 is the invariant of the first principal stress, J 2 is the second deviatoric stress invariant, φ is the internal friction angle of the soil, and θ σ is the stress Rhode angle. e yield surface of the Mohr-Coulomb yield condition is an irregular pyramid in three-dimensional stress space. e yield surface in the π-plane was a hexagon. e 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. is 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. e 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. e 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.

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. e 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.  Shock and Vibration Given that ground motion is a very complex natural phenomenon, it has a high degree of uncertainty. erefore, the selection of appropriate ground motion for dynamic time history analysis of a structure is necessary to achieve reliable calculated results. e 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. ese records include all types of earthquakes, including 11 SMSR records, 10 SMLR records, 25 NEAR records, 22 LMSR records, and 12 LMLR records. erefore, the uncertainty of the earthquake was considered in the simulation. e 80 real ground motions were input to analyze the nonlinear seismic response of RACSR for different types of soil materials. e 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. e corresponding displacement control ratios compared to the rawsoil 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. e 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. is shows that PGA is one of the main factors that affect the ISDA of rawsoil 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.
is proves that the improvement and modification of raw soil can enhance the rigidity of the RACSR.

Probabilistic Seismic Demand
Model. e 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   Figure 7: e first three modal shapes of the RACSR for the LC3 soil. Stereogram and plan view of (a) the first modal shape, (b) the second modal shape, and (c) the third modal shape. e 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 (IM i -EDP i ) under different conditions was obtained. Cornell et al. [36] proposed that the relationship between EDP and IM should satisfy the following logarithmic linear relation: e PSD model under each condition was obtained by fitting the results from the nonlinear dynamic time history analysis of the structure. e logarithmic standard deviation of the seismic demand can be expressed as follows: where σ D|IM is the logarithmic standard deviation of earthquake demand, N is the number of sample points in the regression analysis, D i is the peak value of the ith seismic demand, IM i is the PGA of the ith ground motion, and a and b are the regression parameters. Figures 10-12 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.

Vulnerability Analytical Method.
e 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 σ D|IM and σ C are the standard deviations corresponding to D and C, respectively. e 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. e building damage was classified into five grades, as shown in Table 7.
e 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. erefore, the definition of the limit state is an important component of structural vulnerability analysis. e 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). e 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. Figure 13 shows a comparison of the RACSR vulnerability curves for different limit states and earthquakes. e vulnerability curves of the raw-soil RACSR, modified soil RACSR, and LC3 soil RACSR model are presented for 80 real ground motions. e 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.

Vulnerability Results and Discussion.
As shown in Figure 13(a), considering the case of PGA � 0.1 g as an example, the damage probabilities of rawsoil 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   As shown in Figure 13 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.

II
Slight damage e bearing wall has no damage or some slight cracks in some cases. e roof and arch are basically intact. Some non-load-bearing members have slight damage or some obvious damage, such as some non-loadbearing walls have slight cracks, some have obvious cracks, the gable slightly flashes outside, the roof tiles slide, etc. e basic function of the structure is not affected, and it can be used with little or no repair.

III
Moderate damage Most of the load-bearing walls have slight cracks, some walls have obvious cracks, some walls have serious cracks, and many cave arches have cracks. ere are obvious cracks in some roofs and vaults. Some nonload-bearing components have obvious damage, such as many parts of the wall plastering fall off, some roof tiles fall off, etc. e basic function of the structure is affected to some extent, and it can be used after repair.

IV Severe damage
Most of the load-bearing walls have obvious cracks, and some of them have serious damage, such as wall dislocation, breakage, internal or external tilt, or local collapse. Uplift or collapse of roof or vault, local collapse, or the entire structure is obviously inclined. e basic function of the structure is seriously affected, and even part of the function is lost, which is difficult to repair or has no repair value.

V Collapse
Most of the walls are seriously broken or collapsed, and the roof or vault is seriously damaged and collapsed. e structure is on the verge of collapse or completely collapsed. e function of the structure no longer exists and there is no possibility of repair.

Conclusion
In this study, the seismic fragility of raw-soil structures that are representative of rural areas of China was investigated via numerical simulations. e effects of raw and modified soil on the fragility of rural structures were investigated in detail, considering the influence of LC3 soil. e 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. e LC3 soil strength with 5% or 10% cement can meet the demands of standard and engineering projects. e 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. is 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. erefore, 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. e 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. e 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. e 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. is work serves as a theoretical reference for seismic design and performance improvement of raw-soil structures Data Availability e data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this study.