Drilling and blasting,characterized by their efficiency,ubiquity,and cost-effectiveness,have emerged as predominant techniques in rock excavation;however,they are accompanied by enormous destructive power.Accurately c...Drilling and blasting,characterized by their efficiency,ubiquity,and cost-effectiveness,have emerged as predominant techniques in rock excavation;however,they are accompanied by enormous destructive power.Accurately controlling the blasting energy and achieving the directional fracture of a rock mass have become common problems in the field.A two-dimensional blasting(2D blasting)technique was proposed that utilizes the characteristic that the tensile strength of a rock mass is significantly lower than its compressive strength.After blasting,only a 2D crack surface is generated along the predetermined direction,eliminating the damage to the reserved rock mass caused by conventional blasting.However,the interior of a natural rock mass is a"black box",and the process of crack propagation is difficult to capture,resulting in an unclear 2D blasting mechanism.To this end,a single-hole polymethyl methacrylate(PMMA)test piece was used to conduct a 2D blasting experiment with the help of a high-speed camera to capture the dynamic crack propagation process and the digital image correlation(DIC)method to analyze the evolution law of surface strain on the test piece.On this basis,a three-dimensional(3D)finite element model was established based on the progressive failure theory to simulate the stress,strain,damage,and displacement evolution process of the model under 2D blasting.The simulation results were consistent with the experimental results.The research results reveal the 2D blasting mechanism and provide theoretical support for the application of 2D blasting technology in the field of rock excavation.展开更多
Swells are critical concerns regarding safety,marine transportation,and coastal engineering construction of coastal countries along the Gulf of Guinea and have been scientific problems due to the lack of systematic th...Swells are critical concerns regarding safety,marine transportation,and coastal engineering construction of coastal countries along the Gulf of Guinea and have been scientific problems due to the lack of systematic theoretical,numerical,and observational research.In this study,a double nesting numerical model was constructed and validated from the Atlantic Ocean to the Gulf of Guinea based on simulating waves nearshore(SWAN)to explore the swell characteristics and source tracing in the Gulf of Guinea in winter and summer seasons from 2020 to 2021.Simulation results reveal that swells are stronger and deflect more to the west in winter than summer,even though they dominate in both seasons in the Gulf of Guinea in the S-SW directional range.Simulated two-dimensional(2D)wave spectral patterns not only clarify wave composition,variation,and propagation properties from the central South Atlantic Ocean to the Gulf of Guinea,but also distinguish swell strength and directional range in winter and summer.The NW wind events induce swells which spread toward the SSE-ESE direction from the North Atlantic Ocean,big wind source generates sustained and stable S-SW swells from the South Atlantic Ocean,and corresponding swell-influenced areas are discussed.The strongest swell event in the Gulf of Guinea during the simulation was used as a case study to trace its source.A strong clockwise wind vortex within the Roaring Forties induced these large swells in the Gulf of Guinea approximately 5.5 days later,and swell propagation formed a regular isoline of peak period distribution from the South Atlantic Ocean to the Gulf of Guinea in the SSW-SW direction.展开更多
0 INTRODUCTION In recent years,modern railways have been actively under construction in the complex mountainous area of Southwest China.However,rockfall poses a significant threat to both construction and operation ph...0 INTRODUCTION In recent years,modern railways have been actively under construction in the complex mountainous area of Southwest China.However,rockfall poses a significant threat to both construction and operation phases of railway projects(Yan et al.,2023;Chen et al.,2022;Fanos and Pradhan,2018).展开更多
The numerical simulation is based on the authors' high-order models with a dissipative term for nonlinear and dispersive wave in water of varying depth. Corresponding finite-difference equations and general condit...The numerical simulation is based on the authors' high-order models with a dissipative term for nonlinear and dispersive wave in water of varying depth. Corresponding finite-difference equations and general conditions for open and fixed natural boundaries with an arbitrary reflection coefficient and phase shift are also given in this paper. The systematical tests of numerical simulation show that the theoretical models, the finite-difference algorithms and the boundary conditions can give good calculation results for the wave propagating in shallow and deep water with an arbitrary slope varying from gentle to steep.展开更多
POM97, an oceanic model, has been used for the first time to the numerical study on the tidal waves of the as Regions around Taiwan. In this paper, we have got the result that the semidiurnal tidal waves of these area...POM97, an oceanic model, has been used for the first time to the numerical study on the tidal waves of the as Regions around Taiwan. In this paper, we have got the result that the semidiurnal tidal waves of these area mainly are the co-operating tides which come from the south of 23'N of the western Pacific. Those semidiurnal tidal waves affecting the Taiwan Straits come respectively from the south and the north entrance of the channel, and the north tidal wave is stronger than the south one. The strongest tidal field is the area from the Meizhou Bay to the Xinhua Bay along the coast of Fujian Province, where the biggest amplitude of the M2 partial tide can reach 240 cm. The strongest tidal cur- rent fields lie in the Penghu watercourse, where the maximum velocity of the M2 partial tide can arrive at 196 m/s. In the horizontal structure of the tidal currets, we have found that there is a stream dot in the north of the channel, besides, there still exist four new ones. As for the vertical structure, it mainly is biassed to the right at the surface, and to the left near the bottom layer.展开更多
Dynamic analysis steps and general flow of fast lagrangian analysis of continua in 3 dimensions(FLAC3D) were discussed. Numerical simulation for influence of excavation and blasting vibration on stability of mined-out...Dynamic analysis steps and general flow of fast lagrangian analysis of continua in 3 dimensions(FLAC3D) were discussed. Numerical simulation for influence of excavation and blasting vibration on stability of mined-out area was carried out with FLAC3D. The whole analytical process was divided into two steps, including the static analysis and the dynamic analysis which were used to simulate the influence of excavation process and blasting vibration respectively. The results show that the shape of right upper boundary is extremely irregular after excavation, and stress concentration occurs at many places and higher tensile stress appears. The maximum tensile stress is higher than the tensile strength of rock mass, and surrounding rock of right roof will be damaged with tension fracture. The maximum displacement of surrounding rock is 4.75 mm after excavation. However, the maximum displacement increases to 5.47 mm after the blasting dynamic load is applied. And the covering area of plastic zones expands obviously, especially at the foot of right upper slope. The analytical results are in basic accordance with the observed results on the whole. Damage and disturbance on surrounding rock to some degree are caused by excavation, while blasting dynamic load increases the possibility of occurrence of dynamic instability and destruction further. So the effective supporting and vibration reducing measures should be taken during mining.展开更多
The developing processes of stress and deformation fields of a protected layer after mining an upper-protective layer with a bow pseudo-incline technique were simulated to locate the protection region. The pressure re...The developing processes of stress and deformation fields of a protected layer after mining an upper-protective layer with a bow pseudo-incline technique were simulated to locate the protection region. The pressure relief of the protected layer was analyzed after mining the upper-protective layer. The pressure relief angle along the strike and incline were located according to the roles of protection of the deformation and stress pressure-relief of the protective layer after mining. This results show that the upper-protective layer with the bow pseudo-incline technique have an upper and downside pressure relief angle of 85 and 68 degrees respectively; the distribution of strike pressure relief angles along the pseudo-incline working face is uneven and their values range from 38.3 to 51 degrees. The pressure relief angle of the inclined middle location was the largest. The distribution of the protection region of the upper-protective layer with the bow pseudo-incline teelmique located by practical tests and numerical simulation is essentially consistent, compared with the results obtained by these methods.展开更多
Design of forming dies and whole process of simulation of cold rolling involutes spline can be realized by using of CAD software of PRO-E and CAE software of DEFORM-3D. Software DEFORM-3D provides an automatic and opt...Design of forming dies and whole process of simulation of cold rolling involutes spline can be realized by using of CAD software of PRO-E and CAE software of DEFORM-3D. Software DEFORM-3D provides an automatic and optimized remeshing function, especially for the large deformation. In order to use this function sufficiently, simulation of cold rolling involutes spline can be implemented indirectly. The relationship between die and workpiece, forming force and characteristic of deformation in the forming process of cold rolling involutes spline are analyzed and researched. Meanwhile, reliable proofs for the design of dies and deforming equipment are provided.展开更多
Based on specific geology and mining conditions of certain coal working face in China, a series of numerical models under different unconsolidated layers thickness were respectively established by employing FLAC3D. Th...Based on specific geology and mining conditions of certain coal working face in China, a series of numerical models under different unconsolidated layers thickness were respectively established by employing FLAC3D. The relationship between the unconsolidated layers thickness and surface movement laws was studied. Maximum surface subsidence, Maximum horizontal displacement and surface subsidence degree were obtained. Contours of surface subsidence/horizontal displacement and curves were drawn. Some laws of surface subsidence/horizontal displacement were analyzed. The role of the unconsolidated layers in surface subsidence was revealed. It is significant to predict surface subsidence of thick unconsolidated layers for coal mine and take effective measures to control surface subsidence.展开更多
In order to verify the validity of finite element numerical simulation method for asphalt mixture, which consists of aggregates, mastic （where mastic is a kind of fine mixture composed of asphalt binder mixed with fi...In order to verify the validity of finite element numerical simulation method for asphalt mixture, which consists of aggregates, mastic （where mastic is a kind of fine mixture composed of asphalt binder mixed with fines and fine aggregates） and air voids, based on three-dimensional （3D） heterogeneous specimen, X-ray computerized tomography （X-ray CT） was used to scan the asphalt specimens to obtain the real internal microstrnctures of asphalt mixture. CT images were reconstructed to build up 3D digital specimen, and the viscoelastic properties of mastic were described with Burgers model The uniaxial creep numerical simulations of three different levels of aggregate gradation were conducted. The simulation results agree well with the experimental results. The numerical simulation of asphalt mixture incorporated with real 3D microstructure based on finite element method is a promising application to conduct research of asphalt concrete. Additionally, this method can increase the mechanistic understanding of global viscoelastic properties of asphalt mixtures by linking the real 3D microstructure.展开更多
In order to simulate field distribution rules,mathematical models for 3-D air flows and gas transfer in the goaf of the coal mining are established,based on theories of permeability and dynamic dispersion through poro...In order to simulate field distribution rules,mathematical models for 3-D air flows and gas transfer in the goaf of the coal mining are established,based on theories of permeability and dynamic dispersion through porous media. A gas dispersion equation in a 3-D field is calculated by use of numerical method on a weighted upstream multi-element balance. Based on data of an example with a U type ventilation mode,surface charts of air pressure distribution and gas concentration are drawn by Graphtool software. Finally,a comparison between actually measured results in the model test and the numerical simulation results is made to proves the numerical implementation feasible.展开更多
A numerical method is presented that simulates 3D explosive field problems. A code MMIC3D using this method can be used to simulate the propagation and reflected effects of all kinds of rigid boundaries to shock waves...A numerical method is presented that simulates 3D explosive field problems. A code MMIC3D using this method can be used to simulate the propagation and reflected effects of all kinds of rigid boundaries to shock waves produced by an explosive source. These numerical results indicate that the code MMIC3D has the ability in computing cases such as 3D shock waves produced by air explosion, vortex region of the shock wave, the Mach wave, and reflected waves behind rigid boundaries.展开更多
The finite element simulation software SYSWELD is used to numerically simulate the temperature field,residual stress field,and welding deformation of Q690D thick plate multi-layer and multi-pass welding under differen...The finite element simulation software SYSWELD is used to numerically simulate the temperature field,residual stress field,and welding deformation of Q690D thick plate multi-layer and multi-pass welding under different welding heat input and groove angles.The simulation results show that as the welding heat input increases,the peak temperature during the welding process is higher,and the residual stress increases,they are all between 330–340 MPa,and the residual stress is concentrated in the area near the weld.The hole-drilling method is used to measure the actual welding residual stress,and the measured data is in good agreement with the simulated value.The type of post-welding deformation is angular deformation,and as the welding heat input increases,the maximum deformation also increases.It shows smaller residual stress and deformation when the groove angle is 40°under the same heat input.In engineering applications,under the premise of guaranteeing welding quality,smaller heat input and 40°groove angle should be used.展开更多
Debris flows are rapid mass movements with a mixture of rock,soil and water.High-intensity rainfall events have triggered multiple debris flows around the globe,making it an important concern from the disaster managem...Debris flows are rapid mass movements with a mixture of rock,soil and water.High-intensity rainfall events have triggered multiple debris flows around the globe,making it an important concern from the disaster management perspective.This study presents a numerical model called debris flow simulation 2D(DFS 2D)and applicability of the proposed model is investigated through the values of the model parameters used for the reproduction of an occurred debris flow at Yindongzi gully in China on 13 August 2010.The model can be used to simulate debris flows using three different rheologies and has a userfriendly interface for providing the inputs.Using DFS 2D,flow parameters can be estimated with respect to space and time.The values of the flow resistance parameters of model,dry-Coulomb and turbulent friction,were calibrated through the back analysis and the values obtained are 0.1 and 1000 m/s^(2),respectively.Two new methods of calibration are proposed in this study,considering the crosssectional area of flow and topographical changes induced by the debris flow.The proposed methods of calibration provide an effective solution to the cumulative errors induced by coarse-resolution digital elevation models(DEMs)in numerical modelling of debris flows.The statistical indices such as Willmott's index of agreement,mean-absolute-error,and normalized-root-mean-square-error of the calibrated model are 0.5,1.02 and 1.44,respectively.The comparison between simulated and observed values of topographic changes indicates that DFS 2D provides satisfactory results and can be used for dynamic modelling of debris flows.展开更多
In this paper,various aspects of the 2D and 3D nonlinear liquid sloshing problems in vertically excited containers have been studied numerically along with the help of a modified-transformation.Based on this new numer...In this paper,various aspects of the 2D and 3D nonlinear liquid sloshing problems in vertically excited containers have been studied numerically along with the help of a modified-transformation.Based on this new numerical algorithm,a numerical study on a regularly and randomly excited container in vertical direction was conducted utilizing four different cases： The first case was performed utilizing a 2D container with regular excitations.The next case examined a regularly excited 3D container with two different initial conditions for the liquid free surface,and finally,3D container with random excitation in the vertical direction.A grid independence study was performed along with a series of validation tests.An iteration error estimation method was used to stop the iterative solver（used for solving the discretized governing equations in the computational domain） upon reaching steady state of results at each time step.In the present case,this method was found to produce quite accurate results and to be more time efficient as compared to other conventional stopping procedures for iterative solvers.The results were validated with benchmark results.The wave elevation time history,phase plane diagram and surface plots represent the wave nonlinearity during its motion.展开更多
Based on the 2D horizontal plane numerical model, a quasi-3D numerical model is established for coastal regions of shallow water. The characteristics of this model are that the velocity profiles;can be obtained at the...Based on the 2D horizontal plane numerical model, a quasi-3D numerical model is established for coastal regions of shallow water. The characteristics of this model are that the velocity profiles;can be obtained at the same time when the equations of the value of difference between the horizontal current velocity and its depth-averaged velocity in the vertical direction are solved and the results obtained are consistent with the results of the 2D, model. The circulating flow in the rectangular area induced by wind is simulated and applied to the tidal flow field of the radial sandbanks in the South Yellow Sea. The computational results from this quasi-3D model are in good agreement with analytical results and observed data. The solution of the finite difference equations has been found to be stable, and the model is simple, effective and practical.展开更多
We report our numerical simulation on the dynamic interference photoelectron spectra for a one-dimensional （1D） He model exposed to intense ultrashort extreme ultraviolet （XUV） laser pulses. The results demonstrat...We report our numerical simulation on the dynamic interference photoelectron spectra for a one-dimensional （1D） He model exposed to intense ultrashort extreme ultraviolet （XUV） laser pulses. The results demonstrate an unambiguous interference feature in the photoelectron spectra, and the interference is unveiled to originate from the dynamic Stark effect. The interference photoelectron spectra are prompted for intense sub-femtosecond XUV laser pulses in double ionization. The stationary phase picture is corroborated qualitatively in the two-electron system. The ability of probing the dynamic Stark effect by the photoelectron spectra in a pragmatic experiment of single-photon double ionization of He may shed light on further investigation on multi-electron atoms and molecules.展开更多
According to the influence of the combination of short-distance coal seam group on mining roadway, using numerical simulation software FLAG2D to draw the abutment pressure distribution ahead the working face and the a...According to the influence of the combination of short-distance coal seam group on mining roadway, using numerical simulation software FLAG2D to draw the abutment pressure distribution ahead the working face and the area of influence in fully-mechanized mining conditions, the variation rules of surrounding rock supporting pressure of floor roadway and the deformation rules were summarized. GYS-300 anchor dynamometer was used to measure the roadway surface displacement, and the conclusions of numerical simulation were verified.展开更多
In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform ...In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform cubic grids, the approach used to establish stable formulas with 2M-order accuracy is discussed in detail, with M being a positive integer, and is illustrated by establishing second order （M=1） recursion formulas. The theoretical results presented in this paper are demonstrated through numerical testing.展开更多
The riverbank soil is a natural purifying agent for the polluted river water(Riverbank filtration, RBF). This is of great importance to groundwater safety along the riverbank. This paper examines the migration and tra...The riverbank soil is a natural purifying agent for the polluted river water(Riverbank filtration, RBF). This is of great importance to groundwater safety along the riverbank. This paper examines the migration and transformation rules of ammonia-nitrogen in three typical types of sand soil using the indoor leaching experiment of soil column, and then makes comparison with the indoor experiment results in combination with the numerical simulation method. The experiment process shows that the change in ammonia-nitrogen concentration goes through three stages including "removal-water saturation-saturation". As the contents of clay particles in soil sample increase, the removal of ammonia-nitrogen from soil sample will take more time and gain higher ratio. During the removal period, the removal ratio of Column 1, Column 2 and Column 3 averages 68.8%(1-12 d), 74.6%(1-22 d) and 91.1%(1-26 d). The ammonia-nitrogen removal ratio shows no noticeable change as the depth of soil columns varies. But it is found that the ammonia-nitrogen removal ratio is the least of the whole experiment when the soil columns are at the depth of 15 cm. It can be preliminary inferred that the natural purifying performance of soil along the river for ammonia-nitrogen in river water mainly depends on the proportion of fine particles in soil. HYDRUS-1D model is used to simulate this experiment process, analyze the change of the bottom observation holes by time and depth in three columns(the tenth day), and make comparison with the experiment result. The coefficients of determination for fitting curves of Column 1, Column 2 and Column 3 are 0.953, 0.909, 0.882 and 0.955, 0.740, 0.980 separately. Besides, this paper examines the contribution of absorption, mineralization and nitrification in the simulation process. In the early removal stage, mineralization plays a dominant role and the maximum contribution rate of mineralization is 99%. As time goes by, absorption starts to function and gradually assumes a dominant position. In the middle and late removal stage, nitrification in Column 1 and Column 2 makes more contribution than mineralization. So the experiment result of the ammonia-nitrogen concentration is 0.6% and 2.4% lower than that in effluent and the maximum contribution ratio of nitrification is -4.53% and -5.10% respectively when only the function of absorption is considered. The mineralization in Column 1 and Column 2 in the middle and late removal stage still plays a more important role than nitrification. So the experiment result is 1.4% higher than that in effluent and the maximum contribution ratio of nitrification is -2.51% when only the function of absorption is considered. Therefore, absorption, mineralization and nitrification make different contributions during different part of the stage. This means that the natural purifying performance of soil along the river for ammonia-nitrogen in river water not only depends on the proportion of fine particles in soil, but depends on the mineralization and nitrification environment. This can offer some insights into the protection and recovery of groundwater along the riverbank.展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.52404155 and 52304111)State Key Laboratory for Geomechanics and Deep Underground Engineering,China University of Mining&Technology,Beijing(Grant No.XD2024006).
文摘Drilling and blasting,characterized by their efficiency,ubiquity,and cost-effectiveness,have emerged as predominant techniques in rock excavation;however,they are accompanied by enormous destructive power.Accurately controlling the blasting energy and achieving the directional fracture of a rock mass have become common problems in the field.A two-dimensional blasting(2D blasting)technique was proposed that utilizes the characteristic that the tensile strength of a rock mass is significantly lower than its compressive strength.After blasting,only a 2D crack surface is generated along the predetermined direction,eliminating the damage to the reserved rock mass caused by conventional blasting.However,the interior of a natural rock mass is a"black box",and the process of crack propagation is difficult to capture,resulting in an unclear 2D blasting mechanism.To this end,a single-hole polymethyl methacrylate(PMMA)test piece was used to conduct a 2D blasting experiment with the help of a high-speed camera to capture the dynamic crack propagation process and the digital image correlation(DIC)method to analyze the evolution law of surface strain on the test piece.On this basis,a three-dimensional(3D)finite element model was established based on the progressive failure theory to simulate the stress,strain,damage,and displacement evolution process of the model under 2D blasting.The simulation results were consistent with the experimental results.The research results reveal the 2D blasting mechanism and provide theoretical support for the application of 2D blasting technology in the field of rock excavation.
基金The National Key R&D Program of China under contract No.2023YFE0126300the National Natural Science Foundation of China under contract Nos 42066002 and U20A2099.
文摘Swells are critical concerns regarding safety,marine transportation,and coastal engineering construction of coastal countries along the Gulf of Guinea and have been scientific problems due to the lack of systematic theoretical,numerical,and observational research.In this study,a double nesting numerical model was constructed and validated from the Atlantic Ocean to the Gulf of Guinea based on simulating waves nearshore(SWAN)to explore the swell characteristics and source tracing in the Gulf of Guinea in winter and summer seasons from 2020 to 2021.Simulation results reveal that swells are stronger and deflect more to the west in winter than summer,even though they dominate in both seasons in the Gulf of Guinea in the S-SW directional range.Simulated two-dimensional(2D)wave spectral patterns not only clarify wave composition,variation,and propagation properties from the central South Atlantic Ocean to the Gulf of Guinea,but also distinguish swell strength and directional range in winter and summer.The NW wind events induce swells which spread toward the SSE-ESE direction from the North Atlantic Ocean,big wind source generates sustained and stable S-SW swells from the South Atlantic Ocean,and corresponding swell-influenced areas are discussed.The strongest swell event in the Gulf of Guinea during the simulation was used as a case study to trace its source.A strong clockwise wind vortex within the Roaring Forties induced these large swells in the Gulf of Guinea approximately 5.5 days later,and swell propagation formed a regular isoline of peak period distribution from the South Atlantic Ocean to the Gulf of Guinea in the SSW-SW direction.
基金supported by the Open Research Fund of Key Laboratory of Geological Hazards on Three Gorges Reservoir Area(China Three Gorges University),Ministry of Education(No.2022KDZ03)the Science and Technology Projects of Yunnan Provincial Science and Technology Department(No.202401AT070328)+1 种基金the Young talents project of“Xingdian Talent Support Program”in Yunnan Province(No.YNWR-QNBJ-2020-019)the Fund Project of China Academy of Railway Sciences Co.,Ltd.(No.2021YJ178)。
文摘0 INTRODUCTION In recent years,modern railways have been actively under construction in the complex mountainous area of Southwest China.However,rockfall poses a significant threat to both construction and operation phases of railway projects(Yan et al.,2023;Chen et al.,2022;Fanos and Pradhan,2018).
文摘The numerical simulation is based on the authors' high-order models with a dissipative term for nonlinear and dispersive wave in water of varying depth. Corresponding finite-difference equations and general conditions for open and fixed natural boundaries with an arbitrary reflection coefficient and phase shift are also given in this paper. The systematical tests of numerical simulation show that the theoretical models, the finite-difference algorithms and the boundary conditions can give good calculation results for the wave propagating in shallow and deep water with an arbitrary slope varying from gentle to steep.
文摘POM97, an oceanic model, has been used for the first time to the numerical study on the tidal waves of the as Regions around Taiwan. In this paper, we have got the result that the semidiurnal tidal waves of these area mainly are the co-operating tides which come from the south of 23'N of the western Pacific. Those semidiurnal tidal waves affecting the Taiwan Straits come respectively from the south and the north entrance of the channel, and the north tidal wave is stronger than the south one. The strongest tidal field is the area from the Meizhou Bay to the Xinhua Bay along the coast of Fujian Province, where the biggest amplitude of the M2 partial tide can reach 240 cm. The strongest tidal cur- rent fields lie in the Penghu watercourse, where the maximum velocity of the M2 partial tide can arrive at 196 m/s. In the horizontal structure of the tidal currets, we have found that there is a stream dot in the north of the channel, besides, there still exist four new ones. As for the vertical structure, it mainly is biassed to the right at the surface, and to the left near the bottom layer.
基金Project (50490272) supported by the National Natural Science Foundation of China project(NCET-05-0687) supportedby Programfor New Century Excellent Talents project (040109) supported bythe Doctor Degree Paper Innovation Engineering of CentralSouth University
文摘Dynamic analysis steps and general flow of fast lagrangian analysis of continua in 3 dimensions(FLAC3D) were discussed. Numerical simulation for influence of excavation and blasting vibration on stability of mined-out area was carried out with FLAC3D. The whole analytical process was divided into two steps, including the static analysis and the dynamic analysis which were used to simulate the influence of excavation process and blasting vibration respectively. The results show that the shape of right upper boundary is extremely irregular after excavation, and stress concentration occurs at many places and higher tensile stress appears. The maximum tensile stress is higher than the tensile strength of rock mass, and surrounding rock of right roof will be damaged with tension fracture. The maximum displacement of surrounding rock is 4.75 mm after excavation. However, the maximum displacement increases to 5.47 mm after the blasting dynamic load is applied. And the covering area of plastic zones expands obviously, especially at the foot of right upper slope. The analytical results are in basic accordance with the observed results on the whole. Damage and disturbance on surrounding rock to some degree are caused by excavation, while blasting dynamic load increases the possibility of occurrence of dynamic instability and destruction further. So the effective supporting and vibration reducing measures should be taken during mining.
基金Projects PLN0610 supported by the Open Fund of State Key Lab of Oil and Gas Reservoir Geology and Exploitation (Southwest Petroleum University)HKLGF200706 by the Opening Project of Henan Key Laboratory of Coal Mine Methane and Fire Prevention+3 种基金50334060, 50474025 and 50774106 by the National Natural Science Foundation of China2005CB221502 by the National Basic Research Program of China50621403 by the Natural Science Innova-tion Group Foundation of ChinaCSTC, 2006BB7147, 2006AA7002 by the Natural Science Foundation of Chongqing
文摘The developing processes of stress and deformation fields of a protected layer after mining an upper-protective layer with a bow pseudo-incline technique were simulated to locate the protection region. The pressure relief of the protected layer was analyzed after mining the upper-protective layer. The pressure relief angle along the strike and incline were located according to the roles of protection of the deformation and stress pressure-relief of the protective layer after mining. This results show that the upper-protective layer with the bow pseudo-incline technique have an upper and downside pressure relief angle of 85 and 68 degrees respectively; the distribution of strike pressure relief angles along the pseudo-incline working face is uneven and their values range from 38.3 to 51 degrees. The pressure relief angle of the inclined middle location was the largest. The distribution of the protection region of the upper-protective layer with the bow pseudo-incline teelmique located by practical tests and numerical simulation is essentially consistent, compared with the results obtained by these methods.
文摘Design of forming dies and whole process of simulation of cold rolling involutes spline can be realized by using of CAD software of PRO-E and CAE software of DEFORM-3D. Software DEFORM-3D provides an automatic and optimized remeshing function, especially for the large deformation. In order to use this function sufficiently, simulation of cold rolling involutes spline can be implemented indirectly. The relationship between die and workpiece, forming force and characteristic of deformation in the forming process of cold rolling involutes spline are analyzed and researched. Meanwhile, reliable proofs for the design of dies and deforming equipment are provided.
基金Project(2007BAK28B03)supported by the National Eleventh-Five Year Research Program of ChinaProject(2010YD05)supported by the Fundamental Research Funds for the Central UniversitiesProject(200911036)supported by the Ministry of Land and Resources Research Special
文摘Based on specific geology and mining conditions of certain coal working face in China, a series of numerical models under different unconsolidated layers thickness were respectively established by employing FLAC3D. The relationship between the unconsolidated layers thickness and surface movement laws was studied. Maximum surface subsidence, Maximum horizontal displacement and surface subsidence degree were obtained. Contours of surface subsidence/horizontal displacement and curves were drawn. Some laws of surface subsidence/horizontal displacement were analyzed. The role of the unconsolidated layers in surface subsidence was revealed. It is significant to predict surface subsidence of thick unconsolidated layers for coal mine and take effective measures to control surface subsidence.
基金Project(51038004) supported by the National Natural Science Foundation of China
文摘In order to verify the validity of finite element numerical simulation method for asphalt mixture, which consists of aggregates, mastic （where mastic is a kind of fine mixture composed of asphalt binder mixed with fines and fine aggregates） and air voids, based on three-dimensional （3D） heterogeneous specimen, X-ray computerized tomography （X-ray CT） was used to scan the asphalt specimens to obtain the real internal microstrnctures of asphalt mixture. CT images were reconstructed to build up 3D digital specimen, and the viscoelastic properties of mastic were described with Burgers model The uniaxial creep numerical simulations of three different levels of aggregate gradation were conducted. The simulation results agree well with the experimental results. The numerical simulation of asphalt mixture incorporated with real 3D microstructure based on finite element method is a promising application to conduct research of asphalt concrete. Additionally, this method can increase the mechanistic understanding of global viscoelastic properties of asphalt mixtures by linking the real 3D microstructure.
基金Projects 50534090 and 50674090 supported by the National Natural Science Foundation of China2005CB221503 by the National Key Basic ResearchDevelopment Program (973 Program)
文摘In order to simulate field distribution rules,mathematical models for 3-D air flows and gas transfer in the goaf of the coal mining are established,based on theories of permeability and dynamic dispersion through porous media. A gas dispersion equation in a 3-D field is calculated by use of numerical method on a weighted upstream multi-element balance. Based on data of an example with a U type ventilation mode,surface charts of air pressure distribution and gas concentration are drawn by Graphtool software. Finally,a comparison between actually measured results in the model test and the numerical simulation results is made to proves the numerical implementation feasible.
文摘A numerical method is presented that simulates 3D explosive field problems. A code MMIC3D using this method can be used to simulate the propagation and reflected effects of all kinds of rigid boundaries to shock waves produced by an explosive source. These numerical results indicate that the code MMIC3D has the ability in computing cases such as 3D shock waves produced by air explosion, vortex region of the shock wave, the Mach wave, and reflected waves behind rigid boundaries.
文摘The finite element simulation software SYSWELD is used to numerically simulate the temperature field,residual stress field,and welding deformation of Q690D thick plate multi-layer and multi-pass welding under different welding heat input and groove angles.The simulation results show that as the welding heat input increases,the peak temperature during the welding process is higher,and the residual stress increases,they are all between 330–340 MPa,and the residual stress is concentrated in the area near the weld.The hole-drilling method is used to measure the actual welding residual stress,and the measured data is in good agreement with the simulated value.The type of post-welding deformation is angular deformation,and as the welding heat input increases,the maximum deformation also increases.It shows smaller residual stress and deformation when the groove angle is 40°under the same heat input.In engineering applications,under the premise of guaranteeing welding quality,smaller heat input and 40°groove angle should be used.
基金financially supported by Department of Space,India(Grant No.ISRO/RES/4/663/18-19)。
文摘Debris flows are rapid mass movements with a mixture of rock,soil and water.High-intensity rainfall events have triggered multiple debris flows around the globe,making it an important concern from the disaster management perspective.This study presents a numerical model called debris flow simulation 2D(DFS 2D)and applicability of the proposed model is investigated through the values of the model parameters used for the reproduction of an occurred debris flow at Yindongzi gully in China on 13 August 2010.The model can be used to simulate debris flows using three different rheologies and has a userfriendly interface for providing the inputs.Using DFS 2D,flow parameters can be estimated with respect to space and time.The values of the flow resistance parameters of model,dry-Coulomb and turbulent friction,were calibrated through the back analysis and the values obtained are 0.1 and 1000 m/s^(2),respectively.Two new methods of calibration are proposed in this study,considering the crosssectional area of flow and topographical changes induced by the debris flow.The proposed methods of calibration provide an effective solution to the cumulative errors induced by coarse-resolution digital elevation models(DEMs)in numerical modelling of debris flows.The statistical indices such as Willmott's index of agreement,mean-absolute-error,and normalized-root-mean-square-error of the calibrated model are 0.5,1.02 and 1.44,respectively.The comparison between simulated and observed values of topographic changes indicates that DFS 2D provides satisfactory results and can be used for dynamic modelling of debris flows.
文摘In this paper,various aspects of the 2D and 3D nonlinear liquid sloshing problems in vertically excited containers have been studied numerically along with the help of a modified-transformation.Based on this new numerical algorithm,a numerical study on a regularly and randomly excited container in vertical direction was conducted utilizing four different cases： The first case was performed utilizing a 2D container with regular excitations.The next case examined a regularly excited 3D container with two different initial conditions for the liquid free surface,and finally,3D container with random excitation in the vertical direction.A grid independence study was performed along with a series of validation tests.An iteration error estimation method was used to stop the iterative solver（used for solving the discretized governing equations in the computational domain） upon reaching steady state of results at each time step.In the present case,this method was found to produce quite accurate results and to be more time efficient as compared to other conventional stopping procedures for iterative solvers.The results were validated with benchmark results.The wave elevation time history,phase plane diagram and surface plots represent the wave nonlinearity during its motion.
基金National Natural Science Foundation of China(Grant No.49236120)
文摘Based on the 2D horizontal plane numerical model, a quasi-3D numerical model is established for coastal regions of shallow water. The characteristics of this model are that the velocity profiles;can be obtained at the same time when the equations of the value of difference between the horizontal current velocity and its depth-averaged velocity in the vertical direction are solved and the results obtained are consistent with the results of the 2D, model. The circulating flow in the rectangular area induced by wind is simulated and applied to the tidal flow field of the radial sandbanks in the South Yellow Sea. The computational results from this quasi-3D model are in good agreement with analytical results and observed data. The solution of the finite difference equations has been found to be stable, and the model is simple, effective and practical.
基金Supported by the National Natural Science Foundation of China under Grant Nos 61178028,11674243 and 11674242the National Basic Research Program of China under Grant No 2015CB755403
文摘We report our numerical simulation on the dynamic interference photoelectron spectra for a one-dimensional （1D） He model exposed to intense ultrashort extreme ultraviolet （XUV） laser pulses. The results demonstrate an unambiguous interference feature in the photoelectron spectra, and the interference is unveiled to originate from the dynamic Stark effect. The interference photoelectron spectra are prompted for intense sub-femtosecond XUV laser pulses in double ionization. The stationary phase picture is corroborated qualitatively in the two-electron system. The ability of probing the dynamic Stark effect by the photoelectron spectra in a pragmatic experiment of single-photon double ionization of He may shed light on further investigation on multi-electron atoms and molecules.
文摘According to the influence of the combination of short-distance coal seam group on mining roadway, using numerical simulation software FLAG2D to draw the abutment pressure distribution ahead the working face and the area of influence in fully-mechanized mining conditions, the variation rules of surrounding rock supporting pressure of floor roadway and the deformation rules were summarized. GYS-300 anchor dynamometer was used to measure the roadway surface displacement, and the conclusions of numerical simulation were verified.
基金China Postdoctoral Science Foundation Under Grant No.20100480321National Basic Research Program of China Under Grant No. 2007CB714200
文摘In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform cubic grids, the approach used to establish stable formulas with 2M-order accuracy is discussed in detail, with M being a positive integer, and is illustrated by establishing second order （M=1） recursion formulas. The theoretical results presented in this paper are demonstrated through numerical testing.
基金supported by Special Scientific Research Expenditure for Public Charity Industry of Ministry of Water Resources(No.201501008)Institute of Resources and Environment of North China University of Water Resources and Electric Power
文摘The riverbank soil is a natural purifying agent for the polluted river water(Riverbank filtration, RBF). This is of great importance to groundwater safety along the riverbank. This paper examines the migration and transformation rules of ammonia-nitrogen in three typical types of sand soil using the indoor leaching experiment of soil column, and then makes comparison with the indoor experiment results in combination with the numerical simulation method. The experiment process shows that the change in ammonia-nitrogen concentration goes through three stages including "removal-water saturation-saturation". As the contents of clay particles in soil sample increase, the removal of ammonia-nitrogen from soil sample will take more time and gain higher ratio. During the removal period, the removal ratio of Column 1, Column 2 and Column 3 averages 68.8%(1-12 d), 74.6%(1-22 d) and 91.1%(1-26 d). The ammonia-nitrogen removal ratio shows no noticeable change as the depth of soil columns varies. But it is found that the ammonia-nitrogen removal ratio is the least of the whole experiment when the soil columns are at the depth of 15 cm. It can be preliminary inferred that the natural purifying performance of soil along the river for ammonia-nitrogen in river water mainly depends on the proportion of fine particles in soil. HYDRUS-1D model is used to simulate this experiment process, analyze the change of the bottom observation holes by time and depth in three columns(the tenth day), and make comparison with the experiment result. The coefficients of determination for fitting curves of Column 1, Column 2 and Column 3 are 0.953, 0.909, 0.882 and 0.955, 0.740, 0.980 separately. Besides, this paper examines the contribution of absorption, mineralization and nitrification in the simulation process. In the early removal stage, mineralization plays a dominant role and the maximum contribution rate of mineralization is 99%. As time goes by, absorption starts to function and gradually assumes a dominant position. In the middle and late removal stage, nitrification in Column 1 and Column 2 makes more contribution than mineralization. So the experiment result of the ammonia-nitrogen concentration is 0.6% and 2.4% lower than that in effluent and the maximum contribution ratio of nitrification is -4.53% and -5.10% respectively when only the function of absorption is considered. The mineralization in Column 1 and Column 2 in the middle and late removal stage still plays a more important role than nitrification. So the experiment result is 1.4% higher than that in effluent and the maximum contribution ratio of nitrification is -2.51% when only the function of absorption is considered. Therefore, absorption, mineralization and nitrification make different contributions during different part of the stage. This means that the natural purifying performance of soil along the river for ammonia-nitrogen in river water not only depends on the proportion of fine particles in soil, but depends on the mineralization and nitrification environment. This can offer some insights into the protection and recovery of groundwater along the riverbank.
