==== Front Sci Rep Sci Rep Scientific Reports 2045-2322 Nature Publishing Group UK London 37865 10.1038/s41598-023-37865-9 Article Study on distribution of sidewall earth pressure on open caissons considering soil arching effect Yang Bai 12 Shi Qingye 1 Zhou Hexiang 1223591918@163.com 3 Qin Chao 1 Xiao Weiwei shawvivi@163.com 1 1 grid.440723.6 0000 0001 0807 124X School of Architecture and Transportation Engineering, Guilin University of Electronic Technology, Guilin, 541004 Guangxi China 2 grid.263901.f 0000 0004 1791 7667 Faculty of Geosciences and Environmental Engineering, Southwest Jiaotong University, Chengdu, 610031 China 3 Sichuan Road and Bridge Group Co., Ltd, Chengdu, 610041 Sichuan China 30 6 2023 30 6 2023 2023 13 1065710 4 2023 28 6 2023 © The Author(s) 2023 https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. Based on the soil arching effect theory, the magnitude and distribution of sidewall earth pressure on open caissons when the embedded depth is large was analyzed by using theory of non-limit state earth pressure theory and horizontal differential element method. The theoretical formula was deduced. The theoretical calculation results are compared with the field test results and centrifugal model test results respectively. The results show that when the embedded depth of the open caisson is large, the distribution of earth pressure on the side wall of the open caisson first increases with the increase of embedded depth, reaches a peak value, and then sharply decreases. The peak point is located at 2/3 ~ 4/5 of the embedded depth. In engineering practice, when the embedded depth of the open caisson is 40 m, the relative error between the field test value and the theoretical calculation value is − 55.8% ~ 1.2%, with an average error of 13.8%. When the equivalent embedded depth of the open caisson in the centrifugal model test is 36 m, the relative error between the centrifugal model test value and the theoretical calculation value is − 20.1% ~ 68.0%, with an average error of 10.6%, The results are consistent well. The results of this article provides reference for the design and construction of open caisson. Subject terms Civil engineering Structural geology issue-copyright-statement© Springer Nature Limited 2023 ==== Body pmcIntroduction Open caisson foundation has the characteristics of high bearing capacity, high stiffness, good integrity, etc., it is widely used in large-scale bridge engineering1,2. How to ensure the safe and steady process of open caisson sinking is a key issue in the construction of the open caisson foundation. In order to solve the above problem, it is necessary to conduct in-depth research on the magnitude and distribution of the sidewall earth pressure of the open caisson. At present, the design of open caisson is mainly calculated by standard3, and the distribution law of sidewall friction is mainly based on the analysis results of sinking mechanism of large diameter pile and small open caisson4. The existing field monitoring and laboratory test results5–13 all indicate that the above calculation method is not applicable to large open caisson foundations. In the field test, Chen XP5 conducted real-time monitoring of the whole sinking process of the open caisson on a certain bridge's main pier, and preliminary analyzed the sinking mechanism and mechanical characteristics of the open caisson, and the magnitude and distribution of sidewall earth pressure at different depths of the open caisson are obtained. On the basis of the magnitude and distribution law, a calculation model for sinking resistance of the open caisson is established. In order to obtain the variation process of soil pressure during the first sinking stage of a large caisson, Jiang F et al.6 used finite element method to simulate the caisson foundation of the Changtai Yangtze River Bridge. Combined with field measured data of soil pressure at the blade foot, the variation law of soil pressure at the blade foot during the soil sampling process was obtained. Guo MW et al.7 proposed a new method for calculating the resistance at the end of the caisson based on the blade foot soil pressure. Chen BG et al.8 obtained the expression equation of caisson movement process by establishing the caisson settlement model, and revealed the kinematics characteristics and influence mechanism of caisson settlement. Lv CJ et al.9 conducted real-time monitoring of the sinking process of large open caissons in clay sand interactive formations, and further studied the magnitude and distribution of the reaction force, sidewall pressure, and lateral friction during the sinking process. In terms of indoor model tests, Wang Jian11 and Mu Bao-gang13 successively carried out model tests on sinking resistance of the open caisson, and studied the change characteristics of sidewall friction resistance, subsidence range of soil surface outside the wall, flow trend of soil particles and so on when the open caisson sank to different depths. At present, less research results in theoretical terms. The soil arching effect is the phenomenon that the stress in the soil is transferred from the yield region to the adjacent unyield region14. In 1943, Terzaghi15 proved the existence of soil arching effect for the first time through a trapdoor test. Later, in the study of earth pressure on earth retaining walls, Handy16, Zhao XY17, Liu H18, Lai F19 and Zhang H20, etc. considered the soil arching effect and obtained the calculation method of earth pressure. Jiang XD21 applied the theory of soil arching effect to the sinking calculation of the open caisson, and derived the coefficient of sidewall earth pressure and the calculation formula of sidewall earth pressure. However, it did not take into account the stress relaxation effect of the soil outside the cutting curb5, and there was no theoretical basis for the assumption that the active and passive limit states are above and below the critical depth, respectively. The non-limiting state earth pressure theory is improved from the Rankine or Coulomb earth pressure theory, which is widely used in the calculation of earth pressure on earth retaining wall22,23. The theory believes that the soil is a gradual process from the static state to ultimate state, and in the calculation of earth pressure, the influence of the displacement mode of retaining wall and displacement on the magnitude and distribution of earth pressure is considered. Similarly, the posture change during sinking process of the open caisson can also be regarded as a displacement of the sidewall relative to the soil. In this paper, based on the soil arching effect theory, the magnitude and distribution of sidewall earth pressure on open caissons when the depth of sinking is large was analyzed by using theory of non-limit state earth pressure theory and horizontal differential element method. The theoretical formula was deduced. The theoretical calculation results are compared with the field test results and centrifugal model test results respectively. The results can provide references for the design and construction of open caissons. Considering the internal and external friction angle of sidewall displacement Taking the open caisson sinking in soil with a horizontal non-cohesive soil surface as an example, in order to simplify the calculation, the working condition of the open caisson in a vertical state is taken for analysis, that is, the sidewall transversals relative to the soil. In the non-limit state, there is a quasi-slip surface in the soil. It is assumed that the quasi-slip surface is straight and the sidewall moves away from the soil, the angle between the quasi-slip surface and the sidewall is β=π4-φ2, and when the sidewall transversals towards the soil, the angle between the quasi-slip surface and the wall is β=π4+φ2. At this time, the internal friction angle φm of soil is between the initial internal friction angle φ0 and the ultimate internal friction angle φ, which can be calculated by Eq. (1)241 φm=arctantanφ0+Kdtanφ-tanφ0. Type: Kd=4arctan(S/Sc)/π; S is the actual translational displacement of the sidewall; Sc is the critical displacement of sidewall when the soil reaches the active or passive limit state. When the sidewall transversals horizontally, the corresponding displacement when the sand reaches the active limit state is about 0.001H ~ 0.005H, and when it reaches the passive limit state is about 0.05H. (H is the buried depth of open caisson). Considering the influence of friction angle δ at the interface between soil and sidewall, φ0 can be solved by Eq. (2)23.2 1K0=1cosφ0+tan2φ0+tanφ0tanδ122 For normally consolidated soil, K0=1-sinφ and to be conservative, to take δ=φ/2. Soil arching effect outside the sidewall of open caisson Stress model of principal stress arch When the open caisson sinks to a certain depth, due to excavation, the height of the soil surface inside the well is lower than that outside the well, and a soil pressure difference is formed inside and outside the cutting curb, which makes the soil outside the cutting curb have a tendency to migrate to the inside and form a pressure relaxation area in the soil, as shown in Fig. 1.Therefore, the soil outside the sidewall can be divided into two parts, non-stress relaxation zone and stress relaxation zone, from top to bottom.Figure 1 Schematic diagram of stress relaxation effect. When the sinking depth of the open caisson is large, the sinking rate of the open caisson sucking sludge is small or in a stagnant state. At this time, the soil in the non-stress relaxation zone moves upward relative to the sidewall and is affected by the downward friction of the sidewall, while the soil mass in the stress relaxation zone migrates to the inside of the open caisson through the cutting curb under the action of sucking sludge in the well, and thus is affected by the upward friction of the sidewall. According to the principle of soil arching effect, under the action of the sidewall friction, the soil stress outside the sidewall deflected, causing the stress in the yielding area to transfer to the adjacent unyielding area, forming the principal stress arch between the quasi-slip fracture surface and the sidewall. If the sidewall moves away from the soil, the soil is in an active state, forming a minor principal stress arch. If the sidewall moves toward the soil, the soil is in a passive state, forming a major principal stress arch. The calculation model is shown in Fig. 2. The sinking depth of the open caisson is H, the height of the non-stress relaxation area is H1, and the height of the stress relaxation area is H2.Figure 2 Soil arching effect outside of sidewall. Stress analysis of principal stress arch The differential element body at point A establishes the horizontal and vertical balance equations. In the passive state:3 σh=σ1sin2θ+σ3cos2θ 4 τw=σ1-σ3sinθcosθ 5 Kp=σ1/σ3 In the active state:6 σh=σ1cos2θ+σ3sin2θ 7 τw=σ1-σ3sinθcosθ 8 Ka=σ3/σ1 where θ is the included angle between σ3 and the horizontal direction; σh is the horizontal stress; τw is the shear stress; Kp is the Rankine passive earth pressure coefficient; Ka is the Rankine active earth pressure coefficient. As can be seen from the Mohr stress circle in Fig. 3, whether in the passive state or the active state, there is:9 σh=σ1+σ3-σv Figure 3 Mohr stress circle of the soil outside of sidewall. Type: σv is the vertical stress. The passive states of simultaneous (3), (5) and (9) can be obtained:10 σhσ3=Kpsin2θ+cos2θ 11 σvσ3=Kpcos2θ+sin2θ The active states of the simultaneous (6), (8) and (9) can be obtained:12 σhσ1=Kasin2θ+cos2θ 13 σvσ1=Kacos2θ+sin2θ Dividing Eq. (10) by Eq. (11) and Eq. (12) by Eq. (13), can be obtained, whether active or passive state, the theoretical lateral pressure coefficient is:14 K=σhσv=Ksin2θ+cos2θKcos2θ+sin2θ Type: K is the earth pressure coefficient, Kp is the passive earth pressure coefficient in the passive state, and Ka is the active earth pressure coefficient in the active state. The friction angle δ of the interface between the soil and the sidewall is known, and it can be derived from Fig. 3 that the value of θ in the passive state is:15 tanδ=Kp-1tanθKp+tan2θ Similarly, the value of θ in the active state is:16 tanδ=1-Katanθ1+Katan2θ Shape of stress arch Although the theoretical soil arch curve has been proved to be a part of catenary25, its expression is more complicated and can generally be simplified as a circular arc soil arch for calculation26,27. Taking the non-stress relaxation area as an example, assuming that the stress arch trace is a circular arc, as shown in Fig. 2, point A is the coordinate origin, the x-axis is positive horizontally to the outside of the wellbore, and the y-axis is positive vertically. The relative coordinates of the center of the circular arc arch are (B1,B1tanθ), then the geometric equation of the circular arc arch is:17 x-B12+y-B1tanθ2=B12cos2θ Type: B1 is the horizontal distance between the center of the circular arc arch and the sidewall. Practical lateral earth pressure coefficient As shown in Fig. 2, it can be known from the geometric relationship:18 cosψB1-x=cosθB1 Type: ψ is the included angle between the center line of the circular arc arch differential element and the horizontal direction. It can be obtained from the Molar stress circle, in the passive state:19 σmhσ3=Kpsin2ψ+cos2ψ 20 σmvσ3=Kpcos2ψ+sin2ψ Type: σmh is the horizontal stress on the circular arc arch differential element; σmv is the vertical stress on the circular arc arch differential element. Dividing Eq. (19) by Eq. (20), the practical lateral pressure coefficient in the passive state can be known as:21 Kpw=σhσ¯v1=σh/σ3σ¯v1/σ3=1σ¯v1/σ3Kpsin2θ+cos2θ Type: σ¯v1 is the average vertical stress on the horizontal differential element, which can be calculated by σ¯v1=∫0B1σmvdxB1. Simultaneous Eqs. (18) and (20) and calculating through integrals can be obtained:22 σ¯v1σ3=Kp-1cos2θ13+1 Substitute into Eq. (21) to get:23 Kpw=1Kp-1cos2θ3+1Kpsin2θ+cos2θ Similarly, the practical lateral pressure coefficient in the active state is:24 Kaw=1Ka-1cos2θ3+1Kasin2θ+cos2θ Comparing Eq. (23) and Eq. (24), it can be seen that regardless of the passive state or the active state, the practical lateral pressure coefficient is:25 Kw=1K-1cos2θ3+1Ksin2θ+cos2θ Earth pressure based on horizontal differential element method Establishment of basic equations Non-stress relaxation area Take a horizontal differential element with thickness dh at the distance h1 from the surface of the non-stress relaxation area, and the force acting on the horizontal differential element is shown in Fig. 4.Figure 4 Stress of horizontal differential element outside of sidewall. From the equilibrium conditions of the horizontal differential element in the vertical direction, it can be obtained:26 H1+H2-h1tanβdσ¯vfdh1-σ¯vfKwtanδ=γH1+H2-h1tanβ Type: γ is the soil weight; σ¯vf is the average vertical stress on the horizontal differential element in the non-stress relaxation area; β is the included angle between the quasi-slip surface and the sidewall. Taking the boundary condition as σ¯vf(h1=0)=0 , the special solution of σ¯vf is obtained by Eq. (26):27 σ¯vf=γA-h1cotβBcotβAA-h1cotβB-1 Type: A=(H2+H1)cotβ; B=1+Kwtanδcotβ. Earth pressure on the sidewall of the non-stress relaxation area:28 p1=KwγA-h1cotβBcotβAA-h1cotβB-1 Stress relaxation area Similarly, the vertical equilibrium equation of the force on the horizontal differential element at the distance h2 from the top surface of the stress relaxation area is as follows:29 (H2-h2)tanβdσ¯vsdh2+σ¯vsKwtanδ=γH2-h2tanβ Type: σ¯vs is the average vertical stress on the horizontal differential element in the stress relaxation area. Taking the boundary condition as σ¯vf(h1=H1)=σ¯vs(h2=0), and the special solution of σ¯vs is obtained from Eq. (22):30 σ¯vs=-γ(H2-h2)C+D(H2-h2H2)1-C Type: C=1-Kwtanδcotβ; D=σ¯vfh1=H1+γH2C. Earth pressure on the sidewall of the stress relaxation area:31 p2=-KwγH2-h2C-DH2-h2H21-C Range of stress relaxation area The height of the open caisson stress relaxation area can be determined by referring to existing field monitoring data. The research results of Chen XP5 showed that the range of stress relaxation area is 1-10 m, and the upper limit is taken when the sinking depth is large. Examples verification On-site monitoring of open caisson sinking of main pier of a bridge The size of the open caisson of the main pier of a bridge is 86.9 m × 58.7 m, the radius of the inverted circle is 7.45 m, and the height is 115 m.The open caisson is built in water and needs to pass through dense silty sand layer, dense fine sand layer and dense coarse sand layer successively when sinking. According to the field geological survey data, the saturation weight density of soil is γ=19.2kN/m3, and the angle of internal friction is φ=36.9∘. To be conservative, the angle of friction between the soil and the sidewall is δ=18.5∘. Seven monitoring sections were set at different heights around the outer wall of the open caisson. The sections were 2 m, 5 m, 19 m, 37 m, 51 m, 60 m and 72 m above the cutting curb, and is divided into section 1–section 7 according to the elevation of the section from high to low. The section layout of the sidewall earth pressure sensor is shown in Fig. 5, and the plane layout of each section is shown in Fig. 6.Figure 5 Monitoring sections of earth pressure sensor in sidewall (unit: cm). Figure 6 Plane layout diagram of earth pressure sensor in sidewall. Due to the limitation construction conditions of the open caisson was built in the water and traversed a huge thick soil layer, many earth pressure sensors were damaged and failed after entering the soil, and only some of the sensors in Sections “Earth pressure based on horizontal differential element method” and “Conclusions and suggestions” were intact. Selecting the working condition of the open caisson embedded depth H≈100 m and the open caisson is vertical to the north–south direction for analysis. At this time, the stress relaxation area h2 of the open caisson is 10 m reaches the passive limit state and the displacement Sc is 200 cm. According to the field measurement, the cumulative translational displacement in the south direction of the open caisson is S=72.3cm, the soil on the south side of the open caisson is in a passive state. The calculation results are compared with the field test results, as shown in Fig. 7. The calculated results of the method in this paper is basically consistent with the field measured results. There are few field test results. From the theoretical calculation curve, it can be seen that the sidewall earth pressure first increases with the increase of embedded depth, reaches a peak value, and then sharply decreases. The peak point is 32 m, which is 80.0% of the embedded depth.Figure 7 Comparison between the calculated results of this paper and the field measured results. The theoretical method was used to calculate the earth pressure of Sections “Earth pressure based on horizontal differential element method” and “Conclusions and suggestions” of the caisson. The calculated values were compared with the field values, and the results are shown in Table 1. The relative error between the field value and the theoretical calculation value is -55.8% ~ 1.2%, with an average error of 13.8%.Table 1 Comparison between theoretical calculation results and field test results. Embedded depth/m Theoretical calculation results/kPa Field test results/kPa Relative error/% 3.0 75.9 85.8 13.1 35.0 902.7 913.8 1.2 38.0 419.2 185.2 -55.8 Centrifugal model test of sinking resistance of open caisson The centrifugal model test of the open caisson sinking is used to study the distribution law of earth pressure on the sidewall of open caisson, the outer wall of open caisson is selected as the research object, and the three-dimensional model of the open caisson is simplified into a two-dimensional model of the sidewall. The sidewall model of 45 cm high, 70 cm wide, the width of the cutting curb tread surface is 0.3 cm, the width of the cutting curb slope is 1.7 cm, and the inclination angle is 45°. The open caisson is connected with a fixed bracket by a slide rail and can slide freely in the vertical direction after installation, as shown in Fig. 8.Figure 8 Layout diagram of centrifuge model test (unit: cm). It is very difficult to dynamically simulate the process of mud suction and sinking of open caisson under the current technical conditions of the centrifugal model test technology. In the test, the working condition of the open caisson embedded depth H≈36 m was selected, and the dynamic analysis was replaced by the quasi-static analysis. In addition, the loading system is used to control the vertical force to be constant to simulate the process that the vertical force balance of the open caisson is broken after a certain mud suction, and the open caisson is sank and reached rebalance by the action of self-weight stress. The centrifugal acceleration was set to 90 g, and the initial position of the cutting curb tread surface is level with the soil surface in the well. Silt sand was selected as the test soil, the saturation weight density of soil γ=19.8kN/m3, and internal friction angle φ=36.2∘. For the sake of conservativeness, take δ=18.1∘. The earth pressure sensor is used to measure the earth pressure on the sidewall of the open caisson. The arrangement of the earth pressure sensor is shown in Fig. 9.Figure 9 Layout diagram of earth pressure sensor in sidewall (unit: cm). Due to the limited machining accuracy, when riveting the open caisson model with the loading system, it is necessary to make a slight adjustment of fine-tune the surface of the open caisson model towards the direction of soil, the adjustment value is about 0.5 cm. After the conversion by similarity ratio, S=45cm, when the displacement reaches the passive limit state, Sc=180cm.The soil is in a passive state. Assuming that the stress relaxation area of the open caisson is 18 m, the calculated results of the method in this paper is compared with the results of centrifugal model test, the results are shown in Fig. 10. The calculation results of the method in this paper are in good agreement with the results of the centrifugal model test. The stress relaxation area of the centrifugal model is greater than 10 m, the reason is that when the test data are collected, the open caisson has not appeared obvious subsidence, resulting in a larger range of the soil with downward displacement relative to the sidewall, that is, a larger range of stress relaxation area. The equivalent embedded depth of the open caisson in the centrifugal model test is 36 m. From the theoretical calculation results curve and centrifugal model test results, it can be seen that the earth pressure shows a distribution characteristic of first increasing with the increase of embedded depth, reaching a peak value and then sharply decreasing. The theoretical calculation peak point is 24 m, and the centrifugal model test peak point is 23.4 m, which are 66.7% and 65.0% of the embedded depth, respectively.Figure 10 Comparison between calculated results of this paper and centrifuge model test results. The theoretical method was used to calculate the earth pressure of the sensor burial depth in the centrifugal model test. The calculated value was compared with the experimental value, and the results are shown in Table 2. The error between the experimental value and the theoretical calculation value is − 20.1% ~ 68.0%, with an average error of 10.6%.Table 2 Comparison between theoretical calculation results and centrifugal model test results. Embedded depth/m Theoretical calculation results/kPa Centrifugal model test results/kPa Relative error/% 5.4 131.1 184.0 40.4 14.4 384.4 465.5 21.1 23.4 537.4 514.2 − 4.3 27.9 527.5 460.0 − 12.8 29.7 495.9 396.1 − 20.1 31.9 454.6 372.4 − 18.1 34.2 204.2 343.0 68.0 Conclusions and suggestions Based on the principle of soil arching effect, the non-limit state earth pressure theory and horizontal differential element method are adopted to analyze the magnitude and distribution of the earth pressure on the sidewall of the open caisson, and the theoretical formula of the earth pressure on the sidewall of the open caisson is derived. The formula is suitable for working conditions where the sinking depth of the open caisson is large, the rate of mud absorption of the open caisson is small, or it is in a stagnant state. Based on the results of this article, selecting an appropriate sidewall friction coefficient of the caisson can calculate the frictional resistance value of the caisson sidewall, that providing an important basis for the design of the caisson and the stability judgment during the sinking process. Through centrifugal model tests and theoretical calculations, it was found that when the embedded depth of the open caisson is large, the distribution of sidewall earth pressure of the open caisson first increases with the increase of embedded depth, reaches a peak value, and then sharply decreases. The peak point is located at 2/3 ~ 4/5 of the embedded depth. The calculation model for the earth pressure on the sidewall of the open caisson established in this paper is simple and convenient with clear mechanical concept. In engineering practice, when the embedded depth of the open caisson is 40 m, the relative error between the field test value and the theoretical calculation value is − 55.8% ~ 1.2%, with an average error of 13.8%. When the equivalent embedded depth of the open caisson in the centrifugal model test is 36 m, the relative error between the centrifugal model test value and the theoretical calculation value is − 20.1% ~ 68.0%, with an average error of 10.6%.The calculation results are in good agreement with field measured results and the centrifugal model test results, which can provide reference for the design and construction of open caisson. The distribution of earth pressure on the side wall of deep and large caisson is affected by the attitude of the caisson, and the change of the attitude of the caisson can be divided into two modes: rotation and translation. This paper only gives the theory of computation of the earth pressure distribution on the side wall when the caisson moves horizontally without tilting displacement during the sinking process. The theory of computation of the earth pressure distribution on the side wall when the caisson rotates and tilts is more complex, and needs further research. Author contributions B.Y., Z.H. Q.C. and X.W. wrote the main manuscript text. S.Q. and Z.H. prepared figures. All authors reviewed the manuscript. Data availability The data used to support the findings of this study are available from the corresponding author upon request. Competing interests The authors declare no competing interests. Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. ==== Refs References 1. Zhang YT Yang Z Li DJ Zeng XT Zhu JT Development and prospects of construction technology for large open caisson foundations of bridges Bridge Constr. 2023 53 02 17 26 2. Shi Z Zhong ML Zhou YC Li SY Feng CB Summary of basic research on large bridge caisson China Railw. Sci. 2022 43 05 11 22 3. China Construction Standardization Association. Specification for structural design of reinforced concrete sinking well of water supply and sewerage engineering. CECS 137: 2015 (China Planning Press, 2015). 4. Dong XC Guo MW Jiang ZX Wang SL Chen ZW Calculation and analysis of resistance at the sinking end of super large open caisson based on soil pressure monitoring values Bridge Constr. 2022 52 05 78 84 5. Chen XP Qian PY Zhang ZY Study on penetration resistance distribution characteristic of sunk shaft foundation Chin. J. Geotech. Eng. 2005 27 02 148 152 6. Jiang F Liu H Yue Q Yang WS Study on the variation law of soil pressure at the foot of the edge of a large open caisson foundation during soil sampling and sinking Geotech. Mech. 2022 43 S2 431 442 7. Guo MW Study on the resistance variation law of the sinking end of the foundation of super large open caisson J Rock Mech. Eng. 2021 40 S1 2976 2985 8. Chen BG He JX Luo RP Zhang GH Gao Q Study on kinematics characteristics and settlement control of open caisson sinking process Geotech. Mech. 2022 43 S2 425 430+453 9. Lv CJ Huang R Ma YG Study on the sinking resistance of open caisson in clay sand interactive strata Bridge Constr. 2022 52 04 61 67 10. Chavda JT Mishra S Dodagoudar GR Experimental evaluation of ultimate bearing capacity of the cutting edge of an open caisson Int. J. Phys. Modell. Geotech. 2020 20 5 281 294 10.1680/jphmg.18.00052 11. Wang J Liu Y Zhang Y Model test on sidewall friction of open caisson Rock Soil Mech. 2013 34 03 659 666 12. Chavda JT Dodagoudar GR Experimental studies on a circular open caisson Int. J. Phys. Modell. Geotech. 2022 22 2 70 87 10.1680/jphmg.20.00050 13. Mu BG Bie Q Zhao XL Gong WM Meso-experiment on caisson load distribution characteristics during sinking China J. Highw. Transp. 2014 27 09 49 56 14. Wang J Xia TD He PF Huang B Analysis of active earth pressure on rigid retaining walls considering soil arching Rock Soil Mech. 2014 35 07 1914 1920 15. Terzaghi K Theoretical Soil Mechanics 1943 Wiley 66 76 16. Handy RL The arch in soil arching J. Geotech. Eng. ASCE 1985 111 3 302 318 10.1061/(ASCE)0733-9410(1985)111:3(302) 17. Zhao X Li K Xiao D A simplified method to analyze the load on composite retaining structures based on a novel soil arch model Bull. Eng. Geol. Env. 2020 79 3483 3496 10.1007/s10064-020-01780-4 18. Liu H Kong D Active earth pressure of finite width soil considering intermediate principal stress and soil arching effects Int. J. Geomech. 2022 22 3 04021294 10.1061/(ASCE)GM.1943-5622.0002298 19. Lai F Zhang N Liu S Yang DY A generalised analytical framework for active earth pressure on retaining walls with narrow soil Géotechnique 2022 10.1680/jgeot.21.00305 20. Zhang H Liu M Liu Z Li XL Analysis of the pile spacing and earth pressure of sheet pile walls based on the spatial soil arching model Acta Geotechnica Slovenica 2022 19 1 2 16 21. Jiang XD Study on critical depth of the maximum sidewall earth pressure on open caisson and its influencing factor Water Resour. Power 2016 34 10 110 113 22. Bang S Active earth pressure behind retaining walls J. Geotech. Eng. 1985 111 3 407 412 10.1061/(ASCE)0733-9410(1985)111:3(407) 23. Fang YS Ishibashi I Static earth pressure with various wall movements J. Geotech. Eng. 1986 112 3 317 333 10.1061/(ASCE)0733-9410(1986)112:3(317) 24. Chang MF Lateral earth pressure behind rotating wall Can. Geotech. J. 1997 34 2 498 509 10.1139/t97-016 25. Kingsley HW Arch in soil arching J. Geotech. Eng. ASCE 1989 115 3 415 419 10.1061/(ASCE)0733-9410(1989)115:3(415) 26. Wang YZ Investigation on the stratigraphic response and plugging effect induced by press-in open caisson in mucky soil KSCE J. Civ. Eng. 2023 10.1007/s12205-023-0294-7 27. Zhu X Li W Fei K Jacking penetration resistance and mechanical characteristics of bucket foundations in sand Bull. Eng. Geol. Env. 2020 79 1 1 10 10.1007/s10064-019-01552-9