CN-121977949-A - Hot water power multi-field monitoring analysis method and system suitable for TEHB test system

CN121977949ACN 121977949 ACN121977949 ACN 121977949ACN-121977949-A

Abstract

The invention provides a thermal power multi-field monitoring and analyzing method and a system suitable for TEHB test systems, wherein the thermal power multi-field monitoring and analyzing method comprises a model frame building step, a steady initial field solving step, a thermal power coupling solving step and a multi-field response data acquisition and recording step, and also provides a thermal power multi-field monitoring and analyzing system suitable for TEHB test systems, wherein a three-dimensional model comprising a six-way waveguide rod and a rock sample is built, a mechanical, seepage and temperature steady initial field before impact loading is obtained, thermal power coupling dynamic calculation is carried out under the three-axis six-way synchronous impact action, and the curves and statistical indexes of key results such as stress, crack, pore pressure, temperature and the like are automatically output. The invention can realize impact process simulation under the condition of hot water coupling loading, and monitor and analyze multi-field response data, thereby being used for revealing the dynamic mechanical behavior of the rock under the action of hot water coupling and providing basis for scientific design of deep engineering.

Inventors

ZHU JIANBO
HE QI
WANG SILI
LIANG LINSHENG
DING WENDUO
Zhang Kongwei
ZHOU XINYANG
WANG KANGYU
XIE HEPING
BAO WEIYUE
ZHOU TAO
LIAO ZHIYI
ZHANG SHIWEI
Cen Zhuo
XIE CHENGCHENG

Assignees

深圳大学

Dates

Publication Date: 20260505
Application Date: 20260407

Claims (10)

1. The hydrothermal force multi-field monitoring and analyzing method suitable for TEHB test system is characterized by comprising the following steps: The method comprises the steps of establishing a three-dimensional model comprising six-directional waveguide rods and rock samples, wherein the waveguide rods are waveguide rods in a triaxial six-directional dynamic true triaxial electromagnetic Hopkinson rod TEHB test system, realizing propagation and loading of stress waves based on a finite element method, generating a discrete element rock sample model by the rock samples based on a discrete element method, and setting physical parameters of the waveguide rods, rock mechanical parameters, rock thermal parameters and fluid physical parameters as a hot water force coupling region to form a coupling model frame for multi-field calculation; A steady-state initial field solving step, namely respectively carrying out mechanical, seepage and thermal field steady-state initialization calculation on the coupling model frame to form a steady-state initial field with mechanical balance, stable pore pressure and flow and constant temperature before impact loading; A thermal hydraulic coupling solving step, namely constructing an impact loading waveform and applying the impact loading waveform to the loading end surface of the six-way waveguide rod, and synchronously propelling the mechanical response, the coupling solving of a temperature field and a seepage field until the set impact duration is reached, so as to obtain a multi-field response result of the whole impact process; and the multi-field response data acquisition and recording step is to extract, sort and output the test results of the stress field, the fracture field, the pore pressure field and the temperature field of the rock sample, so as to realize the observable, quantifiable and comparable analysis of the multi-field response of the rock test under the loading of the hydrothermal coupling.
2. The method for monitoring and analyzing the thermal force multiple fields of the TEHB test system according to claim 1, wherein the model frame establishing step includes a model frame preliminary establishing step A1 and a mechanical parameter setting step A2, wherein the mechanical parameter setting step A2 includes the following sub-steps: setting material parameters for the waveguide rod in the three-dimensional model, wherein the material parameters comprise density, elastic modulus and poisson ratio and are used for describing the mechanical response characteristic of the waveguide rod in the three-dimensional model when the waveguide rod is used as a continuous medium; setting rock mechanical parameters for a discrete element rock sample model, wherein the rock mechanical parameters are used for representing basic mechanical behaviors of the rock sample under the condition of true triaxial stress; Setting temperature-related parameters for the rock sample, wherein the temperature-related parameters comprise heat conductivity, specific heat capacity and thermal expansion coefficient and are used for describing the influence of a temperature field on the mechanical behavior of the rock test; Setting seepage related parameters for the rock sample, wherein the seepage related parameters comprise porosity and permeability coefficient and are used for describing the flow characteristics of a seepage field in the rock sample and the coupling effect between the seepage field and the mechanical behavior.
3. The method for monitoring and analyzing the hydrothermal force of a TEHB test system according to claim 2, wherein the rock sample is cubic, the rock sample is discretely divided by adopting complete triangulation, six-way waveguide rods are respectively arranged on the outer sides of 6 outer surfaces of the rock sample, the waveguide rods adopt elastic continuous body models, the rock sample adopts elastic models, the contact adopts a mole-coulomb model and sets contact rigidity, cohesive force, tensile strength and friction angle, and meanwhile, fluid parameters, rock heat conductivity, specific heat and thermal expansion coefficient and fluid-rock heat exchange coefficient are configured, and then a coupling model frame capable of being used for calculation of the hydrothermal force coupling is formed through a continuous body-discrete body coupling interface.
4. The method for monitoring and analyzing the hydrothermal force of TEHB in a test system according to claim 1, wherein the step of solving the steady-state initial field includes a step of solving the mechanical steady-state initial field, a step of solving the seepage steady-state initial field, and a step of solving the temperature steady-state initial field, B3, The processing method of the mechanical steady-state initial field solving step B1 comprises the steps of applying displacement constraint to the six-way waveguide rod in the non-loading direction, enabling the six-way waveguide rod to only allow displacement along the corresponding impact loading direction, applying a gravity field to the three-dimensional model, and obtaining a mechanical balance state before impact loading through iterative solution; Applying constant temperature boundary conditions to six boundary surfaces of a rock sample, performing heat conduction evolution on the rock sample in a preset constant temperature environment, and performing thermal field calculation by loop iteration until a temperature field reaches a stable state, so as to obtain a temperature steady-state initial field before impact loading; Setting an initial pore pressure field of a seepage medium in a rock sample, respectively applying pore pressure boundary conditions at a water inlet and a water outlet to form a seepage field driven by a pressure gradient, and then starting a seepage module to calculate until the seepage field is stable, so as to obtain an initial pore pressure field and a flow field before impact loading.
5. The method for monitoring and analyzing the hydrothermal force of TEHB in the test system according to claim 4, wherein in the step B2 of solving the seepage steady-state initial field, the calculation formula of the pore pressure field is as follows: Wherein, the For the hole pressure in the last calculation step, For the bulk modulus of the fluid, In order for the net inflow to be a net, In order to calculate the step size, In order to calculate the volume of the cell, In order to provide a volume change amount, As volume average value, when pore pressure When the pore pressure field tends to be stable, the pore pressure field is indicated to be stable, so that an initial pore pressure field is obtained; In the temperature steady-state initial field solving step B3, the calculation formula of the temperature field is as follows: Wherein, the As the temperature at the time t is the temperature, For the moment of time Is used for the temperature control of the liquid crystal display device, In order to calculate the step of time, In order to increase the temperature of the product, Is a heat-conducting item, and is a heat-conducting item, As an additional heat source item, the heat source, For the temperature to be updated by the coefficient, Is the total heat capacity of the area, To calculate the step size, when the formula value Approaching 0, it is shown that the temperature has reached a steady state, thus obtaining a temperature steady state initial field.
6. The method for monitoring and analyzing the thermal power multi-field suitable for TEHB test systems according to claim 1, wherein the thermal power coupling solving step includes an energy-absorbing boundary condition setting step C1, an impact waveform constructing and synchronous impact loading step C2, a heat source coupling updating step C3 and a full coupling propulsion solving step C4, The energy absorption boundary condition setting step C1 is used for setting an absorption boundary to inhibit reflected waves; The step C2 of constructing and synchronously loading the shock wave is used for constructing the shock wave and loading the shock wave to the loading end face of the waveguide rod, so that the stress wave in the waveguide rod is transferred to the rock sample, and the dynamic shock loading under the true triaxial condition is realized; the heat source coupling updating step C3 is used for realizing a coupling link of mechanical dissipation, heat source input and temperature rise, and the reaction coupling of a temperature field to mechanical response and crack evolution; And the full-coupling propulsion solving step C4 is used for starting seepage quick solving and thermal field implicit solving to realize multi-field coupling solving in the impact short-time strong nonlinear process, and then propelling the thermal-hydraulic full-coupling analysis to a set duration, and outputting stress wave response, temperature rise evolution, seepage field response and crack evolution test results in the impact whole process.
7. The method for monitoring and analyzing the hydrothermal force of TEHB in a test system according to claim 6, wherein the heat source coupling updating step C3 is performed by: C301, counting contact friction energy and calculating increment, namely, carrying out energy statistics on contact behaviors inside a rock sample in the power stepping process, obtaining contact friction energy and calculating friction energy increment of adjacent time steps The calculation formula is as follows: In the formula, Representing the initial time period of time that is required, Representing the time of the break of the contact, Representing frictional energy generated at the contact; For the contact friction energy at the initial time, The contact friction energy is the contact breaking time; And C302, equivalently using friction energy increment as a heat source and performing space distribution, namely using the friction energy increment released by contact damage as a local heat source, distributing the local heat source to adjacent area blocks according to a space neighborhood rule, and driving a temperature field to update, thereby realizing a coupling link of mechanical dissipation, heat source input and temperature rise, wherein a specific calculation formula is as follows: Wherein, the Representing the heat transfer efficiency of the heat transfer device, Representing the total heat capacity of the zone, Representing the initial temperature of the product, Is the updated temperature; And C303, temperature updating and thermal expansion feedback, namely calculating thermal strain and thermal expansion effect according to the updated temperature field, and feeding back the thermal strain and the thermal expansion effect to a mechanical solving process to update the deformation and contact state of the block, so as to realize reaction coupling of the temperature field to mechanical response and crack evolution, wherein the calculation formula of the thermal strain is as follows: wherein: In terms of the coefficient of linear thermal expansion, As the value of the change in temperature, For the kronecker function, 1 when i=j, or 0, i, j=1, 2,3 respectively represent the spatial coordinate direction, Is the increment of thermal strain.
8. The method for monitoring and analyzing the hydrothermal force of TEHB in a test system according to claim 1, wherein the multi-field response data acquisition and recording step includes a stress field data processing step D1, a fracture monitoring and statistics step D2, a pore pressure data processing step D3 and a temperature data processing step D4, The stress field data processing step D1 is used for outputting stress and strain time-course data and automatically generating a stress-strain curve and a characteristic value based on the stress and strain time-course data; the fracture monitoring and statistics step D2 is used for monitoring fracture states, carrying out statistics on fracture indexes, and outputting curves or summarized results of the fracture indexes changing along with time; the pore pressure data processing step D3 is used for recording the data of the change of the pore pressure of the representative position in the rock along with time, processing the data and outputting the data; And the temperature data processing step D4 is used for recording the time-varying data of the temperature of the sample center and the representative position, processing the data and outputting the data.
9. The method for monitoring and analyzing the hydrothermal force of TEHB in a test system according to claim 8, wherein in the step D2 of monitoring and counting the cracks, the crack index includes the crack opening, the crack type and the number of cracks, The judging and processing formula of the crack type is as follows: Wherein, the In order to increase the amount of shear stress, In order to increase the amount of shear displacement, In order for the shear stress to be a high shear stress, In order to achieve a tangential contact stiffness, In order to contact the cohesive force of the adhesive, In order to contact the friction angle of the friction, For the ultimate shear strength, the shear strength, As an increment of the normal stress, As a result of the normal stress, For the normal contact stiffness, In order to contact the normal displacement increment, For ultimate tensile strength, when exposed to shear forces Reaching the ultimate shear strength After that, marked as shear crack, when contacted with normal force Exceeding the ultimate tensile strength And then marked as tension cracks.
10. A hydrothermal force multi-field monitoring system adapted for use in TEHB test systems for implementing the hydrothermal force multi-field monitoring analysis method of any one of claims 1-9 adapted for use in TEHB test systems, comprising: The model frame building module is used for building a three-dimensional model comprising six-direction waveguide rods and rock samples, wherein the waveguide rods are waveguide rods in a triaxial six-direction dynamic true triaxial electromagnetic Hopkinson rod TEHB test system, stress wave propagation and loading are realized based on a finite element method, the rock samples are used for generating discrete element rock sample models based on a discrete element method and are used as hot water force coupling areas, and waveguide rod physical parameters, rock mechanical parameters, rock thermal parameters and fluid physical parameters are set to form a coupling model frame for multi-field calculation; The steady-state initial field solving module is used for respectively carrying out steady-state initialization calculation on the mechanical, seepage and thermal fields of the coupling model frame to form a steady-state initial field with stable mechanical balance, pore pressure and flow and constant temperature before impact loading; The thermal hydraulic coupling solving module is used for constructing an impact loading waveform and applying the impact loading waveform to the loading end surface of the six-way waveguide rod, and synchronously propelling the mechanical response, the coupling solving of the temperature field and the seepage field until the set impact duration is reached, so as to obtain a multi-field response result of the whole impact process; And the multi-field response data acquisition recording module is used for extracting, sorting and outputting test results of stress fields, fracture fields, pore pressure fields and temperature fields of the rock test sample, and realizing observable, quantifiable and comparable analysis of multi-field response of the rock test under the condition of hydrothermal coupling loading.

Description

Hot water power multi-field monitoring analysis method and system suitable for TEHB test system Technical Field The invention relates to the technical field of rock impact test simulation, in particular to a hydrothermal force multi-field monitoring and analyzing method and system suitable for TEHB test systems. Background The environment where the deep rock mass is located has typical characteristics of high ground stress, high ground temperature (heat), high osmotic pressure (water) and the like, and frequently bears transient impact load in the process of engineering exploitation and construction activities, thus forming typical 'three-high-one disturbance' characteristics. Aiming at the research of rock dynamics characteristics in the complex environment, the method has important engineering application value and scientific significance. The dynamic true triaxial electromagnetic Hopkinson bar (True Triaxial Electromagnetic Hopkinson Bar, TEHB) test system can realize multiaxial and multidirectional synchronous impact loading, and can truly simulate the stress condition of a deep rock mass under the condition of dynamic disturbance. However, when a dynamic true triaxial six-directional impact test is performed under a hot water coupling condition, the test is usually performed in a closed constraint environment in order to ensure stability of confining pressure, pore pressure and temperature control boundary conditions and realize synchronous control of a loading process. The method is limited by a sealing structure form and test space arrangement conditions, and the multi-axis stress field in the sample and the coupling relation between the multi-axis stress field and the pore pressure field, the temperature field and the crack evolution process are difficult to realize visual and continuous monitoring. The information of multiple fields is often obtained by means of indirect measurement, so that deep knowledge of microscopic crack evolution paths, stress wave propagation characteristics and action mechanisms between the microscopic crack evolution paths, stress wave propagation characteristics and the relation between the stress wave propagation characteristics and seepage processes and temperature feedback is restricted to a certain extent. In addition, limited by the sensor layout mode and sampling conditions, the observed quantity which can be obtained through experiments is limited, and multi-field response data are difficult to obtain continuously under the conditions of larger time scale, space scale and richer working conditions. In order to further reveal the dynamic characteristics of the rock under the thermal hydraulic multi-field coupling, and expand the observable and comparable data dimension in the time scale, the space scale and the working condition range, it is necessary to construct a set of reliable numerical simulation methods and analysis systems for important extension and effective supplement of TEHB tests. Disclosure of Invention In order to solve the problem that multi-field response data are difficult to obtain continuously under larger time scale, space scale and richer working conditions in the prior art, the invention provides a hydrothermal force multi-field monitoring and analyzing method suitable for a TEHB test system, and also provides a hydrothermal force multi-field monitoring and analyzing system for realizing the hydrothermal force multi-field monitoring and analyzing method suitable for a TEHB test system, and realizing observable, quantifiable and comparable analysis of multi-field response of a rock sample under the condition of hydrothermal force coupling loading. The invention is suitable for a hydrothermal force multi-field monitoring analysis method of TEHB test systems, which comprises the following steps: The method comprises the steps of establishing a three-dimensional model comprising six-directional waveguide rods and rock samples, wherein the waveguide rods are waveguide rods in a triaxial six-directional dynamic true triaxial electromagnetic Hopkinson rod TEHB test system, realizing propagation and loading of stress waves based on a finite element method, generating a discrete element rock sample model by the rock samples based on a discrete element method, and setting physical parameters of the waveguide rods, rock mechanical parameters, rock thermal parameters and fluid physical parameters as a hot water force coupling region to form a coupling model frame for multi-field calculation; A steady-state initial field solving step, namely respectively carrying out mechanical, seepage and thermal field steady-state initialization calculation on the coupling model frame to form a steady-state initial field with mechanical balance, stable pore pressure and flow and constant temperature before impact loading; A thermal hydraulic coupling solving step, namely constructing an impact loading waveform and applying the impact loading waveform to the loading end