Analysis method and device of hole stress field, computer equipment, readable storage medium and program product

Abstract

The application relates to a hole stress field analysis method and device, computer equipment, a readable storage medium and a program product. The method comprises the following steps: for each hole in a target area, determining plane mapping data based on a plane conformal mapping algorithm, anisotropy characteristic data and hole shape parameters, combining hole edge stress data to construct a uniform stress distribution model; determining first stress data of a to-be-measured point under the action of unit load of each target arc based on the model; determining second stress data of the to-be-measured point based on remote stress field data; processing the first and second stress data through a force analysis algorithm to determine interaction force data; and analyzing the remote stress field data, the first stress data and the interaction force data based on a stress analysis algorithm to obtain target stress field data of the to-be-measured point. The method can improve the efficiency and accuracy of stress field calculation in anisotropic rock layers.

Description

Technical Field

[0001] This application relates to the field of fracture mechanics technology, and in particular to a method, apparatus, computer equipment, readable storage medium, and program product for analyzing the stress field of a cavity. Background Technology

[0002] With the continuous advancement of urban infrastructure construction such as subways, highways, and water conservancy projects, the twin-tunnel structure has been widely applied in engineering practice. Compared with single tunnels, twin-tunnel structures have significant advantages in saving underground space and improving transportation efficiency. However, underground rock strata generally exhibit anisotropy and heterogeneity, especially in weak rock masses such as shale, slate, and interbedded sandstone. These rock strata often have quasi-orthogonal anisotropy characteristics and contain natural joints, fissures, and other structural surfaces. During tunnel excavation, these geological discontinuities can cause a redistribution of the original stress field. Especially in the presence of multiple or adjacent pores, stress concentration can easily be induced in local areas, leading to crack initiation and propagation. In severe cases, this can cause surrounding rock instability or support structure failure, seriously affecting the safety and durability of the tunnel project. Furthermore, in deep tunnel underground engineering, the surrounding rock commonly contains natural or construction-induced pore and fissure structures. Under the combined action of far-field stress and arbitrarily distributed seepage pressure within the pores, high-stress zones can easily form around the pores, becoming potential sources of damage.

[0003] Actual monitoring results show that, in addition to geostress, the dominant far-field factor, the pores often contain non-uniformly distributed pore pressure, seepage disturbances, and construction loads, among other near-field influencing factors. The superposition of these loads significantly alters the stress state and stability of the surrounding rock mass, increasing the uncertainty and risk of the project. When dealing with arbitrary loads in orthotropic porous media, the hypersingular boundary element method is commonly used in related technologies, resulting in low efficiency in analyzing the tunnel stress field. Summary of the Invention

[0004] Therefore, it is necessary to provide a method, apparatus, computer equipment, readable storage medium, and program product for analyzing pore stress fields in anisotropic rock strata, which can improve the efficiency and accuracy of stress field calculation in anisotropic rock strata.

[0005] Firstly, this application provides a method for analyzing the stress field of a cavity, including:

[0006] For each hole in the target area, based on the planar conformal mapping algorithm, anisotropic feature data and the shape parameters of the hole, the planar mapping data corresponding to the hole is determined, and based on the planar mapping data and the hole edge stress data, the uniform stress distribution model of the hole is determined.

[0007] Based on the uniform stress distribution model, the first stress data of the test point under each segment of the target arc in the hole under a unit load is determined; based on the far-end stress field data, the second stress data of the test point under the far-end stress is determined.

[0008] Based on the force analysis algorithm, the first stress data and the second stress data are processed to determine the interaction force data of the point to be measured;

[0009] Based on the stress analysis algorithm, the remote stress field data, the first stress data, and the interaction force data are analyzed to obtain the target stress field data of the test point.

[0010] In one embodiment, determining the uniform stress distribution model of the hole based on the planar mapping data and the hole edge stress data includes:

[0011] Based on the stress data at the edge of the hole, the load data of the hole is determined; based on the load data and the anisotropic characteristic data, a first coefficient is calculated; and based on the first coefficient and the plane mapping data, a first potential function component of the hole is determined.

[0012] Based on the load data, the mapping position data of each target arc segment, and the shape parameters of the hole, the load boundary data of each target arc segment of the hole is determined; based on the first coefficient and the load boundary data, the second potential function components of each target arc segment of the hole are determined.

[0013] Based on the first potential function component and the second potential function component, a uniform stress distribution model for the hole is determined.

[0014] In one embodiment, determining the first stress data of the target arc segment of the hole under a unit load based on the uniform stress distribution model includes:

[0015] For each target arc, the uniform stress distribution data of the target arc is obtained by analyzing the unit load data and the position data of each target arc based on the uniform stress distribution model.

[0016] Based on the anisotropic characteristic data and the uniform stress distribution data, the first stress data of the target arc under the test point under unit load is calculated.

[0017] In one embodiment, the process of processing each of the first stress data and each of the second stress data based on the force analysis algorithm to determine the interaction force data of the test point includes:

[0018] Based on the stress data at the edge of the hole and the second stress data, the external force data of each target arc segment of the hole is determined;

[0019] Based on the force balance algorithm, the external force data, the first stress data of each target arc of each hole, and the shape parameters of each target arc, the interaction force data of the point to be measured is determined.

[0020] In one embodiment, the step of analyzing the remote stress field data, the first stress data, and the interaction force data based on a stress analysis algorithm to obtain the target stress field data of the test point includes:

[0021] Based on the remote stress field data and the position coordinate data of the point to be measured, the first stress field component of the point to be measured under the action of the remote stress is determined.

[0022] Based on the first stress data corresponding to each target arc of each hole and the interaction force data of each target arc, the second stress field component of the test point under the action of each target arc of each hole is determined.

[0023] Based on the first stress field component and the second stress field component, the target stress field data of the point to be measured is determined.

[0024] In one embodiment, the method further includes:

[0025] Based on the mechanical response data of the rock strata in the target area, the equivalent compliance coefficient is determined;

[0026] The equivalent orthogonal flexibility coefficient is calculated based on the equivalent flexibility coefficient, and the anisotropic direction characteristic parameters are determined based on the equivalent orthogonal flexibility coefficient.

[0027] Secondly, this application also provides an analysis device for the stress field of a cavity, comprising:

[0028] The first determining module is used to determine the plane mapping data corresponding to each hole in the target area based on the plane conformal mapping algorithm, anisotropic feature data and the shape parameters of the hole, and to determine the uniform stress distribution model of the hole based on the plane mapping data and the hole edge stress data.

[0029] The second determining module is used to determine, based on the uniform stress distribution model, the first stress data of the target arc segment of the hole under unit load at the point to be measured; and to determine the second stress data of the point to be measured under remote stress based on the remote stress field data.

[0030] The third determining module is used to process the first stress data and the second stress data based on the force analysis algorithm to determine the interaction force data of the point to be measured.

[0031] The analysis module is used to analyze the remote stress field data, the first stress data, and the interaction force data based on the stress analysis algorithm to obtain the target stress field data of the test point.

[0032] Thirdly, this application also provides a computer device, including a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to perform the following steps:

[0033] For each hole in the target area, based on the planar conformal mapping algorithm, anisotropic feature data and the shape parameters of the hole, the planar mapping data corresponding to the hole is determined, and based on the planar mapping data and the hole edge stress data, the uniform stress distribution model of the hole is determined.

[0034] Based on the uniform stress distribution model, the first stress data of the test point under each segment of the target arc in the hole under a unit load is determined; based on the far-end stress field data, the second stress data of the test point under the far-end stress is determined.

[0035] Based on the force analysis algorithm, the first stress data and the second stress data are processed to determine the interaction force data of the point to be measured;

[0036] Based on the stress analysis algorithm, the remote stress field data, the first stress data, and the interaction force data are analyzed to obtain the target stress field data of the test point.

[0037] Fourthly, this application also provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, performs the following steps:

[0038] For each hole in the target area, based on the planar conformal mapping algorithm, anisotropic feature data and the shape parameters of the hole, the planar mapping data corresponding to the hole is determined, and based on the planar mapping data and the hole edge stress data, the uniform stress distribution model of the hole is determined.

[0039] Based on the uniform stress distribution model, the first stress data of the test point under each segment of the target arc in the hole under a unit load is determined; based on the far-end stress field data, the second stress data of the test point under the far-end stress is determined.

[0040] Based on the force analysis algorithm, the first stress data and the second stress data are processed to determine the interaction force data of the point to be measured;

[0041] Based on the stress analysis algorithm, the remote stress field data, the first stress data, and the interaction force data are analyzed to obtain the target stress field data of the test point.

[0042] Fifthly, this application also provides a computer program product, including a computer program that, when executed by a processor, performs the following steps:

[0043] For each hole in the target area, based on the planar conformal mapping algorithm, anisotropic feature data and the shape parameters of the hole, the planar mapping data corresponding to the hole is determined, and based on the planar mapping data and the hole edge stress data, the uniform stress distribution model of the hole is determined.

[0044] Based on the uniform stress distribution model, the first stress data of the test point under each segment of the target arc in the hole under a unit load is determined; based on the far-end stress field data, the second stress data of the test point under the far-end stress is determined.

[0045] Based on the force analysis algorithm, the first stress data and the second stress data are processed to determine the interaction force data of the point to be measured;

[0046] Based on the stress analysis algorithm, the remote stress field data, the first stress data, and the interaction force data are analyzed to obtain the target stress field data of the test point.

[0047] The aforementioned analysis method, apparatus, computer equipment, readable storage medium, and program product for the stress field of holes, for each hole in the target area, determines the corresponding plane mapping data of the hole based on a plane conformal mapping algorithm, anisotropic characteristic data, and the shape parameters of the hole, and determines the uniform stress distribution model of the hole based on the plane mapping data and the hole edge stress data; it realizes the transformation of the anisotropic porous stress field problem, so that the complex stress boundary conditions originally applied on the curved surface of the circular hole are transformed into standard boundary conditions on a unit circumference, improving the efficiency of boundary condition solution; based on the uniform stress distribution model, it determines the first stress number of the target arc segment of the hole at the test point under a unit load. Based on the remote stress field data, the second stress data of the test point under the action of remote stress is determined; a bridge is established from the stress response of a single hole to the analysis of multi-hole interactions; based on the force analysis algorithm, the first stress data and the second stress data are processed to determine the interaction force data of the test point; the accurate direct solution of multi-hole interaction forces is realized; based on the stress analysis algorithm, the remote stress field data, the first stress data and the interaction force data are analyzed to obtain the target stress field data of the test point; the efficiency and accuracy of solving the stress field at any position in the solution domain are realized, the applicable range of multi-hole interactions in anisotropic materials is expanded, and the stability of safety assessment and analysis of multi-hole tunnels is improved. Attached Figure Description

[0048] To more clearly illustrate the technical solutions in the embodiments of this application or related technologies, the drawings used in the description of the embodiments of this application or related technologies will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other related drawings can be obtained based on these drawings without creative effort.

[0049] Figure 1 This is a flowchart illustrating the method for analyzing the stress field of a hole in one embodiment;

[0050] Figure 2 This is a flowchart illustrating step 101 in one embodiment;

[0051] Figure 3 This is a flowchart illustrating step 104 in one embodiment;

[0052] Figure 4 This is a schematic diagram of a portion of an anisotropic rock stratum containing a single circular hole subjected to a uniformly distributed force at the edge of the hole in one embodiment.

[0053] Figure 5 This is a schematic diagram of the superposition of various interaction forces at point A in one embodiment;

[0054] Figure 6 This is a schematic diagram of an anisotropic rock layer containing two unequal-sized circular holes under far-field biaxial compressive stress and surface non-uniform stress in one embodiment.

[0055] Figure 7 This is a structural block diagram of an analysis device for the stress field of a hole in one embodiment;

[0056] Figure 8 This is an internal structural diagram of a computer device in one embodiment. Detailed Implementation

[0057] To make the objectives, technical solutions, and advantages of this application clearer, the following detailed description is provided in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the scope of this application.

[0058] In one exemplary embodiment, such as Figure 1 As shown, a method for analyzing the stress field of a hole is provided. This embodiment illustrates the application of this method to a terminal. It is understood that this method can also be applied to a server, and to a system including both a terminal and a server, and can be implemented through interaction between the terminal and the server. In this embodiment, the method includes the following steps:

[0059] Step 101: For each hole in the target area, based on the planar conformal mapping algorithm, anisotropic feature data and the shape parameters of the hole, determine the planar mapping data corresponding to the hole, and based on the planar mapping data and the hole edge stress data, determine the uniform stress distribution model of the hole.

[0060] The target region is the area where the stress field is to be analyzed. This target region includes multiple rock layers characterized by anisotropy and heterogeneity. Optionally, these rock layers can be underground rock layers, or soft rock masses such as shale, slate, and interbedded sandstone. The cavities are hollow structures present in the target region. Each cavity is located within the multiple rock layers of the target region. It should be understood that in this embodiment, the stress field analysis is performed on the plane stress surface of the target region, analyzing the stress field data at each point on that plane. Optionally, the rock layers include multiple different loading directions, with stress and strain existing in each loading direction. The loading directions can include directions parallel to bedding planes and directions perpendicular to bedding planes. Bedding planes refer to the parallel or near-parallel layered structures formed during sedimentation and diagenesis. Optionally, the parallel bedding plane direction is set as the x-direction, and the perpendicular bedding plane direction is set as... The y-direction is defined as follows: the plane containing the x and y directions is the xy plane; the plane conformal mapping algorithm is used to map holes on the plane stress surface to a unit circle of a specified complex plane; optionally, the plane conformal mapping algorithm is the conformal mapping function from the complex variable plane to the unit circle plane in the complex variable function method of anisotropic elasticity; anisotropic characteristic data is a set of parameters of different mechanical properties of rock strata in different directions, used to quantify the anisotropic characteristics of rock strata; shape parameters are data used to describe the geometry of the hole, for example, the shape parameters may include the radius of the rock stratum plane where the hole is located; plane mapping data refers to the mapping function data mapped to a specified complex plane; the hole edge is the edge formed by the hole and the rock strata, and the hole edge stress data is the load data of the hole under construction; optionally, due to the anisotropy of the rock strata, its mechanical properties change with the direction of force. The uniform stress distribution model can be the potential function of a single hole in the rock strata (anisotropic medium) under uniform stress at the hole edge, which is a mathematical model used to describe the stress field distribution law of the hole under partially uniform stress. The plane mapping data can be the complex variable plane in this anisotropic elasticity.

[0061] Specifically, for each hole in the target area, the terminal calculates the anisotropic feature data and the shape parameters of the hole based on the planar conformal mapping algorithm to obtain the planar mapping data corresponding to the hole; based on the hole edge stress data, the load data of the hole is determined; based on the load data and anisotropic feature data, a first coefficient is calculated; and based on the first coefficient and the planar mapping data, the first potential function component of the hole is determined; based on the load data, the mapping position data of each target arc segment, and the shape parameters of the hole, the load boundary data of each target arc segment of the hole is determined; based on the first coefficient and the load boundary data, the second potential function component of each target arc segment of the hole is determined; and based on the first potential function component and the second potential function component, the uniform stress distribution model of the hole is determined.

[0062] Optionally, the specific expression for calculating the plane mapping data based on the plane function mapping algorithm can be:

[0063]

[0064] in, It is the unit circle plane after the j-th mapping, satisfying The boundary of the circular hole corresponding to the original plane; R j This is the scaling factor of the mapping, which is directly related to the radius of the circular hole corresponding to the hole, and determines the scale of the geometric dimensions after mapping, j=1,2; m j It controls the shape correction factor, which is related to the material anisotropy parameters. Adjusting the mapped shape ensures that the unit circle boundary precisely corresponds to the circular hole boundary on the original plane. Optionally, , u j is the j-th anisotropic feature data, and a is the shape parameter, i.e., the radius of the hole.

[0065] Step 102: Based on the uniform stress distribution model, determine the first stress data of the target arc segment of the hole under unit load; based on the far-end stress field data, determine the second stress data of the target point under far-end stress.

[0066] The first stress data refers to the stress induced at the test point by a uniform force in a circular hole formed by a borehole, including radial, circumferential, and tangential stresses. The test point is any point in the space containing the rock strata; it can be a point on the edge of the borehole or any point in the rock strata not on the edge of the borehole. The target arc is the arc formed by dividing the edge of the borehole. The first stress data is the stress response data experienced by the test point under a unit load. The far-field stress field data refers to the far-field stress of the surrounding rock at an infinitely far (or sufficiently far) distance from the borehole in an anisotropic medium without boreholes. The second stress data refers to the force generated at the test point under the action of the far-field stress field, including tangential and radial forces experienced by the test point.

[0067] Specifically, the terminal can determine the first stress data of the test point under a unit load on each segment of the target arc of the hole based on a uniform stress distribution model; and acquire the remote stress field data; and calculate the second stress data of the test point under the remote stress based on the local polar angle of the test point in the corresponding coordinate system and the remote stress field data. Optionally, when the test point is any point on the edge of a hole, the corresponding coordinate system is a coordinate system (global coordinate system) established with the center of the hole as the origin; when the test point is any point on the rock layer other than any point on the edge of the hole, the corresponding coordinate system is a coordinate system established with the center of the rock layer as the center.

[0068] For example, the calculation expression for the second stress data can be:

[0069]

[0070] in, The point to be measured The radial force acting on the i-th point (on the k-th hole). It is the tangential force acting on the point to be measured. It is the local polar angle of the point to be measured; It is the horizontal normal stress in the global coordinate system; It is a normal stress in the vertical direction; It is far-field shear stress.

[0071] Step 103: Based on the force analysis algorithm, process the first stress data and the second stress data to determine the interaction force data of the point to be measured.

[0072] Among them, the force analysis algorithm is an algorithm for analyzing the internal and external forces acting on the test point; the interaction force data is the additional surface force data generated on each test point due to the presence of other holes when multiple holes coexist, which can be the radial interaction force and tangential interaction force of the test point.

[0073] Specifically, the terminal determines the external force data of each target arc segment of the hole based on the hole edge stress data and each second stress data. The external force data includes the equivalent force generated by the far-field stress at the hole edge and the combined effect of the known stress at the hole edge; based on the force balance algorithm, the external force data, the first stress data of each target arc segment of each hole, and the shape parameters of each target arc, the interaction force data of the measured point is determined.

[0074] Step 104: Based on the stress analysis algorithm, analyze the remote stress field data, the first stress data, and the interaction force data to obtain the target stress field data of the test point.

[0075] Among them, the stress analysis algorithm is used to calculate the total stress data at any point in the porous structure; the target stress field data is the data of multiple stress components experienced by the test point in the target region, including radial stress, circumferential stress and shear stress.

[0076] Specifically, based on the far-end stress field data and the position coordinate data of the test point, the first stress field component of the test point under the action of far-end stress is determined. In a homogeneous anisotropic medium without pores, the first stress field component is equal to the value of the far-end stress at the test point. Based on the first stress data corresponding to each target arc of each pore and the interaction force data of each target arc, the second stress field component of the test point under the action of each target arc of each pore is determined. Specifically, the interaction force of each target arc is multiplied by the corresponding kernel function, and then multiplied by the arc length of that arc to obtain the stress contribution generated by that arc at the test point. The stress contributions of all target arcs are summed to obtain the total second stress field component. Based on the first and second stress field components, the target stress field data of the test point is determined.

[0077] The aforementioned method for analyzing the stress field of holes, targeting each hole in the target area, determines the corresponding plane mapping data for the hole based on a planar conformal mapping algorithm, anisotropic characteristic data, and the shape parameters of the hole. Based on the plane mapping data and the hole edge stress data, a uniform stress distribution model for the hole is determined. This transforms the anisotropic porous stress field problem, converting the complex stress boundary conditions originally applied to the curved surface of the circular hole into standard boundary conditions on a unit circle, thus improving the efficiency of boundary condition solving. Based on the uniform stress distribution model, the method determines the first stress data experienced by the measured point on each segment of the target arc within the hole under a unit load. Based on the far-end stress field data... This study determined the second stress data of the test point under the action of far-end stress; established a bridge from single-hole stress response to multi-hole interaction analysis; processed the first and second stress data based on the force analysis algorithm to determine the interaction force data of the test point; realized the accurate direct solution of multi-hole interaction force; and analyzed the far-end stress field data, first stress data, and interaction force data based on the stress analysis algorithm to obtain the target stress field data of the test point. It improved the efficiency and accuracy of solving the stress field at any location within the solution domain, expanded the applicability of multi-hole interaction in anisotropic materials, and improved the stability of safety assessment and analysis of multi-hole tunnels.

[0078] In one exemplary embodiment, the method for analyzing the stress field of a hole further includes:

[0079] Based on the mechanical response data of the rock strata in the target area, the equivalent flexibility coefficient is determined; the equivalent orthogonal flexibility coefficient is calculated based on the equivalent flexibility coefficient; and the anisotropic direction characteristic parameters are determined based on the equivalent orthogonal flexibility coefficient.

[0080] Among them, the mechanical response data of the rock strata refers to the deformation feedback data generated under stress loading in different directions. This deformation feedback data can be strain data. The mechanical response data includes stress data and corresponding strain data under different loading directions. Stress data refers to the external force data borne by the rock strata, and strain data is the deformation feedback data of the rock strata in response to the force. The equivalent compliance coefficient characterizes the ease of deformation of the rock strata. The equivalent orthogonal compliance coefficient is the compliance coefficient after removing the coupling terms of shear and tension. The anisotropic direction characteristic parameters characterize the degree of anisotropy of the rock strata.

[0081] Specifically, the terminal can acquire mechanical response data of the rock strata in the target area obtained through multi-directional uniaxial compression experiments, and calculate the equivalent compliance coefficient based on the mechanical response data. Optionally, the mechanical response data may include multiple sets of strain and stress data. Each set of strain and stress data includes normal stress-strain data in the x-direction, stress-strain data in the y-direction, and tangential force-tangential strain data in the xy-plane. The expression for calculating the equivalent compliance coefficient can be:

[0082]

[0083] in, These are strain data in the x-direction, y-direction, and xy-plane. These are stress data in the x-direction, y-direction, and xy-plane. It is the equivalent compliance coefficient under stress in the xy plane.

[0084] The terminal can calculate the equivalent orthogonal compliance coefficient based on the relationship between stress and strain data, as well as the equivalent compliance coefficient. The equivalent orthogonal compliance coefficient can be represented as follows:

[0085]

[0086] in, These are the compliance coefficients in different directions.

[0087] Based on the expressions for equivalent orthogonality coefficients and equivalent compliance coefficients, an expression for calculating the anisotropic directional characteristic parameters is constructed. This expression can be:

[0088]

[0089] Solving the above calculation expression, we get:

[0090]

[0091] Wherein, coefficient u k , (k=1,2) are the roots of the fourth-order characteristic equation, and are the anisotropic direction characteristic parameters.

[0092] In this embodiment, the anisotropic directional characteristic parameters are determined by using the mechanical response data of the rock strata in the target area, thereby enabling the analysis of the anisotropic characteristics of the rock strata, extracting the directional characteristics of the rock strata's mechanical properties, and achieving accurate identification of the rock strata.

[0093] In one exemplary embodiment, such as Figure 2 As shown, the specific implementation process of "determining the uniform stress distribution model of the hole based on the plane mapping data and the hole edge stress data" in step 101 may include the following steps:

[0094] Step 201: Based on the stress data at the edge of the hole, determine the load data of the hole; based on the load data and anisotropic characteristic data, calculate the first coefficient; and based on the first coefficient and plane mapping data, determine the first potential function component of the hole.

[0095] Among them, the stress data around the hole refers to the stress distribution in the local area around the hole in a solid structure containing holes (such as rock strata, plates, and rods) when subjected to external forces. Load data refers to the external forces acting on the hole. Load data can include the composite vector of radial and tangential forces in the x and y directions. The first coefficient is related to the load data and the anisotropic characteristics of the material, and is the core source of stress concentration around the hole. The first potential function component reflects the disturbance of the hole to the surrounding stress field.

[0096] Specifically, the terminal determines the load data of the hole based on the stress data at the edge of the hole, including radial force and tangential force; calculates the resultant force in the x and y directions based on the radial force and tangential force; solves for the first coefficient based on the anisotropic characteristic data through the resultant force equilibrium condition; and determines the first potential function component of the hole based on the first coefficient and the plane mapping data.

[0097] The expression for calculating the first coefficient can be:

[0098]

[0099] Among them, A1, A2, , The first coefficient; u1, u2, , Characteristic roots are used to determine the elastic constants of anisotropic materials, reflecting the differences in the material's mechanical response in different directions. X is the radial force, and Y is the tangential force.

[0100] The product of the first coefficient and the logarithmic term of the plane mapping data is determined as the first potential function component.

[0101] Step 202: Based on the load data, the mapping position data of each target arc segment, and the shape parameters of the hole, determine the load boundary data of each target arc segment of the hole; based on the first coefficient and the load boundary data, determine the second potential function components of each target arc segment of the hole.

[0102] Among them, the mapped position data are the corresponding coordinates of the two endpoints of the target arc on the ζ-plane after conformal mapping. The load boundary data are the force boundary conditions on each arc segment of the hole edge on the ζ-plane; the second potential function component is used to describe the uniform stress field in the hole-free region, ensuring unique single-valuedness after eliminating singular terms. The load data refers to the radial and tangential surface forces at the hole edge.

[0103] Specifically, for each target arc, the corresponding position data of the two endpoints on the ζ-plane through the mapping function are calculated. The load data is integrated on the target arc to obtain the load boundary data in the x and y directions (load boundary data of each target arc of the hole). Using the Schwarz formula, the load boundary data is integrated on the unit circle with the load boundary data as the integrand. The integration path starts from the corresponding position data of the two endpoints on the ζ-plane and ends to obtain the second potential function component.

[0104] Optionally, the load data P (normal) and Q (tangential) at the hole edge are integrated along the hole wall to obtain two real functions f1(θ) and f2(θ). f1(θ) is the negative resultant force in the normal (y-direction) direction of the force on the arc segment corresponding to angle θ, integrated from each starting point along the hole wall. f2(θ) is the corresponding tangential resultant force. The specific boundary function expression can be:

[0105]

[0106] Where Yn is the component of the surface force at the hole edge in the normal direction, and P is the radial surface force (normal load) applied at the hole edge; Xn and Yn are the composite vectors of radial and tangential forces (P, Q) in different segments in the x and y directions; a is the radius of the circular hole; σ is a point on the unit circle of the auxiliary plane (ζ plane), σ=e^(iθ), and θ is the polar angle of the current integration point; σ1 is the coordinate of the starting point z1 of the integration arc segment on the unit circle, σ1=e^(iθ1); σ2 is the coordinate of the ending point z2 of the integration arc segment on the unit circle, σ2 = e^(iθ2); s is the arc length parameter along the hole wall.

[0107] Construct parameters m1+n1i=A1, m2+n2i=A2. According to complex function theory, Φj0(zj) has the following form. Construct a system of boundary condition equations. The expressions for this system of boundary condition equations can be:

[0108]

[0109] in, It is the modified normal boundary function, which serves as the known input for the Schwarz formula integral; This is the corrected tangential boundary function. The specific calculation expression can be:

[0110]

[0111]

[0112] Integrating the above equation using Schwarz's formula yields expressions for two holomorphic functions:

[0113]

[0114] in, It is the holomorphic part of the first complex potential function, defined outside the unit circle (| | > 1), It is the independent variable; ln(σ) A - ) is a logarithmic singular term, when →σ A The time tends towards infinity, representing the stress concentration effect at the endpoint A of the load arc; ln(σ B - ) is a logarithmic singular term, representing the stress concentration effect at the endpoint B of the load arc segment; b is a geometric parameter related to the mapping function, and the relationship between b and a in the standard circular mapping depends on the specific mapping form.

[0115] It should be understood that It is a logarithmic fundamental term related to the overall angular span of the arc segment, describing the far-field effect of the load arc segment as a whole.

[0116] It describes the local stress response generated at the starting point A of the load arc segment; It describes the local stress response generated at the end point B of the load arc segment.

[0117]

[0118] in, It is the holomorphic part of the second complex potential function.

[0119] The expression for calculating the potential function Φj(zj) is as follows:

[0120]

[0121]

[0122] in, , It is the second potential function component; z1 is the first generalized physical coordinate, z1 = x + μ1y, where x and y are physical rectangular coordinates, and μ1 and μ2 are anisotropic eigenvalues; ζ(z1) is the inverse function of the conformal mapping, mapping points in the z1 plane back to... The plane. z2 is the second generalized physical coordinate.

[0123] Step 203: Determine the uniform stress distribution model of the hole based on the first potential function component and the second potential function component.

[0124] Specifically, the sum of the first and second potential function components is defined as the uniform stress distribution model of the hole, and the specific expression can be:

[0125]

[0126] In this embodiment, by using conformal mapping and the Schwarz formula, the complex stress boundary conditions of circular holes in anisotropic media are transformed into an integral problem on a standard unit circle, resulting in an analytical potential function with a clear structure and high degree of modularity. This supports stress conditions with arbitrary distribution along the hole edges, providing a reliable mathematical foundation for subsequent kernel function construction and multi-hole superposition calculations.

[0127] In an exemplary embodiment, the specific implementation process of step 102, "based on a uniform stress distribution model, determining the first stress data of the target arc segment of the hole under a unit load at the test point," may include:

[0128] For each target arc, the uniform stress distribution data of the target arc is obtained by analyzing the unit load data and the position data of each target arc based on the uniform stress distribution model; based on the anisotropic characteristic data and the uniform stress distribution data, the first stress data of the target arc at the test point under the unit load is calculated.

[0129] The target arc position data consists of polar coordinates describing the target arc's location on the hole boundary. The uniform stress distribution data is the stress distribution data when the target arc acts as a stress source under a unit load.

[0130] Specifically, for each target arc segment, the uniform stress distribution data is obtained by analyzing the uniform stress distribution model, unit load data, and position data of each target arc. Two unit load conditions are applied: the first condition sets the unit radial force P=1 and the unit tangential force Q=0; the second condition sets the unit radial force P=0 and the unit tangential force Q=1. For each condition, the stress response data of the test point is calculated using the uniform stress distribution model. Based on the anisotropic characteristic data and the aforementioned uniform stress distribution data, the first stress data of the test point under the unit load on each target arc segment is calculated. The first stress data includes six kernel functions, which describe the effects of the unit radial force and unit tangential force on the radial stress, circumferential stress, and shear stress of the test point.

[0131] Optionally, the stress σ induced at any point by a single circular hole under partially uniform stress can be calculated. r , σ θ ,τ rθ and far-field stress , and The stress σ induced at any point by far-field stress in a non-porous anisotropic medium is calculated. r , σ θ , τ rθ The expression can be:

[0132]

[0133] Where, σ r , σ θ , τ rθ These represent the radial, circumferential, and tangential stresses at any given point.

[0134] Alternatively, setting P=1 and Q=0 yields the kernel function (G) of the interaction of the unit radial force with respect to any point. rt (z1z2,z),G rs (z1z2,z),G rr (z1z2,z)); Setting P=0 and Q=1, we can obtain the kernel function (G) of the interaction of the unit tangential force with respect to any point. st (z1z2,z),G ss (z1z2,z),G sr (z1z2,z)).

[0135] In this embodiment, six kernel functions containing material anisotropy information were constructed using the unit load method, establishing a complete linear mapping relationship between the stress source at the hole edge and the stress response at any point in space. This can be reused under different load conditions, improving the computational efficiency of multi-condition analysis and realizing a systematic extension from single-hole solutions to multi-hole interaction problems.

[0136] In an exemplary embodiment, the specific implementation process of step 103, "processing the first stress data and the second stress data based on the force analysis algorithm to determine the interaction force data of the point to be measured," may include:

[0137] Based on the stress data at the edge of the hole and the second stress data of each hole, the external force data of each target arc segment of the hole is determined; based on the force balance algorithm, the external force data, the first stress data of each target arc segment of each hole and the shape parameters of each target arc, the interaction force data of the point to be measured is determined.

[0138] Among them, the external force data are the external load data acting on the target arc, including the equivalent force of the far-field stress generated at the hole edge and the known stress at the hole edge; the force balance algorithm is an algorithm determined based on the principle of static equilibrium. For example, for each target arc, the total internal force on it must be equal to the total external force.

[0139] Specifically, based on the stress data at the hole edge and the data of each secondary stress, the external force data of each target arc segment of the hole are determined. The external force data includes the equivalent force generated by the far-field stress at the hole edge and the combined effect of the known stresses at the hole edge. The calculation expression for the stress caused by the far-field stress at any point in a non-porous anisotropic medium can be:

[0140]

[0141] in, and These are the radial and tangential forces at the edge of the k-th hole, respectively, and the far-field stress ( , and They are the same size but opposite in direction.

[0142] For each target arc segment, establish a static equilibrium equation. The interaction force on that segment equals the external force acting on that segment plus the stress influence of all other target arcs on that segment through the kernel function. Combine the equilibrium equations for all target arcs of all holes to form a system of linear algebraic equations.

[0143] The system of linear algebraic equations comprises 2KM equations, where K is the number of holes and M is the number of discrete segments per hole. The unknowns in the system include the radial and tangential interaction forces on each target arc segment. Solving this system numerically yields all the interaction force data. The computational expression for the linear equations can be:

[0144]

[0145] in, M is the number of edge segments of the hole, K is the number of holes, and δ m Let be the arc length of each line segment of the m-th hole. and These are the radial and tangential interaction forces of the i-th arc segment of the k-th hole, respectively.

[0146] Setting P=1 and Q=0 respectively, we can obtain the kernel function (G) of the interaction of the unit radial force with respect to any point. rt (z1z2,z),G rs (z1z2,z) and G rr (z1z2,z)); Setting P=0 and Q=1, we can obtain the kernel function (G) of the interaction of the unit tangential force with respect to any point. st (z1z2,z), G ss (z1z2,z) and G sr (z1z2,z)), its expression is as follows:

[0147]

[0148]

[0149] Among them, G rt (z1z2,z), G rs (z1z2,z) and G rr (z1z2, z) represent the circumferential stress, shear stress, and radial stress at point z caused by a unit uniform radial stress (P) on a portion of the hole edge (z1Mz2), respectively. G st (z1z2,z), G ss (z1z2,z) and G sr (z1z2,z) represent the circumferential stress, shear stress, and radial stress at point z caused by a unit uniform tangential stress (Q) on a portion of the hole edge (z1Mz2), respectively.

[0150] In this embodiment, the problem of porous interaction is transformed into the direct solution of a system of linear algebraic equations by using a force balance algorithm, thereby achieving accurate quantification of the stress interference effect between porous structures.

[0151] In one exemplary embodiment, such as Figure 3 As shown, the specific implementation process of step 104, "Based on the stress analysis algorithm, analyze the remote stress field data, the first stress data of each target arc of each hole, and the data of each interaction force to obtain the target stress field data of the point to be measured," may include the following steps:

[0152] Step 301: Based on the remote stress field data and the position coordinate data of the point to be measured, determine the first stress field component of the point to be measured under the action of the remote stress.

[0153] Among them, the position coordinate data is the polar coordinate data of the point to be measured in the global coordinate system; the first stress field component refers to the stress generated at the point to be measured only by the far-end stress field.

[0154] Specifically, the terminal determines the first stress field component of the test point under the action of the far-end stress based on the far-end stress field data and the position coordinate data of the test point. In a pore-free, homogeneous, anisotropic medium, the first stress field component is equal to the value of the far-end stress at the test point.

[0155] Step 302: Based on the first stress data corresponding to each target arc of each hole and the interaction force data of each target arc, determine the second stress field component of the test point under the action of each target arc of each hole.

[0156] The second stress field component is the additional stress generated at the test point by all the holes and their interactions.

[0157] Specifically, based on the first stress data corresponding to each target arc in each hole and the interaction force data of each target arc, the terminal determines the second stress field component of the test point under the action of each target arc in each hole. Specifically, the interaction force of each target arc segment is multiplied by the corresponding kernel function, and then multiplied by the arc length of that segment to obtain the stress contribution generated by that segment at the test point. The total second stress field component is obtained by summing the stress contributions of all target arcs.

[0158] Step 303: Determine the target stress field data of the point to be measured based on the first stress field component and the second stress field component.

[0159] Specifically, based on the first and second stress field components, the calculation expression for determining the target stress field data of the point to be measured can be:

[0160]

[0161] in, , and These represent any points generated by the interaction of any number of loaded circular holes under far-field stress. Radial stress, circumferential stress, and tangential stress at point ) Let A represent the radial stress caused at any point A by a unit radial force along the j-th arc of the m-th hole. Let A represent the radial stress caused at any point A by a unit tangential force along the j-th arc of the m-th hole. Let A represent the circumferential stress caused at any point A by a unit radial traction force along the j-th arc of the m-th hole. Let A represent the circumferential stress caused at any point A by the unit tangential traction force along the j-th arc of the m-th hole. Let A represent the tangential stress caused at any point A by the unit radial traction force along the j-th arc of the m-th hole. Let A represent the tangential stress caused at any point A by the unit tangential traction force along the j-th arc of the m-th hole.

[0162] In one example, such as Figure 4 , 5 As shown, Figure 4 This is a schematic diagram of a single circular hole in an anisotropic rock stratum subjected to a uniformly distributed force along the hole edge, as shown below. Figure 5 As shown, Figure 5 It is the superposition of the various interaction forces at point A. It should be understood that this explanation is only for illustrative purposes and does not constitute a specific limitation.

[0163] In this embodiment, the stress field is modularly decomposed into far-field components and void disturbance components through the superposition principle, thereby realizing modular calculation of the stress field. This enables efficient acquisition of the calculation of any test point within the solution domain, improving the accuracy and efficiency of the analysis of the stress received at any test point.

[0164] In one embodiment, such as Figure 6 As shown, Figure 6 This embodiment analyzes the stress field of an anisotropic rock stratum containing two unequal-sized circular holes under far-field biaxial compressive stress and surface non-uniform stress. The analysis is based on the measured far-field stress of the surrounding rock. , and and the hole edge stress p n Establish such Figure 6 The model shown is an anisotropic elastic rock stratum containing two parallel unequal circular holes under far-field compressive stress and non-uniform surface compressive stress. The radii of the two holes are a1 = 0.5 m and a2 = 1 m, respectively. The far-field stress is... = =1MPa, =0 and hole edge stress p n=cos2θ1MPa, where t is the distance between the two circles (0.1-20m). The stress field analysis is performed using the method described in the above embodiment. Tables 1 and 2 show the stress field analysis of holes 1 and 2 at different distances t (0.1, 0.3, 0.5, 1, 2, 3, 5, 8, 12, 20m) and at different points on the two holes (…). =0°, 90°, 180°) under far-field stress (far-field stress) = =1MPa, =0) and surface uneven stress (the surface of the first circular hole is subjected to p) n A comparison of the analytical solution and the finite element method (Abaqus-based) solution for the tangential stress of an anisotropic rock stratum containing two unequal-sized circular holes (cos2θ = 1 MPa). The finite element method is used to simulate a 200m x 200m plate. It can be seen that the tangential stress on the two circular holes first decreases rapidly with increasing spacing t and then remains stable. Since the surface stress is greatest along the x-axis on the first circular hole, when... When the angle is 0°, the stress values ​​of the two circles are at their maximum. Furthermore, the error between the solution calculated in this embodiment and the numerical solution does not exceed 10%. It should be understood that this embodiment is for illustrative purposes only and does not constitute a specific limitation.

[0165] Table 1

[0166]

[0167] Table 2

[0168]

[0169] It should be understood that although the steps in the flowcharts of the embodiments described above are shown sequentially according to the arrows, these steps are not necessarily executed in the order indicated by the arrows. Unless explicitly stated herein, there is no strict order restriction on the execution of these steps, and they can be executed in other orders. Moreover, at least some steps in the flowcharts of the embodiments described above may include multiple steps or multiple stages. These steps or stages are not necessarily completed at the same time, but can be executed at different times. The execution order of these steps or stages is not necessarily sequential, but can be performed alternately or in turn with other steps or at least some of the steps or stages of other steps.

[0170] Based on the same inventive concept, this application also provides an apparatus for analyzing the pore stress field to implement the above-described method for analyzing the pore stress field. The solution provided by this apparatus is similar to the solution described in the above-described method. Therefore, the specific limitations of one or more embodiments of the apparatus for analyzing the pore stress field provided below can be found in the limitations of the pore stress field analysis method described above, and will not be repeated here.

[0171] In one exemplary embodiment, such as Figure 7 As shown, an analysis device 70 for a pore stress field is provided, comprising: a first determination module 71, a second determination module 72, a third determination module 73, and an analysis module 74, wherein:

[0172] The first determining module 71 is used to determine the plane mapping data corresponding to each hole in the target area based on the plane conformal mapping algorithm, anisotropic feature data and the shape parameters of the hole, and to determine the uniform stress distribution model of the hole based on the plane mapping data and the hole edge stress data.

[0173] The second determining module 72 is used to determine, based on the uniform stress distribution model, the first stress data of the target arc segment of the hole under unit load at the point to be measured; and to determine the second stress data of the point to be measured under remote stress based on the remote stress field data.

[0174] The third determining module 73 is used to process the first stress data and the second stress data based on the force analysis algorithm to determine the interaction force data of the point to be measured.

[0175] Analysis module 74 is used to analyze the remote stress field data, the first stress data and the interaction force data based on the stress analysis algorithm to obtain the target stress field data of the test point.

[0176] In one embodiment, the first determining module 71 is configured to determine the load data of the hole based on the hole edge stress data; calculate a first coefficient based on the load data and the anisotropic characteristic data; and determine a first potential function component of the hole based on the first coefficient and the plane mapping data.

[0177] Based on the load data, the mapping position data of each target arc segment, and the shape parameters of the hole, the load boundary data of each target arc segment of the hole is determined; based on the first coefficient and the load boundary data, the second potential function components of each target arc segment of the hole are determined.

[0178] Based on the first potential function component and the second potential function component, a uniform stress distribution model for the hole is determined.

[0179] In one embodiment, the second determining module 72 is used to analyze, for each of the target arcs, the uniform stress distribution model, the unit load data, and the position data of each target arc to obtain the uniform stress distribution data of the target arc.

[0180] Based on the anisotropic characteristic data and the uniform stress distribution data, the first stress data of the target arc under the test point under unit load is calculated.

[0181] In one embodiment, the third determining module 73 is used to determine the external force data of each target arc segment of the hole based on the hole edge stress data and each of the second stress data;

[0182] Based on the force balance algorithm, the external force data, the first stress data of each target arc of each hole, and the shape parameters of each target arc, the interaction force data of the point to be measured is determined.

[0183] In one embodiment, the analysis module 74 is used to determine the first stress field component of the test point under the action of the remote stress based on the remote stress field data and the position coordinate data of the test point.

[0184] Based on the first stress data corresponding to each target arc of each hole and the interaction force data of each target arc, the second stress field component of the test point under the action of each target arc of each hole is determined.

[0185] Based on the first stress field component and the second stress field component, the target stress field data of the point to be measured is determined.

[0186] In one embodiment, the device further includes:

[0187] The fourth determining module is used to determine the equivalent compliance coefficient based on the mechanical response data of the rock strata in the target area;

[0188] The equivalent orthogonal flexibility coefficient is calculated based on the equivalent flexibility coefficient, and the anisotropic direction characteristic parameters are determined based on the equivalent orthogonal flexibility coefficient.

[0189] Each module in the aforementioned pore stress field analysis device can be implemented entirely or partially through software, hardware, or a combination thereof. These modules can be embedded in the processor of a computer device in hardware form or independent of it, or stored in the memory of a computer device in software form, so that the processor can call and execute the operations corresponding to each module.

[0190] In one exemplary embodiment, a computer device is provided, which may be a terminal, and its internal structure diagram may be as follows: Figure 8 As shown, the computer device includes a processor, memory, input / output interface, communication interface, display unit, and input device. The processor, memory, and input / output interface are connected via a system bus, and the communication interface, display unit, and input device are also connected to the system bus via the input / output interface. The processor provides computational and control capabilities. The memory includes non-volatile storage media and internal memory. The non-volatile storage media stores the operating system and computer programs. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage media. The input / output interface is used for exchanging information between the processor and external devices. The communication interface is used for wired or wireless communication with external terminals; wireless communication can be achieved through Wi-Fi, mobile cellular networks, Near Field Communication (NFC), or other technologies. When executed by the processor, the computer program implements a method for analyzing the stress field of a hole. The display unit is used to form a visually visible image and can be a display screen, projection device, or virtual reality imaging device. The display screen can be an LCD screen or an e-ink screen. The input device of the computer device can be a touch layer covering the display screen, or buttons, trackballs, or touchpads set on the casing of the computer device, or external keyboards, touchpads, or mice, etc.

[0191] Those skilled in the art will understand that Figure 8 The structure shown is merely a block diagram of a portion of the structure related to the present application and does not constitute a limitation on the computer device to which the present application is applied. Specific computer devices may include more or fewer components than those shown in the figure, or combine certain components, or have different component arrangements.

[0192] In one embodiment, a computer device is provided, including a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to implement the steps in the above-described method embodiments.

[0193] In one embodiment, a computer-readable storage medium is provided having a computer program stored thereon that, when executed by a processor, implements the steps in the above method embodiments.

[0194] In one embodiment, a computer program product is provided, including a computer program that, when executed by a processor, implements the steps in the above method embodiments.

[0195] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, data stored, data displayed, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of the relevant data must comply with relevant regulations.

[0196] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium, and when executed, it can include the processes of the embodiments of the above methods. Any references to memory, databases, or other media used in the embodiments provided in this application can include at least one of non-volatile memory and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM can take many forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM). The databases involved in the embodiments provided in this application may include at least one type of relational database and non-relational database. Non-relational databases may include, but are not limited to, blockchain-based distributed databases. The processors involved in the embodiments provided in this application may be general-purpose processors, central processing units, graphics processing units, digital signal processors, programmable logic devices, quantum computing-based data processing logic devices, artificial intelligence (AI) processors, etc., and are not limited to these.

[0197] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this application.

[0198] The embodiments described above are merely illustrative of several implementation methods of this application, and while the descriptions are specific and detailed, they should not be construed as limiting the scope of this patent application. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of this application, and these all fall within the protection scope of this application. Therefore, the protection scope of this application should be determined by the appended claims.

Claims

1. A method for analyzing the stress field of a cavity, characterized in that, The method includes: For each hole in the target area, based on the planar conformal mapping algorithm, anisotropic feature data and the shape parameters of the hole, the planar mapping data corresponding to the hole is determined, and based on the planar mapping data and the hole edge stress data, the uniform stress distribution model of the hole is determined. Based on the uniform stress distribution model, the first stress data of the test point under each segment of the target arc in the hole under a unit load is determined; based on the far-end stress field data, the second stress data of the test point under the far-end stress is determined. Based on the force analysis algorithm, the first stress data and the second stress data are processed to determine the interaction force data of the point to be measured; Based on the stress analysis algorithm, the remote stress field data, the first stress data, and the interaction force data are analyzed to obtain the target stress field data of the test point.

2. The method according to claim 1, characterized in that, The step of determining the uniform stress distribution model of the hole based on the planar mapping data and the hole edge stress data includes: Based on the stress data at the edge of the hole, the load data of the hole is determined; based on the load data and the anisotropic characteristic data, a first coefficient is calculated; and based on the first coefficient and the plane mapping data, a first potential function component of the hole is determined. Based on the load data, the mapping position data of each target arc segment, and the shape parameters of the hole, the load boundary data of each target arc segment of the hole is determined; based on the first coefficient and the load boundary data, the second potential function components of each target arc segment of the hole are determined. Based on the first potential function component and the second potential function component, a uniform stress distribution model for the hole is determined.

3. The method according to claim 1, characterized in that, Based on the uniform stress distribution model, the determination of the first stress data of the target arc segment of the hole under unit load at the test point includes: For each target arc, the uniform stress distribution data of the target arc is obtained by analyzing the unit load data and the position data of each target arc based on the uniform stress distribution model. Based on the anisotropic characteristic data and the uniform stress distribution data, the first stress data of the target arc under the test point under unit load is calculated.

4. The method according to claim 1, characterized in that, The force analysis algorithm is used to process the first stress data and the second stress data to determine the interaction force data of the point to be measured, including: Based on the stress data at the edge of the hole and the second stress data, the external force data of each target arc segment of the hole is determined; Based on the force balance algorithm, the external force data, the first stress data of each target arc of each hole, and the shape parameters of each target arc, the interaction force data of the point to be measured is determined.

5. The method according to claim 1, characterized in that, The stress analysis algorithm is used to analyze the remote stress field data, the first stress data, and the interaction force data to obtain the target stress field data of the test point, including: Based on the remote stress field data and the position coordinate data of the point to be measured, the first stress field component of the point to be measured under the action of the remote stress is determined. Based on the first stress data corresponding to each target arc of each hole and the interaction force data of each target arc, the second stress field component of the test point under the action of each target arc of each hole is determined. Based on the first stress field component and the second stress field component, the target stress field data of the point to be measured is determined.

6. The method according to claim 1, characterized in that, The method further includes: Based on the mechanical response data of the rock strata in the target area, the equivalent compliance coefficient is determined; The equivalent orthogonal flexibility coefficient is calculated based on the equivalent flexibility coefficient, and the anisotropic direction characteristic parameters are determined based on the equivalent orthogonal flexibility coefficient.

7. An analytical device for analyzing the stress field of a cavity, characterized in that, The device includes: The first determining module is used to determine the plane mapping data corresponding to each hole in the target area based on the plane conformal mapping algorithm, anisotropic feature data and the shape parameters of the hole, and to determine the uniform stress distribution model of the hole based on the plane mapping data and the hole edge stress data. The second determining module is used to determine, based on the uniform stress distribution model, the first stress data of the target arc segment of the hole under unit load at the point to be measured; and to determine the second stress data of the point to be measured under remote stress based on the remote stress field data. The third determining module is used to process the first stress data and the second stress data based on the force analysis algorithm to determine the interaction force data of the point to be measured. The analysis module is used to analyze the remote stress field data, the first stress data, and the interaction force data based on the stress analysis algorithm to obtain the target stress field data of the test point.

8. A computer device comprising a memory and a processor, wherein the memory stores a computer program, characterized in that, When the processor executes the computer program, it implements the steps of the method according to any one of claims 1 to 6.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it implements the steps of the method according to any one of claims 1 to 6.

10. A computer program product, comprising a computer program, characterized in that, When the computer program is executed by a processor, it implements the steps of the method according to any one of claims 1 to 6.

Images