This study aims to develop a chloride diffusion simulation method that considers the hydration microstructure and pore solution properties during the hydration of tricalcium silicate(C3S).The method combines the hydra...This study aims to develop a chloride diffusion simulation method that considers the hydration microstructure and pore solution properties during the hydration of tricalcium silicate(C3S).The method combines the hydration simulation,thermodynamic calculation,and finite element analysis to examine the effects of pore solution,including effect of electrochemical potential,effect of chemical activity,and effect of mechanical interactions between ions,on the chloride effective diffusion coefficient of hydrated C3S paste.The results indicate that the effect of electrochemical potential on chloride diffusion becomes stronger with increasing hydration age due to the increase in the content of hydrated calcium silicate;as the hydration age increases,the effect of chemical activity on chloride diffusion weakens when the number of diffusible elements decreases;the effect of mechanical interactions between ions on chloride diffusion decreases with the increase of hydration age.展开更多
Without knowing the emittance value, it is difficult to optimize ion beam optics for minimum beam loss during transmission, especially considering the very high emittance values of electron cyclotron resonance(ECR) io...Without knowing the emittance value, it is difficult to optimize ion beam optics for minimum beam loss during transmission, especially considering the very high emittance values of electron cyclotron resonance(ECR) ion sources.With this in mind, to measure the emittance of the ion beams produced by the mVINIS ECR, which is part of the FAMA facility at the Vin?a Institute of Nuclear Sciences, we have developed a pepper-pot scintillator screen system combined with a CMOS camera. The application, developed on the Lab VIEW platform, allows us to control the camera's main attribute settings, such as the shutter speed and the gain, record the images in the region of interest, and process and filter the images in real time. To analyze the data from the obtained image, we have developed an algorithm called measurement and analysis of ion beam luminosity(MAIBL) to reconstruct the four-dimensional(4D) beam profile and calculate the root mean square(RMS) emittance. Before measuring emittance, we performed a simulated experiment using the pepper-pot simulation(PPS) program. An exported file(PPS) gives a numerically generated raw image(mock image) of a beam with a predefined emittance value after it has passed through a pepper-pot mask. By analyzing data from mock images instead of the image obtained by the camera and putting it into the MAIBL algorithm, we can compare the calculated emittance with PPS's initial emittance value. In this paper, we present our computational tools and explain the method for verifying the correctness of the calculated emittance values.展开更多
The initial shape of the secondary arc considerably influences its subsequent shape.To establish the model for the arcing time of the secondary arc and modify the single-phase reclosing sequence,theoretical and experi...The initial shape of the secondary arc considerably influences its subsequent shape.To establish the model for the arcing time of the secondary arc and modify the single-phase reclosing sequence,theoretical and experimental analysis of the evolution process of the short-circuit arc to the secondary arc is critical.In this study,an improved charge simulation method was used to develop the internal-space electric-field model of the short-circuit arc.The intensity of the electric field was used as an independent variable to describe the initial shape of the secondary arc.A secondary arc evolution model was developed based on this model.Moreover,the accuracy of the model was evaluated by comparison with physical experimental results.When the secondary arc current increased,the arcing time and dispersion increased.There is an overall trend of increasing arc length with increasing arcing time.Nevertheless,there is a reduction in arc length during arc ignition due to short circuits between the arc columns.Furthermore,the arcing time decreased in the range of 0°-90°as the angle between the wind direction and the x-axis increased.This work investigated the method by which short-circuit arcs evolve into secondary arcs.The results can be used to develop the secondary arc evolution model and to provide both a technical and theoretical basis for secondary arc suppression.展开更多
This article gives a mechanical model, in which the layers of lithosphere are assumed to be the creep materials, to study the coupling mechanism of a syn-basin-mountain system quantitatively by using the numerical s...This article gives a mechanical model, in which the layers of lithosphere are assumed to be the creep materials, to study the coupling mechanism of a syn-basin-mountain system quantitatively by using the numerical simulating method. A geological dynamic extensional mode given by some geologists is theoretically discussed and verified. The study shows that lithosphere thickening or thinning is closely related to the thermal activity, or in other words, thermal convection beneath the lithosphere. It is one of the important factors affecting the formation of the basin-mountain coupling system. As an essential condition, only the upward buoyant force and the horizontal dragging force caused by the thermal convection jointly act on the bottom of the lithosphere, the stress and strain states in rock's layers are advantageous to forming the tectonic-landforms of the basin-mountain coupling system. A study on the creep features of the lithosphere shows that the stress and strain in the rock's layers vary with time when the lasting forces act on the boundary. They increase rapidly at initial stage and decrease steadily after reaching the peak value. Phenomena of stress relaxation are significant for studying the tectonic evolution.展开更多
The dynamic performance of high-speed trains is significantly influenced by sudden changes in aerodynamic loads(ADLs)when exiting a tunnel in a windy environment.Focusing on a double-track tunnel under construction in...The dynamic performance of high-speed trains is significantly influenced by sudden changes in aerodynamic loads(ADLs)when exiting a tunnel in a windy environment.Focusing on a double-track tunnel under construction in a mountain railway,we established an aerodynamic model involving a train exiting the tunnel,and verified it in the Fluent environment.Overset mesh technology was adopted to characterize the train’s movement.The flow field involving the train,tunnel,and crosswinds was simulated using the Reynolds-averaged turbulence model.Then,we built a comprehensive train-track coupled dynamic model considering the influences of ADLs,to investigate the vehicles’dynamic responses.The aerodynamics and dynamic behaviors of the train when affected by crosswinds with different velocities and directions are analyzed and discussed.The results show that the near-wall side crosswind leads to sharper variations in ADLs than the far-wall side crosswind.The leading vehicle suffers from more severe ADLs than other vehicles,which worsens the wheel-rail interaction and causes low-frequency vibration of the car body.When the crosswind velocity exceeds 20 m/s,significant wheel-rail impacts occur,and the running safety of the train worsens rapidly.展开更多
The fatigue of heavy-haul railway bridges is considered a key concern due to high stress levels and cyclic loading.The evaluation of fatigue reliability is required to include factor correlations.A major challenge is ...The fatigue of heavy-haul railway bridges is considered a key concern due to high stress levels and cyclic loading.The evaluation of fatigue reliability is required to include factor correlations.A major challenge is presented by the construction of the cumulative distribution function(CDF)and the description of correlations between random variables.In this study,the copula function is used to analyze the fatigue failure probability of the Shuohuang heavy-haul railway bridge.A C-vine copula(CVC)-based joint probability density function(JPDF)is derived with eight correlated parameters.To enhance efficiency in small failure probability calculations,the subset simulation and most probable point(MPP)Monte Carlo importance sampling are introduced based on the Rosenblatt transform and C-vine model.Comparisons with traditional Monte Carlo methods confirm that high accuracy and efficiency are achieved.The results show that when parameter correlations are ignored,failure probability is underestimated,increasing safety risks in bridge assessments.展开更多
Grouting is a widely applied technique for reinforcing fractured zones in deep soft rock tunnels.By infiltrating rock fissures,slurry materials enhance structural integrity and improve the overall stability of the sur...Grouting is a widely applied technique for reinforcing fractured zones in deep soft rock tunnels.By infiltrating rock fissures,slurry materials enhance structural integrity and improve the overall stability of the surrounding rock.The performance of grouting is primarily governed by the flow behavior and diffusion extent of the slurry.This review considers recent advances in the theory and methodology of slurry flow and diffusion in fractured rock.It examines commonly used grout materials,including cement-based,chemical,and composite formulations,each offering distinct advantages for specific geological conditions.Themechanisms of reinforcement vary significantly across materials,requiring tailored application strategies.The rheological properties of grouting slurries,particularly cement-based types,have been widely modeled using classical constitutive approaches.However,the influence of time-and space-dependent viscosity evolution on slurry behavior remains underexplored.Experimental studies have provided valuable insights into slurry diffusion,yet further research is needed to capture real-time behavior under multi-scale and multi-physics coupling conditions.Similarly,current numerical simulations are largely limited to twoand three-dimensional models of single-fracture flow.These models often neglect the complexity of fracture networks and geological heterogeneity,highlighting a need for more realistic and integrated simulation frameworks.Future research should focus on:(1)fine-scale modeling of slurry hydration and mechanical reinforcement processes;(2)cross-scale analysis of slurry flow under coupled thermal,hydraulic,andmechanical fields;and(3)development of realtime,three-dimensional dynamic simulation tools to capture the full grouting process.These efforts will strengthen the theoretical foundation and practical effectiveness of grouting in complex underground environments.展开更多
A 3-D digital core describes the pore space microstructure of rocks. An X-ray micro CT scan is the most accurate and direct but costly method to obtain a 3-D digital core. In this study, we propose a hybrid method whi...A 3-D digital core describes the pore space microstructure of rocks. An X-ray micro CT scan is the most accurate and direct but costly method to obtain a 3-D digital core. In this study, we propose a hybrid method which combines sedimentation simulation and simulated annealing （SA） method to generate 3-D digital cores based on 2-D images of rocks. The method starts with the sedimentation simulation to build a 3-D digital core, which is the initial configuration for the SA method. We update the initial digital core using the SA method to match the auto-correlation function of the 2-D rock image and eventually build the final 3-D digital core. Compared with the typical SA method, the hybrid method has significantly reduced the computation time. Local porosity theory is applied to quantitatively compare the reconstructed 3-D digital cores with the X-ray micro CT 3-D images. The results indicate that the 3-D digital cores reconstructed by the hybrid method have homogeneity and geometric connectivity similar to those of the X-ray micro CT image. The formation factors and permeabilities of the reconstructed 3-D digital cores are estimated using the finite element method （FEM） and lattice Boltzmann method （LBM）, respectively. The simulated results are in good agreement with the experimental measurements. Comparison of the simulation results suggests that the digital cores reconstructed by the hybrid method more closely reflect the true transport properties than the typical SA method alone.展开更多
It is difficult to temporally and spatially track and characterize the slurry viscosity in flowing water during grouting simulation.In this study,a sequential flow and solidification(SFS)method considering the spatial...It is difficult to temporally and spatially track and characterize the slurry viscosity in flowing water during grouting simulation.In this study,a sequential flow and solidification(SFS)method considering the spatial-temporal evolution of slurry viscosity in flowing water in karst conduit is proposed.First,a time-dependent model for the threshold function of slurry viscosity is established.During the grouting process,the spatial-temporal evolution of slurry viscosity is revealed by tracking the diffusion behavior of the slurry injected at different times.This method is capable of describing the gradual solidification process of the slurry during grouting.Furthermore,a physical model of grouting in a karst conduit is developed.Second,the effectiveness of the SFS method in grouting simulation is verified by the experiment of grouting conduit in flowing water.The SFS method enables real-time monitoring of fluid velocity and pressure during grouting in flowing water and provides a feasible calculation method for revealing the grouting plugging mechanism in complex karst conduits at different engineering scales.In addition,it can be used to guide the design of grouting tests in flowing water,improve cost efficiency,and provide theoretical basis for optimizing grouting design and slurry selection.展开更多
The particle simulation method is used to study the effects of loading waveforms （i.e. square, sinusoidal and triangle waveforms） on rock damage at mesoscopic scale. Then some influencing factors on rock damage at t...The particle simulation method is used to study the effects of loading waveforms （i.e. square, sinusoidal and triangle waveforms） on rock damage at mesoscopic scale. Then some influencing factors on rock damage at the mesoscopic scale, such as loading frequency, stress amplitude, mean stress, confining pressure and loading sequence, are also investigated with sinusoidal waveform in detail. The related numerical results have demonstrated that： 1） the loading waveform has a certain effect on rock failure processes. The square waveform has the most damage within these waveforms, while the triangle waveform has less damage than sinusoidal waveform. In each cycle, the number of microscopic cracks increases in the loading stage, while it keeps nearly constant in the unloading stage. 2） The loading frequency, stress amplitude, mean stress, confining pressure and loading sequence have considerable effects on rock damage subjected to cyclic loading. The higher the loading frequency, stress amplitude and mean stress, the greater the damage the rock accumulated; in contrast, the lower the confining pressure, the greater the damage the rock has accumulated. 3） There is a threshold value of mean stress and stress amplitude, below which no further damage accumulated after the first few cycle loadings. 4） The high-to-low loading sequence has more damage than the low-to-high loading sequence, suggesting that the rock damage is loading-path dependent.展开更多
The control rod drive mechanism(CRDM)is an essential part of the control and safety protection system of pressurized water reactors.Current CRDM simulations are mostly performed collectively using a single method,igno...The control rod drive mechanism(CRDM)is an essential part of the control and safety protection system of pressurized water reactors.Current CRDM simulations are mostly performed collectively using a single method,ignoring the influence of multiple motion units and the differences in various features among them,which strongly affect the efficiency and accuracy of the simulations.In this study,we constructed a flow field fusion simulation method based on model features by combining key motion unit analysis and various simulation methods and then applied the method to the CRDM simulation process.CRDM performs motion unit decomposition through the structural hierarchy of function-movement-action method,and the key meta-actions are identified as the nodes in the flow field simulation.We established a fused feature-based multimethod simulation process and processed the simulation methods and data according to the features of the fluid domain space and the structural complexity to obtain the fusion simulation results.Compared to traditional simulation methods and real measurements,the simulation method provides advantages in terms of simulation efficiency and accuracy.展开更多
A passive simulation method based on the six degrees of freedom(6-DOF)model and dynamic mesh is proposed according to the working principle to study the dynamic characteristics of the turbine flow sensors.This simulat...A passive simulation method based on the six degrees of freedom(6-DOF)model and dynamic mesh is proposed according to the working principle to study the dynamic characteristics of the turbine flow sensors.This simulation method controls the six degrees of freedom of the impeller using the user-defined functions(UDF)program so that it can only rotate under the impact of fluid.The impeller speed can be calculated in real-time,and the inlet speed can be set with time to obtain the dynamic performance of the turbine flow sensors.Based on this simulation method,three turbine flow sensors with different diameters were simulated,and the reliability of the simulation method was verified by both steady-state and unsteady-state experiments.The results show that the trend of meter factor with flow rate acquired from the simulation is close to the experimental results.The deviation between the simulation and experiment results is low,with a maximum deviation of 2.88%.In the unsteady simulation study,the impeller speed changed with the inlet velocity of the turbine flow sensor,showing good tracking performance.The passive simulation method can be used to predict the dynamic performance of the turbine flow sensor.展开更多
The investigation was carried out on the technical problems of finishing the inner surface of elbow parts and the action mechanism of particles in elbow precision machining by abrasive flow.This work was analyzed and ...The investigation was carried out on the technical problems of finishing the inner surface of elbow parts and the action mechanism of particles in elbow precision machining by abrasive flow.This work was analyzed and researched by combining theory,numerical and experimental methods.The direct simulation Monte Carlo(DSMC)method and the finite element analysis method were combined to reveal the random collision of particles during the precision machining of abrasive flow.Under different inlet velocity,volume fraction and abrasive particle size,the dynamic pressure and turbulence flow energy of abrasive flow in elbow were analyzed,and the machining mechanism of particles on the wall and the influence of different machining parameters on the precision machining quality of abrasive flow were obtained.The test results show the order of the influence of different parameters on the quality of abrasive flow precision machining and establish the optimal process parameters.The results of the surface morphology before and after the precision machining of the inner surface of the elbow are discussed,and the surface roughness Ra value is reduced from 1.125μm to 0.295μm after the precision machining of the abrasive flow.The application of DSMC method provides special insights for the development of abrasive flow technology.展开更多
To overcome the default of single search tendency, the ants in the colony are divided into several sub-groups. The ants in different subgroups have different trail information and expectation coefficients. The simulat...To overcome the default of single search tendency, the ants in the colony are divided into several sub-groups. The ants in different subgroups have different trail information and expectation coefficients. The simulated annealing method is introduced to the algorithm. Through setting the temperature changing with the iterations, after each turn of tours, the solution set obtained by the ants is taken as the candidate set. The update set is obtained by adding the solutions in the candidate set to the previous update set with the probability determined by the temperature. The solutions in the candidate set are used to update the trail information. In each turn of updating, the current best solution is also used to enhance the trail information on the current best route. The trail information is reset when the algorithm is in stagnation state. The computer experiments demonstrate that the proposed algorithm has higher stability and convergence speed.展开更多
One of the most critical and complicated steps in mine design is a selection of suitable mining method based upon geological,geotechnical,geographical,safety and economical parameters.The aim of this study is developi...One of the most critical and complicated steps in mine design is a selection of suitable mining method based upon geological,geotechnical,geographical,safety and economical parameters.The aim of this study is developing a Monte Carlo simulation to selection the optimum mining method by using effective and major criteria and at the same time,taking subjective judgments of decision makers into consideration.Proposed approach is based on the combination of Monte Carlo simulation with conventional Analytic Hierarchy Process(AHP).Monte Carlo simulation is used to determine the confdence level of each alternative’s score,is calculated by AHP,with the respect to the variance of decision makers’opinion.The proposed method is applied for Jajarm Bauxite Mine in Iran and eventually the most appropriate mining methods for this mine are ranked.展开更多
This paper presents a rapid and simple risk calculation method for large and complex engineering systems, the simulated maximum entropy method （SMEM）, which is based on integration of the advantages of the Monte Car...This paper presents a rapid and simple risk calculation method for large and complex engineering systems, the simulated maximum entropy method （SMEM）, which is based on integration of the advantages of the Monte Carlo and maximum entropy methods, thus avoiding the shortcoming of the slow convergence rate of the Monte Carlo method in risk calculation. Application of SMEM in the calculation of reservoir flood discharge risk shows that this method can make full use of the known information under the same conditions and obtain the corresponding probability distribution and the risk value. It not only greatly improves the speed, compared with the Monte Carlo method, but also provides a new approach for the risk calculation in large and complex engineering systems.展开更多
In this paper .fatigue strength of T-type tubularjoints subjected to in-planebending or out-of-plane bending load is investigated By, considering material constantsand initial crack sizes as random variables and appl...In this paper .fatigue strength of T-type tubularjoints subjected to in-planebending or out-of-plane bending load is investigated By, considering material constantsand initial crack sizes as random variables and applying Monte Carlo simulation method, we have giyen a statistical analvsis on fatigue life. Simultaneously, linearregression analyses of computed resuldts are performed and compared with that of theknown experimental data.展开更多
For the sake of a more accurate shell boundary and calculation of radiation heat transfer in the Directional Solidification(DS) process, a radiation heat transfer model based on the Finite Element Method(FEM)is develo...For the sake of a more accurate shell boundary and calculation of radiation heat transfer in the Directional Solidification(DS) process, a radiation heat transfer model based on the Finite Element Method(FEM)is developed in this study. Key technologies, such as distinguishing boundaries automatically, local matrix and lumped heat capacity matrix, are also stated. In order to analyze the effect of withdrawing rate on DS process,the solidification processes of a complex superalloy turbine blade in the High Rate Solidification(HRS) process with different withdrawing rates are simulated; and by comparing the simulation results, it is found that the most suitable withdrawing rate is determined to be 5.0 mm·min^(-1). Finally, the accuracy and reliability of the radiation heat transfer model are verified, because of the accordance of simulation results with practical process.展开更多
Each joint of hydraulic drive quadruped robot is driven by the hydraulic drive unit（HDU）, and the contacting between the robot foot end and the ground is complex and variable, which increases the difficulty of force...Each joint of hydraulic drive quadruped robot is driven by the hydraulic drive unit（HDU）, and the contacting between the robot foot end and the ground is complex and variable, which increases the difficulty of force control inevitably. In the recent years, although many scholars researched some control methods such as disturbance rejection control, parameter self-adaptive control, impedance control and so on, to improve the force control performance of HDU, the robustness of the force control still needs improving. Therefore, how to simulate the complex and variable load characteristics of the environment structure and how to ensure HDU having excellent force control performance with the complex and variable load characteristics are key issues to be solved in this paper. The force control system mathematic model of HDU is established by the mechanism modeling method, and the theoretical models of a novel force control compensation method and a load characteristics simulation method under different environment structures are derived, considering the dynamic characteristics of the load stiffness and the load damping under different environment structures. Then, simulation effects of the variable load stiffness and load damping under the step and sinusoidal load force are analyzed experimentally on the HDU force control performance test platform, which provides the foundation for the force control compensation experiment research. In addition, the optimized PID control parameters are designed to make the HDU have better force control performance with suitable load stiffness and load damping, under which the force control compensation method is introduced, and the robustness of the force control system with several constant load characteristics and the variable load characteristics respectively are comparatively analyzed by experiment. The research results indicate that if the load characteristics are known, the force control compensation method presented in this paper has positive compensation effects on the load characteristics variation, i.e., this method decreases the effects of the load characteristics variation on the force control performance and enhances the force control system robustness with the constant PID parameters, thereby, the online PID parameters tuning control method which is complex needs not be adopted. All the above research provides theoretical and experimental foundation for the force control method of the quadruped robot joints with high robustness.展开更多
An element-free Galerkin method(EFGM) is used to solve the two-dimensional(2D) ground penetrating radar(GPR)modelling problems, due to its simple pre-processing, the absence of elements and high accuracy. Different fr...An element-free Galerkin method(EFGM) is used to solve the two-dimensional(2D) ground penetrating radar(GPR)modelling problems, due to its simple pre-processing, the absence of elements and high accuracy. Different from element-based numerical methods, this approach makes nodes free from the elemental restraint and avoids the explicit mesh discretization. First, we derived the boundary value problem for the 2D GPR simulation problems. Second, a penalty function approach and a boundary condition truncated method were used to enforce the essential and the absorbing boundary conditions, respectively. A three-layered GPR model was used to verify our element-free approach. The numerical solutions show that our solutions have an excellent agreement with solutions of a finite element method(FEM). Then, we used the EFGM to simulate one more complex model to show its capability and limitations. Simulation results show that one obvious advantage of EFGM is the absence of element mesh, which makes the method very flexible. Due to the use of MLS fitting, a key feature of EFM, is that both the dependent variable and its gradient are continuous and have high precision.展开更多
基金Funded by the Natural Science Foundation of Jiangsu Province(No.BK20241529)China Postdoctoral Science Foundation(No.2024M750736)。
文摘This study aims to develop a chloride diffusion simulation method that considers the hydration microstructure and pore solution properties during the hydration of tricalcium silicate(C3S).The method combines the hydration simulation,thermodynamic calculation,and finite element analysis to examine the effects of pore solution,including effect of electrochemical potential,effect of chemical activity,and effect of mechanical interactions between ions,on the chloride effective diffusion coefficient of hydrated C3S paste.The results indicate that the effect of electrochemical potential on chloride diffusion becomes stronger with increasing hydration age due to the increase in the content of hydrated calcium silicate;as the hydration age increases,the effect of chemical activity on chloride diffusion weakens when the number of diffusible elements decreases;the effect of mechanical interactions between ions on chloride diffusion decreases with the increase of hydration age.
基金funded by the Ministry of Education, Science and Technological Development of the Republic of Serbia via the FAMA project (research topic “Physics and chemistry with ion beams”)。
文摘Without knowing the emittance value, it is difficult to optimize ion beam optics for minimum beam loss during transmission, especially considering the very high emittance values of electron cyclotron resonance(ECR) ion sources.With this in mind, to measure the emittance of the ion beams produced by the mVINIS ECR, which is part of the FAMA facility at the Vin?a Institute of Nuclear Sciences, we have developed a pepper-pot scintillator screen system combined with a CMOS camera. The application, developed on the Lab VIEW platform, allows us to control the camera's main attribute settings, such as the shutter speed and the gain, record the images in the region of interest, and process and filter the images in real time. To analyze the data from the obtained image, we have developed an algorithm called measurement and analysis of ion beam luminosity(MAIBL) to reconstruct the four-dimensional(4D) beam profile and calculate the root mean square(RMS) emittance. Before measuring emittance, we performed a simulated experiment using the pepper-pot simulation(PPS) program. An exported file(PPS) gives a numerically generated raw image(mock image) of a beam with a predefined emittance value after it has passed through a pepper-pot mask. By analyzing data from mock images instead of the image obtained by the camera and putting it into the MAIBL algorithm, we can compare the calculated emittance with PPS's initial emittance value. In this paper, we present our computational tools and explain the method for verifying the correctness of the calculated emittance values.
基金supported by National Natural Science Foundation of China(Nos.92066108 and 51277061)。
文摘The initial shape of the secondary arc considerably influences its subsequent shape.To establish the model for the arcing time of the secondary arc and modify the single-phase reclosing sequence,theoretical and experimental analysis of the evolution process of the short-circuit arc to the secondary arc is critical.In this study,an improved charge simulation method was used to develop the internal-space electric-field model of the short-circuit arc.The intensity of the electric field was used as an independent variable to describe the initial shape of the secondary arc.A secondary arc evolution model was developed based on this model.Moreover,the accuracy of the model was evaluated by comparison with physical experimental results.When the secondary arc current increased,the arcing time and dispersion increased.There is an overall trend of increasing arc length with increasing arcing time.Nevertheless,there is a reduction in arc length during arc ignition due to short circuits between the arc columns.Furthermore,the arcing time decreased in the range of 0°-90°as the angle between the wind direction and the x-axis increased.This work investigated the method by which short-circuit arcs evolve into secondary arcs.The results can be used to develop the secondary arc evolution model and to provide both a technical and theoretical basis for secondary arc suppression.
基金ThispaperissupportedbytheNationalNaturalScienceFoundationofChina (No .1 9972 0 64)andtheMinistryofLandandResources (No.991 0 0 5)
文摘This article gives a mechanical model, in which the layers of lithosphere are assumed to be the creep materials, to study the coupling mechanism of a syn-basin-mountain system quantitatively by using the numerical simulating method. A geological dynamic extensional mode given by some geologists is theoretically discussed and verified. The study shows that lithosphere thickening or thinning is closely related to the thermal activity, or in other words, thermal convection beneath the lithosphere. It is one of the important factors affecting the formation of the basin-mountain coupling system. As an essential condition, only the upward buoyant force and the horizontal dragging force caused by the thermal convection jointly act on the bottom of the lithosphere, the stress and strain states in rock's layers are advantageous to forming the tectonic-landforms of the basin-mountain coupling system. A study on the creep features of the lithosphere shows that the stress and strain in the rock's layers vary with time when the lasting forces act on the boundary. They increase rapidly at initial stage and decrease steadily after reaching the peak value. Phenomena of stress relaxation are significant for studying the tectonic evolution.
基金National Natural Science Foundation of China(No.52388102)New Cornerstone Science Foundation through the Xplorer Prize.
文摘The dynamic performance of high-speed trains is significantly influenced by sudden changes in aerodynamic loads(ADLs)when exiting a tunnel in a windy environment.Focusing on a double-track tunnel under construction in a mountain railway,we established an aerodynamic model involving a train exiting the tunnel,and verified it in the Fluent environment.Overset mesh technology was adopted to characterize the train’s movement.The flow field involving the train,tunnel,and crosswinds was simulated using the Reynolds-averaged turbulence model.Then,we built a comprehensive train-track coupled dynamic model considering the influences of ADLs,to investigate the vehicles’dynamic responses.The aerodynamics and dynamic behaviors of the train when affected by crosswinds with different velocities and directions are analyzed and discussed.The results show that the near-wall side crosswind leads to sharper variations in ADLs than the far-wall side crosswind.The leading vehicle suffers from more severe ADLs than other vehicles,which worsens the wheel-rail interaction and causes low-frequency vibration of the car body.When the crosswind velocity exceeds 20 m/s,significant wheel-rail impacts occur,and the running safety of the train worsens rapidly.
基金Project supported by the National Natural Science Foundation of China(No.52278180)。
文摘The fatigue of heavy-haul railway bridges is considered a key concern due to high stress levels and cyclic loading.The evaluation of fatigue reliability is required to include factor correlations.A major challenge is presented by the construction of the cumulative distribution function(CDF)and the description of correlations between random variables.In this study,the copula function is used to analyze the fatigue failure probability of the Shuohuang heavy-haul railway bridge.A C-vine copula(CVC)-based joint probability density function(JPDF)is derived with eight correlated parameters.To enhance efficiency in small failure probability calculations,the subset simulation and most probable point(MPP)Monte Carlo importance sampling are introduced based on the Rosenblatt transform and C-vine model.Comparisons with traditional Monte Carlo methods confirm that high accuracy and efficiency are achieved.The results show that when parameter correlations are ignored,failure probability is underestimated,increasing safety risks in bridge assessments.
基金funded by the National Natural Science Foundation of China[U22A20234]Hubei Province key research and development project[2023BCB121]Wuhan innovation supporting projects[2023020201010079].
文摘Grouting is a widely applied technique for reinforcing fractured zones in deep soft rock tunnels.By infiltrating rock fissures,slurry materials enhance structural integrity and improve the overall stability of the surrounding rock.The performance of grouting is primarily governed by the flow behavior and diffusion extent of the slurry.This review considers recent advances in the theory and methodology of slurry flow and diffusion in fractured rock.It examines commonly used grout materials,including cement-based,chemical,and composite formulations,each offering distinct advantages for specific geological conditions.Themechanisms of reinforcement vary significantly across materials,requiring tailored application strategies.The rheological properties of grouting slurries,particularly cement-based types,have been widely modeled using classical constitutive approaches.However,the influence of time-and space-dependent viscosity evolution on slurry behavior remains underexplored.Experimental studies have provided valuable insights into slurry diffusion,yet further research is needed to capture real-time behavior under multi-scale and multi-physics coupling conditions.Similarly,current numerical simulations are largely limited to twoand three-dimensional models of single-fracture flow.These models often neglect the complexity of fracture networks and geological heterogeneity,highlighting a need for more realistic and integrated simulation frameworks.Future research should focus on:(1)fine-scale modeling of slurry hydration and mechanical reinforcement processes;(2)cross-scale analysis of slurry flow under coupled thermal,hydraulic,andmechanical fields;and(3)development of realtime,three-dimensional dynamic simulation tools to capture the full grouting process.These efforts will strengthen the theoretical foundation and practical effectiveness of grouting in complex underground environments.
基金sponsored by NSFC(Grant No.40574030)CNPC Research Project(Grant No.06A30102)
文摘A 3-D digital core describes the pore space microstructure of rocks. An X-ray micro CT scan is the most accurate and direct but costly method to obtain a 3-D digital core. In this study, we propose a hybrid method which combines sedimentation simulation and simulated annealing （SA） method to generate 3-D digital cores based on 2-D images of rocks. The method starts with the sedimentation simulation to build a 3-D digital core, which is the initial configuration for the SA method. We update the initial digital core using the SA method to match the auto-correlation function of the 2-D rock image and eventually build the final 3-D digital core. Compared with the typical SA method, the hybrid method has significantly reduced the computation time. Local porosity theory is applied to quantitatively compare the reconstructed 3-D digital cores with the X-ray micro CT 3-D images. The results indicate that the 3-D digital cores reconstructed by the hybrid method have homogeneity and geometric connectivity similar to those of the X-ray micro CT image. The formation factors and permeabilities of the reconstructed 3-D digital cores are estimated using the finite element method （FEM） and lattice Boltzmann method （LBM）, respectively. The simulated results are in good agreement with the experimental measurements. Comparison of the simulation results suggests that the digital cores reconstructed by the hybrid method more closely reflect the true transport properties than the typical SA method alone.
基金financial support from the National Natural Science Foundation of China(Grant Nos.52022053 and 51879153)the China National Postdoctoral Program for Innovative Talents(Grant No.BX2021172)。
文摘It is difficult to temporally and spatially track and characterize the slurry viscosity in flowing water during grouting simulation.In this study,a sequential flow and solidification(SFS)method considering the spatial-temporal evolution of slurry viscosity in flowing water in karst conduit is proposed.First,a time-dependent model for the threshold function of slurry viscosity is established.During the grouting process,the spatial-temporal evolution of slurry viscosity is revealed by tracking the diffusion behavior of the slurry injected at different times.This method is capable of describing the gradual solidification process of the slurry during grouting.Furthermore,a physical model of grouting in a karst conduit is developed.Second,the effectiveness of the SFS method in grouting simulation is verified by the experiment of grouting conduit in flowing water.The SFS method enables real-time monitoring of fluid velocity and pressure during grouting in flowing water and provides a feasible calculation method for revealing the grouting plugging mechanism in complex karst conduits at different engineering scales.In addition,it can be used to guide the design of grouting tests in flowing water,improve cost efficiency,and provide theoretical basis for optimizing grouting design and slurry selection.
基金Projects(11702235,51641905,41472269) supported by the National Natural Science Foundation of ChinaProject(2017JJ3290) supported by the Natural Science Foundation of Hunan Province,China+1 种基金Project(17C1540) supported by the Scientific Research Foundation of Education Department of Hunan Province,ChinaProject(16GES07) supported by the Open Research Fund of Hunan Key Laboratory of Geomechanics and Engineering Safety,China
文摘The particle simulation method is used to study the effects of loading waveforms （i.e. square, sinusoidal and triangle waveforms） on rock damage at mesoscopic scale. Then some influencing factors on rock damage at the mesoscopic scale, such as loading frequency, stress amplitude, mean stress, confining pressure and loading sequence, are also investigated with sinusoidal waveform in detail. The related numerical results have demonstrated that： 1） the loading waveform has a certain effect on rock failure processes. The square waveform has the most damage within these waveforms, while the triangle waveform has less damage than sinusoidal waveform. In each cycle, the number of microscopic cracks increases in the loading stage, while it keeps nearly constant in the unloading stage. 2） The loading frequency, stress amplitude, mean stress, confining pressure and loading sequence have considerable effects on rock damage subjected to cyclic loading. The higher the loading frequency, stress amplitude and mean stress, the greater the damage the rock accumulated; in contrast, the lower the confining pressure, the greater the damage the rock has accumulated. 3） There is a threshold value of mean stress and stress amplitude, below which no further damage accumulated after the first few cycle loadings. 4） The high-to-low loading sequence has more damage than the low-to-high loading sequence, suggesting that the rock damage is loading-path dependent.
基金supported by the National Natural Science Foundation of China (No. 52075350)the Special City School Strategic Cooperation Project of Sichuan University and Zigong (No.2021CDZG-3)
文摘The control rod drive mechanism(CRDM)is an essential part of the control and safety protection system of pressurized water reactors.Current CRDM simulations are mostly performed collectively using a single method,ignoring the influence of multiple motion units and the differences in various features among them,which strongly affect the efficiency and accuracy of the simulations.In this study,we constructed a flow field fusion simulation method based on model features by combining key motion unit analysis and various simulation methods and then applied the method to the CRDM simulation process.CRDM performs motion unit decomposition through the structural hierarchy of function-movement-action method,and the key meta-actions are identified as the nodes in the flow field simulation.We established a fused feature-based multimethod simulation process and processed the simulation methods and data according to the features of the fluid domain space and the structural complexity to obtain the fusion simulation results.Compared to traditional simulation methods and real measurements,the simulation method provides advantages in terms of simulation efficiency and accuracy.
基金The National Natural Science Foundation of China(No.62173122)the Hebei Key Project of Natural Science Foundation(No.F2021201031)。
文摘A passive simulation method based on the six degrees of freedom(6-DOF)model and dynamic mesh is proposed according to the working principle to study the dynamic characteristics of the turbine flow sensors.This simulation method controls the six degrees of freedom of the impeller using the user-defined functions(UDF)program so that it can only rotate under the impact of fluid.The impeller speed can be calculated in real-time,and the inlet speed can be set with time to obtain the dynamic performance of the turbine flow sensors.Based on this simulation method,three turbine flow sensors with different diameters were simulated,and the reliability of the simulation method was verified by both steady-state and unsteady-state experiments.The results show that the trend of meter factor with flow rate acquired from the simulation is close to the experimental results.The deviation between the simulation and experiment results is low,with a maximum deviation of 2.88%.In the unsteady simulation study,the impeller speed changed with the inlet velocity of the turbine flow sensor,showing good tracking performance.The passive simulation method can be used to predict the dynamic performance of the turbine flow sensor.
基金Projects(51206011,U1937201)supported by the National Natural Science Foundation of ChinaProject(20200301040RQ)supported by the Science and Technology Development Program of Jilin Province,China+1 种基金Project(JJKH20190541KJ)supported by the Education Department of Jilin Province,ChinaProject(18DY017)supported by Changchun Science and Technology Program of Changchun City,China。
文摘The investigation was carried out on the technical problems of finishing the inner surface of elbow parts and the action mechanism of particles in elbow precision machining by abrasive flow.This work was analyzed and researched by combining theory,numerical and experimental methods.The direct simulation Monte Carlo(DSMC)method and the finite element analysis method were combined to reveal the random collision of particles during the precision machining of abrasive flow.Under different inlet velocity,volume fraction and abrasive particle size,the dynamic pressure and turbulence flow energy of abrasive flow in elbow were analyzed,and the machining mechanism of particles on the wall and the influence of different machining parameters on the precision machining quality of abrasive flow were obtained.The test results show the order of the influence of different parameters on the quality of abrasive flow precision machining and establish the optimal process parameters.The results of the surface morphology before and after the precision machining of the inner surface of the elbow are discussed,and the surface roughness Ra value is reduced from 1.125μm to 0.295μm after the precision machining of the abrasive flow.The application of DSMC method provides special insights for the development of abrasive flow technology.
基金Project supported by the National Natural Science Foundation of China (Grant No.50608069)
文摘To overcome the default of single search tendency, the ants in the colony are divided into several sub-groups. The ants in different subgroups have different trail information and expectation coefficients. The simulated annealing method is introduced to the algorithm. Through setting the temperature changing with the iterations, after each turn of tours, the solution set obtained by the ants is taken as the candidate set. The update set is obtained by adding the solutions in the candidate set to the previous update set with the probability determined by the temperature. The solutions in the candidate set are used to update the trail information. In each turn of updating, the current best solution is also used to enhance the trail information on the current best route. The trail information is reset when the algorithm is in stagnation state. The computer experiments demonstrate that the proposed algorithm has higher stability and convergence speed.
文摘One of the most critical and complicated steps in mine design is a selection of suitable mining method based upon geological,geotechnical,geographical,safety and economical parameters.The aim of this study is developing a Monte Carlo simulation to selection the optimum mining method by using effective and major criteria and at the same time,taking subjective judgments of decision makers into consideration.Proposed approach is based on the combination of Monte Carlo simulation with conventional Analytic Hierarchy Process(AHP).Monte Carlo simulation is used to determine the confdence level of each alternative’s score,is calculated by AHP,with the respect to the variance of decision makers’opinion.The proposed method is applied for Jajarm Bauxite Mine in Iran and eventually the most appropriate mining methods for this mine are ranked.
基金supported by the National Water Pollution Control and Management Technology Major Projects(Grant No. 2009ZX07423-001)the National Natural Science Foundation of China (Grants No.51179069and 40971300)the Fundamental Research Funds for the Central Universities (Grants No.10QX43,09MG16,and 10QG23)
文摘This paper presents a rapid and simple risk calculation method for large and complex engineering systems, the simulated maximum entropy method （SMEM）, which is based on integration of the advantages of the Monte Carlo and maximum entropy methods, thus avoiding the shortcoming of the slow convergence rate of the Monte Carlo method in risk calculation. Application of SMEM in the calculation of reservoir flood discharge risk shows that this method can make full use of the known information under the same conditions and obtain the corresponding probability distribution and the risk value. It not only greatly improves the speed, compared with the Monte Carlo method, but also provides a new approach for the risk calculation in large and complex engineering systems.
文摘In this paper .fatigue strength of T-type tubularjoints subjected to in-planebending or out-of-plane bending load is investigated By, considering material constantsand initial crack sizes as random variables and applying Monte Carlo simulation method, we have giyen a statistical analvsis on fatigue life. Simultaneously, linearregression analyses of computed resuldts are performed and compared with that of theknown experimental data.
基金financially supported by the Program for New Century Excellent Talents in University(No.NCET-13-0229,NCET-09-0396)the National Science & Technology Key Projects of Numerical Control(No.2012ZX04010-031,2012ZX0412-011)the National High Technology Research and Development Program("863"Program)of China(No.2013031003)
文摘For the sake of a more accurate shell boundary and calculation of radiation heat transfer in the Directional Solidification(DS) process, a radiation heat transfer model based on the Finite Element Method(FEM)is developed in this study. Key technologies, such as distinguishing boundaries automatically, local matrix and lumped heat capacity matrix, are also stated. In order to analyze the effect of withdrawing rate on DS process,the solidification processes of a complex superalloy turbine blade in the High Rate Solidification(HRS) process with different withdrawing rates are simulated; and by comparing the simulation results, it is found that the most suitable withdrawing rate is determined to be 5.0 mm·min^(-1). Finally, the accuracy and reliability of the radiation heat transfer model are verified, because of the accordance of simulation results with practical process.
基金Supported by National Key Basic Research Program of China(973 Program,Grant No.2014CB046405)State Key Laboratory of Fluid Power and Mechatronic Systems(Zhejiang University)Open Fund Project(Grant No.GZKF-201502)Hebei Military and Civilian Industry Development Funds Projects of China(Grant No.2015B060)
文摘Each joint of hydraulic drive quadruped robot is driven by the hydraulic drive unit（HDU）, and the contacting between the robot foot end and the ground is complex and variable, which increases the difficulty of force control inevitably. In the recent years, although many scholars researched some control methods such as disturbance rejection control, parameter self-adaptive control, impedance control and so on, to improve the force control performance of HDU, the robustness of the force control still needs improving. Therefore, how to simulate the complex and variable load characteristics of the environment structure and how to ensure HDU having excellent force control performance with the complex and variable load characteristics are key issues to be solved in this paper. The force control system mathematic model of HDU is established by the mechanism modeling method, and the theoretical models of a novel force control compensation method and a load characteristics simulation method under different environment structures are derived, considering the dynamic characteristics of the load stiffness and the load damping under different environment structures. Then, simulation effects of the variable load stiffness and load damping under the step and sinusoidal load force are analyzed experimentally on the HDU force control performance test platform, which provides the foundation for the force control compensation experiment research. In addition, the optimized PID control parameters are designed to make the HDU have better force control performance with suitable load stiffness and load damping, under which the force control compensation method is introduced, and the robustness of the force control system with several constant load characteristics and the variable load characteristics respectively are comparatively analyzed by experiment. The research results indicate that if the load characteristics are known, the force control compensation method presented in this paper has positive compensation effects on the load characteristics variation, i.e., this method decreases the effects of the load characteristics variation on the force control performance and enhances the force control system robustness with the constant PID parameters, thereby, the online PID parameters tuning control method which is complex needs not be adopted. All the above research provides theoretical and experimental foundation for the force control method of the quadruped robot joints with high robustness.
基金Project(41074085)supported by the National Natural Science Foundation of ChinaProject(NCET-12-0551)supported by the Funds for New Century Excellent Talents in University,ChinaProject supported by Shenghua Yuying Program of Central South University,China
文摘An element-free Galerkin method(EFGM) is used to solve the two-dimensional(2D) ground penetrating radar(GPR)modelling problems, due to its simple pre-processing, the absence of elements and high accuracy. Different from element-based numerical methods, this approach makes nodes free from the elemental restraint and avoids the explicit mesh discretization. First, we derived the boundary value problem for the 2D GPR simulation problems. Second, a penalty function approach and a boundary condition truncated method were used to enforce the essential and the absorbing boundary conditions, respectively. A three-layered GPR model was used to verify our element-free approach. The numerical solutions show that our solutions have an excellent agreement with solutions of a finite element method(FEM). Then, we used the EFGM to simulate one more complex model to show its capability and limitations. Simulation results show that one obvious advantage of EFGM is the absence of element mesh, which makes the method very flexible. Due to the use of MLS fitting, a key feature of EFM, is that both the dependent variable and its gradient are continuous and have high precision.
