Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The t...Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The traditional thermal elastic-plastic finite element method(TEP-FEM)can accurately predict welding deformation.However,its efficiency is low because of the complex nonlinear transient computation,making it difficult to meet the needs of rapid engineering evaluation.To address this challenge,this study proposes an efficient prediction method for welding deformation in marine thin plate butt welds.This method is based on the coupled temperature gradient-thermal strain method(TG-TSM)that integrates inherent strain theory with a shell element finite element model.The proposed method first extracts the distribution pattern and characteristic value of welding-induced inherent strain through TEP-FEM analysis.This strain is then converted into the equivalent thermal load applied to the shell element model for rapid computation.The proposed method-particularly,the gradual temperature gradient-thermal strain method(GTG-TSM)-achieved improved computational efficiency and consistent precision.Furthermore,the proposed method required much less computation time than the traditional TEP-FEM.Thus,this study lays the foundation for future prediction of welding deformation in more complex marine thin plates.展开更多
As a pyrometallurgical process,circulating fluidized bed(CFB) roasting has good potential for application in desulfurization of high-sulfur bauxite.The gas-solid distribution and reaction during CFB roasting of high-s...As a pyrometallurgical process,circulating fluidized bed(CFB) roasting has good potential for application in desulfurization of high-sulfur bauxite.The gas-solid distribution and reaction during CFB roasting of high-sulfur bauxite were simulated using the computational particle fluid dynamics(CPFD) method.The effect of primary air flow velocity on particle velocity,particle volume distribution,furnace temperature distribution and pressure distribution were investigated.Under the condition of the same total flow of natural gas,the impact of the number of inlets on the desulfurization efficiency,atmosphere mass fraction distribution and temperature distribution in the furnace was further investigated.展开更多
After a long period of water flooding development,the oilfield has entered the middle and high water cut stage.The physical properties of reservoirs are changed by water erosion,which directly impacts reservoir develo...After a long period of water flooding development,the oilfield has entered the middle and high water cut stage.The physical properties of reservoirs are changed by water erosion,which directly impacts reservoir development.Conventional numerical reservoir simulation methodologies typically employ static assumptions for model construction,presuming invariant reservoir geological parameters throughout the development process while neglecting the reservoir’s temporal evolution characteristics.Although such simplifications reduce computational complexity,they introduce substantial descriptive inaccuracies.Therefore,this paper proposes a meshless numerical simulation method for reservoirs that considers time-varying characteristics.This method avoids the meshing in traditional numerical simulation methods.From the fluid flow perspective,the reservoir’s computational domain is discretized into a series of connection units.An influence domain with a certain radius centered on the nodes is selected,and one-dimensional connection units are established between the nodes to achieve the characterization of the flow topology structure of the reservoir.In order to reflect the dynamic evolution of the reservoir’s physical properties during the water injection development process,the time-varying characteristics are incorporated into the formula of the seepage characteristic parameters in the meshless calculation.The change relationship of the permeability under different surface fluxes is considered to update the calculated connection conductivity in real time.By combining with the seepage control equation for solution,a time-varying meshless numerical simulation method is formed.The results show that compared with the numerical simulationmethod of the connection elementmethod(CEM)that only considers static parameters,this method has higher simulation accuracy and can better simulate the real migration and distribution of oil and water in the reservoir.Thismethod improves the accuracy of reservoir numerical simulation and the development effect of oilfields,providing a scientific basis for optimizing the water injection strategy,adjusting the production plan,and extending the effective production cycle of the oilfield.展开更多
Deepwater drilling riser is the key equipment connecting the subsea wellhead and floating drilling platform.Due to complex marine environment,vortex-induced vibration(ViV)will be generated on riser,which will induce f...Deepwater drilling riser is the key equipment connecting the subsea wellhead and floating drilling platform.Due to complex marine environment,vortex-induced vibration(ViV)will be generated on riser,which will induce fatigue failure and even cause unpredictable drilling accidents.Therefore,it is important to study the ViV characteristics of deepwater drilling riser and reveal the main controlling factors for ensuring the safe and efficient operation of deepwater drilling engineering.In this paper,the ViV of deepwater drilling riser is numerically simulated in time domain based on the discrete vortex method(DvM).A hydrodynamic analysis model and governing equation of VIV is proposed with considering the effect of riser motion using DVM and slice method,where the governing equation is solved by Runge-Kutta method.Model validation is performed,which verified the correctness and accuracy of the mechanical model and the solution method.On this basis,the influence of the number of control points,current velocity,riser outer diameter,shear flow and top tension on the ViV characteristics of deepwater drilling risers are discussed in detail.The results show that with the increase of current velocity,the vibration amplitude of deepwater drilling riser decreases obviously,while the vibration frequency increases gradually.However,if the outer diameter of riser increases,the vibration amplitude increases,while the vibration frequency decreases gradually.The top tension also has great influence on the VIV of riser.When the top tension is 1.25 G,the VIV is suppressed to a certain extent.This study has guiding significance for optimal design and engineering control of deepwater drilling riser.展开更多
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).展开更多
In order to predict the long-term rutting of asphalt pavement, the effective temperature for pavement rutting is calculated using the numerical simulation method. The transient temperature field of asphalt pavement wa...In order to predict the long-term rutting of asphalt pavement, the effective temperature for pavement rutting is calculated using the numerical simulation method. The transient temperature field of asphalt pavement was simulated based on actual meteorological data of Nanjing. 24-hour rutting development under a transient temperature field was calculated in each month. The rutting depth accumulated under the static temperature field was also estimated and the relationship between constant temperature parameters was analyzed. Then the effective temperature for pavement rutting was determined based on the rutting equivalence principle. The results show that the monthly effective temperature is above 40 t in July and August, while in June and September it ranges from 30 to 40 Rutting development can be ignored when the monthly effective temperature is less than 30 t. The yearly effective temperature for rutting in Nanjing is around 38. 5 t. The long-term rutting prediction model based on the effective temperature can reflect the influences of meteorological factors and traffic time distribution.展开更多
LS-DYNA program and the principle of ALE method were introduced, and the target features of the reinforced concrete penetration were analyzed by using the D material model and the ALE method. A numerical simulation ha...LS-DYNA program and the principle of ALE method were introduced, and the target features of the reinforced concrete penetration were analyzed by using the D material model and the ALE method. A numerical simulation has been done to show the penetration visually and veritably. The simulation results are analyzed carefully and explicitly prove their significance to the research of reinforced concrete penetration.展开更多
Shale gas production involves complex gas-water two-phase flow,with flow patterns in proppant-filled fractures playing a critical role in determining production efficiency.In this study,3D geometric models of 40/70 me...Shale gas production involves complex gas-water two-phase flow,with flow patterns in proppant-filled fractures playing a critical role in determining production efficiency.In this study,3D geometric models of 40/70 mesh ceramic particles and quartz sand proppant clusters were elaborated using computed tomography(CT)scanning.These models were used to develop a numerical simulation framework based on the lattice Boltzmann method(LBM),enabling the investigation of gas-water flow behavior within proppant-filled fractures under varying driving forces and surface tensions.Simulation results at a closure pressure of 15 MPa have revealed that ceramic particles exhibit a simpler and more porous internal structure than quartz sand of the same size.Under identical flow conditions,ceramic proppants demonstrate higher fluid replacement efficiency.Replacement efficiency increases with higher porosity,greater driving force,and lower surface tension.Furthermore,fluid displacement is strongly influenced by pore geometry:flow is faster in straighter and wider channels,with preferential movement through larger pores forming dominant flow paths.The replacement velocity exhibits a characteristic time evolution,initially rapid,then gradually decreasing,correlating positively with the development of these dominant channels.展开更多
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.展开更多
Considering the limitation of computational capacity, a new finite element solution is used to simulate the welding deformation of the side sill of railroad car＇ s bogie frame based on the local-global method. Firstl...Considering the limitation of computational capacity, a new finite element solution is used to simulate the welding deformation of the side sill of railroad car＇ s bogie frame based on the local-global method. Firstly, a volumetric heat source defined by a double ellipsoid is adopted to simulate the thermal distributions of the arc welding process. And then, the local models extracted from the global model are computed with refined meshes. On these bases, the global distortions of the subject studied are ascertained by transferring the inner forces of computed local models to the global model. It indicates that the local-global method is feasible for simulating the large welded structures by comparing the computed results with the corresponding actual measured values. The work provides basis for optimizing the welding sequence and clamping conditions, and has theoretical values and engineering significance in the integral design, manufacturing technique selection of the bogie frame, as well as other kinds of large welded structures.展开更多
The mechanism of stress generation and propagation by detonation loading in five separate independent advance of ore breaking patterns is discussed in the paper. An elastic numerical model was developed using AN- SYS/...The mechanism of stress generation and propagation by detonation loading in five separate independent advance of ore breaking patterns is discussed in the paper. An elastic numerical model was developed using AN- SYS/LS-DYNA 3D Nonlinear Dynamic Finite Element Software. In this package ANSYS is the preprocessor and LS-DYNA is the postprocessor. Numerical models in the paper to actual were l：10 and the element mesh was dissected in scanning mode utilizing the symmetry characteristics of the numerical model. Five different advance rates were studied. Parameters, such as the time required to maximum stress, the action time of the available stress, the maximum velocity of the nodes, the stress penetration time, the magnitude of the stress peak and the time duration for high stress were numerically simulated. The 2.2 m advance appeared optimum from an analysis of the simulation results. The results from numerical simulation have been validated by tests with physical models.展开更多
In this paper, a method to develop a hierarchy of explicit recursion formulas for numerical simulation in an irregular grid for scalar wave equations is presented and its accuracy is illustrated via 2-D and 1-D models...In this paper, a method to develop a hierarchy of explicit recursion formulas for numerical simulation in an irregular grid for scalar wave equations is presented and its accuracy is illustrated via 2-D and 1-D models. Approaches to develop the stable formulas which are of 2M-order accuracy in both time and space with Mbeing a positive integer for regular grids are discussed and illustrated by constructing the second order （M= 1） and the fourth order （M = 2） recursion formulas.展开更多
Non-pillar continuous mining(NPCM) is regarded as a high-efficient, high-level and one-step mining technology, which can be divided into two substopes. Back fill stability status in substope I, which directly influenc...Non-pillar continuous mining(NPCM) is regarded as a high-efficient, high-level and one-step mining technology, which can be divided into two substopes. Back fill stability status in substope I, which directly influence the loss rate and dilution rate, etc, will determine whether the experimental research is successful or not. By employing energy method of limit analysis and finite element numerical simulation method, the critical backfill height was determined under the prerequisite condition of its stability, which put forward theoretical basis for reasonable and correct selection of backfill’s parameters. The result showed that the first backfill could not keep stable for NPCM, while the other was able to.展开更多
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.展开更多
The electrophysical property of saturated rocks is very important for reservoir identification and evaluation. In this paper, the lattice Boltzmann method （LBM） was used to study the electrophysical property of rock...The electrophysical property of saturated rocks is very important for reservoir identification and evaluation. In this paper, the lattice Boltzmann method （LBM） was used to study the electrophysical property of rock saturated with fluid because of its advantages over conventional numerical approaches in handling complex pore geometry and boundary conditions. The digital core model was constructed through the accumulation of matrix grains based on their radius distribution obtained by the measurements of core samples. The flow of electrical current through the core model saturated with oil and water was simulated on the mesoscopic scale to reveal the non-Archie relationship between resistivity index and water saturation （I-Sw）. The results from LBM simulation and laboratory measurements demonstrated that the I-Sw relation in the range of low water saturation was generally not a straight line in the log-log coordinates as described by the Archie equation. We thus developed a new equation based on numerical simulation and physical experiments. This new equation was used to fit the data from laboratory core measurements and previously published data. Determination of fluid saturation and reservoir evaluation could be significantly improved by using the new equation.展开更多
The numerical simulation model of single roller melt spinning for rapid quenching process of NdFeB alloy was built,and the vacuum chamber,cooling roller and sample were taken into account as a system.The existing matu...The numerical simulation model of single roller melt spinning for rapid quenching process of NdFeB alloy was built,and the vacuum chamber,cooling roller and sample were taken into account as a system.The existing mature technology was in order to verify the correctness of simulation.The rapid quenching ribbons with different roll speeds were used as the simulation objects.The results of the numerical simulation and experiments show that the validity of the model has been testified and the reasons of the formation of complete quenching ribbons and by-product have been explained.The experimental thickness of the ribbons is proportional to the theoretical thickness.In the same spray condition,with the roll speed increasing,the thickness decreases linearly.At the speed range of25-30 m·s^(-1),the simulated calculation date is close to the experimental date,which can be considered as an ideal technological parameter.展开更多
Numerical wave tanks are widely-acknowledged tools in studying waves and wave-structure interactions. They can generate waves under realistic scales and offers more information on the fluid field. However, most numeri...Numerical wave tanks are widely-acknowledged tools in studying waves and wave-structure interactions. They can generate waves under realistic scales and offers more information on the fluid field. However, most numerical wave tanks suffer from issues known as the numerical dissipation and numerical dispersion. The former causes wave energy to be slowly dissipated and the latter shifts wave frequencies during wave propagation. This paper proposes a simple method of depressing numerical dissipation effects on the basis of solving Euler equations using the finite difference method(FDM). The wave propagation solutions are solved analytically taking into account the influence of the damping terms. The main idea of the method is to append a source term to the momentum equation, whose strength is determined by how strong the numerical damping effect is. The method is verified by successfully depressing numerical effects during the simulation of regular linear waves, Stokes waves and irregular waves. By applying the method, wave energy is able to be close to its initial value after long distance of travel.展开更多
Wire arc additive manufacture(WAAM) is a new technique to fabricate large-scale complex aluminum alloy components.However, the performance of the parts is critically influenced by residual stresses and deformation. A ...Wire arc additive manufacture(WAAM) is a new technique to fabricate large-scale complex aluminum alloy components.However, the performance of the parts is critically influenced by residual stresses and deformation. A sequentially thermal-mechanical coupled model of residual stress and deformation for aluminum alloy WAAM parts was established based on commercial FE software ABAQUS. The temperature field was calculated by the moving heat source(MHS) method. The temperature function was obtained according to the distribution of the peak temperature. Furthermore, the MHS method and segmented temperature function(STF) method were used to calculate the residual stress and deformation. The results show that the STF method satisfies both the efficiency and accuracy requirements. 1-segment, 3-segment, and 5-segment STF methods can shorten the time for mechanical analysis by 91%, 79%, 63%, respectively.The error of the residual stress and deformation are all less than 20%. STF method provides an effective way to predict the residual stress and deformation of large-scale WAAM parts.展开更多
Hydraulic fracturing （HF） technique has been extensively used for the exploitation of unconventional oiland gas reservoirs. HF enhances the connectivity of less permeable oil and gas-bearing rock formationsby fluid ...Hydraulic fracturing （HF） technique has been extensively used for the exploitation of unconventional oiland gas reservoirs. HF enhances the connectivity of less permeable oil and gas-bearing rock formationsby fluid injection, which creates an interconnected fracture network and increases the hydrocarbonproduction. Meanwhile, microseismic （MS） monitoring is one of the most effective approaches to evaluatesuch stimulation process. In this paper, the combined finite-discrete element method （FDEM） isadopted to numerically simulate HF and associated MS. Several post-processing tools, includingfrequency-magnitude distribution （b-value）, fractal dimension （D-value）, and seismic events clustering,are utilized to interpret numerical results. A non-parametric clustering algorithm designed specificallyfor FDEM is used to reduce the mesh dependency and extract more realistic seismic information.Simulation results indicated that at the local scale, the HF process tends to propagate following the rockmass discontinuities; while at the reservoir scale, it tends to develop in the direction parallel to themaximum in-situ stress. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.展开更多
In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibi...In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibility and computational efficiency of wavelet multi-resolution method with easy implementation of the finite-difference method. The orthogonal wavelet basis provides a natural framework, which adapt spatial grids to local wavefield properties. Numerical results show usefulness of the approach as an accurate and stable tool for simulation of wave propagation in fluid-saturated porous media.展开更多
基金Supported by the National Natural Science Foundation of China under Grant No.51975138the High-Tech Ship Scientific Research Project from the Ministry of Industry and Information Technology under Grant No.CJ05N20the National Defense Basic Research Project under Grant No.JCKY2023604C006.
文摘Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The traditional thermal elastic-plastic finite element method(TEP-FEM)can accurately predict welding deformation.However,its efficiency is low because of the complex nonlinear transient computation,making it difficult to meet the needs of rapid engineering evaluation.To address this challenge,this study proposes an efficient prediction method for welding deformation in marine thin plate butt welds.This method is based on the coupled temperature gradient-thermal strain method(TG-TSM)that integrates inherent strain theory with a shell element finite element model.The proposed method first extracts the distribution pattern and characteristic value of welding-induced inherent strain through TEP-FEM analysis.This strain is then converted into the equivalent thermal load applied to the shell element model for rapid computation.The proposed method-particularly,the gradual temperature gradient-thermal strain method(GTG-TSM)-achieved improved computational efficiency and consistent precision.Furthermore,the proposed method required much less computation time than the traditional TEP-FEM.Thus,this study lays the foundation for future prediction of welding deformation in more complex marine thin plates.
基金supported by the National Key Research and Development Program of China(2022YFC2904400)Guangxi Science and Technology Major Project(Gui Ke AA23023033)。
文摘As a pyrometallurgical process,circulating fluidized bed(CFB) roasting has good potential for application in desulfurization of high-sulfur bauxite.The gas-solid distribution and reaction during CFB roasting of high-sulfur bauxite were simulated using the computational particle fluid dynamics(CPFD) method.The effect of primary air flow velocity on particle velocity,particle volume distribution,furnace temperature distribution and pressure distribution were investigated.Under the condition of the same total flow of natural gas,the impact of the number of inlets on the desulfurization efficiency,atmosphere mass fraction distribution and temperature distribution in the furnace was further investigated.
基金funded by the 14th Five-Year Plan Major Science and Technology Project of CNOOC project number KJGG2021-0506.
文摘After a long period of water flooding development,the oilfield has entered the middle and high water cut stage.The physical properties of reservoirs are changed by water erosion,which directly impacts reservoir development.Conventional numerical reservoir simulation methodologies typically employ static assumptions for model construction,presuming invariant reservoir geological parameters throughout the development process while neglecting the reservoir’s temporal evolution characteristics.Although such simplifications reduce computational complexity,they introduce substantial descriptive inaccuracies.Therefore,this paper proposes a meshless numerical simulation method for reservoirs that considers time-varying characteristics.This method avoids the meshing in traditional numerical simulation methods.From the fluid flow perspective,the reservoir’s computational domain is discretized into a series of connection units.An influence domain with a certain radius centered on the nodes is selected,and one-dimensional connection units are established between the nodes to achieve the characterization of the flow topology structure of the reservoir.In order to reflect the dynamic evolution of the reservoir’s physical properties during the water injection development process,the time-varying characteristics are incorporated into the formula of the seepage characteristic parameters in the meshless calculation.The change relationship of the permeability under different surface fluxes is considered to update the calculated connection conductivity in real time.By combining with the seepage control equation for solution,a time-varying meshless numerical simulation method is formed.The results show that compared with the numerical simulationmethod of the connection elementmethod(CEM)that only considers static parameters,this method has higher simulation accuracy and can better simulate the real migration and distribution of oil and water in the reservoir.Thismethod improves the accuracy of reservoir numerical simulation and the development effect of oilfields,providing a scientific basis for optimizing the water injection strategy,adjusting the production plan,and extending the effective production cycle of the oilfield.
基金the financial support from National Key R&D Program of China(Grant number:2024YFC2815100)Natural Science Foundation of China(Grant number:52322110)Beijing Nova Program(Grant number:20230484341).
文摘Deepwater drilling riser is the key equipment connecting the subsea wellhead and floating drilling platform.Due to complex marine environment,vortex-induced vibration(ViV)will be generated on riser,which will induce fatigue failure and even cause unpredictable drilling accidents.Therefore,it is important to study the ViV characteristics of deepwater drilling riser and reveal the main controlling factors for ensuring the safe and efficient operation of deepwater drilling engineering.In this paper,the ViV of deepwater drilling riser is numerically simulated in time domain based on the discrete vortex method(DvM).A hydrodynamic analysis model and governing equation of VIV is proposed with considering the effect of riser motion using DVM and slice method,where the governing equation is solved by Runge-Kutta method.Model validation is performed,which verified the correctness and accuracy of the mechanical model and the solution method.On this basis,the influence of the number of control points,current velocity,riser outer diameter,shear flow and top tension on the ViV characteristics of deepwater drilling risers are discussed in detail.The results show that with the increase of current velocity,the vibration amplitude of deepwater drilling riser decreases obviously,while the vibration frequency increases gradually.However,if the outer diameter of riser increases,the vibration amplitude increases,while the vibration frequency decreases gradually.The top tension also has great influence on the VIV of riser.When the top tension is 1.25 G,the VIV is suppressed to a certain extent.This study has guiding significance for optimal design and engineering control of deepwater drilling riser.
基金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 National Natural Science Foundation of China(No.51378121)the Fok Ying Tung Education Foundation(No.141076)the Scientific Innovation Research of College Graduates in Jiangsu Province(No.KYLX_0164)
文摘In order to predict the long-term rutting of asphalt pavement, the effective temperature for pavement rutting is calculated using the numerical simulation method. The transient temperature field of asphalt pavement was simulated based on actual meteorological data of Nanjing. 24-hour rutting development under a transient temperature field was calculated in each month. The rutting depth accumulated under the static temperature field was also estimated and the relationship between constant temperature parameters was analyzed. Then the effective temperature for pavement rutting was determined based on the rutting equivalence principle. The results show that the monthly effective temperature is above 40 t in July and August, while in June and September it ranges from 30 to 40 Rutting development can be ignored when the monthly effective temperature is less than 30 t. The yearly effective temperature for rutting in Nanjing is around 38. 5 t. The long-term rutting prediction model based on the effective temperature can reflect the influences of meteorological factors and traffic time distribution.
文摘LS-DYNA program and the principle of ALE method were introduced, and the target features of the reinforced concrete penetration were analyzed by using the D material model and the ALE method. A numerical simulation has been done to show the penetration visually and veritably. The simulation results are analyzed carefully and explicitly prove their significance to the research of reinforced concrete penetration.
文摘Shale gas production involves complex gas-water two-phase flow,with flow patterns in proppant-filled fractures playing a critical role in determining production efficiency.In this study,3D geometric models of 40/70 mesh ceramic particles and quartz sand proppant clusters were elaborated using computed tomography(CT)scanning.These models were used to develop a numerical simulation framework based on the lattice Boltzmann method(LBM),enabling the investigation of gas-water flow behavior within proppant-filled fractures under varying driving forces and surface tensions.Simulation results at a closure pressure of 15 MPa have revealed that ceramic particles exhibit a simpler and more porous internal structure than quartz sand of the same size.Under identical flow conditions,ceramic proppants demonstrate higher fluid replacement efficiency.Replacement efficiency increases with higher porosity,greater driving force,and lower surface tension.Furthermore,fluid displacement is strongly influenced by pore geometry:flow is faster in straighter and wider channels,with preferential movement through larger pores forming dominant flow paths.The replacement velocity exhibits a characteristic time evolution,initially rapid,then gradually decreasing,correlating positively with the development of these dominant channels.
基金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.
文摘Considering the limitation of computational capacity, a new finite element solution is used to simulate the welding deformation of the side sill of railroad car＇ s bogie frame based on the local-global method. Firstly, a volumetric heat source defined by a double ellipsoid is adopted to simulate the thermal distributions of the arc welding process. And then, the local models extracted from the global model are computed with refined meshes. On these bases, the global distortions of the subject studied are ascertained by transferring the inner forces of computed local models to the global model. It indicates that the local-global method is feasible for simulating the large welded structures by comparing the computed results with the corresponding actual measured values. The work provides basis for optimizing the welding sequence and clamping conditions, and has theoretical values and engineering significance in the integral design, manufacturing technique selection of the bogie frame, as well as other kinds of large welded structures.
文摘The mechanism of stress generation and propagation by detonation loading in five separate independent advance of ore breaking patterns is discussed in the paper. An elastic numerical model was developed using AN- SYS/LS-DYNA 3D Nonlinear Dynamic Finite Element Software. In this package ANSYS is the preprocessor and LS-DYNA is the postprocessor. Numerical models in the paper to actual were l：10 and the element mesh was dissected in scanning mode utilizing the symmetry characteristics of the numerical model. Five different advance rates were studied. Parameters, such as the time required to maximum stress, the action time of the available stress, the maximum velocity of the nodes, the stress penetration time, the magnitude of the stress peak and the time duration for high stress were numerically simulated. The 2.2 m advance appeared optimum from an analysis of the simulation results. The results from numerical simulation have been validated by tests with physical models.
基金National Basic Research Program of China Under Grant No. 2007CB714200National Natural Science Foundation of China Under Grant No. 90715038
文摘In this paper, a method to develop a hierarchy of explicit recursion formulas for numerical simulation in an irregular grid for scalar wave equations is presented and its accuracy is illustrated via 2-D and 1-D models. Approaches to develop the stable formulas which are of 2M-order accuracy in both time and space with Mbeing a positive integer for regular grids are discussed and illustrated by constructing the second order （M= 1） and the fourth order （M = 2） recursion formulas.
文摘Non-pillar continuous mining(NPCM) is regarded as a high-efficient, high-level and one-step mining technology, which can be divided into two substopes. Back fill stability status in substope I, which directly influence the loss rate and dilution rate, etc, will determine whether the experimental research is successful or not. By employing energy method of limit analysis and finite element numerical simulation method, the critical backfill height was determined under the prerequisite condition of its stability, which put forward theoretical basis for reasonable and correct selection of backfill’s parameters. The result showed that the first backfill could not keep stable for NPCM, while the other was able to.
基金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.
基金sponsored by the project No.50404001 from the National Natural Science Foundation of Chinathe National Key Fundamental Research & Development Project(Grant No.2007CB209601)+1 种基金the China National PetroleumCorporation Fundamental Research Program (Grant No.06A30102)the China Postdoctoral Science Foundation(Project No.2004035350)
文摘The electrophysical property of saturated rocks is very important for reservoir identification and evaluation. In this paper, the lattice Boltzmann method （LBM） was used to study the electrophysical property of rock saturated with fluid because of its advantages over conventional numerical approaches in handling complex pore geometry and boundary conditions. The digital core model was constructed through the accumulation of matrix grains based on their radius distribution obtained by the measurements of core samples. The flow of electrical current through the core model saturated with oil and water was simulated on the mesoscopic scale to reveal the non-Archie relationship between resistivity index and water saturation （I-Sw）. The results from LBM simulation and laboratory measurements demonstrated that the I-Sw relation in the range of low water saturation was generally not a straight line in the log-log coordinates as described by the Archie equation. We thus developed a new equation based on numerical simulation and physical experiments. This new equation was used to fit the data from laboratory core measurements and previously published data. Determination of fluid saturation and reservoir evaluation could be significantly improved by using the new equation.
基金financially supported by the National Natural Science Foundation of China(No.51571064)the National Basic Research Program of China(No.2014CB643701)。
文摘The numerical simulation model of single roller melt spinning for rapid quenching process of NdFeB alloy was built,and the vacuum chamber,cooling roller and sample were taken into account as a system.The existing mature technology was in order to verify the correctness of simulation.The rapid quenching ribbons with different roll speeds were used as the simulation objects.The results of the numerical simulation and experiments show that the validity of the model has been testified and the reasons of the formation of complete quenching ribbons and by-product have been explained.The experimental thickness of the ribbons is proportional to the theoretical thickness.In the same spray condition,with the roll speed increasing,the thickness decreases linearly.At the speed range of25-30 m·s^(-1),the simulated calculation date is close to the experimental date,which can be considered as an ideal technological parameter.
基金The National Natural Science Foundation of China under contract No.51609101 and 51909103the Natural Science Foundation of Fujian Province of China under contract Nos 2017J01701,2017J05085 and 2018J05090the Outstanding Young University Scientific Research Talents Cultivation Plan of Fujian Province of China
文摘Numerical wave tanks are widely-acknowledged tools in studying waves and wave-structure interactions. They can generate waves under realistic scales and offers more information on the fluid field. However, most numerical wave tanks suffer from issues known as the numerical dissipation and numerical dispersion. The former causes wave energy to be slowly dissipated and the latter shifts wave frequencies during wave propagation. This paper proposes a simple method of depressing numerical dissipation effects on the basis of solving Euler equations using the finite difference method(FDM). The wave propagation solutions are solved analytically taking into account the influence of the damping terms. The main idea of the method is to append a source term to the momentum equation, whose strength is determined by how strong the numerical damping effect is. The method is verified by successfully depressing numerical effects during the simulation of regular linear waves, Stokes waves and irregular waves. By applying the method, wave energy is able to be close to its initial value after long distance of travel.
基金supported by the National Key Technologies R&D Program (Grant No. 2018YFB1106000)Innovation Funds of China Academy of Launch Vehicle Technology (CALT) for Universities (Grant No.CALT201709)Tsinghua Grants for Autonomous Research。
文摘Wire arc additive manufacture(WAAM) is a new technique to fabricate large-scale complex aluminum alloy components.However, the performance of the parts is critically influenced by residual stresses and deformation. A sequentially thermal-mechanical coupled model of residual stress and deformation for aluminum alloy WAAM parts was established based on commercial FE software ABAQUS. The temperature field was calculated by the moving heat source(MHS) method. The temperature function was obtained according to the distribution of the peak temperature. Furthermore, the MHS method and segmented temperature function(STF) method were used to calculate the residual stress and deformation. The results show that the STF method satisfies both the efficiency and accuracy requirements. 1-segment, 3-segment, and 5-segment STF methods can shorten the time for mechanical analysis by 91%, 79%, 63%, respectively.The error of the residual stress and deformation are all less than 20%. STF method provides an effective way to predict the residual stress and deformation of large-scale WAAM parts.
基金supported by the Natural Sciences and Engineering Research Council of Canada through Discovery Grant 341275 (G. Grasselli) and Engage EGP 461019-13
文摘Hydraulic fracturing （HF） technique has been extensively used for the exploitation of unconventional oiland gas reservoirs. HF enhances the connectivity of less permeable oil and gas-bearing rock formationsby fluid injection, which creates an interconnected fracture network and increases the hydrocarbonproduction. Meanwhile, microseismic （MS） monitoring is one of the most effective approaches to evaluatesuch stimulation process. In this paper, the combined finite-discrete element method （FDEM） isadopted to numerically simulate HF and associated MS. Several post-processing tools, includingfrequency-magnitude distribution （b-value）, fractal dimension （D-value）, and seismic events clustering,are utilized to interpret numerical results. A non-parametric clustering algorithm designed specificallyfor FDEM is used to reduce the mesh dependency and extract more realistic seismic information.Simulation results indicated that at the local scale, the HF process tends to propagate following the rockmass discontinuities; while at the reservoir scale, it tends to develop in the direction parallel to themaximum in-situ stress. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.
基金the National Natural Science Foundation of China(No.40774056)Program of Excellent Team in Harbin Institute of Technology
文摘In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibility and computational efficiency of wavelet multi-resolution method with easy implementation of the finite-difference method. The orthogonal wavelet basis provides a natural framework, which adapt spatial grids to local wavefield properties. Numerical results show usefulness of the approach as an accurate and stable tool for simulation of wave propagation in fluid-saturated porous media.
