A new hydraulic system of a novel automatic transmission (AT) was designed. The dimension and structure of valves and cylinders were designed by theoretical calculation. The dynamic simulation model of hydraulic syste...A new hydraulic system of a novel automatic transmission (AT) was designed. The dimension and structure of valves and cylinders were designed by theoretical calculation. The dynamic simulation model of hydraulic system of AT was established by ITI-SimulationX. Simulation results and theoretical design results were compared to confirm the simulation model. Based on the confirmed simulation model, the simulation results of pressure and flow of the hydraulic system were analyzed. The dynamic simulation method is very helpful for designing and analyzing the performance of hydraulic system and further optimization design. The theoretical design method and dynamic simulation model are feasible for the real industrial applications. The research results can be used in hydraulic system design and optimization.展开更多
With comprehensive considerations of the operational safety and collection efficiency for the tracked miner collecting the seafloor poly-metallic nodules, two new improved mining paths for the miner on the deep seaflo...With comprehensive considerations of the operational safety and collection efficiency for the tracked miner collecting the seafloor poly-metallic nodules, two new improved mining paths for the miner on the deep seafloor were proposed. Compared to the conventional mining path, the design principles and superiorities of the two new paths are that the miner turning with relative long radius should avoid large sinkage and high slip, so as to ensure its operational safety, while the space between its straight-line trajectories before and after the turning is optimum, which is designed as the total width of the miner, and collect nodules as more as possible, so as to ensure its collection efficiency. To realize the new mining paths, theoretical designs and quantitative calculations were carried out to determine the exact positions for the speed controls of the miner during its whole operation process. With the new dynamic model of the miner, and through regulations of the speeds of the left and right tracks of the miner on the exact motion positions according to the theoretical calculations, the two new improved mining paths for the miner on the seafloor were successfully simulated, thus the turning radius of the miner in the simulation is about 21.8 m, while the distance between the straight-line trajectories before and after the turning is about 5.2 m. The dynamic simulation results preliminarily prove the feasibility of these two new mining paths, and further can provide important theoretical guidance and useful technical reference for the practical tracked miner operation and control on the seafloor.展开更多
Accurate 3-D fracture network model for rock mass in dam foundation is of vital importance for stability,grouting and seepage analysis of dam foundation.With the aim of reducing deviation between fracture network mode...Accurate 3-D fracture network model for rock mass in dam foundation is of vital importance for stability,grouting and seepage analysis of dam foundation.With the aim of reducing deviation between fracture network model and measured data,a 3-D fracture network dynamic modeling method based on error analysis was proposed.Firstly,errors of four fracture volume density estimation methods(proposed by ODA,KULATILAKE,MAULDON,and SONG)and that of four fracture size estimation methods(proposed by EINSTEIN,SONG and TONON)were respectively compared,and the optimal methods were determined.Additionally,error index representing the deviation between fracture network model and measured data was established with integrated use of fractal dimension and relative absolute error(RAE).On this basis,the downhill simplex method was used to build the dynamic modeling method,which takes the minimum of error index as objective function and dynamically adjusts the fracture density and size parameters to correct the error index.Finally,the 3-D fracture network model could be obtained which meets the requirements.The proposed method was applied for 3-D fractures simulation in Miao Wei hydropower project in China for feasibility verification and the error index reduced from 2.618 to 0.337.展开更多
An improved whole model of beam pumping system was built. In the detail, for surface transmission system(STS), a new mathematical model was established considering the influence of some factors on the STS's torsio...An improved whole model of beam pumping system was built. In the detail, for surface transmission system(STS), a new mathematical model was established considering the influence of some factors on the STS's torsional vibration, such as the time variation characteristic of equivalent stiffness of belt and equivalent rotational inertia of crank. For the sucker rod string(SRS), an improved mathematical model was built considering the influence of some parameters on the SRS's longitudinal vibration, such as the nonlinear friction of plunger, hydraulic loss of pump and clearance leakage. The dynamic response and system efficiency of whole system were analyzed. The results show that there is a jumping phenomenon in the amplitude frequency curve, and the system efficiency is sensitive to motor power, pump diameter, stroke number, ratio of gas and oil, and submergence depth. The simulation results have important significance for improving the efficiency of beam pumping system.展开更多
Based on the theory of complex adaptive system(CAS),the optimal allocation model of water resources in sewage irrigation areas was established,which provided new ideas and application value for the rational utilizatio...Based on the theory of complex adaptive system(CAS),the optimal allocation model of water resources in sewage irrigation areas was established,which provided new ideas and application value for the rational utilization of agricultural production and waste water resources.The results demonstrated that the difference of crop energy capture mainly depended on the development stage.Waste water with a certain concentration was able to promote crop growth,while excessive concentration inhibited crop growth.The correlation between water absorption rate and leaf area index was close(R=0.9498,p<0.01).The amount of bad seeds increased at a speed of 34.7·d^-1,when system irrigated randomly in the seedling stage,while it tended to remain stable at a speed of 0.3·d^-1 after plants entering the mature stage which impacted the total yields of crops.展开更多
A deep understanding of explosive sensitivities and their factors is important for safe and reliable applications.However,quantitative prediction of the sensitivities is difficult.Here,reactive molecular dynamics simu...A deep understanding of explosive sensitivities and their factors is important for safe and reliable applications.However,quantitative prediction of the sensitivities is difficult.Here,reactive molecular dynamics simulation models for high-speed piston impacts on explosive supercells were established.Simulations were also performed to investigate shock-induced reactions of various high-energy explosives.The fraction of reacted explosive molecules in an initial supercell changed linearly with the propagation distance of the shock-wave front.The corresponding slope could be used as a reaction rate for a specific shock-loading velocity.Reaction rates that varied with the shock-loading pressure exhibited two-stage linearities with different slopes.The two inflection points corresponded to the initial and accelerated reactions,which respectively correlated to the thresholds of shock-induced ignition and detonation.Therefore,the ignition and detonation critical pressures could be determined.The sensitivity could then be a quantitative prediction of the critical pressure.The accuracies of the quantitative shock sensitivity predictions were verified by comparing the impact and shock sensitivities of common explosives and the characteristics of anisotropic shock-induced reactions.Molecular dynamics simulations quantitatively predict and rank shock sensitivities by using only crystal structures of the explosives.Overall,this method will enable the design and safe use of explosives.展开更多
According to the recently developed single-trough floating machine with the world's largest volume(inflatable mechanical agitation flotation machine with volume of 320 m3) in China, the gas-fluid two-phase flow in...According to the recently developed single-trough floating machine with the world's largest volume(inflatable mechanical agitation flotation machine with volume of 320 m3) in China, the gas-fluid two-phase flow in flotation cell was simulated using computational fluid dynamics method. It is shown that hexahedral mesh scheme is more suitable for the complex structure of the flotation cell than tetrahedral mesh scheme, and a mesh quality ranging from 0.7 to 1.0 is obtained. Comparative studies of the standard k-ε, k-ω and realizable k-ε turbulence models were carried out. It is indicated that the standard k-ε turbulence model could give a result relatively close to the practice and the liquid phase flow field is well characterized. In addition, two obvious recirculation zones are formed in the mixing zones, and the pressure on the rotor and stator is well characterized. Furthermore, the simulation results using improved standard k-ε turbulence model show that surface tension coefficient of 0.072, drag model of Grace and coefficient of 4, and lift coefficient of 0.001 can be achieved. The research results suggest that gas-fluid two-phase flow in large flotation cell can be well simulated using computational fluid dynamics method.展开更多
The effects of tensile temperatures ranging from 100 K to 900 K on the phase transition of hexagonal close-packed(HCP)zirconium were investigated by molecular dynamics simulations,which were combined with experimental...The effects of tensile temperatures ranging from 100 K to 900 K on the phase transition of hexagonal close-packed(HCP)zirconium were investigated by molecular dynamics simulations,which were combined with experimental observation under high resolution transmission electron microscopy.The results show that externally applied loading first induced the HCP to body-centered cubic(BCC)phase transition in the Pitsch-Schrader(PS)orientation relationship(OR).Then,the face-centered cubic(FCC)structure transformed from the BCC phase in the Bain path.However,the HCP-to-BCC transition was incomplete at 100 K and 300 K,resulting in a prismatic-type OR between the FCC and original HCP phase.Additionally,at the temperature ranging from 100 K to 600 K,the inverse BCC-to-HCP transition occurred locally following other variants of the PS OR,resulting in a basal-type relation between the newly generated HCP and FCC phases.A higher tensile temperature promoted the amount of FCC phase transforming into the BCC phase when the strain exceeded 45%.Besides,the crystal stretched at lower temperatures exhibits relatively higher strength but by the compromise of plasticity.This study reveals the deformation mechanisms in HCP-Zr at different temperatures,which may provide a better understanding of the deformation mechanism of zirconium alloys under different application environments.展开更多
To investigate the effect of void defects on the shock response of hexanitrohexaazaisowurtzitane(CL-20)co-crystals,shock responses of CL-20 co-crystals with energetic materials ligands trinitrotoluene(TNT),1,3-dinitro...To investigate the effect of void defects on the shock response of hexanitrohexaazaisowurtzitane(CL-20)co-crystals,shock responses of CL-20 co-crystals with energetic materials ligands trinitrotoluene(TNT),1,3-dinitrobenzene(DNB),solvents ligands dimethyl carbonate(DMC) and gamma-butyrolactone(GBL)with void were simulated,using molecular dynamics method and reactive force field.It is found that the CL-20 co-crystals with void defects will form hot spots when impacted,significantly affecting the decomposition of molecules around the void.The degree of molecular fragmentation is relatively low under the reflection velocity of 2 km/s,and the main reactions are the formation of dimer and the shedding of nitro groups.The existence of voids reduces the safety of CL-20 co-crystals,which induced the sensitivity of energetic co-crystals CL-20/TNT and CL-20/DNB to increase more significantly.Detonation has occurred under the reflection velocity of 4 km/s,energetic co-crystals are easier to polymerize than solvent co-crystals,and are not obviously affected by voids.The results show that the energy of the wave decreases after sweeping over the void,which reduces the chemical reaction frequency downstream of the void and affects the detonation performance,especially the solvent co-crystals.展开更多
In this study,circular dichroism(CD)and molecular dynamics(MD)simulation were used to investigate the thermal unfolding pathway of staphylococcal enterotoxin B(SEB)at temperatures of 298–371 and 298–500 K,and the re...In this study,circular dichroism(CD)and molecular dynamics(MD)simulation were used to investigate the thermal unfolding pathway of staphylococcal enterotoxin B(SEB)at temperatures of 298–371 and 298–500 K,and the relationship between the experimental and simulation results were explored.Our computational findings on the secondary structure of SEB showed that at room temperature,the CD spectroscopic results were highly consistent with the MD results.Moreover,under heating conditions,the changing trends of helix,sheet and random coil obtained by CD spectral fitting were highly consistent with those obtained by MD.In order to gain a deeper understanding of the thermal stability mechanism of SEB,the MD trajectories were analyzed in terms of root mean square deviation(RMSD),secondary structure assignment(SSA),radius of gyration(R_(g)),free energy surfaces(FES),solvent-accessible surface area(SASA),hydrogen bonds and salt bridges.The results showed that at low heating temperature,domain Ⅰ without loops(omitting the mobile loop region)mainly relied on hydrophobic interaction to maintain its thermal stability,whereas the thermal stability of domain Ⅱ was mainly controlled by salt bridges and hydrogen bonds.Under high heating temperature conditions,the hydrophobic interactions in domain Ⅰ without loops were destroyed and the secondary structure was almost completely lost,while domain Ⅱ could still rely on salt bridges as molecular staples to barely maintain the stability of the secondary structure.These results help us to understand the thermodynamic and kinetic mechanisms that maintain the thermal stability of SEB at the molecular level,and provide a direction for establishing safer and more effective food sterilization processes.展开更多
Molecular dynamics simulations are performed to observe the evolutions of 512 and 51262 cage-like water clusters filled with or without a methane molecule immersed in bulk liquid water at 250 K and 230 K. The lifetime...Molecular dynamics simulations are performed to observe the evolutions of 512 and 51262 cage-like water clusters filled with or without a methane molecule immersed in bulk liquid water at 250 K and 230 K. The lifetimes of these clusters are calculated according to their Lindemann index δ (t) using the criteria of δ≥0.07. For both the filled and empty clusters, we find the dynamics of bulk water determines the lifetimes of cage-like water clusters, and that the lifetime of 512 62 cage-like cluster is the same as that of 512 cage-like cluster. Although the methane molecule indeed makes the filled cage-like cluster more stable than the empty one, the empty cage-like cluster still has chance to be long-lived compared with the filled clusters. These observations support the labile cluster hypothesis on the formation mechanisms of gas hydrates.展开更多
The microstructures and thermodynamic properties of mixed systems comprising pyridinium ionic liquid[HPy][BF_(4)]and acetonitrile at different mole fractions were studied using molecular dynamics simulation in this wo...The microstructures and thermodynamic properties of mixed systems comprising pyridinium ionic liquid[HPy][BF_(4)]and acetonitrile at different mole fractions were studied using molecular dynamics simulation in this work.The following properties were determined:density,self-diffusion coefficient,excess molar volume,and radial distribution function.The results show that with an increase in the mole fraction of[HPy][BF_(4)],the self-diffusion coefficient decreases.Additionally,the excess molar volume initially decreases,reaches a minimum,and then increases.The rules of radial distribution functions(RDFs)of characteristic atoms are different.With increasing the mole fraction of[HPy][BF_(4)],the first peak of the RDFs of HA1-F decreases,while that of CT6-CT6 rises at first and then decreases.This indicates that the solvent molecules affect the polar and non-polar regions of[HPy][BF_(4)]differently.展开更多
The dynamic model of cold rolling mill based on strip flatness and thickness integrated control was proposed,containing the following sub-models:the rolling process model,the dynamic model of rolls along axial directi...The dynamic model of cold rolling mill based on strip flatness and thickness integrated control was proposed,containing the following sub-models:the rolling process model,the dynamic model of rolls along axial direction,and the compensation model.Based on the rule of volume flow rate,the dynamic rolling process model was built.The work roll and backup roll were taken as elastic continuous bodies,the effect of shear and moment of inertia were taken into consideration,and then the dynamic model of rolls was built.The two models were coupled together,and the dynamic model of rolling mill was built.In the dynamic model,the thermal expansion of the rolls,the wear of the rolls and other related parameters can not be considered.In order to compensate the dynamic model,the coupled static model of rolls and strip was applied.Then,according to the inner relationship of these models,the dynamic model and the compensation model were coupled,and the dynamic model of rolling mill based on the strip flatness and thickness integrated control was built.The dynamic simulation of the rolling process was made,and the dynamic thickness and the dynamic flatness information were obtained.This model not only provides a theory basis for the virtual rolling,but also provides a platform for the application of advanced control theory.展开更多
Switching expansion reduction(SER)uses a switch valve instead of the throttle valve to realize electronically controlled pressure reduction for high pressure pneumatics.A comprehensive and interactive pneumatic simula...Switching expansion reduction(SER)uses a switch valve instead of the throttle valve to realize electronically controlled pressure reduction for high pressure pneumatics.A comprehensive and interactive pneumatic simulation model according to the experimental setup of SER has been built.The mathematical model considers heat exchanges,source air pressure and temperature,environmental temperatures and heat transfer coefficients variations.In addition,the compensation for real gas effect is used in the model building.The comparison between experiments and simulations of SER indicates that,to compensate the real gas effect in high pressure discharging process,the thermal capacity of air supply container in simulation should be less than the actual value.The higher the pressure range,the greater the deviation.Simulated and experimental results are highly consistent within pressure reduction ratios ranging from 1.4 to 20 and output air mass flow rates ranging from 3.5 to 132 g/s,which verifies the high adaptability of SER and the validity of the mathematic model and the compensation method.展开更多
The characteristic of easy sintering of aluminum nanoparticle(ANP)limits its application in solid propellants.Coating ANP with fluoropolymer could effectively improve its combustion performance.To find out how the coa...The characteristic of easy sintering of aluminum nanoparticle(ANP)limits its application in solid propellants.Coating ANP with fluoropolymer could effectively improve its combustion performance.To find out how the coating layer inhibits sintering and promotes complete combustion of particles from an atomic view,a comparative study has been done for bare ANP and PTFE coated ANP by using reactive molecular dynamics simulations.The sintering process is quantified by shrinkage ratio and gyration radius.Our results show that,at the same heating rate and combustion temperatures,bare ANPs are sintered together after the temperature exceeds the melting point of aluminum but the decomposition of PTFE coating layer pushes particles away and increases reaction surface area by producing small Al-F clusters.The sintering of ANPs which are heated in PTFE is alleviated compared with particles heated in oxygen,but particles still sinter together due to the lack of intimate contact between PTFE and alumina surface.The effect of temperature on the combustion of PTFE coated ANPs is also studied from 1000 to3500 K.The number density analysis shows the particles will not be sintered at any temperature.Aluminum fluoride prefers diffusing to the external space and the remained particles are mainly composed of Al,C and O.Fast ignition simulations are performed by adopting micro canonical ensemble.With the expansion of aluminum core and the melting of alumina shell,bare ANPs are sintered into a liquid particle directly.For PTFE coated ANPs,the volatilization of gaseous aluminum fluoride products continually endows particles opposite momentum.展开更多
The interface properties between hydrated cement paste(hcp)and aggregates largely determine the various performances of concrete.In this work,molecular dynamics simulations were employed to explore the atomistic inter...The interface properties between hydrated cement paste(hcp)and aggregates largely determine the various performances of concrete.In this work,molecular dynamics simulations were employed to explore the atomistic interaction mechanisms between the commonly used aggregate phase calcite/silica and calcium silicate hydrates(C-S-H),as well as the effect of moisture.The results suggest that the C-S-H/calcite interface is relatively strong and stable under both dry and moist conditions,which is caused by the high-strength interfacial connections formed between calcium ions from calcite and high-polarity non-bridging oxygen atoms from the C-S-H surface.Silica can be also adsorbed on the dry C-S-H surface by the H-bonds;however,the presence of water molecules on the interface may substantially decrease the affinities.Furthermore,the dynamics interface separation tests of C-S-H/aggregates were also implemented by molecular dynamics.The shape of the calculated stress-separation distance curves obeys the quasi-static cohesive law obtained experimentally.The moisture conditions and strain rates were found to affect the separation process of C-S-H/silica.A wetter interface and smaller loading rate may lead to a lower adhesion strength.The mechanisms interpreted here may shed new lights on the understandings of hcp/aggregate interactions at a nano-length scale and creation of high performance cementitious materials.展开更多
Nanograins are characterized by a typical grain size from 1 to 100 nm. Molecular dynamics simulations have been carried out for the nanograin sphere with the diameters from 1.45 to 10.12 nm. We study the influence of ...Nanograins are characterized by a typical grain size from 1 to 100 nm. Molecular dynamics simulations have been carried out for the nanograin sphere with the diameters from 1.45 to 10.12 nm. We study the influence of grain size on structure and diffusion properties of the nanograins. The results reveal that as the grain size is reduced, the fraction of grain surface increases significantly, and the surface width is approximately constant; the mean atomic energy of the surface increases distinctly, but that of the grain interior varies insignificantly; the diffusion coefficient is increased sharply, and the relation of the diffusion coefficient and the grain size is close to exponential relation below 10 nm.展开更多
The distribution of Al （j） and the structural units distribution of Qi T in calcium aluminosilicate melts were studied by means of molecular dynamics simulation. The results show that provided there exists lower-fie...The distribution of Al （j） and the structural units distribution of Qi T in calcium aluminosilicate melts were studied by means of molecular dynamics simulation. The results show that provided there exists lower-field strength cation relative to Al3+, such as alkaline and alkaline earth metals, Al will be four-coordinated but not six-coordinated. Meanwhile, if there exist a large number of higher-field strength cations such as Si4+ and little lower-field strength cation, six-coordinated aluminum will be formed. The relation of structural units distribution of Qi T with chemical composition shift was also extracted, showing that as Ca2+ exists, the distributions of Qi Si, Qi Al or Qi T have the similar changing trend with the variation of component. Because of high-temperature effect, the Al-tetrahedral units in melts are greatly active and unstable and there exist dynamic transforming equilibria of Al（3）Al（4） and （Al（5））Al（4）. The three-coordinated oxygen and charge-compensated bridging oxygen are proposed to explain phenomena of the negative charge redundancy of AlO4 and location of network modifier with charge-compensated function in aluminosilicate melts.展开更多
The pre-crack blast technology has been used to control the induction caving area in the roof. The key is to form the pre-crack seam and predict the effect of the seam. The H-J-C blast model was built in the roof. Bas...The pre-crack blast technology has been used to control the induction caving area in the roof. The key is to form the pre-crack seam and predict the effect of the seam. The H-J-C blast model was built in the roof. Based on the theories of dynamic strength and failure criterion of dynamic rock, the rock dynamic damage and the evolution of pre-crack seam were simulated by the tensile damage and shear failure of the model. According to the actual situation of No. 92 ore body test stope at Tongkeng Mine, the formation process of the pre-crack blast seam was simulated by Ansys/Ls-dyna software, the pre-crack seam was inspected by a system of digital panoramic borehole camera. The pre-crack seam was inspected by the system of digital panoramic borehole in the roof. The results of the numerical simulation and inspection show that in the line of centers of pre-hole, the minimum of the tensile stress reaches 20 MPa, which is much larger than 13.7 MPa of the dynamic tensile strength of rock. The minimum particle vibration velocity reaches 50 cm/s, which is greater than 30-40 cm/s of the allowable vibration velocity. It is demonstrated that the rock is destroyed near the center line and the pre-crack is successfully formed by the large diameters and large distances pre-crack holes in the roof.展开更多
An efficient and real-time simulation method is proposed for the dynamic electromagnetic characteristics of cluster targets to meet the requirements of engineering practical applications.First,the coordinate transform...An efficient and real-time simulation method is proposed for the dynamic electromagnetic characteristics of cluster targets to meet the requirements of engineering practical applications.First,the coordinate transformation method is used to establish a geometric model of the observation scene,which is described by the azimuth angles and elevation angles of the radar in the target reference frame and the attitude angles of the target in the radar reference frame.Then,an approach for dynamic electromagnetic scattering simulation is proposed.Finally,a fast-computing method based on sparsity in the time domain,space domain,and frequency domain is proposed.The method analyzes the sparsity-based dynamic scattering characteristic of the typical cluster targets.The error between the sparsity-based method and the benchmark is small,proving the effectiveness of the proposed method.展开更多
基金Project(911901204) supported by Youth Innovation Foundation of Beijing University of Aeronautics and Astronautics
文摘A new hydraulic system of a novel automatic transmission (AT) was designed. The dimension and structure of valves and cylinders were designed by theoretical calculation. The dynamic simulation model of hydraulic system of AT was established by ITI-SimulationX. Simulation results and theoretical design results were compared to confirm the simulation model. Based on the confirmed simulation model, the simulation results of pressure and flow of the hydraulic system were analyzed. The dynamic simulation method is very helpful for designing and analyzing the performance of hydraulic system and further optimization design. The theoretical design method and dynamic simulation model are feasible for the real industrial applications. The research results can be used in hydraulic system design and optimization.
基金Project(DYXM-115-04-02-01) supported by the National Deep-sea Technology Project of Development and Research, ChinaProject(2011QNZT058) supported by the Fundamental Research Funds for the Central Universities, ChinaProject(51105386) supported by the National Natural Science Foundation of China
文摘With comprehensive considerations of the operational safety and collection efficiency for the tracked miner collecting the seafloor poly-metallic nodules, two new improved mining paths for the miner on the deep seafloor were proposed. Compared to the conventional mining path, the design principles and superiorities of the two new paths are that the miner turning with relative long radius should avoid large sinkage and high slip, so as to ensure its operational safety, while the space between its straight-line trajectories before and after the turning is optimum, which is designed as the total width of the miner, and collect nodules as more as possible, so as to ensure its collection efficiency. To realize the new mining paths, theoretical designs and quantitative calculations were carried out to determine the exact positions for the speed controls of the miner during its whole operation process. With the new dynamic model of the miner, and through regulations of the speeds of the left and right tracks of the miner on the exact motion positions according to the theoretical calculations, the two new improved mining paths for the miner on the seafloor were successfully simulated, thus the turning radius of the miner in the simulation is about 21.8 m, while the distance between the straight-line trajectories before and after the turning is about 5.2 m. The dynamic simulation results preliminarily prove the feasibility of these two new mining paths, and further can provide important theoretical guidance and useful technical reference for the practical tracked miner operation and control on the seafloor.
基金Project(51321065)supported by the Innovative Research Groups of the National Natural Science Foundation of ChinaProject(2013CB035904)supported by the National Basic Research Program of China(973 Program)Project(51439005)supported by the National Natural Science Foundation of China
文摘Accurate 3-D fracture network model for rock mass in dam foundation is of vital importance for stability,grouting and seepage analysis of dam foundation.With the aim of reducing deviation between fracture network model and measured data,a 3-D fracture network dynamic modeling method based on error analysis was proposed.Firstly,errors of four fracture volume density estimation methods(proposed by ODA,KULATILAKE,MAULDON,and SONG)and that of four fracture size estimation methods(proposed by EINSTEIN,SONG and TONON)were respectively compared,and the optimal methods were determined.Additionally,error index representing the deviation between fracture network model and measured data was established with integrated use of fractal dimension and relative absolute error(RAE).On this basis,the downhill simplex method was used to build the dynamic modeling method,which takes the minimum of error index as objective function and dynamically adjusts the fracture density and size parameters to correct the error index.Finally,the 3-D fracture network model could be obtained which meets the requirements.The proposed method was applied for 3-D fractures simulation in Miao Wei hydropower project in China for feasibility verification and the error index reduced from 2.618 to 0.337.
基金Projects(50974108,51174175) supported by the National Natural Science Foundation of China
文摘An improved whole model of beam pumping system was built. In the detail, for surface transmission system(STS), a new mathematical model was established considering the influence of some factors on the STS's torsional vibration, such as the time variation characteristic of equivalent stiffness of belt and equivalent rotational inertia of crank. For the sucker rod string(SRS), an improved mathematical model was built considering the influence of some parameters on the SRS's longitudinal vibration, such as the nonlinear friction of plunger, hydraulic loss of pump and clearance leakage. The dynamic response and system efficiency of whole system were analyzed. The results show that there is a jumping phenomenon in the amplitude frequency curve, and the system efficiency is sensitive to motor power, pump diameter, stroke number, ratio of gas and oil, and submergence depth. The simulation results have important significance for improving the efficiency of beam pumping system.
基金Supported by the Science and Technology Research Project of the Ministry of Education(14YJCZH017)the Major State Basic Research Development Program of China(973 Program)(2017YFC0404503)+1 种基金Key Cultivation Project of Lingnan Normal University in 2019(LZ1903)Lingnan Normal University Special Talent Program(ZL2007)
文摘Based on the theory of complex adaptive system(CAS),the optimal allocation model of water resources in sewage irrigation areas was established,which provided new ideas and application value for the rational utilization of agricultural production and waste water resources.The results demonstrated that the difference of crop energy capture mainly depended on the development stage.Waste water with a certain concentration was able to promote crop growth,while excessive concentration inhibited crop growth.The correlation between water absorption rate and leaf area index was close(R=0.9498,p<0.01).The amount of bad seeds increased at a speed of 34.7·d^-1,when system irrigated randomly in the seedling stage,while it tended to remain stable at a speed of 0.3·d^-1 after plants entering the mature stage which impacted the total yields of crops.
基金supported by the National Natural Science Foundation of China(Grant No.11832006).
文摘A deep understanding of explosive sensitivities and their factors is important for safe and reliable applications.However,quantitative prediction of the sensitivities is difficult.Here,reactive molecular dynamics simulation models for high-speed piston impacts on explosive supercells were established.Simulations were also performed to investigate shock-induced reactions of various high-energy explosives.The fraction of reacted explosive molecules in an initial supercell changed linearly with the propagation distance of the shock-wave front.The corresponding slope could be used as a reaction rate for a specific shock-loading velocity.Reaction rates that varied with the shock-loading pressure exhibited two-stage linearities with different slopes.The two inflection points corresponded to the initial and accelerated reactions,which respectively correlated to the thresholds of shock-induced ignition and detonation.Therefore,the ignition and detonation critical pressures could be determined.The sensitivity could then be a quantitative prediction of the critical pressure.The accuracies of the quantitative shock sensitivity predictions were verified by comparing the impact and shock sensitivities of common explosives and the characteristics of anisotropic shock-induced reactions.Molecular dynamics simulations quantitatively predict and rank shock sensitivities by using only crystal structures of the explosives.Overall,this method will enable the design and safe use of explosives.
基金Project(51074027)supported by the National Natural Science Foundation of China
文摘According to the recently developed single-trough floating machine with the world's largest volume(inflatable mechanical agitation flotation machine with volume of 320 m3) in China, the gas-fluid two-phase flow in flotation cell was simulated using computational fluid dynamics method. It is shown that hexahedral mesh scheme is more suitable for the complex structure of the flotation cell than tetrahedral mesh scheme, and a mesh quality ranging from 0.7 to 1.0 is obtained. Comparative studies of the standard k-ε, k-ω and realizable k-ε turbulence models were carried out. It is indicated that the standard k-ε turbulence model could give a result relatively close to the practice and the liquid phase flow field is well characterized. In addition, two obvious recirculation zones are formed in the mixing zones, and the pressure on the rotor and stator is well characterized. Furthermore, the simulation results using improved standard k-ε turbulence model show that surface tension coefficient of 0.072, drag model of Grace and coefficient of 4, and lift coefficient of 0.001 can be achieved. The research results suggest that gas-fluid two-phase flow in large flotation cell can be well simulated using computational fluid dynamics method.
基金Projects(51901248,51828102)supported by the National Natural Science Foundation of ChinaProject(2018JJ3649)supported by the Natural Science Foundation of Hunan Province,ChinaProject(2019CX026)supported by the Innovation-driven Plan in Central South University,China。
文摘The effects of tensile temperatures ranging from 100 K to 900 K on the phase transition of hexagonal close-packed(HCP)zirconium were investigated by molecular dynamics simulations,which were combined with experimental observation under high resolution transmission electron microscopy.The results show that externally applied loading first induced the HCP to body-centered cubic(BCC)phase transition in the Pitsch-Schrader(PS)orientation relationship(OR).Then,the face-centered cubic(FCC)structure transformed from the BCC phase in the Bain path.However,the HCP-to-BCC transition was incomplete at 100 K and 300 K,resulting in a prismatic-type OR between the FCC and original HCP phase.Additionally,at the temperature ranging from 100 K to 600 K,the inverse BCC-to-HCP transition occurred locally following other variants of the PS OR,resulting in a basal-type relation between the newly generated HCP and FCC phases.A higher tensile temperature promoted the amount of FCC phase transforming into the BCC phase when the strain exceeded 45%.Besides,the crystal stretched at lower temperatures exhibits relatively higher strength but by the compromise of plasticity.This study reveals the deformation mechanisms in HCP-Zr at different temperatures,which may provide a better understanding of the deformation mechanism of zirconium alloys under different application environments.
基金supported by the National Natural Science Foundation of China (22275018)the Project of State Key Laboratory of Explosion Science and Technology (Beijing Institute of Technology)(Grant No.QNKT20-04)。
文摘To investigate the effect of void defects on the shock response of hexanitrohexaazaisowurtzitane(CL-20)co-crystals,shock responses of CL-20 co-crystals with energetic materials ligands trinitrotoluene(TNT),1,3-dinitrobenzene(DNB),solvents ligands dimethyl carbonate(DMC) and gamma-butyrolactone(GBL)with void were simulated,using molecular dynamics method and reactive force field.It is found that the CL-20 co-crystals with void defects will form hot spots when impacted,significantly affecting the decomposition of molecules around the void.The degree of molecular fragmentation is relatively low under the reflection velocity of 2 km/s,and the main reactions are the formation of dimer and the shedding of nitro groups.The existence of voids reduces the safety of CL-20 co-crystals,which induced the sensitivity of energetic co-crystals CL-20/TNT and CL-20/DNB to increase more significantly.Detonation has occurred under the reflection velocity of 4 km/s,energetic co-crystals are easier to polymerize than solvent co-crystals,and are not obviously affected by voids.The results show that the energy of the wave decreases after sweeping over the void,which reduces the chemical reaction frequency downstream of the void and affects the detonation performance,especially the solvent co-crystals.
文摘In this study,circular dichroism(CD)and molecular dynamics(MD)simulation were used to investigate the thermal unfolding pathway of staphylococcal enterotoxin B(SEB)at temperatures of 298–371 and 298–500 K,and the relationship between the experimental and simulation results were explored.Our computational findings on the secondary structure of SEB showed that at room temperature,the CD spectroscopic results were highly consistent with the MD results.Moreover,under heating conditions,the changing trends of helix,sheet and random coil obtained by CD spectral fitting were highly consistent with those obtained by MD.In order to gain a deeper understanding of the thermal stability mechanism of SEB,the MD trajectories were analyzed in terms of root mean square deviation(RMSD),secondary structure assignment(SSA),radius of gyration(R_(g)),free energy surfaces(FES),solvent-accessible surface area(SASA),hydrogen bonds and salt bridges.The results showed that at low heating temperature,domain Ⅰ without loops(omitting the mobile loop region)mainly relied on hydrophobic interaction to maintain its thermal stability,whereas the thermal stability of domain Ⅱ was mainly controlled by salt bridges and hydrogen bonds.Under high heating temperature conditions,the hydrophobic interactions in domain Ⅰ without loops were destroyed and the secondary structure was almost completely lost,while domain Ⅱ could still rely on salt bridges as molecular staples to barely maintain the stability of the secondary structure.These results help us to understand the thermodynamic and kinetic mechanisms that maintain the thermal stability of SEB at the molecular level,and provide a direction for establishing safer and more effective food sterilization processes.
基金supported by the National Natural Science Foundation of China(Grant No.40102005 and No.49725205).
文摘Molecular dynamics simulations are performed to observe the evolutions of 512 and 51262 cage-like water clusters filled with or without a methane molecule immersed in bulk liquid water at 250 K and 230 K. The lifetimes of these clusters are calculated according to their Lindemann index δ (t) using the criteria of δ≥0.07. For both the filled and empty clusters, we find the dynamics of bulk water determines the lifetimes of cage-like water clusters, and that the lifetime of 512 62 cage-like cluster is the same as that of 512 cage-like cluster. Although the methane molecule indeed makes the filled cage-like cluster more stable than the empty one, the empty cage-like cluster still has chance to be long-lived compared with the filled clusters. These observations support the labile cluster hypothesis on the formation mechanisms of gas hydrates.
文摘The microstructures and thermodynamic properties of mixed systems comprising pyridinium ionic liquid[HPy][BF_(4)]and acetonitrile at different mole fractions were studied using molecular dynamics simulation in this work.The following properties were determined:density,self-diffusion coefficient,excess molar volume,and radial distribution function.The results show that with an increase in the mole fraction of[HPy][BF_(4)],the self-diffusion coefficient decreases.Additionally,the excess molar volume initially decreases,reaches a minimum,and then increases.The rules of radial distribution functions(RDFs)of characteristic atoms are different.With increasing the mole fraction of[HPy][BF_(4)],the first peak of the RDFs of HA1-F decreases,while that of CT6-CT6 rises at first and then decreases.This indicates that the solvent molecules affect the polar and non-polar regions of[HPy][BF_(4)]differently.
基金Project(E2012203177)supported by the Natural Science Foundation of Hebei Province,ChinaProject(2011BAF15B01)supported by the National Science and Technology Support Plan of China+1 种基金Project(E2006001038)supported by Great Natural Science Foundation of Hebei Province,ChinaProject(NECSR-201202)supported by Open Project Program of National Engineering Research Center for Equipment and Technology of Cold Strip Rolling,China
文摘The dynamic model of cold rolling mill based on strip flatness and thickness integrated control was proposed,containing the following sub-models:the rolling process model,the dynamic model of rolls along axial direction,and the compensation model.Based on the rule of volume flow rate,the dynamic rolling process model was built.The work roll and backup roll were taken as elastic continuous bodies,the effect of shear and moment of inertia were taken into consideration,and then the dynamic model of rolls was built.The two models were coupled together,and the dynamic model of rolling mill was built.In the dynamic model,the thermal expansion of the rolls,the wear of the rolls and other related parameters can not be considered.In order to compensate the dynamic model,the coupled static model of rolls and strip was applied.Then,according to the inner relationship of these models,the dynamic model and the compensation model were coupled,and the dynamic model of rolling mill based on the strip flatness and thickness integrated control was built.The dynamic simulation of the rolling process was made,and the dynamic thickness and the dynamic flatness information were obtained.This model not only provides a theory basis for the virtual rolling,but also provides a platform for the application of advanced control theory.
基金Project(51205421)supported by the National Natural Science Foundation of ChinaProject(2012M521647)supported by the Postdoctoral Science Foundation of China
文摘Switching expansion reduction(SER)uses a switch valve instead of the throttle valve to realize electronically controlled pressure reduction for high pressure pneumatics.A comprehensive and interactive pneumatic simulation model according to the experimental setup of SER has been built.The mathematical model considers heat exchanges,source air pressure and temperature,environmental temperatures and heat transfer coefficients variations.In addition,the compensation for real gas effect is used in the model building.The comparison between experiments and simulations of SER indicates that,to compensate the real gas effect in high pressure discharging process,the thermal capacity of air supply container in simulation should be less than the actual value.The higher the pressure range,the greater the deviation.Simulated and experimental results are highly consistent within pressure reduction ratios ranging from 1.4 to 20 and output air mass flow rates ranging from 3.5 to 132 g/s,which verifies the high adaptability of SER and the validity of the mathematic model and the compensation method.
基金supported by the fellowship of China Postdoctoral Science Foundation(Grant No.2021TQ0267)。
文摘The characteristic of easy sintering of aluminum nanoparticle(ANP)limits its application in solid propellants.Coating ANP with fluoropolymer could effectively improve its combustion performance.To find out how the coating layer inhibits sintering and promotes complete combustion of particles from an atomic view,a comparative study has been done for bare ANP and PTFE coated ANP by using reactive molecular dynamics simulations.The sintering process is quantified by shrinkage ratio and gyration radius.Our results show that,at the same heating rate and combustion temperatures,bare ANPs are sintered together after the temperature exceeds the melting point of aluminum but the decomposition of PTFE coating layer pushes particles away and increases reaction surface area by producing small Al-F clusters.The sintering of ANPs which are heated in PTFE is alleviated compared with particles heated in oxygen,but particles still sinter together due to the lack of intimate contact between PTFE and alumina surface.The effect of temperature on the combustion of PTFE coated ANPs is also studied from 1000 to3500 K.The number density analysis shows the particles will not be sintered at any temperature.Aluminum fluoride prefers diffusing to the external space and the remained particles are mainly composed of Al,C and O.Fast ignition simulations are performed by adopting micro canonical ensemble.With the expansion of aluminum core and the melting of alumina shell,bare ANPs are sintered into a liquid particle directly.For PTFE coated ANPs,the volatilization of gaseous aluminum fluoride products continually endows particles opposite momentum.
基金Projects(6512009004A,51908119,U1706222)supported by the National Natural Science Foundation of ChinaProject(BK20190367)supported by the Natural Science Foundation of Jiangsu Province,China。
文摘The interface properties between hydrated cement paste(hcp)and aggregates largely determine the various performances of concrete.In this work,molecular dynamics simulations were employed to explore the atomistic interaction mechanisms between the commonly used aggregate phase calcite/silica and calcium silicate hydrates(C-S-H),as well as the effect of moisture.The results suggest that the C-S-H/calcite interface is relatively strong and stable under both dry and moist conditions,which is caused by the high-strength interfacial connections formed between calcium ions from calcite and high-polarity non-bridging oxygen atoms from the C-S-H surface.Silica can be also adsorbed on the dry C-S-H surface by the H-bonds;however,the presence of water molecules on the interface may substantially decrease the affinities.Furthermore,the dynamics interface separation tests of C-S-H/aggregates were also implemented by molecular dynamics.The shape of the calculated stress-separation distance curves obeys the quasi-static cohesive law obtained experimentally.The moisture conditions and strain rates were found to affect the separation process of C-S-H/silica.A wetter interface and smaller loading rate may lead to a lower adhesion strength.The mechanisms interpreted here may shed new lights on the understandings of hcp/aggregate interactions at a nano-length scale and creation of high performance cementitious materials.
基金supported by the National Naturl Science Foundation of China(No.10172088)the Potdoctoral Saience Foundation of China
文摘Nanograins are characterized by a typical grain size from 1 to 100 nm. Molecular dynamics simulations have been carried out for the nanograin sphere with the diameters from 1.45 to 10.12 nm. We study the influence of grain size on structure and diffusion properties of the nanograins. The results reveal that as the grain size is reduced, the fraction of grain surface increases significantly, and the surface width is approximately constant; the mean atomic energy of the surface increases distinctly, but that of the grain interior varies insignificantly; the diffusion coefficient is increased sharply, and the relation of the diffusion coefficient and the grain size is close to exponential relation below 10 nm.
文摘The distribution of Al （j） and the structural units distribution of Qi T in calcium aluminosilicate melts were studied by means of molecular dynamics simulation. The results show that provided there exists lower-field strength cation relative to Al3+, such as alkaline and alkaline earth metals, Al will be four-coordinated but not six-coordinated. Meanwhile, if there exist a large number of higher-field strength cations such as Si4+ and little lower-field strength cation, six-coordinated aluminum will be formed. The relation of structural units distribution of Qi T with chemical composition shift was also extracted, showing that as Ca2+ exists, the distributions of Qi Si, Qi Al or Qi T have the similar changing trend with the variation of component. Because of high-temperature effect, the Al-tetrahedral units in melts are greatly active and unstable and there exist dynamic transforming equilibria of Al（3）Al（4） and （Al（5））Al（4）. The three-coordinated oxygen and charge-compensated bridging oxygen are proposed to explain phenomena of the negative charge redundancy of AlO4 and location of network modifier with charge-compensated function in aluminosilicate melts.
基金Project(2006BAB02A02)supported by the National Key Technology R&D Program for the 11th Five-year Plan of ChinaProject(09JJ4025)supported by the National Natural Science Foundation of Hunan Province,ChinaProject(51074178)supported by the National Natural Science Foundation of China
文摘The pre-crack blast technology has been used to control the induction caving area in the roof. The key is to form the pre-crack seam and predict the effect of the seam. The H-J-C blast model was built in the roof. Based on the theories of dynamic strength and failure criterion of dynamic rock, the rock dynamic damage and the evolution of pre-crack seam were simulated by the tensile damage and shear failure of the model. According to the actual situation of No. 92 ore body test stope at Tongkeng Mine, the formation process of the pre-crack blast seam was simulated by Ansys/Ls-dyna software, the pre-crack seam was inspected by a system of digital panoramic borehole camera. The pre-crack seam was inspected by the system of digital panoramic borehole in the roof. The results of the numerical simulation and inspection show that in the line of centers of pre-hole, the minimum of the tensile stress reaches 20 MPa, which is much larger than 13.7 MPa of the dynamic tensile strength of rock. The minimum particle vibration velocity reaches 50 cm/s, which is greater than 30-40 cm/s of the allowable vibration velocity. It is demonstrated that the rock is destroyed near the center line and the pre-crack is successfully formed by the large diameters and large distances pre-crack holes in the roof.
文摘An efficient and real-time simulation method is proposed for the dynamic electromagnetic characteristics of cluster targets to meet the requirements of engineering practical applications.First,the coordinate transformation method is used to establish a geometric model of the observation scene,which is described by the azimuth angles and elevation angles of the radar in the target reference frame and the attitude angles of the target in the radar reference frame.Then,an approach for dynamic electromagnetic scattering simulation is proposed.Finally,a fast-computing method based on sparsity in the time domain,space domain,and frequency domain is proposed.The method analyzes the sparsity-based dynamic scattering characteristic of the typical cluster targets.The error between the sparsity-based method and the benchmark is small,proving the effectiveness of the proposed method.
