Numerical Analysis on the Load Sharing Performance of Long-Short CFG Pile Composite Foundation Subjected to Rotation of Adjacent Retaining Wall

*e stress and displacement boundary conditions of excavation retaining structures affect the deformation mechanism and movement of the retained soil mass. *e soil movement disturbs the load sharing performance and structural integrity of cementfly ash-gravel (CFG) pile composite foundations existing in the vicinity, which merits considerable research work. *is article presents results from 3D finite element analyses performed to study the influence of retaining wall rotation on the load sharing characteristics of adjacent CFG pile composite foundation comprising long and short piles. To verify the numerical model, a relatively large-scale 1 g physical model test was conducted. It is revealed that to arrive at a new static equilibrium during progression of wall rotation, the percentage load sharing ratios of the long and short piles change increasingly while the load proportion carried by the upper soil reduces remarkably. *e percentage load sharing characteristics of CFG pile composite foundation are more affected in immediate proximity to the wall than those located at far distance. For the foundation having 3× 3 long and short piles placed at 3.0–15.0m away from the wall, the location resulted in a reduction of soil bearing capacity ranging between 1.4 and 7.5% of the total imposed load while the corresponding increase in the % load borne by the long and short pile range was 0.83–4.15 and 0.59–3.36%, respectively. For the other parameters considered in this article viz. pile spacing, subsoil stiffness, cushion stiffness and thickness, and applied working load, the increment in % load sharing of the long and short pile range was 3.45–4.15, 1.3–5.79, 1.48–3.36, 4.15–4.79, and 3.67–4.15% and 3.36–4.67, 1.43–4.99, 1.48–3.36, 3.36–3.64, and 1.38–3.36% of the imposed load, respectively. Moreover, the long piles’ load sharing proportion was higher than that of short piles, and peripheral piles received larger load proportion.


Introduction
Escalating urban land demands and growing population size in emerging cities have promoted the quest for possible alternative solutions to prevailing geotechnical problems, notable among which is enhancing poor ground conditions to make them suitable for construction. Engineers have to choose between the best available foundation treatment method and/or in combination with special foundation system that will allow overlaying structures to satisfy admissible technical requirements. CFG (cement-fly ashgravel) pile composite foundation with mechanical characteristics intermediate between conventional shallow and deep foundations is generally a disconnected type of piled raft foundation that serves to achieve the performance requirements in almost all kinds of weak soils to support high-rise buildings and other infrastructures [1]. e concept is more popular across Europe and Asia where serviceability requirement for mat foundation calls for the use of settlement reducing piles [2]. e CFG piles crowned with a cushion layer receive the foundation load indirectly; and owing to the deformability characteristics of the cushion, piles' head pricks into the mattress layer, leading to relative pile/soil settlement that provokes negative side friction (NSF) to develop on the upper 15-50% of the pile length [3][4][5][6][7]. In this regard, researchers such as Halder and Manna [3], Rui et al. [4], and Fioravante and Giretti [8] conducted centrifuge and static geotechnical model tests, and Tradigo et al. [6] and Ata et al. [9] performed numerical analyses while Jiang and Liu [5] and Wu et al. [7] established differential equation based on typical unit calculation model to deduce the calculation formula of the depth of neutral plane and pile-soil stress ratio.
In order to meet further economic advantages of CFG pile reinforced ground, a combination of different piles with varying stiffness [10,11] and/or length [12,13] has been used in practice. Earlier studies have disclosed that the natural unfavourable shallow soil would be strengthened by the short/flexible piles whereas the long piles play a major role in reducing the foundation settlement. Moreover, the influence of various components such as cushion's property and geometry, subgrade material stiffness properties, geometric and stiffness properties of piles, and raft thickness on the vertical load sharing proportion between the piles and soil and overall foundation settlement characteristics have been discussed by a number of authors [14][15][16]. e recent experimental study of Guo et al. [17] on the pile-soil synergistic mechanism of such multipile foundation system revealed a significant influence of long piles relative to the short ones on overall settlement and stiffness of the foundation under vertical working load. In reality, however, a foundation system is not only subject to vertical loading but also to lateral loadings [18][19][20]. Ma et al. [18] conducted experimental and numerical study using ANSYS software to investigate the bearing characteristics of long-short pile composite foundation subjected to horizontal forces while Zhu et al. [20] studied the horizontal bearing characteristics of disconnected piled rafts with equal length piles subjected to lateral loads. It can be seen from the results of these studies that under both vertical and lateral loads, the interposed layer plays a significant role as it adjusts the vertical load sharing while offering horizontal resistance to the lateral load by the friction developing at its interface with the other contacting materials. Nonetheless, studies on horizontal load bearing behaviour are still limited, probably because of the complex soil structure interaction.
In recent years, the safety of conventional pile foundation subjected to ground movements arising from nearby excavations and/or other sources has particularly garnered the attention of scholars globally [21][22][23][24]. is is because the risk associated with foundation pit excavation-induced phenomenon could bring about causalities and huge economic losses as a result of catastrophic collapse of infrastructures [25][26][27]. To avoid such undesirable events, stringent performance requirements are nowadays put as risk control measures for the safety of neighbouring facilities in congested municipalities. Along with this, scholars [28,29] have suggested that lessons taken from previous foundation pit excavation failures could potentially be beneficial in preventing future failure events. According to Poulos [19], the mechanism of loading arising from ground movements is different in action to that of lateral loadings acting on the pile head and even if the ultimate bearing capacity may not be affected, the load from such hidden sources has a significant influence on adjacent piles. In the same manner, the lateral ground disturbances associated with the process of excavation affect the load sharing mechanism of nearby CFG pile composite foundations [30]. e schematic representation in Figure 1 illustrates a typical noncontact long-short piled raft foundation situated at a distance x 0 away from an excavation of depth H e retained by a retaining wall embedded to a depth H p below the dredge line. e common pile to pile spacing in composite foundation usually ranges from 2 to 5 times the pile diameter; and the outward projection of the mat edge from the outer pile is 1.5 to 2 times the pile diameter [31]. Unless the excavation-induced deformation of the support structure is kept under control within a very small range, the resulting soil movement would influence the performance of such foundations in immediate vicinity, especially if they exist within excavation influence zone. Nevertheless, in recent years, many efforts have been made toward development and distribution of lateral earth pressure acting on retaining wall neighbouring composite foundations [32,33].
is is because of the lateral reinforcement effect of CFG piles that cannot be ignored and has to be acknowledged to optimize the rigidity and working mechanism of foundation pit support structures [34]. For instance, when it comes to the active earth pressure acting on the wall, studies have shown that the calculation results obtained from the conventional theories are normally larger than experimental observations [33,35] and numerical findings [36] due the shielding effect of the CFG piles. According to Wei [37], the conventional algorithm overestimates the active earth pressure distribution on the wall in the upper negative friction zone of the CFG piles.
Even though CFG pile composite foundation relatively improves the safety of the foundation pit excavation, the load redistribution and transfer mechanism in adjacent piles and the soil between piles need to be reconsidered to better understand the soil-structure interaction synergy. According to Wang and Yang [38], the load sharing phenomenon in CFG pile composite ground continues to alter throughout the process of neighbouring excavation work. e influence of excavation sequence was also studied by Fu and Li [39] using the centrifuge test, and the result indicated that the later excavation stage highly affected the load sharing characteristics of CFG piles. Physical model test results of Li et al. [40] on the load distribution of rigid-pile composite foundation revealed the same characteristics under the condition of soil displacement arising from retaining wall rotation about the base. Even though such studies highlighted the effects of soil displacement on the load bearing behaviour, they were limited to equal length piles. e current study is aimed to investigate the response of long-short pile composite foundation in immediate vicinity to retaining wall rotating about its base using numerical simulation. A large-scale indoor 1 g model test was carried out to verify the finiteelement model used. 2 Advances in Civil Engineering

Finite-Element
Modelling. e commercial FE software ABAQUS was used to develop the finite-element model. e model comprised of structural (pile, raft, and wall) and geotechnical parts (soil and cushion) that were modelled with the elastic and elastic ideal plastic Mohr-Coulomb constitutive model, respectively. Figure 2 depicts the typical meshed geometry employed in the model. e retaining wall in each case was 10.0 m deep, to avoid boundary effects, and the homogeneously treated soil medium was made to be at least five times as wide as the breadth B r of the raft in the lateral extent and more than three times as deep as the length of the long pile vertically.
Considering computational expenses without compromising accuracy, the mesh density was made finer in the neighbourhood of structural members with gradual change to coarser mesh sizes for farther locations. Each part was discretized using an eight-node linear brick element with reduced-integration and hourglass control (C3D8R). e surrounding vertical boundaries were modelled to restrain lateral displacements while the bottom was set to constrain both horizontal and vertical displacements. e top surface boundary was set to be free in each direction. e analysis was initiated by generating K0 geostatic stress field within the soil mass that produced very small deformation yet in equilibrium with gravity loading. is initial condition was followed by introducing sequentially piles, cushion layer, and raft by establishing surface-to-surface contacts along interacting interfaces. e face of stiffer contacting part was set as a master surface contact candidate while the face of the remaining part was described as a slave surface. Penalty friction formulation was enforced for both tangential and normal behaviour assuming 0.005 m limiting displacement for elastic slip to develop full contact friction, which overcomes overestimation of pile bending moment as suggested by Miao et al. [41].

Verification of the Finite-Element Model.
e FE model employed in this study was compared with the results from indoor experiment conducted using a model box of 1.6 m × 1.6 m × 2.5 m inner dimension, which was sufficient to ignore boundary effects. e apparatus was designed and developed by the Research Institute of Geotechnical Engineering at Zhengzhou University. e model test setup is shown schematically in Figure 3 and comprised of stationary side walls braced with square section RHS frame; movable front wall equipped with bottom roller and slideway; and the rear wall was formed from dismountable ten split-steel plates connected to the frame with high-strength bolts. Displacement control screw rods were symmetrically located at top, middle, and lower level of the right and left edges of the front wall to permit rotation/translation in a similar way used by other researchers such as Yang et al. [42]. Rectangular RHS steel reaction frame was positioned beside sidewalls. e test was conducted on 2 × 2 (2 long-2 short) pile composite foundation. Instrumented hollow aluminium close-ended tubes of length 2100 mm and 1000 mm were used as long and short piles, respectively, considering the length similarity ratio (λ L ) of 3 to account for scale effect. Young's modulus of each pile was 72 GPa and was not simulated for scaling. e piles' outer diameter and thickness were 100 mm and 3 mm, respectively. e roughness of the pile shaft surface was established with knurling treatment and effected an interface angle of friction 27.3°. e spacing of the piles in both model and prototype was three times the pile diameter so that the area replacement ratio remains 0.087.
A plain carbon structural steel plate (quality Q235B) of dimension 600 mm × 600 mm × 50 mm was used as a rigid raft according to the stiffness requirement suggested by Horikoshi and Randolph [43]. Fine dry river sand of Zhengzhou City, China, and standard sand with the grain size distribution presented in Figure 4 were used as foundation soil and a cushion layer, respectively, in accordance with the test requirement of the Chinese specification (GB/T  Advances in Civil Engineering 50783-2012). e scale effect of particle size on ultimate skin friction was also cross checked to ensure the ratio of pile diameter (d) to the average particle size of the soil (D 50 ) was in the range 30-50 in agreement with Fioravante's [44] suggestion. e maximum and minimum dry density of the analogous fine sand was 1.795 g/cm 3 and 1.592 g/cm 3 , respectively. e cushion layer had a thickness of 50 mm. e bearing layer was prepared by pouring and tamping technique with the sand being freely dropped from the height of the model box by a conveyor belt. Once the fill was made to the level of the prescribed pile toe, the pile was positioned vertically and the next layer was filled up to the top surface. After the top sand surface was levelled, soil pressure cells were placed at specified locations before laying the cushion layer. After the bounding steel frame was laid, the cushion layer was placed. en, the loading plate was stationed on top of the cushion. e vertical working load was applied using a hydraulic jack pressed against the top beam of the reaction frame, and a load cell was mounted between the jack and beam. e model ground was given at least a week before the loading process begun in order to balance the stress state of the soil. e loading application was performed according to the Chinese Technical specifications for building foundation inspection (JGJ340-2015), Technical code for testing of building foundation piles (JGJ106-2014) and Technical code for ground treatment of buildings (JGJ 79-2012). e settlement of raft surface during the loading process was obtained from the reading taken by YHD series electronic strain type displacement sensors. All the strain gauges attached to the pile body, displacement transducers, and pressure sensors placed at the top soil surface between piles and on the head and tip of pile were all connected to static digital data acquisition system to collect real-time measurements. Each earth pressure-measuring cell and load sensor were verified for repeatability through indoor calibration procedure preceding the experiment in order to correct calibration coefficients and stress hysteresis effect.
Based on in-house existing research works conducted with the same model box [33,37], the top of the wall was displaced laterally an amount of Δ � 10-15 mm (0.47-0.75% wall height), to attain active limit state by the rotation of the wall about the base. Accordingly, this study used 10 m horizontal movement of the wall top away from the foundation in 11 stages, with 0.5 mm at the beginning and 1 mm for each of the remaining. Consequently, the displacement of the screw rods on the top were controlled manually to obtain 0.5 or 1.0 mm at each stage according to the calibration of screw rode rotation to yield the same. After the data in each stage were stable, the next rotation was carried out. Since settlement occurred as wall rotation advanced, the applied load on the foundation was manually sustained to maintain the hydraulic jack pressure unchanged throughout the process.
In order to compare experimental and numerical results, an FEA was performed on the same geometric configuration but with the depth of the sand twice the long piles' length. Table 1 summarizes material parameters used to back analyse the experiment, while Figure 5 presents the load settlement before wall rotation and the load sharing characteristics during wall rotation. It can be seen that the prediction from the numerical simulation was in agreement with the results from the experiment in terms of both magnitude and trend. us, the behaviour of long-short pile composite foundation can reasonably be predicted following similar numerical simulation procedure for succeeding analyses. However, there is limitation of computational results owing to the commonly applied way of contact formulation in ABAQUS software through the zero-thickness algorithm, the use of reduced-integration to take advantage of computational cost, and employing Mohr-Coulomb constitutive relationship to approximately capture the nonlinear soil behaviour. Despite these

Parametric Studies
e chosen influencing factors in Table 2 were considered based on the literature on load sharing behaviour and pile responses to neighbouring excavations.

Effect of Number of Piles.
e particular position of the pile in a group affects the response of individual piles owing to shadowing and reinforcing effect [47,48]. e pile-soil interaction and load carrying mechanism of internal (central), lateral (edge), and corner piles are influenced by the overlapping zones of shearing resistance, which also garnered attention to use dissimilar piles for practical purposes of reducing differential settlement even in connected pile rafts [49]. Under the influence of soil movement generated because of rotating rigid retaining wall about its base, the presence of piles not only hinders the free-field soil movement to a certain level but also changes the stress state of the soil close to the piles from its already stressed condition.
As depicted in Figure 6, when the composite foundation was subjected to vertical loading of 180 kPa, increasing pile number from 4 to 25 with the same cover area ratio of 0.087 increased the load borne by long and short piles from 18.81 to 29.04% and 12.17 to 15.80%, respectively, before commencing wall rotation. As a consequence of retaining wall rotation, the long piles' load sharing proportion increased from the earlier 18.81 to 29.04% to 22.94 to 30.34% increments at completion of wall rotation, respectively. Similarly, the percentage of load sharing for short piles has increased, but the relative increment due to number of piles was very little at the end of wall rotation, i.e., an increment of ≤1%. It can also be observed that the % load borne by soil between piles during progression of wall movement decreased at a higher rate for the group having four piles as compared to the other cases, indicating the behaviour of predominant pile load sharing proportion with an increasing number of piles under vertical loading results to negligible adjustment during wall rotation. Figure 7 indicates the proportion of load carried by the central, inner, edge, and corner individual piles. Before wall rotation, the corner pile carried the highest pile head load compared to the central pile, which is in tandem with previous studies [12,14,50]. However, at the final rotation of the wall, the pile head load computed for shorter piles close to the wall became relatively larger and increased at a higher rate than the others while the increment for the long pile head load increased with increasing the wall rotation to a certain amount, e.g., 20-25 × 10 −4 rad, and then due to pilepile interaction effect, it afterward started decreasing with the corner pile (pile A) attaining a 4% increment but the edge (pile C) 1% decrement at the end of wall rotation as compared to their initial values.
For short piles, the closer the pile (pile B) is to the wall, the greater the pile is loaded compared to the lateral pile (pile 2), inner pile (pile D), and/or rear pile (pile J). at is, at the end of wall rotation, piles B, 2, D, and J have borne 26, 22.3, 22.6, and 2% more pile head load, respectively. is is mainly attributed to the effect of reduction in load bearing contribution from the soil between piles due to wall rotationinduced soil settlement.
Moreover, the stress on long piles' top is greater than that on the short piles throughout wall rotation process ( Figure 6). Owing to the group effect, the load shared by the front and inner short piles keeps increasing with an increasing number of piles while the edge short piles gradually stabilize their share in a steady manner to accommodate the soil bearing capacity reduction. In other words, this increase in the piles' load sharing can be interpreted as the piles' bearing capacity loss because the piles act as if they are subjected to a new loading condition violating the settlement criteria due to the effect of the induced soil movement. Such "apparent loss of pile capacity" happened due to external soil displacement, and the induced extra pile settlement can be correlated to describe the bearing capacity loss of existing piles in a similar concept as the routinely used settlement criteria for pile failure (settlement due to degradation of shear stiffness rather than "bulb stress"), which however does not necessarily mean it affects the ultimate geotechnical load bearing capacity of piles. erefore, due diligence should be given to avoid or reduce the adverse effect on the structural integrity [19,51].   Advances in Civil Engineering

Effect of Pile Spacing.
For passive piles that are arranged in a group with closer pile spacing, the increase in spacing increases the lateral soil movement-induced forces acting on them due to soil arching phenomenon. Beyond the spacing of 8 times the pile diameter, the group effect and arching would not be felt, and the failure mode of the soil surrounding the piles would be the same as that of a single pile [48]. In composite pile foundation, pile spacing directly reflects one of the basic load sharing parameters, the pile area replacement ratio, m. Apart from changing the spacing, the area replacement ratio can also be varied by increasing the diameter of the pile (or its head) [50,52,53] and/or provision of pile cap [54,55]. From the mechanical characteristics of CFG pile composite foundation, it is well known that the smaller the clearance between neighbouring piles, the larger the cover ratio of piles and the smaller the load sharing proportion of the soil, meanwhile reducing the foundation settlement. For piles of diameter d arranged in a group with a square configuration of spacing s, it is common to describe the area replacement ratio using m � πd 2 /4 s 2 . (1) As predictable from the characteristics of composite foundation under vertical loading, the lower the area replacement ratio, the lower the soil load sharing proportion in the initial stage (see Figure 8). Provoking wall rotation only increases the pile load sharing proportion, keeping the trend almost similar, in an attempt to equilibrate the additional load impinged by soil stress release from imposed soil movement.
Comparing the effect of replacement ratio on the load sharing proportion of soil between piles during the progression of wall movement, it is observed that the percentage load transmitted to the soil beneath the cushion layer reduces at a higher rate for the group with larger pile spacing than the configuration with smaller spacing between piles. What is understandable also is that as the area replacement ratio is increased (and gradually approaches the optimal pile area cover), the more favourably it contributes to the settlement and bearing capacity enhancement of the whole composite ground; consequently, it makes the composite ground stiffer so that its load sharing mechanism is less disturbed. e larger the number of piles that are stiffer than the surrounding soil, the better the performance in reducing the upcoming decrease in bearing capacity from wall movement-induced soil displacements.

Effect of Pile Tip Location.
Assessing earlier research indicates that increasing the ratio of long to short pile length for a given length of long pile increases the pile-soil stress ratio of the short pile while reducing the same for the long piles [56]. On the other hand, the formation level of excavation with respect to the depth of pile tip affects the response of pile to neighbouring excavation [51], and for short piles as compared to long ones, excavation-induced pile settlement is likely to be a significant factor than lateral pile deflections [21]. Akin to the stress release attributable to excavation, wall movement-induced deformation of soil leads to reduction in pile-soil stiffness, which in turn brings pile settlement owing to degradation of pile-soil interface stiffness and changes in pile-soil load sharing ratio. Alongside it, when floating piles are subjected to both axial load on pile head and a surcharge on top of the surrounding soil, the relative pile-soil settlement characteristics and the down-drag load on the piles are affected by the location of pile tip relative to the stable and/or firm stratum [57]. In this section, in order to demonstrate the influence of pile tip location, the length of short piles is varied whereby the pile toe is positioned above, at and below the depth of axis of wall rotation. Figure 9(a) shows the lateral deflection of selected piles. It demonstrates that the horizontal displacement and deflection of short pile positioned above the depth of axis of wall rotation are greater than the computed results for piles located at and below it. In addition, the longer the piles, the more slender they become and the pattern of deflection profiles follows a similar response like the relatively flexible semirigid conventional piles. Moreover, the longer the pile, the larger its lateral head displacement. However, for short piles, the lateral pile toe displacement is far larger than the long pile's tip horizontal displacement, and it is more pronounced for the pile tip located above the depth of axis of wall rotation. At first glance, it is apparent that both the head and tip of short piles undergo a certain amount of translation and rotation while the base of the 20 m long piles encounter restraint boundary condition from the surrounding soil. It is also important to note that rear pile experiences smaller translation and deflection than front pile along the same axis.
At the end of the wall rotation, as can be observed from Figure 9(b), irrespective of the pile tip location, the maximum positive bending moment induced in the front corner pile close to the wall is greater than that of the central and rear corner piles. e axial strains at the outermost nodes of the pile are used to calculate the bending moment distribution, M(z), at any depth z from the pile top according to  10 Advances in Civil Engineering where E p is Young's modulus of the pile, I p is the moment of inertia of the pile, ε 1 /ε 2 are the axial strains along the compression and tension face of pile shaft at the depth z on both sides of neutral axis, and r is the radius of the pile. For convenience, the bending moment is considered positive for the induced tensile stress facing towards the retaining wall. It can also be demonstrated from Figure 9 that letting the length of short piles to be equal to that of the long piles (i.e., both piles assume the same length: 15 m) decreases the maximum bending moment induced in the front peripheral pile of the long-short pile (long 15 m and short 10 m long) placed at the same location. Moreover, owing to the influence of the surrounding reinforcing piles (also affected by distance from wall), the group interaction effect makes the edge and middle piles experience less exposure to disturbance from soil movement and hence the smaller the induced bending moment along them as compared to the peripheral ones but with similar trend. For instance, if the group having 15 m long piles and 10 m short piles is taken, the shielding effect resulted in a 53% and 109% higher bending moment in the front long pile (pile A) as compared to that of the bending moment developed in the central (pile F) and rear corner pile (pile K), respectively, whereas the rear short pile (pile I) sustained a bending moment 53% less than that of the short pile adjacent to the wall (pile C).
Consistent with the difference in displacement of short piles for each case and further pile settlement to mobilize additional pile tip and shaft resistance, the change in pile load sharing ratio of the group of piles with short piles at and above the depth of wall rotation axis is different from that located below it. Figure 10 shows that the rate at which the % load proportion carried by the soil between piles reduces is greater for short piles placed at closer depth to the slip surface (i.e., short piles at and above the depth of axis of wall rotation) than that at a deeper position because of the resultant soil disturbance acting on larger portion of the pile body. Accordingly, the percentage load increment received by long piles increases at higher rate in the composite foundation with short piles of length 6 and 10 m than that of 15 m. Observable from Figure 10 is that the larger the length of piles, the more the load transferred.

Effect of Subsoil Stiffness.
Increase in soil stiffness decreases excavation-induced bending moment and horizontal deformation of nearby piles [58]. It also affects pile's response under vertical loading [53,59]. In addition, it is apparent that the compressibility characteristics of the subsoil depends on its stiffness and affects the pile-soil relative settlement, hence influencing the load sharing characteristics of a composite foundation.
As illustrated in Figure 11, piles embedded in a soil of lower Young's modulus take larger proportion of load and the % load carried by the soil does not drop as rapidly as observed in higher stiffness soils. On the other hand, the stiffer the interpile soil, the smaller the imposed load proportion of the piles; that is, it exhibits an inverse relationship, and as the wall movement progresses, the % load received by the piles increases at a comparatively higher rate until a wall rotation value of 30 × 10 −4 rad and begins to be relatively moderate. e higher the modulus of the subsoil between piles, the larger the load percentage the soil shares prior to wall rotation, and the mobilized bearing capacity reduction in the process of wall rotation increases the concentration of stresses on the CFG piles at a higher rate than the increment for piles installed in soils with lower modulus. e larger the value of subsoil Young's modulus, the smaller the soil displacement due to wall rotation (also explained by [60]) and the larger the stress borne by the soil between the piles. In spite of the small soil displacement in the earlier wall rotation, there is a significant reduction in the load taken by the shallow soil between piles.

Effect of Stiffness and ickness of Cushion.
e variation of load sharing ratio for different cushion thickness and stiffness values during progression of wall rotation is depicted in Figure 12. In the interest of foundation load finetuning ability of the interposed layer, a higher modulus of elasticity of the cushion changes the percentage of the load transferred to piles appreciably in the earlier stage of wall rotation (i.e., up to 15 × 10 −4 rad). Subsequently, in the mentioned wall rotation, the load sharing percentage of the soil decreases at larger rate, and beyond that point both piles and soil attain almost uniform load redistribution except for the case of the 15 MPa cushion. Nevertheless, from the load sharing ratio curves of varying cushion thickness in Figure 12(b), it is evident that there is a higher % load sharing ratio of soil between piles for lower cushion thickness which decreases with increase in thickness until completion of wall rotation. Wall movement (10 -4 rad) Figure 10: Load sharing proportion during wall rotation for short pile tip located above, at, and below the axis of wall rotation.

Effect of Location Away from the Wall.
e pile position in a group and its perpendicular distance from the retaining wall are among the important factors affecting the response of pile groups to imposed soil displacements [45,61]. erefore, it was sought to give insight on the influence of the same in CFG pile composite foundations. Besides, with the increase of retaining wall rotation, the soil behind the retaining wall tends to move, and then the soil between the piles settles so that the upper load is gradually transferred to the pile top. e triggered soil displacement field behind the wall diminishes as the distance from the supporting structure increases, and so it does the soil settlement along the depth. In line with the largest ground displacement near the wall, the expected maximum lateral pile deflection occurs around its head [62]. In this regard, it is common to consider the tolerance to soil displacement of neighbouring structures located within the "Apparent excavation Influence Range (AIR)" [63,64]. Unsurprisingly, Figure 13 shows that the incident of proceeding wall movement affects very little the load sharing proportion of CFG pile composite foundation placed at far distance as compared to those which are closer to the wall. In fact, owing to the cushion's stress diffusive mechanism that produces uneven stress distribution between the pile and soil under vertical loading condition, bending moments are generated at different places of the pile body before wall rotation, but the bending moment value is mostly small. Figure 14 shows the bending moment profile of short and long piles at various distance x 0 away from the retaining wall for a group of 9 piles. Interestingly, it was observed that increasing x 0 from 3.0 m to 15.0 m decreases the maximum bending moment developing both in long and short piles in a similar pattern, except for the cases where curvature deviates for the piles at far distance whose considerable pile body lies outside the theoretical failure zone, such as Rankine's slip surface. In addition, the farther the group is located away from the retaining structure, the farther its centroid shifts away from the back face of the wall and experiences lesser wall rotation-induced soil displacements, and in turn a lesser pile deflection of the group as a whole.
is is consistent with the decrease in soil deformation zone for locations away from the retaining wall so that the decrease in shear strain magnitude will not weaken the majority of the surrounding soil confining the rear piles, whereby the smaller the pile proportion within the disturbed soil mass, the lower the bending moment. Wall movement (10 -4 rad) Figure 11: Influence subsoil stiffness on load sharing ratio of 9 long-short pile composite foundation configuration.
Having in mind the reduction of pile axial load carrying capacity as excavation-induced bending moment increases [58,65], sudden structural failure by bending can be prevented by assessing whether or not the pile's section yield moment is exceeded. Moreover, since CFG piles are usually unreinforced, the possibility of developing sever cracks due to the induced bending moment may cause lateral bearing capacity loss, which would affect the cooperation and integrity of the CFG piles in the group. In the conventional pile group back-analysis using excavation-induced free-filed soil movement, it is recommended to take into account pile moment of inertia deterioration to avoid overconservative structural pile responses corresponding to the pile cracking bending moment for large soil displacements at fully mobilized active lateral soil stresses [66,67]. In this regard, it is often the horizontal component of excavation-induced soil movement and the associated large lateral pile displacement recognized to be more crucial for adjacent concrete piles [68,69]. One can also notice at this point that the further the rear piles are located from the wall, the lesser the soil movement experienced by the piles and consequently the lesser the pile deflection. is puts the piles in a relatively more stable situation as compared to the soil dragging action impinged on the front row piles. e peripheral CFG piles, especially those close to the wall, experience larger bending moment as compared to the middle piles in each row, which may need special attention to confirm whether or not the section's capacity is exceeded under a given working load.

Effect of Applied Working Load.
e intensity of foundation load has been reported to have negligible impact on the lateral reaction of conventional pile to excavation-induced soil movements [58,65] but significantly affects its settlement behaviour [70,71]; hence, owing to the distinct feature of relative pile-soil settlement in composite foundation, observation was sought to assist shed more light on the effect of higher stress level encountered prior to introducing wall rotation-induced soil movement. Figure 15 shows the increase in % load sharing ratio of piles with increase in the applied foundation load preceding wall movement, which is a trend in agreement with the pilesoil stress ratio characteristics of composite foundations. e bigger the applied working load, the larger the load carried by the piles and the smaller the proportion of load borne by the upper soil. As the wall movement proceeds, the reduction of % load share in the soil and the increment in both long and short piles occur following a moderate trend but with fairly higher change for lower working load applied on the raft. It can also be seen from the axial force distribution along the Load sharing ratio (%)  Advances in Civil Engineering pile length that the neutral plane where the maximum normal force happens to exist shifts downwards due to the soil movement-induced relative pile-soil settlement.

Conclusion
In this article, 3D finite element modelling was established based on 1 g large-scale experiment to investigate the influence of retaining wall rotation on the load sharing behaviour of adjacent long-short CFG pile composite foundation. e following conclusions were drawn: (i) e process of wall rotation about the base increased the displacement of the soil near the wall, and the anticipated soil movement decreased with increase in depth and distance away from the wall. e prompted soil settlement associated with progression of wall movement resulted in continuous reduction of the mobilized bearing capacity of the soil between piles. e soil bearing capacity loss for the composite foundation having 9 long and short piles range was 1.4-7.5% of the total imposed load depending on how far the foundation was situated away from the wall.
Consequently, the percentage load sharing of piles increased, ranging between 0.83-4.15 and 0.59-3.36%, respectively, for the long and short piles depending on the location of the foundation behind the wall, as the proportion of the subsequent decrease in stress borne by the soil between piles transferred to the piles. (ii) During wall rotation, the 3D soil movement behind the wall increased the vertical and lateral soil strain near the wall. As the soil elements near the wall were subjected to larger principal strain increments with near-lateral extension, it led to larger stress relief experienced by front piles in comparison to the piles far away from the wall. On the other hand, the lateral displacement experienced by the soil between piles with higher elastic modulus was smaller than the soil with low stiffness but as it bore higher load proportion prior to wall rotation, significant reduction in the load taken by the soil was observed more in the earlier wall rotation, i.e., up to 15 × 10 −4 rad. For the composite foundation having 9 long and short piles placed at 3.0 m behind the wall, the soil bearing capacity drop observed was 10.77 and 2.73% of the imposed load for the subsoil elastic modulus of 18.5 and 5 MPa, respectively. (iii) Depending on the pile tip location with respect to the failure surface, the pile lateral response varied. e proportion of the pile body intersected by the slip line and its embedded depth in the stable layer influenced the reaction of the pile to the soil movement.
(iv) During wall rotation progression, the increment in proportion of the load taken by each pile varied depending on the position of the piles in the group and the adjacent wall. e contribution of short piles increased during the course of wall rotation. For the composite foundation having 25 long and short piles placed at 3 m behind the wall, the front long pile attained 4% pile head increment while the short pile bore as much as 26% more head load as compared to the 2% rear short pile's head load increment. Peripheral piles and those close to the wall received the highest load.
(v) For each key factor chosen, long piles' load sharing proportion was bigger than that of the short piles. e larger the area replacement ratio, the more it improved the vertical bearing capacity and the lesser the disturbance in percentage of load shared by the piles and the soil during wall rotation compared to the group with smaller area replacement ratio.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.