Rock mass is a fractured porous medium usually subjected to complex geostress and fluid pressure simultaneously.Moreover,the properties of rock mass change in time and space due to mining-induced fractures.Therefore,i...Rock mass is a fractured porous medium usually subjected to complex geostress and fluid pressure simultaneously.Moreover,the properties of rock mass change in time and space due to mining-induced fractures.Therefore,it is always challenging to accurately measure rock mass properties.In this study,a three-dimensional(3D)microseismic(MS)data-driven damage model for jointed rock mass under hydro-mechanical coupling conditions is proposed.It is a 3D finite element model that takes seepage,damage and stress field effects into account jointly.Multiple factors(i.e.joints,water and microseismicity)are used to optimize the rock mass mechanical parameters at different scales.The model is applied in Shirengou iron mine to study the damage evolution of rock mass and assess the crown pillar stability during the transition from open-pit to underground mining.It is found that the damage pattern is mostly controlled by the structure,water and rock mass parameters.The damage pattern is evidently different from the two-dimensional result and is more consistent with the field observations.This difference is caused by the MS-derived damage acting on the rock mass.MS data are responsible for gradually correcting the damage zone,changing the direction in which it expands,and promoting it to evolve close to reality.For the crown pillar,the proposed model yields a more trustworthy safety factor.In order to guarantee the stability of the pillar,it is suggested to take waterproof and reinforcement measures in areas with a high degree of damage.展开更多
The joint of clay core-wall and concrete cut-off wall is one of the weakest parts in high earth and rockftll dams.A kind of highly plastic clay is always fixed on the joint to fit the large shear deformation between c...The joint of clay core-wall and concrete cut-off wall is one of the weakest parts in high earth and rockftll dams.A kind of highly plastic clay is always fixed on the joint to fit the large shear deformation between clay core-wall and concrete cut-offwall,so the hydro-mechanical coupling mechanisms on the joint under high stress,high hydraulic gradient,and large shear deformation are of great importance for the evaluation of dam safety.The hydro-mechanical coupling characteristics of the joint of the highly plastic clay and the concrete cut-off wall in a high earth and rockfill dam in China were studied by using a newly designed soil-structure contact erosion apparatus.The experimental results indicate that：1） Shear failure on the joint is due to the hydro-mechanical coupling effect of stress and seepage failure.The seepage failure will induce the final shear failure when the ratio of deviatoric stress to confining pressure is within 1.0-1.2; 2） A negative exponential permeability empirical model for the joint denoted by a newly defined principal stress function,which considers the coupling effect of confining pressure and axial pressure on the permeability,is established based on hydro-mechanical coupling experiments.3） The variation of the settlement before and after seepage failure is very different.The settlement before seepage failure changes very slowly,while it increases significantly after the seepage failure.4） The stress-strain relationship is of a strain softening type.5） Flow along the joint still follows Darcian flow rule.The results will provide an important theoretical basis for the further evaluation on the safety of the high earth and rockfill dam.展开更多
At the Mont Terri Underground Research Laboratory （Switzerland）, a field-scale investigation has been conducted in order to investigate the hydro-mechanical and chemical perturbations induced in the argilla- ceous f...At the Mont Terri Underground Research Laboratory （Switzerland）, a field-scale investigation has been conducted in order to investigate the hydro-mechanical and chemical perturbations induced in the argilla- ceous formation by forced ventilation through a tunnel. This experiment has been selected to be used for processing model development and validation in the international project DECOVALEX-2011. The con- ceptual and mathematical representation of the engineered void, which itself forms a major part of the experiment and is not simply a boundary condition, is the subject of this paper. A variety of approaches have been examined by the contributors to DECOVALEX and a summary of their findings is presented here. Two major aspects are discussed. Firstly, the approaches for the treatment of the surface condition at the porous media/tunnel interface are examined, with two equivalent but differing formulations successfully demonstrated. Secondly, approaches for representing the tunnel with associated air and water vapour movement, when coupled with the hydro-mechanical （HM） representation of the porous medium, are also examined. It is clearly demonstrated that, for the experimental conditions of the ventilation experiment （VE） that abstracted physical and empirical models of the tunnel, can be used successfully to represent the hydraulic behaviour of the tunnel and the hydraulic interaction between the tunnel and the surrounding rock mass.展开更多
Mechanical failure of materials adjacent to the production cavity and material disaggregation caused by fluid drag are considered as the most important parameters that affect sand production.In light of such factors,t...Mechanical failure of materials adjacent to the production cavity and material disaggregation caused by fluid drag are considered as the most important parameters that affect sand production.In light of such factors,the coupling of two mechanisms-mechanical instability and hydrodynamic erosion-is indispensable in order to model this phenomenon successfully.This paper examines the applicability of a coupled hydro-mechanical erosion criterion for simulating sand production using the finite element method.The porous medium was considered fully saturated.The onset of sanding and production of sand were predicted by coupling mechanical failure and subsequent erosion of the grain particles utilizing a sanding model.To consider the erosion process,the Papamichos and Stavropoulou(1998)’s sand erosion criterion was incorporated into the finite element code.Arbitrary Lagrangian-Eulerian(ALE)adaptive mesh approach was used to account for large amounts of erosive material loss.Besides,in order to address the problem of severe mesh distortion,the“mesh mapping technique”was employed.Sand production in a horizontal wellbore and in a field case was simulated to demonstrate capabilities of the proposed model.In addition,principal parameters affecting sand production,including in situ stresses,cohesion,perforation orientation,and drawdown were examined.The results indicated the efficiency of the model used in evaluation of sanding in the field.Parametric studies indicated that in situ stresses and formation cohesion could be considered as dominant factors affecting the amount of sand production.展开更多
The objective of this paper is to develop a methodology for calibration of a discrete element grain-based model(GBM)to replicate the hydro-mechanical properties of a brittle rock measured in the laboratory,and to appl...The objective of this paper is to develop a methodology for calibration of a discrete element grain-based model(GBM)to replicate the hydro-mechanical properties of a brittle rock measured in the laboratory,and to apply the calibrated model to simulating the formation of excavation damage zone(EDZ)around underground excavations.Firstly,a new cohesive crack model is implemented into the universal distinct element code(UDEC)to control the fracturing behaviour of materials under various loading modes.Next,a methodology for calibration of the components of the UDEC-Voronoi model is discussed.The role of connectivity of induced microcracks on increasing the permeability of laboratory-scale samples is investigated.The calibrated samples are used to investigate the influence of pore fluid pressure on weakening the drained strength of the laboratory-scale rock.The validity of the Terzaghi’s effective stress law for the drained peak strength of low-porosity rock is tested by performing a series of biaxial compression test simulations.Finally,the evolution of damage and pore pressure around two unsupported circular tunnels in crystalline granitic rock is studied.展开更多
The construction of moral education in primary and secondary schools is facingmany dilemmas,such as formalism,passivity,task-type,Tacitus trap and new media interaction.There is a high degree of coupling between the p...The construction of moral education in primary and secondary schools is facingmany dilemmas,such as formalism,passivity,task-type,Tacitus trap and new media interaction.There is a high degree of coupling between the progressive,stepped and interlocking supply of Problem Chain mode and the demand of moral education construction in primary and secondary schools.Problem Chain mode is helpful to solve the problems of teaching content,teaching effects and educational object in moral education construction in primary and secondary schools.展开更多
To explore the variation of permeability and deformation behaviors of a fractured rock mass in high water pressure,a high pressure permeability test（HPPT）,including measuring sensors of pore water pressure and displ...To explore the variation of permeability and deformation behaviors of a fractured rock mass in high water pressure,a high pressure permeability test（HPPT）,including measuring sensors of pore water pressure and displacement of the rock mass,was designed according to the hydrogeological condition of Heimifeng pumped storage power station.With the assumption of radial water flow pattern in the rock mass during the HPPT,a theoretical formula was presented to estimate the coefficient of permeability of the rock mass using water pressures in injection and measuring boreholes.The variation in permeability of the rock mass with the injected water pressure was studied according to the suggested formula.By fitting the relationship between the coefficient of permeability and the injected water pressure,a mathematical expression was obtained and used in the numerical simulations.For a better understanding of the relationship between the pore water pressure and the displacement of the rock mass,a 3D numerical method based on a coupled hydro-mechanical theory was employed to simulate the response of the rock mass during the test.By comparison of the calculated and measured data of pore water pressure and displacement,the deformation behaviors of the rock mass were analyzed.It is shown that the variation of displacement in the fractured rock mass is caused by water flow passing through it under high water pressure,and the rock deformation during the test could be calculated by using the coupled hydro-mechanical model.展开更多
Sandstone“injectite”intrusions are generally developed by the fluidization of weakly cemented sandstones and their subsequent injection into fractured reservoirs.In this work,a continuum coupled hydromechanical mode...Sandstone“injectite”intrusions are generally developed by the fluidization of weakly cemented sandstones and their subsequent injection into fractured reservoirs.In this work,a continuum coupled hydromechanical model TOUGH-FLAC3D is applied to simulate the discrete fracture network in large-scale sand injectite complexes.A sand production constitutive model is incorporated to consider the formation of sand through plastic deformation and its influence on evolution of fracture permeability.Overpressures in the fluidized sand slurry drives the injection with sand dikes intruded upwards,typically into previously low permeability“tight”mudstone formations.The contrast in poroelastic properties of the underlying weak sandstone and overlying injectite receptor directly affects the evolution of fracture aperture both during and after intrusion.Fluid drainage into the unconsolidated matrix may reduce the extent of fracture aperture growth,through the formation of shear bands.The results of this work have broad implications related to the emplacement of sandstone intrusions and subsequent hydrocarbon accumulation,maturation and then production.展开更多
Callovo-Oxfordian （COx） argillite obtained from the excavation of high-level radioactive waste geological disposal has been evaluated as an alternative sealing/backfill material in France. This paper presents an exp...Callovo-Oxfordian （COx） argillite obtained from the excavation of high-level radioactive waste geological disposal has been evaluated as an alternative sealing/backfill material in France. This paper presents an experimental investigation into the hydro-mechanical behaviour of compacted crushed COx argillite. A series of oedorneter compressive tests including various loading-unloading cycles were conducted on COx argillite powders at different initial water contents. After reaching the desired dry density （2.0 Mg/m^3）, the vertical stress was reduced to different levels （7.0 and 0.5 MPa） and the compacted sample was then flooded under constant volume conditions while measuring the changes in the vertical stress. It was found that the initial water content significantly affects the compressive behaviour. The measured saturated hydraulic conductivity is less than 1×10^-10m/s.展开更多
A methodology for identifying and calibrating the material parameters for a coupled hydro-mechanical problem is presented in this pape r.For validation purpose,a laboratory-based water infiltration test was numericall...A methodology for identifying and calibrating the material parameters for a coupled hydro-mechanical problem is presented in this pape r.For validation purpose,a laboratory-based water infiltration test was numerically simulated using finite element method(FEM).The test was conducted using a self-designed column-type experimental device,which mimicked the wetting process of a candidate backfill material in a nuclear waste repository.The real-time measurements of key state variables(e.g.water content,relative humidity,temperature,and total stresses)were performed with the monitoring sensors along the height of cylindrical soil sample.For numerical simulation,the modified Barcelona Basic Model(BBM)along with soil-water retention model for compacted bentonite was used.It shows that the identified model parameters successfully captured the moisture migration process under an applied hydraulic gradient in a bentonite-based compacted soil sample.A comparison between the measured and predicted values of total stresses both in axial and lateral directions along with other state variables revealed that heterogeneous moisture content was distributed along the hydration-path,resulting in non-uniform stress-deformation characteristics of soil.展开更多
Explicit solution techniques have been widely used in geotechnical engineering for simulating the coupled hydro-mechanical(H-M) interaction of fluid flow and deformation induced by structures built above and under sat...Explicit solution techniques have been widely used in geotechnical engineering for simulating the coupled hydro-mechanical(H-M) interaction of fluid flow and deformation induced by structures built above and under saturated ground, i.e. circular footing and deep tunnel. However, the technique is only conditionally stable and requires small time steps, portending its inefficiency for simulating large-scale H-M problems. To improve its efficiency, the unconditionally stable alternating direction explicit(ADE)scheme could be used to solve the flow problem. The standard ADE scheme, however, is only moderately accurate and is restricted to uniform grids and plane strain flow conditions. This paper aims to remove these drawbacks by developing a novel high-order ADE scheme capable of solving flow problems in nonuniform grids and under axisymmetric conditions. The new scheme is derived by performing a fourthorder finite difference(FD) approximation to the spatial derivatives of the axisymmetric fluid-diffusion equation in a non-uniform grid configuration. The implicit Crank-Nicolson technique is then applied to the resulting approximation, and the subsequent equation is split into two alternating direction sweeps,giving rise to a new axisymmetric ADE scheme. The pore pressure solutions from the new scheme are then sequentially coupled with an existing geomechanical simulator in the computer code fast Lagrangian analysis of continua(FLAC). This coupling procedure is called the sequentially-explicit coupling technique based on the fourth-order axisymmetric ADE scheme or SEA-4-AXI. Application of SEA-4-AXI for solving axisymmetric consolidation of a circular footing and of advancing tunnel in deep saturated ground shows that SEA-4-AXI reduces computer runtime up to 42%-50% that of FLAC’s basic scheme without numerical instability. In addition, it produces high numerical accuracy of the H-M solutions with average percentage difference of only 0.5%-1.8%.展开更多
The present work deals with the calculation of transition probability between two diabatic potentials coupled by any arbitrary coupling. The method presented in this work is applicable to any type of coupling while fo...The present work deals with the calculation of transition probability between two diabatic potentials coupled by any arbitrary coupling. The method presented in this work is applicable to any type of coupling while for numerical calculations we have assumed the arbitrary coupling as Gaussian coupling. This arbitrary coupling is expressed as a collection of Dirac delta functions and by the use of the transfer matrix technique the transition probability from one diabatic potential to another diabatic potential is calculated. We examine our approach by considering the case of two constant potentials coupled by Gaussian coupling as an arbitrary coupling.展开更多
The material point method(MPM)has been gaining increasing popularity as an appropriate approach to the solution of coupled hydro-mechanical problems involving large deformation.In this paper,we survey the current stat...The material point method(MPM)has been gaining increasing popularity as an appropriate approach to the solution of coupled hydro-mechanical problems involving large deformation.In this paper,we survey the current state-of-the-art in the MPM simulation of hydro-mechanical behaviour in two-phase porous geomaterials.The review covers the recent advances and developments in the MPM and their extensions to capture the coupled hydro-mechanical problems involving large deformations.The focus of this review is aiming at providing a clear picture of what has or has not been developed or implemented for simulating two-phase coupled large deformation problems,which will provide some direct reference for both practitioners and researchers.展开更多
Significant progress has been made in mixed boundary-value problems associated with three-dimensional（3D） crack and contact analyses of advanced materials featuring more complexities compared to the conventional iso...Significant progress has been made in mixed boundary-value problems associated with three-dimensional（3D） crack and contact analyses of advanced materials featuring more complexities compared to the conventional isotropic elastic materials.These include material anisotropy and multifield coupling,two typical characteristics of most current multifunctional materials.In this paper we try to present a state-of-the-art description of 3D exact/analytical solutions derived for crack and contact problems of elastic solids with both transverse isotropy and multifield coupling in the latest decade by the potential theory method in the spirit of V.I.Fabrikant.whose ingenious breakthrough brings new vigor and vitality to the old research subject of classical potential theory.We are particularly interested in crack and contact problems with certain nonlinear features.Emphasis is also placed on the coupling between the temperature field（or the like） and other physical fields（e.g.,elastic,electric,and magnetic fields）.We further highlight the practical significance of 3D contact solutions,in particular in applications related to modern scanning probe microscopes.展开更多
A new isospectral problem is firstly presented, then we derive integrable system of soliton hierarchy. Also we obtain new integrable couplings of the generalized Kaup-Newell soliton equations hierarchy and its Hamilto...A new isospectral problem is firstly presented, then we derive integrable system of soliton hierarchy. Also we obtain new integrable couplings of the generalized Kaup-Newell soliton equations hierarchy and its Hamiltonian structures by using Tu scheme and the quadratic-form identity. The method can be generalized to other soliton hierarchy.展开更多
This paper derives new discrete integrable system based on discrete isospectral problem. It shows that the hierarchy is completely integrable in the Liouville sense and possesses bi-Hamiltonian structure. Finally, int...This paper derives new discrete integrable system based on discrete isospectral problem. It shows that the hierarchy is completely integrable in the Liouville sense and possesses bi-Hamiltonian structure. Finally, integrable couplings of the obtained system is given by means of semi-direct sums of Lie algebras.展开更多
In order to investigate the hydro-mechanical （HM） and chemical perturbations induced in an argillaceous formation by forced ventilation during the operational period of a nuclear waste repository, a specific experim...In order to investigate the hydro-mechanical （HM） and chemical perturbations induced in an argillaceous formation by forced ventilation during the operational period of a nuclear waste repository, a specific experiment has been performed in a tunnel, at Mont Terri Underground Research Laboratory （URL） in Switzerland. This experiment has been selected in the international project DECOVALEX for model vali- dation and the numerical simulation of this ventilation experiment （VE） is the object of the present paper. Since the argillaceous rock exhibits anisotropic properties, particular attention is given to the evaluation of the effects of various anisotropic features on the predicted results. In situ measurements such as relative humidity （RH）, global water mass extracted, pore water pressure, water content, and relative displace- ments are compared to predictions using both isotropic and anisotropic parameters. Water permeability anisotropy is shown to be the most influencing parameter by far, whereas in situ stress anisotropy has an effect only during the excavation phase. The anisotropy for mechanical parameterization has also some influence, in particular through HM couplings. These HM couplings have the potential to be very significant in terms of providing confidence in describing the experimental observation, and should be considered for further investigation.展开更多
This paper reports a multiscale analysis method to predict the thermomechanical coupling performance of composite structures with quasi-periodic properties.In these material structures,the configurations are periodic,...This paper reports a multiscale analysis method to predict the thermomechanical coupling performance of composite structures with quasi-periodic properties.In these material structures,the configurations are periodic,and the material coefficients are quasi-periodic,i.e.,they depend not only on the microscale information but also on the macro location.Also,a mutual interaction between displacement and temperature fields is considered in the problem,which is our particular interest in this study.The multiscale asymptotic expansions of the temperature and displacement fields are constructed and associated error estimation in nearly pointwise sense is presented.Then,a finite element-difference algorithm based on the multiscale analysis method is brought forward in detail.Finally,some numerical examples are given.And the numerical results show that the multiscale method presented in this paper is effective and reliable to study the nonlinear thermo-mechanical coupling problem of composite structures with quasiperiodic properties.展开更多
This paper presents an overview of the recent progress of potential theory method in the analysis of mixed boundary value problems mainly stemming from three-dimensional crack or contact problems of multi-field couple...This paper presents an overview of the recent progress of potential theory method in the analysis of mixed boundary value problems mainly stemming from three-dimensional crack or contact problems of multi-field coupled media. This method was used to derive a series of exact three dimensional solutions which should be of great theoretical significance because most of them usually cannot be derived by other methods such as the transform method and the trial-and-error method. Further, many solutions are obtained in terms of elementary functions that enable us to treat more complicated problems easily. It is pointed out here that the method is usually only applicable to media characterizing transverse isotropy, from which, however, the results for the isotropic case can be readily obtained.展开更多
We propose a multi-field implicit finite element method for analyzing the electromechanical behavior of dielectric elastomers. This method is based on a four-field variational principle, which includes displacement an...We propose a multi-field implicit finite element method for analyzing the electromechanical behavior of dielectric elastomers. This method is based on a four-field variational principle, which includes displacement and electric potential for the electromechanical coupling analysis, and additional independent fields to address the incompressible constraint of the hyperelastic material. Linearization of the variational form and finite element discretization are adopted for the numerical implementation. A general FEM program framework is devel- oped using C＋＋ based on the open-source finite element library deal.II to implement this proposed algorithm. Numerical examples demonstrate the accuracy, convergence properties, mesh-independence properties, and scalability of this method. We also use the method for eigenvalue analysis of a dielectric elastomer actuator subject to electromechanical loadings. Our finite element implementation is available as an online supplementary material.展开更多
基金We acknowledge the combined support from the National Natural Science Foundation of China(Grant Nos.52039007 and 42102325)Tiandi Science and Technology Co.,Ltd.(Grant No.2022-2-TD-MS012).
文摘Rock mass is a fractured porous medium usually subjected to complex geostress and fluid pressure simultaneously.Moreover,the properties of rock mass change in time and space due to mining-induced fractures.Therefore,it is always challenging to accurately measure rock mass properties.In this study,a three-dimensional(3D)microseismic(MS)data-driven damage model for jointed rock mass under hydro-mechanical coupling conditions is proposed.It is a 3D finite element model that takes seepage,damage and stress field effects into account jointly.Multiple factors(i.e.joints,water and microseismicity)are used to optimize the rock mass mechanical parameters at different scales.The model is applied in Shirengou iron mine to study the damage evolution of rock mass and assess the crown pillar stability during the transition from open-pit to underground mining.It is found that the damage pattern is mostly controlled by the structure,water and rock mass parameters.The damage pattern is evidently different from the two-dimensional result and is more consistent with the field observations.This difference is caused by the MS-derived damage acting on the rock mass.MS data are responsible for gradually correcting the damage zone,changing the direction in which it expands,and promoting it to evolve close to reality.For the crown pillar,the proposed model yields a more trustworthy safety factor.In order to guarantee the stability of the pillar,it is suggested to take waterproof and reinforcement measures in areas with a high degree of damage.
基金Projects(51009053,51079039)supported by the National Natural Science Foundation of ChinaProject(20100094120004)supported by the Doctoral Program of Higher Education of China
文摘The joint of clay core-wall and concrete cut-off wall is one of the weakest parts in high earth and rockftll dams.A kind of highly plastic clay is always fixed on the joint to fit the large shear deformation between clay core-wall and concrete cut-offwall,so the hydro-mechanical coupling mechanisms on the joint under high stress,high hydraulic gradient,and large shear deformation are of great importance for the evaluation of dam safety.The hydro-mechanical coupling characteristics of the joint of the highly plastic clay and the concrete cut-off wall in a high earth and rockfill dam in China were studied by using a newly designed soil-structure contact erosion apparatus.The experimental results indicate that：1） Shear failure on the joint is due to the hydro-mechanical coupling effect of stress and seepage failure.The seepage failure will induce the final shear failure when the ratio of deviatoric stress to confining pressure is within 1.0-1.2; 2） A negative exponential permeability empirical model for the joint denoted by a newly defined principal stress function,which considers the coupling effect of confining pressure and axial pressure on the permeability,is established based on hydro-mechanical coupling experiments.3） The variation of the settlement before and after seepage failure is very different.The settlement before seepage failure changes very slowly,while it increases significantly after the seepage failure.4） The stress-strain relationship is of a strain softening type.5） Flow along the joint still follows Darcian flow rule.The results will provide an important theoretical basis for the further evaluation on the safety of the high earth and rockfill dam.
基金the context of the international DECOVALEX Project (DEmonstration of COupled models and their VALidation against EXperiments)EC project NF-PRO (Contract number FI6W-CT-2003-02389) under the coordination of ENRESA (Empresa Nacional de Residuos Radiactivos)
文摘At the Mont Terri Underground Research Laboratory （Switzerland）, a field-scale investigation has been conducted in order to investigate the hydro-mechanical and chemical perturbations induced in the argilla- ceous formation by forced ventilation through a tunnel. This experiment has been selected to be used for processing model development and validation in the international project DECOVALEX-2011. The con- ceptual and mathematical representation of the engineered void, which itself forms a major part of the experiment and is not simply a boundary condition, is the subject of this paper. A variety of approaches have been examined by the contributors to DECOVALEX and a summary of their findings is presented here. Two major aspects are discussed. Firstly, the approaches for the treatment of the surface condition at the porous media/tunnel interface are examined, with two equivalent but differing formulations successfully demonstrated. Secondly, approaches for representing the tunnel with associated air and water vapour movement, when coupled with the hydro-mechanical （HM） representation of the porous medium, are also examined. It is clearly demonstrated that, for the experimental conditions of the ventilation experiment （VE） that abstracted physical and empirical models of the tunnel, can be used successfully to represent the hydraulic behaviour of the tunnel and the hydraulic interaction between the tunnel and the surrounding rock mass.
基金funded by the Iran National Science Foundation(INSF)(Grant No.96001589).
文摘Mechanical failure of materials adjacent to the production cavity and material disaggregation caused by fluid drag are considered as the most important parameters that affect sand production.In light of such factors,the coupling of two mechanisms-mechanical instability and hydrodynamic erosion-is indispensable in order to model this phenomenon successfully.This paper examines the applicability of a coupled hydro-mechanical erosion criterion for simulating sand production using the finite element method.The porous medium was considered fully saturated.The onset of sanding and production of sand were predicted by coupling mechanical failure and subsequent erosion of the grain particles utilizing a sanding model.To consider the erosion process,the Papamichos and Stavropoulou(1998)’s sand erosion criterion was incorporated into the finite element code.Arbitrary Lagrangian-Eulerian(ALE)adaptive mesh approach was used to account for large amounts of erosive material loss.Besides,in order to address the problem of severe mesh distortion,the“mesh mapping technique”was employed.Sand production in a horizontal wellbore and in a field case was simulated to demonstrate capabilities of the proposed model.In addition,principal parameters affecting sand production,including in situ stresses,cohesion,perforation orientation,and drawdown were examined.The results indicated the efficiency of the model used in evaluation of sanding in the field.Parametric studies indicated that in situ stresses and formation cohesion could be considered as dominant factors affecting the amount of sand production.
文摘The objective of this paper is to develop a methodology for calibration of a discrete element grain-based model(GBM)to replicate the hydro-mechanical properties of a brittle rock measured in the laboratory,and to apply the calibrated model to simulating the formation of excavation damage zone(EDZ)around underground excavations.Firstly,a new cohesive crack model is implemented into the universal distinct element code(UDEC)to control the fracturing behaviour of materials under various loading modes.Next,a methodology for calibration of the components of the UDEC-Voronoi model is discussed.The role of connectivity of induced microcracks on increasing the permeability of laboratory-scale samples is investigated.The calibrated samples are used to investigate the influence of pore fluid pressure on weakening the drained strength of the laboratory-scale rock.The validity of the Terzaghi’s effective stress law for the drained peak strength of low-porosity rock is tested by performing a series of biaxial compression test simulations.Finally,the evolution of damage and pore pressure around two unsupported circular tunnels in crystalline granitic rock is studied.
文摘The construction of moral education in primary and secondary schools is facingmany dilemmas,such as formalism,passivity,task-type,Tacitus trap and new media interaction.There is a high degree of coupling between the progressive,stepped and interlocking supply of Problem Chain mode and the demand of moral education construction in primary and secondary schools.Problem Chain mode is helpful to solve the problems of teaching content,teaching effects and educational object in moral education construction in primary and secondary schools.
文摘To explore the variation of permeability and deformation behaviors of a fractured rock mass in high water pressure,a high pressure permeability test（HPPT）,including measuring sensors of pore water pressure and displacement of the rock mass,was designed according to the hydrogeological condition of Heimifeng pumped storage power station.With the assumption of radial water flow pattern in the rock mass during the HPPT,a theoretical formula was presented to estimate the coefficient of permeability of the rock mass using water pressures in injection and measuring boreholes.The variation in permeability of the rock mass with the injected water pressure was studied according to the suggested formula.By fitting the relationship between the coefficient of permeability and the injected water pressure,a mathematical expression was obtained and used in the numerical simulations.For a better understanding of the relationship between the pore water pressure and the displacement of the rock mass,a 3D numerical method based on a coupled hydro-mechanical theory was employed to simulate the response of the rock mass during the test.By comparison of the calculated and measured data of pore water pressure and displacement,the deformation behaviors of the rock mass were analyzed.It is shown that the variation of displacement in the fractured rock mass is caused by water flow passing through it under high water pressure,and the rock deformation during the test could be calculated by using the coupled hydro-mechanical model.
基金the financial support from the Laboratory of Coal Resources and Safe Mining(China University of Mining and Technology,Beijing)(Grant No.SKLCRSM16KFC01)。
文摘Sandstone“injectite”intrusions are generally developed by the fluidization of weakly cemented sandstones and their subsequent injection into fractured reservoirs.In this work,a continuum coupled hydromechanical model TOUGH-FLAC3D is applied to simulate the discrete fracture network in large-scale sand injectite complexes.A sand production constitutive model is incorporated to consider the formation of sand through plastic deformation and its influence on evolution of fracture permeability.Overpressures in the fluidized sand slurry drives the injection with sand dikes intruded upwards,typically into previously low permeability“tight”mudstone formations.The contrast in poroelastic properties of the underlying weak sandstone and overlying injectite receptor directly affects the evolution of fracture aperture both during and after intrusion.Fluid drainage into the unconsolidated matrix may reduce the extent of fracture aperture growth,through the formation of shear bands.The results of this work have broad implications related to the emplacement of sandstone intrusions and subsequent hydrocarbon accumulation,maturation and then production.
基金Supported by the French National Radioactive Waste Management Agency
文摘Callovo-Oxfordian （COx） argillite obtained from the excavation of high-level radioactive waste geological disposal has been evaluated as an alternative sealing/backfill material in France. This paper presents an experimental investigation into the hydro-mechanical behaviour of compacted crushed COx argillite. A series of oedorneter compressive tests including various loading-unloading cycles were conducted on COx argillite powders at different initial water contents. After reaching the desired dry density （2.0 Mg/m^3）, the vertical stress was reduced to different levels （7.0 and 0.5 MPa） and the compacted sample was then flooded under constant volume conditions while measuring the changes in the vertical stress. It was found that the initial water content significantly affects the compressive behaviour. The measured saturated hydraulic conductivity is less than 1×10^-10m/s.
基金the German Research Foundation(DFG)for the financial support(Grant No.SCHA 675/17-1)。
文摘A methodology for identifying and calibrating the material parameters for a coupled hydro-mechanical problem is presented in this pape r.For validation purpose,a laboratory-based water infiltration test was numerically simulated using finite element method(FEM).The test was conducted using a self-designed column-type experimental device,which mimicked the wetting process of a candidate backfill material in a nuclear waste repository.The real-time measurements of key state variables(e.g.water content,relative humidity,temperature,and total stresses)were performed with the monitoring sensors along the height of cylindrical soil sample.For numerical simulation,the modified Barcelona Basic Model(BBM)along with soil-water retention model for compacted bentonite was used.It shows that the identified model parameters successfully captured the moisture migration process under an applied hydraulic gradient in a bentonite-based compacted soil sample.A comparison between the measured and predicted values of total stresses both in axial and lateral directions along with other state variables revealed that heterogeneous moisture content was distributed along the hydration-path,resulting in non-uniform stress-deformation characteristics of soil.
基金the support from the University Transportation Center for Underground Transportation Infrastructure at the Colorado School of Mines for partially funding this research under Grant No. 69A3551747118 of the Fixing America's Surface Transportation Act (FAST Act) of U.S. DoT FY2016
文摘Explicit solution techniques have been widely used in geotechnical engineering for simulating the coupled hydro-mechanical(H-M) interaction of fluid flow and deformation induced by structures built above and under saturated ground, i.e. circular footing and deep tunnel. However, the technique is only conditionally stable and requires small time steps, portending its inefficiency for simulating large-scale H-M problems. To improve its efficiency, the unconditionally stable alternating direction explicit(ADE)scheme could be used to solve the flow problem. The standard ADE scheme, however, is only moderately accurate and is restricted to uniform grids and plane strain flow conditions. This paper aims to remove these drawbacks by developing a novel high-order ADE scheme capable of solving flow problems in nonuniform grids and under axisymmetric conditions. The new scheme is derived by performing a fourthorder finite difference(FD) approximation to the spatial derivatives of the axisymmetric fluid-diffusion equation in a non-uniform grid configuration. The implicit Crank-Nicolson technique is then applied to the resulting approximation, and the subsequent equation is split into two alternating direction sweeps,giving rise to a new axisymmetric ADE scheme. The pore pressure solutions from the new scheme are then sequentially coupled with an existing geomechanical simulator in the computer code fast Lagrangian analysis of continua(FLAC). This coupling procedure is called the sequentially-explicit coupling technique based on the fourth-order axisymmetric ADE scheme or SEA-4-AXI. Application of SEA-4-AXI for solving axisymmetric consolidation of a circular footing and of advancing tunnel in deep saturated ground shows that SEA-4-AXI reduces computer runtime up to 42%-50% that of FLAC’s basic scheme without numerical instability. In addition, it produces high numerical accuracy of the H-M solutions with average percentage difference of only 0.5%-1.8%.
文摘The present work deals with the calculation of transition probability between two diabatic potentials coupled by any arbitrary coupling. The method presented in this work is applicable to any type of coupling while for numerical calculations we have assumed the arbitrary coupling as Gaussian coupling. This arbitrary coupling is expressed as a collection of Dirac delta functions and by the use of the transfer matrix technique the transition probability from one diabatic potential to another diabatic potential is calculated. We examine our approach by considering the case of two constant potentials coupled by Gaussian coupling as an arbitrary coupling.
基金The financial supports from National Outstanding Youth Science Fund Project of National Natural Science Foundation of China(Grant No.52022112)the International Postdoctoral Exchange Fellowship Program(Talent-Introduction Program,Grant No.YJ20220219)。
文摘The material point method(MPM)has been gaining increasing popularity as an appropriate approach to the solution of coupled hydro-mechanical problems involving large deformation.In this paper,we survey the current state-of-the-art in the MPM simulation of hydro-mechanical behaviour in two-phase porous geomaterials.The review covers the recent advances and developments in the MPM and their extensions to capture the coupled hydro-mechanical problems involving large deformations.The focus of this review is aiming at providing a clear picture of what has or has not been developed or implemented for simulating two-phase coupled large deformation problems,which will provide some direct reference for both practitioners and researchers.
基金supported by the National Natural Science Foundation of China(Grant 11321202)the Specialized Research Fund for the Doctoral Program of Higher Education(Grant 20130101110120)
文摘Significant progress has been made in mixed boundary-value problems associated with three-dimensional（3D） crack and contact analyses of advanced materials featuring more complexities compared to the conventional isotropic elastic materials.These include material anisotropy and multifield coupling,two typical characteristics of most current multifunctional materials.In this paper we try to present a state-of-the-art description of 3D exact/analytical solutions derived for crack and contact problems of elastic solids with both transverse isotropy and multifield coupling in the latest decade by the potential theory method in the spirit of V.I.Fabrikant.whose ingenious breakthrough brings new vigor and vitality to the old research subject of classical potential theory.We are particularly interested in crack and contact problems with certain nonlinear features.Emphasis is also placed on the coupling between the temperature field（or the like） and other physical fields（e.g.,elastic,electric,and magnetic fields）.We further highlight the practical significance of 3D contact solutions,in particular in applications related to modern scanning probe microscopes.
基金Supported by the Natural Science Foundation of China under Grant Nos. 61072147, 11071159, and 10971031by the Natural Science Foundation of Shanghai and Zhejiang Province under Grant Nos. 09ZR1410800 and Y6100791+1 种基金the Shanghai Shuguang Tracking Project under Grant No. 08GG01the Shanghai Leading Academic Discipline Project under Grant No. J50101
文摘A new isospectral problem is firstly presented, then we derive integrable system of soliton hierarchy. Also we obtain new integrable couplings of the generalized Kaup-Newell soliton equations hierarchy and its Hamiltonian structures by using Tu scheme and the quadratic-form identity. The method can be generalized to other soliton hierarchy.
文摘This paper derives new discrete integrable system based on discrete isospectral problem. It shows that the hierarchy is completely integrable in the Liouville sense and possesses bi-Hamiltonian structure. Finally, integrable couplings of the obtained system is given by means of semi-direct sums of Lie algebras.
基金the context of the international DECOVALEX Project (DEmonstration of Coupled models and their VALidation against EXperiments)Quintessa Ltd. and University of Edinburgh were supported by the Nuclear Decommissioning Authority (NDA), UK+2 种基金CEA was supported by Institut de Radioprotection et de Sreté Nucléaire(IRSN)The Japan Atomic Energy Agency (JAEA) and the Institute of Rock and Soil Mechanics, Chinese Academy of Sciences (CAS)funded DECOVALEX and participated in the workEC project NF-PRO (Contract number FI6W-CT-2003-02389) under the coordination of ENRESA (Empresa Nacional de Residuos Radiactivos)
文摘In order to investigate the hydro-mechanical （HM） and chemical perturbations induced in an argillaceous formation by forced ventilation during the operational period of a nuclear waste repository, a specific experiment has been performed in a tunnel, at Mont Terri Underground Research Laboratory （URL） in Switzerland. This experiment has been selected in the international project DECOVALEX for model vali- dation and the numerical simulation of this ventilation experiment （VE） is the object of the present paper. Since the argillaceous rock exhibits anisotropic properties, particular attention is given to the evaluation of the effects of various anisotropic features on the predicted results. In situ measurements such as relative humidity （RH）, global water mass extracted, pore water pressure, water content, and relative displace- ments are compared to predictions using both isotropic and anisotropic parameters. Water permeability anisotropy is shown to be the most influencing parameter by far, whereas in situ stress anisotropy has an effect only during the excavation phase. The anisotropy for mechanical parameterization has also some influence, in particular through HM couplings. These HM couplings have the potential to be very significant in terms of providing confidence in describing the experimental observation, and should be considered for further investigation.
基金financially supported by the National Natural Science Foundation of China(11501449)the Fundamental Research Funds for the Central Universities(3102017zy043)+2 种基金the China Postdoctoral Science Foundation(2016T91019)the fund of the State Key Laboratory of Solidification Processing in NWPU(SKLSP201628)the Scientific Research Program Funded by Shaanxi Provincial Education Department(14JK1353).
文摘This paper reports a multiscale analysis method to predict the thermomechanical coupling performance of composite structures with quasi-periodic properties.In these material structures,the configurations are periodic,and the material coefficients are quasi-periodic,i.e.,they depend not only on the microscale information but also on the macro location.Also,a mutual interaction between displacement and temperature fields is considered in the problem,which is our particular interest in this study.The multiscale asymptotic expansions of the temperature and displacement fields are constructed and associated error estimation in nearly pointwise sense is presented.Then,a finite element-difference algorithm based on the multiscale analysis method is brought forward in detail.Finally,some numerical examples are given.And the numerical results show that the multiscale method presented in this paper is effective and reliable to study the nonlinear thermo-mechanical coupling problem of composite structures with quasiperiodic properties.
基金Project (No. 10372088) supported by the National Natural Science Foundation of China
文摘This paper presents an overview of the recent progress of potential theory method in the analysis of mixed boundary value problems mainly stemming from three-dimensional crack or contact problems of multi-field coupled media. This method was used to derive a series of exact three dimensional solutions which should be of great theoretical significance because most of them usually cannot be derived by other methods such as the transform method and the trial-and-error method. Further, many solutions are obtained in terms of elementary functions that enable us to treat more complicated problems easily. It is pointed out here that the method is usually only applicable to media characterizing transverse isotropy, from which, however, the results for the isotropic case can be readily obtained.
基金the support under A*STAR SERC grant (132-183-0025)
文摘We propose a multi-field implicit finite element method for analyzing the electromechanical behavior of dielectric elastomers. This method is based on a four-field variational principle, which includes displacement and electric potential for the electromechanical coupling analysis, and additional independent fields to address the incompressible constraint of the hyperelastic material. Linearization of the variational form and finite element discretization are adopted for the numerical implementation. A general FEM program framework is devel- oped using C＋＋ based on the open-source finite element library deal.II to implement this proposed algorithm. Numerical examples demonstrate the accuracy, convergence properties, mesh-independence properties, and scalability of this method. We also use the method for eigenvalue analysis of a dielectric elastomer actuator subject to electromechanical loadings. Our finite element implementation is available as an online supplementary material.
