To enhance the computational efficiency of spatio-temporally discretized phase-field models,we present a high-speed solver specifically designed for the Poisson equations,a component frequently used in the numerical c...To enhance the computational efficiency of spatio-temporally discretized phase-field models,we present a high-speed solver specifically designed for the Poisson equations,a component frequently used in the numerical computation of such models.This efficient solver employs algorithms based on discrete cosine transformations(DCT)or discrete sine transformations(DST)and is not restricted by any spatio-temporal schemes.Our proposed methodology is appropriate for a variety of phase-field models and is especially efficient when combined with flow field systems.Meanwhile,this study has conducted an extensive numerical comparison and found that employing DCT and DST techniques not only yields results comparable to those obtained via the Multigrid(MG)method,a conventional approach used in the resolution of the Poisson equations,but also enhances computational efficiency by over 90%.展开更多
In polar regions, floating ice exhibits distinct characteristics across a range of spatial scales. It is well recognized that the irregular geometry of these ice formations markedly influences their dynamic behavior. ...In polar regions, floating ice exhibits distinct characteristics across a range of spatial scales. It is well recognized that the irregular geometry of these ice formations markedly influences their dynamic behavior. This study introduces a polyhedral Discrete Element Method (DEM) tailored for polar ice, incorporating the Gilbert-Johnson-Keerthi (GJK) and Expanding Polytope Algorithm (EPA) for contact detection. This approach facilitates the simulation of the drift and collision processes of floating ice, effectively capturing its freezing and fragmentation. Subsequently, the stability and reli ability of this model are validated by uniaxial compression on level ice fields, focusing specifically on the influence of compression strength on deformation resistance. Additionally, clusters of ice floes nav igating through narrow channels are simulated. These studies have qualitatively assessed the effects of Floe Size Distribution (FSD), initial concentration, and circularity on their flow dynamics. The higher power-law exponent values in the FSD, increased circularity, and decreased concentration are each as sociated with accelerated flow in ice floe fields. The simulation results distinctly demonstrate the con siderable impact of sea ice geometry on the movement of clusters, offering valuable insights into the complexities of polar ice dynamics.展开更多
Signal filtering and differential acquisition are classic yet challenging issues in control engineering.The discrete-time optimal control(DTOC)based on classic tracking differentiator(TD)can effectively extract differ...Signal filtering and differential acquisition are classic yet challenging issues in control engineering.The discrete-time optimal control(DTOC)based on classic tracking differentiator(TD)can effectively extract differentiation signals and filter signals,while eliminating the chattering problem that arises during the discretization of the continuous solution.However,under external disturbance,the convergence mode may change,leading to overshoot and noise amplification.In this paper,a dual-switching strategy is proposed,which can alternate between the base double-integral system and its dual system according to the quadrant of the system’s state.And a novel linearized control law is also introduced,deriving a novel dual-switch tracking differentiator.Further analysis of system convergence and time optimality is provided.Simulation results show that the application of this dual-switching strategy notably reduces overshoot in both tracking and differential signals while enhancing noise filtering performance.Moreover,experiments conducted on a permanent magnet synchronous motor(PMSM)platform,where the proposed TD acts as a filter in the speed feedback loop,demonstrate that the standard deviation between the reference speed and the target speed(at a constant speed of 378 r/min)decreased from 5.63 r/min to 4.93 r/min,compared to the moving average algorithm.展开更多
Outwash deposit is a unique type of geological materials, and its features such as heterogeneity, discontinuity and nonlinearity determine the complexity of mechanical characteristics and failure mechanism. In this wo...Outwash deposit is a unique type of geological materials, and its features such as heterogeneity, discontinuity and nonlinearity determine the complexity of mechanical characteristics and failure mechanism. In this work, random meso-structure of outwash deposits was constructed by the technique of computer random simulation based on characteristics of its meso-structure in the statistical sense and some simplifications, and a series of large direct shear tests on numerical samples of outwash deposits with stone contents of 15%, 30%, 45% and 60% were conducted using the discrete element method to further investigate its mechanical characteristics and failure mechanism under external load. The results show that the deformation characteristics and shear strength of outwash deposits are to some extent improved with the increase of stone content, and the shear stress–shear displacement curves of outwash deposits show great differences at the post-peak stage due to the random spatial distribution and content of stones. From the mesoscopic view, normal directions of contacts between "soil" and "stone" particles undergo apparent deflection as the shear displacement continues during the shearing process, accompanying redistribution of the magnitude of contact forces during the shearing process. For outwash deposits, the shear zone formed after shear failure is an irregular stripe due to the movements of stones near the shear zone, and it expands gradually with the increase of stone content. In addition, there is an approximately linear relation between the mean increment of internal friction angle and the stone content lying between 30% and 60%, and a concave nonlinear relation between the mean increment of cohesion and stone content, which are in good agreement with the existing research results.展开更多
With the development of high-speed and heavy-haul railway in China, problems like insufficient thickness of ballast bed and overlarge track stiffness are obvious. Ballast may break into small particles and their conta...With the development of high-speed and heavy-haul railway in China, problems like insufficient thickness of ballast bed and overlarge track stiffness are obvious. Ballast may break into small particles and their contact status will deteriorate under cyclic loading, resulting in ballast degradation. Discrete element method(DEM) was used to research improved performance of ballast bed using elastic sleeper. Clusters were generated by bonding spheres to model real ballasts, while broken bonds were utilized to distinguish breakage. Two kinds of ballast beds with elastic sleeper and conventional sleeper were established, respectively. After applying cyclic loading to the models, differences of mechanical properties between two models were analyzed by contrasting their dynamic behavior indexes, such as particle contact force, sleeper settlement, vibration velocity and acceleration, breakage characteristic. The results illustrate that compared with conventional sleeper, elastic sleeper increases sleeper settlement, while reduces ballast vibration and contact force between particles, which could depress ballast breakage.展开更多
The objective of this work is to model the microstructure of asphalt mixture and build virtual test for asphalt mixture by using Particle Flow Code in three dimensions(PFC^(3D))based on three-dimensional discrete elem...The objective of this work is to model the microstructure of asphalt mixture and build virtual test for asphalt mixture by using Particle Flow Code in three dimensions(PFC^(3D))based on three-dimensional discrete element method.A randomly generating algorithm was proposed to capture the three-dimensional irregular shape of coarse aggregate.And then,modeling algorithm and method for graded aggregates were built.Based on the combination of modeling of coarse aggregates,asphalt mastic and air voids,three-dimensional virtual sample of asphalt mixture was modeled by using PFC^(3D).Virtual tests for penetration test of aggregate and uniaxial creep test of asphalt mixture were built and conducted by using PFC^(3D).By comparison of the testing results between virtual tests and actual laboratory tests,the validity of the microstructure modeling and virtual test built in this study was verified.Additionally,compared with laboratory test,the virtual test is easier to conduct and has less variability.It is proved that microstructure modeling and virtual test based on three-dimensional discrete element method is a promising way to conduct research of asphalt mixture.展开更多
A user-defined micromechanical model was developed to investigate the fracture mechanism of asphalt concrete （AC） using the discrete element method （DEM）. A three-dimensional （3D） AC beam was built using the ＂F...A user-defined micromechanical model was developed to investigate the fracture mechanism of asphalt concrete （AC） using the discrete element method （DEM）. A three-dimensional （3D） AC beam was built using the ＂Fish＂ language provided by PFC3D and was employed to simulate the three-point bending beam test at two temperature levels： -10 ℃ and 15℃. The AC beam was modeled with the consideration of the microstructural features of asphalt mixtures. Uniaxial complex modulus test and indirect tensile strength test were conducted to obtain material input parameters for numerical modeling. The 3D predictions were validated using laboratory experimental measurements of AC beams prepared by the same mixture design. Effects of mastic stiffness, cohesive and adhesive strength on AC fracture behavior were investigated using the DEM model. The results show that the 3D DEM fracture model can accurately predict the fracture patterns of asphalt concrete. The ratio of stress at interfaces to the stress in mastics increases as the mastic stiffness decreases; however, the increase in the cohesive strength or adhesive strength shows no significant influence on the tensile strength.展开更多
The behavior of the active layer of material bed within rotary kilns plays a key role in industrial applications.To obtain its influences on industrial process,different regimes of particle motion have been simulated ...The behavior of the active layer of material bed within rotary kilns plays a key role in industrial applications.To obtain its influences on industrial process,different regimes of particle motion have been simulated by discrete element method(DEM) in three dimensions under variant rotation speeds,filling degree,based on the background of induration process of iron ore pellets.The influences of the mentioned factors on the maximum thickness of the active layer and the average velocity of particles have been investigated.The average velocity of particles increases with Froude number following the power function over a wide range,and the maximum thickness rises with increasing rotation speed in a way of logarithm.The influence of the filling degree f on the maximum thickness exhibits a good linearity under two classic regimes,but the increasing of the average velocity of the active layer is limited at f=0.4.This basic research highlights the impact of the active layer within rotary kilns,and lays a good foundation for the further investigation in mixing and heat transfer within the particle bed inside rotary kilns.展开更多
The mechanical properties of outwash deposits which are taken as unconsolidated geo-materials with the characteristics of non-uniformity, heterogeneity and multiphase have attracted much attention in engineering. Acco...The mechanical properties of outwash deposits which are taken as unconsolidated geo-materials with the characteristics of non-uniformity, heterogeneity and multiphase have attracted much attention in engineering. According to the results of laboratory direct shear test on the remolded samples, the soil particle parameters of numerical model based on in-situ particle size cumulative curves and 3D granular discrete element method were determined. Then, numerical experiments on different lithology, stone content and gradation composition were conducted. The results show that it is not a flat surface but a shear band that yields in the sample. The curve of particle velocity vs distance from the designed shear surface of test model that is taken as a datum plane in the vertical section of sample shows in "S" shape. The shear disturbance area is about twice the maximum diameter of stone blocks. The greater the stiffness of stone is, the rougher the shear surface is. The shear strength of outwash deposits is largely controlled by lithology and stone content, and the bite force between stone blocks is the root reason of larger friction angle. It is also shown that strain hardening and low shear dilatancy occur under high confining pressure as well as possibility of shear shrinkage. But it is easy to behave shear dilatation and strain softening under low confining pressure. The relationship between particle frictional coefficient and stone content presents an approximately quadratic parabola increase. The strain energy first increases and then drops with the increase of frictional energy. The cohesion increases with soil stiffness increasing but decreases with stone stiffness increasing. Numerical results are consistent with the laboratory test results of remolded samples, which indicate that this method can be a beneficial supplement to determine the parameters of engineering deposit bodies.展开更多
The lateral resistance of sleeper plays an important role in ensuring the stability of a railway track, which may change in the operation of railway, due to the fouling in the ballast bed. In this work, discrete eleme...The lateral resistance of sleeper plays an important role in ensuring the stability of a railway track, which may change in the operation of railway, due to the fouling in the ballast bed. In this work, discrete element method was adopted to investigate the effect of fouling on the lateral resistance of sleeper. The shape information of ballast was captured by method of three-dimensional vision reconstruction. In order to calibrate the mechanical parameters and verify the models, a lateral resistance field test was carried out by using a custom-made device. The contact force distributions in the different parts of sleeper as well as the interaction between ballast and sleeper were discussed in depth. The results show that fouling of ballast bed evidently reduces the lateral resistance of sleeper and the decreasing degree is also related to the fouled position of ballast bed, in the order of shoulder > bottom > side.Therefore, the effect of fouling, especially the fouling in the ballast shoulder, on the lateral resistance of sleeper, should be taken into account in ballast track maintenance work.展开更多
The adhesion between the mining machine and the deep-sea sediments will significantly affect the driving performance of the mining machine in the deep-sea environment.When the mining machine and the deep-sea sediment ...The adhesion between the mining machine and the deep-sea sediments will significantly affect the driving performance of the mining machine in the deep-sea environment.When the mining machine and the deep-sea sediment interaction simulation was carried out,the accuracy of the particle interaction parameters will directly affect the simulation results.This study proposed a method to systematically calibrate the interaction parameters between deep-sea sediment and grouser through the combination of experiment and simulation.The uniaxial compression test and macro adhesion test and corresponding discrete element numerical simulation were carried out,modifying the contact parameters until the simulation results are close to the experimental results.Then the micro-parameters of the JKR adhesion contact model were back calibrated with the test results,and the contact parameters between soil particle-soil particle and soil particle-metal are calibrated.Besides,the adhesion test shows that the adhesion forces were ranked in the order of 5052<STi80<TA2<TC4 under the same surface roughness,which indicates the aluminum alloy 5052 has the best anti-adhesion performance.The relationship between surface adhesion force and microscopic contact parameters was studied by discrete element numerical simulation,and the result shows that the coefficient of static friction and the coefficient of rolling friction has little effect on adhesion force.While it is mainly affected by the coefficient of restitution and surface energy,the surface adhesion force tends to decrease with the increase of the coefficient of restitution and increase with the growth of surface energy.The obtained parameters of soil particle to soil particle and soil particle to metal affecting the adhesion can contribute to the design optimization for the grouser of mining machines to decrease surface adhesion and enhance its movability and mining efficiency.展开更多
This study investigated the physicochemical properties,enzyme activities,volatile flavor components,microbial communities,and sensory evaluation of high-temperature Daqu(HTD)during the maturation process,and a standar...This study investigated the physicochemical properties,enzyme activities,volatile flavor components,microbial communities,and sensory evaluation of high-temperature Daqu(HTD)during the maturation process,and a standard system was established for comprehensive quality evaluation of HTD.There were obvious changes in the physicochemical properties,enzyme activities,and volatile flavor components at different storage periods,which affected the sensory evaluation of HTD to a certain extent.The results of high-throughput sequencing revealed significant microbial diversity,and showed that the bacterial community changed significantly more than did the fungal community.During the storage process,the dominant bacterial genera were Kroppenstedtia and Thermoascus.The correlation between dominant microorganisms and quality indicators highlighted their role in HTD quality.Lactococcus,Candida,Pichia,Paecilomyces,and protease activity played a crucial role in the formation of isovaleraldehyde.Acidic protease activity had the greatest impact on the microbial community.Moisture promoted isobutyric acid generation.Furthermore,the comprehensive quality evaluation standard system was established by the entropy weight method combined with multi-factor fuzzy mathematics.Consequently,this study provides innovative insights for comprehensive quality evaluation of HTD during storage and establishes a groundwork for scientific and rational storage of HTD and quality control of sauce-flavor Baijiu.展开更多
This paper discuss stability of the full discrete nonlinear Galerkin method based on the approximation inertial manifold method for some nonlinear evolution equation, for example, some nonlinear reactor equation and N...This paper discuss stability of the full discrete nonlinear Galerkin method based on the approximation inertial manifold method for some nonlinear evolution equation, for example, some nonlinear reactor equation and Navier-Stokes Equation. In the paper we provide some necessary and sufficient conditions of stability.展开更多
Optimization of design features of reinforced sheet is investigated. Initially, equations governing composite structures are extracted based on Kirchhoff sheet model under bending using Hamilton's principal. Then,...Optimization of design features of reinforced sheet is investigated. Initially, equations governing composite structures are extracted based on Kirchhoff sheet model under bending using Hamilton's principal. Then, design parameters for the composite structure are extracted with simple supportive boundary conditions from proposed solution. Next, optimization is achieved by determining dimensions of a reinforced sheet specimen. Weight optimization of reinforced sheet structure has been obtained based on variations in thickness and number of longitudinal and transverse reinforcements. Buckling static characteristic is utilized in optimization process. To solve the extracted equations, semi-analytical method of CS-DSG3 has been applied. Results are presented in graphs that show variation of design parameters by changing the geometric parameters. ABAQUS software has been used for design verification. The results show that an increase in thickness of 3 mm skip value tends to be zero. Also, there is a change in the amount of deflection for sheets with a minimum thickness of 3 mm by increasing the number of longitudinal and transverse reinforcement. There is a good agreement between the numerical method of finite elements and the method X-FEM-DSG3.展开更多
Shock wave caused by a sudden release of high-energy,such as explosion and blast,usually affects a significant range of areas.The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise re...Shock wave caused by a sudden release of high-energy,such as explosion and blast,usually affects a significant range of areas.The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise results is inefficient in terms of computational resource.This is particularly evident when large-scale fluid field simulations are conducted with significant differences in computational domain size.In this work,a variable-domain-size adaptive mesh enlargement(vAME)method is developed based on the proposed adaptive mesh enlargement(AME)method for modeling multi-explosives explosion problems.The vAME method reduces the division of numerous empty areas or unnecessary computational domains by adaptively suspending enlargement operation in one or two directions,rather than in all directions as in AME method.A series of numerical tests via AME and vAME with varying nonintegral enlargement ratios and different mesh numbers are simulated to verify the efficiency and order of accuracy.An estimate of speedup ratio is analyzed for further efficiency comparison.Several large-scale near-ground explosion experiments with single/multiple explosives are performed to analyze the shock wave superposition formed by the incident wave,reflected wave,and Mach wave.Additionally,the vAME method is employed to validate the accuracy,as well as to investigate the performance of the fluid field and shock wave propagation,considering explosive quantities ranging from 1 to 5 while maintaining a constant total mass.The results show a satisfactory correlation between the overpressure versus time curves for experiments and numerical simulations.The vAME method yields a competitive efficiency,increasing the computational speed to 3.0 and approximately 120,000 times in comparison to AME and the fully fine mesh method,respectively.It indicates that the vAME method reduces the computational cost with minimal impact on the results for such large-scale high-energy release problems with significant differences in computational domain size.展开更多
In this paper,a new technique is introduced to construct higher-order iterative methods for solving nonlinear systems.The order of convergence of some iterative methods can be improved by three at the cost of introduc...In this paper,a new technique is introduced to construct higher-order iterative methods for solving nonlinear systems.The order of convergence of some iterative methods can be improved by three at the cost of introducing only one additional evaluation of the function in each step.Furthermore,some new efficient methods with a higher-order of convergence are obtained by using only a single matrix inversion in each iteration.Analyses of convergence properties and computational efficiency of these new methods are made and testified by several numerical problems.By comparison,the new schemes are more efficient than the corresponding existing ones,particularly for large problem sizes.展开更多
To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based sim...To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based simulation(NNS)method with higher accuracy and better efficiency was proposed.The NNS method consisted of three main steps.First,the parameters of blast loads,including the peak pressures and impulses of cylindrical charges with different aspect ratios(L/D)at different stand-off distances and incident angles were obtained by two-dimensional numerical simulations.Subsequently,incident shape factors of cylindrical charges with arbitrary aspect ratios were predicted by a neural network.Finally,reflected shape factors were derived and implemented into the subroutine of the ABAQUS code to modify the CONWEP model,including modifications of impulse and overpressure.The reliability of the proposed NNS method was verified by related experimental results.Remarkable accuracy improvement was acquired by the proposed NNS method compared with the unmodified CONWEP model.Moreover,huge efficiency superiority was obtained by the proposed NNS method compared with the CEL method.The proposed NNS method showed good accuracy when the scaled distance was greater than 0.2 m/kg^(1/3).It should be noted that there is no need to generate a new dataset again since the blast loads satisfy the similarity law,and the proposed NNS method can be directly used to simulate the blast loads generated by different cylindrical charges.The proposed NNS method with high efficiency and accuracy can be used as an effective method to analyze the dynamic response of structures under blast loads,and it has significant application prospects in designing protective structures.展开更多
Red-green-blue(RGB)beam combiners are widely used in scenarios such as augmented reality/virtual reality(AR/VR),laser projection,biochemical detection,and other fields.Optical waveguide combiners have attracted extens...Red-green-blue(RGB)beam combiners are widely used in scenarios such as augmented reality/virtual reality(AR/VR),laser projection,biochemical detection,and other fields.Optical waveguide combiners have attracted extensive attention due to their advantages of small size,high multiplexing efficiency,convenient mass production,and low cost.An RGB beam combiner based on directional couplers is designed,with a core-cladding relative refractive index difference of 0.75%.The RGB beam combiner is optimized from the perspective of parameter optimization.Using the beam propagation method(BPM),the relationship between the performance of the RGB beam combiner and individual parameters is studied,achieving preliminary optimization of the device’s performance.The key parameters of the RGB beam combiner are optimized using the entropy weight-technique for order preference by similarity to an ideal solution TOPSIS method,establishing the optimal parameter scheme and further improving the device’s performance indicators.The results show that after optimization,the multiplexing efficiencies for red,green,and blue lights,as well as the average multiplexing efficiency,reached 99.17%,99.76%,96.63%and 98.52%,respectively.The size of the RGB beam combiner is 4.768 mm×0.062 mm.展开更多
The surrounding rock is prone to large-scale loosening and failure after the excavation of shallow large-span caverns because of the thin overlying strata and large cross-section span.The rational design of bolt suppo...The surrounding rock is prone to large-scale loosening and failure after the excavation of shallow large-span caverns because of the thin overlying strata and large cross-section span.The rational design of bolt support is very important to the safety control of surrounding rock as a common support means.The control mechanism and design method of bolt support for shallow-buried large-span caverns is carried out.The calculation method of bolt prestress and length based on arched failure and collapsed failure mode is established.The influence mechanism of different influencing factors on the bolt prestress and length is clarified.At the same time,the constant resistance energy-absorbing bolt with high strength and high toughness is developed,and the comparative test of mechanical properties is carried out.On this basis,the design method of high prestressed bolt support for shallow-buried large-span caverns is put forward,and the field test is carried out in Qingdao metro station in China.The monitoring results show that the maximum roof settlement is 6.8 mm after the new design method is adopted,and the effective control of the shallow-buried large-span caverns is realized.The research results can provide theoretical and technical support for the safety control of shallow-buried large-span caverns.展开更多
To study the influence of silicon(Si)on 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane(CL-20),NC/CL-20 composite explosives and Si/NC/CL-20 composite explosives were prepared by the electrostatic spraying ...To study the influence of silicon(Si)on 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane(CL-20),NC/CL-20 composite explosives and Si/NC/CL-20 composite explosives were prepared by the electrostatic spraying method.The morphology,structure and thermal decomposition properties of the samples were analyzed using scanning electron microscopy(SEM),X-ray energy spectroscopy(EDS),infrared spectroscopy(FT-IR),and simultaneous thermal analyzer(TG-DSC).Additionally,the combustion process of the samples was tested using a high-speed camera.The results show that the addition of nano-Si contributes to the formation of composite explosives with regular morphology and smaller particle size.The Si/NC/CL-20 composite explosive has better and more uniform sphericity,with an average particle size of 73.4 nm,compared to the NC/CL-20 composite explosive.The Si/NC/CL-20 composite explosive which produced by the electrostatic spraying method,achieves physically uniform distribution of the components including NC,CL-20,Si.The addition of Si promotes the thermal decomposition of CL-20.In comparison to the NC/CL-20 composite explosive,the activation energy of the Si/NC/CL-20 composite explosive decreases by 16.78 kJ/mol,and the self-accelerated decomposition temperature and the critical temperature of thermal explosion decreases by 3.12 K and 2.61 K,respectively.Furthermore,Si/NC/CL-20 composite explosive has shorter ignition delay time and faster combustion rate compared to the NC/CL-20 composite explosive,which shows that Si can improve the combustion performance of CL-20.展开更多
基金Supported by Shanxi Province Natural Science Research(202203021212249)Special/Youth Foundation of Taiyuan University of Technology(2022QN101)+3 种基金National Natural Science Foundation of China(12301556)Research Project Supported by Shanxi Scholarship Council of China(2021-029)International Cooperation Base and Platform Project of Shanxi Province(202104041101019)Basic Research Plan of Shanxi Province(202203021211129)。
文摘To enhance the computational efficiency of spatio-temporally discretized phase-field models,we present a high-speed solver specifically designed for the Poisson equations,a component frequently used in the numerical computation of such models.This efficient solver employs algorithms based on discrete cosine transformations(DCT)or discrete sine transformations(DST)and is not restricted by any spatio-temporal schemes.Our proposed methodology is appropriate for a variety of phase-field models and is especially efficient when combined with flow field systems.Meanwhile,this study has conducted an extensive numerical comparison and found that employing DCT and DST techniques not only yields results comparable to those obtained via the Multigrid(MG)method,a conventional approach used in the resolution of the Poisson equations,but also enhances computational efficiency by over 90%.
文摘In polar regions, floating ice exhibits distinct characteristics across a range of spatial scales. It is well recognized that the irregular geometry of these ice formations markedly influences their dynamic behavior. This study introduces a polyhedral Discrete Element Method (DEM) tailored for polar ice, incorporating the Gilbert-Johnson-Keerthi (GJK) and Expanding Polytope Algorithm (EPA) for contact detection. This approach facilitates the simulation of the drift and collision processes of floating ice, effectively capturing its freezing and fragmentation. Subsequently, the stability and reli ability of this model are validated by uniaxial compression on level ice fields, focusing specifically on the influence of compression strength on deformation resistance. Additionally, clusters of ice floes nav igating through narrow channels are simulated. These studies have qualitatively assessed the effects of Floe Size Distribution (FSD), initial concentration, and circularity on their flow dynamics. The higher power-law exponent values in the FSD, increased circularity, and decreased concentration are each as sociated with accelerated flow in ice floe fields. The simulation results distinctly demonstrate the con siderable impact of sea ice geometry on the movement of clusters, offering valuable insights into the complexities of polar ice dynamics.
基金Project(QZKFKT2023-012)supported by the State Key Laboratory of Heavy-duty and Express High-power Electric Locomotive,China。
文摘Signal filtering and differential acquisition are classic yet challenging issues in control engineering.The discrete-time optimal control(DTOC)based on classic tracking differentiator(TD)can effectively extract differentiation signals and filter signals,while eliminating the chattering problem that arises during the discretization of the continuous solution.However,under external disturbance,the convergence mode may change,leading to overshoot and noise amplification.In this paper,a dual-switching strategy is proposed,which can alternate between the base double-integral system and its dual system according to the quadrant of the system’s state.And a novel linearized control law is also introduced,deriving a novel dual-switch tracking differentiator.Further analysis of system convergence and time optimality is provided.Simulation results show that the application of this dual-switching strategy notably reduces overshoot in both tracking and differential signals while enhancing noise filtering performance.Moreover,experiments conducted on a permanent magnet synchronous motor(PMSM)platform,where the proposed TD acts as a filter in the speed feedback loop,demonstrate that the standard deviation between the reference speed and the target speed(at a constant speed of 378 r/min)decreased from 5.63 r/min to 4.93 r/min,compared to the moving average algorithm.
基金Project(2011CB013504) supported by the National Basic Research Program(973 Program)of ChinaProject(2013BAB06B01) supported by the National Science&Technology Pillar Program during the Twelfth Five-year Plan Period+2 种基金Projects(11772118,51479049,51709282) supported by the National Natural Science Foundation of ChinaProject(2017M620838) supported by the Postdoctoral Science Foundation of ChinaProject(487237) supported by the Natural Sciences and Engineering Research Council of Canada
文摘Outwash deposit is a unique type of geological materials, and its features such as heterogeneity, discontinuity and nonlinearity determine the complexity of mechanical characteristics and failure mechanism. In this work, random meso-structure of outwash deposits was constructed by the technique of computer random simulation based on characteristics of its meso-structure in the statistical sense and some simplifications, and a series of large direct shear tests on numerical samples of outwash deposits with stone contents of 15%, 30%, 45% and 60% were conducted using the discrete element method to further investigate its mechanical characteristics and failure mechanism under external load. The results show that the deformation characteristics and shear strength of outwash deposits are to some extent improved with the increase of stone content, and the shear stress–shear displacement curves of outwash deposits show great differences at the post-peak stage due to the random spatial distribution and content of stones. From the mesoscopic view, normal directions of contacts between "soil" and "stone" particles undergo apparent deflection as the shear displacement continues during the shearing process, accompanying redistribution of the magnitude of contact forces during the shearing process. For outwash deposits, the shear zone formed after shear failure is an irregular stripe due to the movements of stones near the shear zone, and it expands gradually with the increase of stone content. In addition, there is an approximately linear relation between the mean increment of internal friction angle and the stone content lying between 30% and 60%, and a concave nonlinear relation between the mean increment of cohesion and stone content, which are in good agreement with the existing research results.
基金Project(U1234211)supported by the National Natural Science Foundation of ChinaProject(2013G009-B)supported by China Railway Corporation
文摘With the development of high-speed and heavy-haul railway in China, problems like insufficient thickness of ballast bed and overlarge track stiffness are obvious. Ballast may break into small particles and their contact status will deteriorate under cyclic loading, resulting in ballast degradation. Discrete element method(DEM) was used to research improved performance of ballast bed using elastic sleeper. Clusters were generated by bonding spheres to model real ballasts, while broken bonds were utilized to distinguish breakage. Two kinds of ballast beds with elastic sleeper and conventional sleeper were established, respectively. After applying cyclic loading to the models, differences of mechanical properties between two models were analyzed by contrasting their dynamic behavior indexes, such as particle contact force, sleeper settlement, vibration velocity and acceleration, breakage characteristic. The results illustrate that compared with conventional sleeper, elastic sleeper increases sleeper settlement, while reduces ballast vibration and contact force between particles, which could depress ballast breakage.
基金Project(51378006) supported by National Natural Science Foundation of ChinaProject(141076) supported by Huoyingdong Foundation of the Ministry of Education of China+1 种基金Project(2242015R30027) supported by Excellent Young Teacher Program of Southeast University,ChinaProject(BK20140109) supported by the Natural Science Foundation of Jiangsu Province,China
文摘The objective of this work is to model the microstructure of asphalt mixture and build virtual test for asphalt mixture by using Particle Flow Code in three dimensions(PFC^(3D))based on three-dimensional discrete element method.A randomly generating algorithm was proposed to capture the three-dimensional irregular shape of coarse aggregate.And then,modeling algorithm and method for graded aggregates were built.Based on the combination of modeling of coarse aggregates,asphalt mastic and air voids,three-dimensional virtual sample of asphalt mixture was modeled by using PFC^(3D).Virtual tests for penetration test of aggregate and uniaxial creep test of asphalt mixture were built and conducted by using PFC^(3D).By comparison of the testing results between virtual tests and actual laboratory tests,the validity of the microstructure modeling and virtual test built in this study was verified.Additionally,compared with laboratory test,the virtual test is easier to conduct and has less variability.It is proved that microstructure modeling and virtual test based on three-dimensional discrete element method is a promising way to conduct research of asphalt mixture.
基金Project(51208178)supported by the National Natural Science Foundation of ChinaProject(2012M520991)supported by China Postdoctoral Science Foundation
文摘A user-defined micromechanical model was developed to investigate the fracture mechanism of asphalt concrete （AC） using the discrete element method （DEM）. A three-dimensional （3D） AC beam was built using the ＂Fish＂ language provided by PFC3D and was employed to simulate the three-point bending beam test at two temperature levels： -10 ℃ and 15℃. The AC beam was modeled with the consideration of the microstructural features of asphalt mixtures. Uniaxial complex modulus test and indirect tensile strength test were conducted to obtain material input parameters for numerical modeling. The 3D predictions were validated using laboratory experimental measurements of AC beams prepared by the same mixture design. Effects of mastic stiffness, cohesive and adhesive strength on AC fracture behavior were investigated using the DEM model. The results show that the 3D DEM fracture model can accurately predict the fracture patterns of asphalt concrete. The ratio of stress at interfaces to the stress in mastics increases as the mastic stiffness decreases; however, the increase in the cohesive strength or adhesive strength shows no significant influence on the tensile strength.
基金Project(FRF-AS-10-0058) supported by the Fundamental Research Funds for the Central Universities,China
文摘The behavior of the active layer of material bed within rotary kilns plays a key role in industrial applications.To obtain its influences on industrial process,different regimes of particle motion have been simulated by discrete element method(DEM) in three dimensions under variant rotation speeds,filling degree,based on the background of induration process of iron ore pellets.The influences of the mentioned factors on the maximum thickness of the active layer and the average velocity of particles have been investigated.The average velocity of particles increases with Froude number following the power function over a wide range,and the maximum thickness rises with increasing rotation speed in a way of logarithm.The influence of the filling degree f on the maximum thickness exhibits a good linearity under two classic regimes,but the increasing of the average velocity of the active layer is limited at f=0.4.This basic research highlights the impact of the active layer within rotary kilns,and lays a good foundation for the further investigation in mixing and heat transfer within the particle bed inside rotary kilns.
基金Project(2011CB013504) supported by the National Basic Research Program of ChinaProjects(50911130366, 11172090) supported by the National Natural Science Foundation of ChinaProject supported by Central University Basic Research Special Fund, China
文摘The mechanical properties of outwash deposits which are taken as unconsolidated geo-materials with the characteristics of non-uniformity, heterogeneity and multiphase have attracted much attention in engineering. According to the results of laboratory direct shear test on the remolded samples, the soil particle parameters of numerical model based on in-situ particle size cumulative curves and 3D granular discrete element method were determined. Then, numerical experiments on different lithology, stone content and gradation composition were conducted. The results show that it is not a flat surface but a shear band that yields in the sample. The curve of particle velocity vs distance from the designed shear surface of test model that is taken as a datum plane in the vertical section of sample shows in "S" shape. The shear disturbance area is about twice the maximum diameter of stone blocks. The greater the stiffness of stone is, the rougher the shear surface is. The shear strength of outwash deposits is largely controlled by lithology and stone content, and the bite force between stone blocks is the root reason of larger friction angle. It is also shown that strain hardening and low shear dilatancy occur under high confining pressure as well as possibility of shear shrinkage. But it is easy to behave shear dilatation and strain softening under low confining pressure. The relationship between particle frictional coefficient and stone content presents an approximately quadratic parabola increase. The strain energy first increases and then drops with the increase of frictional energy. The cohesion increases with soil stiffness increasing but decreases with stone stiffness increasing. Numerical results are consistent with the laboratory test results of remolded samples, which indicate that this method can be a beneficial supplement to determine the parameters of engineering deposit bodies.
基金Projects(U1234211,61472029,51208034)supported by the National Natural Science Foundation of China
文摘The lateral resistance of sleeper plays an important role in ensuring the stability of a railway track, which may change in the operation of railway, due to the fouling in the ballast bed. In this work, discrete element method was adopted to investigate the effect of fouling on the lateral resistance of sleeper. The shape information of ballast was captured by method of three-dimensional vision reconstruction. In order to calibrate the mechanical parameters and verify the models, a lateral resistance field test was carried out by using a custom-made device. The contact force distributions in the different parts of sleeper as well as the interaction between ballast and sleeper were discussed in depth. The results show that fouling of ballast bed evidently reduces the lateral resistance of sleeper and the decreasing degree is also related to the fouled position of ballast bed, in the order of shoulder > bottom > side.Therefore, the effect of fouling, especially the fouling in the ballast shoulder, on the lateral resistance of sleeper, should be taken into account in ballast track maintenance work.
基金Project(12072309)supported by the National Natural Science Foundation of ChinaProject(19B546)supported by the Education Department Foundation of Hunan Province,ChinaProject(2019RS1059)supported by the Hunan Innovative Province Construction Project,China。
文摘The adhesion between the mining machine and the deep-sea sediments will significantly affect the driving performance of the mining machine in the deep-sea environment.When the mining machine and the deep-sea sediment interaction simulation was carried out,the accuracy of the particle interaction parameters will directly affect the simulation results.This study proposed a method to systematically calibrate the interaction parameters between deep-sea sediment and grouser through the combination of experiment and simulation.The uniaxial compression test and macro adhesion test and corresponding discrete element numerical simulation were carried out,modifying the contact parameters until the simulation results are close to the experimental results.Then the micro-parameters of the JKR adhesion contact model were back calibrated with the test results,and the contact parameters between soil particle-soil particle and soil particle-metal are calibrated.Besides,the adhesion test shows that the adhesion forces were ranked in the order of 5052<STi80<TA2<TC4 under the same surface roughness,which indicates the aluminum alloy 5052 has the best anti-adhesion performance.The relationship between surface adhesion force and microscopic contact parameters was studied by discrete element numerical simulation,and the result shows that the coefficient of static friction and the coefficient of rolling friction has little effect on adhesion force.While it is mainly affected by the coefficient of restitution and surface energy,the surface adhesion force tends to decrease with the increase of the coefficient of restitution and increase with the growth of surface energy.The obtained parameters of soil particle to soil particle and soil particle to metal affecting the adhesion can contribute to the design optimization for the grouser of mining machines to decrease surface adhesion and enhance its movability and mining efficiency.
文摘This study investigated the physicochemical properties,enzyme activities,volatile flavor components,microbial communities,and sensory evaluation of high-temperature Daqu(HTD)during the maturation process,and a standard system was established for comprehensive quality evaluation of HTD.There were obvious changes in the physicochemical properties,enzyme activities,and volatile flavor components at different storage periods,which affected the sensory evaluation of HTD to a certain extent.The results of high-throughput sequencing revealed significant microbial diversity,and showed that the bacterial community changed significantly more than did the fungal community.During the storage process,the dominant bacterial genera were Kroppenstedtia and Thermoascus.The correlation between dominant microorganisms and quality indicators highlighted their role in HTD quality.Lactococcus,Candida,Pichia,Paecilomyces,and protease activity played a crucial role in the formation of isovaleraldehyde.Acidic protease activity had the greatest impact on the microbial community.Moisture promoted isobutyric acid generation.Furthermore,the comprehensive quality evaluation standard system was established by the entropy weight method combined with multi-factor fuzzy mathematics.Consequently,this study provides innovative insights for comprehensive quality evaluation of HTD during storage and establishes a groundwork for scientific and rational storage of HTD and quality control of sauce-flavor Baijiu.
文摘This paper discuss stability of the full discrete nonlinear Galerkin method based on the approximation inertial manifold method for some nonlinear evolution equation, for example, some nonlinear reactor equation and Navier-Stokes Equation. In the paper we provide some necessary and sufficient conditions of stability.
文摘Optimization of design features of reinforced sheet is investigated. Initially, equations governing composite structures are extracted based on Kirchhoff sheet model under bending using Hamilton's principal. Then, design parameters for the composite structure are extracted with simple supportive boundary conditions from proposed solution. Next, optimization is achieved by determining dimensions of a reinforced sheet specimen. Weight optimization of reinforced sheet structure has been obtained based on variations in thickness and number of longitudinal and transverse reinforcements. Buckling static characteristic is utilized in optimization process. To solve the extracted equations, semi-analytical method of CS-DSG3 has been applied. Results are presented in graphs that show variation of design parameters by changing the geometric parameters. ABAQUS software has been used for design verification. The results show that an increase in thickness of 3 mm skip value tends to be zero. Also, there is a change in the amount of deflection for sheets with a minimum thickness of 3 mm by increasing the number of longitudinal and transverse reinforcement. There is a good agreement between the numerical method of finite elements and the method X-FEM-DSG3.
基金supported by the National Natural Science Foundation of China(Grant Nos.12302435 and 12221002)。
文摘Shock wave caused by a sudden release of high-energy,such as explosion and blast,usually affects a significant range of areas.The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise results is inefficient in terms of computational resource.This is particularly evident when large-scale fluid field simulations are conducted with significant differences in computational domain size.In this work,a variable-domain-size adaptive mesh enlargement(vAME)method is developed based on the proposed adaptive mesh enlargement(AME)method for modeling multi-explosives explosion problems.The vAME method reduces the division of numerous empty areas or unnecessary computational domains by adaptively suspending enlargement operation in one or two directions,rather than in all directions as in AME method.A series of numerical tests via AME and vAME with varying nonintegral enlargement ratios and different mesh numbers are simulated to verify the efficiency and order of accuracy.An estimate of speedup ratio is analyzed for further efficiency comparison.Several large-scale near-ground explosion experiments with single/multiple explosives are performed to analyze the shock wave superposition formed by the incident wave,reflected wave,and Mach wave.Additionally,the vAME method is employed to validate the accuracy,as well as to investigate the performance of the fluid field and shock wave propagation,considering explosive quantities ranging from 1 to 5 while maintaining a constant total mass.The results show a satisfactory correlation between the overpressure versus time curves for experiments and numerical simulations.The vAME method yields a competitive efficiency,increasing the computational speed to 3.0 and approximately 120,000 times in comparison to AME and the fully fine mesh method,respectively.It indicates that the vAME method reduces the computational cost with minimal impact on the results for such large-scale high-energy release problems with significant differences in computational domain size.
基金Supported by the National Natural Science Foundation of China(12061048)NSF of Jiangxi Province(20232BAB201026,20232BAB201018)。
文摘In this paper,a new technique is introduced to construct higher-order iterative methods for solving nonlinear systems.The order of convergence of some iterative methods can be improved by three at the cost of introducing only one additional evaluation of the function in each step.Furthermore,some new efficient methods with a higher-order of convergence are obtained by using only a single matrix inversion in each iteration.Analyses of convergence properties and computational efficiency of these new methods are made and testified by several numerical problems.By comparison,the new schemes are more efficient than the corresponding existing ones,particularly for large problem sizes.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.52271317 and 52071149)the Fundamental Research Funds for the Central Universities(HUST:2019kfy XJJS007)。
文摘To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based simulation(NNS)method with higher accuracy and better efficiency was proposed.The NNS method consisted of three main steps.First,the parameters of blast loads,including the peak pressures and impulses of cylindrical charges with different aspect ratios(L/D)at different stand-off distances and incident angles were obtained by two-dimensional numerical simulations.Subsequently,incident shape factors of cylindrical charges with arbitrary aspect ratios were predicted by a neural network.Finally,reflected shape factors were derived and implemented into the subroutine of the ABAQUS code to modify the CONWEP model,including modifications of impulse and overpressure.The reliability of the proposed NNS method was verified by related experimental results.Remarkable accuracy improvement was acquired by the proposed NNS method compared with the unmodified CONWEP model.Moreover,huge efficiency superiority was obtained by the proposed NNS method compared with the CEL method.The proposed NNS method showed good accuracy when the scaled distance was greater than 0.2 m/kg^(1/3).It should be noted that there is no need to generate a new dataset again since the blast loads satisfy the similarity law,and the proposed NNS method can be directly used to simulate the blast loads generated by different cylindrical charges.The proposed NNS method with high efficiency and accuracy can be used as an effective method to analyze the dynamic response of structures under blast loads,and it has significant application prospects in designing protective structures.
基金Project(52175445)supported by the National Natural Science Foundation of ChinaProject(2022JJ30743)supported by the Natural Science Foundation of Hunan Province,China+1 种基金Project(2023GK2024)supported by the Key Research and Development Program of Hunan Province,ChinaProject(2023ZZTS0391)supported by the Fundamental Research Funds for the Central Universities of China。
文摘Red-green-blue(RGB)beam combiners are widely used in scenarios such as augmented reality/virtual reality(AR/VR),laser projection,biochemical detection,and other fields.Optical waveguide combiners have attracted extensive attention due to their advantages of small size,high multiplexing efficiency,convenient mass production,and low cost.An RGB beam combiner based on directional couplers is designed,with a core-cladding relative refractive index difference of 0.75%.The RGB beam combiner is optimized from the perspective of parameter optimization.Using the beam propagation method(BPM),the relationship between the performance of the RGB beam combiner and individual parameters is studied,achieving preliminary optimization of the device’s performance.The key parameters of the RGB beam combiner are optimized using the entropy weight-technique for order preference by similarity to an ideal solution TOPSIS method,establishing the optimal parameter scheme and further improving the device’s performance indicators.The results show that after optimization,the multiplexing efficiencies for red,green,and blue lights,as well as the average multiplexing efficiency,reached 99.17%,99.76%,96.63%and 98.52%,respectively.The size of the RGB beam combiner is 4.768 mm×0.062 mm.
基金Project(2023YFC3805700) supported by the National Key Research and Development Program of ChinaProjects(42477166,42277174) supported by the National Natural Science Foundation of China+2 种基金Project(2024JCCXSB01) supported by the Fundamental Research Funds for the Central Universities,ChinaProject(KFJJ24-01M) supported by the State Key Laboratory of Explosion Science and Safety Protection,Beijing Institute of Technology,ChinaProject(HLCX-2024-04) supported by the Open Foundation of Collaborative Innovation Center of Green Development and Ecological Restoration of Mineral Resources,China。
文摘The surrounding rock is prone to large-scale loosening and failure after the excavation of shallow large-span caverns because of the thin overlying strata and large cross-section span.The rational design of bolt support is very important to the safety control of surrounding rock as a common support means.The control mechanism and design method of bolt support for shallow-buried large-span caverns is carried out.The calculation method of bolt prestress and length based on arched failure and collapsed failure mode is established.The influence mechanism of different influencing factors on the bolt prestress and length is clarified.At the same time,the constant resistance energy-absorbing bolt with high strength and high toughness is developed,and the comparative test of mechanical properties is carried out.On this basis,the design method of high prestressed bolt support for shallow-buried large-span caverns is put forward,and the field test is carried out in Qingdao metro station in China.The monitoring results show that the maximum roof settlement is 6.8 mm after the new design method is adopted,and the effective control of the shallow-buried large-span caverns is realized.The research results can provide theoretical and technical support for the safety control of shallow-buried large-span caverns.
基金National Natural Science Foundation of China(No.22275150)。
文摘To study the influence of silicon(Si)on 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane(CL-20),NC/CL-20 composite explosives and Si/NC/CL-20 composite explosives were prepared by the electrostatic spraying method.The morphology,structure and thermal decomposition properties of the samples were analyzed using scanning electron microscopy(SEM),X-ray energy spectroscopy(EDS),infrared spectroscopy(FT-IR),and simultaneous thermal analyzer(TG-DSC).Additionally,the combustion process of the samples was tested using a high-speed camera.The results show that the addition of nano-Si contributes to the formation of composite explosives with regular morphology and smaller particle size.The Si/NC/CL-20 composite explosive has better and more uniform sphericity,with an average particle size of 73.4 nm,compared to the NC/CL-20 composite explosive.The Si/NC/CL-20 composite explosive which produced by the electrostatic spraying method,achieves physically uniform distribution of the components including NC,CL-20,Si.The addition of Si promotes the thermal decomposition of CL-20.In comparison to the NC/CL-20 composite explosive,the activation energy of the Si/NC/CL-20 composite explosive decreases by 16.78 kJ/mol,and the self-accelerated decomposition temperature and the critical temperature of thermal explosion decreases by 3.12 K and 2.61 K,respectively.Furthermore,Si/NC/CL-20 composite explosive has shorter ignition delay time and faster combustion rate compared to the NC/CL-20 composite explosive,which shows that Si can improve the combustion performance of CL-20.
