Method and system for calculating pile driving and soil squeezing effect of shaped pile considering stress anisotropy

Abstract

The present application relates to the technical field of pile foundation engineering and geotechnical mechanics, in particular to a calculation method and system for the pile sinking and soil squeezing effect of a special-shaped pile considering stress anisotropy. The method comprises: obtaining relevant parameters of the special-shaped pile, and setting boundary stress and assumption conditions; constructing a reaming model according to the relevant parameters, the boundary stress and the assumption conditions; calculating the like-plane radial velocity and the like-plane hoop velocity based on the fluid mechanics theory and by using the pre-obtained conformal transformation function, the reaming model and the assumption conditions; calculating the like-plane soil body displacement and the like-plane stress distribution based on the reaming model, the like-plane radial velocity, the like-plane hoop velocity, the conformal transformation function, the assumption conditions and the pre-obtained K0 anisotropy correction Cambridge constitutive model; and calculating the stress field and the displacement field of the soil body around the special-shaped pile by using the assumption conditions, the like-plane soil body displacement, the like-plane stress distribution and the conformal transformation function. The present application solves the problem of inaccurate calculation of the pile sinking and soil squeezing effect of a special-shaped pile in the prior art.

Description

Technical Field

[0001] This invention relates to the fields of pile foundation engineering and geotechnical mechanics, specifically to a method and system for calculating the soil squeezing effect of irregularly shaped piles under stress anisotropy. Background Technology

[0002] Major infrastructure projects are placing increasingly higher demands on pile foundation performance. Reducing material consumption, minimizing construction energy consumption, and improving pile foundation performance have become bottlenecks restricting the development of pile foundation technology. In recent years, various new types of irregularly shaped piles (such as rectangular piles, X-shaped piles, and other irregularly shaped cross-section piles) have been continuously developed and applied. Due to their special cross-sectional design, they significantly improve the bearing characteristics of pile foundations and achieve a substantial increase in bearing capacity. They possess numerous advantages such as "high bearing capacity," "good economy," and "high construction efficiency," and are widely used in major infrastructure projects such as highways and high-speed railways.

[0003] To date, although irregularly shaped pile foundations have a very broad development and application prospect, their design and calculation system still lags behind engineering practice. In particular, research on the soil squeezing effect of irregularly shaped pile driving, considering stress anisotropy, is still not comprehensive or in-depth, and a mature and reliable design and calculation method is lacking. During the driving of irregularly shaped piles, the initial stress of the natural soil often exhibits anisotropic characteristics due to sedimentation and consolidation environments, leading to significant changes in the displacement and stress fields of the surrounding soil. This, in turn, affects the bearing capacity of the pile foundation itself and may even have adverse effects on adjacent buildings and structures. However, existing methods for calculating disturbances during irregularly shaped pile driving all assume that the soil stress is isotropic or do not comprehensively consider stress anisotropy, resulting in inaccurate calculations of the soil squeezing effect of irregularly shaped pile driving and limitations in engineering applications.

[0004] Therefore, there is an urgent need for a method to calculate the soil squeezing effect of irregularly shaped piles that take into account stress anisotropy. Summary of the Invention

[0005] To address the shortcomings of existing technologies, this invention provides a method and system for calculating the soil squeezing effect of irregularly shaped piles during pile driving that takes into account stress anisotropy, thus solving the problem that the calculation of the soil squeezing effect of irregularly shaped piles during pile driving is not accurate enough in existing technologies.

[0006] To achieve the above objectives, one aspect of the present invention provides a method for calculating the soil squeezing effect of irregularly shaped piles under stress anisotropy. The method includes: acquiring relevant parameters of the irregularly shaped pile and setting boundary stresses and assumptions; constructing a borehole expansion model based on the relevant parameters, the boundary stresses, and the assumptions; calculating the radial velocity and circumferential velocity of the image plane using a pre-acquired conformal transformation function, the borehole expansion model, and the assumptions based on fluid mechanics theory; calculating the soil displacement and stress distribution of the image plane based on the borehole expansion model, the radial velocity, the circumferential velocity, the conformal transformation function, the assumptions, and a pre-acquired K0 anisotropic modified Cambridge constitutive model; and calculating the stress field and displacement field of the soil under the irregularly shaped pile under driving using the assumptions, the soil displacement, the stress distribution, and the conformal transformation function.

[0007] This invention acquires parameters of irregularly shaped piles and sets boundary conditions to build a hole enlargement model, accurately reproducing the physical scenario of pile driving and hole enlargement for irregularly shaped piles. It avoids the distortion problem of adapting a circular hole enlargement model to irregular cross-sections. Combining fluid mechanics and conformal transformation functions to solve the image plane velocity, it maps irregular irregular holes to a standard unit circle to simplify the solution difficulty. Relying on the K0 anisotropic correction Cambridge constitutive model, it solves the soil displacement and stress in the image plane, fully restoring the true mechanical properties of the anisotropic initial stress of natural soil. Then, by combining conformal transformation and reverse conversion, it obtains the actual soil stress and displacement field, significantly reducing the calculation error of the soil squeezing effect of irregularly shaped piles and improving the calculation accuracy of the soil squeezing effect of pile driving for irregularly shaped piles.

[0008] Optionally, the calculation of the radial velocity and circumferential velocity of the image plane based on fluid dynamics theory using a pre-obtained conformal transformation function, the borehole enlargement model, and the assumed conditions includes: determining the Laplace equation for the velocity of the soil around the borehole based on fluid dynamics theory using the assumed conditions; and calculating the radial velocity and circumferential velocity of the image plane based on the Laplace equation for the velocity of the soil around the borehole, the pre-obtained conformal transformation function, and the borehole enlargement model.

[0009] This invention derives the Laplace equation for soil velocity based on fluid mechanics theory and presupposed conditions, accurately characterizing the irrotational and incompressible motion of soil during pile driving. Then, it solves the radial and circumferential velocities in the image plane by simultaneously applying the Laplace equation, conformal transformation function, and borehole enlargement model. This transforms the physical plane irregular hole into a regular unit circle in the image plane, eliminating the obstacle of solving irregular cross-sections, simplifying the solution process of complex partial differential equations, and improving the efficiency and computational stability of solving the velocity field of soil around irregular pile holes.

[0010] Optionally, the calculation of soil displacement and stress distribution in the image plane based on the borehole enlargement model, the radial velocity of the image plane, the circumferential velocity of the image plane, the conformal transformation function, the assumptions, and the pre-acquired K0 anisotropic modified Cambridge constitutive model includes: calculating the complex velocity variables of the physical plane and the complex strain rate variables of the image plane based on the borehole enlargement model, the radial velocity of the image plane, the circumferential velocity of the image plane, the conformal transformation function, and the assumptions; and calculating soil displacement and stress distribution in the image plane using the complex velocity variables of the physical plane, the complex strain rate variables of the image plane, and the pre-acquired K0 anisotropic modified Cambridge constitutive model.

[0011] This invention solves for the complex variables of physical plane velocity and image plane strain rate by combining the borehole expansion model, image plane velocity and conformal transformation function. The complex variable function uniformly represents the global velocity and strain information, simplifies the coupled calculation of multi-component mechanical quantities under irregular cross sections, and then combines the K0 anisotropic modified Cambridge constitutive model to solve for soil displacement and stress. It fully considers the anisotropic elastic-plastic deformation characteristics of natural soil under initial stress, and improves the calculation fit and accuracy of image plane soil displacement and stress distribution.

[0012] Optionally, the calculation of the complex velocity variables of the physical plane and the complex strain rate variables of the image plane based on the aperture expansion model, the radial velocity of the image plane, the circumferential velocity of the image plane, the conformal transformation function, and the assumptions includes: calculating the radial velocity of the physical plane using the conformal transformation function and the aperture expansion model; calculating the radial velocity on the unit circle of the image plane using the conformal transformation function, the radial velocity of the physical plane, and the assumptions; calculating the complex velocity variables of the physical plane based on the radial velocity on the unit circle of the image plane, the radial velocity of the image plane, the circumferential velocity of the image plane, and the conformal transformation function; and calculating the complex strain rate variables of the image plane using the complex velocity variables of the physical plane based on complex variable function theory.

[0013] This invention solves for the radial velocity of the physical plane by relying on conformal transformation functions and hole expansion models, accurately matching the expansion motion characteristics of irregular hole boundaries. Then, it combines transformation functions and the radial velocity of the physical plane to solve for the radial velocity of the unit circle in the image plane, completing the equivalent transformation from the irregular boundary to the velocity boundary of the regular unit circle. By combining the radial and circumferential velocities of the unit circle with the transformation function, the complex variable of the physical plane velocity is obtained, and the velocity components in the plane are unified and integrated. Finally, the complex variable of strain rate is derived from the complex variable of velocity using complex variable function theory, simplifying the solution process of multi-dimensional strain component separation and improving the adaptability and calculation accuracy of the velocity and strain field conversion of irregular pile soil.

[0014] Optionally, the complex velocity variable of the physical plane satisfies the following formula: , in, Let the velocity of the physical plane be a complex variable. For complex variables in the image plane The conjugate of complex numbers, The complex conjugate of the first derivative of the conformal transformation function. The coefficients in the potential function obtained using the method of separation of variables, For simplification symbols.

[0015] This invention constructs a complex variable expression for the velocity in the physical plane, establishes a velocity mapping relationship from the image plane to the physical plane based on the complex variables, constant terms, and series coefficients related to conformal transformation, and adopts an infinite series form to adapt to the boundary contour of any irregular pile section, uniformly representing the velocity components in the plane, avoiding the problem that the velocity field is difficult to solve directly under irregular boundary conditions, and improving the versatility and calculation accuracy of the soil velocity field solution during the pile driving process of irregular piles.

[0016] Optionally, the complex strain rate variable of the image plane satisfies the following formula: , , in, Let the strain rate be a complex variable on the image plane. Let the velocity of the physical plane be a complex variable. For complex variables in the image plane The conjugate of complex numbers, For complex variables in the image plane, The complex conjugate of the first derivative of the conformal transformation function. The complex conjugate of the second derivative of the conformal transformation function. The coefficients in the potential function obtained using the method of separation of variables, For simplification symbols.

[0017] This invention establishes an explicit calculation formula for the complex variable of the image plane strain rate, and completes the complex variable differentiation operation based on the complex variable of the physical plane velocity and the derivatives of each order of conformal transformation. It directly realizes the transformation from the velocity field to the strain rate field without having to solve multiple strain components separately. It unifies the complex variable form to adapt to the conformal transformation mapping relationship of irregular cross sections, avoids the problem of difficult decoupling of strain components under irregular boundaries, and improves the calculation efficiency and solution stability of the strain rate field of the soil around the irregular pile.

[0018] Optionally, the calculation of soil displacement and stress distribution on the image plane using the complex velocity variable of the physical plane, the complex strain rate variable of the image plane, and the pre-acquired K0 anisotropic modified Cambridge constitutive model includes: determining the kinematic equation of the soil on the image plane using the complex strain rate variable of the image plane and the complex velocity variable of the physical plane; calculating the elastoplastic constitutive matrix of the soil on the image plane using the yield function of the pre-acquired K0 anisotropic modified Cambridge constitutive model; obtaining the constitutive equation of the image plane based on the elastoplastic constitutive matrix of the soil on the image plane; and calculating the soil displacement and stress distribution on the image plane using the kinematic equation and the constitutive equation of the image plane.

[0019] This invention establishes the image plane soil kinematic equation by combining the complex variables of image plane strain rate and physical plane velocity, thus fully characterizing the soil motion and deformation relationship. Based on Hooke's law, superposition principle, related flow law and hardening law, and combined with the yield function of the K0 anisotropic modified Cambridge constitutive model, the elastoplastic constitutive matrix is ​​solved, and the image plane constitutive equation is constructed. The kinematic equation and constitutive equation are solved simultaneously to solve the soil displacement and stress distribution. It fully considers the anisotropy of the initial stress and the elastoplastic deformation characteristics of the soil, and improves the consistency between the image plane soil displacement and stress calculation results and the actual soil mechanical response under the condition of irregular pile driving.

[0020] Optionally, calculating the stress field and displacement field of the soil body for the irregularly shaped pile driving using the assumed conditions, the image plane soil displacement, the image plane stress distribution, and the conformal transformation function includes: calculating the physical plane soil displacement and physical plane stress distribution using the conformal transformation function, the image plane soil displacement, and the image plane stress distribution; calculating the physical plane pore water pressure based on the stress balance equation, using the assumed conditions, the conformal transformation function, and the image plane stress distribution; and calculating the stress field and displacement field of the soil body for the irregularly shaped pile driving based on the physical plane pore water pressure, the physical plane soil displacement, and the physical plane stress distribution.

[0021] This invention combines conformal transformation functions and image plane soil displacement and stress distribution conversion to obtain physical plane soil displacement and stress distribution, realizing the mapping of image plane solution results to the real physical region. Relying on stress balance equations, combined with assumptions, conformal transformation functions, and image plane stress distribution, the physical plane pore water pressure is solved, improving the saturated soil mechanical variable system. The final field quantity results are obtained by integrating physical plane pore water pressure, soil displacement, and stress distribution, completely restoring the coupling effect between soil skeleton and pore water during pile driving, and improving the completeness and accuracy of predicting soil stress field and displacement field for irregular pile driving.

[0022] Optionally, the physical plane pore water pressure includes: , , , , in, For physical plane soil pore water pressure, For example, the pore water pressure of a planar soil mass, The initial pore water pressure, Let the radial radius of the image plane be . For the effective radial stress in the image plane, For plane shear stress, The polar angle of the image plane, For the effective circumferential stress in the image plane, These are the coefficients in the conformal transformation function that control the size of the hole. The first derivative of the conformal transformation function, The complex conjugate of the first derivative of the conformal transformation function. Let be the real part of the complex number. The second derivative of the conformal transformation function, The imaginary unit, Let be a point on the unit circle of the image plane. , and All are simplified symbols.

[0023] This invention constructs an integral calculation formula for pore water pressure in the physical plane, defines intermediate simplification coefficients based on the derivatives of conformal transformation, and establishes a quantitative mapping relationship between effective stress and shear stress in the image plane and pore water pressure in the physical plane. This formula can adapt to the irregular cross-sectional boundaries of irregular piles, fully consider the combined contribution of radial, circumferential, and shear stresses to pore water pressure, and accurately convert the image plane solution results to the real physical plane, thereby improving the completeness and accuracy of pore water pressure calculation in saturated soil during the pile driving process of irregular piles.

[0024] Another aspect of the present invention provides a calculation system for the soil squeezing effect of irregularly shaped pile driving considering stress anisotropy, comprising: a processor, an input device, an output device, and a memory, wherein the processor, the input device, the output device, and the memory are interconnected, wherein the memory is used to store a computer program, the computer program including program instructions, and the processor is configured to call the program instructions to execute the calculation method for the soil squeezing effect of irregularly shaped pile driving considering stress anisotropy as described in any of the preceding aspects of the present invention.

[0025] The soil squeezing effect calculation system for irregularly shaped pile driving considering stress anisotropy of the present invention has a compact structure, stable performance, high integration and simple composition. It can stably execute the soil squeezing effect calculation method for irregularly shaped pile driving considering stress anisotropy provided in the first aspect of the present invention, further improving the overall applicability and practical application capability of the present invention. Attached Figure Description

[0026] Figure 1 This is a flowchart of the method for calculating the soil squeezing effect of irregularly shaped piles considering stress anisotropy in an embodiment of the present invention. Figure 2 This is a schematic diagram of a soil borehole enlargement model for driving irregularly shaped piles considering stress anisotropy, as described in an embodiment of the present invention. Figure 3 This is a schematic diagram of the physical plane and image plane in an embodiment of the present invention; Figure 4 This refers to the normalized displacement of the soil mass during the driving of irregularly shaped piles in this embodiment of the invention. Figure 5 The normalized effective stress of the soil around the pile hole of the irregularly shaped pile in this embodiment of the invention; Figure 6 This is a normalized shear stress distribution cloud map of the soil mass during the driving of irregularly shaped piles, according to an embodiment of the present invention. Figure 7 This is a normalized pore water pressure distribution cloud map of the soil mass during the driving of irregularly shaped piles, according to an embodiment of the present invention. Figure 8 This is a schematic diagram of the structure of the calculation system for the soil squeezing effect of irregularly shaped piles considering stress anisotropy, according to an embodiment of the present invention. Detailed Implementation

[0027] Specific embodiments of the present invention will now be described in detail. It should be noted that the embodiments described herein are for illustrative purposes only and are not intended to limit the invention. In the following description, numerous specific details are set forth in order to provide a thorough understanding of the invention. However, it will be apparent to those skilled in the art that these specific details are not necessary to practice the invention. In other instances, well-known circuits, software, or methods have not been specifically described to avoid obscuring the invention.

[0028] Throughout this specification, references to "an embodiment," "an embodiment," "an example," or "an example" mean that a particular feature, structure, or characteristic described in connection with that embodiment or example is included in at least one embodiment of the invention. Therefore, the phrases "in an embodiment," "in an embodiment," "an example," or "an example" appearing in various places throughout the specification do not necessarily refer to the same embodiment or example. Furthermore, specific features, structures, or characteristics can be combined in one or more embodiments or examples in any suitable combination and / or sub-combination. Moreover, those skilled in the art will understand that the illustrations provided herein are for illustrative purposes and are not necessarily drawn to scale.

[0029] Please see Figure 1 In order to solve the problems in the prior art, in an alternative embodiment, such as Figure 1 The method for calculating the soil displacement effect of irregularly shaped piles considering stress anisotropy, as shown, includes the following steps: Step S1: Obtain the relevant parameters of the irregular pile and set the boundary stress and assumptions.

[0030] In this embodiment, the irregular pile can be any one of a rectangular pile, an X-shaped pile, or an elliptical pile. The relevant parameters of the irregular pile include geometric parameters, soil initial stress parameters, and soil material parameters.

[0031] Geometric parameters include the first set of parameters Second set parameters When the irregularly shaped pile is a rectangular pile, the longer side is... The shorter side is When the irregularly shaped pile is an X-shaped pile, the distance between opposite sides is... The length of the straight boundary is When the irregularly shaped pile is an elliptical pile, the major axis is... The minor axis is .

[0032] Initial stress parameters of soil include initial horizontal stress. Initial vertical stress Initial pore water pressure Initial effective horizontal stress Initial effective vertical stress Initial shear stress and soil undrained shear strength The initial effective horizontal stress and the initial effective vertical stress respectively satisfy the following conditions: , Initial shear stress .

[0033] Soil material parameters include initial void ratio Compression parameters springback parameters Critical stress ratio Poisson's ratio and shear modulus Soil material parameters are used to describe the elastoplastic mechanical behavior of soil.

[0034] Boundary stresses include the far-field horizontal stress, far-field vertical stress, and far-field pore water pressure of the soil. The far-field horizontal stress is set as the initial horizontal stress. The far-field vertical stress of the soil is set as the initial vertical stress. The far-field pore water pressure of the soil is set as the initial pore water pressure. Boundary stress is used to set the stress boundary conditions of the hole expansion model at infinity.

[0035] Assumptions are satisfied: During the expansion of irregular pores, the soil exists in an elastic zone, a plastic zone, and a critical state zone. The elastic-plastic deformation of the soil is described by the K0 anisotropic modified Cambridge constitutive model (the definition of K0 is given below). The vertical strain of the soil is 0. Both the soil skeleton and pore water are incompressible; The boundary of the hole is a smooth boundary.

[0036] Step S2: Construct a hole expansion model based on the relevant parameters, the boundary stress, and the assumptions.

[0037] In this embodiment, as Figure 2 As shown, the pile driving process of an irregularly shaped pile can be viewed as the process of enlarging an irregularly shaped hole in an infinite domain of soil. The initial size of this irregularly shaped hole is 0, and it expands radially until it reaches a size identical to the cross-sectional size of the irregularly shaped pile. During the enlargement process, the shape of the irregularly shaped hole remains consistent with the cross-sectional shape of the irregularly shaped pile. The radius of the non-circular hole (irregularly shaped hole) is... ,in It is the polar angle of the physical plane. The radial velocity of the physical plane. Defined as the radius of a non-circular hole For hole enlargement time The derivative of For the normal velocity of the boundary of the non-circular hole, The tangential velocity at the boundary of the non-circular hole. The initial horizontal stress (i.e., the far-field horizontal stress of the soil) is ), The initial vertical stress is given by the initial pore water pressure of the soil. (i.e., the far-field pore water pressure of the soil is) Therefore, the initial effective horizontal stress Initial effective vertical stress Initial shear stress The shear strength of the undrained soil is 0. The soil displacement is The initial coordinates of the soil particles are , Let be the initial coordinates of the soil particles in the physical plane, where The initial horizontal coordinates of the physical plane. The initial vertical coordinates are given for the physical plane. This invention provides three typical irregularly shaped piles with different cross-sections, such as: Figure 2 (E) in the diagram is a schematic diagram of the cross-section of a rectangular irregular pile; Figure 2 (F) in the diagram is a schematic diagram of the cross-section of the X-shaped pile; Figure 2 (J) in the diagram represents a schematic cross-section of an elliptical irregular pile. For the three typical irregularly shaped piles with different cross-sections presented in this invention, geometric parameters a and b are used to describe the geometric characteristics of the irregular piles. Thus, the above physical description, kinematic definition, boundary stress conditions, and material parameter settings together constitute the borehole enlargement model for driving irregularly shaped piles into soil considering stress anisotropy.

[0038] Step S3: Based on fluid dynamics theory, the radial velocity and circumferential velocity of the image plane are calculated using the pre-acquired conformal transformation function, the aperture expansion model, and the assumed conditions.

[0039] The calculation of the radial velocity and circumferential velocity of the image plane based on fluid dynamics theory using a pre-obtained conformal transformation function, the aperture expansion model, and the assumed conditions specifically includes the following sub-steps: Step S301: Based on fluid mechanics theory, determine the Laplace equation for the velocity of the soil around the hole using the aforementioned assumptions.

[0040] In this embodiment, since the vertical strain of the soil is 0, the volumetric strain rate of the soil should satisfy the following formula: , in, The soil strain rate. and These represent the strain rates in the horizontal x-coordinate and vertical y-coordinate of the physical plane, respectively.

[0041] Under the condition that both the soil skeleton and pore water are incompressible, the soil strain rate is 0, therefore: , Based on the relationship between strain rate and velocity, the above equation can be transformed into: , in, The velocity in the x-direction of the physical plane's horizontal coordinate. The velocity in the vertical y-direction of the physical plane.

[0042] The theoretical framework based on the strain path method states that the motion of the soil is irrotational. and The following relationship should be satisfied: , Based on the above equations and fluid mechanics theory, there exists a potential function φ that satisfies the following equation: , , in, It is a potential function.

[0043] Substituting the above equation into the relationship between strain rate and velocity, we can obtain the Laplace equation for the velocity of the soil around the hole during the driving of irregularly shaped piles: , Step S302: Calculate the radial velocity and circumferential velocity of the image plane based on the Laplace equation of the soil velocity around the borehole, the pre-obtained conformal transformation function, and the borehole enlargement model.

[0044] In this embodiment, for the Laplace equation, since the hole is non-circular, the radial radius in polar coordinates... and physical plane polar angle Since direct separation is not possible, the traditional method of separation of variables is no longer applicable, and a conformal transformation function is needed to solve the problem. The conformal transformation function maps the non-circular hole onto the image plane, transforming the non-circular hole in the original plane (physical plane) into a unit circular hole in the image plane, thereby solving the Laplace equation. Figure 3 (a1) in the diagram is a physical plane diagram. Figure 3 (a2) in the diagram is a schematic diagram of the image plane, showing the initial coordinates of the soil particles on the physical plane. Initial coordinates of soil particles mapped to the image plane , The normal velocity of the non-circular hole in the physical plane (original plane) The tangential velocity at the boundary of the non-circular hole in the physical plane (original plane) is... The radial velocity of the physical plane (original plane) The polar angle of the physical plane. The angle between the normal to the physical plane hole boundary and the horizontal. For the radial velocity of the image plane, For the circumferential velocity of the image plane, The polar angle of the image plane. Represents the unit circle in the image plane. The conformal transformation function used in the solution is as follows: , , , in, For conformal transformation function, These are the coefficients in the conformal transformation function that control the size of the hole. These are the coefficients controlling the shape of the hole in the conformal transformation function. and Determined by the geometric parameters of the non-circular hole. For complex variables in the image plane, For the physical plane horizontal coordinates, For the image physical plane vertical coordinates, Let the radial radius of the image plane be . The polar angle of the image plane, For physical plane complex variables, The horizontal coordinates of the physical plane The vertical coordinates of the physical plane. The radial radius of the physical plane. The polar angle of the physical plane. , It is a natural constant. The value represents the sequence number and ranges from [1, ∞).

[0045] This invention provides conformal transformation function coefficients for rectangular, X-shaped, and elliptical holes, as shown in the table below.

[0046] The coefficient table of conformal transformation function for rectangular holes is as follows:

[0047] In the table of coefficients for the conformal transformation function of the rectangular hole, a and b represent the long side and short side of the rectangle, respectively (see dimensioning). Figure 2 Furthermore, for rectangular holes, an infinite series is used. The first three terms of the polynomial ( , , ,in, , and The coefficients of the first, second, and third polynomials, respectively, are sufficient to meet the required accuracy for the project.

[0048] The coefficient table of the conformal transformation function for the X-shaped hole is as follows:

[0049] In the conformal transformation function coefficient table for the X-shaped hole, 'a' represents the distance between the straight boundaries of the two opposite corners of the X-shape, and 'b' represents the length of the straight boundary of the X-shape (see dimension annotations). Figure 2 Furthermore, for X-shaped holes, an infinite series is used. The first six terms of the polynomial ( , , , , , in, , , , , and The coefficients of the first, second, third, fourth, fifth, and sixth polynomials, respectively, are sufficient to meet the required accuracy for the project.

[0050] The coefficient table of conformal transformation function for elliptical holes is as follows:

[0051] In the conformal transformation function coefficient table for the elliptical hole, a and b represent the major and minor axes of the ellipse, respectively (see dimensioning). Figure 2 Furthermore, for elliptical holes, an infinite series is used. The first term polynomial ( ,in, (where is the coefficient of the first polynomial), which is sufficient to meet the accuracy required for the project.

[0052] The Laplace equation and potential function are transformed into the following form using conformal transformation functions: , , , in, Let ξ be the velocity in the direction of ξ. Let η be the velocity in the direction of η. It is a potential function.

[0053] Because the soil is in a state of motion, that is The above formula can be simplified to: , Transforming the above equation into polar coordinates of the image plane, we obtain the following equation: , Solving the above equation using the method of separation of variables yields the potential function: , in, Let be the potential function. , , and All of these are coefficients in the potential function obtained using the method of separation of variables, and the parameters are different for each of the three hole shapes given in this invention. The value represents the sequence number and ranges from [1, ∞).

[0054] Based on the relationship between the potential function and the velocity, the radial velocity and circumferential velocity of the image plane are obtained: , , in, The expression for the radial velocity in the image plane, Let be the potential function. Let the radial radius of the image plane be . It is the expression for the circumferential velocity in the image plane. The polar angle of the image plane.

[0055] Step S4: Calculate the soil displacement and stress distribution on the image plane based on the hole enlargement model, the radial velocity of the image plane, the circumferential velocity of the image plane, the conformal transformation function, the assumptions, and the pre-acquired K0 anisotropic modified Cambridge constitutive model.

[0056] The calculation of soil displacement and stress distribution in the image plane based on the enlarged borehole model, the radial velocity of the image plane, the circumferential velocity of the image plane, the conformal transformation function, the assumptions, and the pre-obtained K0 anisotropic modified Cambridge constitutive model specifically includes the following sub-steps: Step S401: Calculate the complex velocity variables of the physical plane and the complex strain rate variables of the image plane based on the hole enlargement model, the radial velocity of the image plane, the circumferential velocity of the image plane, the conformal transformation function, and the assumptions.

[0057] The calculation of the complex velocity variables of the physical plane and the complex strain rate variables of the image plane based on the aperture expansion model, the radial velocity of the image plane, the circumferential velocity of the image plane, the conformal transformation function, and the assumptions specifically includes the following sub-steps: Step S40101: Calculate the radial velocity of the physical plane using the conformal transformation function and the hole enlargement model.

[0058] In this embodiment, a conformal transformation function is used to map the velocity boundary of the physical plane onto the image plane (the velocity boundaries of the physical plane and the image plane are as follows). Figure 3 As shown), the following equation holds: , , in, It is a natural constant. The imaginary unit, The angle between the radial and horizontal sides of the physical plane hole boundary (physical plane polar angle). The angle between the normal to the physical plane hole boundary and the horizontal. The complex conjugate of the first derivative of the conformal transformation function. The first derivative of the conformal transformation function, Let be the magnitude of the first derivative of the conformal transformation function. Conformal transformation function The length of the module.

[0059] This invention defines Represents a point on the unit circle of the image plane ( ).

[0060] According to the conformal transformation function, the radius of the non-circular hole in the physical plane is obtained as follows: , in, Let be the radius of a non-circular hole in a physical plane. for The length of the mold, for The conformal transformation function value.

[0061] According to the above equation, the radial velocity of the physical plane is obtained. : , in, It is the radial velocity of the physical plane. This is the coefficient that controls the size of the hole in the conformal transformation function.

[0062] Step S40102: Calculate the radial velocity on the unit circle of the image plane using the conformal transformation function, the radial velocity of the physical plane, and the assumed conditions.

[0063] In this embodiment, when the hole boundary is a smooth boundary, the normal velocity of the physical plane hole boundary is equal to the projection of the physical plane radial velocity onto the boundary normal, that is: , , in, Let be the normal velocity at the boundary of the physical plane hole. It is the cosine of the angle between the normal and radial directions of the physical plane hole boundary. The radial velocity of the physical plane. It is the physical plane normal velocity (i.e., the normal velocity of the boundary of a non-circular hole). The angle between the normal of the physical plane hole boundary and the horizontal.

[0064] Will , and Substituting the expression into the above equation, we obtain the physical plane normal velocity. The expression: , According to the conformal transformation function, the velocities of the physical plane and the image plane should satisfy the following relationship: , in, This is the first derivative of the conformal transformation function.

[0065] According to the physical plane coordinate transformation relationship, the physical plane velocity must satisfy the following relationship: , Combining the above equations, we obtain the following expression: , Will Substituting the expression into the above formula, we obtain the radial velocity on the unit circle of the image plane: , in, for The conformal transformation function value, for The first derivative, for The conjugate of complex numbers, Let the radial velocity be on the unit circle of the image plane. Let be the real part of the complex number. These are the coefficients in the conformal transformation function that control the size of the hole. For conformal transformation function, It is a natural constant. The imaginary unit, For complex variables in the image plane, Let be a point on the unit circle of the image plane. The velocity in the x-direction of the physical plane's horizontal coordinate. The velocity in the vertical y-direction of the physical plane. The tangential velocity at the boundary of the non-circular hole. The radial velocity is the image plane.

[0066] Step S40103: Calculate the complex velocity variables of the physical plane based on the radial velocity on the unit circle of the image plane, the radial velocity of the image plane, the circumferential velocity of the image plane, and the conformal transformation function.

[0067] Will Substitute the radial velocity of the image plane The radial velocity on the unit circle of the image plane is obtained as follows: , in, Let be the radial velocity on the unit circle of the image plane.

[0068] because and If both are expressed as radial velocities on the unit circle of the image plane, then the following equation holds: , Using the above equation, the coefficients can be obtained through integration. and The expression: , , in, and These are all symbols defined to simplify mathematical expressions.

[0069] because It is about an even function, then .

[0070] Will and Substitute the radial velocity of the image plane and circumferential velocity have to: , According to the conformal transformation function, the image plane and the physical plane have the following velocity relationship: , Combining the above equations, we obtain the following expression: , To facilitate subsequent calculations, we define the velocity complex variable of the physical plane. The definition is as follows: , in, Let the velocity of the physical plane be a complex variable. For complex variables in the image plane The conjugate of complex numbers, The complex conjugate of the first derivative of the conformal transformation function. The coefficients in the potential function obtained using the method of separation of variables, To simplify symbols, Let the radial radius of the image plane be . The polar angle of the image plane, For the radial velocity of the image plane, For the circumferential velocity of the image plane, Let be a point on the unit circle of the image plane. for The conformal transformation function value, for The conjugate of complex numbers, For conformal transformation functions in The first derivative value at that point, These are the coefficients in the conformal transformation function that control the size of the hole. Let be the real part of the complex number. The velocity in the x-direction of the physical plane's horizontal coordinate. The velocity in the vertical y-direction of the physical plane.

[0071] Step S40104: Based on the complex variable function theory, the complex variable of the velocity of the physical plane is used to calculate the complex variable of the strain rate of the image plane.

[0072] According to the definition of strain rate, we have: , , , in, It is the shear strain rate of the physical plane.

[0073] To facilitate subsequent calculations, a new complex variable is defined. and The definition is as follows: , , in, It is the complex strain rate variable of the physical plane. It is a complex strain rate variable like a plane. It is the strain rate in the radial direction of the image plane. It is the strain rate in the circumferential direction of the image plane. It is the shear strain rate of the plane.

[0074] Based on the theory of differentiation of multivariable functions and the theory of complex functions, the following differential expression is obtained: , , Will , , Substituting the above differential expression into From the definition, we get: , According to the conformal transformation function, the strain rate complex variable and The following relationship must be satisfied: , Will Substituting the expression into the above relationship, we can obtain the complex variable of strain rate of the image plane: , in, It is a complex strain rate variable like a plane. The expression is as follows: , in, Let the strain rate be a complex variable on the image plane. Let the velocity of the physical plane be a complex variable. For complex variables in the image plane The conjugate of complex numbers, For complex variables in the image plane, The complex conjugate of the first derivative of the conformal transformation function. The complex conjugate of the second derivative of the conformal transformation function. The coefficients in the potential function obtained using the method of separation of variables, For simplification symbols.

[0075] Step S402: Calculate the soil displacement and stress distribution on the image plane using the complex velocity variable of the physical plane, the complex strain rate variable of the image plane, and the pre-acquired K0 anisotropic modified Cambridge constitutive model.

[0076] The calculation of soil displacement and stress distribution on the image plane using the complex velocity variable of the physical plane, the complex strain rate variable of the image plane, and the pre-acquired K0 anisotropic modified Cambridge constitutive model specifically includes the following sub-steps: Step S40201: Using the complex strain rate variable of the image plane and the complex velocity variable of the physical plane, determine the kinematic equation of the soil in the image plane.

[0077] Soil displacement and strain can be obtained by integrating over velocity and strain rate, respectively, as shown in the following integral expressions: , , , , in, It is like a plane strain complex variable. yes The conjugate of complex numbers, It is the increment of the complex variable in the physical plane. yes The conjugate of complex numbers, It is like the increment of a complex plane strain variable. yes The conjugate of complex numbers, The velocity complex variable of the physical plane The conjugate of complex numbers.

[0078] To account for the large deformation effect in the plastic region, the above integral expression is transformed into the following differential expression: , , , , in, Let z be the differential of the physical plane complex variable z with respect to R. yes The conjugate of complex numbers, For plane strain complex variables Differential with respect to R, yes The conjugate of complex numbers.

[0079] According to the theory of differentiation of multivariable functions, the following equation holds: , , Substituting the differential expression into the equation above, we obtain the following equation: , , , , in, For complex variables in the image plane Differential with respect to R, yes The conjugate of complex numbers, , , and These are all simplified symbols for differential equations.

[0080] The above equations are transformed into the following form to obtain the kinematic equations of the image plane soil: , Step S40202: Calculate the elastoplastic constitutive matrix of the image plane soil using the yield function of the pre-acquired K0 anisotropic modified Cambridge constitutive model.

[0081] In this embodiment, the stress-strain matrix of the elastic region in the image plane is obtained according to Hooke's law, and its expression is as follows: , in, Shear modulus Poisson's ratio, Radial strain in the elastic region of the image plane right The differential, For the circumferential strain of the elastic region of the image plane right The differential, For the strain in the vertical direction of the elastic region of the image plane right The differential, For the shear strain of the elastic region of the image plane right The differential, Effective radial stress in the image plane right The differential, Effective circumferential stress in the image plane right The differential, Effective vertical stress in the image plane right The differential, For plane shear stress right The differential.

[0082] The elastoplastic deformation is described by the K0 anisotropic modified Cambridge constitutive model, and the yield function of the K0 anisotropic modified Cambridge constitutive model is as follows: , Where λ is the compression parameter and κ is the springback parameter. The initial void ratio, For the average effective stress, The initial average effective stress, For plastic strain, The relative stress ratio The relative stress ratio under critical conditions is expressed as follows: , , , , , , , , , , , , Where M is the critical stress ratio. , and These are all symbols defined to simplify mathematical expressions. , , and They are respectively The radial, circumferential, vertical, and shear components, , , and They are respectively The radial, circumferential, vertical, and shear components, The symbol for Kroneck. The initial effective stress ratio, For the effective radial stress in the image plane, For the effective circumferential stress in the image plane, For the effective vertical stress of the image plane, The initial effective radial stress of the image plane, The initial effective circumferential stress of the image plane, The initial effective vertical stress of the image plane, For plane shear stress, For the initial shear stress in the image plane, For effective stress, This represents the initial effective stress.

[0083] Based on the yield function, hardening law, and associated flow rule of the K0 anisotropic modified Cambridge constitutive model, the following equation is obtained: , , , , in, For the strain increment in the plastic region, For the radial strain increment in the plastic region of the image plane, For the circumferential strain increment in the plastic region of the image plane, For the vertical strain increment in the plastic region of the image plane, For the shear strain increment in the plastic region of the image plane, For the effective stress increment, For the effective radial stress increment in the image plane, For the effective circumferential stress increment in the image plane, For the effective vertical stress increment in the image plane, For the shear stress increment in the image plane, Symbols defined to simplify mathematical expressions.

[0084] The above expression is simplified to the following form, yielding the stress-strain increment matrix of the plastic region in the image plane: , in, Radial strain in the plastic region of the image plane right The differential, Circumferential strain in the plastic region of the image plane right The differential, For the vertical strain of the plastic region of the image plane right The differential, For the shear strain of the plastic region of the image plane right The differential, , , , and These are all symbols defined to simplify mathematical expressions, as shown below: , , , , , Using the superposition principle, the stress-strain matrix of the elastic region and the stress-strain increment matrix of the plastic region in the image plane are superimposed to obtain the elastoplastic constitutive matrix of the soil in the image plane: , , , , , , , , , , , in, , , , , , , , , , , , , , , and These are all symbols defined to simplify mathematical expressions. For image plane radial strain Differential with respect to R, For the circumferential strain of the image plane Differential with respect to R, For image plane perpendicular strain Differential with respect to R, For image plane shear strain Differential with respect to R.

[0085] To maintain consistency with the kinematic equations, the strain in the elastoplastic matrix needs to be transformed into a complex variable form. The transformed elastoplastic constitutive matrix of the image plane soil is as follows: , , Here, D is a matrix defined to simplify the equation.

[0086] Step S40203: Obtain the constitutive equation of the image plane based on the elastic-plastic constitutive matrix of the soil in the image plane.

[0087] In this embodiment, utilizing the properties of the inverse matrix, the elastoplastic constitutive matrix of the image plane soil can be transformed into the following form: , make The following equation is obtained: , Here, V is a column vector with four rows and one column.

[0088] To maintain consistency, the above equation is rewritten in the following form to obtain the constitutive equation of the image plane. , in, , , and These are all simplified symbols for differential equations, specifically: The expression representing the first line of V; The expression representing the second line of V; The expression in the third line representing V; The expression in the fourth line representing V.

[0089] Step S40204: Calculate the soil displacement and stress distribution in the image plane using the kinematic equations and constitutive equations of the image plane.

[0090] In this embodiment, the kinematic equations of the image plane soil and the constitutive equations of the image plane soil are combined to obtain the soil governing differential equations: , Based on the initial coordinates of the soil particles in the image plane , to obtain a plane complex variable Initial values: , in, Let ζ be the initial value of the complex variable ζ in the image plane. Since the soil deforms to zero in its initial state, the following equation can be obtained: , , , , , in, Let these be the initial values ​​for the complex variables of plane strain. Let the initial value be the effective radial stress in the image plane. Let be the initial value of the effective circumferential stress in the image plane. Let this be the initial value of the effective perpendicular stress in the image plane. This is the initial value of the shear stress in the image plane.

[0091] Based on the above equation, the initial values ​​of the soil governing differential equation are obtained: , in, for The conjugate of complex numbers, for The conjugate of complex numbers.

[0092] Substituting the initial values ​​into the soil governing differential equation and solving iteratively, we obtain the soil displacement and stress distribution in the image plane.

[0093] Step S5: Calculate the stress field and displacement field of the soil body for the irregular pile driving using the assumed conditions, the soil displacement in the image plane, the stress distribution in the image plane, and the conformal transformation function.

[0094] The calculation of the stress field and displacement field of the soil for driving irregularly shaped piles using the aforementioned assumptions, the soil displacement in the image plane, the stress distribution in the image plane, and the conformal transformation function specifically includes the following sub-steps: Step S501: Calculate the physical plane soil displacement and physical plane stress distribution using the conformal transformation function, the image plane soil displacement, and the image plane stress distribution.

[0095] In this embodiment, according to the conformal transformation function, the stress-strain transformation relationship between the image plane and the physical plane is as follows: , , , , , , in, For the effective radial stress of the physical plane, For the effective circumferential stress of the physical plane, For physical plane shear stress, For physical plane radial strain, For physical plane circumferential strain For physical plane shear strain, For stress or strain on a physical plane, For stress or strain on a plane, It is the physical planar displacement of soil. This represents the displacement of the soil mass in the image plane.

[0096] Substituting the image plane soil displacement and stress distribution into the above transformation relationship, we obtain the physical plane soil displacement and physical plane stress distribution.

[0097] By utilizing the physical plane soil displacement, the normalized soil displacement of irregularly shaped holes is obtained.

[0098] By utilizing the physical plane stress distribution, the normalized effective stress of the soil in irregularly shaped holes is obtained.

[0099] By utilizing the physical plane stress distribution, the normalized shear stress of the soil in irregularly shaped holes is obtained.

[0100] This invention provides normalized displacement diagrams, normalized effective stress distribution diagrams, and normalized shear stress diagrams for rectangular holes, X-shaped holes, and elliptical holes (since the rectangular holes, X-shaped holes, and elliptical holes are all axisymmetric figures, images of 1 / 4 of the holes are given here).

[0101] Figure 4 (a3) in the figure represents the normalized displacement of the soil mass in the rectangular borehole (in the figure). For physical planar displacement, Let be the radius of the irregular hole, and at this time, Take the special value 0, Radial radius For normalized displacement, For normalized radial distance, Figure 4 The same applies to the other diagrams). Figure 4 (a4) in the figure is the normalized displacement diagram of the soil in the X-shaped hole; Figure 4 (a5) in the figure represents the normalized displacement of the soil mass in the elliptical borehole. From... Figure 4 It can be seen that the normalized displacement of the X-shaped and elliptical holes has a similar variation pattern, with a faster initial descent rate. However, the normalized displacement of the rectangular hole decreases more slowly in the initial stage, and the normalized displacement of all three irregular holes tends to 0 as the normalized radial distance increases.

[0102] Figure 5 (a6) in the diagram is the normalized effective stress distribution of the soil around the rectangular hole; Figure 5 (a7) in the figure is the normalized effective stress distribution diagram of the soil for the X-shaped hole (in the figure) The polar angle of the physical plane. , and These are the effective radial, circumferential, and vertical stresses in the physical plane, respectively. For effective stress, The initial average effective stress, To normalize the effective stress, Figure 5 The same applies to the other diagrams). Figure 5 (a8) in the diagram represents the normalized effective stress distribution of the soil around the elliptical hole. Figure 5 As can be seen, the effective stress curves of rectangular and X-shaped holes are more complex, while the trend of elliptical curves is simpler. This is mainly due to the complex geometry of rectangular and X-shaped holes.

[0103] Figure 6 (a9) in the figure is the normalized shear stress distribution cloud map of the soil in the rectangular hole (in the figure) For physical plane shear stress, For undrained shear strength, To normalize the shear stress, Figure 6 The same applies to the other diagrams). Figure 6 (a) 10 ( ) is a normalized shear stress distribution cloud map of soil in an X-shaped hole; Figure 6 (a) 11 (This is a normalized shear stress distribution cloud map of the soil mass within an elliptical borehole.) Figure 6 It can be seen that both the X-shaped and rectangular shapes exhibit significant stress concentration at their corners, and both tend to... ( With the polar angle as the axis of symmetry, the cloud map of the ellipse is wavy and has no negative values.

[0104] Step S502: Based on the stress balance equation, the physical plane pore water pressure is calculated using the assumptions, the conformal transformation function, and the image plane stress distribution.

[0105] In this embodiment, based on the fact that the vertical strain of the soil is 0 and combined with the conformal transformation function, the soil stress balance equation is as follows: , in, For radial stress in the image plane, For the circumferential stress in the image plane, H is a symbol defined to simplify the mathematical expression. Let H be the derivative of ρ. Let H be the derivative of ω, expressed as follows: , , , According to the effective stress principle: , , Where uphase is the pore water pressure in the image plane soil.

[0106] Substituting the above equation into the stress balance equation, we get: , Integrating the above equation along the ρ direction yields the expression for the pore water pressure in the image plane: , in, This represents the initial pore water pressure.

[0107] Substituting the stress distribution in the image plane into the pore water pressure expression, the pore water pressure of the soil in the image plane is obtained.

[0108] The relationship between pore water pressure in the physical plane and pore water pressure in the image plane is as follows: , Where u is the physical plane pore water pressure.

[0109] The physical plane pore water pressure is obtained based on the above relationship.

[0110] This invention provides normalized pore water pressure distribution cloud maps of soil for rectangular, X-shaped, and elliptical holes. (Since the given rectangular, X-shaped, and elliptical holes are all axisymmetric figures, images of 1 / 4 of the hole are shown here).

[0111] Figure 7 (a) 12 (This is a normalized pore water pressure distribution cloud map of a rectangular pore soil mass.) For physical plane pore water pressure, For undrained shear strength, To normalize the pore water pressure, Figure 7 The same applies to the other diagrams); Figure 7 (a) 13 This is a cloud map showing the normalized pore water pressure distribution in soil with X-shaped pores. Figure 7 (a) 14 (This is a contour map of normalized pore water pressure distribution in elliptical pore soil.) Figure 7 It can be seen that the extreme values ​​of the pore pressure of the rectangular and X-shaped holes appear at the corners, while the extreme values ​​of the elliptical holes appear around the endpoints of the major axis. However, similarly, negative pore pressure values ​​appear in the middle area of ​​all three figures.

[0112] The physical plane pore water pressure includes: , , , , in, For physical plane soil pore water pressure, For example, the pore water pressure of a planar soil mass, The initial pore water pressure, Let the radial radius of the image plane be . For the effective radial stress in the image plane, For plane shear stress, The polar angle of the image plane, For the effective circumferential stress in the image plane, These are the coefficients in the conformal transformation function that control the size of the hole. The first derivative of the conformal transformation function, The complex conjugate of the first derivative of the conformal transformation function. Let be the real part of the complex number. The second derivative of the conformal transformation function, The imaginary unit, Let be a point on the unit circle of the image plane. , and All are simplified symbols.

[0113] Step S503: Calculate the stress field and displacement field of the soil body for the irregular pile driving based on the pore water pressure in the physical plane, the soil displacement in the physical plane, and the stress distribution in the physical plane.

[0114] In this embodiment, the physical plane pore water pressure, the physical plane soil displacement, and the physical plane stress distribution are combined to obtain the stress field and displacement field of the soil for the irregular pile driving.

[0115] like Figure 8 As shown, in another aspect, the present invention also provides a calculation system for the soil squeezing effect of irregularly shaped pile driving considering stress anisotropy, comprising: a processor, an input device, an output device, and a memory, wherein the processor, the input device, the output device, and the memory are interconnected, wherein the memory is used to store a computer program, the computer program including program instructions, and the processor is configured to call the program instructions to execute the relevant steps of the relevant embodiments of the calculation method for the soil squeezing effect of irregularly shaped pile driving considering stress anisotropy of the present invention.

[0116] The soil squeezing effect calculation system for irregularly shaped pile driving considering stress anisotropy provided by this invention allows for the integration of functional components into a single processing unit, or for each component to exist independently, or for two or more components to be integrated into a single unit. These integrated components can be implemented in hardware or software.

[0117] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some or all of the technical features therein. Such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of the present invention, and they should all be covered within the scope of the claims and specification of the present invention.

Claims

1. A method for calculating the soil squeezing effect of irregularly shaped piles considering stress anisotropy, characterized in that, The method includes: Obtain the relevant parameters of the irregular pile and set the boundary stress and assumptions; A hole enlargement model is constructed based on the relevant parameters, the boundary stress, and the assumptions. Based on fluid dynamics theory, the radial velocity and circumferential velocity of the image plane are calculated using the pre-obtained conformal transformation function, the aperture expansion model, and the assumed conditions. Based on the hole enlargement model, the radial velocity of the image plane, the circumferential velocity of the image plane, the conformal transformation function, the assumptions, and the pre-obtained K0 anisotropic modified Cambridge constitutive model, the soil displacement and stress distribution in the image plane are calculated. The stress field and displacement field of the soil for driving irregularly shaped piles are calculated using the aforementioned assumptions, the soil displacement in the image plane, the stress distribution in the image plane, and the conformal transformation function.

2. The method for calculating the soil squeezing effect of irregularly shaped piles considering stress anisotropy as described in claim 1, characterized in that, The calculation of the radial velocity and circumferential velocity of the image plane based on fluid dynamics theory using a pre-obtained conformal transformation function, the aperture enlargement model, and the assumed conditions includes: Based on fluid mechanics theory, the Laplace equation for the velocity of the soil around the hole is determined using the aforementioned assumptions. The radial velocity and circumferential velocity of the image plane are calculated based on the Laplace equation of the soil velocity around the borehole, the pre-obtained conformal transformation function, and the borehole enlargement model.

3. The method for calculating the soil squeezing effect of irregularly shaped piles considering stress anisotropy as described in claim 1, characterized in that, The calculation of soil displacement and stress distribution in the image plane based on the enlarged borehole model, the radial velocity of the image plane, the circumferential velocity of the image plane, the conformal transformation function, the assumptions, and the pre-obtained K0 anisotropic modified Cambridge constitutive model includes: Based on the hole enlargement model, the radial velocity of the image plane, the circumferential velocity of the image plane, the conformal transformation function, and the assumptions, calculate the complex velocity variables of the physical plane and the complex strain rate variables of the image plane; The soil displacement and stress distribution on the image plane are calculated using the complex velocity variable of the physical plane, the complex strain rate variable of the image plane, and the pre-acquired K0 anisotropic modified Cambridge constitutive model.

4. The method for calculating the soil squeezing effect of irregularly shaped piles considering stress anisotropy as described in claim 3, characterized in that, The calculation of the complex velocity variables of the physical plane and the complex strain rate variables of the image plane based on the aperture expansion model, the radial velocity of the image plane, the circumferential velocity of the image plane, the conformal transformation function, and the assumptions includes: Using the conformal transformation function and the hole enlargement model, the radial velocity of the physical plane is calculated; Using the conformal transformation function, the radial velocity of the physical plane, and the assumptions, calculate the radial velocity on the unit circle of the image plane; The velocity complex variables of the physical plane are calculated based on the radial velocity on the unit circle of the image plane, the radial velocity of the image plane, the circumferential velocity of the image plane, and the conformal transformation function. Based on complex variable function theory, the complex variable of the velocity of the physical plane is used to calculate the complex variable of the strain rate of the image plane.

5. The method for calculating the soil squeezing effect of irregularly shaped piles considering stress anisotropy as described in claim 4, characterized in that, The complex velocity variable of the physical plane satisfies the following formula: ， in, Let the velocity of the physical plane be a complex variable. For complex variables in the image plane The conjugate of complex numbers, The complex conjugate of the first derivative of the conformal transformation function. The coefficients in the potential function obtained using the method of separation of variables, To simplify symbols, For serial numbers.

6. The method for calculating the soil squeezing effect of irregularly shaped piles considering stress anisotropy as described in claim 4, characterized in that, The complex strain rate variable of the image plane satisfies the following formula: ， ， in, Let the strain rate be a complex variable on the image plane. Let the velocity of the physical plane be a complex variable. For complex variables in the image plane The conjugate of complex numbers, The first derivative of the conformal transformation function, The complex conjugate of the first derivative of the conformal transformation function. The complex conjugate of the second derivative of the conformal transformation function. The coefficients in the potential function obtained using the method of separation of variables, To simplify symbols, For serial numbers.

7. The method for calculating the soil squeezing effect of irregularly shaped piles considering stress anisotropy as described in claim 3, characterized in that, The calculation of soil displacement and stress distribution on the image plane using the complex velocity variable of the physical plane, the complex strain rate variable of the image plane, and the pre-acquired K0 anisotropic modified Cambridge constitutive model includes: The kinematic equations of the soil in the image plane are determined using the complex strain rate variable of the image plane and the complex velocity variable of the physical plane. The elastoplastic constitutive matrix of the image plane soil is calculated using the yield function of the pre-acquired K0 anisotropic modified Cambridge constitutive model. Based on the elastoplastic constitutive matrix of the soil in the image plane, the constitutive equation of the image plane is obtained; Using the kinematic equations and constitutive equations of the image plane soil, the displacement and stress distribution of the image plane soil are calculated.

8. The method for calculating the soil squeezing effect of irregularly shaped piles considering stress anisotropy as described in claim 1, characterized in that, The calculation of the stress field and displacement field of the soil for irregularly shaped pile driving using the assumed conditions, the soil displacement in the image plane, the stress distribution in the image plane, and the conformal transformation function includes: Using the conformal transformation function, the image plane soil displacement, and the image plane stress distribution, the physical plane soil displacement and physical plane stress distribution are calculated; Based on the stress balance equation, the physical plane pore water pressure is calculated using the aforementioned assumptions, the conformal transformation function, and the image plane stress distribution. The stress field and displacement field of the soil body for the irregular pile driving are calculated based on the pore water pressure in the physical plane, the soil displacement in the physical plane, and the stress distribution in the physical plane.

9. The method for calculating the soil squeezing effect of irregularly shaped piles considering stress anisotropy as described in claim 8, characterized in that, The physical plane pore water pressure includes: ， ， ， ， in, For physical plane soil pore water pressure, For example, the pore water pressure of a planar soil mass, The initial pore water pressure, Let the radial radius of the image plane be . For the effective radial stress in the image plane, For plane shear stress, The polar angle of the image plane, For the effective circumferential stress in the image plane, These are the coefficients in the conformal transformation function that control the size of the hole. The first derivative of the conformal transformation function, The complex conjugate of the first derivative of the conformal transformation function. Let be the real part of the complex number. The second derivative of the conformal transformation function, The imaginary unit, Let be a point on the unit circle of the image plane. , and All are simplified symbols.

10. A calculation system for the soil squeezing effect of irregularly shaped piles considering stress anisotropy, characterized in that, include: The system includes a processor, an input device, an output device, and a memory, all interconnected, wherein the memory stores a computer program comprising program instructions, and the processor is configured to invoke the program instructions to execute the calculation method for the soil squeezing effect of irregularly shaped pile driving considering stress anisotropy as described in any one of claims 1 to 9.

Images