In geotechnical engineering, hydraulic engineering, coastal and deep-water engineering, foundations are usually subjected to vertical loading with a long-period acting and simultaneous or cyclic loading components including horizontal loading and moment loading. The vertical loading V, horizontal loading H, and moment loading M can be transferred to the underground soil layers through the foundation beneath the building. Such a loading mode is defined as the combined loading mode. The prediction of the ultimate load capacity of such foundations carrying vertical, horizontal and moment loads is a fundamental engineering problem. In the recent two decades, many researchers attempted to find the closed form of the failure envelope locus and to derive the ultimate bearing capacity of the deep-water foundation under combined vertical, horizontal and moment loadings based on experimental laboratory small scale model tests and finite element analysis programs[1-6]. The failure envelope is defined as a convex curved face in vertical, horizontal and moment loading space when the general shear failure mode occurs.

The failure envelope of shallow foundations undergoing combined vertical, horizontal and moment loadings has become a fundamental soil mechanics research problem. Researchers recently analyzed the combined loading responses of shallow foundations by means of plasticity methods. The failure envelope in V-H-M load space and yield surface equations has been deduced for different foundation systems. A general review of different proposed three-dimensional yield surface equations of foundations installed in different soil types and subjected to combined V-H-M loadings is presented in the following.

For any foundation, there is a three-dimensional yield surface equation containing all combinations of loads V, H, and M, which represents a failure envelope for the foundation. Any combination of loads inside the failure envelope equation can be regarded as a safe load combination for the foundation. This idea has been supported at the model scale by a large number of experiments on shallow strip footings on the sand by Butterfield and Gottardi[7]. They showed a simple three-dimensional failure envelope in V-H-M load space and defined the failure envelope as the location of the end-points of all conceivable load paths to the failure of footing. The failure locus on the V-H plane or the V-M plane is a parabolic curve whereas on the H-M plane, at a vertical load ratio V/Vu=0.5, the failure locus is an ellipse rotated through an angle of 13° from the H axis toward the negative direction of the M axis.

A finite element study of the failure locus for strip foundations on non-homogeneous undrained clay under combined V, H, and M loading was proposed by Bransby and Randolph[8]. New information has been obtained about the shape of the yield locus and the soil deformation mechanisms at yield from the finite element analyses. The shape of the yield locus was found to be similar to that predicted by previous workers in V-H and V-M load spaces but differed significantly in H-M load space. The soil deformation mechanisms calculated in the finite element analyses were supported using upper-bound plasticity mechanisms. A three-dimensional equation that approximates the shape of the yield locus is suggested, with a simplifying transformation for the yield locus in H-M load space. The transformation of the yield locus in H-M load space is obtained by the function of the height of the rotation center above the base when foundations become unstable under the pure moment.

Taiebat and Carter[9] presented a series of three-dimensional finite element analyses in order to investigate the failure envelope of circular foundations on the surface of homogeneous, purely cohesive soil subjected to combined V-H-M loadings. They showed two and three-dimensional failure envelopes in the non-dimensional space for foundations under combined loads and moment, based on numerical predictions. A three-dimensional equation that approximates the shape of the failure locus was also suggested to deal with the bearing capacity of these foundations. In the formulation of this equation, the moment capacity of the foundation is related to the horizontal load acting simultaneously on the foundation.

Loukidis et al.[10] presented a series of two-dimensional finite element methods to determine the collapse load of a rigid strip footing placed on a uniform layer of purely frictional soil subjected to eccentric and inclined loading. Two series of analyses were performed, one using an associated flow rule and one using a non-associated flow rule. Both probe-type analyses and swipe-type analyses were undertaken. Three different equations are proposed to predict the three-dimensional failure envelope of these foundations subjected to combined V-H-M loadings. The failure locus on the V-H plane or the V-M plane is a parabolic curve whereas on the H-M plane, at different vertical load ratios, the failure locus is a rotated ellipse. The elliptical equation parameters H|e=0 and M|α=0 are equal to the H and M values yielded by the equations used to predict V-H and V-M cross-sections, respectively.

In the present study, a series of three-dimensional finite element analyses using the ABAQUS program are performed to investigate the shape and locus of the failure envelope of a caisson foundation installed into homogenous sand soil, where the foundation embedment ratio D/B=1.0, and the soil caisson interface condition friction coefficient μ=0.3. The capacities of caisson foundations under combined V-H, V-M and H-M load spaces are studied and expressed by failure envelopes. Different failure mode mechanisms that occurred during the system (soil-foundation) failure are investigated and presented in this work. Finally, equations are proposed to predict different cross-sections of the failure envelope of a caisson foundations under combined V-H-M loadings.

1 Numerical Modelling

1.1 System soil-foundation geometries

The simplified analysis of the numerical modeling results provided a basis for developing the relationships between different applied vertical, horizontal and moment loads in the form of a complete failure surface in V-H-M space.

To investigate the failure envelope of a caisson foundation under combined vertical, horizontal and moment loadings, a three-dimensional finite element simulation was carried out using the ABAQUS program on a circular concrete caisson foundation embedded in homogeneous sandy soil. The caisson foundation embedment ratio is D/B=1, and the soil-caisson interface friction coefficient is μ=0.3. Fig.1 shows the modeled caisson foundation. The depth of the caisson foundation is considered to be constant and equals 9.0 m. The outside diameter of the caisson is 9 m, and the wall thickness is 0.25 m. The inside diameters vary from 4.8 to 5.2 m. The thickness of the lower part (tip-face) of the caisson is 0.4 m. The geometrical dimensions of the surrounding sand soil (see Fig.2) depend on the boundaries conditions which are placed far enough to avoid its influence on the different soil deformation mechanisms mobilized during the system’s failure and defined as length×width×depth (40 m×40 m×40 m).

1.2 Sand soil and caisson foundation properties

In the analysis, the sandy soil cohesion c=0.5 kPa; the frictional angle φ=20°; elastic Young’s modulus E and Poisson’s ratio ν are 80 MPa and 0.3, respectively. The mass density and the dilatancy angle of sand soil are 1 800 kg/m3 and 0°, respectively. In the concrete caisson foundation, the mass density is 2 500 kg/m3, elastic Young’s modulus E=30 GPa and Poisson’s ratio ν=0.2.

1.3 Contact modeling

The interface conditions between deep foundations and the surrounding soil are defined by interaction properties in normal and tangential behaviors. However, the caisson is connected to the soil with special contact surfaces (using the Master-Slave option), allowing for realistic simulation of possible detachment and sliding at the soil-caisson interfaces. The normal behavior or the surface-to-surface contact interaction was modeled by hard contact pressure-overclosure formulation, and the hard contact formulation shows that no contact pressure is transmitted without contact between the surfaces. Also, there is no limit to the contact pressure that can be transmitted when the surfaces are in contact. The hard contact formulation can be used along with the augmented lagrange, penalty and direct standards; i.e., the constraint enforcement method is selected as default in this study, and the separation contact between the respective surfaces is allowed all the time. The tangential behavior or the slippage of the caisson along the interfaces is governed by Coulomb’s friction law by appropriately assigning friction coefficient μ.

1.4 The finite element mesh

The finite element mesh used in the analysis is shown in Fig.2. A thin layer of continuum elements was used in the region of the soil-caisson interface to improve the prediction of the response of the caisson with the surrounding soil. The mesh comprises of 13 989 hex elements for both soil and caisson instances with the element type of C3D8R: an 8-node linear brick, reduced integration, and hourglass control.

1.5 Load paths

Three different steps are followed in this analysis to simulate the actual field conditions. In the first step, the numerical simulation soil undergoes geostatic loading. In the second step, part of the sand soil will be removed and replaced by the caisson foundation.The third step is defined as the load step, in which the caisson foundation will be subjected to a vertical load applied at the load reference point Pref, as shown in Fig.3.

In this work, the failure modes of the foundation are investigated by the load-controlled method. The loading was started with a single vertical load applied at the reference point on the top face of the caisson foundation and was increased gradually until failure was reached. Afterward, multiplying the vertical load by different ratios varying from 0.8 to 0 and applying another horizontal load or bending moment at the same load reference point Pref was done to obtain a vertical horizontal V-H or a vertical moment V-M plane, respectively. When the vertical load ratio equals zero, the foundation is then under pure horizontal or pure moment loading.

For the vertical, horizontal and moment load V-H-M space, the applied vertical failure load was reduced and multiplied by a vertical load ratio of 80%. Then a moment loading was applied at the same load reference point and increased gradually until failure was reached. This failure moment loading was multiplied by different ratios varying from 0.8 to 0 (for positive load combinations), and varying from zero to more than one (for negative load combinations). Then, another horizontal load was applied at the same reference point Pref and increased gradually until failure occured. The shape of the V-H-M failure envelope was identified when the vertical failure load was reduced to more ratios varying from 80% to 0%. When the applied vertical load reduced to 0%, the foundation is then under H-M loadings.

2 Results and Analysis

2.1 Vertical bearing capacity

The load-settlement curves were used to determine the vertical bearing capacity of the caissons.Fig.4(a) shows the load-settlement results obtained by the finite element analysis of the caisson under pure vertical loads plotted against the vertical displacements, where the depth to breadth ratio D/B=1, and the interface friction coefficient between the caisson and surrounding soil is μ=0.3.

The present results show that the load-settlement curve of the caisson foundation deflected gradually without the inflection point. This failure mode corresponds to a punching failure mode, and the vertical ultimate bearing capacity vu depends on the foundation settlement. From Fig.4(a), we can see that the vertical bearing capacity cannot be reached even when the caisson displacement is larger than 0.3B. Wang and Jin[11] studied the failure loci of suction caisson foundations under combined loading conditions and defined the vertical load capacities as the cut-off values, or the load corresponding to a vertical settlement of 0.05B. Liu et al.[12] defined the vertical load capacity where the vertical settlement is equal to 0.06B, and B is the diameter of the caisson foundation. In order to ensure the normal function of the caisson foundation, the ultimate vertical loading corresponding to 0.05B was taken as the ultimate bearing capacity vu.

Fig.4(b) shows the undeformed shape of the failure mechanism of this caisson foundation under pure vertical loading by incremental displacement vectors at the point of failure for the caisson foundation. The failure mechanism of a caisson foundation penetrating into sand soil is significantly different from shallow footings resting on sandy soil. Caissons translate vertical and lateral loadings from the top of the caisson to the surrounding lateral soil in front of the caisson and the underlying soil at the tip surface level. Shallow footings resting on sand soil translate vertical and lateral loadings from the top of the foundation to the underlying soil at the tip surface level only.

2.2 Horizontal bearing capacity

Fig.5(a) shows the behavior of load-displacement results obtained by the finite element analysis of the caisson under pure horizontal loads plotted against the horizontal displacements. Fig.5(b) shows that the failure mechanism of the caisson foundation is under pure horizontal loading by incremental displacement vectors at the point of failure for the caisson foundation. Based on the terminology introduced by Gerolymos et al.[13], the failure mechanism observed for the caisson foundation under pure horizontal loading is classified as a scoop slide mechanism. But in this study, the failure mode is similar to the scoop mechanism failure mode observed for the case when the foundation is under pure moment loading. The main reason for this is when the horizontal load is applied to the foundation, the rotation is not constrained. Thus, this failure mode that occurred can be classified as a scoop mechanism.

2.3 Moment bearing capacity

Fig.6(a) shows the behavior of the moment-rotational angle θ° curve and the ultimate moment capacity mu of the caisson foundation obtained by finite element analysis, where the caisson foundation bearing capacity under pure moment increases linearly. Fig.6(b) shows a snapshot of the failure mechanism of this caisson foundation under pure moment loading by incremental displacement vectors at the point of failure for the caisson foundation. According to the previous terminologies, the category of the failure mechanism observed under pure moment loading is classified as a scoop mechanism.

2.4 The failure envelope in VHM space

2.4.1 The failure envelope on V-H plane at M=0

Fig.7 shows the cross-section of the failure envelope on the V-H plane when the caisson foundation is under combined vertical and horizontal loads. The horizontal axis represents the vertical failure load components v normalized by the ultimate bearing capacity of the caisson foundation under pure vertical load vu, while the vertical axis represents the horizontal failure load components h normalized by the ultimate bearing capacity of the caisson foundation under pure horizontal load hu. The vertical v and the horizontal h are failure load components, which are obtained as the ultimate vertical vu and horizontal hu bearing capacities, respectively. In this figure, we can see that as the ratio of the vertical load v/vu increases, the ratio of the horizontal load h/hu increases until it reaches a peak (h/hu)max =1.6. This maximum horizontal load coincides with a vertical load ratio v/vu=0.54. However, as the ratio of the vertical loads (v/vu) becomes significant, the ratio values of the horizontal loads h/hu reverse and decrease until they become zero. In other words, as the vertical components of the failure load v approaches the ultimate bearing capacity of the foundation under pure vertical load, the horizontal load capacity declines sharply until it is nullified.

2.4.2 The failure envelope on V-M plane at H=0

Fig.8 shows the failure envelope of vertical and moment load combinations (V-M). The horizontal axis represents a ratio of the vertical failure load components v over the ultimate bearing capacity of the caisson foundation under pure vertical load vu, while the vertical axis represents the ratio of the moment failure load components m over the ultimate bearing capacity of the caisson foundation under pure moment loading mu. The vertical v and moment m failure load components are estimated as the ultimate vertical vu and moment mu bearing capacities, respectively. In Fig.8, we can also see that the maximum ratio of the normalized moment loading (m/mu)max=1.61 corresponds to a normalized vertical loading equal to 0.54. Also, when the values of the vertical load ratio vary from 0 to 0.54, the moment load ratios increase until they reach the maximum load (m/mu)max, while for the case when the vertical load ratios vary from 0.54 to 1, the moment load ratio decreases until it becomes nullified.

2.4.3 The failure envelope on H-M plane at V=0

A series of numerical analysis on the H-M plane with V/Vu=0 was performed. The failure envelope obtained from the analysis of the horizontal and moment loading is illustrated in Fig.9(a). We can see that the H-M failure envelope has a shape of half an oblique ellipse. In the case when the horizontal and moment act in a positive direction, the values of horizontal and moment load ratios v/vu and h/hu are less than or equal to one. While for the case when the horizontal and moment loads act in the opposite direction, the applied overturning moment M counter balancing the applied horizontal load H confirms the efficiency of negative load combinations in increasing and enhancing the ultimate bearing capacity of the caisson foundation. For example, the horizontal and moment failure loads can reach their maximum values hmax and mmax, where hmax is 2.77 times higher than the bearing capacity of the foundation under pure horizontal load, and mmax is 2.44 times higher than the foundation capacity under pure moment loading.

Fig.9(b) shows the normalized H-M failure envelope corresponding to different vertical load ratios V/Vu≥0. It is obvious that the failure envelope has a shape of a half oblique ellipse although there is influence from the V/Vu ratio on the expansion of the failure envelope shown for the case when V/Vu=0, i.e. the vertical load ratios increase the failure envelopes’ expanding until h/hu and m/mu obtain their maximum values, where the vertical load ratio V/Vu=0.2. However, the failure envelope starts to reduce and shrink back when the vertical load ratio V/Vu

≥ 0.2. We can also see that the maximum horizontal and moment failure loads corresponding to a vertical load ratio V/Vu=0 are smaller compared to the enhanced horizontal and moment failure loads when the vertical load ratio V/Vu=0.2. Thus, to investigate the failure mode mechanisms of this foundation, we take into account the horizontal-moment H-M cross-sections corresponding to vertical load ratios V/Vu of 0 and 0.2, respectively.

2.5 Deformation mechanisms

Zafeirakos and Gerolymos[14] used a finite element analysis program to study the soil deformation mechanisms of a caisson foundation in frictional soil subjected to combined V-H-M loadings, with a D/B ratio varying from 1 to 3 and the interface friction coefficient μ=1. This study is similar to that presented by Gerolymos et al.[13] on caisson foundations in cohesive soil subjected to static and dynamic loadings. Gerolymos et al.[13] presented the contours of resultant displacement at failure and defined the location of the rotation pole Zp (Zp is the distance from the load reference point to the center of the rotation of contours) to analyze different failure mode mechanisms corresponding to specific points on the horizontal-moment H-M cross-sections (see Fig.9(a)).

To investigate different failure mode mechanisms in this study, we first focus on the horizontal-moment (H-M) cross-section corresponding to the vertical load ratios V/Vu=0. Also, the terminology introduced by Gerolymos et al.[13-14] is used to classify different failure mode mechanisms.

Fig.10 shows snapshots of the failure mode mechanisms corresponding to a vertical load ratio V/Vu=0. At point a, a scoop mechanism is observed when the foundation is subjected to a pure moment loading as shown in Fig.10(a). Scoop mechanisms are developed between points a and c, similar to the case when the foundation is subjected to a pure horizontal loading (see Fig.10(b)). The rotation point of the vector displacement of different failure mode mechanisms was found to move from the lower-left part (under pure moment loading) to the lower-right part of the caisson (under pure horizontal loading).

At point c, the location of the vector displacement rotation point appeared on the tip face of the caisson foundation as illustrated in Fig.10(c). This failure mechanism is called the inverted pendulum mechanism. At point d, where the maximum horizontal load ratio h/hu is attained, the location of the vector displacement rotation point Zp is defined by Zafeirakos and Gerolymos[14] to be infinity (see Fig.10(d)). At point e, where the maximum moment ratio m/mu is attained, the vector displacement rotation point Zp is located at a depth corresponding to a half embedment of the caisson foundation as illustrated in Fig.10(e). And finally, a reverse scoop mechanism appeared at point f corresponding to the remaining part of the H-M cross section (see Fig.10(f)).

Figs.11 shows the snapshots of the failure mode mechanisms corresponding to a vertical load ratio V/Vu=0.2 by the incremental displacement vectors. Feld[15] studied the failure mode mechanism for combined V-H-M loading of the bucket foundation subjected to large vertical loads and defined three different failure mode mechanisms (shallow rotational failure, intermediate rotational failure, and deep rotational failure). The illustrated snapshot of the caisson foundation subjected to vertical and moment loadings corresponds to the intermediate rotational failure mode presented by Feld[15](see Fig.11(a)). The failure mode mechanism presented in this case is similar to that presented for the caisson foundation under pure moment loading, as shown in Fig.10(a). However, the rotation point of the vector displacement presented in Fig.11(a) seems slightly pushed down by the presence of the vertical load, i.e. the presence of the vertical load enhanced the bearing capacity of the moment loading up to about 70%, as shown in Fig.9(b).

A scoop slide mechanism is observed between points a and c, which is similar to the case that the foundation is subjected to a pure horizontal loading as illustrated in Fig.10(b). However, the rotation point of the vector displacement presented in Fig.11(b) was found to have moved from the right side to the left side of the caisson and being slightly pushed down by the presence of the vertical load. The presence of the vertical load enhanced the bearing capacity of the horizontal loading to about 60%, as shown in Fig.9 (b).

Fig.11(c) illustrates a failure mode mechanism, which is similar to the deep rotational failure presented by Feld[15]. In this case, the location of the rotation point of vector displacement appearing on the tip surface of the caisson (see Fig.10(c)) is moved from the right side to the left side by the presence of the vertical load. Fig.11(d) shows the failure mode mechanism corresponding to the maximum horizontal load of a caisson foundation under combined V-H-M loading with a depth to breadth ratio equaling 1, and the vector displacement rotation point seems to be infinity.

Fig.11(e) shows a maximum moment ratio (m/mu) of a caisson foundation under combined (vertical, horizontal and moment) loadings with a depth to breadth ratio equaling 1. In this case the vector displacement rotation point Zp is located at a depth corresponding to a half embedment of the caisson foundation. This failure mode mechanism is similar to the shallow rotational failure mode presented by Feld[15]. Finally, a reverse scoop mechanism at the remaining part of the failure envelope is presented in Fig.11(f).

2.6 General failure equations

This study aims to develop new formulas to predict the locus of the failure envelope of caisson foundations subjected to combined vertical-horizontal V-H, vertical-moment V-M and horizontal-moment H-M loadings.

The V-H, V-M and H-M cross-sections of the failure envelopes are obtained after normalization of the different vertical, horizontal and moment failure loads by the corresponding ultimate vertical, horizontal and moment loads, respectively. However, to better predict the results of the numerical analysis, a new equation is proposed to describe the shape of the V-H and V-M cross-sections of the failure envelopes(see Figs.12(a) and (b)). This equation can be written as

(1)

where v/vu, h/hu and m/mu are the vertical, horizontal and moment failure loadings normalized by their vertical ultimate bearing capacities, respectively.The values of the vertical ratio v/vu vary from 0 to 1.The parameters α1, α2, β1, β2 and β3 take the values of 0.5,15,1,0.730 6 and 1.730 6 for the V-H cross-section and 0.7,10, 1,0.907 3 and 1.907 3 for the V-M cross-section, respectively.

After replacing different parameters α1, α2, β1, β2 and β3 in Eq.(1), the equations used to predict the shape of the failure envelopes in V-H and V-M planes can be respectively rewritten as

(2)

(3)

Fig.12(c) illustrates the results of the H-M cross-sections of the failure envelopes obtained after renormalization of horizontal h and moment m loads by horizontal h* and moment m* loads corresponding to a fixed vertical load ratio V/Vu, respectively. The proposed equation to predict the locus of the H-M cross-sections can be written as

(4)

where h* is the horizontal failure load corresponding to a fixed vertical load ratio V/Vu in the V-H plane; m* is the moment failure load corresponding to a fixed vertical load ratio V/Vu in V-M plane; and the normalized horizontal and moment loadings h/h* and m/m* can be calculated from Eqs.(2) and (3), respectively.

The approximating half elliptical expression describes the inner boundary of the failure envelopes at each level of the vertical load ratio v/vu, where the values of the vertical load ratio vary from 0 to 1. The results of the failure envelope corresponding to the vertical load ratio V/Vu=0.6 seems slightly eccentric along the major axis of the half oblique ellipse. Although the approximate half elliptical expression is proposed to give a good fitting of different load combinations scenarios, it seems slightly conservative when the failure envelope approaches the maximum values h/h* and m/m*.

3 Conclusions

1) The response of caisson foundations under pure vertical, horizontal and moment loadings, as well as the results of two-dimensional cross-sections of the failure envelopes of caisson foundations in vertical-horizontal V-H, vertical-moment V-M and horizontal-moment H-M corresponding to different vertical load ratios V/Vu were presented and investigated. It is found that the presence of the vertical load at V/Vu=0.2 enhanced the horizontal and moment bearing capacities corresponding to the negative load combination plane to 60% and 70%, respectively.

2) The failure mode mechanisms of the caisson foundations under pure vertical, horizontal and moment loadings, as well as the accompanying failure mode mechanisms of different horizontal-moment H-M load combinations were investigated and presented, where the vertical load ratios V/Vu are of 0 and 0.2, respectively. Also, different failure mode mechanisms were defined according to the location of the rotation point of the incremental displacement vectors center Zp.

3) Three proposed formulas to describe the shape of the failure envelope of both vertical-horizontal V-H and vertical-moment V-M cross-sections of caisson foundations, as well as a half elliptical equation to describe the shape of the horizontal-moment H-M cross-sections corresponding to different vertical load ratios V/Vu, were presented. These proposed formulas allow us to predict the three-dimensional failure envelope of a caisson foundation subjected to combined vertical, horizontal and moment loadings.

4) The conclusions drawn in this work are valid for concrete caisson foundations in homogenous sand soil subjected to combined vertical, horizontal and moment loads, where the soil-caisson interface friction coefficient μ=0.3, and the depth to breadth ratio D/B=1.

References

[1] Murff J D. Limit analysis of multi-footing foundation systems [C]//Proceedings of the 18th International Conference on Computer Methods and Advances in Geomechanics. Rotterdam, the Netherland, 1994: 233-244.

[2] Gourvenec S, Barnett S. Undrained failure envelope for skirted foundations under general loading[J]. Géotechnique, 2011, 61(3): 263-270. DOI:10.1680/geot.9.t.027.

[3] Erik S. Bearing capacity failure envelope of suction caissons subjected to combined loading [D]. Norway: Norwegian University of Science and Technology, 2013.

[4] Mehravar M, Harireche O, Faramarzi A. Evaluation of undrained failure envelopes of caisson foundations under combined loading[J]. Applied Ocean Research, 2016, 59: 129-137. DOI:10.1016/j.apor.2016.05.001.

[5] Vulpe C, Gourvenec S M, Cornelius A F. Effect of embedment on consolidated undrained capacity of skirted circular foundations in soft clay under planar loading[J]. Canadian Geotechnical Journal, 2017, 54(2): 158-172. DOI:10.1139/cgj-2016-0265.

[6] Bagheri P, Son S W, Kim J M. Investigation of the load-bearing capacity of suction caissons used for offshore wind turbines [J]. Applied Ocean Research, 2017, 67: 148-161. DOI:10.1016/j.apor.2017.07.002.

[7] Butterfield R,Gottardi G. A complete three-dimensional failure envelope for shallow footings on sand[J]. Géotechnique, 1994, 44(1): 181-184. DOI:10.1680/geot.1994.44.1.181.

[8] Bransby M F, Randolph M F. Combined loading of skirted foundations[J]. Géotechnique, 1998, 48(5): 637-655. DOI:10.1680/geot.1998.48.5.637.

[9] Taiebat H A, Carter J P. Numerical studies of the bearing capacity of shallow foundations on cohesive soil subjected to combined loading[J]. Géotechnique, 2000, 50(4): 409-418. DOI:10.1680/geot.2000.50.4.409.

[10] Loukidis D, Chakraborty T, Salgado R. Bearing capacity of strip footings on purely frictional soil under eccentric and inclined loads[J]. Canadian Geotechnical Journal, 2008, 45(6): 768-787. DOI:10.1139/t08-015.

[11] Wang D, Jin X. Failure loci of suction caisson foundations under combined loading conditions[J]. China Ocean Engineering, 2008, 22(3): 455-464.

[12] Liu M M, Yang M, Wang H J. Bearing behavior of wide-shallow bucket foundation for offshore wind turbines in drained silty sand[J]. Ocean Engineering, 2014, 82: 169-179. DOI:10.1016/j.oceaneng.2014.02.034.

[13] Gerolymos N, Zafeirakos A, Karapiperis K. Generalized failure envelope for caisson foundations in cohesive soil: Static and dynamic loading[J]. Soil Dynamics and Earthquake Engineering, 2015, 78: 154-174.DOI:10.1016/j.soildyn.2015.07.012.

[14] Zafeirakos A, Gerolymos N. Bearing strength surface for bridge caisson foundations in frictional soil under combined loading[J]. Acta Geotechnica, 2016, 11(5): 1189-1208. DOI:10.1007/s11440-015-0431-7.

[15] Feld T. Suction Buckets, a new innovative foundation concept, applied to offshore wind turbines [D]. Denmark: Aalborg University, 2001.

竖向、水平与弯矩荷载联合作用下沉井基础破坏包络线研究

Abdellatif Boucheloukh 龚维明 戴国亮 竺明星

(东南大学土木工程学院， 南京 211189)

摘要：为研究竖向、水平与弯矩荷载联合作用下沉井基础承载特性，开展了砂土中圆形沉井基础在组合荷载作用下的有限元分析.沉井模型的高宽比为1，沉井-土体界面摩擦系数μ为0.3.首先开展了沉井基础单独受竖向荷载、水平荷载以及弯矩荷载作用下的沉井响应研究.其次基于荷载控制法，分别开展了竖向-水平联合荷载(V-H)、竖向-弯矩联合荷载(V-M)和水平-弯矩联合荷载(H-M)作用下的沉井基础归一化破坏包络线研究.针对不同竖向荷载比(V/Vu)情况，分析了竖向-水平-弯矩联合荷载(V-H-M)作用下的沉井承载特性，并进一步阐明了单向荷载及组合荷载下的沉井基础承载失效动态机理.最后，根据破坏包络三维形状特征建立了3个方程，通过这些方程可便利地得出单向荷载及竖向、水平与弯矩荷载联合作用下沉井基础承载力.

关键词：沉井基础；破坏包络线；有限元分析；承载力；砂性土

DOI:10.3969/j.issn.1003-7985.2020.02.009

Received 2019-10-08，Revised 2019-12-20．

Biographies：Abdellatif Boucheloukh (1982—), male, Ph.D. candidate; Gong Weiming (corresponding author), male, doctor, professor, wmgong@seu.edu.cn.

Foundation items：The National Natural Science Foundation of China (No.51808112, 51878160, 51678145), the Natural Science Foundation of Jiangsu Province (No.BK20180155).

Citation：Abdellatif Boucheloukh, Gong Weiming, Dai Guoliang, et al．Failure envelope of a caisson foundation under combined vertical, horizontal and moment loadings[J]．Journal of Southeast University (English Edition),2020,36(2):188-197．DOI:10.3969/j.issn.1003-7985.2020.02.009.

中图分类号：TU473.2