CN110895276A - Expansion evolution considered method and device for simulating hard gypsum rock tunnel - Google Patents
Expansion evolution considered method and device for simulating hard gypsum rock tunnel Download PDFInfo
- Publication number
- CN110895276A CN110895276A CN201910723608.XA CN201910723608A CN110895276A CN 110895276 A CN110895276 A CN 110895276A CN 201910723608 A CN201910723608 A CN 201910723608A CN 110895276 A CN110895276 A CN 110895276A
- Authority
- CN
- China
- Prior art keywords
- expansion
- water
- temperature
- water absorption
- establishing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000011435 rock Substances 0.000 title claims abstract description 95
- 238000000034 method Methods 0.000 title claims abstract description 73
- 229910052602 gypsum Inorganic materials 0.000 title claims abstract description 30
- 239000010440 gypsum Substances 0.000 title claims abstract description 30
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 253
- 238000004364 calculation method Methods 0.000 claims abstract description 55
- 238000010521 absorption reaction Methods 0.000 claims abstract description 54
- 230000006870 function Effects 0.000 claims abstract description 47
- 238000004088 simulation Methods 0.000 claims abstract description 46
- 238000002076 thermal analysis method Methods 0.000 claims abstract description 43
- 230000008569 process Effects 0.000 claims abstract description 36
- 238000012360 testing method Methods 0.000 claims abstract description 31
- 229910052925 anhydrite Inorganic materials 0.000 claims abstract description 30
- OSGAYBCDTDRGGQ-UHFFFAOYSA-L calcium sulfate Chemical compound [Ca+2].[O-]S([O-])(=O)=O OSGAYBCDTDRGGQ-UHFFFAOYSA-L 0.000 claims abstract description 30
- 230000009471 action Effects 0.000 claims abstract description 28
- 239000000463 material Substances 0.000 claims abstract description 25
- 230000014509 gene expression Effects 0.000 claims abstract description 15
- 230000004907 flux Effects 0.000 claims abstract description 8
- 230000008859 change Effects 0.000 claims abstract description 7
- 230000001052 transient effect Effects 0.000 claims description 23
- 238000004458 analytical method Methods 0.000 claims description 18
- 238000009412 basement excavation Methods 0.000 claims description 11
- 238000006243 chemical reaction Methods 0.000 claims description 11
- 230000000694 effects Effects 0.000 claims description 9
- 238000010586 diagram Methods 0.000 description 10
- 238000003860 storage Methods 0.000 description 8
- 238000004590 computer program Methods 0.000 description 7
- 238000012545 processing Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 230000008961 swelling Effects 0.000 description 3
- 230000036962 time dependent Effects 0.000 description 3
- 238000009826 distribution Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000012916 structural analysis Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- XLYOFNOQVPJJNP-ZSJDYOACSA-N heavy water Substances [2H]O[2H] XLYOFNOQVPJJNP-ZSJDYOACSA-N 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/082—Investigating permeability by forcing a fluid through a sample
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/16—Investigating or analyzing materials by the use of thermal means by investigating thermal coefficient of expansion
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/20—Investigating or analyzing materials by the use of thermal means by investigating the development of heat, i.e. calorimetry, e.g. by measuring specific heat, by measuring thermal conductivity
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Pathology (AREA)
- Immunology (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Medicinal Chemistry (AREA)
- Food Science & Technology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Dispersion Chemistry (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a method and a device for simulating a hard gypsum tunnel by considering expansion evolution, wherein the method comprises the following steps: carrying out an expansion test under the action of pressure water to obtain a first derivative function of water absorption and water pressure; introducing the first derivative function into a seepage differential equation, and establishing a water absorption-expansion evolution model in the seepage field; establishing a corresponding relation between a thermal analysis model and a water absorption-expansion evolution model in a seepage field, and converting specific values of a thermal conductivity coefficient, a boundary temperature or a heat flux in the thermal analysis model, and expressions of specific heat capacity and a thermal expansion coefficient by using specific water capacity, relative density of rock relative to water, expansion stress, elastic modulus of a material, expansion modulus of the anhydrite, an expansion coefficient in the thermal analysis model and an actual expansion coefficient of the anhydrite; and establishing a general calculation flow, simulating the process of the deformation of the surrounding rock and the change of the supporting stress along with time under the action of seepage and water absorption-expansion evolution by using the general calculation flow, and outputting a simulation result.
Description
Technical Field
The invention relates to the technical field of tunnel simulation, in particular to a method and a device for simulating a hard gypsum rock tunnel by considering expansion evolution.
Background
The seepage movement has obvious influence on the stability and deformation characteristics of the expansive rock tunnel. Compared with other expansive rocks, the seepage process has more obvious influence on the hard gypsum tunnel and is reflected in two typical characteristics of the hard gypsum tunnel: the swelling time is long (including the seepage time and the water absorption time) and the swelling pressure is not obviously stable (seepage causes the water absorption to change along with the time). In engineering, the expansion force of a rock test piece is usually obtained through laboratory tests and is used as a fixed input parameter for simulating the deformation characteristic of an expansion rock tunnel.
However, in the existing engineering, the expansive force of a rock test piece is usually obtained through indoor tests and is used as an input parameter for simulating the deformation characteristic of an expansive rock tunnel, so that the influence of seepage on the expansive rock cannot be reflected, and the expansion evolution and spatial distribution characteristics of the tunnel anhydrite surrounding rock cannot be reflected.
Disclosure of Invention
The present invention aims to provide a method and apparatus for simulating a hard gypsum tunnel taking into account dilatational evolution which overcomes one of the problems mentioned above or at least partially solves any of the problems mentioned above.
In order to achieve the purpose, the technical scheme of the invention is realized as follows:
one aspect of the present invention provides a method for simulating a hard gypsum tunnel considering expansive evolution, including: carrying out an expansion test under the action of pressure water to obtain a first derivative function of water absorption and water pressure; introducing the first derivative function into a seepage differential equation, and establishing a water absorption-expansion evolution model in the seepage field; establishing a corresponding relation between a thermal analysis model and a water absorption-expansion evolution model in a seepage field, and converting specific values of a thermal conductivity coefficient, a boundary temperature or a heat flux in the thermal analysis model, and expressions of specific heat capacity and a thermal expansion coefficient by using specific water capacity, relative density of rock relative to water, expansion stress, elastic modulus of a material, expansion modulus of the anhydrite, an expansion coefficient in the thermal analysis model and an actual expansion coefficient of the anhydrite; and establishing a general calculation flow, simulating the process of the deformation of the surrounding rock and the change of the supporting stress along with time under the action of seepage and water absorption-expansion evolution by using the general calculation flow, and outputting a simulation result.
Wherein, the expansion test under the action of pressure water is carried out, and the first derivative function of the water absorption and the water pressure is obtained by the following steps: carrying out expansion tests on the anhydrite under different water pressure conditions, and fitting the relation between the water pressure and the water absorption rate; determining the relation between the pressure water head and the water content according to the relation between the water pressure and the water absorption; and obtaining a first derivative function of the water absorption and the water pressure according to the relation between the pressure water head and the water content.
The method for building the water absorption-expansion evolution model in the seepage field by using the first derivative function as a reference in the seepage differential equation comprises the following steps: determining the value of the specific water capacity of the rock according to a first derivative function of the water absorption and the water pressure; and establishing a water absorption-expansion evolution model in the seepage field according to the value of the specific water capacity of the rock.
Wherein, the general calculation process includes: establishing a geometric model division network; carrying out parameter assigning and boundary condition and heat load applying on the geometric model; carrying out heat conduction steady state calculation to obtain an initial state; carrying out heat conduction transient calculation to obtain an instantaneous state; extracting the temperature of the expansion surrounding rock unit; updating the specific heat capacity and the thermal expansion coefficient of the material; storing the thermal expansion coefficient and the initial temperature; judging whether the ending time is reached, if not, increasing the time and returning to execute the transient calculation of heat conduction to obtain the transient state; if yes, saving the calculation result and entering a structure analysis process, wherein the structure analysis process comprises the following steps: the method comprises the steps of simulating ground stress, simulating excavation and support, extracting a thermal expansion coefficient to endow a unit, converting temperature into load to be applied to surrounding rock, and simulating the expansion effect of the surrounding rock.
Wherein, the specific heat capacity and the thermal expansion coefficient of the updating material comprise: setting the thermal expansion coefficient to be zero when the temperature of a certain unit is lower than the temperature of the previous step; when the temperature of one unit is higher than the maximum expansion temperature, setting a thermal expansion coefficient according to the ratio of the maximum expansion temperature to the temperature of the other unit; and when the temperature of a certain unit is greater than or equal to the temperature of the previous step and less than the maximum expansion temperature, setting the thermal expansion coefficient by utilizing the conversion relation between the specific water capacity and the specific heat capacity.
In another aspect, the present invention provides a device for simulating a hard gypsum rock tunnel considering expansion evolution, including: the calculation module is used for performing an expansion test under the action of pressure water to obtain a first derivative function of water absorption and water pressure; the establishing module is used for introducing the first derivative function into a seepage differential equation and establishing a water absorption-expansion evolution model in the seepage field; the conversion module is used for establishing a corresponding relation between the thermal analysis model and a water absorption-expansion evolution model in the seepage field, and converting specific values of heat conductivity coefficient, boundary temperature or heat flux in the thermal analysis model, and specific heat capacity and thermal expansion coefficient expressions by using specific water capacity, relative density of rock relative to water, expansion stress, elastic modulus of material and expansion modulus of the hard gypsum as well as expansion coefficient in the thermal analysis model and actual expansion coefficient of the hard gypsum; and the simulation module is used for establishing a general calculation flow, simulating the process of time-varying deformation and supporting stress of the surrounding rock under the action of seepage and water absorption-expansion evolution by using the general calculation flow, and outputting a simulation result.
The calculation module performs an expansion test under the action of pressure water in the following way to obtain a first derivative function of water absorption and water pressure: the calculation module is specifically used for performing expansion tests on the anhydrite under different water pressure conditions and fitting the relation between the water pressure and the water absorption rate; determining the relation between the pressure water head and the water content according to the relation between the water pressure and the water absorption; and obtaining a first derivative function of the water absorption and the water pressure according to the relation between the pressure water head and the water content.
The establishing module is used for establishing a water absorption-expansion evolution model in the seepage field by introducing a first derivative function into a seepage differential equation in the following way: the establishing module is specifically used for determining the value of the specific water capacity of the rock according to a first derivative function of the water absorption and the water pressure; and establishing a water absorption-expansion evolution model in the seepage field according to the value of the specific water capacity of the rock.
The simulation module establishes a general calculation flow in the following way: the simulation module is specifically used for establishing a geometric model division network; carrying out parameter assigning and boundary condition and heat load applying on the geometric model; carrying out heat conduction steady state calculation to obtain an initial state; carrying out heat conduction transient calculation to obtain an instantaneous state; extracting the temperature of the expansion surrounding rock unit; updating the specific heat capacity and the thermal expansion coefficient of the material; storing the thermal expansion coefficient and the initial temperature; judging whether the ending time is reached, if not, increasing the time and returning to execute the transient calculation of heat conduction to obtain the transient state; if yes, saving the calculation result and entering a structure analysis process, wherein the structure analysis process comprises the following steps: the method comprises the steps of simulating ground stress, simulating excavation and support, extracting a thermal expansion coefficient to endow a unit, converting temperature into load to be applied to surrounding rock, and simulating the expansion effect of the surrounding rock.
Wherein, the simulation module updates the specific heat capacity and the thermal expansion coefficient of the material by the following method: the simulation module is specifically used for setting the thermal expansion coefficient to be zero when the temperature of a certain unit is lower than the temperature of the previous step; when the temperature of one unit is higher than the maximum expansion temperature, setting a thermal expansion coefficient according to the ratio of the maximum expansion temperature to the temperature of the other unit; and when the temperature of a certain unit is greater than or equal to the temperature of the previous step and less than the maximum expansion temperature, setting the thermal expansion coefficient by utilizing the conversion relation between the specific water capacity and the specific heat capacity.
Therefore, the method and the device for simulating the hard gypsum rock tunnel considering the expansion evolution provided by the embodiment of the invention are based on the expansion test under the action of pressure water, and obtain the first derivative function of the water absorption rate to the water pressure from the aspect of phenomenology. The derivative function is introduced into a seepage differential equation, and a water absorption-expansion evolution model in the seepage field is established. Based on the consistency of the model and the thermal analysis model on a mathematical expression, a general calculation flow is established by means of an ANSYS thermal analysis model, and the process of time-dependent changes of surrounding rock deformation and supporting stress under the action of seepage and water absorption-expansion evolution is simulated. By adopting a water absorption-expansion evolution model in the seepage field, the purpose of simulating the water absorption-expansion evolution of the anhydrite rock can be realized, namely the dynamic changes of the seepage of the effluent, the rock water absorption and the rock expansion along with the time can be simulated, so that the simulation result is closer to the real situation; by adopting the parameter corresponding relation between the water absorption-expansion evolution model in the seepage field and the ANSYS thermal analysis model, a simulation result can be simply and conveniently obtained and is used for analyzing the influence of expansion evolution on the tunnel.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the description of the embodiments are briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on the drawings without creative efforts.
Fig. 1 is a flowchart of a method for simulating a hard gypsum tunnel considering expansive evolution according to an embodiment of the present invention;
fig. 2 is a flowchart of a concrete example of a simulation method for a hard gypsum tunnel considering expansion evolution according to an embodiment of the present invention;
fig. 3 is a schematic structural diagram of the apparatus for simulating a hard gypsum tunnel considering expansion evolution according to the embodiment of the present invention.
Detailed Description
Exemplary embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.
Fig. 1 shows a flowchart of a method for simulating a hard gypsum tunnel considering dilatational evolution according to an embodiment of the present invention, and referring to fig. 1, the method for simulating a hard gypsum tunnel considering dilatational evolution according to an embodiment of the present invention includes:
and S1, performing an expansion test under the action of pressure water to obtain a first derivative function of water absorption and water pressure.
And S2, introducing the first derivative function into a seepage differential equation, and establishing a water absorption-expansion evolution model in the seepage field.
Specifically, as an alternative implementation of the embodiment of the present invention, the performing of the expansion test under the action of the pressurized water to obtain the first derivative function of the water absorption and the water pressure includes: carrying out expansion tests on the anhydrite under different water pressure conditions, and fitting the relation between the water pressure and the water absorption rate; determining the relation between the pressure water head and the water content according to the relation between the water pressure and the water absorption; and obtaining a first derivative function of the water absorption and the water pressure according to the relation between the pressure water head and the water content.
As an optional implementation manner of the embodiment of the present invention, the building of the water absorption-expansion evolution model in the seepage field by referring the first derivative function to the seepage differential equation includes: determining the value of the specific water capacity of the rock according to a first derivative function of the water absorption and the water pressure; and establishing a water absorption-expansion evolution model in the seepage field according to the value of the specific water capacity of the rock.
In practical application, firstly, the expansion test of the anhydrite rock under the action of pressure water is carried out, and the method specifically comprises the following steps:
expansion tests of the anhydrite rock were carried out under different hydraulic pressure conditions. The water absorption rate A in the expansion process of the anhydrite rock is obtained through teststAnd fitting the relation between the water pressure and the water absorption according to the formula (1):
wherein Delta A is a parameter reflecting the water absorption capacity, Ap0Is the maximum water absorption of the rock under a water pressure of 0MPa, ap0Is the water absorption coefficient under the water pressure of 0MPa, delta a is a parameter reflecting the change of the water absorption coefficient along with the water pressure, p is the water pressure, and t is the time.
The water pressure p in formula (1) can be replaced by a pressure head H, AtAvailable mass water content wmInstead, after conversion, formula (1) becomes formula (2):
the above equation can be derived for the variable H:
definition of C in the present methodmFor rock specific water capacity, it means the mass (kg) of water that needs to be absorbed per 1m head added to a rock mass (kg), which is expressed in m-1. From the definition of mass water cut, C can be determinedmThe values of (A) are as follows:
c obtained by experimentmThe method can reflect the change relation of the rock water absorption along with the water pressure and the time, and lays a foundation for further simulating the water absorption-expansion evolution of the anhydrite.
S3, establishing a corresponding relation between the thermal analysis model and the water absorption-expansion evolution model in the seepage field, and converting specific values of heat conductivity coefficient, boundary temperature or heat flux in the thermal analysis model, and specific heat capacity and thermal expansion coefficient expressions by using specific water capacity, relative density of rock relative to water, expansion stress, elastic modulus of material and expansion modulus of hard gypsum, and expansion coefficient in the thermal analysis model and actual expansion coefficient of hard gypsum.
Specifically, as an optional implementation manner of the embodiment of the present invention, the general computing process includes: establishing a geometric model division network; carrying out parameter assigning and boundary condition and heat load applying on the geometric model; carrying out heat conduction steady state calculation to obtain an initial state; carrying out heat conduction transient calculation to obtain an instantaneous state; extracting the temperature of the expansion surrounding rock unit; updating the specific heat capacity and the thermal expansion coefficient of the material; storing the thermal expansion coefficient and the initial temperature; judging whether the ending time is reached, if not, increasing the time and returning to execute the transient calculation of heat conduction to obtain the transient state; if yes, saving the calculation result and entering a structure analysis process, wherein the structure analysis process comprises the following steps: the method comprises the steps of simulating ground stress, simulating excavation and support, extracting a thermal expansion coefficient to endow a unit, converting temperature into load to be applied to surrounding rock, and simulating the expansion effect of the surrounding rock.
As an optional implementation manner of the embodiment of the present invention, the updating of the specific heat capacity and the thermal expansion coefficient of the material includes: setting the thermal expansion coefficient to be zero when the temperature of a certain unit is lower than the temperature of the previous step; when the temperature of one unit is higher than the maximum expansion temperature, setting a thermal expansion coefficient according to the ratio of the maximum expansion temperature to the temperature of the other unit; and when the temperature of a certain unit is greater than or equal to the temperature of the previous step and less than the maximum expansion temperature, setting the thermal expansion coefficient by utilizing the conversion relation between the specific water capacity and the specific heat capacity.
In practical application, a parameter value method in an ANSYS thermal analysis model needs to be established, the value method can have a corresponding relation with a water absorption-expansion evolution model in a seepage field, and parameters in the thermal analysis model are converted by using the corresponding relation.
The invention establishes the corresponding relation between a thermal analysis system and a water absorption-expansion evolution model in a seepage field, and the corresponding relation is shown in table 1. It should be noted that: the subscript i in the table indicates the direction of propagation; x, y and z respectively correspond to three directions in a space coordinate system; cmIs specific water, determined by the results of the swelling test in the above step; rhoreIs the relative density of rock relative to water; sigmasIs an expansion stress; e and Esα, the modulus of elasticity of the material and the modulus of expansion of the anhydrite rock, respectivelyTAnd α are the actual coefficients of expansion in the thermal analysis model and in the anhydrite rock, respectively;
TABLE 1 corresponding relationship between thermal analysis model and water absorption-expansion evolution model in seepage field
The parameter value taking method comprises the following steps:
the temperature is in equal proportion to the water head, namely 1 ℃ corresponds to the water head of 1 m; the absorbed heat and the water absorption mass are in equal proportion, namely 1J corresponds to 1kg of water; the mass and density of the heat transfer material are made to be the same as those of the anhydrite; other key parameters take the following values:
λ＝k (5)
CT＝Cm/1000 (6)
αT＝(Es/E)αCmΔH/ΔT＝(Es/E)αCm(7)
and S4, establishing a general calculation process, simulating the process of time-varying surrounding rock deformation and supporting stress under the action of seepage and water absorption-expansion evolution by using the general calculation process, and outputting a simulation result.
The water absorption-expansion evolution characteristic of the hard gypsum in the seepage field is simulated by means of an ANSYS thermal analysis model, and the influence of the characteristic on the tunnel engineering can be analyzed in combination with the tunnel excavation process. The general simulation flow is shown in fig. 2, and the specific steps are discussed as follows:
4.1 establishing a geometric model and dividing meshes. The geometric model comprises surrounding rocks, a supporting structure and a tunnel outline; the surrounding rock considering seepage and water absorption-expansion characteristics adopts a thermal analysis unit; the supporting structure can comprise anchor rods, primary lining concrete and secondary lining concrete.
4.2 parameters are assigned, boundary conditions are applied and thermal loads are applied. The invention establishes the corresponding relation between the thermal analysis system and the water absorption-expansion equation system in the seepage field, and the corresponding relation is shown in table 1 and table 1. From this correspondence, specific values of the thermal conductivity, the temperature or the heat flux at the boundary, and the specific heat capacity (expression (3)), the thermal expansion coefficient expression are converted. The concrete values of the parameters and the mechanical parameters of the surrounding rock and the supporting structure are directly used as the input of the model, and the expressions of specific heat capacity and thermal expansion coefficient are reserved for the next analysis. It should be noted that: the thermal expansion coefficient of the surrounding rock without taking the water-swelling property into consideration is zero.
4.3 heat conduction steady state and transient calculations. Transient heat transfer analysis is used because the analysis process involves dynamic changes in seepage fields and water absorption-expansion. It should be noted that, when the initial state of the seepage field is unknown or the position head needs to be considered, a small time step, such as 0.001, may be set first, and the time integrator is turned off, that is, the initial state of the transient analysis is obtained by the steady-state analysis method, and then the time integrator is turned on for the transient analysis.
4.4 extracting the temperature of the expansion surrounding rock unit. In transient thermal analysis, a small time step (which can be set to 1 day) needs to be set, and the temperature (corresponding to the water pressure) of the expansion surrounding rock unit in the calculation result of each step is extracted and reserved for the next step of calculation.
4.5 updating the specific heat capacity and thermal expansion coefficient of the material. The unit temperature and time are substituted into formula (3), and the specific heat capacity is calculated and given to the corresponding unit. This step updates the thermal expansion coefficient of the material, considering three cases, which are as follows:
4.5.1 when the temperature of a certain unit is lower than the temperature of the previous step, the thermal expansion coefficient is zero;
4.5.2 when the CaSO4 in the interior of the anhydrite is completely converted into CaSO 4.2H2O, the anhydrite reaches the theoretical maximum expansion amount, and the cell reaching the maximum expansion amount is defined as the maximum expansion cell. The temperature corresponding to the water head can be determined according to the water absorption equation, and the temperature is defined as the maximum expansion temperature delta Tmax. When a certain cell temperature is higher than the maximum expansion temperature, the thermal expansion coefficient is given as the following equation (8).
Wherein the subscript j is the cell number and indicates the variable corresponding to the jth cell
4.5.3 for cells with a temperature greater than or equal to the temperature of the previous step and less than the maximum expansion temperature, an assignment is made according to the parameter conversion relationship mentioned earlier, i.e., equation (3).
4.6 store the coefficient of thermal expansion and the initial temperature. And storing the thermal expansion coefficient of each expansion surrounding rock unit at each time step and the temperature of each unit at the initial state, and reserving the thermal expansion coefficient and the temperature for structural analysis.
4.7 simulated ground stress. And when the preset simulation time is reached, the structural analysis is carried out. The ground stress is simulated by applying a downward acceleration to the unit.
4.8 simulation excavation and support. And simulating the excavation and supporting process of the tunnel by using an ANSYS unit life and death technology. In the actual process, the ground stress after excavation is released in a time-based manner, the node stress of the node of the excavation surface can be extracted, the node stress in the opposite direction is applied according to engineering experience, and the process of slowly releasing the stress of the surrounding rock in the actual engineering is simulated. In addition, a creep material can be added to the surrounding rock, and the time-varying rule of the engineering structure body, namely the releasing process of the stress of the surrounding rock, is reflected through the creep characteristic of the material.
4.9 simulation of the expansion effect of the surrounding rock. Firstly, the initial temperature of each expansion surrounding rock unit stored in 3.6 is taken as a reference temperature and is respectively given to each unit; secondly, extracting the thermal expansion coefficient stored in 3.6 at the corresponding moment and endowing the thermal expansion coefficient to a corresponding unit for obtaining the stress and displacement distribution of the surrounding rock at a certain moment; finally, the LDEAD command in ANSYS is adopted to convert the temperature into the load applied to the unit.
Therefore, the method for simulating the hard gypsum rock tunnel considering the expansion evolution provided by the embodiment of the invention is based on the expansion test under the action of pressure water, and obtains the first derivative function of the water absorption rate to the water pressure from the aspect of phenomenology. The derivative function is introduced into a seepage differential equation, and a water absorption-expansion evolution model in the seepage field is established. Based on the consistency of the model and the thermal analysis model on a mathematical expression, a general calculation flow is established by means of an ANSYS thermal analysis model, and the process of time-dependent changes of surrounding rock deformation and supporting stress under the action of seepage and water absorption-expansion evolution is simulated. By adopting a water absorption-expansion evolution model in the seepage field, the purpose of simulating the water absorption-expansion evolution of the anhydrite rock can be realized, namely the dynamic changes of the seepage of the effluent, the rock water absorption and the rock expansion along with the time can be simulated, so that the simulation result is closer to the real situation; by adopting the parameter corresponding relation between the water absorption-expansion evolution model in the seepage field and the ANSYS thermal analysis model, a simulation result can be simply and conveniently obtained and is used for analyzing the influence of expansion evolution on the tunnel.
Therefore, the invention has the following effects:
1. the specific water capacity is obtained through an expansion test under the action of pressure water and is used as a key input parameter of a water absorption-expansion evolution model in the seepage field, so that the purpose of simulating expansion evolution is realized, and a simulation result is closer to a real situation.
2. Based on the consistency of the new model and the thermal analysis simulation on the mathematical expression, the aim of quickly and simply obtaining the simulation result is fulfilled by means of an ANSYS thermal analysis model.
Fig. 3 shows a schematic structural diagram of the expansion evolution considered anhydrite tunnel simulation device according to the embodiment of the present invention, where the expansion evolution considered anhydrite tunnel simulation device is applied to the expansion evolution considered anhydrite tunnel simulation method, and only a brief description is given below of a structure of the expansion evolution considered anhydrite tunnel simulation device, but other things are not considered to the utmost, please refer to the relevant description of the expansion evolution considered anhydrite tunnel simulation method, and are not repeated here. Referring to fig. 3, an embodiment of the present invention provides a device for simulating a hard gypsum rock tunnel considering expansion evolution, including:
the calculation module is used for performing an expansion test under the action of pressure water to obtain a first derivative function of water absorption and water pressure;
the establishing module is used for introducing the first derivative function into a seepage differential equation and establishing a water absorption-expansion evolution model in the seepage field;
the conversion module is used for establishing a corresponding relation between the thermal analysis model and a water absorption-expansion evolution model in the seepage field, and converting specific values of heat conductivity coefficient, boundary temperature or heat flux in the thermal analysis model, and specific heat capacity and thermal expansion coefficient expressions by using specific water capacity, relative density of rock relative to water, expansion stress, elastic modulus of material and expansion modulus of the hard gypsum as well as expansion coefficient in the thermal analysis model and actual expansion coefficient of the hard gypsum;
and the simulation module is used for establishing a general calculation flow, simulating the process of time-varying deformation and supporting stress of the surrounding rock under the action of seepage and water absorption-expansion evolution by using the general calculation flow, and outputting a simulation result.
As an optional implementation manner of the embodiment of the present invention, the calculation module performs an expansion test under the action of pressurized water in the following manner to obtain a first derivative function of water absorption and water pressure: the calculation module is specifically used for performing expansion tests on the anhydrite under different water pressure conditions and fitting the relation between the water pressure and the water absorption rate; determining the relation between the pressure water head and the water content according to the relation between the water pressure and the water absorption; and obtaining a first derivative function of the water absorption and the water pressure according to the relation between the pressure water head and the water content.
As an optional implementation manner of the embodiment of the present invention, the establishing module refers the first derivative function to the seepage differential equation in the following manner, and establishes the water absorption-expansion evolution model in the seepage field: the establishing module is specifically used for determining the value of the specific water capacity of the rock according to a first derivative function of the water absorption and the water pressure; and establishing a water absorption-expansion evolution model in the seepage field according to the value of the specific water capacity of the rock.
As an optional implementation of the embodiment of the present invention, the simulation module establishes a general computation flow by: the simulation module is specifically used for establishing a geometric model division network; carrying out parameter assigning and boundary condition and heat load applying on the geometric model; carrying out heat conduction steady state calculation to obtain an initial state; carrying out heat conduction transient calculation to obtain an instantaneous state; extracting the temperature of the expansion surrounding rock unit; updating the specific heat capacity and the thermal expansion coefficient of the material; storing the thermal expansion coefficient and the initial temperature; judging whether the ending time is reached, if not, increasing the time and returning to execute the transient calculation of heat conduction to obtain the transient state; if yes, saving the calculation result and entering a structure analysis process, wherein the structure analysis process comprises the following steps: the method comprises the steps of simulating ground stress, simulating excavation and support, extracting a thermal expansion coefficient to endow a unit, converting temperature into load to be applied to surrounding rock, and simulating the expansion effect of the surrounding rock.
As an optional implementation of the embodiment of the present invention, the simulation module updates the specific heat capacity and the thermal expansion coefficient of the material by: the simulation module is specifically used for setting the thermal expansion coefficient to be zero when the temperature of a certain unit is lower than the temperature of the previous step; when the temperature of one unit is higher than the maximum expansion temperature, setting a thermal expansion coefficient according to the ratio of the maximum expansion temperature to the temperature of the other unit; and when the temperature of a certain unit is greater than or equal to the temperature of the previous step and less than the maximum expansion temperature, setting the thermal expansion coefficient by utilizing the conversion relation between the specific water capacity and the specific heat capacity.
Therefore, the device for simulating the hard gypsum rock tunnel considering expansion evolution provided by the embodiment of the invention obtains the first derivative function of water absorption rate to water pressure from the aspect of phenomenology based on the expansion test under the action of pressure water. The derivative function is introduced into a seepage differential equation, and a water absorption-expansion evolution model in the seepage field is established. Based on the consistency of the model and the thermal analysis model on a mathematical expression, a general calculation flow is established by means of an ANSYS thermal analysis model, and the process of time-dependent changes of surrounding rock deformation and supporting stress under the action of seepage and water absorption-expansion evolution is simulated. By adopting a water absorption-expansion evolution model in the seepage field, the purpose of simulating the water absorption-expansion evolution of the anhydrite rock can be realized, namely the dynamic changes of the seepage of the effluent, the rock water absorption and the rock expansion along with the time can be simulated, so that the simulation result is closer to the real situation; by adopting the parameter corresponding relation between the water absorption-expansion evolution model in the seepage field and the ANSYS thermal analysis model, a simulation result can be simply and conveniently obtained and is used for analyzing the influence of expansion evolution on the tunnel.
Therefore, the invention has the following effects:
1. the specific water capacity is obtained through an expansion test under the action of pressure water and is used as a key input parameter of a water absorption-expansion evolution model in the seepage field, so that the purpose of simulating expansion evolution is realized, and a simulation result is closer to a real situation.
2. Based on the consistency of the new model and the thermal analysis simulation on the mathematical expression, the aim of quickly and simply obtaining the simulation result is fulfilled by means of an ANSYS thermal analysis model.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
In a typical configuration, a computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.
The memory may include forms of volatile memory in a computer readable medium, Random Access Memory (RAM) and/or non-volatile memory, such as Read Only Memory (ROM) or flash memory (flash RAM). The memory is an example of a computer-readable medium.
Computer-readable media, including both non-transitory and non-transitory, removable and non-removable media, may implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), Static Random Access Memory (SRAM), Dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), Read Only Memory (ROM), Electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), Digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information that can be accessed by a computing device. As defined herein, a computer readable medium does not include a transitory computer readable medium such as a modulated data signal and a carrier wave.
The above are merely examples of the present application and are not intended to limit the present application. Various modifications and changes may occur to those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the scope of the claims of the present application.
Claims (10)
1. A method for simulating a hard gypsum tunnel by considering expansion evolution is characterized by comprising the following steps:
carrying out an expansion test under the action of pressure water to obtain a first derivative function of water absorption and water pressure;
the first derivative function is referred to a seepage differential equation, and a water absorption-expansion evolution model in a seepage field is established;
establishing a corresponding relation between a thermal analysis model and a water absorption-expansion evolution model in the seepage field, and converting specific values of a thermal conductivity coefficient, a boundary temperature or a heat flux in the thermal analysis model, and expressions of specific heat capacity and a thermal expansion coefficient by using specific water capacity, relative density of rock relative to water, expansion stress, elastic modulus of a material, expansion modulus of the hard gypsum, expansion coefficient in the thermal analysis model and actual expansion coefficient of the hard gypsum;
and establishing a general calculation flow, simulating the process of the deformation of the surrounding rock and the change of the supporting stress along with time under the action of seepage and water absorption-expansion evolution by utilizing the general calculation flow, and outputting a simulation result.
2. The method of claim 1, wherein said performing a test of expansion under pressurized water to obtain a first derivative of water absorption and water pressure comprises:
carrying out expansion tests on the anhydrite under different water pressure conditions, and fitting the relation between the water pressure and the water absorption rate;
determining the relation between a pressure water head and the water content according to the relation between the water pressure and the water absorption;
and obtaining a first derivative function of the water absorption and the water pressure according to the relation between the pressure water head and the water content.
3. The method of claim 1 or 2, wherein the first derivative function is referred to a seepage differential equation, and establishing a water absorption-expansion evolution model in a seepage field comprises:
determining the value of the specific water capacity of the rock according to the first derivative function of the water absorption and the water pressure;
and establishing a water absorption-expansion evolution model in the seepage field according to the value of the specific water capacity of the rock.
4. The method of claim 1, wherein the general computing process comprises:
establishing a geometric model division network;
carrying out parameter assigning, boundary condition applying and heat load on the geometric model;
carrying out heat conduction steady state calculation to obtain an initial state;
carrying out heat conduction transient calculation to obtain an instantaneous state;
extracting the temperature of the expansion surrounding rock unit;
updating the specific heat capacity and the thermal expansion coefficient of the material;
storing the thermal expansion coefficient and the initial temperature;
judging whether the ending time is reached, if not, increasing the time and returning to execute the transient calculation of heat conduction to obtain the transient state; if so, saving the calculation result and entering a structure analysis process, wherein the structure analysis process comprises the following steps: the method comprises the steps of simulating ground stress, simulating excavation and support, extracting a thermal expansion coefficient to endow a unit, converting temperature into load to be applied to surrounding rock, and simulating the expansion effect of the surrounding rock.
5. The method of claim 4, wherein the updating the specific heat capacity and coefficient of thermal expansion of the material comprises:
setting the thermal expansion coefficient to zero when the temperature of a certain unit is lower than the temperature of the previous step;
when the temperature of one unit is higher than the maximum expansion temperature, setting the thermal expansion coefficient according to the ratio of the maximum expansion temperature to the temperature of the other unit;
and when the temperature of a certain unit is greater than or equal to the temperature of the previous step and less than the maximum expansion temperature, setting the thermal expansion coefficient by utilizing the conversion relation between the specific water capacity and the specific heat capacity.
6. A anhydrite tunnel simulation device considering expansion evolution is characterized by comprising:
the calculation module is used for performing an expansion test under the action of pressure water to obtain a first derivative function of water absorption and water pressure;
the establishing module is used for introducing the first derivative function into a seepage differential equation and establishing a water absorption-expansion evolution model in a seepage field;
the conversion module is used for establishing a corresponding relation between a thermal analysis model and a water absorption-expansion evolution model in the seepage field, and converting specific values of heat conductivity coefficient, boundary temperature or heat flux in the thermal analysis model, and specific heat capacity and thermal expansion coefficient expressions by using specific water capacity, relative density of rock relative to water, expansion stress, elastic modulus of material and expansion modulus of the gypsum as well as expansion coefficient in the thermal analysis model and actual expansion coefficient of the gypsum;
and the simulation module is used for establishing a general calculation flow, simulating the process of time-varying surrounding rock deformation and supporting stress under the action of seepage and water absorption-expansion evolution by utilizing the general calculation flow, and outputting a simulation result.
7. The device of claim 6, wherein the calculation module performs the expansion test under the action of the pressurized water by obtaining a first derivative function of the water absorption and the water pressure as follows:
the calculation module is specifically used for carrying out expansion tests on the anhydrite under different water pressure conditions and fitting the relation between the water pressure and the water absorption rate; determining the relation between a pressure water head and the water content according to the relation between the water pressure and the water absorption; and obtaining a first derivative function of the water absorption and the water pressure according to the relation between the pressure water head and the water content.
8. The apparatus of claim 6 or 7, wherein the establishing module establishes the water uptake-expansion evolution model in the seepage field by referencing the first derivative function to a seepage differential equation as follows:
the establishing module is specifically used for determining a value of the specific water capacity of the rock according to the first derivative function of the water absorption and the water pressure; and establishing a water absorption-expansion evolution model in the seepage field according to the value of the specific water capacity of the rock.
9. The apparatus of claim 6, wherein the simulation module establishes the general computational flow by:
the simulation module is specifically used for establishing a geometric model division network; carrying out parameter assigning, boundary condition applying and heat load on the geometric model; carrying out heat conduction steady state calculation to obtain an initial state; carrying out heat conduction transient calculation to obtain an instantaneous state; extracting the temperature of the expansion surrounding rock unit; updating the specific heat capacity and the thermal expansion coefficient of the material; storing the thermal expansion coefficient and the initial temperature; judging whether the ending time is reached, if not, increasing the time and returning to execute the transient calculation of heat conduction to obtain the transient state; if so, saving the calculation result and entering a structure analysis process, wherein the structure analysis process comprises the following steps: the method comprises the steps of simulating ground stress, simulating excavation and support, extracting a thermal expansion coefficient to endow a unit, converting temperature into load to be applied to surrounding rock, and simulating the expansion effect of the surrounding rock.
10. The apparatus of claim 9, wherein the simulation module updates the specific heat capacity and coefficient of thermal expansion of the material by:
the simulation module is specifically used for setting the thermal expansion coefficient to be zero when the temperature of a certain unit is lower than the temperature of the previous step; when the temperature of one unit is higher than the maximum expansion temperature, setting the thermal expansion coefficient according to the ratio of the maximum expansion temperature to the temperature of the other unit; and when the temperature of a certain unit is greater than or equal to the temperature of the previous step and less than the maximum expansion temperature, setting the thermal expansion coefficient by utilizing the conversion relation between the specific water capacity and the specific heat capacity.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910723608.XA CN110895276B (en) | 2019-08-07 | 2019-08-07 | Expansion evolution considered method and device for simulating hard gypsum rock tunnel |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910723608.XA CN110895276B (en) | 2019-08-07 | 2019-08-07 | Expansion evolution considered method and device for simulating hard gypsum rock tunnel |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110895276A true CN110895276A (en) | 2020-03-20 |
| CN110895276B CN110895276B (en) | 2022-03-08 |
Family
ID=69785748
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910723608.XA Active CN110895276B (en) | 2019-08-07 | 2019-08-07 | Expansion evolution considered method and device for simulating hard gypsum rock tunnel |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110895276B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111458273A (en) * | 2020-04-16 | 2020-07-28 | 四川大学 | Model and method for discriminating hydraulic pressure lifting of porous medium under chemical reaction |
| CN113092720A (en) * | 2021-04-02 | 2021-07-09 | 交通运输部公路科学研究所 | Rock lateral confinement expansion constitutive relation analysis method |
| CN113515803A (en) * | 2021-09-14 | 2021-10-19 | 中国矿业大学(北京) | Method and system for evaluating heat transfer and storage of surrounding rocks of subway under action of periodic surface temperature and wind temperature |
| CN114739883A (en) * | 2022-03-31 | 2022-07-12 | 中铁十一局集团第五工程有限公司 | Generalized analysis method for rock expansion |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104634729A (en) * | 2015-02-04 | 2015-05-20 | 中电建路桥集团有限公司 | Multi-factor gyprock expansion testing method |
| CN106202980A (en) * | 2016-08-24 | 2016-12-07 | 山西省交通科学研究院 | A kind of swelled ground is humidified under Condition of Rainfall Infiltration and expands method for numerical simulation |
| CN206161649U (en) * | 2016-09-21 | 2017-05-10 | 中国矿业大学(北京) | Tunnelling simulation modeling experiment device |
| CN106951612A (en) * | 2017-03-06 | 2017-07-14 | 河海大学 | Dynamic water storage capacity Runoff calculation method in freeze-thawing process of soil |
| CN107729645A (en) * | 2017-10-11 | 2018-02-23 | 重庆大学 | A kind of decision method of Swelling Rock tunnel inverted arch supporting stability |
| CN108038305A (en) * | 2017-12-08 | 2018-05-15 | 重庆大学 | A kind of anhydrite tunnel reliability analysis method for considering expansion and developing |
| CN108344675A (en) * | 2018-02-08 | 2018-07-31 | 四川大学 | Coal body adopts the test method of permeation fluid mechanics rule under the conditions of simulation protective coat extracted |
| CN108959758A (en) * | 2018-06-27 | 2018-12-07 | 大连海大恒信土木工程技术有限公司 | A kind of tunnel Local security analysis method considering seepage flow creep |
| CN109522588A (en) * | 2018-09-27 | 2019-03-26 | 天津大学 | A kind of foundation and calculation method of water-heat of oil-well cement ring-mechanics multi- scenarios method model |
| CN110018290A (en) * | 2019-04-16 | 2019-07-16 | 中国矿业大学(北京) | One kind being based on two-dimensional analog testing stand overlying strata crack evolution Fractal Methods |
-
2019
- 2019-08-07 CN CN201910723608.XA patent/CN110895276B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104634729A (en) * | 2015-02-04 | 2015-05-20 | 中电建路桥集团有限公司 | Multi-factor gyprock expansion testing method |
| CN106202980A (en) * | 2016-08-24 | 2016-12-07 | 山西省交通科学研究院 | A kind of swelled ground is humidified under Condition of Rainfall Infiltration and expands method for numerical simulation |
| CN206161649U (en) * | 2016-09-21 | 2017-05-10 | 中国矿业大学(北京) | Tunnelling simulation modeling experiment device |
| CN106951612A (en) * | 2017-03-06 | 2017-07-14 | 河海大学 | Dynamic water storage capacity Runoff calculation method in freeze-thawing process of soil |
| CN107729645A (en) * | 2017-10-11 | 2018-02-23 | 重庆大学 | A kind of decision method of Swelling Rock tunnel inverted arch supporting stability |
| CN108038305A (en) * | 2017-12-08 | 2018-05-15 | 重庆大学 | A kind of anhydrite tunnel reliability analysis method for considering expansion and developing |
| CN108344675A (en) * | 2018-02-08 | 2018-07-31 | 四川大学 | Coal body adopts the test method of permeation fluid mechanics rule under the conditions of simulation protective coat extracted |
| CN108959758A (en) * | 2018-06-27 | 2018-12-07 | 大连海大恒信土木工程技术有限公司 | A kind of tunnel Local security analysis method considering seepage flow creep |
| CN109522588A (en) * | 2018-09-27 | 2019-03-26 | 天津大学 | A kind of foundation and calculation method of water-heat of oil-well cement ring-mechanics multi- scenarios method model |
| CN110018290A (en) * | 2019-04-16 | 2019-07-16 | 中国矿业大学(北京) | One kind being based on two-dimensional analog testing stand overlying strata crack evolution Fractal Methods |
Non-Patent Citations (5)
| Title |
|---|
| 任松 等: "压力水作用下硬石膏岩膨胀性研究", 《岩土力学》 * |
| 吴顺川 等: "基于膨胀本构的石膏岩隧道", 《岩土力学》 * |
| 李强 等: "硬石膏岩膨胀特性试验及隧道抗膨胀衬砌设计", 《地下空间与工程学报》 * |
| 王辅圣: "渗水条件下膨胀土隧道围岩湿度场", 《铁道勘察》 * |
| 陈钒 等: "基于湿度应力场理论的硬石膏岩膨胀试验研究", 《岩土力学》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111458273A (en) * | 2020-04-16 | 2020-07-28 | 四川大学 | Model and method for discriminating hydraulic pressure lifting of porous medium under chemical reaction |
| CN113092720A (en) * | 2021-04-02 | 2021-07-09 | 交通运输部公路科学研究所 | Rock lateral confinement expansion constitutive relation analysis method |
| CN113515803A (en) * | 2021-09-14 | 2021-10-19 | 中国矿业大学(北京) | Method and system for evaluating heat transfer and storage of surrounding rocks of subway under action of periodic surface temperature and wind temperature |
| CN113515803B (en) * | 2021-09-14 | 2021-12-03 | 中国矿业大学(北京) | Method and system for evaluating heat transfer and storage of surrounding rocks of subway under action of periodic surface temperature and wind temperature |
| CN114739883A (en) * | 2022-03-31 | 2022-07-12 | 中铁十一局集团第五工程有限公司 | Generalized analysis method for rock expansion |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110895276B (en) | 2022-03-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110895276B (en) | Expansion evolution considered method and device for simulating hard gypsum rock tunnel | |
| Peng et al. | LOQUAT: an open-source GPU-accelerated SPH solver for geotechnical modeling | |
| Lisjak et al. | Acceleration of a 2D/3D finite-discrete element code for geomechanical simulations using General Purpose GPU computing | |
| CN109992900B (en) | Multi-field real-time online collaborative intelligent simulation method and system for mass concrete | |
| CN108038305B (en) | Anhydrite tunnel reliability analysis method considering expansion evolution | |
| Alves et al. | An object-oriented approach to the generalized finite element method | |
| CN111125963B (en) | Numerical simulation system and method based on Lagrange integral point finite element | |
| Tian et al. | A challenging dam structural analysis: large-scale implicit thermo-mechanical coupled contact simulation on Tianhe-II | |
| Khoshghalb et al. | Does the upper bound solution property of the Node-based Smoothed Point Interpolation Methods (NSPIMs) hold true in coupled flow-deformation problems of porous media? | |
| Ceyhan et al. | Towards a faster and improved ADCIRC (ADvanced Multi-Dimensional CIRCulation) model | |
| Rathnayaka et al. | A three-dimensional (3-D) meshfree-based computational model to investigate stress-strain-time relationships of plant cells during drying | |
| Zheng et al. | A mathematical programming approach for frictional contact problems with the extended finite element method | |
| Pons et al. | Adaptive mesh refinement method applied to shallow water model: a mass conservative projection | |
| CN103279685B (en) | Based on the Reservoir behavior analogy method of reverse automatic differential | |
| Li et al. | Memory access optimization of molecular dynamics simulation software crystal-md on sunway taihulight | |
| Witteveen et al. | On the projection of a flexible bodies modal coordinates onto another finite element model with local modifications | |
| CN113051790B (en) | Steam load loading method, system, equipment and medium for finite element simulation | |
| Skalko et al. | Application of a fundamental solution to a problem operator for modeling vibrations of elastically supported load-bearing structures | |
| Desai et al. | Damage model and implementation in nonlinear dynamic problems | |
| Dong et al. | Assessment of applicability of the material point method in offshore geotechnical engineering | |
| Si et al. | Development of simulation program for temperature field of mass concrete structures | |
| Shafee et al. | Augmenting low‐order finite element method with partial nodal strain smoothing for flow‐deformation analysis of geomechanical problems | |
| Labib | Utilizing High Performance Computing to Improve the Application of Machine Learning for Time-Efficient Prediction of Buildings’ Daylighting Performance | |
| Iizuka et al. | Parallel simulation system for earthquake generation: fault analysis modules and parallel coupling analysis | |
| Lim et al. | Using adaptive control for dynamic substructuring tests |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |