To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based sim...To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based simulation(NNS)method with higher accuracy and better efficiency was proposed.The NNS method consisted of three main steps.First,the parameters of blast loads,including the peak pressures and impulses of cylindrical charges with different aspect ratios(L/D)at different stand-off distances and incident angles were obtained by two-dimensional numerical simulations.Subsequently,incident shape factors of cylindrical charges with arbitrary aspect ratios were predicted by a neural network.Finally,reflected shape factors were derived and implemented into the subroutine of the ABAQUS code to modify the CONWEP model,including modifications of impulse and overpressure.The reliability of the proposed NNS method was verified by related experimental results.Remarkable accuracy improvement was acquired by the proposed NNS method compared with the unmodified CONWEP model.Moreover,huge efficiency superiority was obtained by the proposed NNS method compared with the CEL method.The proposed NNS method showed good accuracy when the scaled distance was greater than 0.2 m/kg^(1/3).It should be noted that there is no need to generate a new dataset again since the blast loads satisfy the similarity law,and the proposed NNS method can be directly used to simulate the blast loads generated by different cylindrical charges.The proposed NNS method with high efficiency and accuracy can be used as an effective method to analyze the dynamic response of structures under blast loads,and it has significant application prospects in designing protective structures.展开更多
In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step ...In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step for the dynamic response simulation of rock mass in a high in-situ stress field.In this paper,stress initialization methods,including their principles and operating procedures for reproducing steady in-situ stress state in LS-DYNA,are first introduced.Then the most popular four methods,i.e.,explicit dynamic relaxation(DR)method,implicit-explicit sequence method,Dynain file method and quasi-static method,are exemplified through a case analysis by using the RHT and plastic hardening rock material models to simulate rock blasting under in-situ stress condition.Based on the simulations,it is concluded that the stress initialization results obtained by implicit-explicit sequence method and dynain file method are closely related to the rock material model,and the explicit DR method has an obvious advantage in solution time when compared to other methods.Besides that,it is recommended to adopt two separate analyses for the whole numerical simulation of rock mass under the combined action of in-situ stress and dynamic disturbance.展开更多
Outwash deposit is a unique type of geological materials, and its features such as heterogeneity, discontinuity and nonlinearity determine the complexity of mechanical characteristics and failure mechanism. In this wo...Outwash deposit is a unique type of geological materials, and its features such as heterogeneity, discontinuity and nonlinearity determine the complexity of mechanical characteristics and failure mechanism. In this work, random meso-structure of outwash deposits was constructed by the technique of computer random simulation based on characteristics of its meso-structure in the statistical sense and some simplifications, and a series of large direct shear tests on numerical samples of outwash deposits with stone contents of 15%, 30%, 45% and 60% were conducted using the discrete element method to further investigate its mechanical characteristics and failure mechanism under external load. The results show that the deformation characteristics and shear strength of outwash deposits are to some extent improved with the increase of stone content, and the shear stress–shear displacement curves of outwash deposits show great differences at the post-peak stage due to the random spatial distribution and content of stones. From the mesoscopic view, normal directions of contacts between "soil" and "stone" particles undergo apparent deflection as the shear displacement continues during the shearing process, accompanying redistribution of the magnitude of contact forces during the shearing process. For outwash deposits, the shear zone formed after shear failure is an irregular stripe due to the movements of stones near the shear zone, and it expands gradually with the increase of stone content. In addition, there is an approximately linear relation between the mean increment of internal friction angle and the stone content lying between 30% and 60%, and a concave nonlinear relation between the mean increment of cohesion and stone content, which are in good agreement with the existing research results.展开更多
Fetr6 is an underground mine using the stope-and-pillar mining method. As there was some evidence regarding pillar failure in this mine, improving works such as roof support and replacing existing pillars with concret...Fetr6 is an underground mine using the stope-and-pillar mining method. As there was some evidence regarding pillar failure in this mine, improving works such as roof support and replacing existing pillars with concrete pillars (CP) were carried out. During the construction of the second CP, in the space between the remaining pillars, one of the pillars failed leading to the progressive failure of other pillars until 4 000 m 2 of mine had collapsed within a few minutes. In this work, this phenomenon is described by applying both numerical and empirical methods and the respective results are compared. The results of numerical modelling are found to be closer to the actual condition than those of the empirical method. Also, a width-to-height (W/H) ratio less than 1, an inadequate support system and the absence of a detailed program for pillar recovery are shown to be the most important causes of the Domino failure in this mine.展开更多
The behavior of single bubble rising in quiescent shear-thinning tlmds was lnvestlgateO numerically by level set metnoa. number of bubbles in a large range of Reynolds number and Eotvos number were investigated includ...The behavior of single bubble rising in quiescent shear-thinning tlmds was lnvestlgateO numerically by level set metnoa. number of bubbles in a large range of Reynolds number and Eotvos number were investigated including spherical, oblate and spherical. The bubble shape and drag coefficient were compared with experimental results. It is observed that the simulated results show good conformity to experimental results over a wide range of Reynolds number. In addition, the detailed flow field based on the reference coordinate system moving with the bubble is obtained, and the relationship among flow field, bubble shape and velocity is discussed.展开更多
Columnar jointed rock mass with unique geometric and geological properties is one spectacular example of geometrical order in nature.Columnar joints are generally accepted to be formed by spatially uniform volume cont...Columnar jointed rock mass with unique geometric and geological properties is one spectacular example of geometrical order in nature.Columnar joints are generally accepted to be formed by spatially uniform volume contraction during cooling.In this paper,substantial field work was performed to study the geological characteristics of irregular columnar jointed basalt on the left bank dam foundation in the Baihetan Hydropower Station,where the columnar jointed rock mass is extensively exposed due to excavation.The quantitative measurements of the sizing of polygonal crack pattern of columnar joints and assessment of their degree of irregularity were summarized.Considering the irregularity of polygonal crack pattern,a modified Voronoi polygon(MVP)method was developed to model the special polygonal crack pattern of columnar joints.The new polygonal pattern obtained by the MVP method consists of a large number of irregular polygons,of which the degree of irregularity is consistent with the field measurement results.This method can reproduce the rapid evolution from an initial ideal regular hexagonal pattern to a final actual irregular polygonal pattern as the degree of irregularity increases.The compression tests of columnar jointed rock mass with different irregularity show that the geometric irregularity has a great influence on its mechanical properties.This numerical construction method provides a reliable way to reconstruct columnar joint structure with specific polygonal crack pattern,which is consistent with onsite columnar jointed basalt.展开更多
Several parameter identification methods of thermal response test were evaluated through numerical and experimental study.A three-dimensional finite-volume numerical model was established under the assumption that the...Several parameter identification methods of thermal response test were evaluated through numerical and experimental study.A three-dimensional finite-volume numerical model was established under the assumption that the soil thermal conductivity had been known in the simulation of thermal response test.The thermal response curve was firstly obtained through numerical calculation.Then,the accuracy of the numerical model was verified with measured data obtained through a thermal response test.Based on the numerical and experimental thermal response curves,the thermal conductivity of the soil was calculated by different parameter identification methods.The calculated results were compared with the assumed value and then the accuracy of these methods was evaluated.Furthermore,the effects of test time,variable data quality,borehole radius,initial ground temperature,and heat injection rate were analyzed.The results show that the method based on cylinder-source model has a low precision and the identified thermal conductivity decreases with an increase in borehole radius.For parameter estimation,the measuring accuracy of the initial temperature of the deep ground soil has greater effect on identified thermal conductivity.展开更多
Partition of unity based numerical manifold method can solve continuous and discontinuous problems in a unified framework with a two-cover system,i.e.,the mathematical cover and physical cover.However,renewal of the t...Partition of unity based numerical manifold method can solve continuous and discontinuous problems in a unified framework with a two-cover system,i.e.,the mathematical cover and physical cover.However,renewal of the topology of the two-cover system poses a challenge for multiple crack propagation problems and there are few references.In this study,a robust and efficient strategy is proposed to update the cover system of the numerical manifold method in simulation of multiple crack propagation problems.The proposed algorithm updates the cover system with a bottom-up process:1)identification of fractured manifold elements according to the previous and latest crack tip position;and 2)local topological update of the manifold elements,physical patches,block boundary loops,and non-persistent joint loops according to the scenario classification of the propagating crack.The proposed crack tracking strategy and classification of the renewal cases promote a robust and efficient cover renewal algorithm for multiple crack propagation analysis.Three crack propagation examples show that the proposed algorithm performs well in updating the cover system.This cover renewal methodology can be extended for numerical manifold method with polygonal mathematical covers.展开更多
Combined with the kinetic model of liquid film spreading, a new numerical method of solid-liquid-gas three-phase flow was developed for the moving of contact line, which was a hybrid method of computational fluid dyna...Combined with the kinetic model of liquid film spreading, a new numerical method of solid-liquid-gas three-phase flow was developed for the moving of contact line, which was a hybrid method of computational fluid dynamics and lattice Boltzmalm method （LBM）. By taking the effect of molecule force in droplet and the wall surface on liquid film into account, the changing law of contact angle with different surface tensions was analyzed on glass and aluminum foil surfaces. Compared with experimental results, the standard deviation by using LBM is less than 0.5°, which validates the feasibility of LBM simulation on the dynamic process of liquid film spreading. In addition, oscillations are discovered both at the initial and end phases. The phenomenon of retraction is also found and the maximum retraction angle is 7.58°. The obtained result shows that the retraction is proved to be correlative with precursor film by tracking the volume change of liquid film contour. Furthermore, non-dimensional coefficient 2 is introduced to measure the liquid film retraction capacity.展开更多
A class of modified parallel combined methods of real-time numerical simulation are presented for a stiff dynamic system. By combining the parallelism across the system with the parallelism across the method, and rela...A class of modified parallel combined methods of real-time numerical simulation are presented for a stiff dynamic system. By combining the parallelism across the system with the parallelism across the method, and relaxing the dependence of stage value computation on sampling time of input function, a class of modified real-time parallel combined methods are constructed. Stiff and nonstiff subsystems are solved in parallel on a parallel computer by a parallel Rosen-brock method and a parallel RK method, respectively. Their order conditions and convergences are discussed. The numerical simulation experiments show that this class of modified algorithms can get high speed and efficiency.展开更多
This paper deals with the numerical solution of initial value problems for systems of differential equations with two delay terms. We investigate the stability of adaptations of the θ-methods in the numerical solutio...This paper deals with the numerical solution of initial value problems for systems of differential equations with two delay terms. We investigate the stability of adaptations of the θ-methods in the numerical solution of test equations u'(t) = a 11 u(t) + a12v(t) + b11 u(t - τ1) + b12v(t-τ2,v'(t) = a21 u(t) + a22 v(t) + b21 u(t -τ1,) + b22 v(t -τ2), t>0,with initial conditionsu(t)=u0(t),v(t) =v0(t), t≤0.where aij, bij∈C, τj >0, i,j = 1,2,, and u0(t), v0(t)are continuous and complex valued. Sufficient conditions for the asymptotic stability of test equation are derived. Furthermore, with respect to an appropriate definition of stability for the numerical method, it is proved that the linear θ-method is stable if and only if 1/2≤θ≤1 and the one-leg θ-method is stable if and only if θ= 1.展开更多
The absorbing boundary is the key in numerical simulation of borehole radar.Perfect match layer(PML) was chosen as the absorbing boundary in numerical simulation of GPR.But CPML(convolutional perfect match layer) appr...The absorbing boundary is the key in numerical simulation of borehole radar.Perfect match layer(PML) was chosen as the absorbing boundary in numerical simulation of GPR.But CPML(convolutional perfect match layer) approach that we have chosen has the advantage of being media independent.Beginning with the Maxwell equations in a two-dimensional structure,numerical formulas of finite-difference time-domain(FDTD) method with CPML boundary condition for transverse electric(TE) or transverse magnetic(TM) wave are presented in details.Also,there are three models for borehole-GPR simulation.By analyzing the simulation results,the features of targets in GPR are obtained,which can provide a better interpretation of real radar data.The results show that CPML is well suited for the simulation of borehole-GPR.展开更多
In order to investigate the process of incremental sheet forming (ISF) through both experimental and numerical approaches, a three-dimensional elasto-plastic finite element model (FEM) was developed to simulate the pr...In order to investigate the process of incremental sheet forming (ISF) through both experimental and numerical approaches, a three-dimensional elasto-plastic finite element model (FEM) was developed to simulate the process and the simulated results were compared with those of experiment. The results of numerical simulations, such as the strain history and distribution, the stress state and distribution, sheet thickness distribution, etc, were discussed in details, and the influences of process parameters on these results were also analyzed. The simulated results of the radial strain and the thickness distribution are in good agreement with experimental results. The simulations reveal that the deformation is localized around the tool and constantly remains close to a plane strain state. With decreasing depth step, increasing tool diameter and wall inclination angle, the axial stress reduces, leading to less thinning and more homogeneous plastic strain and thickness distribution. During ISF, the plastic strain increases stepwise under the action of the tool. Each increase in plastic strain is accompanied by hydrostatic pressure, which explains why obtainable deformation using ISF exceeds the forming limits of conventional sheet forming.展开更多
In order to obtain satisfactory mechanical properties for the cam used in high-power ship diesel engines, a new quenching technology was proposed by designing a two-stage quenching process with an alkaline bath as the...In order to obtain satisfactory mechanical properties for the cam used in high-power ship diesel engines, a new quenching technology was proposed by designing a two-stage quenching process with an alkaline bath as the quenching medium. To demonstrate the effectiveness of the proposed new quenching technology, both numerical analysis and experimental study were performed. The new quenching technology was analyzed using finite element method. The combined effects of the temperature, stress and microstructure fields were investigated considering nonlinear material properties. Finally, an experimental study was performed to verify the effectiveness of the proposed new quenching technology. The numerical results show that internal stress is affected by both thermal stress and transformation stress. In addition, the direction of the internal stress is changed several times due to thermal interaction and microstructure evolution during the quenching process. The experimental results show that the proposed new quenching technology significantly improves the mechanical properties and microstructures of the cam. The tensile strength, the impact resistance and the hardness value of the cam by the proposed new quenching technology are improved by 4.3%, 8.9% and 3.5% compared with those by the traditional quenching technology. Moreover, the residual stress and cam shape deformation are reduced by 40.0% and 48.9% respectively for the cam manufactured by the new quenching technology.展开更多
The instantaneous variations of the hydraulic characteristics take place in centrifugal pumps during their start-up,shutdown and other variable speed operations.In this paper,the variable speed method was proposed to ...The instantaneous variations of the hydraulic characteristics take place in centrifugal pumps during their start-up,shutdown and other variable speed operations.In this paper,the variable speed method was proposed to simulate the transient internal flow field and the external performance of the pump during starting and stopping periods.The terms of accelerations due to variable speeds in the flow governing equations were analyzed in a multiple reference of frame(MRF).A transient CFD simulation was performed for a typical centrifugal pump by using ANSYS-CFX with the standard k-εturbulence model.The entire simulation process was composed of four stages:start-up,normal run,shutdown and post-shutdown.The function of rotating speed with regard to time was set by CEL language directly into the impeller domain in the pre-processor of the software to conduct variable speed simulation.The variations of the flow field in the centrifugal pump were obtained from the transient simulation.The changing laws of flow rate,head and other performance parameters over time were also analyzed and summarized.展开更多
Being aimed at the inside wall wrinkling and sinking phenomenon of palladium-yttrium alloy thin wall spiral tube used for preparation of high purity hydrogen, extraction of hydrogen isotope, and purification and separ...Being aimed at the inside wall wrinkling and sinking phenomenon of palladium-yttrium alloy thin wall spiral tube used for preparation of high purity hydrogen, extraction of hydrogen isotope, and purification and separation of hydrogen in the winding process, this article analyzed the reasons for above phenomena, established a numerical simulation model of winding process of above tube, using elastic-plastic Finite Element method analyzed the max. tensile stress and max. compression stress and their locations, thereby provides a theory base for the control of working forming course of thin wall spiral tube.展开更多
In this paper, the internal fluid motion of a jet system is described by the Navier Stokes mechanics equations. For the simulation of the motion, the penalty function finite element method is used, and the velocity ve...In this paper, the internal fluid motion of a jet system is described by the Navier Stokes mechanics equations. For the simulation of the motion, the penalty function finite element method is used, and the velocity vectors and stream function curves are obtained. Using the Prandtl theory, this paper derives the free jet velocity and the jet bunch width in a half-space, the latter of which is amended by experiment. The results obtained in this paper are applied to micro-type high pressure water jet cleaner and the ejector of rocket engine.展开更多
In order to investigate the temperature distribution, a three-dimensional finite element model （FEM） was developed to simulate the temperature regime in the channels of double-loop inductor, and the simulated result...In order to investigate the temperature distribution, a three-dimensional finite element model （FEM） was developed to simulate the temperature regime in the channels of double-loop inductor, and the simulated results were compared with experimental data from low load trials of a 400 kW inductor. The results of numerical simulations, such as the temperature and Joule heating rate, show reasonable correlation with experimental data. The results indicate that Joule heating rate and the temperature reach the maximum at the comers and the minimum at the centre of the cross-section area. The temperature difference between the inlet and outlet is in an inverse proportion to mass transport. Joule heating rate and the temperature are directly proportional to power frequency. It is concluded that mass transport and power frequency play a critical role in determining the temperature regime and Joule heating rate, the relative permeability of the magnetic core shows no significant influence on temperature regime and Joule heating rate, when the relative permeability varies from 5 000 to 10 000.展开更多
基金financially supported by the National Natural Science Foundation of China(Grant Nos.52271317 and 52071149)the Fundamental Research Funds for the Central Universities(HUST:2019kfy XJJS007)。
文摘To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based simulation(NNS)method with higher accuracy and better efficiency was proposed.The NNS method consisted of three main steps.First,the parameters of blast loads,including the peak pressures and impulses of cylindrical charges with different aspect ratios(L/D)at different stand-off distances and incident angles were obtained by two-dimensional numerical simulations.Subsequently,incident shape factors of cylindrical charges with arbitrary aspect ratios were predicted by a neural network.Finally,reflected shape factors were derived and implemented into the subroutine of the ABAQUS code to modify the CONWEP model,including modifications of impulse and overpressure.The reliability of the proposed NNS method was verified by related experimental results.Remarkable accuracy improvement was acquired by the proposed NNS method compared with the unmodified CONWEP model.Moreover,huge efficiency superiority was obtained by the proposed NNS method compared with the CEL method.The proposed NNS method showed good accuracy when the scaled distance was greater than 0.2 m/kg^(1/3).It should be noted that there is no need to generate a new dataset again since the blast loads satisfy the similarity law,and the proposed NNS method can be directly used to simulate the blast loads generated by different cylindrical charges.The proposed NNS method with high efficiency and accuracy can be used as an effective method to analyze the dynamic response of structures under blast loads,and it has significant application prospects in designing protective structures.
基金Project(41630642)supported by the Key Project of National Natural Science Foundation of ChinaProject(51974360)supported by the National Natural Science Foundation of ChinaProject(2018JJ3656)supported by the Natural Science Foundation of Hunan Province,China。
文摘In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step for the dynamic response simulation of rock mass in a high in-situ stress field.In this paper,stress initialization methods,including their principles and operating procedures for reproducing steady in-situ stress state in LS-DYNA,are first introduced.Then the most popular four methods,i.e.,explicit dynamic relaxation(DR)method,implicit-explicit sequence method,Dynain file method and quasi-static method,are exemplified through a case analysis by using the RHT and plastic hardening rock material models to simulate rock blasting under in-situ stress condition.Based on the simulations,it is concluded that the stress initialization results obtained by implicit-explicit sequence method and dynain file method are closely related to the rock material model,and the explicit DR method has an obvious advantage in solution time when compared to other methods.Besides that,it is recommended to adopt two separate analyses for the whole numerical simulation of rock mass under the combined action of in-situ stress and dynamic disturbance.
基金Project(2011CB013504) supported by the National Basic Research Program(973 Program)of ChinaProject(2013BAB06B01) supported by the National Science&Technology Pillar Program during the Twelfth Five-year Plan Period+2 种基金Projects(11772118,51479049,51709282) supported by the National Natural Science Foundation of ChinaProject(2017M620838) supported by the Postdoctoral Science Foundation of ChinaProject(487237) supported by the Natural Sciences and Engineering Research Council of Canada
文摘Outwash deposit is a unique type of geological materials, and its features such as heterogeneity, discontinuity and nonlinearity determine the complexity of mechanical characteristics and failure mechanism. In this work, random meso-structure of outwash deposits was constructed by the technique of computer random simulation based on characteristics of its meso-structure in the statistical sense and some simplifications, and a series of large direct shear tests on numerical samples of outwash deposits with stone contents of 15%, 30%, 45% and 60% were conducted using the discrete element method to further investigate its mechanical characteristics and failure mechanism under external load. The results show that the deformation characteristics and shear strength of outwash deposits are to some extent improved with the increase of stone content, and the shear stress–shear displacement curves of outwash deposits show great differences at the post-peak stage due to the random spatial distribution and content of stones. From the mesoscopic view, normal directions of contacts between "soil" and "stone" particles undergo apparent deflection as the shear displacement continues during the shearing process, accompanying redistribution of the magnitude of contact forces during the shearing process. For outwash deposits, the shear zone formed after shear failure is an irregular stripe due to the movements of stones near the shear zone, and it expands gradually with the increase of stone content. In addition, there is an approximately linear relation between the mean increment of internal friction angle and the stone content lying between 30% and 60%, and a concave nonlinear relation between the mean increment of cohesion and stone content, which are in good agreement with the existing research results.
文摘Fetr6 is an underground mine using the stope-and-pillar mining method. As there was some evidence regarding pillar failure in this mine, improving works such as roof support and replacing existing pillars with concrete pillars (CP) were carried out. During the construction of the second CP, in the space between the remaining pillars, one of the pillars failed leading to the progressive failure of other pillars until 4 000 m 2 of mine had collapsed within a few minutes. In this work, this phenomenon is described by applying both numerical and empirical methods and the respective results are compared. The results of numerical modelling are found to be closer to the actual condition than those of the empirical method. Also, a width-to-height (W/H) ratio less than 1, an inadequate support system and the absence of a detailed program for pillar recovery are shown to be the most important causes of the Domino failure in this mine.
基金Project(21406141)supported by the National Natural Science Foundation of ChinaProject(20141078)supported by the Scientific Research Starting Foundation for Doctors of Liaoning Province,China+1 种基金Project(L2014060)supported by the Foundation of Department of Education of Liaoning Province,ChinaProject(157B21)supported by the Scientific Research Starting Foundation for Doctors of Shenyang Aerospace University,China
文摘The behavior of single bubble rising in quiescent shear-thinning tlmds was lnvestlgateO numerically by level set metnoa. number of bubbles in a large range of Reynolds number and Eotvos number were investigated including spherical, oblate and spherical. The bubble shape and drag coefficient were compared with experimental results. It is observed that the simulated results show good conformity to experimental results over a wide range of Reynolds number. In addition, the detailed flow field based on the reference coordinate system moving with the bubble is obtained, and the relationship among flow field, bubble shape and velocity is discussed.
基金Projects(51621006,51779251)supported by the National Natural Science Foundation of China。
文摘Columnar jointed rock mass with unique geometric and geological properties is one spectacular example of geometrical order in nature.Columnar joints are generally accepted to be formed by spatially uniform volume contraction during cooling.In this paper,substantial field work was performed to study the geological characteristics of irregular columnar jointed basalt on the left bank dam foundation in the Baihetan Hydropower Station,where the columnar jointed rock mass is extensively exposed due to excavation.The quantitative measurements of the sizing of polygonal crack pattern of columnar joints and assessment of their degree of irregularity were summarized.Considering the irregularity of polygonal crack pattern,a modified Voronoi polygon(MVP)method was developed to model the special polygonal crack pattern of columnar joints.The new polygonal pattern obtained by the MVP method consists of a large number of irregular polygons,of which the degree of irregularity is consistent with the field measurement results.This method can reproduce the rapid evolution from an initial ideal regular hexagonal pattern to a final actual irregular polygonal pattern as the degree of irregularity increases.The compression tests of columnar jointed rock mass with different irregularity show that the geometric irregularity has a great influence on its mechanical properties.This numerical construction method provides a reliable way to reconstruct columnar joint structure with specific polygonal crack pattern,which is consistent with onsite columnar jointed basalt.
基金Project(xjj20100078) supported by the Fundamental Research Funds for the Central Universities in China
文摘Several parameter identification methods of thermal response test were evaluated through numerical and experimental study.A three-dimensional finite-volume numerical model was established under the assumption that the soil thermal conductivity had been known in the simulation of thermal response test.The thermal response curve was firstly obtained through numerical calculation.Then,the accuracy of the numerical model was verified with measured data obtained through a thermal response test.Based on the numerical and experimental thermal response curves,the thermal conductivity of the soil was calculated by different parameter identification methods.The calculated results were compared with the assumed value and then the accuracy of these methods was evaluated.Furthermore,the effects of test time,variable data quality,borehole radius,initial ground temperature,and heat injection rate were analyzed.The results show that the method based on cylinder-source model has a low precision and the identified thermal conductivity decreases with an increase in borehole radius.For parameter estimation,the measuring accuracy of the initial temperature of the deep ground soil has greater effect on identified thermal conductivity.
基金Project(51321065,51479191,11672360)supported by the National Natural Science Foundation of China。
文摘Partition of unity based numerical manifold method can solve continuous and discontinuous problems in a unified framework with a two-cover system,i.e.,the mathematical cover and physical cover.However,renewal of the topology of the two-cover system poses a challenge for multiple crack propagation problems and there are few references.In this study,a robust and efficient strategy is proposed to update the cover system of the numerical manifold method in simulation of multiple crack propagation problems.The proposed algorithm updates the cover system with a bottom-up process:1)identification of fractured manifold elements according to the previous and latest crack tip position;and 2)local topological update of the manifold elements,physical patches,block boundary loops,and non-persistent joint loops according to the scenario classification of the propagating crack.The proposed crack tracking strategy and classification of the renewal cases promote a robust and efficient cover renewal algorithm for multiple crack propagation analysis.Three crack propagation examples show that the proposed algorithm performs well in updating the cover system.This cover renewal methodology can be extended for numerical manifold method with polygonal mathematical covers.
基金Project(U1261107)supported by the National Natural Science Foundation of China
文摘Combined with the kinetic model of liquid film spreading, a new numerical method of solid-liquid-gas three-phase flow was developed for the moving of contact line, which was a hybrid method of computational fluid dynamics and lattice Boltzmalm method （LBM）. By taking the effect of molecule force in droplet and the wall surface on liquid film into account, the changing law of contact angle with different surface tensions was analyzed on glass and aluminum foil surfaces. Compared with experimental results, the standard deviation by using LBM is less than 0.5°, which validates the feasibility of LBM simulation on the dynamic process of liquid film spreading. In addition, oscillations are discovered both at the initial and end phases. The phenomenon of retraction is also found and the maximum retraction angle is 7.58°. The obtained result shows that the retraction is proved to be correlative with precursor film by tracking the volume change of liquid film contour. Furthermore, non-dimensional coefficient 2 is introduced to measure the liquid film retraction capacity.
基金This project was supported by the National Natural Science Foundation of China (19871080).
文摘A class of modified parallel combined methods of real-time numerical simulation are presented for a stiff dynamic system. By combining the parallelism across the system with the parallelism across the method, and relaxing the dependence of stage value computation on sampling time of input function, a class of modified real-time parallel combined methods are constructed. Stiff and nonstiff subsystems are solved in parallel on a parallel computer by a parallel Rosen-brock method and a parallel RK method, respectively. Their order conditions and convergences are discussed. The numerical simulation experiments show that this class of modified algorithms can get high speed and efficiency.
文摘This paper deals with the numerical solution of initial value problems for systems of differential equations with two delay terms. We investigate the stability of adaptations of the θ-methods in the numerical solution of test equations u'(t) = a 11 u(t) + a12v(t) + b11 u(t - τ1) + b12v(t-τ2,v'(t) = a21 u(t) + a22 v(t) + b21 u(t -τ1,) + b22 v(t -τ2), t>0,with initial conditionsu(t)=u0(t),v(t) =v0(t), t≤0.where aij, bij∈C, τj >0, i,j = 1,2,, and u0(t), v0(t)are continuous and complex valued. Sufficient conditions for the asymptotic stability of test equation are derived. Furthermore, with respect to an appropriate definition of stability for the numerical method, it is proved that the linear θ-method is stable if and only if 1/2≤θ≤1 and the one-leg θ-method is stable if and only if θ= 1.
基金Project(41174061) supported by the National Natural Science Foundation of ChinaProject(2011QNZT011) supported by the Free Exploration Program of Central South University,China
文摘The absorbing boundary is the key in numerical simulation of borehole radar.Perfect match layer(PML) was chosen as the absorbing boundary in numerical simulation of GPR.But CPML(convolutional perfect match layer) approach that we have chosen has the advantage of being media independent.Beginning with the Maxwell equations in a two-dimensional structure,numerical formulas of finite-difference time-domain(FDTD) method with CPML boundary condition for transverse electric(TE) or transverse magnetic(TM) wave are presented in details.Also,there are three models for borehole-GPR simulation.By analyzing the simulation results,the features of targets in GPR are obtained,which can provide a better interpretation of real radar data.The results show that CPML is well suited for the simulation of borehole-GPR.
基金Project(50175034) supported by the National Natural Science Foundation of China
文摘In order to investigate the process of incremental sheet forming (ISF) through both experimental and numerical approaches, a three-dimensional elasto-plastic finite element model (FEM) was developed to simulate the process and the simulated results were compared with those of experiment. The results of numerical simulations, such as the strain history and distribution, the stress state and distribution, sheet thickness distribution, etc, were discussed in details, and the influences of process parameters on these results were also analyzed. The simulated results of the radial strain and the thickness distribution are in good agreement with experimental results. The simulations reveal that the deformation is localized around the tool and constantly remains close to a plane strain state. With decreasing depth step, increasing tool diameter and wall inclination angle, the axial stress reduces, leading to less thinning and more homogeneous plastic strain and thickness distribution. During ISF, the plastic strain increases stepwise under the action of the tool. Each increase in plastic strain is accompanied by hydrostatic pressure, which explains why obtainable deformation using ISF exceeds the forming limits of conventional sheet forming.
基金Project(50875268) supported by the National Natural Science Foundation of China Project(CSTC2008AB3057) supported by Foundation of Chongqing Science and Technology Commission, China+1 种基金 Project(108107) supported by the Key Project of Ministry of Education of China Project(50925518) supported by the National Science Fund for Distinguished Young Scholars
文摘In order to obtain satisfactory mechanical properties for the cam used in high-power ship diesel engines, a new quenching technology was proposed by designing a two-stage quenching process with an alkaline bath as the quenching medium. To demonstrate the effectiveness of the proposed new quenching technology, both numerical analysis and experimental study were performed. The new quenching technology was analyzed using finite element method. The combined effects of the temperature, stress and microstructure fields were investigated considering nonlinear material properties. Finally, an experimental study was performed to verify the effectiveness of the proposed new quenching technology. The numerical results show that internal stress is affected by both thermal stress and transformation stress. In addition, the direction of the internal stress is changed several times due to thermal interaction and microstructure evolution during the quenching process. The experimental results show that the proposed new quenching technology significantly improves the mechanical properties and microstructures of the cam. The tensile strength, the impact resistance and the hardness value of the cam by the proposed new quenching technology are improved by 4.3%, 8.9% and 3.5% compared with those by the traditional quenching technology. Moreover, the residual stress and cam shape deformation are reduced by 40.0% and 48.9% respectively for the cam manufactured by the new quenching technology.
文摘The instantaneous variations of the hydraulic characteristics take place in centrifugal pumps during their start-up,shutdown and other variable speed operations.In this paper,the variable speed method was proposed to simulate the transient internal flow field and the external performance of the pump during starting and stopping periods.The terms of accelerations due to variable speeds in the flow governing equations were analyzed in a multiple reference of frame(MRF).A transient CFD simulation was performed for a typical centrifugal pump by using ANSYS-CFX with the standard k-εturbulence model.The entire simulation process was composed of four stages:start-up,normal run,shutdown and post-shutdown.The function of rotating speed with regard to time was set by CEL language directly into the impeller domain in the pre-processor of the software to conduct variable speed simulation.The variations of the flow field in the centrifugal pump were obtained from the transient simulation.The changing laws of flow rate,head and other performance parameters over time were also analyzed and summarized.
文摘Being aimed at the inside wall wrinkling and sinking phenomenon of palladium-yttrium alloy thin wall spiral tube used for preparation of high purity hydrogen, extraction of hydrogen isotope, and purification and separation of hydrogen in the winding process, this article analyzed the reasons for above phenomena, established a numerical simulation model of winding process of above tube, using elastic-plastic Finite Element method analyzed the max. tensile stress and max. compression stress and their locations, thereby provides a theory base for the control of working forming course of thin wall spiral tube.
文摘In this paper, the internal fluid motion of a jet system is described by the Navier Stokes mechanics equations. For the simulation of the motion, the penalty function finite element method is used, and the velocity vectors and stream function curves are obtained. Using the Prandtl theory, this paper derives the free jet velocity and the jet bunch width in a half-space, the latter of which is amended by experiment. The results obtained in this paper are applied to micro-type high pressure water jet cleaner and the ejector of rocket engine.
基金Project(50876116) supported by the National Natural Science Foundation of ChinaProject(2007CK3077) supported by the Innovative Program of Hunan Science and Technology Agency, ChinaProject(1343-77225) supported by the Graduate School of Central South University, China
文摘In order to investigate the temperature distribution, a three-dimensional finite element model （FEM） was developed to simulate the temperature regime in the channels of double-loop inductor, and the simulated results were compared with experimental data from low load trials of a 400 kW inductor. The results of numerical simulations, such as the temperature and Joule heating rate, show reasonable correlation with experimental data. The results indicate that Joule heating rate and the temperature reach the maximum at the comers and the minimum at the centre of the cross-section area. The temperature difference between the inlet and outlet is in an inverse proportion to mass transport. Joule heating rate and the temperature are directly proportional to power frequency. It is concluded that mass transport and power frequency play a critical role in determining the temperature regime and Joule heating rate, the relative permeability of the magnetic core shows no significant influence on temperature regime and Joule heating rate, when the relative permeability varies from 5 000 to 10 000.
