Numerical simulation method and device for simulating square pile hole-forming effect

Abstract

The application provides a numerical simulation method and device for square pile hole forming effect simulation, and relates to the technical field of geotechnical engineering and pile foundation construction, which comprises the following steps: after grid division is performed on an integrated model, the integrated model is imported into numerical analysis software, material parameters are given to a slope in the integrated model, and boundary conditions are set; a particle flow model based on discrete elements is embedded in a preset drilling area of the integrated model to obtain a drill bit-stratum coupling numerical model; numerical simulation of rig drilling and square pile hole forming is performed on the drill bit-stratum coupling numerical model, and force data of the drill bit and deformation response data of a hole wall are obtained in real time; multi-parameter coupling analysis is performed on the force data of the drill bit and the deformation response data of the hole wall, and a drill bit cutting stress-hole wall displacement coupling relationship is established. Through establishment of a quantitative coupling relationship between drill bit cutting stress and hole wall displacement, the application can accurately predict the hole forming stability and adaptability of the rig, and improves the construction reliability of square pile hole forming.

Description

Technical Field

[0001] This invention relates to the fields of geotechnical engineering and pile foundation construction technology, and more specifically, to a numerical simulation method and apparatus for simulating the hole-forming effect of square piles. Background Technology

[0002] In pile foundation engineering, manually excavated piles and circular drilled piles are two mainstream construction methods. However, the former suffers from low construction efficiency and high safety risks, while the latter is prone to borehole instability and poor borehole shape due to uneven distribution of cutting force from the rotating drill bit. For square piles, which are better matched with rectangular pile caps, mechanized drilling technology is still immature and lacks effective theoretical guidance and design tools. Existing numerical simulation methods are mostly based on the finite element theory of continuous media, which makes it difficult to realistically simulate the discontinuous deformation behavior of soil such as particle-level crushing and spalling during the cutting process of a square drill bit. This results in a large deviation between the simulation results and the actual situation, and cannot provide a reliable basis for the optimized design and construction control of square pile drilling rigs. Summary of the Invention

[0003] The purpose of this invention is to provide a numerical simulation method and apparatus for simulating the hole-forming effect of square piles, thereby improving the aforementioned problems. To achieve the above objective, the technical solution adopted by this invention is as follows:

[0004] Firstly, this application provides a numerical simulation method for simulating the hole-forming effect of square piles, including:

[0005] Geological data of the target slope is acquired, a three-dimensional model of the slope including the preset borehole area is established based on the geological data, and the three-dimensional model of the slope is spatially registered with the independently constructed drill bit model to obtain an integrated model;

[0006] After meshing the integrated model, it is imported into numerical analysis software to assign material parameters and set boundary conditions for the slopes in the integrated model.

[0007] A discrete element-based particle flow model is embedded within the pre-defined borehole area of ​​the integrated model, and the particle flow model is brought to an initial stress equilibrium state through servo control to obtain a drill bit-formation coupled numerical model.

[0008] Numerical simulations of drilling rig drilling and square pile hole formation were performed on the drill bit-formation coupled numerical model, and the stress data of the drill bit and the deformation response data of the hole wall were acquired in real time.

[0009] Multi-parameter coupling analysis was performed on the stress data of the drill bit and the deformation response data of the borehole wall to establish the coupling relationship between the drilling rig cutting stress and the borehole wall displacement.

[0010] Secondly, this application also provides a numerical simulation device for simulating the hole-forming effect of square piles, comprising:

[0011] The geological data acquisition module is used to acquire geological data of the target slope, establish a three-dimensional model of the slope including the preset borehole area based on the geological data, and spatially register the three-dimensional model of the slope with the independently constructed drill bit model to obtain an integrated model.

[0012] The mesh processing module is used to mesh the integrated model and import it into the numerical analysis software, assign material parameters to the slopes in the integrated model and set boundary conditions.

[0013] The particle flow modeling module is used to embed a discrete element-based particle flow model into a preset borehole area of ​​the integrated model, and to make the particle flow model reach an initial stress equilibrium state through servo control, thereby obtaining a drill bit-formation coupled numerical model.

[0014] The drilling simulation module is used to perform numerical simulation of drilling rig drilling and square pile hole formation on the drill bit-formation coupled numerical model, and to acquire the stress data of the drill bit and the deformation response data of the hole wall in real time.

[0015] The data analysis module is used to perform multi-parameter coupling analysis on the stress data of the drill bit and the deformation response data of the borehole wall to establish the coupling relationship between the drilling rig cutting stress and the borehole wall displacement.

[0016] Thirdly, this application also provides a numerical simulation device for simulating the hole-forming effect of square piles, comprising:

[0017] Memory, used to store computer programs;

[0018] A processor is used to execute the computer program to implement the steps of a numerical simulation method for simulating the hole-forming effect of the square pile.

[0019] Fourthly, this application also provides a readable storage medium storing a computer program, which, when executed by a processor, implements the steps of the numerical simulation method based on the simulation of the hole-forming effect of square piles described above.

[0020] The beneficial effects of this invention are as follows:

[0021] This invention integrates multi-source geological data, 3D modeling, continuous-discrete coupled numerical calculation, and multi-parameter collaborative analysis to construct a complete simulation and evaluation system for square pile drilling effects. This system achieves high-precision reconstruction of the drilling rig-soil interaction mechanism during square pile drilling in complex strata. By establishing a quantitative coupling relationship between drilling rig cutting stress and borehole wall displacement, the system can accurately predict the drilling stability and adaptability of the drilling rig, thereby improving the construction reliability of square pile drilling.

[0022] Other features and advantages of the invention will be set forth in the following description, and will be apparent in part from the description, or may be learned by practicing embodiments of the invention. The objects and other advantages of the invention may be realized and obtained by means of the structures particularly pointed out in the written description, claims, and drawings. Attached Figure Description

[0023] To more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings used in the embodiments will be briefly introduced below. It should be understood that the following drawings only show some embodiments of the present invention and should not be regarded as a limitation on the scope. For those skilled in the art, other related drawings can be obtained based on these drawings without creative effort.

[0024] Figure 1 This is a schematic diagram of the numerical simulation method for simulating the hole-forming effect of square piles as described in this embodiment of the invention;

[0025] Figure 2 This is the integrated model diagram described in the embodiments of the present invention;

[0026] Figure 3 This is a schematic diagram illustrating the hole-forming effect of the square pile as described in an embodiment of the present invention;

[0027] Figure 4 This is the drilling rig cutting stress cloud diagram described in the embodiments of the present invention;

[0028] Figure 5 This is the stress cloud diagram of the drilling rig described in the embodiments of the present invention;

[0029] Figure 6 This is a schematic diagram of the numerical simulation equipment used to simulate the hole-forming effect of square piles as described in this embodiment of the invention.

[0030] Marked in the image:

[0031] 800. Numerical simulation equipment for simulating the hole-forming effect of square piles; 801. Processor; 802. Memory; 803. Multimedia components; 804. I / O interface; 805. Communication components. Detailed Implementation

[0032] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of the present invention. The components of the embodiments of the present invention described and shown in the accompanying drawings can generally be arranged and designed in various different configurations. Therefore, the following detailed description of the embodiments of the present invention provided in the accompanying drawings is not intended to limit the scope of the claimed invention, but merely to illustrate selected embodiments of the invention. All other embodiments obtained by those skilled in the art based on the embodiments of the present invention without inventive effort are within the scope of protection of the present invention.

[0033] It should be noted that similar reference numerals and letters in the following figures indicate similar items; therefore, once an item is defined in one figure, it does not need to be further defined and explained in subsequent figures. Furthermore, in the description of this invention, terms such as "first," "second," etc., are used only to distinguish descriptions and should not be construed as indicating or implying relative importance.

[0034] Example 1:

[0035] This embodiment provides a numerical simulation method for simulating the hole-forming effect of square piles.

[0036] See Figure 1 , Figure 2 The figure shows that this method includes:

[0037] S1. Obtain geological data of the target slope, establish a three-dimensional model of the slope including the preset drilling area based on the geological data, and spatially register the three-dimensional model of the slope with the independently constructed drill bit model to obtain an integrated model;

[0038] Specifically, step S1 includes:

[0039] S11. Use lidar to scan the target slope and obtain three-dimensional terrain point cloud data;

[0040] S12. The three-dimensional terrain point cloud data is digitally processed to extract the terrain morphology, slope, aspect and elevation information of the target slope;

[0041] Specifically, the 3D terrain point cloud data is imported into Context Capture software for denoising, filtering, and stitching to generate a high-precision digital elevation model. Subsequently, in Global Mapper software, the surface analysis toolset is used to automatically extract information such as slope gradient, aspect, strike, and dip angle from the digital elevation model.

[0042] S13. Obtain the layering information of underground rock and soil masses through drilling sampling and geological exploration in order to determine the thickness of each stratum and the data of rock and soil mechanical parameters;

[0043] Specifically, at least three geological boreholes should be drilled in the target slope area, and rotary core sampling should be performed using an XY-100 drilling rig, with a drilling depth of no less than 5 meters below the potential slip surface. The extracted core samples should be logged to determine the burial depth, thickness, lithology, density, elastic modulus, Poisson's ratio, cohesion, internal friction angle, and tensile strength of each soil and rock layer.

[0044] Specifically, step S1 further includes:

[0045] S14. Based on the surface modeling method, a three-dimensional model of the slope is established according to the topographic morphology, slope, aspect and elevation information of the target slope;

[0046] Specifically, the geological data is imported into Rhinoceros 7.0 software, and the NURBS surface modeling method is used to generate topographic surfaces based on surface elevation points. Then, based on the contour lines of each stratigraphic interface, the corresponding three-dimensional slope model is constructed.

[0047] S15. Based on the thickness data of each stratum, generate the stratigraphic interfaces representing different soil and rock layers in the three-dimensional model of the slope;

[0048] S16. Cut out the preset borehole area in the three-dimensional model of the slope of the integrated stratigraphic interface to obtain the final three-dimensional model of the slope; specifically, the preset borehole area has a square cross section of 400mm×400mm.

[0049] Based on the above embodiments, this method further includes:

[0050] S2. After meshing the integrated model, import it into numerical analysis software, assign material parameters to the slopes in the integrated model, and set boundary conditions.

[0051] Specifically, step S2 includes:

[0052] S21. In the numerical analysis software, the slope and drill bit models of the integrated model are grouped, and each soil and rock layer is divided according to the stratigraphic interface.

[0053] Specifically, the integrated model, which has already been spatially registered and uniformly meshed in Rhinoceros, is exported as a .f3grid format mesh file that can be recognized by FLAC3D through the Griddle plugin; in FLAC3D, different components in the overall model are grouped, preferably with the geological body named slope and the drill bit named cutter.

[0054] S22. Define an elastic constitutive model for each soil and rock layer, and assign the corresponding elastic modulus and Poisson's ratio parameters respectively;

[0055] S23. Apply fixed constraints to the bottom boundary of the integrated model and roll constraints to the four lateral boundaries of the integrated model.

[0056] Specifically, fixed constraints are applied to all mesh nodes at the bottom of the integrated model, and rolling constraints are applied to the lateral boundaries around the model.

[0057] Based on the above embodiments, this method further includes:

[0058] S3. Embed a discrete element-based particle flow model within the preset borehole area of ​​the integrated model, and use a servo control method to make the particle flow model reach an initial stress equilibrium state to obtain a drill bit-formation coupled numerical model.

[0059] Specifically, step S3 includes:

[0060] S31. A wall is generated at the boundary of the preset drilling area to constrain the particles;

[0061] Specifically, for the preset drilling area, a square cylindrical wall matching the shape of the hole wall is generated, and the wall is used to define the spatial boundary of the particles. At the same time, a horizontal boundary wall is generated at the bottom of the preset drilling area.

[0062] S32. Generate a particle flow model composed of spherical particles within the drilling area, and preset the particle size distribution, density, and contact stiffness parameters of the spherical particles.

[0063] Specifically, random spherical particles are generated within the borehole area surrounded by a wall, and the particle size is set to follow a uniform distribution, ranging from 5 mm to 15 mm, to simulate the non-uniformity of the soil. Simultaneously, basic physical parameters of the spherical particles are set, including void ratio, maximum and minimum particle diameter, density, damping ratio, normal stiffness, and tangential stiffness. Among these, the tangential stiffness... .

[0064] S33. Define the interaction between spherical particles using a linear contact model;

[0065] S34. Obtain the actual stress of the wall and calculate the difference between the actual stress of the wall structure and the preset target stress:

[0066] ;

[0067] In the formula, Indicates the stress difference. Indicates the preset target stress. Represents the actual stress, where, , Indicates the area of ​​the wall subjected to force;

[0068] when At that time, the wall moves inward, achieving compaction loading; when At that time, the wall moves outward, achieving relaxation and unloading.

[0069] S35. Dynamically adjust the running speed of the wall based on the difference, and repeatedly calculate the difference between the actual stress of the wall and the preset target stress:

[0070] ;

[0071] In the formula, Indicates running speed.

[0072] The actual stress of the wall is updated through iterative calculation:

[0073] ;

[0074] In the formula, This represents the servo proportional gain coefficient, which controls the adjustment rate. This indicates the control stress to be applied at the next moment;

[0075] By repeatedly acquiring the difference between the current wall stress and the preset target stress, the running speed of the wall is dynamically adjusted to achieve automatic compaction and balance of the particle system.

[0076] In this embodiment, six walls are generated in the model: top and bottom, left and right, front and back, which serve as control boundaries.

[0077] Define a preset target stress for each wall, for example: vertical target stress: 200 kPa; horizontal template stress: 100 kPa.

[0078] S36. Until the difference is less than the preset stress tolerance, it is determined that the particle flow model has reached the initial stress equilibrium state;

[0079] Specifically, when When the stress tolerance is less than the preset tolerance and the system kinetic energy approaches zero, it indicates that the granular flow model has reached a servo equilibrium state. In the formula, This represents the stress tolerance. At this point, the particle density is most uniform, the porosity is stable, and the calculation model has the compactness characteristics of a real stratum.

[0080] After the granular flow model reaches stress equilibrium, the following steps are also included:

[0081] S37. Apply a parallel cementation model to the particle flow model that has reached the initial stress equilibrium state, and use the parallel cementation model to assign compressive strength, tensile strength and cementation stiffness parameters to the spherical particles. Specifically, the cementation stiffness parameters include normal cementation stiffness, tangential cementation stiffness and cementation radius multiplier.

[0082] S38. Perform equilibrium calculations on the particle flow model again to make the particle flow model reach a new stress equilibrium state in the cemented state.

[0083] S39. Reset the velocity and displacement of the spherical particles to zero as the initial state for the drilling simulation, ensuring that when the drilling simulation begins, the entire particle system is in an ideal initial state that is static, free from initial kinetic energy disturbance, and has real in-situ stress and intrinsic strength.

[0084] Based on the above embodiments, this method further includes:

[0085] S4. Perform numerical simulation of drilling rig drilling and square pile hole formation on the drill bit-formation coupled numerical model, and obtain the stress data of the drill bit and the deformation response data of the hole wall in real time.

[0086] Specifically, step S4 includes:

[0087] S41. Using the center point of the square pile as the geometric reference, delete the particle units within the range tangent to the outer circle of the square pile to form an initial pilot hole;

[0088] Specifically, in PFC3D, the radius of the circumcircle of the pre-defined drilling area is calculated. Taking the ground projection point of the square pile design center line as the center and the radius of the circumcircle as the range of action, all spherical particles in the cylindrical area are deleted.

[0089] S42. Multiple monitoring points are set up on the surface of the borehole wall along the pile axis direction of the square pile. Specifically, a monitoring section is set up every 0.5 meters along the depth direction on the surface of the initial pilot hole wall. On each section, a monitoring point is set up at each of the four corners of the square borehole wall and at the midpoint of each side to obtain the displacement, stress and inclination changes of the borehole wall in real time during the excavation process.

[0090] S43. Assign wall properties to the drill bit model, giving it contact characteristics that interact with particles, and control the drill bit model to rotate around its central axis and advance axially, resulting in a square drill hole as shown. Figure 3 As shown;

[0091] Specifically, motion control scripts are written using the FISH language to achieve synchronous control of the drilling rig's rotation around its central axis and axial thrust. By setting the rotational angular velocity and axial thrust velocity, the process of continuous cutting and fracturing by the drilling rig in the formation is simulated. Preferably, an elastic body model is used as the material model for the drilling rig.

[0092] S44. Real-time acquisition of contact stress distribution and torque data on the drill bit model surface, as well as displacement data at each monitoring point, such as... Figure 4 As shown;

[0093] Based on the above embodiments, this method further includes:

[0094] S5. Perform multi-parameter coupling analysis on the stress data of the drill bit and the deformation response data of the borehole wall to establish the coupling relationship between drilling rig cutting stress and borehole wall displacement;

[0095] ;

[0096] In the formula, This indicates the total radial displacement of the borehole wall. Indicates the initial displacement. Indicates the cutting stress of the drilling rig. All represent formation response parameters, obtained by fitting simulation data using the least squares method.

[0097] The overall radial displacement of the borehole wall is set with a critical threshold for borehole stability. When the overall radial displacement of the borehole wall is detected to be close to or exceed the critical threshold, the risk of borehole wall instability can be warned, thereby accurately judging the degree of influence of drill bit cutting stress on borehole wall stability.

[0098] Among them, drilling rig cutting stress The calculation method is as follows:

[0099] ;

[0100] In the formula, This indicates the axial thrust of the drill bit. This represents the cross-sectional area of ​​the drill bit. Indicates drill bit torque. Indicates the equivalent radius of the drill bit. This represents the lateral surface area of ​​the drill bit. This represents the torque contribution coefficient.

[0101] Among them, the comprehensive radial displacement of the borehole wall The calculation method is as follows:

[0102] ;

[0103] In the formula, This indicates the displacement of the hole wall in the X direction. This indicates the displacement of the hole wall in the Y direction.

[0104] Example 2:

[0105] like Figure 5As shown in the figure, this embodiment provides a numerical simulation device for simulating the hole-forming effect of square piles. The device includes:

[0106] The geological data acquisition module is used to acquire geological data of the target slope, establish a three-dimensional model of the slope including the preset borehole area based on the geological data, and spatially register the three-dimensional model of the slope with the independently constructed drill bit model to obtain an integrated model.

[0107] The mesh processing module is used to mesh the integrated model and import it into the numerical analysis software, assign material parameters to the slopes in the integrated model and set boundary conditions.

[0108] The particle flow modeling module is used to embed a discrete element-based particle flow model into a preset borehole area of ​​the integrated model, and to make the particle flow model reach an initial stress equilibrium state through servo control, thereby obtaining a drill bit-formation coupled numerical model.

[0109] The drilling simulation module is used to perform numerical simulation of drilling rig drilling and square pile hole formation on the drill bit-formation coupled numerical model, and to acquire the stress data of the drill bit and the deformation response data of the hole wall in real time.

[0110] The data analysis module is used to perform multi-parameter coupling analysis on the stress data of the drill bit and the deformation response data of the borehole wall to establish the coupling relationship between the drilling rig cutting stress and the borehole wall displacement.

[0111] Based on the above embodiments, the geological data acquisition module includes:

[0112] The laser scanning unit is used to scan the target slope using lidar to acquire three-dimensional terrain point cloud data;

[0113] The terrain processing unit is used to digitize the three-dimensional terrain point cloud data to extract the terrain morphology, slope, aspect and elevation information of the target slope.

[0114] A geological exploration unit is used to obtain layered information of underground rock and soil masses through drilling and geological exploration in order to determine the thickness data of each stratum.

[0115] Based on the above embodiments, the geological data acquisition module further includes:

[0116] The surface modeling unit is used to establish a three-dimensional model of the target slope based on the surface modeling method and the topographic morphology, slope, aspect and elevation information of the slope.

[0117] The stratigraphic interface generation unit is used to generate stratigraphic interfaces representing different soil and rock layers in the three-dimensional model of the slope based on the thickness data of each stratum.

[0118] The borehole area setting unit is used to cut out the preset borehole area in the slope 3D model of the integrated stratigraphic interface to obtain the final slope 3D model.

[0119] Based on the above embodiments, the mesh processing module includes:

[0120] The model grouping unit is used to group the slope and drill bit models of the integrated model in the numerical analysis software, and to divide each soil and rock layer according to the stratigraphic interface.

[0121] Material definition unit, used to define elastic constitutive models for each soil and rock layer, and assign corresponding elastic modulus and Poisson's ratio parameters respectively;

[0122] Boundary constraint elements are used to apply fixed constraints to the bottom boundary of the integrated model and roll constraints to the four lateral boundaries of the integrated model.

[0123] Based on the above embodiments, the granular flow modeling module includes:

[0124] A wall generation unit is used to generate a wall at the boundary of the preset drilling area to constrain particles;

[0125] The particle generation unit is used to generate a particle flow model composed of spherical particles in the borehole area, and preset the particle size distribution, density and contact stiffness parameters of the spherical particles.

[0126] Contact model unit, used to define the interaction between spherical particles using a linear contact model;

[0127] A stress monitoring unit is used to acquire the actual stress of the wall and calculate the difference between the actual stress of the wall structure and the preset target stress.

[0128] A servo control unit is used to dynamically adjust the running speed of the wall based on the difference and repeatedly calculate the difference between the actual stress of the wall and the preset target stress.

[0129] The equilibrium determination unit is used to determine that the particle flow model has reached the initial stress equilibrium state when the difference is less than the preset stress tolerance.

[0130] Based on the above embodiments, the granular flow modeling module further includes:

[0131] The cementation model unit is used to apply a parallel cementation model to the particle flow model that has reached the initial stress equilibrium state, and to assign compressive strength, tensile strength and cementation stiffness parameters to the spherical particles using the parallel cementation model.

[0132] The secondary balancing unit is used to perform balancing calculations on the particle flow model again, so that the particle flow model reaches a new stress equilibrium state in the cemented state.

[0133] The state reset unit is used to reset the velocity and displacement of the spherical particles to zero, which is the initial state of the drilling simulation.

[0134] Based on the above embodiments, the drilling simulation module includes:

[0135] The pilot hole forming unit is used to form an initial pilot hole by deleting the particle units within the range tangent to the outer circle of the square pile, with the center point of the square pile as the geometric reference.

[0136] The monitoring point deployment unit is used to deploy multiple monitoring points on the surface of the borehole wall along the pile axis direction of the square pile.

[0137] The drill bit control unit is used to assign wall properties to the drill bit model and control the rotation and axial advance of the drill bit model around the central axis.

[0138] The data acquisition unit is used to acquire in real time the contact stress distribution, contact force change and torque data of the drill bit model surface, as well as the displacement and stress time history data of each monitoring point.

[0139] It should be noted that the specific manner in which each module performs its operation in the apparatus described in the above embodiments has been described in detail in the embodiments of the method, and will not be elaborated here.

[0140] Example 3:

[0141] Corresponding to the above method embodiments, this embodiment also provides a numerical simulation device for simulating the hole-forming effect of square piles. The numerical simulation device for simulating the hole-forming effect of square piles described below and the numerical simulation method for simulating the hole-forming effect of square piles described above can be referred to in correspondence with each other.

[0142] Figure 6 This is a block diagram of a numerical simulation device 800 for simulating the drilling effect of square piles, according to an exemplary embodiment. Figure 6 As shown, the numerical simulation device 800 for simulating the hole-forming effect of square piles may include: a processor 801 and a memory 802. The numerical simulation device 800 may also include one or more of a multimedia component 803, an I / O interface 804, and a communication component 805.

[0143] The processor 801 controls the overall operation of the numerical simulation device 800 for simulating the hole-forming effect of the square pile, in order to complete all or part of the steps in the numerical simulation method for simulating the hole-forming effect of the square pile described above. The memory 802 stores various types of data to support the operation of the numerical simulation device 800 for simulating the hole-forming effect of the square pile. This data may include, for example, instructions for any application or method operating on the numerical simulation device 800 for simulating the hole-forming effect of the square pile, as well as application-related data, such as contact data, sent and received messages, images, audio, video, etc. The memory 802 can be implemented using any type of volatile or non-volatile storage device or a combination thereof, such as Static Random Access Memory (SRAM), Electrically Erasable Programmable Read-Only Memory (EEPROM), Erasable Programmable Read-Only Memory (EPROM), Programmable Read-Only Memory (PROM), Read-Only Memory (ROM), magnetic storage, flash memory, magnetic disk, or optical disk. The multimedia component 803 may include a screen and an audio component. The screen may be, for example, a touchscreen, and the audio component is used to output and / or input audio signals. For example, the audio component may include a microphone for receiving external audio signals. The received audio signals may be further stored in the memory 802 or transmitted via the communication component 805. The audio component also includes at least one speaker for outputting audio signals. I / O interface 804 provides an interface between processor 801 and other interface modules, such as keyboards, mice, and buttons. These buttons can be virtual or physical. Communication component 805 is used for wired or wireless communication between the numerical simulation device 800 simulating the hole-forming effect of the square pile and other devices. Wireless communication includes, for example, Wi-Fi, Bluetooth, Near Field Communication (NFC), 2G, 3G, or 4G, or a combination thereof. Therefore, the corresponding communication component 805 may include a Wi-Fi module, a Bluetooth module, or an NFC module.

[0144] In an exemplary embodiment, the numerical simulation device 800 for simulating the hole-forming effect of square piles may be implemented by one or more application-specific integrated circuits (ASICs), digital signal processors (DSPs), digital signal processing devices (DSPDs), programmable logic devices (PLDs), field-programmable gate arrays (FPGAs), controllers, microcontrollers, microprocessors, or other electronic components to perform the numerical simulation method for simulating the hole-forming effect of square piles described above.

[0145] In another exemplary embodiment, a computer-readable storage medium including program instructions is also provided, which, when executed by a processor, implement the steps of the numerical simulation method for simulating the hole-forming effect of square piles described above. For example, the computer-readable storage medium may be the memory 802 including program instructions described above, which may be executed by the processor 801 of the numerical simulation device 800 for simulating the hole-forming effect of square piles to complete the numerical simulation method for simulating the hole-forming effect of square piles described above.

[0146] Example 4:

[0147] Corresponding to the above method embodiments, this embodiment also provides a readable storage medium. The readable storage medium described below and the numerical simulation method for simulating the hole-forming effect of square piles described above can be referred to in correspondence.

[0148] A readable storage medium storing a computer program, wherein when the computer program is executed by a processor, the steps of the numerical simulation method for simulating the hole-forming effect of square piles in the above-described method embodiments are implemented.

[0149] Specifically, the readable storage medium can be a USB flash drive, a portable hard drive, a read-only memory (ROM), a random access memory (RAM), a magnetic disk, or an optical disk, or any other readable storage medium capable of storing program code.

[0150] The above description is merely a preferred embodiment of the present invention and is not intended to limit the invention. Various modifications and variations can be made to the present invention by those skilled in the art. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.

[0151] The above description is merely a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the technical scope disclosed in the present invention should be included within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be determined by the scope of the claims.

Claims

1. A numerical simulation method for simulating the hole-forming effect of square piles, characterized in that, include: Geological data of the target slope is acquired, a three-dimensional model of the slope including the preset borehole area is established based on the geological data, and the three-dimensional model of the slope is spatially registered with the independently constructed drill bit model to obtain an integrated model; After meshing the integrated model, it is imported into numerical analysis software to assign material parameters and set boundary conditions for the slopes in the integrated model. A discrete element-based particle flow model is embedded within the pre-defined borehole area of ​​the integrated model, and the particle flow model is brought to an initial stress equilibrium state through servo control to obtain a drill bit-formation coupled numerical model. Numerical simulations of drilling rig drilling and square pile hole formation were performed on the drill bit-formation coupled numerical model, and the stress data of the drill bit and the deformation response data of the hole wall were acquired in real time. Multi-parameter coupling analysis was performed on the stress data of the drill bit and the deformation response data of the borehole wall to establish the coupling relationship between the drilling rig cutting stress and the borehole wall displacement.

2. The numerical simulation method for simulating the hole-forming effect of square piles according to claim 1, characterized in that, The acquisition of geological data for the target slope includes: The target slope is scanned using lidar to obtain three-dimensional terrain point cloud data; The three-dimensional terrain point cloud data is digitally processed to extract the terrain morphology, slope, aspect and elevation information of the target slope; By drilling and sampling and geological exploration, we can obtain the layering information of underground rock and soil masses to determine the thickness data of each stratum.

3. The numerical simulation method for simulating the hole-forming effect of square piles according to claim 2, characterized in that, The process of establishing a three-dimensional slope model containing a pre-defined borehole area based on the geological data includes: Based on the surface modeling method, a three-dimensional model of the target slope is established according to the topographic morphology, slope, aspect and elevation information of the slope. Based on the thickness data of each stratum, stratigraphic interfaces representing different soil and rock layers are generated in the three-dimensional model of the slope. The pre-defined drilling area is cut out in the three-dimensional slope model of the integrated stratigraphic interface to obtain the final three-dimensional slope model.

4. The numerical simulation method for simulating the hole-forming effect of square piles according to claim 1, characterized in that, After meshing the integrated model and importing it into numerical analysis software, material parameters are assigned to the slopes of the integrated model and boundary conditions are set, including: In the numerical analysis software, the slope and drill bit models of the integrated model are grouped, and each soil and rock layer is divided according to the stratigraphic interface. Define an elastic constitutive model for each soil and rock layer, and assign corresponding elastic modulus and Poisson's ratio parameters respectively; Apply fixed constraints to the bottom boundary of the integrated model and roll constraints to the four lateral boundaries of the integrated model.

5. The numerical simulation method for simulating the hole-forming effect of square piles according to claim 1, characterized in that, The step of embedding a discrete element-based particle flow model within a preset borehole region of the integrated model, and using a servo control method to bring the particle flow model to an initial stress equilibrium state, includes: A wall is generated at the boundary of the preset drilling area to constrain the particles; A particle flow model composed of spherical particles is generated within the drilling area, and the particle size distribution, density, and contact stiffness parameters of the spherical particles are preset. The interaction between spherical particles is defined using a linear contact model; Obtain the actual stress of the wall and calculate the difference between the actual stress of the wall structure and the preset target stress; The running speed of the wall is dynamically adjusted based on the difference, and the difference between the actual stress of the wall and the preset target stress is repeatedly calculated. The particle flow model is considered to have reached an initial stress equilibrium state when the difference is less than the preset stress tolerance.

6. A numerical simulation device for simulating the hole-forming effect of square piles, characterized in that, include: The geological data acquisition module is used to acquire geological data of the target slope, establish a three-dimensional model of the slope including the preset borehole area based on the geological data, and spatially register the three-dimensional model of the slope with the independently constructed drill bit model to obtain an integrated model. The mesh processing module is used to mesh the integrated model and import it into the numerical analysis software, assign material parameters to the slopes in the integrated model and set boundary conditions. The particle flow modeling module is used to embed a discrete element-based particle flow model into a preset borehole area of ​​the integrated model, and to make the particle flow model reach an initial stress equilibrium state through servo control, thereby obtaining a drill bit-formation coupled numerical model. The drilling simulation module is used to perform numerical simulation of drilling rig drilling and square pile hole formation on the drill bit-formation coupled numerical model, and to acquire the stress data of the drill bit and the deformation response data of the hole wall in real time. The data analysis module is used to perform multi-parameter coupling analysis on the stress data of the drill bit and the deformation response data of the borehole wall to establish the coupling relationship between the drilling rig cutting stress and the borehole wall displacement.

7. The numerical simulation device for simulating the hole-forming effect of square piles according to claim 6, characterized in that, The geological data acquisition module includes: The laser scanning unit is used to scan the target slope using lidar to acquire three-dimensional terrain point cloud data; The terrain processing unit is used to digitize the three-dimensional terrain point cloud data to extract the terrain morphology, slope, aspect and elevation information of the target slope. A geological exploration unit is used to obtain layered information of underground rock and soil masses through drilling and geological exploration in order to determine the thickness data of each stratum.

8. The numerical simulation device for simulating the hole-forming effect of square piles according to claim 7, characterized in that, The geological data acquisition module also includes: The surface modeling unit is used to establish a three-dimensional model of the target slope based on the surface modeling method and the topographic morphology, slope, aspect and elevation information of the slope. The stratigraphic interface generation unit is used to generate stratigraphic interfaces representing different soil and rock layers in the three-dimensional model of the slope based on the thickness data of each stratum. The borehole area setting unit is used to cut out the preset borehole area in the slope 3D model of the integrated stratigraphic interface to obtain the final slope 3D model.

9. The numerical simulation device for simulating the hole-forming effect of square piles according to claim 6, characterized in that, The mesh processing module includes: The model grouping unit is used to group the slope and drill bit models of the integrated model in the numerical analysis software, and to divide each soil and rock layer according to the stratigraphic interface. Material definition unit, used to define elastic constitutive models for each soil and rock layer, and assign corresponding elastic modulus and Poisson's ratio parameters respectively; Boundary constraint elements are used to apply fixed constraints to the bottom boundary of the integrated model and roll constraints to the four lateral boundaries of the integrated model.

10. The numerical simulation device for simulating the hole-forming effect of square piles according to claim 6, characterized in that, The granular flow modeling module includes: A wall generation unit is used to generate a wall at the boundary of the preset drilling area to constrain particles; The particle generation unit is used to generate a particle flow model composed of spherical particles in the borehole area, and preset the particle size distribution, density and contact stiffness parameters of the spherical particles. Contact model unit, used to define the interaction between spherical particles using a linear contact model; A stress monitoring unit is used to acquire the actual stress of the wall and calculate the difference between the actual stress of the wall structure and the preset target stress. A servo control unit is used to dynamically adjust the running speed of the wall based on the difference and repeatedly calculate the difference between the actual stress of the wall and the preset target stress. The equilibrium determination unit is used to determine that the particle flow model has reached the initial stress equilibrium state when the difference is less than the preset stress tolerance.

Images