Understanding the rock mass response to excavation and thermal loading and improving the capability of the numerical models for simulating the progressive failure process of brittle rocks are important for safety asse...Understanding the rock mass response to excavation and thermal loading and improving the capability of the numerical models for simulating the progressive failure process of brittle rocks are important for safety assessment and optimization design of nuclear waste repositories.The international cooperative DECOVALEX-2011 project provides a platform for development,validation and comparison of numerical models,in which the sp pillar stability experiment(APSE) was selected as the modeling target for Task B.This paper presents the modeling results of Wuhan University(WHU) team for stages 1 and 2 of Task B by using a coupled thermo-mechanical model within the framework of continuum mechanics.The rock mass response to excavation is modeled with linear elastic,elastoplastic and brittle-plastic models,while the response to heating is modeled with a coupled thermo-elastic model.The capabilities and limitations of the model for representation of the thermo-mechanical responses of the rock pillar are discussed by comparing the modeling results with experimental observations.The results may provide a helpful reference for the stability and safety assessment of the hard granite host rock in China's Beishan preselected area for high-level radioactive waste disposal.展开更多
The shear failure of intact rock under thermo-mechanical(TM)coupling conditions is common,such as in enhanced geothermal mining and deep mine construction.Under the effect of a continuous engineering disturbance,shear...The shear failure of intact rock under thermo-mechanical(TM)coupling conditions is common,such as in enhanced geothermal mining and deep mine construction.Under the effect of a continuous engineering disturbance,shear-formed fractures are prone to secondary instability,posing a severe threat to deep engineering.Although numerous studies regarding three-dimensional(3D)morphologies of fracture surfaces have been conducted,the understanding of shear-formed fractures under TM coupling conditions is limited.In this study,direct shear tests of intact granite under various TM coupling conditions were conducted,followed by 3D laser scanning tests of shear-formed fractures.Test results demonstrated that the peak shear strength of intact granite is positively correlated with the normal stress,whereas it is negatively correlated with the temperature.The internal friction angle and cohesion of intact granite significantly decrease with an increase in the temperature.The anisotropy,roughness value,and height of the asperities on the fracture surfaces are reduced as the normal stress increases,whereas their variation trends are the opposite as the temperature increases.The macroscopic failure mode of intact granite under TM coupling conditions is dominated by mixed tensileeshear and shear failures.As the normal stress increases,intragranular fractures are developed ranging from a local to a global distribution,and the macroscopic failure mode of intact granite changes from mixed tensileeshear to shear failure.Finally,3D morphological characteristics of the asperities on the shear-formed fracture surfaces were analyzed,and a quadrangular pyramid conceptual model representing these asperities was proposed and sufficiently verified.展开更多
This paper, with a finite element method, studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations. With this fully coupled model, the...This paper, with a finite element method, studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations. With this fully coupled model, the rigid structure is taken as "fictitious" fluid with zero strain rate. Both fluid and structure are described by velocity and pressure. The whole domain, including fluid region and structure region, is modeled by the incompressible Navier-Stokes equations which are discretized with fixed Eulerian mesh. However, to keep the structure' s rigid body shape and behavior, a rigid body constraint is enforced on the "fictitious" fluid domain by use of the Distributed Lagrange Multipher/Fictitious Domain (DLM/ FD) method which is originally introduced to solve particulate flow problems by Glowinski et al. For the verification of the model presented herein, a 2D numerical wave tank is established to simulate small amplitude wave propagations, and then numerical results are compared with analytical solutions. Finally, a 2D example of fluid-structure interaction under wave dynamic forces provides convincing evidences for the method excellent solution quality and fidelity.展开更多
Tension Leg Platform(TLP)in deepwater oil and gas field development usually consists of a hull,tendons,and top tension risers(TTRs).To maintain its top tension,each TTR is connected with a tensioner system to the hull...Tension Leg Platform(TLP)in deepwater oil and gas field development usually consists of a hull,tendons,and top tension risers(TTRs).To maintain its top tension,each TTR is connected with a tensioner system to the hull.Owing to the complicated configuration of the tensioners,the hull and TTRs form a strong coupled system.Traditionally,some simplified tensioner models are applied to analyze the TLP structures.There is a large discrepancy between their analysis results and the actual mechanism behaviors of a tensioner.It is very necessary to develop a more detailed tensioner model to consider the coupling effects between TLP and TTRs.In the present study,a fully coupled TLP hull-TTR system for hydrodynamic numerical simulation is established.A specific hydraulic pneumatic tensioner is modeled by considering 4 cylinders.The production TTR model is stacked up by specific riser joints.The simulation is also extended to analyze an array of TTRs.Different regular and irregular waves are considered.The behaviors of different cylinders are presented.The results show that it is important to consider the specific configurations of the tensioner and TTRs,which may lead to obviously different response behaviors,compared with those from a simplified model.展开更多
This paper developed a statistical damage constitutive model for deep rock by considering the effects of external load and thermal treatment temperature based on the distortion energy.The model parameters were determi...This paper developed a statistical damage constitutive model for deep rock by considering the effects of external load and thermal treatment temperature based on the distortion energy.The model parameters were determined through the extremum features of stress−strain curve.Subsequently,the model predictions were compared with experimental results of marble samples.It is found that when the treatment temperature rises,the coupling damage evolution curve shows an S-shape and the slope of ascending branch gradually decreases during the coupling damage evolution process.At a constant temperature,confining pressure can suppress the expansion of micro-fractures.As the confining pressure increases the rock exhibits ductility characteristics,and the shape of coupling damage curve changes from an S-shape into a quasi-parabolic shape.This model can well characterize the influence of high temperature on the mechanical properties of deep rock and its brittleness-ductility transition characteristics under confining pressure.Also,it is suitable for sandstone and granite,especially in predicting the pre-peak stage and peak stress of stress−strain curve under the coupling action of confining pressure and high temperature.The relevant results can provide a reference for further research on the constitutive relationship of rock-like materials and their engineering applications.展开更多
The process of plane strain compression of 06Cr19Ni9NbN steel was carried out in the temperature range of 1000-1200℃ and the reduction ratio range of 10-50%.Combining the finite element numerical simulation,a new con...The process of plane strain compression of 06Cr19Ni9NbN steel was carried out in the temperature range of 1000-1200℃ and the reduction ratio range of 10-50%.Combining the finite element numerical simulation,a new constitutive model of thermal deformation was established, which provides the theoretical basis to optimize the plarie strain compression process of the steel.The temperature and grain size at different regions were achieved by experiment and simulation,respectively.According to the results,the mathematical models of stress and temperature during the plane strain compression were established by mathematical analysis method.The new temperature models were established in three regions,respectively,and the stress models took account of the variation of temperature and'st^rain rate.Finally,by comparing the results of calculation,numerical simulation and experiment,the accuracy and validity of these mathematical models were verified.展开更多
Multi-fidelity simulations incorporate computational fluid dynamics(CFD) models into a thermodynamic model,enabling the simulation of the overall performance of an entire gas turbine with high-fidelity components.Trad...Multi-fidelity simulations incorporate computational fluid dynamics(CFD) models into a thermodynamic model,enabling the simulation of the overall performance of an entire gas turbine with high-fidelity components.Traditional iterative coupled methods rely on characteristic maps,while fully coupled methods directly incorporate high-fidelity simulations.However,fully coupled methods face challenges in simulating rotating components,including weak convergence and complex implementation.To address these challenges,a fully coupled method with logarithmic transformations was developed to directly integrate high-fidelity CFD models of multiple rotating components.The developed fully coupled method was then applied to evaluate the overall performance of a KJ66 micro gas turbine across various off-design simulations.The developed fully coupled method was also compared with the traditional iterative coupled method.Furthermore,experimental data from ground tests were conducted to verify its effectiveness.The convergence history indicated that the proposed fully coupled method exhibited stable convergence,even under far-off-design simulations.The experimental verification demonstrated that the multi-fidelity simulation with the fully coupled method achieved high accuracy in off-design conditions.Further analysis revealed inherent differences in the coupling methods of CFD models between the developed fully coupled and traditional iterative coupled methods.These inherent differences provide valuable insights for reducing errors between the component-level model and CFD models in different coupling methods.The developed fully coupled method,introducing logarithmic transformations,offers more realistic support for the detailed and optimal design of high-fidelity rotating components within the overall performance platform of gas turbines.展开更多
Based on the elastic theory of porous media,embedded discrete fracture model and finite volume method,and considering the micro-seepage mechanism of shale gas,a fully coupled seepage-geomechanical model suitable for f...Based on the elastic theory of porous media,embedded discrete fracture model and finite volume method,and considering the micro-seepage mechanism of shale gas,a fully coupled seepage-geomechanical model suitable for fractured shale gas reservoirs is established,the optimization method of refracturing timing is proposed,and the influencing factors of refracturing timing are analyzed based on the data from shale gas well in Fuling of Sichuan Basin.The results show that due to the depletion of formation pressure,the percentage of the maximum horizontal principal stress reversal area in the total area increases and then decreases with time.The closer the area is to the hydraulic fracture,the shorter the time for the peak of the stress reversal area percentage curve to appear,and the shorter the time for the final zero return(to the initial state).The optimum time of refracturing is affected by matrix permeability,initial stress difference and natural fracture approach angle.The larger the matrix permeability and initial stress difference is,the shorter the time for stress reversal area percentage curve to reach peak and return to the initial state,and the earlier the time to take refracturing measures.The larger the natural fracture approach angle is,the more difficult it is for stress reversal to occur near the fracture,and the earlier the optimum refracturing time is.The more likely the stress reversal occurs at the far end of the artificial fracture,the later the optimal time of refracturing is.Reservoirs with low matrix permeability have a rapid decrease in single well productivity.To ensure economic efficiency,measures such as shut-in or gas injection can be taken to restore the stress,and refracturing can be implemented in advance.展开更多
Numerical simulation is a useful tool in investigating the loading performance of rock bolts.The cable structural elements(cableSELs)in FLAC3D are commonly adopted to simulate rock bolts to solve geotechnical issues.I...Numerical simulation is a useful tool in investigating the loading performance of rock bolts.The cable structural elements(cableSELs)in FLAC3D are commonly adopted to simulate rock bolts to solve geotechnical issues.In this study,the bonding performance of the interface between the rock bolt and the grout material was simulated with a two-stage shearing coupling model.Furthermore,the FISH language was used to incorporate this two-stage shear coupling model into FLAC3D to modify the current cableSELs.Comparison was performed between numerical and experimental results to confirm that the numerical approach can properly simulate the loading performance of rock bolts.Based on the modified cableSELs,the influence of the bolt diameter on the performance of rock bolts and the shear stress propagation along the interface between the bolt and the grout were studied.The simulation results indicated that the load transfer capacity of rock bolts rose with the rock bolt diameter apparently.With the bolt diameter increasing,the performance of the rock bolting system was likely to change from the ductile behaviour to the brittle behaviour.Moreover,after the rock bolt was loaded,the position where the maximum shear stress occurred was variable.Specifically,with the continuous loading,it shifted from the rock bolt loaded end to the other end.展开更多
To keep coal workers away from the hazardous area with frequent accidents such as the roof fall and rib spalling in an underground coalmine,we put forward the solution with robotized self-moving anchor-supporting unit...To keep coal workers away from the hazardous area with frequent accidents such as the roof fall and rib spalling in an underground coalmine,we put forward the solution with robotized self-moving anchor-supporting unit.The existing research shows that the surrounding rock of the roadway has self-stability,and the early or late support is not conducive to the safe and reliable support of the roadway,so there is a problem of support opportunity.In order to study the supporting effect and the optimal supporting time of the above solution,we established the mechanical coupling model of surrounding rock and advance support,and investigated the surrounding rock deformation and advance support pressure distribution under different reserved roof subsidence by using the numerical simulation software FLAC3D.The results show that the deformation of surrounding rock increases and finally tends to a stable level with the increase of pre settlement of roadway roof,and when the pre settlement of roof is between 8-15 mm,the vertical pressure of the top beam of advance support reaches the minimum value,about 0.58 MPa.Based on the above research,we put forward the optimum supporting time in roadway excavation,and summarized the evaluation method based on the mechanical coupling model of surrounding rock-advance support.展开更多
基金Supported by the National Natural Science Foundation of China(51079107,50839004)the Program for New Century Excellent Talents in University (NCET-09-0610)
文摘Understanding the rock mass response to excavation and thermal loading and improving the capability of the numerical models for simulating the progressive failure process of brittle rocks are important for safety assessment and optimization design of nuclear waste repositories.The international cooperative DECOVALEX-2011 project provides a platform for development,validation and comparison of numerical models,in which the sp pillar stability experiment(APSE) was selected as the modeling target for Task B.This paper presents the modeling results of Wuhan University(WHU) team for stages 1 and 2 of Task B by using a coupled thermo-mechanical model within the framework of continuum mechanics.The rock mass response to excavation is modeled with linear elastic,elastoplastic and brittle-plastic models,while the response to heating is modeled with a coupled thermo-elastic model.The capabilities and limitations of the model for representation of the thermo-mechanical responses of the rock pillar are discussed by comparing the modeling results with experimental observations.The results may provide a helpful reference for the stability and safety assessment of the hard granite host rock in China's Beishan preselected area for high-level radioactive waste disposal.
基金supported by the National Natural Science Foundation of China(Grant No.51974173)the Natural Science Foundation of Shandong Province,China(Grant No.ZR2020QD122).
文摘The shear failure of intact rock under thermo-mechanical(TM)coupling conditions is common,such as in enhanced geothermal mining and deep mine construction.Under the effect of a continuous engineering disturbance,shear-formed fractures are prone to secondary instability,posing a severe threat to deep engineering.Although numerous studies regarding three-dimensional(3D)morphologies of fracture surfaces have been conducted,the understanding of shear-formed fractures under TM coupling conditions is limited.In this study,direct shear tests of intact granite under various TM coupling conditions were conducted,followed by 3D laser scanning tests of shear-formed fractures.Test results demonstrated that the peak shear strength of intact granite is positively correlated with the normal stress,whereas it is negatively correlated with the temperature.The internal friction angle and cohesion of intact granite significantly decrease with an increase in the temperature.The anisotropy,roughness value,and height of the asperities on the fracture surfaces are reduced as the normal stress increases,whereas their variation trends are the opposite as the temperature increases.The macroscopic failure mode of intact granite under TM coupling conditions is dominated by mixed tensileeshear and shear failures.As the normal stress increases,intragranular fractures are developed ranging from a local to a global distribution,and the macroscopic failure mode of intact granite changes from mixed tensileeshear to shear failure.Finally,3D morphological characteristics of the asperities on the shear-formed fracture surfaces were analyzed,and a quadrangular pyramid conceptual model representing these asperities was proposed and sufficiently verified.
基金This study is supported by the National Natural Science Foundation of China (Grant No50579046) the Science Foundation of Tianjin Municipal Commission of Science and Technology (Grant No043114711)
文摘This paper, with a finite element method, studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations. With this fully coupled model, the rigid structure is taken as "fictitious" fluid with zero strain rate. Both fluid and structure are described by velocity and pressure. The whole domain, including fluid region and structure region, is modeled by the incompressible Navier-Stokes equations which are discretized with fixed Eulerian mesh. However, to keep the structure' s rigid body shape and behavior, a rigid body constraint is enforced on the "fictitious" fluid domain by use of the Distributed Lagrange Multipher/Fictitious Domain (DLM/ FD) method which is originally introduced to solve particulate flow problems by Glowinski et al. For the verification of the model presented herein, a 2D numerical wave tank is established to simulate small amplitude wave propagations, and then numerical results are compared with analytical solutions. Finally, a 2D example of fluid-structure interaction under wave dynamic forces provides convincing evidences for the method excellent solution quality and fidelity.
基金The research was financially supported by the National Natural Science Foundation of China for Youth(Grant No.51609169)Guangxi Science and Technology Major Project(Grant No.Guike AA17292007)+2 种基金the National Key R&D Program of China(Grant No.2018YFC0310502)National Natural Science Foundation of China(Grant No.51779173)China Scholarship Council(CSC).
文摘Tension Leg Platform(TLP)in deepwater oil and gas field development usually consists of a hull,tendons,and top tension risers(TTRs).To maintain its top tension,each TTR is connected with a tensioner system to the hull.Owing to the complicated configuration of the tensioners,the hull and TTRs form a strong coupled system.Traditionally,some simplified tensioner models are applied to analyze the TLP structures.There is a large discrepancy between their analysis results and the actual mechanism behaviors of a tensioner.It is very necessary to develop a more detailed tensioner model to consider the coupling effects between TLP and TTRs.In the present study,a fully coupled TLP hull-TTR system for hydrodynamic numerical simulation is established.A specific hydraulic pneumatic tensioner is modeled by considering 4 cylinders.The production TTR model is stacked up by specific riser joints.The simulation is also extended to analyze an array of TTRs.Different regular and irregular waves are considered.The behaviors of different cylinders are presented.The results show that it is important to consider the specific configurations of the tensioner and TTRs,which may lead to obviously different response behaviors,compared with those from a simplified model.
基金Project(11272119)supported by the National Natural Science Foundation of China。
文摘This paper developed a statistical damage constitutive model for deep rock by considering the effects of external load and thermal treatment temperature based on the distortion energy.The model parameters were determined through the extremum features of stress−strain curve.Subsequently,the model predictions were compared with experimental results of marble samples.It is found that when the treatment temperature rises,the coupling damage evolution curve shows an S-shape and the slope of ascending branch gradually decreases during the coupling damage evolution process.At a constant temperature,confining pressure can suppress the expansion of micro-fractures.As the confining pressure increases the rock exhibits ductility characteristics,and the shape of coupling damage curve changes from an S-shape into a quasi-parabolic shape.This model can well characterize the influence of high temperature on the mechanical properties of deep rock and its brittleness-ductility transition characteristics under confining pressure.Also,it is suitable for sandstone and granite,especially in predicting the pre-peak stage and peak stress of stress−strain curve under the coupling action of confining pressure and high temperature.The relevant results can provide a reference for further research on the constitutive relationship of rock-like materials and their engineering applications.
基金National Natural Science Foundation of China (51275330,51775361)Natural Science Foundation of Shanxi Province (2014011015-5).
文摘The process of plane strain compression of 06Cr19Ni9NbN steel was carried out in the temperature range of 1000-1200℃ and the reduction ratio range of 10-50%.Combining the finite element numerical simulation,a new constitutive model of thermal deformation was established, which provides the theoretical basis to optimize the plarie strain compression process of the steel.The temperature and grain size at different regions were achieved by experiment and simulation,respectively.According to the results,the mathematical models of stress and temperature during the plane strain compression were established by mathematical analysis method.The new temperature models were established in three regions,respectively,and the stress models took account of the variation of temperature and'st^rain rate.Finally,by comparing the results of calculation,numerical simulation and experiment,the accuracy and validity of these mathematical models were verified.
基金funded by the Science and Technology Innovation Committee Foundation of Shenzhen,Grant No.JCYJ20200109141403840 and Grant No.ZDSYS20220527171405012the National Natural Science Foundation of China (NSFC),Grant No.52106045。
文摘Multi-fidelity simulations incorporate computational fluid dynamics(CFD) models into a thermodynamic model,enabling the simulation of the overall performance of an entire gas turbine with high-fidelity components.Traditional iterative coupled methods rely on characteristic maps,while fully coupled methods directly incorporate high-fidelity simulations.However,fully coupled methods face challenges in simulating rotating components,including weak convergence and complex implementation.To address these challenges,a fully coupled method with logarithmic transformations was developed to directly integrate high-fidelity CFD models of multiple rotating components.The developed fully coupled method was then applied to evaluate the overall performance of a KJ66 micro gas turbine across various off-design simulations.The developed fully coupled method was also compared with the traditional iterative coupled method.Furthermore,experimental data from ground tests were conducted to verify its effectiveness.The convergence history indicated that the proposed fully coupled method exhibited stable convergence,even under far-off-design simulations.The experimental verification demonstrated that the multi-fidelity simulation with the fully coupled method achieved high accuracy in off-design conditions.Further analysis revealed inherent differences in the coupling methods of CFD models between the developed fully coupled and traditional iterative coupled methods.These inherent differences provide valuable insights for reducing errors between the component-level model and CFD models in different coupling methods.The developed fully coupled method,introducing logarithmic transformations,offers more realistic support for the detailed and optimal design of high-fidelity rotating components within the overall performance platform of gas turbines.
基金Supported by National Natural Science Foundation Joint Fund Project(U21B2071)National Natural Science Foundation of China(52174033)National Natural Science Youth Foundation of China(52304041).
文摘Based on the elastic theory of porous media,embedded discrete fracture model and finite volume method,and considering the micro-seepage mechanism of shale gas,a fully coupled seepage-geomechanical model suitable for fractured shale gas reservoirs is established,the optimization method of refracturing timing is proposed,and the influencing factors of refracturing timing are analyzed based on the data from shale gas well in Fuling of Sichuan Basin.The results show that due to the depletion of formation pressure,the percentage of the maximum horizontal principal stress reversal area in the total area increases and then decreases with time.The closer the area is to the hydraulic fracture,the shorter the time for the peak of the stress reversal area percentage curve to appear,and the shorter the time for the final zero return(to the initial state).The optimum time of refracturing is affected by matrix permeability,initial stress difference and natural fracture approach angle.The larger the matrix permeability and initial stress difference is,the shorter the time for stress reversal area percentage curve to reach peak and return to the initial state,and the earlier the time to take refracturing measures.The larger the natural fracture approach angle is,the more difficult it is for stress reversal to occur near the fracture,and the earlier the optimum refracturing time is.The more likely the stress reversal occurs at the far end of the artificial fracture,the later the optimal time of refracturing is.Reservoirs with low matrix permeability have a rapid decrease in single well productivity.To ensure economic efficiency,measures such as shut-in or gas injection can be taken to restore the stress,and refracturing can be implemented in advance.
基金This paper was funded by the following:National Natural Science Foundation of China(51974317,51904302,52034009)Yue Qi Distinguished Scholar Project(800015Z1179,800015Z1138)China University of Mining and Technology(Beijing)and the Fundamental Research Funds for the Central Universities(2020YQNY06).
文摘Numerical simulation is a useful tool in investigating the loading performance of rock bolts.The cable structural elements(cableSELs)in FLAC3D are commonly adopted to simulate rock bolts to solve geotechnical issues.In this study,the bonding performance of the interface between the rock bolt and the grout material was simulated with a two-stage shearing coupling model.Furthermore,the FISH language was used to incorporate this two-stage shear coupling model into FLAC3D to modify the current cableSELs.Comparison was performed between numerical and experimental results to confirm that the numerical approach can properly simulate the loading performance of rock bolts.Based on the modified cableSELs,the influence of the bolt diameter on the performance of rock bolts and the shear stress propagation along the interface between the bolt and the grout were studied.The simulation results indicated that the load transfer capacity of rock bolts rose with the rock bolt diameter apparently.With the bolt diameter increasing,the performance of the rock bolting system was likely to change from the ductile behaviour to the brittle behaviour.Moreover,after the rock bolt was loaded,the position where the maximum shear stress occurred was variable.Specifically,with the continuous loading,it shifted from the rock bolt loaded end to the other end.
基金National Key Basic Research and Development Program Fund project(Grant No.2014CB046306)the Central University Funding Project for Basic Scientific Research Operations(Grant No.2009QJ16)
文摘To keep coal workers away from the hazardous area with frequent accidents such as the roof fall and rib spalling in an underground coalmine,we put forward the solution with robotized self-moving anchor-supporting unit.The existing research shows that the surrounding rock of the roadway has self-stability,and the early or late support is not conducive to the safe and reliable support of the roadway,so there is a problem of support opportunity.In order to study the supporting effect and the optimal supporting time of the above solution,we established the mechanical coupling model of surrounding rock and advance support,and investigated the surrounding rock deformation and advance support pressure distribution under different reserved roof subsidence by using the numerical simulation software FLAC3D.The results show that the deformation of surrounding rock increases and finally tends to a stable level with the increase of pre settlement of roadway roof,and when the pre settlement of roof is between 8-15 mm,the vertical pressure of the top beam of advance support reaches the minimum value,about 0.58 MPa.Based on the above research,we put forward the optimum supporting time in roadway excavation,and summarized the evaluation method based on the mechanical coupling model of surrounding rock-advance support.