Calculation model of the equivalent spiral shear stress of conditioned sand

Soil is the medium that balances the soil and water pressure on the excavation surface of an earth pressure balance (EPB) shield machine, and the soil flow plasticity directly affects the working performance of the shield machine. The spiral shear stress is the main determinant of the soil flow plasticity. During the progress of the EPB shield screw conveyor normal operation, the shear effects between the muck in the spiral tube and the screw shaft, as well as the flight, are different from the shear action, such as the direct shear or the vane shear. The direct or vane shear test is now widely used to evaluate soil flow plasticity. To better investigate the shear properties between the soil and the shield conveyor casing, a model screw conveyor was developed. Based on the analysis of the stress conditions and the movement of the conditioned soil plug, the theoretical calculation model of the equivalent helical shear stress of the EPB shield screw conveyor was deduced based on the steady-state equilibrium conditions of moment of momentum. The conditioned standard sand was used in the indoor tests. The results of the indoor test of the model machine verified the rationality and effectiveness of the theoretical model. The model provides an important theoretical basis and references for the design and optimization of soil conditioning during shield tunneling.


Introduction
Research on the static and dynamic characteristics of the interface between muck and an earth pressure balance (EPB) shield screw conveyor structure is concerned with frontier issues, such as nonlinearity, large creep and discontinuity. It is currently one of the key research subjects related to the physical and mechanical properties of conditioned soil. The model machine test is a basic method that studies the mechanical properties of the interface between muck and the EPB structure. For example, the main factors affecting the mechanical properties of the contact surface and basic laws have been obtained via experimental studies, and a constitutive model describing the mechanical properties of the contact surface has been established by Zhang and Min [1].
Under the influence of a dynamic load from the working chamber pressure and rotation of the shield screw conveyor, the conditioned soil in the spiral groove is subjected to additional stress in the pure shear state and in the spiral rotation of the principal stress axis. The direct or vane shear test can't effectively model this type of stress environment. In this study, the shield a1111111111 a1111111111 a1111111111 a1111111111 a1111111111 screw model machine simulates the stress state of the muck during normal shield tunneling and then the shear stress and creep characteristics between the conditioned soil and the spiral structure in the construction environment are analyzed to guide and solve soil conditioning during the construction process. The screw conveyor model machine performs the spiral shear test on the conditioned soil under the action of the soil pressure and the rotating spiral thrust, where the loading condition is the same as that of the muck in the spiral groove of the construction site.
In light of the importance of muck properties for spewing control of the shield screw conveyors and regulation of the working chamber pressure, the addition of conditioning agents for soil improvement is an effective technical means that has been confirmed by many engineering practices. There are three main types of commonly used additives: foaming agents, bentonite slurry and macromolecular polymers. Scholars have conducted extensive research on the mechanical properties of improved muck.
The slump test, which determines the workability of concrete, is widely used to evaluate the flow characteristics of improved muck. Many researchers have reported results in this area [2][3][4][5][6][7][8]. The slump test is simple and economical, and the test results generally reflect the flow characteristics of the muck. Williamson et al. [3]and Jancsecz et al. [4] studied the improvement effect of foam on sand, suggesting a reasonable slump range of 200-250 mm. The improvement effect of pebble stratum with a foam and bentonite slurry was indicated to be between 150-200mm by Jiang et al. [9]. Gharahbagh et al. [10] combined soil conditioning with the tool wear test, suggesting a reasonable range of slump values of 100-250mm.
From the research results of the above studies, we can see that the slump used to measure the flow characteristics of muck is an important indicator of the conditioning effect. However, shear stress is the main influencing factor of the slump of conditioned soil. Therefore, it is necessary to study the calculation method of the spiral shear stress of conditioned soil. The shield screw conveyor is installed between the bottom of the front shield and the erector at the tail of the shield, wherein the midline rises at a certain inclination angle. When the screw conveyor is normally drained, the shaft and the spiral flight, which extend into the working chamber, are rotated by a hydraulic motor; moreover, the muck is transported and raised along a certain spiral angle under the joint action of the spiral and the outer shell and is discharged to the soil outlet. The EPB shield screw conveyor is shown in Fig 1. The helical shear stress occurs The spiral shear stress between the muck and the spiral structure. The friction and wear problems between the muck and the spiral structure and the particularity of the spiral structure make it difficult to directly measure the helical shear stress.
In recent years, direct shears (Ehsan and Adam [20]; Psomas [21]; Houlsby and Psomas [22]) or modified vane shears (Zumstegt et al. [23]; Messerklinger et al. [24]; Zumsteg et al. [25]; Mori et al. [26]) have been used to study the effects of the loading rate and confining pressure on the mechanical properties of conditioned soils, and specially designed rheometers have been used to test the characteristic parameters (such as yield stress and viscosity) (Karmakar [27]; Meng et al. [28]) of muck. In general, the results of this type of test method can simulate the influence of the pressure environment of the working chamber on the strength of muck to some extent and can characterize the viscoplastic parameters of muck.
The above results have advanced the understanding of the mechanical properties of conditioned soil. To more accurately simulate the mechanical environment and stress state of the muck in the screw conveyor of the EPB shield, this paper developed a shield screw conveyor model. By analyzing the forces and motion of the soil microbody according to the steady-state equilibrium conditions of moment of momentum, the equivalent spiral shear stress calculation model of the EPB shield screw conveyor is derived. The test results confirm the validity and rationality of the model.

Calculation model for the equivalent helical shear stress of conditioned soil
For the convenience of analysis, the symbols used in the paper are listed in S1 Table.

Description of the spiral structure
The parameters of the spiral blade and screw shaft of an EPB shield screw conveyor are shown in Fig 2. The angle between the tangent of any point on the spiral blade and its horizontal projection is the helix angle. Unfolding the spiral of a pitch, the spiral is the oblique side of the rightangled triangle, the pitch is the cathetus of the right triangle, and the spiral is another cathetuswhen it is projected in the horizontal direction. The unfolding result is shown in Fig 3. From this figure, the helix angle φ of any point on the spiral blade can be derived.
The spiral angle of the outer edge of the spiral blade is defined by Eq (2), and the definition of the helix angle at the spiral axis is given by Eq (3).

Stress and motion analysis of the soil microbody
It is assumed that the conditioned sandy soil is a homogeneous fluid body and is isotropic. As shown in Fig 4, the soil microbody in the spiral groove is taken as the research object. The size of the microbody is defined in Fig 5. Eqs (4) to (10) are derived from the relationships among the sizes of the microbody.
During the conveying of the muck by the shield screw conveyor, the soil in the spiral groove is subjected to seven forces: resultants of shear stresses (shear stress of the upper spiral blade (τ f ) t , shear stress of the lower spiral blade (τ f ) b , shear stress of the helical axis τ s , and shear stress of the spiral shell τ c ), the combined pressing force Q n of the upper and lower spiral blades, the pressure P in the direction of the spiral groove and the weight G of the soil microbody. Fig 6 shows the forces acted on the soil microbody. The two-axis system is defined in Fig 6; the (x, y) axis is perpendicular and parallel to the screw shaft, and the (l,u) axis is parallel and perpendicular to the screw blades.
During the normal tunneling process of the shield, due to the relative movement between muck and both the spiral groove and the spiral shell, friction shear stress is generated between the soil in the spiral groove and the spiral casing, as well as the spiral shaft and the spiral blade. Under the action of the helical axis and the shear stress of the spiral blade, the soil microbody has a linear velocity v m following the rotation of the spiral blade v s , and the angle between the velocity v m and the direction of the vertical spiral axis is θ, as shown in  The rotational speed component of the spiral blade tangent to the spiral casing is The speed of the spiral blade along the axis is Eqs (13) and (14) can be derived from   Using Eqs (12) and (14), the conveying efficiency η of the screw conveyor is Referring to the geometry of the screw conveyor, the maximum volume delivery in theory is The actual delivery volume is The following equations can be derived from Eqs (15) and (17) As a result, the calculation formula of the flow direction of the muck can be obtained as follows:

Derivation of the equivalent spiral shear stress formula
When the shield machine is smoothly boring through soil, the screw conveyor quantitatively discharges the conditioned soil. As the muck is spirally moved upward in the spiral groove, a reverse helical drive torque is generated under the action of the force in the direction of the opposite spiral rotation. The equivalent spiral shear stress calculation model can be derived from the equilibrium of moment of momentum in the direction of the vertical spiral axis. It can be seen from the motion and the forces of the research object (Figs 5-7) that only the friction force between the muck and the casing and the weight of the soil have the components in the opposite direction of the rotation. According to the momentum equilibrium condition of the research object along the vertical spiral axis direction, the following formula can be written: The negative sign on the right side of Eq (20) indicates that the weight of the muck's microbody and the shear stress of the spiral shell perpendicular to the direction of the spiral axis are opposite to the direction of the rotational torque. Substituting Eqs (4), (5), (8) and (9) into Eq (20), the following equation can be obtained (to simplify the analysis, the thickness e of the spiral blade is ignored): Ignoring the difference between the helix angle φ of the outer edge of the spiral blade and the average helix angle φ m of the spiral blade and using Eq (2), Eq (21) can be further simplified to Eq (22) is integrated along the direction of the helix axis, the length of which is L. Assuming the muck filling rate is to be 100% (this is in line with reality), the equivalent spiral shear stress formula is available According to Eq (23), the calculated shear stress is affected by the helical geometric parameters, the density of the conditioned muck, the angle of the screw installation, and the parameters of the operating state. Through this formula, the shear stress at the interface between the conditioned soil and the conveyor casing can be calculated for given torques and flow directions of the muck. When θ is 0˚, the screw conveyor is clogged, and the muck does not move in the direction of the screw axis. The shear stress when the muck is clogged can be calculated from Eq (23); thus, the quantitative analysis and suggested calculation parameters can be provided for preventing the blockage of the shield screw conveyor. Under the premise of determining the geometry and installation inclination angle of the shield screw conveyor machine, the shear stress of the conditioned soil is proportional to the screw machine driving torque and the flow angle of the muck. In other words, the more rigid the muck is, the greater the required driving torque and the higher the shear stress at the same soil flow angle.
It is seen that the muck density is of great importance for the stress created at the interface between the conveyor casing and the conditioned soil. Then how to get the desired density is to be concerned. From the difference of the density between the soil and the conditioned soil (the muck), the conditioned muck relates with the specialized additives and the injection amount of the specialized additives, both of which can be decided quantitatively in the laboratory.

Indoor model machine test
To test the validity of the model for calculating the equivalent shear stress of the conditioned soil, a series of muck conveyance tests were conducted on the shield screw conveyor model machine, and the difference between the measured values of the model test and the theoretical calculation results was compared and analyzed. The influence of different factors on the shear stress of conditioned soil was discussed.

Shield screw conveyor model machine
The indoor model test device is shown in Fig 9. The model machine system consists of a soil sample box, a tilted spiral tube with a screw inside, a coupling, a torque sensor, a geared motor, a frequency converter and a horizontally applied pressure system in the soil sample box. Load cells, stress sensors and pore sensors are installed outside the spiral tube for measuring the normal stress, the pore water pressure, and the shear stress in the direction of the parallel and vertical spiral axes. Motor drive torque and speed are collected by the torque sensor. According to the soil transportation situation in the soil box, the muck transportation rate is calculated. The specific technical parameters of each component are shown in Table 1. The conveyor model was modeled after the screw conveyors installed in the shield machines widely used in Beijing tunnel construction. For convenience of easy operation and manufacture, the minification ratio of 1:8 was adopted, resulting in the parameters listed in Table 1.

Experiment material
Standard sand was selected for the test, which was manufactured by Xiamen Iso Standard Sand Co., Ltd. After screening with a particle size of 2 mm, the particle sizes were between 0.075 and 2.0 mm, the silica (SiO 2 ) content was greater than 96%, and the mud content (including soluble salts) was not more than 0.2%. The physical and mechanical indicators are shown in Table 2, which are provided by the supplier. A series of slump tests were conducted using this standard sand to determine the workability of the muck and to determine reasonable additives and ratios.
The foaming agent solution can enhance the fluidity of the muck, and the high-quality bentonite slurry has a good cohesive force. When too much slurry is added, the fluidity of the sand will also increase. According to this feature, a sodium-based bentonite slurry with a mass ratio of 1:5 is selected, with a slurry volume injection ratio of 15%. At this time, the sand has a  The spiral shear stress certain adhesive force, but it is still loose and brittle, and the fluidity is poor. Subsequently, different amounts of foam with a concentration of 7% are added, increasing the blending ratio by 10% each time until the desired fluidity effect is achieved (the optimal slump of sandy soil is in the range of 100-200mm [29]).The test results of standard sand conditioned by bentonite slurry and foam with different injection amount are shown in Fig 10. It can be seen from the test results in Fig 10 that, when only 15% bentonite slurry is added, the sand is in a relatively loose state. The soil sample can reach a relatively good flow state after adding 10% foam. The variation of the conditioned sandy soil slump with different foam solution injection ratios is shown in Fig 11. Based on the above test results, the soil conditioning parameters used in the model machine test are a concentration of 7% foaming agent solution with a blending ratio of 10% and a sodium-based bentonite slurry with a water-soil quality of 5:1 and a blending ratio of 15%. The conditioned muck density was 1452 kg/m 3 .

Test results
The helical shear stress of each measuring point is synthesized by the horizontal component and the vertical component, and is calculated according to Eq (24).  The spiral shear stress The laboratory test was repeated three times under the same test conditions. The data measured below are the average of the three measurements. After the working state of the shield screw conveyor model machine is stable, the volume transportation rate of the conditioned soil becomes constant. According to the volume of the soil sample in the sand sample box and the volume of the muck discharged from the outlet, the average transporting rate of the muck is calculated to be 0.18 (m 3 /h); thus, the average transport efficiency of the model machine is 79%. The driving torque is shown in Fig 12. After the model machine runs for approximately 30 seconds, the driving torque enters the steady-state range. The average torque is 40.95 N � m. According to Eq (19), the flow angle of the muck is calculated to be 51˚. From Eq (23), the helical shear stress at the interface between the modified soil and the conveyor casing can be calculated to be 4.15 KPa.
Twelve shear stress measuring points are arranged on the conveyor casing, and the shear stress of each measuring point changes with time, as shown in Figs 13-18.
Based on Figs 13-18, which show the shear stress-time history curves of the measured points, the following results can be summarized: 1. The shear stress of measuring point 1 is slightly larger than the value of measuring point 2 because measuring point 1 is close to the soil sample box; thus, the compressive stress is higher than that of measuring point 2, and under the action of the compressive stress, the shear stress will increase slightly. When the model screw conveyor is working normally, the linear fitting residual of the shear stress at each measuring point indicates that the shear stress of the spiral casing is approximately constant, as shown in Fig 19. The mechanism of the change of spiral shear stress with compressive stress will be the next problem studied.  The spiral shear stress From the comparison results in Fig 20, the mean shear stress of each measuring point is close to the theoretical calculation result, that is, 4.15 KPa. It is shown that the equivalent spiral shear stress calculation model derived from the equilibrium condition of moment of momentum can accurately predict the shear behavior of the conditioned soil. Moreover, the calculation results are related to the physical and mechanical parameters of the conditioned muck, the operating conditions, and the geometric structure of the shield screw conveyor.

Conclusions
It is seen that the driving torque of a screw conveyor is closely related with the shear stress created at the interface between the conveyor casing and the conditioned soil. The driving torque is easy to measure, whereas the shear stress is difficult to get due to its occurrence especially in the conveyor casing of a real shield machine. From the screw conveyor model machine tests, the shear stress along the casing is measured. The measured results coincide well the calculated results by Eq (23), which also validate the derived relationship between the driving torque and the shear stress.