Orthogonal Test on the True Triaxial Mechanical Properties of Frozen Calcareous Clay and Analysis of Inﬂuencing Factors

: In the Huainan and Huaibei mining areas, a layer of calcareous clay is buried deep in the surface soil layer (at approximately 400 m). This layer is in a high-stress state and is prone to freezing pipe fractures in the freezing method. To obtain the true triaxial mechanical properties of this clay in its frozen state, this study conducted a cross test (L 16 (4 5 )) to explore the change law of the strength of frozen calcareous clay under the inﬂuence of multiple factors. The results showed that the true triaxial stress–strain curve of frozen calcareous clay was divided into three stages: the strain within 0.5% showed linear elasticity. Under compressive stress, ice crystals and their cements were damaged or melted and shrank. At approximately 5%, they showed plastic hardening. The soil particles and ice crystals in the frozen soil recombined and became denser, resulting in irreversible deformation. As the compression progressed, cracks bred and swelled. The failure stage was manifested as strain hardening due to the test loading conditions. As the deformation increased, the stress also slightly increased. The consistent strength-inﬂuencing factors could be obtained through range and hierarchy analyses. The primary and secondary order of inﬂuence of σ 1 was the conﬁning pressure σ 3 , water content ω , temperature T , Bishop parameter b , and salt content ψ . The inﬂuence weight of each factor was quantitatively calculated. In the signiﬁcance analysis, when the interaction was not considered, the effects of the pressure and moisture content on the strength were always signiﬁcant. The effect of temperature was signiﬁcant only when the signiﬁcance level Ω > 0.05. The salt content and b value had no signiﬁcant inﬂuence on the strength, and the signiﬁcance of each factor followed the order of the results of the range analysis method and analytic hierarchy process; when considering the interaction, the interaction factors had different effects on the strength. When Ω > 0.01, the inﬂuence of factor A (temperature T ) × B (water content ω ) on the strength showed signiﬁcance, even exceeding that of temperature. This demonstrated that when studying the strength characteristics of frozen soil, it is necessary to comprehensively consider the various factors and their interaction to more accurately characterize the mechanical behavior of frozen solids.


Introduction
The geotechnical strength index is an important survey parameter for the development and utilization of underground space.Under normal circumstances, conventional soil is considered a two-(saturated) or three-phase system, and laboratory testing and research have been conducted on such soils.For shallow rock and soil, the unconfined uniaxial compressive strength under different factors have been considered [1][2][3].Such soils meet the engineering requirements; however, when the deep rock and soil are in a state of high stress, the contribution of the confining pressure to the strength cannot be ignored.To simulate a real stress environment, it was necessary to conduct false or true triaxial indoor tests [4][5][6].The factors influencing the strength in such tests include the water content (or saturation), soil properties, loading rate, confining pressure, and dry density.The strength of conventional soil is generally low, and the influence of various factors on the strength was not significant.Frozen soil is a four-phase system.Artificial ground freezing curtains are used as temporary support measures in major underground projects such as coal mine shaft excavation, subway tunnels, and special underground spaces.The influence of various factors on the properties of frozen structures, such as thermal conductivity [7,8] and strength, is greatly enhanced and more complex.Cai [9] conducted an unconfined compressive strength test on frozen clay under different water and Na 2 SO 4 contents.The test results showed that the deformation of frozen soil was significantly affected by the change in the water and salt contents.Its strength increased to the maximum value first and then decreased with the increase of water content and salt content.Through a series of triaxial tests, Niu [10] found that the deformation and strength of frozen soil were significantly affected by temperature, moisture, and pressure.When the moisture content was low, and the confining pressure increased, a silty clay transitioned from strain softening to strain hardening.With the increase in the water content, the initial tangent modulus gradually changed from a linear distribution to a parabolic one with the confining pressure.Lai [11] and Du [12] analyzed the change characteristics of sand strength with the moisture content at different temperatures and the dependence of sand strength on the temperature under the same moisture content; the influence of moisture content on the strength of frozen soil did not purely depend on the moisture content.For a certain type of soil, there is a fixed moisture content (optimal moisture content).When it is greater than this value, the strength of frozen soil decreases with the increase in the moisture content [9,[13][14][15][16][17].The temperature, water content, salt content, and confining pressure significantly affect the strength of frozen soil.In terms of true triaxiality, the influence of Bishop parameter b on the strength of frozen soil, particularly the second principal stress, cannot be ignored.Hence, scholars have conducted true triaxial tests that are more representative of the real stratum stress state.For example, Dai [18] discussed the effects of the water content, consolidated confining pressure, and intermediate principal stress ratio on the strength of branched clay and reached a reliable conclusion.As reported by Chen [19], when the consolidated confining pressure was less than the structural strength of the structural loess, the stress-strain curve was a softening type; otherwise, it was a hardening type.Ma [20] and Zhang [21] conducted a true triaxial test on frozen sand under the influence of Bishop parameter b and accordingly performed an analysis.
Many factors affect the strength of frozen soil, not just the ones mentioned in the above studies.Most scholars considered only two or three of the many factors, and research on the interaction effects of the various factors is lacking.Hence, in this study, we considered a layer of calcareous clay in the deep surface of the Huainan and Huaibei mining area (buried at approximately 400 m) as the research object.This layer is in a high-stress state and highly prone to freezing pipe fractures in the freezing method.To obtain the true triaxial mechanical properties of this clay in its frozen state, we conducted true triaxial orthogonal tests on calcareous clay under the influence of five factors: temperature, water content, salt content, confining pressure, and Bishop parameter b.The primary and secondary order of influence of the various factors on the strength of the frozen soil was analyzed, and the weight of the influence of each factor on the peak strength was quantitatively calculated through the analytic hierarchy process.Finally, the interaction and non-interaction between the factors were considered, and statistical theory was applied to elaborate the influence of each factor.The degree of significance provided a certain reference value for the design and construction of underground structures in frozen soil.

Orthogonal Experimental Design Scheme
Frozen soil is used as a curtain for coal mine shaft excavation, and its strength is an important design index.Studies have shown that many factors affect the strength of frozen soil, including water and salt content, confining pressure, temperature, and indoor test loading rate.In this study, we used a calcareous clay exhibiting strong plasticity as the test object.Table 1 presents the basic physical properties of this clay.Considering the temperature, water and salt content, confining pressure, and intermediate principal stress coefficient b, a true triaxial test was conducted on the frozen calcareous clay.Because of the many factors considered and the large number of tests required, an orthogonal test method was used to explore the influence of the various factors on the strength of this frozen soil.The orthogonal test design is presented in Table 2.A five-factor four-level (L 16 (4 5 )) orthogonal test was conducted.The corresponding levels of each factor were:   [23]).

Sample preparation:
The dimensions of the test block sample were 100 mm × 100 mm × 100 mm.After the undisturbed soil sample was dried, ground, and sieved, it was sprayed with different concentrations of CaCl 2 solution in layers, mixed evenly, placed in a Ziplock bag, and allowed to stand for 24 h.The sample was then controlled using the dry density and volume.The mixed soil samples of corresponding quality were packed in a grinding tool in layers and pressed to prepare standard samples for use, as shown in Figure 1.Before the test, the standard sample was frozen at the required low temperature for more than 48 h, and a true triaxial test was conducted.Loading mode: For the test, we adopted a ZSZ-2000 microcomputer-controlled true triaxial frozen soil test platform.Figure 2 shows its loading device.The frozen soil triaxial test platform could realize independent sample loading in the three principal stress directions.The sample could be tested along different stress paths in the 3D stress space.The sample was loaded along the three principal stress directions using rigid plates.As shown Loading mode: For the test, we adopted a ZSZ-2000 microcomputer-controlled true triaxial frozen soil test platform.Figure 2 shows its loading device.The frozen soil triaxial test platform could realize independent sample loading in the three principal stress directions.The sample could be tested along different stress paths in the 3D stress space.The sample was loaded along the three principal stress directions using rigid plates.As shown in the three views in Figure 3, the F6 direction was the fixed end, which gave the reaction force in the F2 direction.Loading in the five directions, F1 to F5 was done independently by the servo cylinder through the force transmission rod.F3 and F5 were the loading directions of the first principal stress; F1 and F4 were the loading directions of the second principal stress; F2 and F6 were the third principal stress loading directions.Six steel plates were placed in a staggered overlap manner to realize the compression process of the sample.This test adopted a loading method with a constant principal stress ratio b in the control, where b = (σ 2 − σ 3 )/(σ 1 − σ 3 ).The parameter σ 3 was a constant.By shifting the term and taking the derivative of this formula, we find that ξ(σ 2 ) = bξ(σ 1 ), and ξ were the loading rate in the direction of the first principal stress, and ξ = 25N/s.The specific loading process of the test is shown in Figure 1.After the sample reached the required temperature for the test and stabilized for 30 min, a three-way synchronous pressure σ 3 was applied to place the specimen in a state of equal stress in the three directions.The specimen was then stabilized and consolidated until its deformation was unchanged.Different loading rates were subsequently set in the σ 1 and σ 2 directions.The ratio of the two loading rates was equal to b to realize the true triaxial test process while keeping b unchanged.The test failure threshold was when the deformation in the direction of the major principal strain reached 15 mm, that is ε 1 = 15%, or when the test monitoring curve shows that the sample was evidently damaged, that is, when the strain in the first principal strain-first principal stress curve grew almost vertically, the corresponding strength was the failure strength.Loading mode: For the test, we adopted a ZSZ-2000 microcomputer-controlled true triaxial frozen soil test platform.Figure 2 shows its loading device.The frozen soil triaxial test platform could realize independent sample loading in the three principal stress directions.The sample could be tested along different stress paths in the 3D stress space.The sample was loaded along the three principal stress directions using rigid plates.As shown in the three views in Figure 3, the F6 direction was the fixed end, which gave the reaction force in the F2 direction.Loading in the five directions, F1 to F5 was done independently by the servo cylinder through the force transmission rod.F3 and F5 were the loading directions of the first principal stress; F1 and F4 were the loading directions of the second principal stress; F2 and F6 were the third principal stress loading directions.Six steel plates were placed in a staggered overlap manner to realize the compression process of the sample.This test adopted a loading method with a constant principal stress ratio b in the control, where ( ) ( ) . The parameter σ3 was a constant.By shifting the term and taking the derivative of this formula, we find that ( ) ( ) and ξ were the loading rate in the direction of the first principal stress, and 25 N s ξ = . The specific loading process of the test is shown in Figure 1.After the sample reached the required temperature for the test and stabilized for 30 min, a three-way synchronous pressure σ3 was applied to place the specimen in a state of equal stress in the three directions.The specimen was then stabilized and consolidated until its deformation was unchanged.Different loading rates were subsequently set in the σ1 and σ2 directions.The ratio of the two loading rates was equal to b to realize the true triaxial test process while keeping b unchanged.The test failure threshold was when the deformation in the direction of the major principal strain reached 15 mm, that is 1 15% ε = , or when the test monitoring curve shows that the sample was evidently damaged, that is, when the strain in the first principal strain-first principal stress curve grew almost vertically, the corresponding strength was the failure strength.

Test Results and Discussion
Considering that the test factors and loading environment of each sample of this true triaxial orthogonal test were different, there were many test curves; therefore, only some test equivalent stress-strain curves were selected for the analysis, as shown in Figure 4, where q was the von Mises stress, see Equation ( 1).The first principal strain ε1 was the sum of the strains in the F3 and F5 directions; the volume strain εv was the sum of the strains in directions F1 to F6, and the sign indicated that the pressure was positive.

Test Results and Discussion
Considering that the test factors and loading environment of each sample of this true triaxial orthogonal test were different, there were many test curves; therefore, only some test equivalent stress-strain curves were selected for the analysis, as shown in Figure 4, where q was the von Mises stress, see Equation (1).The first principal strain ε 1 was the sum of the strains in the F3 and F5 directions; the volume strain ε v was the sum of the strains in directions F1 to F6, and the sign indicated that the pressure was positive.
Figure 3. Three views of lapped rigid plates for frozen soil samples.

Test Results and Discussion
Considering that the test factors and loading environment of each sample of this true triaxial orthogonal test were different, there were many test curves; therefore, only some test equivalent stress-strain curves were selected for the analysis, as shown in Figure 4, where q was the von Mises stress, see Equation (1).The first principal strain ε1 was the sum of the strains in the F3 and F5 directions; the volume strain εv was the sum of the strains in directions F1 to F6, and the sign indicated that the pressure was positive.

(
) ( ) ( ) Figure 4 shows that the equivalent stress-strain relationship of the frozen calcareous clay was divided into three stages.The first was the linear elastic stage.Due to the high water content in the three selected samples, the samples were consolidated and frozen for a long duration.When the first principal strain was less than 0.5%, there was a certain instantaneous strength between the soil particles and cementation between the soil particles and ice crystals.Frozen soil samples exhibited good integrity and resistance to small deformations and rapid stress growth.The second was the plastic hardening stage.During this period, due to the increase in stress, the phenomenon of ice crystal pressure melting occurs.The soil and ice crystal particles, and pores reappeared in the frozen soil sample.The sample became denser, the bearing capacity was further improved, and plastic deformation from the first principal strain occurred, reaching approximately 5%.The third was the strain hardening failure stage.The occurrence of this stage depended on the test loading conditions.As mentioned above, a six-direction independent rigid loading mode was applied to control the load applied to the frozen soil samples.Frozen soil typically has certain residual stress after failure; if the frozen soil sample does not completely Figure 4 shows that the equivalent stress-strain relationship of the frozen calcareous clay was divided into three stages.The first was the linear elastic stage.Due to the high water content in the three selected samples, the samples were consolidated and frozen for a long duration.When the first principal strain was less than 0.5%, there was a certain instantaneous strength between the soil particles and cementation between the soil particles and ice crystals.Frozen soil samples exhibited good integrity and resistance to small deformations and rapid stress growth.The second was the plastic hardening stage.During this period, due to the increase in stress, the phenomenon of ice crystal pressure melting occurs.The soil and ice crystal particles, and pores reappeared in the frozen soil sample.The sample became denser, the bearing capacity was further improved, and plastic deformation from the first principal strain occurred, reaching approximately 5%.The third was the strain hardening failure stage.The occurrence of this stage depended on the test loading conditions.As mentioned above, a six-direction independent rigid loading mode was applied to control the load applied to the frozen soil samples.Frozen soil typically has certain residual stress after failure; if the frozen soil sample does not completely lose its bearing capacity or the test process was terminated, the force acting on the sample will continue to increase.As shown in Figure 4, the reason for the gradual growth in the deviator stress in the later stage was the rapid growth and expansion of the cracks in the sample under the action of the high stress applied (close to the ultimate bearing capacity).The samples expanded first and then shrank during the entire loading process.The expansion increased the gap between the rigid plates.When the specimen was placed between the rigid plates became incompressible, it overflowed from the gap between the plates, resulting in large deformation and unloading of the frozen soil specimens.This was also the reason for the evident shrinkage of the body.
Table 3 presents the results of the true triaxial orthogonal test conducted on the saline calcareous clay.The range analysis of the data listed in Table 3 was performed.Tables 4 and 5 list the analysis results.From a comparison of the R-values in Table 5, the primary and secondary order of influence affecting the peak strength σ 1 of frozen calcareous clay was as follows: pressure σ 3 , water content ω, temperature T, Bishop parameter b, and salt content ψ.The primary and secondary order affecting σ 2 is shown in Table 6: the Bishop parameter b, confining pressure σ 3 , water content ω, temperature T, and salt content ψ.Comparing the value of L, factor C (that is, the salt content ψ) had a certain effect on the peak strength and second principal stress of frozen calcareous clay compared with the other four factors, but the degree of influence was weaker, and the four levels show a trend of first increasing and then decreasing.The Bishop parameter reflected the magnitude of the second principal stress.Therefore, the influence of this factor on the second principal stress was the most significant.This was self-evident; the sample was more compact, the existence of an external force increased the friction on the sliding surface, the cohesive force increased, and the peak strength increased.Existing research results show that [24,25] temperature was the main factor affecting the strength of frozen soil; this was because the moisture in the sample freezes in a low-temperature environment.The lower the temperature, the greater the volume percentage of the ice crystal cement.The lower the unfrozen water content, which increased the strength of the cementation between the soil particles, the higher the strength.However, the presence of salt changes the concentration of the aqueous solution in the sample, thereby changing its freezing temperature, and the water was difficult to freeze.Therefore, the greater the concentration of the salt solution, the greater the decrease in the peak intensity.To systematically study the relationship between the peak stress σ 1 and the intermediate principal stress σ 2 of frozen saline calcareous clay and the confining pressure σ 3 , water content ω, temperature T, Bishop parameter b, and salt content ψ, according to the dimension analysis theory, an exponential function was used to compare the test data in Table 4.The calculation formula could be obtained by fitting in Origin as follows: To more intuitively reflect the significant effects of the confining pressure and water content on the peak stress and the second principal stress in the range analysis, the confining pressure σ 3 and water content ω were drawn for a temperature T = −10 • C and salt content ψ = 1%.Figures 5 and 6 show the 3D trend charts of the peak stress σ 1 and second principal stress σ 2 under different b values, respectively.The 3D trend chart of σ2 under different b values, shown in Figure 6, was the same as that of the peak stress σ1 because σ2 and σ1 were dependent on each other.When the stress in the σ1 direction reached the peak, the sample was broken, and there was no room for growth in σ2.As the confining pressure increased, the sample was sufficiently dense, and   The 3D trend chart of σ2 under different b values, shown in Figure 6, was the same as that of the peak stress σ1 because σ2 and σ1 were dependent on each other.When the stress in the σ1 direction reached the peak, the sample was broken, and there was no room for Figure 5 shows that as the value of b increases, the slope of the curved surface becomes steeper, indicating that an increase in b directly increased the stress in the σ 2 direction, the average principal stress increased, and the frozen soil was compacted.Although the second principal stress also increased, it accelerated the melting effect of ice crystals in the frozen soil and increased the restraint and occlusion of soil particles.Therefore, the strength significantly increased when b = 0 or b = 0, while the peak stress between b = 0.33 and 1 gradually decreased.This was because although the stress in the σ 2 direction increased, the occlusion between the soil particles was increased to a certain extent, and the strength was increased; however, the σ 3 direction constraint remained unchanged.The principal stress σ 1 significantly increased, and the strength of the sample decreased as the moisture content increased; the change in the moisture content cannot significantly affect the first principal stress σ 1 .For example, when the confining pressure was 1 MPa, the first principal stress was 4.895 MPa for 15% water content, and the first principal stress was 2.972 MPa for 22.5% water content, which was reduced by 1.923 MPa, and the rate of change was 39.3%.The confining pressure had a significant influence on the first principal stress σ 1 .For example, when the water content was 15%, the first principal stress with a confining pressure of 1 MPa was 4.895 MPa, and the first principal stress with a confining pressure of 4 MPa was as high as 11.193 MPa, which was an increase of 6.298 MPa, and the rate of change was 56.3%.
The 3D trend chart of σ 2 under different b values, shown in Figure 6, was the same as that of the peak stress σ 1 because σ 2 and σ 1 were dependent on each other.When the stress in the σ 1 direction reached the peak, the sample was broken, and there was no room for growth in σ 2 .As the confining pressure increased, the sample was sufficiently dense, and it was difficult for dilatancy to occur.The restraint and compaction effect of σ 2 became insignificant; therefore, the distance between the curved surfaces became smaller as the value of b increased; the water content increased under freezing conditions.The most intuitive effect of the enhancement in the peak stress of the soil was the volume content of the ice crystals in the frozen soil and the cementation of the ice crystals on the soil particles.However, the increase in the confining pressure also caused compression and melting and weakened the contribution of the increase in the water content to the increase in the peak stress.Therefore, the confining pressure was the main significant factor for the increase in the peak stress and the second principal stress, whereas, for the stress in the σ 2 direction, the b value was the most significant factor.For example, when the confining pressure was 4 MPa, the Bishop parameter b increased from 0 to 1, and the corresponding second principal stress increased from 4 to 9.128 MPa, which was an evident increase.

Levels of Analysis
The range analysis method was conducted to obtain the order of priority among the influencing factors without considering the interaction of the various factors.However, it was not possible to obtain the influence of the factors and their levels on the test results' influence weight.Therefore, the analytic hierarchy process was introduced to explore the weights of the above five factors on the strength of the frozen calcareous clay.

Hierarchical Analysis Theory
Similar to Table 2, considering its more general situation, we supposed that there were k influencing factors, namely A (1) , A (2) , . . ., A (k) , and the number of levels of each factor was n 1 , n 2 , . . ., n k .As shown in Figure 7, the AHP model of the orthogonal experiment had a three-layer structure comprising the index, factor, and level layers.

Hierarchical Analysis Theory
Similar to Table 2, considering its more general situation, we supposed that there were k influencing factors, namely A (1) , A (2) ,…, A (k) , and the number of levels of each factor was n1, n2, …, nk.As shown in Figure 7, the AHP model of the orthogonal experiment had a three-layer structure comprising the index, factor, and level layers.3), and S is expressed in Equation ( 4).
Therefore, matrix A right multiplied by matrix S was the normalization process of each column of matrix A, and matrix AS was called the level standard influence effect matrix.The range of the factors was i R (i = 1, 2, …, m), which was called the influence of Therefore, the sum of the test data under the j-th level of the defined factors were called the influence effect of the j-th level of the factor on the test, where i = 1, 2 . . ., k; j = 1, 2 . . ., n k .If the value of the inspection test index was high, let M ij = K ij ; otherwise, let M ij = 1/K ij .Among them, the effect matrix of the horizontal layer on the experiment is shown in Equation ( 3), and S is expressed in Equation ( 4). where Therefore, matrix A right multiplied by matrix S was the normalization process of each column of matrix A, and matrix AS was called the level standard influence effect matrix.The range of the factors was R i (i = 1, 2, . . ., m), which was called the influence of factors on the experiment.The weight matrix of the influence of factors on the experiment is expressed in Equation ( 5): Finally, the influence weight of each factor level on the index = ASC T was obtained by the analytic hierarchy process.was (n 1+ n 2+ . . .n k+ ) × 1 vector, which indicated the influence of each level of the factor on the experiment.

Discussion on the AHP Result
The above method yields the following: C = 0.155 0.262 0.077 0.398 0.108 (6) Appl.Sci.2022, 12, 8712 Therefore, the degree of influence of the various factors on the peak intensity σ 1 was as follows: From the calculation result obtained using Equation ( 9) without considering the interaction, the four levels of each factor have the highest influence weight on the peak stress σ 1 : factor A was A 4 = 0.0431; factor B was B 1 = 0.0826; factor C was C 2 = 0.0206; factor D was D 4 = 0.1338; and factor E was E 4 = 0.0289.Thus, the optimal level of the experiment was A 4 B 1 C 2 D 4 E 4 , and the primary and secondary order of influence of each factor was determined by the sum of the weights at each level of the five factors.The order of the primary and secondary influence of each factor was D > B > A > E > C. Similarly, the degree of influence of each factor on the intermediate principal stress σ 2 was obtained.As expressed in Equation (10), the weights with the highest influence were A 2 = 0.397; B 1 = 0.0541; C 2 = 0.0166; D 4 = 0.0896; and E 4 = 1.2239.The optimal level of σ 2 was E > D > B > A > C. Both these results were consistent with the range analysis results.Moreover, the influence weight of each level of the factor on the peak intensity was obtained through the analytic hierarchy process.

Significance Analysis Method
The magnitude of the variation was characterized by the sum of the squares of the variation, which was defined as: where, The estimated value of the deviation SS i was expressed as: where f was the degree of freedom, and f T = N − 1, N was the total number of trials, Based on the above deviation definition method, the calculation process of the significance level of the five-factor four-level orthogonal test reported in this paper was as follows: The sum of the squared deviations S T was: From the data listed in Table 3, we obtained η= 645.63; then, the sum of the squared deviations of each factor S A , S B , . . ., S E was expressed as: In the formula: T Ai , T Bi , • • •, T Ei were the sum of the test values of each factor at a certain level, e.g., T Ai was the sum of the four groups of test values at the level of factor 1.
The sum of squares of the error effect deviation S e was: The sum of the squared deviations of the error effect reflected the fluctuation value of the test error in this set of data.From Equations ( 13)- (15), the sum of the squared deviations corresponding to each factor was obtained.For example, the calculation process of S A was: Table 7 presents the rest of the calculation results.Referring to the corresponding degrees of freedom of the factors according to Equation (12), the estimated value of the deviation SS i of each factor was determined.Considering the ratio of deviation estimates: For a given significance level Ω, if The influence of factor i on the test index was significant; otherwise, it was not significant.

Significance Analysis without Considering Interaction
The test data listed in Table 3 were calculated and tested using the above calculation method, and the critical values of the F test listed in Table 6 were compared.Table 7 lists the test results.
The results of the significance study of the peak stress σ 1 showed that under complex stress conditions (under true triaxial conditions), the confining pressure and water content had significant effects on the strength, and the temperature effect was significant when Ω > 0.05.The influences of the salt content and b value on the intensity were not very significant.The order of significance of each factor was consistent with the results of the range analysis method and analytic hierarchy process.

Significance Analysis Considering Interaction
In the real state, the influence of frozen soil strength was not the result of single factors that were independently affected and then simply superimposed.The interactive influence of various factors could not be ignored.For example, the water content was the factor determining the volumetric ice content of frozen soil.However, an increase in the salt content reduced the freezing temperature of the water in the soil, making it difficult to freeze, and the volumetric ice content naturally decreased, thereby reducing the strength.At the same time, the application of a confining pressure caused the compression-thaw effect in the frozen soil, and the ice crystal cementation and melting affected the intensity.There were complex interactions between the various factors affecting the intensity.The above analysis results that did not consider the interaction between the factors could only be qualitatively analyzed to a certain extent.Hence, it was necessary to make the following significance analysis considering these interactions.
According to the basic idea of the orthogonal experiment design, the L 32 (4 9 ) orthogonal experiment table was used to design the experiment.The orthogonal experiment table head design was implemented in strict accordance with the experimental procedures.Due to many factors, the table head was too large, and the table head design is not presented herein.This study only considered the interaction between factors A, B, and C. Some of the 32 groups of effective tests were consistent with the test conditions listed in Table 3, and the test results could be used.Most of the remaining test results were calculated from Equation (2).The data obtained were analyzed for significance using the method described in Section 4.1.Table 8 presents the results.The test results showed that for the significance levels listed in Table 6, the effects of the confining pressure, water content, and intermediate principal stress coefficient b on the peak stress were all significant; the interactions A × B, A × C, and B × C had different degrees of influence on σ 1 , and the significance of the three interactions exceeded the degree of influence of the salt content.When Ω > 0.01, the influence of A × B on the equivalent stress was significant, and its influence even exceeded that of the temperature.From the results listed in Table 8, we concluded that the order of significance of the influence of various factors on the peak stress was as follows: confining pressure D, water content B, intermediate principal stress coefficient b, interaction A × B, temperature A, interaction A × C, interaction B × C, and salt content C. Hence, when studying the strength of frozen soil, it was necessary to comprehensively consider the various factors and their interactions to more accurately characterize the mechanical behavior of frozen soil.

Conclusions
In this study, a true triaxial orthogonal test was conducted on frozen calcareous clay to explore its strength change characteristics under the influence of multiple factors, including the degree of influence of each factor.The following conclusions were drawn from the results: (1) The characteristics of the stress-strain curve of frozen calcareous clay could be divided into three stages: the strain was within 0.5%, showing linear elasticity, and under the effect of pressure, the ice crystals, and their cements were damaged or compressed, and they shrink; at approximately 5%, they showed plastic hardening.The soil particles and ice crystals in the frozen soil recombined and became denser, resulting in irreversible deformation.With the compression process, cracks bred and expanded.
The failure stage was manifested as strain hardening due to the test loading conditions.As the deformation increased, the stress also slightly increased.In this stage, the soil sample was squeezed out, and the volume deformation could not represent its true change.(2) Through analytic hierarchy process, the peak stress σ 1 of each factor under four levels was quantitatively calculated, and the maximum value of influence weight of σ 1 was A 4 = 0.0431; B 1 = 0.0826; C 2 = 0.0206; and D 4 = 0.1338; E 4 = 0.0289, the maximum and minimum weight was 0.1338 and 0.0197, respectively, and the optimal level was A 4 B 1 C 2 D 4 E 4 .By comparing the sum of the weights of the four levels of each factor, the order of primary and secondary influence consistent with the range method was obtained, which was D > B > A > E > C. Similarly, the maximum and minimum weight of σ 2 was 1.2239 and 0.014, respectively, and the order of σ 2 was E > D >B > A > C. (3) When the interaction was not considered under true triaxial conditions, the effects of the confining pressure and water content on the strength were always significant.The temperature effect was significant only when the significance level Ω > 0.05.The salinity and b value had a significant effect on the strength.The impact was not significant, and the ranking of the significance of each factor was consistent with the results of the range analysis method and the analytic hierarchy process.When considering the interaction, the interaction factors had different effects on the intensity.
When Ω > 0.01, A × B was equivalent, and the effect of stress was significant, even exceeding that of the temperature.When studying the strength of frozen soil, it was necessary to comprehensively consider the various factors and the interaction between them, to more accurately characterize the mechanical behavior of frozen soil.

Figure 2 .
Figure 2. True triaxial test platform for frozen soil.Figure 2. True triaxial test platform for frozen soil.

Figure 3 .
Figure 3. Three views of lapped rigid plates for frozen soil samples.

Figure 4 .
Figure 4. True triaxial stress-strain and volume deformation curves of frozen calcareous clay.

Figure 4 .
Figure 4. True triaxial stress-strain and volume deformation curves of frozen calcareous clay.
Appl.Sci.2022, 12, x FOR PEER REVIEW 8 of 16 pressure of 4 MPa was as high as 11.193 MPa, which was an increase of 6.298 MPa, and the rate of change was 56.3%.

First princ ipal stres s σ 1 /Figure 5 .
Figure 5. Three-dimensional trend chart of water content, confining pressure, and high principal stress under different b values [23].

Figure 6 .
Figure 6.Three-dimensional trend chart of water content, confining pressure, and intermediate principal stress under different b values [23].

Figure 5 .
Figure 5. Three-dimensional trend chart of water content, confining pressure, and high principal stress under different b values [23].

Figure 5 .
Figure 5. Three-dimensional trend chart of water content, confining pressure, and high principal stress under different b values [23].

Figure 6 .
Figure 6.Three-dimensional trend chart of water content, confining pressure, and intermediate principal stress under different b values [23].

Figure 6 .
Figure 6.Three-dimensional trend chart of water content, confining pressure, and intermediate principal stress under different b values [23].
Appl.Sci.2022, 12, 8712 9 of 15 weights of the above five factors on the strength of the frozen calcareous clay.

Figure 7 ..
Figure 7. Orthogonal test AHP model.Therefore, the sum of the test data under the j-th level of the defined factors were called the influence effect of the j-th level of the factor on the test, where i = 1, 2…, k; j = 1, 2…, nk.If the value of the inspection test index was high, let ij ij M K =

Table 2 .
True triaxial test factors of calcareous clay and their corresponding levels.(Note: The salt content in this table was the salt content in 100 mL of pure water, that is, the mass percentage of anhydrous CaCl 2

Table 4 .
[23]ysis of the degree of influence of various factors on the major principal stress σ 1[23].

Table 5 .
[23]ysis of the degree of influence of various factors on the intermediate principal stress σ 2[23].

Table 7 .
Significance test on influencing factors of peak intensity σ 1 .

Table 8 .
Significance test of factors affecting the peak intensity σ 1 when the interaction was considered.