Abstract

The self-stability height of the foundation pit sidewall is an important criterion for evaluating the safety degree and designing the supporting structure. The strength reduction elastic-plastic finite element numerical calculation method has been adopted in this paper. Based on comparative analysis of the stability characteristics for deep foundation pit in binary strata of upper soil and lower rock under multiple working conditions, the potential fracture surface of deep foundation pit and the evolution law of corresponding safety factor have been revealed under different Hs and H. A new idea that the vertical soil sidewall height (Hs) and the vertical rock sidewall height (Hr) are used as two independent evaluation indexes, respectively, for deep foundation pit stability in binary strata of upper soil and lower rock has been put forward. The distribution characteristics and variation law of Hs0 and Hr0 under different Hs and different H have been revealed, respectively. The spatial distribution map of the self-stabilizing height for deep foundation pit vertical sidewall in upper soil and lower rock binary stratum has been constructed, and the mathematical fitting equation between Hr0 and Hs has been obtained. Finally, combined with the implementation effect of the deep foundation pit project of Ningxia Road Station for Qingdao Metro Line 3, the rationality of the conclusions is verified. The research results provide theoretical basis for quickly determining the self-stability characteristics of foundation pit vertical sidewall.

1. Introduction

With the accelerating process of urbanization in China, a type of deep foundation pit whose bottom is above the interface between soil stratum and rock stratum comes into being. The unique feature of this deep foundation pit type is that the strata within its excavation depth are composed of soil and rock from top to bottom, and the difference of the main mechanical parameters such as cohesion and elastic modulus between the two strata is almost thousands of times, such as Wusi Square Station foundation pit of Qingdao Metro Line 3 [1], Anshan Road Station foundation pit of Qingdao Metro Line 4 [2, 3], and Chashan station foundation pit of Shenzhen Metro Line 7 [4]. In this paper, this type of deep foundation pit is uniformly called deep foundation pit in the binary strata of upper soil and lower rock.

The self-stability height of foundation pit sidewall is an important criterion for evaluating the safety degree and designing the supporting structure. As early as the 1860s, Culann [5] deduced the formula for calculating the vertical sidewall critical instability height according to the balance condition of force. After that, Terzaghi et al. [610] modified the Culann method and formed a series of theoretical formulas for the vertical sidewall critical instability height. Pufahl [11] deduced the self-stable critical height of the foundation pit vertical excavation slope by calculating and analyzing the condition that the sliding wedge is in the limit equilibrium state. Zhang et al. [12] adopted the beam column mechanical model under the action of self- weight and obtained the critical instability height of vertical layered rock slope. Liu et al. [13, 14] constructed the instability calculation model of vertical layered rock slope based on statistical damage model and obtained the critical height, based on Euler compression bar theory and capacity method. Based on the upper bound theorem of limit analysis method, Chen [15] and Wang et al. [16] deduced the calculation formula of double-layer slope critical height at any layered position under different failure mechanisms and obtained the linear variation law of double-layer slope critical height with layered depth under different conditions. Li et al. [17, 18] studied the overall stability and failure mode of deep foundation pit slope in binary stratum of soil and fully weathered rock, and soil and strongly weathered rock, respectively.

Strength reduction method (SRM) was first proposed by Zienkiewicz et al. in 1975 and used for slope stability analysis [19]. And then, many scholars have done a lot of work on the applicability and feasibility of SRM in slope stability analysis, such as Ugai [20], Matsui and San [21], Griffiths and Lane [22], and Dawson [23] et al. For the whole safety storage of the slope based on the double safety factors, nonproportional correlative reduction finite element method with the contribution of the shear strength parameter to the sliding resistance force as weight is proposed by Xue et al. [24]. A parallel-local strength reduction method is proposed by Zhang et al. [25], which is based on the idea of the local strength reduction method. For the heterogeneous slope, a novel hierarchical multiscale strength reduction method is proposed by Meng et al. [26]. For the stability evaluation of rock slope, the strength reduction method considering the ubiquitous joint model is proposed by Liu et al. [27]. Based on the strength reduction method, the changing positions of the slope sliding surface during the filling process are studied by Yuan et al. [2830].

The self-stability of deep foundation pit sidewall has been taken under serious consideration, and the SRM method has been used to evaluate it as well. However, the vertical sidewall of deep foundation pit in the binary strata of upper soil and lower rock as the independent factor is less considered. There are few reports on the research based on the self-stable height of vertical sidewall for deep foundation pit. In this paper, the strength reduction elastic-plastic finite element numerical analysis method is adopted to calculate and analyze the potential fracture position and FOS of the vertical sidewall for deep foundation pit in the binary strata of upper soil and lower rock. And then, the distribution characteristics and variation law of Hs0 and Hr0 under different Hs and different H have been revealed, respectively. The spatial distribution map of the self-stabilizing height for deep foundation pit vertical sidewall in upper soil and lower rock binary stratum has been constructed, and the mathematical fitting equation between Hr0 and Hs has been obtained. Finally, combined with the implementation effect of the deep foundation pit project of Ningxia Road Station of Qingdao Metro Line 3, the rationality of the conclusions is verified.

2. Methodology

2.1. Strength Reduction Method (SRM)

The strength reduction method is a method combining strength reduction technology, limit equilibrium principle, and elastic-plastic finite element calculation principle. Firstly, the stress and deformation state under the condition of original material parameters is calculated. Then, the material strength parameters c and tanφ are reduced simultaneously according to formula, and a (1) set of new strength parameters and are obtained and calculated as new material strength parameters. Finally, the critical fracture surface is obtained by continuously adjusting the reduction factor k until the material is in the limit equilibrium state. At this time, the reduction factor k of the material is the stability FOS.where k is the reduction factor.

The calculation principle of SRM [31] is shown in Figure 1.

2.2. Implementation Process and Example Analysis

The numerical calculation of deep foundation pit stability analysis adopts Midas GTS analysis software and is modeled according to 3D solid element. The upper boundary of the calculation model is taken to the surface. The lower boundary is taken to 3 times of the excavation depth below the base, and the horizontal boundary is more than 3 times of the excavation depth from the sidewall of the foundation pit. The thickness is taken as 1.0 m. Gravity constraints and model boundary constraints are applied. A uniform ground load of 20 kPa is applied to the upper surface of the calculation model without considering the effects of groundwater and earthquake. The grid size of the foundation pit excavation part is set as 0.5 m, and the other grid sizes are set as 2.0 m. The improved Mohr–Coulomb constitutive model is adopted for rock and soil mass, and the effects of rock and soil self-weight, ground load, and boundary constraints are considered in the initial stress. The calculation model is shown in Figure 2.

The implementation process is completed by six key steps. Now, take the deep foundation pit in the binary strata of upper soil and lower rock with the upper soil stratum thickness of 10 m and the excavation depth of 15 m as an example.Step 1: import the drawn calculation model into Midas GTS software, as shown in Figure 3(a).Step 2: expand the whole part of the model by 1.0 m along the z-axis, and then control the size of the foundation pit location entity and other parts, as shown in Figure 3(b).Step 3: divide the grid into 4 parts: soil stratum excavation part, rock stratum excavation part, soil stratum unexcavated part, and rock stratum unexcavated part. The grid size of the excavation part is 0.5 m × 0.5 m, and the other parts are 2.0 m × 2.0 m, as shown in Figure 3(c).Step 4: apply self-weight, boundary constraints, and 20 kPa ground uniformly distributed load, as shown in Figure 3(d).Step 5: carry out analysis and calculation. The solution type is set as SRM calculation of slope stability. The solution type is set as the SRM calculation model of slope stability. Passivate the grid distributed in the excavated soil stratum and rock stratum. Activate the grid of the unexcavated soil stratum and rock stratum. At the same time, activate the self-weight of rock and soil mass, boundary constraint, and uniformly distributed load.Step 6: run the calculation, and get the results. Get the foundation pit stability FOS and the corresponding distribution diagram of potential fracture surface, as shown in Figure 3(e).

3. Study on Stability Characteristics of Deep Foundation Pit

3.1. Analysis Model and Calculation Parameters

The deep foundation pit in binary strata of upper soil and lower rock is that the strata within their excavation depth are composed of soil and rock from top to bottom, and the difference of the main mechanical parameters such as cohesion and elastic modulus between the two strata is almost thousands of times. For the convenience of description, the part of foundation pit vertical sidewall in the soil stratum is called foundation pit vertical soil sidewall, and its height is called the vertical soil sidewall height (Hs). Correspondingly, the part of the foundation pit vertical sidewall in the rock stratum is called the foundation pit vertical rock sidewall, and its height is called the vertical rock sidewall height (Hr), as shown in Figure 4.

According to the statistical analysis of indoor test data of rock and soil mass at the site of deep foundation pit project of Ningxia Road Station of Qingdao metro, combined with the latest research results, we perform the stability analysis of vertical sidewall of deep foundation pit in upper soil and lower rock binary stratum, and the values of physical and mechanical parameters of rock and soil mass are shown in Table 1.

3.2. Calculation Results and Analysis

The thickness of upper soil stratum for deep foundation pit in binary strata of upper soil and lower rock is taken as 4 types, which are 5 m, 10 m, 15 m, and 20 m in turn. The foundation pit excavation depth of each stratum type is taken as 5 working conditions, which are 5 m, 10 m, 15 m, 20 m, and 30 m in turn. The distribution of foundation pit potential fracture surface under different working conditions is shown in Figure 5, and the calculation results of FOS are shown in Table 2 and Figure 6.

The following conclusions can be drawn from Figures 5 and 6: (1) the potential fracture surfaces of deep foundation pit in binary strata of upper soil and lower rock are located within the upper soil stratum excavation depth. The potential fracture surface penetrates from the top of foundation pit to the interface between soil stratum and rock stratum, which is approximately in the circular arc. The failure mechanism of foundation pit is basically consistent with that of the homogeneous soil stratum slope. (2) When the foundation pit is excavated in the upper soil stratum, FOS decreases with the increase of excavation depth; when the excavation depth completely enters the lower rock stratum, FOS is approximately unchanged. (3) The stability is mainly controlled by the excavation height of upper soil stratum and is less affected by the excavation height of lower rock stratum. As the traditional index to measure the stability of deep foundation pit, excavation depth is not suitable for deep foundation pit in binary strata of upper soil and lower rock. It has important theoretical value and practical significance to study vertical soil sidewall height and vertical rock sidewall height as the independent evaluation indexes of deep foundation pit in binary strata of upper soil and lower rock.

4. Study on Self-Stability Height of Vertical Sidewall for Deep Foundation Pit

4.1. Evaluation Criteria

FOS is an important index to measure the stability of deep foundation pit, which has been widely used in geotechnical engineering [32]. Due to the complexity of foundation pit itself, the uncertainty of rock and soil characteristics, the limitation of human cognitive level, and other factors, FOS with certain reserve is generally used for stability evaluation of deep foundation pit. According to the national standard of the People’s Republic of China “Technical code for building slope engineering” (GB 50330-2013), the slope is divided into two types: temporary slope and permanent slope. According to the stability coefficient, it is divided into 4 states: stable, basically stable, basically unstable, and unstable (as shown in Table 3). The slope stability safety factor (FST) is taken as shown in Table 4.

In this paper, the deep foundation pit vertical sidewall is considered as the conservative situation of Class I permanent slope, and FOS = 1.35 is taken as the foundation pit slope self-stability criterion. That is, when FOS of foundation pit vertical soil sidewall (or rock sidewall) is greater than or equal to 1.35, it is the foundation pit vertical soil sidewall (or rock sidewall) to meet the self-stability requirements. On the contrary, when FOS is less than 1.35, it does not meet the self-stability requirements. In this paper, the critical height when the foundation pit vertical soil wall meets self-stable state is called self-stable critical height of vertical soil wall (Hs0); similarly, the critical height when the foundation pit vertical rock wall meets the self-stable state is called the self-stable critical height of vertical rock wall (Hr0).

4.2. Study on the Self-Stabilizing Height of Vertical Soil Sidewall

The thickness of upper soil stratum is taken as 11 types, which are 5 m, 6 m, 7 m, 8 m, 9 m, 10 m, 11 m, 12 m, 13 m, 14 m, and 15 m in turn. The excavation depth of the foundation pit is the same as the thickness of the upper soil stratum, and the calculation parameters and calculation process are the same as those in Section 3.1. The calculation results of FOS under different Hs are shown in Table 5 and Figure 7.

The following conclusions can be drawn from Table 5 and Figure 7: (1) FOS of foundation pit vertical soil sidewall decreases with the increase of soil wall excavation height. According to the criterion of FST = 1.35, Hs0 is 6.92 m. (2) According to the relationship between Hs and Hs0, the stable state of the vertical soil wall can be quickly determined: when Hs ≤ Hs0, soil wall excavation meets the self-stability requirements. On the contrary, when Hs > Hs0, soil wall excavation cannot meet the requirements of self-stability.

4.3. Study on Self-Stabilizing Height of Vertical Rock Sidewall
4.3.1. Computational Assumptions

It can be seen from Section 3 that it is impossible to directly obtain the potential fracture surface and its corresponding FOS of the foundation pit vertical rock sidewall. Therefore, the improved calculation model is adopted. Firstly, the action of upper soil stratum is converted into an equivalent uniformly distributed load, and the load is applied to the lower rock stratum. Then, it is calculated according to the conventional finite element strength reduction method, and the potential fracture surface of the vertical rock sidewall and its corresponding FOS are obtained.

4.3.2. Calculation Results and Analysis

The thickness of upper soil stratum in the binary strata of upper soil and lower rock is taken as 5 types, which are 0 m, 5 m, 10 m, 15 m, and 20 m in turn. The vertical rock sidewall of each stratum type is taken as 12 working conditions, which are 5 m, 10 m, 15 m, 20 m, 30 m, 50 m, 70, 90 m, 110 m, 130, 150 m,170, and 190 m in turn. The calculation parameters and calculation process are the same as those in Section 3.1. The calculation results of FOS under different working conditions are shown in Table 6 and Figure 8.

The following conclusions can be drawn from Figure 8: (1) When Hs is the same, FOS of vertical rock sidewall decreases with the increase of Hr, but the decreased amplitude decreases gradually. When Hr is the same, FOS of vertical rock wall decreases with the increase of Hs. (2) If we take FST = 1.35 as the critical point, calculated by interpolation method, Hr0 under different Hs is shown in Table 7 and Figure 9.

The following conclusions can be drawn from Figure 9: Hr0 is large as a whole, in the range of 100∼125 m, and Hr0 decreases approximately linearly with the increase of Hs.

The mathematical fitting equation between Hr0 and Hs is shown as follows:

According to the relationship between Hr and Hr0, the stable state of the vertical rock wall can be quickly determined:(1)When Hr ≤ Hr0, the rock sidewall excavation of deep foundation pit can meet the requirements of self-stability(2)when Hr > Hr0, the rock sidewall excavation of deep foundation pit cannot meet the requirements of self-stability

4.3.3. Confirmatory Analysis

The confirmatory analysis model of deep foundation pit in binary strata of upper soil and lower rock still converts the action of upper soil stratum into equivalent uniform load and applies the load directly to the lower rock stratum. At the same time, it is assumed that the sliding body is a rigid body when the foundation pit rock wall is damaged, and the sliding surface is a plane. Thus, the analysis model of vertical rock wall stability is constructed, as shown in Figure 10.

It is assumed that the foundation pit sliding body is ABC, the length of the sliding surface AC is s, and the included angle between the sliding surface and the horizontal plane is α. The sliding weight is γ. The cohesion is c and the internal friction angle is φ. The sum of the ground load and the equivalent uniformly distributed load of the upper soil stratum is q.

For ease of description, the sliding body self-weight and the overlying load are combined into force G′. The normal component of the force on the sliding surface is . The tangential component is .

When the sliding body slides, the tangential component is equal to shear strength of sliding body, as shown in the following equation:

According to the mechanical equilibrium relationship, equations (4) and (5) can be obtained.where and s are, respectively, as shown in (6) and (7).

Substituting equations (3), (6), and (7) into (4) and (5), the critical depth can be obtained when the deep foundation pit vertical rock wall is in limit equilibrium state, as shown in

In order to obtain the most dangerous sliding surface critical height of the deep foundation pit vertical rock wall, let , and (9) can be obtained

By introducing (9) into equation (8), the following equation can be obtained:

When the upper soil stratum thickness is 0 m, 5 m, 10 m, 15 m, and 20 m, can be obtained by SRM, and theoretical calculations are adopted, respectively, as shown in Table 8 and Figure 11.

It can be seen from Table 8 and Figure 11 that, compared with the theoretical calculation, the calculation results of the obtained by SRM have the same change trend. The calculation results of SRM are slightly larger overall, but there are few differences, all in the range of 5.2%∼7.4%, which verifies the reliability of the calculation results for SRM.

5. Self-Stable Spatial Distribution Characteristics and Discussion

If we take FST = 1.35 as the critical point, the deep foundation pit vertical sidewall in binary strata of upper soil and lower rock is divided into two types of spatial distribution areas: self-stable spatial distribution area and non-self-stable spatial distribution area, which are called E area and D area, respectively, in this paper. According to the stratum distribution type, it can be further divided into the self-stable height spatial distribution area of vertical soil sidewall and the non-self-stable height spatial distribution area of vertical soil sidewall (called E1 area and D1 area, respectively), and the self-stable height spatial distribution area of vertical rock sidewall and the non-self-stable height spatial distribution area of vertical rock sidewall (called E2 area and D2 area, respectively), as shown in Figure 12.

In the actual project, due to the limited excavation depth, the deep foundation pit generally will not enter the D2 area. When the thickness of the over soil stratum is less than or equal to Hs0, the D1 area disappears. That is, the foundation pit directly enters E2 area from E1 area, and the foundation pit is in a stable state as a whole. When D1 area exists, it shall be key protected during the excavation of foundation pit.

The significance of dividing the spatial distribution characteristics of self-stable height is that the unique spatial distribution characteristics of deep foundation pit self-stable area in the binary strata of upper soil and lower rock are clarified, which provides a theoretical basis for determining its reasonable support measures.

6. Typical Project

Ningxia Road Station of Qingdao Metro Line 3 is located in the south of the intersection of Nanjing Road and Ningxia road, Shibei District, Qingdao, and in the prosperous urban area along the Nanjing road from south to north. The surrounding environment of the station is complex, with an urban overpass in the north, a demolished residential community in the East, and 6-storey business building in the West. Station mileage from start to end is K9 + 362.229∼K9 + 516.979. The total length of the station main foundation pit is 154.75 m, and the excavation depth is 18.5∼20.6 m. The excavation width of the foundation pit standard section is 20.60 m, and that of the hanging section outside the equipment room with K9 + 473.279∼K9 + 516.979 is 22.60 m, as shown in Figure 13.

Overall, the project site is high in the north and low in the south, with the maximum elevation difference of 5.26 m. The geomorphic type is Piedmont erosion accumulation sloping fields. The strata distribution types from the surface top to bottom are Quaternary soil stratum (Q4), late Yanshanian intrusive granite (γ53), and lamprophyre stratum (x53) in turn. Lamprophyre is vein interspersed in granite, and cataclastic rocks are seen in contact zones of different lithology. The bedrock can be divided into strongly weathered stratum, moderately weathered stratum, and slightly (non)weathered stratum according to the weathering degree. According to the statistical analysis of the test data of undisturbed soil physical property test, rock and soil compression test, and rapid shear test, the physical and mechanical parameters of each layer of the project site are shown in Table 9.

It can be seen from Table 9 that the difference of main strength and deformation parameters such as cohesion and elastic modulus between Quaternary soil stratum and moderately weathered and slightly weathered rock stratum is thousands of times. The original rock structure of the strongly weathered zone has been damaged, and the rock core is in the state of half rock and half soil or breccia. Before disturbance, it is relatively dense and has high bearing capacity. After disturbance, the strength decreases rapidly, softens, and disintegrates, and the self-stability is poor, but its engineering properties are still slightly better than that of the Quaternary soil stratum as a whole. In order to simplify the analysis, the Quaternary soil stratum and strongly weathered rock stratum are collectively referred to as soft soil stratum. The moderately weathered rock stratum and slightly (non)weathered rock stratum are collectively referred to as hard rock stratum. Thus, the binary strata of upper soil and lower rock are constructed [33].

A total of 17 geological exploration holes are arranged in the project site, including 6 control holes and 11 general holes. The controlled drilling hole enters about 10 m below the foundation pit bottom. If there is medium∼slightly weathered rock stratum below the foundation pit bottom, it shall be drilled 3∼5 m. The general drilling hole enters about 6 m below the foundation pit bottom. If there is medium∼slightly weathered rock stratum below the foundation pit bottom, it shall be drilled 3∼5 m. The borehole located at the fault fracture zone shall open the upper and lower walls so as to reveal the basic situation of the fault. The statistical analysis results of geological survey hole survey data are shown in Table 10.

It can be seen from Table 10 that the upper soft soil stratum thickness on the project site is uneven. The soil stratum thickness of the external equipment room at K9 + 473.279∼K9 + 516.979 is about 7.80∼12.30 m, and the lower rock stratum thickness is about 5.30∼10.80 m.

Foundation pit retaining structure adopts φ800 mm@1200 mm reinforced concrete bored pile, and the soil between piles is protected by C20 net shotcrete with the thickness of 20 cm. The first internal support adopts φ690 mm@16 mm supported steel pipe, and its lower part is supported by prestressed anchor cable. The pile bottom of the retaining pile is driven into the moderately weathered rock stratum with the depth of 2.0 m. The rock shoulder with the width of 1.0 m is reserved at the pile bottom, and the prestressed anchor cable is set at the pile foot to ensure that the stability of the pile meets the requirements. In order to protect the reserved rock shoulder, φ168 mm@1200 mm steel pipe piles are set on the foundation pit rock wall. The steel pipe pile is embedded 1.0 m below the foundation pit bottom, and M30 cement mortar is poured into the pipe. The typical supporting structure section of deep foundation pit in binary strata of upper soil and lower rock of Ningxia road station of Qingdao Metro is shown in Figure 14.

According to the statistical analysis of the field measured data of 11 surface settlement monitoring points (No. DC12∼dc18, DC24∼dc27, and dc36), 6 pile top vertical displacement monitoring points (No. QC06∼QC11), 6 pile top horizontal displacement monitoring points (No. QW06∼QW11), and 2 pile horizontal displacement monitoring points (No. CX02, CX04) are arranged by the third-party monitoring unit around the station foundation pit. The maximum surface settlement around the foundation pit is 19.25 mm, the minimum value is 4.16 mm, and the average value is 12.83 mm. The maximum vertical displacement of pile top is 8.72 mm, the minimum is 1.10 mm, and the average is 3.47 mm. The maximum horizontal displacement of pile top is 5.72 mm, the minimum is 0.24 mm, and the average is 2.27 mm. The maximum value of pile body is 18.63 mm, the minimum value is 11.25 mm, and the average value is 14.78 mm, as shown in Table 11. The layout of monitoring points is shown in Figure 15. Typical time history curves of monitoring data are shown in Figure 16.

It can be seen from Table 11 that the ground settlement around the foundation pit and the deformation of retaining pile support structure caused by the foundation pit excavation of Ningxia Road Station of Qingdao metro are within the allowable range. It shows that the stability of the foundation pit is good, and the rationality of the conclusion of the foundation pit vertical sidewall self-stable height is verified.

7. Conclusions

(1)Based on comparative analysis of the stability characteristics for deep foundation pit in binary strata of upper soil and lower rock under multiple working conditions, the potential fracture surface of deep foundation pit and the evolution law of corresponding stability FOS are revealed under different upper soil stratum thickness (Hs) and excavation depth (H). A new idea that the vertical rock sidewall height and the vertical rock sidewall height are used as two independent evaluation indexes, respectively, for deep foundation pit stability in the binary strata of upper soil and lower rock is put forward.(2)The concepts of critical self-stable height of vertical soil sidewall (Hs0) and critical self-stable height of vertical rock sidewall (Hr0) for deep foundation pit in the binary strata of upper soil and lower rock are proposed. The distribution characteristics and variation law of Hs0 and Hr0 under different Hs and different H are revealed, respectively. The spatial distribution map of the self-stabilizing height for deep foundation pit vertical sidewall in upper soil and lower rock binary stratum is constructed, and the mathematical fitting equation between Hr0 and Hs is obtained.(3)Combined with the implementation effect of the deep foundation pit project of Ningxia Road Station for Qingdao Metro Line 3, the rationality of the conclusions is verified.

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.

Authors’ Contributions

Ziguang Zhang is in charge of the paper and completed the construction of thesis framework system and thesis writing. Jiesheng Zhang and Tao Xu completed the collection, sorting, and analysis of relevant project case data in the paper. Yihang Li and Guangyong Cao assisted in completing the related work of numerical calculation in the paper. YiKang Xu completed the related work of the language check and editable picture design in the paper.

Acknowledgments

This study was funded by Scientific Research Project of Anhui Jianzhu University (no. 2019QDZ24), Science and Technology Plan of Housing and Urban Cural Construction in Anhui Province (no. 202053), and Nantong Science and Technology Plan Project (no. JC2021169). The authors appreciatively acknowledge the financial support of the abovementioned agencies.