Abstract

In this study, a detailed analysis was conducted to evaluate the impacts of the deviatoric stress component and spherical stress component on the stability of surrounding rocks in the roadway via the theoretical analysis and calculation and numerical simulation. Based on the analysis, the distribution laws guiding the main stress differences, plastic zone, convergence of surrounding rocks, and third invariant of stress under various conditions (such as equal spherical stress and unequal deviatoric stress and equal deviatoric stress and unequal spherical stress) were developed, providing an optimization scheme for roadway support misunderstanding under the conditions of high spherical stress field and high deviator stress field. The study further reveals that under the circumstance of the constant spherical stress, the greater the deviatoric stress, the plastic zone range of the surrounding rock of the roadway, the range of tensile deformation of the surrounding rock, the amount of convergence of the surrounding rock, the probability of separation of the roof and floor of the roadway, and the principal stress difference and the main stress, the greater the concentration range of the maximum stress difference is, and the maximum principal stress difference is mainly concentrated in the roof and floor rocks of the roadway, and the greater the deviatoric stress, the greater the probability that the roof and floor rocks of the roadway will be separated, and the maximum principal stress difference is mainly concentrated in the roof and floor rocks of the roadway, the greater the deviator stress, the greater the concentration range of the maximum value of the principal stress difference and the principal stress difference; when the deviator stress is constant, the range of the plastic zone and the maximum principal stress difference concentration range of the surrounding rock of the roadway decrease with the increase of the ball stress, and the principal stress difference, the amount of convergence of the surrounding rock, and the range of tensile deformation increase with the increase of the ball stress. The maximum principal stress difference is mainly concentrated in the roof and floor rocks of the roadway. The principal stress difference increases with the increase of the spherical stress, and the maximum concentration range of the principal stress difference decreases with the increase of the spherical stress. After the method proposed in this paper optimizes the actual roadway support on site, the surrounding rock deformation of the roadway is small and the control is relatively ideal, which basically meets the engineering needs.

1. Introduction

Ground stress has been considered as the fundamental force causing deformations and damages in various underground excavation projects [1, 2] and one of the important bases for designing the support and protection system of the underground projects. Ground stress has been divided into two categories including self-weight stress and horizontal tectonic stress. Therefore, scholars domestically and abroad have invested substantial time and energy in researching the impacts of self-weight stress [36] and horizontal tectonic stress [710] on the stability of the roadway surrounding rocks.

Based on plastic mechanics, the stress of surrounding rocks can be categorized into spherical stress tensor and deviatoric stress tensor, with the former one deciding the shape changes of the rocks and the latter one dominating the volume changes [11]. Substantial studies have been conducted to explore the distribution laws of spherical stress of roadway surrounding rocks and plastic zone. Ma et al. researched the deviatoric stress field and the distribution law of plastic zone when the roadway surrounding rocks are exposed to uneven stresses [12]. Yu et al. investigated the impacts of deviatoric stress on the plastic zone distribution of roadway surrounding rocks, resulting in the identification of an instability model of the roadway under various lateral pressure coefficients [13]. Xie et al. studied the variance law of deviatoric stress of deep roadway surrounding rocks and proposed asymmetric control technology of surrounding rocks based on various sections [14]. He et al. focused their studies on the structural stress of deep roadway surrounding rocks at high elevation including the damage and the distribution of deviatoric stress, along with corresponding controlling measures [15]. Xu et al. took the superelevation section of a coal mine track in Shanxi as the research object, used UDEC to simulate the deviatoric stress distribution of the surrounding rock of the roadway at different roadway heights, compared the degree of deviatoric stress changes between the roof and floor and the two sides of the roadway, and finally proposed targeted support technology [16]. Luo et al. thoroughly studied the influence of the intermediate principal stress and the rheological properties of the surrounding rock on the displacement of the surrounding rock of the roadway and the plastic zone of the surrounding rock and finally found that ignoring the rheological properties of the rock would overestimate the lithology of the surrounding rock. In control roadway deformation and plastic zone expansion [17], Zhang et al. used the D-P yield criterion to calculate the analytical solutions for the elastoplasticity, plastic zone radius, and displacement of the surrounding rock under bidirectional isobaric conditions and discovered the importance of the intermediate principal stress to the stress distribution of the surrounding rock [18]. According to the D-P yield criterion and the nonassociated flow rule, Chen et al. derived a closed analytical solution for the stress, deformation, and radius of the plastic zone of the surrounding rock of a deep circular roadway under hydrostatic pressure [19].

Scholars have conducted sufficient studies on the deviatoric stress distribution of roadway surrounding rocks and the distribution law of the plastic zone after the roadway excavation. However, limited studies have been performed on the impacts of both spherical stress tensor and deviatoric stress tensor before the roadway excavation on the roadway stability after excavation. Therefore, in this study, a detailed analysis was conducted to evaluate the impacts of deviatoric stress component and spherical stress component on the stability of surrounding rocks in roadway via the theoretical analysis and calculation, numerical simulation, and two sets of loading tests.

2. Loading Plan and Numerical Model

2.1. Loading Plan and the Analysis of Principal Stress Difference of Roadway Surrounding Rock

A plan including two sets of loading tests was proposed to study the impacts of the deviatoric stress component and spherical stress component on the stability of surrounding rocks in the roadway. The detailed loading plan is provided below.

A loading test was performed at a constant spherical stress component () and various deviatoric stress components (, 7.2 MPa, 12 MPa, and 16.8 MPa). In other words, during this loading test, constant spherical stress was imposed on the model with different deviatoric stresses.

A loading test was performed at a constant deviatoric stress component () and various spherical stress components (, 12 MPa, 15 MPa, and 18 MPa). In other words, during this loading test, constant deviatoric stress was imposed on the model with different spherical stresses.

Equation (3) can be obtained based on equations (1) and (2) listed above.

where refers to spherical stress, indicates deviatoric stress, and is the loading angle.

A stress space diagram was obtained from equation (3) and illustrated in Figure 1.

The strain increment of the deep surrounding rocks under the plastic strain state demonstrated pure shear deformation with the maximum shear stress dominating the generation of the stratum plastic zone and development [20, 21]. The principal stress difference can reflect the shear stress distribution and the damage in the surrounding rocks [22, 23]. The principal stress difference can be illustrated in

Based on the elasticity theory, the stress of the circular hole in the two-way stress infinite plate is shown in

where refers to the vertical stress in MPa, indicates the horizontal stress in MPa, is the radius of the circular hole in m, and represents the polar angle.

Under the state of the plane strain, the principal stress can be calculated following

where refers to the maximum principal stress in MPa and indicates the minimum principal stress in MPa. Equations (6) and (7) can be converted into equations (8) and (9) under polar coordinates.

Considering the roadway as the plane strain allows the overlook of the principal stress, based on equations (3) and (5), equation (10) can be obtained.

Under various deviatoric stresses, -10° load was selected for load while -20° load was selected under various spherical stresses. Based on equation (4) to equation (10), the surrounding principal stress differences of the circular hole under various conditions were obtained and presented in Figure 2.

Based on Figure 2, the following results can be obtained. (1) Under constant spherical stress, as the deviatoric stress increases, the principal stress difference on the vertical direction tends to increase accordingly, with the principal stress difference of shallow surrounding rocks greater than that of deep surrounding rocks. Horizontally, the principal stress difference of shallow surrounding rocks decreases as the deviatoric stress increases while the principal stress difference of deep surrounding rocks increases. Overall, the vertical principal stress difference is greater than the horizontal one. (2) Under constant deviatoric stress, as the spherical stress increases, the vertical principal stress difference increases, with the principal stress difference of shallow surrounding rocks greater than that of deep surrounding rocks. Horizontally, the principal stress difference of shallow surrounding rocks increases as the spherical stress increases while the principal stress difference of deep surrounding rocks decreases. Overall, the vertical principal stress difference is greater than the horizontal one.

2.2. Numerical Model and Parameters

A numerical simulation was performed to analyze the impacts of deviatoric stress and spherical stress on the stability of surrounding rocks in the roadway. In the simulation model, the cross section of the roadway was designed in a semicircular arch, with a height of 3500 mm and a width of 4000 mm. The calculation range was set at . The sides and the bottom of the model limited the horizontal and vertical displacements. The rock mass mechanical parameters of each rock stratum are shown in Table 1. Stresses were loaded to the model following the loading plan listed in Section 2.1.

3. The Test Results under a Constant Spherical Stress

3.1. The Distribution Law of the Principal Stress Difference of the Roadway Surrounding Rocks

With FLAC3D, a cloud diagram shown in Figure 3 was drafted to demonstrate the principal stress difference of the roadway surrounding rocks under constant spherical stress and various deviatoric stresses.

From Figure 3, the following results can be obtained. (1) When , the maximum principal stress difference concentrated within 2 to 5 m range of the roadway surrounding rocks, with maximum principal stress difference observed in the lower stratum of the surrounding rocks. As the deviatoric stress increased, the maximum principal stress difference continued to increase, which is consistent with the trend demonstrated in Figure 2(a), suggesting a higher shear resistance needed in the anchor bolt (cable) installed in the roadway roof in a high deviatoric stress field. (2) As the deviatoric stress increased, the concentration zone of the maximum principal stress in the roadway surrounding rock continued to expand and migrate to the deep layers of rocks, which suggests that longer anchor cables are needed to keep anchor cables fixed to the roof for the purpose of support protection. (3) The deviatoric stress had limited impacts on the principal stress difference of the roadway sides.

3.2. The Variance Law of the Plastic Zone of the Roadway Surrounding Rocks

The distribution of the plastic zone contributes to the identification of the supported depth. The distribution of plastic zone under various deviatoric stresses is demonstrated in Figure 4. The damage depth of the surrounding rock corresponding to the deviatoric stress is shown in Figure 5.

Based on Figures 4 and 5, the following results can be drawn. (1) As the deviatoric stress increases, the plastic zone in the roadway roof grows large, suggesting that longer supporting plates are needed to enhance the stability of the roadway in a high deviatoric stress field, which is consistent with Figure 3. (2) When , the increase of deviatoric stress tends to have limited impacts on the damage of roadway sides. When , the damage of roadway sides increases drastically. Taking Figure 3 into consideration, the maximum principal stress difference concentrates in the lower stratum of the roadway roof. In addition, higher deviatoric stress leads to a larger maximum principal stress difference concentration zone. Therefore, the expansion of the maximum principal stress difference concentration zone is the main cause of the drastic increases in the damage of roadway sides, indicating that proper support and protection are critical to prevent the roadway instability caused by the damage of roadway sides under high deviatoric stress.

3.3. The Deformation Law of the Surrounding Rocks

The deformation degree of the surrounding rocks can reflect the need for ductility of the support system. Under various deviatoric stresses, the deformation of the roadway surrounding rocks from the surface to the depth is demonstrated in Figure 6.

According to Figure 6, (1) the displacement of surrounding rock from shallow to deep decreases in a “logarithmic” pattern, until reaching stability. (2) When , negative displacement was observed at the 3 m of roadway floor in-depth and 2 m of roadway roof, which suggests that in the floor, the deep rock stratum submerges while the shallow ones rise. On the other hand, the deep rock stratum rises while the shallow ones submerge in the roof. The observations suggest that higher deviatoric stress tends to increase the probability of separation in the roadway roof and floor, resulting in a higher possibility of roof collapse and bottom drum. (3) The surface displacement of the roadway can reflect the maximum deformation of the roadway surrounding rocks, the damage of the shallow surrounding rocks, and separation. According to Figure 6, as the deviatoric stress increases, the surface displacement of the roadway increases in a linear pattern. When increased above 2 MPa, the surface displacement of the roadway decreases, which suggests that when is higher than 12 MPa, the deep stratum in the roadway roof starts to rise substantially, resulting in a reduced convergence of the surrounding rocks.

3.4. The Distribution Law of the Third Invariant of Deviatoric Stress in Surrounding Rocks

The third invariant of deviatoric stress can be used to identify the deformation type of the surrounding rocks. When , the deformation is categorized as compressive deformation. When , the deformation is considered as planar deformation. When , the deformation is rated as tensile deformation [8]. The third invariant of deviatoric stress can be a comprehensive index of combining the maximum principal stress, minimum principal stress, and medium principal stress, indicating the deformation type of the surrounding rocks. The third invariant of deviatoric stress can be expressed in

A cloud diagram shown in Figure 7 was obtained to demonstrate the third invariant under various deviatoric stresses. According to Figure 7, when , a tensile stress zone was observed within 0.3 to 2.5 m range of the roadway surrounding rocks. As the deviatoric stress increases, the tensile stress zone expands, resulting in a larger tensile stress concentration zone in roadway sides and bottom than the roof. Due to the poor tensile strength of the rocks, high tensile stress can jeopardize the stability of the roadway surrounding rocks. Therefore, a high pretorque value should be adopted to the supporting system of the roadway surrounding rocks in a high deviatoric stress field to enhance the load capacity of the surrounding the rocks in roadway.

4. The Test Results under a Constant Deviatoric Stress

4.1. The Distribution Law of the Principal Stress Difference of the Roadway Surrounding Rocks

A cloud diagram shown in Figure 8 was drafted to demonstrate the principal stress difference of the roadway surrounding rocks under constant deviatoric stress and various spherical stresses.

According to Figure 8, (1) as the spherical stress increases, the maximum principal stress concentration range in the deep surrounding rocks tends to decrease, suggesting that in a high spherical stress field, the layout and installation of the anchor bolts and cables should avoid the principal stress concentration zone for better supporting and stability. (2) As the spherical stress increases, the principal stress difference increases in the surrounding rocks of the roadway bottom, which is consistent with Figure 2(b), demanding a higher shear resistance in the anchor bolts and cables. (3) The spherical stress exerts some influences on the principal stress difference of the shallow surrounding rocks in the roadway sides. However, such influence is limited to the deep surrounding rocks, which is also demonstrated in Figure 2(b), resulting in a lower shear resistance requirement to the supporting and protection system of the roadway sides.

4.2. The Deformation Law of the Surrounding Rocks

The distribution of plastic zone under various spherical stresses is demonstrated in Figure 9. The damage depth of the surrounding rock corresponding to the spherical stress is shown in Figure 10.

According to Figures 9 and 10, (1) as the spherical stress increases, the plastic zone of the roadway surrounding rocks tends to decrease due to the increased surrounding pressure on the roadway surrounding rocks under higher spherical stress. In other words, a higher surrounding pressure leads to a higher strength of the surrounding rock and a smaller plastic zone. (2) Under various spherical stresses, the plastic zone of the roadway surrounding rock demonstrates a butterfly shape. The spherical stress only impacts the range of the plastic zone instead of the shape. Comparing with Figure 5 has revealed that the deviatoric stress tends to exert a greater impact on the surrounding rocks of the roadway. (3) The damage depth of the roadway surrounding rocks demonstrates a higher sensitivity to the spherical stress while the damage of the roadway sides is less sensitive.

4.3. The Deformation Law of the Surrounding Rocks

The deformation of the roadway surrounding rocks under various spherical stresses is demonstrated in Figure 11.

According to Figure 11, (1) a threshold depth has been observed in rating the impact of spherical stress on the roadway surrounding rocks. When the spherical stress is higher than the threshold depth, the deformation is limited with a low degree of dispersion associated with the surrounding rock deformation. When the spherical stress is lower than the threshold depth, the deformation is severe with a high degree of dispersion associated with the surrounding rock deformation. (2) Under low spherical stress, the convergences in the deep and shallow surrounding rocks of roadway are negative and positive, respectively. As the spherical stress increases, the convergences from the deep to the shallow surrounding rocks become positive, suggesting that the increase of the spherical stress can lower the probability of the stratum separation. (3) As the spherical stress increases, the convergence value of the roadway surface surrounding rocks follows a linear growth.

4.4. The Distribution Law of the Third Invariant of Deviatoric Stress in Surrounding Rocks

A cloud shown in 12 was developed to demonstrate the distribution of the third invariant of deviatoric stress in the roadway surrounding rocks under various spherical stresses.

According to Figure 12, when , the roadway surrounding rocks are under pressure as a whole. When , a tensile stress zone was spotted around the roadway shallow rocks. As the spherical stress increased, multiple tensile stress zones were observed in the roadway sides and surrounding rocks of the roof. In other words, as the spherical stress increases, the tensile stress zone expands in the roadway surrounding rocks, reducing the loading capacity of the surrounding rocks. Therefore, in a high spherical stress environment, a high pretorque value should be adopted to the supporting system of the roadway surrounding rocks to enhance the load capacity of the surrounding rocks in roadway.

5. The Instability Analysis of the Roadway Surrounding Rocks

5.1. The Analysis of the Impact of Deviatoric Stress on the Instability of the Roadway Surrounding Rocks

Taking equation (2) into consideration, the actual deviatoric stress of the roadway surrounding rocks depends on the ratio of three principal stresses instead of the individual value of three principal stresses, which is known as the lateral pressure coefficient. In other words, the deviatoric stress of the roadway surrounding rocks in a deep roadway under high ground stress does not necessarily lead to high deviatoric stress. The deviatoric stress of the roadway surrounding rocks in a shallow roadway under high ground stress does not necessarily lead to low deviatoric stress either. According to the analysis in Section 3, the deviatoric stress tends to play a critical role in influencing the stability of the roadway surrounding rocks. In fields, many shallow roadways under low stresses tend to experience roof collapse and ground drums, such as the roadway located in the mines of the Shendong area. Therefore, according to the analysis in Section 3 along with the field investigation, some common misunderstandings have been identified to the support and protection system in the shallow surrounding rocks under low stress.

Firstly, according to Figures 4 and 5, the deviatoric stress tends to have a significant impact on the damage range and depth of the roadway surrounding rocks. However, regarding the shallow roadway under high deviatoric stress, the impacts of the stress on the stability of the roadway surrounding rocks are often overlooked during the design of the support and protection system, assuming that due to the shallow location of the roadway, no additional support system is needed, leading to a design of support and protection system shown in Figure 13.

As a result, the support and protection system lacks solid anchor points, accompanied by anchor bolts and cables following the movements of the surrounding rocks.

According to Figure 6 and the analysis in Section 3.3, as the deviatoric stress increases, the convergence of the roadway surrounding rocks increases accordingly, along with potential separation in the roadway roof. Without a strong support system such as a low elongation in anchor bolts and cables, a low support density, and a low strength in the metal protection net, the protection system could risk a total failure, resulting in severe deformation in the roadway sides including roof collapse and floor drums.

According to Figure 7, as the deviatoric stress increases, the tensile deformation continues to expand in the roadway surrounding rocks. Due to the poor tensile strength of the rocks, the expansion of the tensile deformation contributes to the increasing of the plastic zone and decreasing of the roadway stability, which is also reflected in Figure 4. With a low rigid support and protection system, such as a low pretorque value and poor timing for the installation of the support and protection system, a large tensile deformation zone can occur in the surrounding rocks, resulting in the continuous expansion of the plastic zone and severe deformation to the surrounding rocks in a high deviatoric stress environment.

5.2. The Analysis of the Impact of Spherical Stress on the Instability of the Roadway Surrounding Rocks

According to equation (1), the spherical stress value mainly depends on three principal stresses. It is commonly known that a greater depth leads to higher spherical stress. Based on Section 4, the spherical stress exerts a significant impact on the damage range of the roadway surrounding rocks, bearing capacity, and the deformation of the surrounding rocks. Therefore, according to the analysis in Section 4 and Figure 14 along with the field investigation, some common misunderstandings have been identified to the support and protection system in the shallow surrounding rocks under low stress including increasing the supported depth of the anchor bolts and cables blindingly and lacking specificity in the design of the support system.

In the fields, engineers and designers are misled to believe that a deeper roadway tends to experience more severe damage, resulting in increasing the supported depth blindingly. Based on the analysis in Sections 3 and 4, the damage range of the roadway surrounding rocks increases along with the increase of deviatoric stress and decrease of the spherical stress. As a result, regarding a roadway in an environment featured with high spherical stress and low deviatoric stress, the damage of the surrounding rocks should be limited, eliminating the need of increasing the supported depth. However, the increase of spherical stress tends to lead to the increase in the deformation of the roadway surrounding rocks, principal stress difference, and tensile stress zone, demanding a support system featured with high strength, shear resistance, and contractibility. Therefore, for the deep roadway in an environment featured with high spherical stress and low deviatoric stress, the support system should be improved by increasing the strength, shear resistance, and contractibility instead of the supported depth.

Due to insufficient theoretic guidance, the understandings of the critical factors of the stability of roadway surrounding rocks are limited. For instance, as demonstrated in Figure 9, the damage of the surrounding rocks concentrated in the roadway shoulders and bottom corners instead of the roadway sides. Therefore, the roadway shoulders and bottom corners should be the main protection targets.

6. Engineering Validation

6.1. Project Overview and Maintenance Measures

A mine located in Inner Mongolia was selected for the purpose of engineering validation. The mining focuses on the #3 coal seam with an average thickness of 5 m, categorized as a horizontal coal seam with an average depth of 275 m. Currently, 310 working face is about 300 m away from the stopping mining line. No. 311 working face is featured with a through open cut. The transportation roadway of the No. 311 working face is 230 m away from the No. 310 goaf. Therefore, the No. 311 transportation roadway is not influenced by the mining process. Due to the depth of the #3 coal seam, the limited ground force, and relatively simple geological structure, supports were only added to the transportation roadway roof and sides, as demonstrated in Figure 15. The support system includes steel anchor rods in , arranged in a triangular layout with a row distance of . Roadway sides were protected with steel anchor rods in , arranged in a flower layout with a row distance of .

After the installation of the support system demonstrated in Figure 15, severe damage was spotted in the roof of the No. 311 transport roadway and roadway sides, with some anchor rod failures, such as migrating with the surrounding rocks, tears in the metal net, and more. Several deep damages concentrated on the roof along with maximum deformation of 200 mm in the bottom, risking roof collapse.

Considering the limited ground force at the location of No. 311 transport roadway and no influence from the mining process, No. 311 transport roadway damage was mainly caused by high deviatoric stress. Following the analysis in Section 3, the deformation of roadway surrounding rocks in a high deviatoric stress field is featured with large damage scale, large roadway surface displacement, potential stratum displacement in the roof, and large tensile stress zone, which were consistent with the deformation of No. 311 transport roadway. In addition, some improper designs were identified in the support plan of No. 311 transport roadway. Therefore, based on the research results, some adjustments were made to the original design plan shown in Figure 16.

Anchor cables can link the load structure of the anchor rods to the deep surrounding rocks, resulting in an interconnected and overlapped network structure in the effective stress zone, resulting in decreasing the tensile stress zone of the surrounding rocks and increasing the stability of the load-bearing structure [24]. Following this concept, the support design of the No. 311 transport roadway was modified as follows. In order to minimize the separation in the roof surrounding rocks and reduce the tensile zone, two additional high-strength and high-elongation prestressed steel strands (dimension of ) were added with a preload no less than120 kN. The distance between two strands was designed at . Meanwhile, in order to control the roadway deformation and tensile deformation and improve the support efficiency of the anchor rods, two additional high-strength prestressed steel strands were added with a row distance of . Besides, an inclined anchor rod was installed at the roadway sides, 15° in contrast with the horizontal direction for the purpose of improving the stability of roadway bottom corners.

6.2. The Analysis of the Modified Support System

Multiple monitoring and detection points were set up to keep track of the roadway surface displacements. The collected data indicated that with the modified support system, the deformation rate of the surrounding rocks was initially maintained at 8~20 mm/d with a higher deformation rate identified in the roadway sides and roof. A total of 4 changes were observed to the roadway sides and roof. The deformation rate demonstrated a downward trend. At day 37, the deformation reached stability, indicating an effective support and protection to the transport roadway.

7. Conclusion

For the principal stress difference distribution of the circular roadway surrounding rocks under various spherical stresses and deviatoric stresses via theoretical analysis, the calculation results were consistent with the numerical simulation. (1)Under constant spherical stress, the plastic zone, the probability of stratum separation of the roof, the convergence of the surrounding rocks, and the tensile zone increase along with the deviatoric stress. The principal stress difference of the surrounding rocks and the concentration zone of the maximum principal stress difference tends to increase as the deviatoric stress increases(2)Under constant deviatoric stress, as the spherical stress increases, the plastic zone decreases while the convergence of the surrounding rocks and the tensile zone increases. In addition, higher spherical stress increases the principal stress difference of the surrounding rocks. Meanwhile, the concentration zone of the maximum principal stress difference tends to decrease as the spherical stress increases(3)Some common misunderstandings about the roadway support in a high spherical stress field and in a high deviatoric stress field were proposed, which was further validated in the field. The field test suggests that the modified support system following the research results can achieve effective protection and support

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was financially supported by the central government guides local science and technology development fund projects (216Z5401G), Guizhou Province High-Level Innovative Talent Training Program Funding Project ([2019]5675), and Scientific Research Projects of Higher Education Institutions in Hebei Province (Z2020124).