The application of nitrogen fertilizers in agricultural fields can lead to the release of nitrogen-containing gases(NCGs),such as NO_(x),NH_(3) and N_(2)O,which can significantly impact regional atmospheric environmen...The application of nitrogen fertilizers in agricultural fields can lead to the release of nitrogen-containing gases(NCGs),such as NO_(x),NH_(3) and N_(2)O,which can significantly impact regional atmospheric environment and con-tribute to global climate change.However,there remain considerable research gaps in the accurate measurement of NCGs emissions from agricultural fields,hindering the development of effective emission reduction strategies.We improved an open-top dynamic chambers(OTDCs)system and evaluated the performance by comparing the measured and given fluxes of the NCGs.The results showed that the measured fluxes of NO,N_(2)O and NH_(3)were 1%,2%and 7%lower than the given fluxes,respectively.For the determination of NH_(3) concentration,we employed a stripping coil-ion chromatograph(SC-IC)analytical technique,which demonstrated an absorption efficiency for atmospheric NH_(3) exceeding 96.1%across sampling durations of 6 to 60 min.In the summer maize season,we utilized the OTDCs system to measure the exchange fluxes of NO,NH_(3),and N_(2)O from the soil in the North China Plain.Substantial emissions of NO,NH_(3) and N_(2)O were recorded following fertilization,with peaks of 107,309,1239 ng N/(m^(2)·s),respectively.Notably,significant NCGs emissions were observed following sus-tained heavy rainfall one month after fertilization,particularly with NH_(3) peak being 4.5 times higher than that observed immediately after fertilization.Our results demonstrate that the OTDCs system accurately reflects the emission characteristics of soil NCGs and meets the requirements for long-term and continuous flux observation.展开更多
The feasibility of using a problem-dependent method to solve systems of second order ODEs is corroborated by an eigen-based theory and a methodology to develop such a numerical method is constructed.The key steps of t...The feasibility of using a problem-dependent method to solve systems of second order ODEs is corroborated by an eigen-based theory and a methodology to develop such a numerical method is constructed.The key steps of this methodology are to decouple a system of ODEs of second order into a set of uncoupled ODEs of second order;next,an eigen-dependent method is proposed to approximate the solution of each uncoupled ODE of second order.It is vital to transform all eigen-dependent methods to a problem-dependent method to bypass an Eigen analysis.The development of an eigen-dependent method plays a key role in this methodology so that slow eigenmodes can be accurately integrated while there is no instability or excessive amplitude growth in fast eigenmodes.This can explain why a problem-dependent method can simultaneously combine the explicitness of each step and A-stability.Consequently,huge computational efforts can be saved for solving nonlinear stiff problems.A new family of problem-dependent methods is developed in this work so that the feasibility of the proposed methodology can be affirmed.It has almost the same performance as that of the HHT-αmethod.However,it can save more than 99.5%of CPU demand in approximating a solution for a system of 1000 nonlinear second order ODEs.展开更多
The complex vibration directly affects the dynamic safety of drill string in ultra-deep wells and extra-deep wells.It is important to understand the dynamic characteristics of drill string to ensure the safety of dril...The complex vibration directly affects the dynamic safety of drill string in ultra-deep wells and extra-deep wells.It is important to understand the dynamic characteristics of drill string to ensure the safety of drill string.Due to the super slenderness ratio of drill string,strong nonlinearity implied in dynamic analysis and the complex load environment,dynamic simulation of drill string faces great challenges.At present,many simulation methods have been developed to analyze drill string dynamics,and node iteration method is one of them.The node iteration method has a unique advantage in dealing with the contact characteristics between drill string and borehole wall,but its drawback is that the calculation consumes a considerable amount of time.This paper presents a dynamic simulation method of drilling string in extra-deep well based on successive over-relaxation node iterative method(SOR node iteration method).Through theoretical analysis and numerical examples,the correctness and validity of this method were verified,and the dynamics characteristics of drill string in extra-deep wells were calculated and analyzed.The results demonstrate that,in contrast to the conventional node iteration method,the SOR node iteration method can increase the computational efficiency by 48.2%while achieving comparable results.And the whirl trajectory of the extra-deep well drill string is extremely complicated,the maximum rotational speed downhole is approximately twice the rotational speed on the ground.The dynamic torque increases rapidly at the position of the bottom stabilizer,and the lateral vibration in the middle and lower parts of drill string is relatively intense.展开更多
The high-speed winding spindle employs a flexible support system incorporating rubber O-rings.By precisely configuring the structural parameters and the number of the O-rings,the spindle can stably surpass its critica...The high-speed winding spindle employs a flexible support system incorporating rubber O-rings.By precisely configuring the structural parameters and the number of the O-rings,the spindle can stably surpass its critical speed points and maintain operational stability across the entire working speed range.However,the support stiffness and damping of rubber O-rings exhibit significant nonlinear frequency dependence.Conventional experimental methods for deriving equivalent stiffness and damping,based on the principle of the forced non-resonance method,require fabricating custom setups for each O-ring specification and conducting vibration tests at varying frequencies,resulting in low efficiency and high costs.This study proposes a hybrid simulation-experimental method for dynamic parameter identification.Firstly,the frequency-dependent dynamic parameters of a specific O-ring support system are experimentally obtained.Subsequently,a corresponding parametric finite element model is established to simulate and solve the equivalent elastic modulus and equivalent stiffness-damping coefficient of this O-ring support system.Ultimately,after iterative simulation,the simulated and experimental results achieve a 99.7%agreement.The parametric finite element model developed herein can directly simulate and inversely estimate frequency-dependent dynamic parameters for O-rings of different specifications but identical elastic modulus.展开更多
This paper proposes a non-intrusive computational method for mechanical dynamic systems involving a large-scale of interval uncertain parameters,aiming to reduce the computational costs and improve accuracy in determi...This paper proposes a non-intrusive computational method for mechanical dynamic systems involving a large-scale of interval uncertain parameters,aiming to reduce the computational costs and improve accuracy in determining bounds of system response.The screening method is firstly used to reduce the scale of active uncertain parameters.The sequential high-order polynomials surrogate models are then used to approximate the dynamic system’s response at each time step.To reduce the sampling cost of constructing surrogate model,the interaction effect among uncertain parameters is gradually added to the surrogate model by sequentially incorporating samples from a candidate set,which is composed of vertices and inner grid points.Finally,the points that may produce the bounds of the system response at each time step are searched using the surrogate models.The optimization algorithm is used to locate extreme points,which contribute to determining the inner points producing system response bounds.Additionally,all vertices are also checked using the surrogate models.A vehicle nonlinear dynamic model with 72 uncertain parameters is presented to demonstrate the accuracy and efficiency of the proposed uncertain computational method.展开更多
In recent years,scholars around the world have shown increasing interest in elastic support structures,leading to significant progress in dynamic modeling techniques for pipeline systems.Although multiple analytical a...In recent years,scholars around the world have shown increasing interest in elastic support structures,leading to significant progress in dynamic modeling techniques for pipeline systems.Although multiple analytical approaches exist,engineers increasingly prioritize computationally efficient,precise low-order models for practical implementation.In order to address this need,this study develops an innovative nonlinear dynamic formulation for pipelines accounting for both foundation and boundary nonlinearities.The proposed solution methodology initiates with global mode extraction using the global mode technique,followed by a detailed implementation procedure.Model validation is conducted through a cantilever pipeline case study featuring nonlinear support conditions,where strong agreement between the proposed model's predictions and finiteelement benchmark solutions demonstrates its reliability.Subsequently,a comprehensive parametric study investigates the combined effects of foundation stiffness,boundary constraints,excitation intensity,and nonlinear interaction terms on the vibrational response of the cantilever pipe.This systematic approach yields critical insights for practical engineering designs and applications.展开更多
In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step ...In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step for the dynamic response simulation of rock mass in a high in-situ stress field.In this paper,stress initialization methods,including their principles and operating procedures for reproducing steady in-situ stress state in LS-DYNA,are first introduced.Then the most popular four methods,i.e.,explicit dynamic relaxation(DR)method,implicit-explicit sequence method,Dynain file method and quasi-static method,are exemplified through a case analysis by using the RHT and plastic hardening rock material models to simulate rock blasting under in-situ stress condition.Based on the simulations,it is concluded that the stress initialization results obtained by implicit-explicit sequence method and dynain file method are closely related to the rock material model,and the explicit DR method has an obvious advantage in solution time when compared to other methods.Besides that,it is recommended to adopt two separate analyses for the whole numerical simulation of rock mass under the combined action of in-situ stress and dynamic disturbance.展开更多
An advanced Actuator Surface Method(ASM)coupled with Computational Fluid Dynamics(CFD)is developed and applied to the complex unsteady aerodynamic simulation of helicopter.By introducing an improved three-dimensional ...An advanced Actuator Surface Method(ASM)coupled with Computational Fluid Dynamics(CFD)is developed and applied to the complex unsteady aerodynamic simulation of helicopter.By introducing an improved three-dimensional anisotropic Gaussian kernel,this method effectively addresses the severe aerodynamic load fluctuations commonly associated with traditional Virtual Blade Method(VBM)due to turbulent flow around blade elements.To manage the issues of regional shape and grid cell quantity variations caused by virtual blade sweeping,a universal hybrid grid generation strategy is established without body-fitted and disk interpolation grids,which enhances the computational stability at both blade elements and blade edges.Aerodynamic numerical simulations of helicopter are performed using this method,focusing on rotor/fuselage interaction dominated by rotor wake motion and fuselage blockage effects,Blade-Vortex Interaction(BVI)induced by tip vortices,and maneuvering flights involving collective pitch ramp increases.The results indicate that the advanced ASM demonstrates reliability and robustness in the simulation of complex unsteady flow fields around helicopter.Under similar computational accuracy,the advanced ASM improves computational efficiency by nearly 40 times compared to the oversetgrid-based full Blade-Resolved(B-R)method,and by 6 times compared to the VBM.It shows significant advantages when applied to complex full-aircraft interaction and maneuvering flight conditions that require substantial computational resources.展开更多
In the leakage detection of reservoir dam bodies,traditional geophysical methods can only achieve one-time detection.Meanwhile,due to the non-uniqueness of geophysicalin version,how to improve the fidelity of geophysi...In the leakage detection of reservoir dam bodies,traditional geophysical methods can only achieve one-time detection.Meanwhile,due to the non-uniqueness of geophysicalin version,how to improve the fidelity of geophysical detection inversion profiles has become a key challenge in the industry.This study aims to construct a long-term real-time monitoring system using the time-lapse resistivity method,reveal the spatiotemporal evolution law of dam leakage,and provide technical support for accurate treatment.By integrating the Internet of Things(IoT),5G technology,and AI technology,real-time data acquisition,real-time transmission,and automatic inversionare realized.Through dynamic imaging analysis of the electrical anomaly characteristics of the leakage area and comparison between corresponding rainfall events,the leakage range and resistivity changes,the reliability and efficiency of dam leakagedete ction are significantly improved.This achieves long-termdynamic monitoring of dam leakage and provides a new perspective for the safe operation and maintenance of reservoirs.展开更多
Steep bedding slopes are widely distributed in Southwestern China’s mountainous regions and have complex seismic responses and instability risks,causing casualties and property losses.Considering the high-seismic-int...Steep bedding slopes are widely distributed in Southwestern China’s mountainous regions and have complex seismic responses and instability risks,causing casualties and property losses.Considering the high-seismic-intensity environment,the dynamic failure evolution and instability mechanism of high-steep bedding slopes are simulated via the discrete element method and shaking table test.The dynamic response characteristics and cumulative failure effects of slopes subjected to continuous ground motion are investigated.The results show that the dynamic response characteristics of slopes under continuous earthquakes are influenced by geological and topographic conditions.Elevation has a distinct impact on both the slope interior and surface,with amplification effects more pronounced on the surface.The weak interlayers have different influences on the dynamic amplification effect of slopes.Weak interlayers have dynamic magnification effects on the slope surface at relative elevations of 0-0.33 and 0.82-1.0 but have weakening effects between 0.33 and 0.82.Moreover,the weak interlayers also have controlling effects on the dynamic instability mode of slopes.The characteristics of intergranular contact failure,fracture propagation,and displacement distribution are analyzed to reveal the dynamic failure evolution and instability mechanism through the discrete-element model.The dynamic instability process of slopes includes three stages:fracture initiation(0-0.2g),fracture expansion(0.2g-0.3g),and sliding instability(0.3g-0.6g).This work can provide a valuable reference for the seismic stability and reinforcement of complex slopes.展开更多
This paper presents a framework for constructing surrogate models for sensitivity analysis of structural dynamics behavior.Physical models involving deformation,such as collisions,vibrations,and penetration,are devel-...This paper presents a framework for constructing surrogate models for sensitivity analysis of structural dynamics behavior.Physical models involving deformation,such as collisions,vibrations,and penetration,are devel-oped using the material point method.To reduce the computational cost of Monte Carlo simulations,response surface models are created as surrogate models for the material point system to approximate its dynamic behavior.An adaptive randomized greedy algorithm is employed to construct a sparse polynomial chaos expansion model with a fixed order,effectively balancing the accuracy and computational efficiency of the surrogate model.Based on the sparse polynomial chaos expansion,sensitivity analysis is conducted using the global finite difference and Sobol methods.Several examples of structural dynamics are provided to demonstrate the effectiveness of the proposed method in addressing structural dynamics problems.展开更多
0 INTRODUCTION In recent years,modern railways have been actively under construction in the complex mountainous area of Southwest China.However,rockfall poses a significant threat to both construction and operation ph...0 INTRODUCTION In recent years,modern railways have been actively under construction in the complex mountainous area of Southwest China.However,rockfall poses a significant threat to both construction and operation phases of railway projects(Yan et al.,2023;Chen et al.,2022;Fanos and Pradhan,2018).展开更多
There has been a growing trend in the development of offshore deep-water ports in China. For such deep sea projects, all-vertical-piled wharves are suitable structures and generally located in open waters, greatly aff...There has been a growing trend in the development of offshore deep-water ports in China. For such deep sea projects, all-vertical-piled wharves are suitable structures and generally located in open waters, greatly affected by wave action. Currently, no systematic studies or simplified numerical methods are available for deriving the dynamic characteristics and dynamic responses of all-vertical-piled wharves under wave cyclic loads. In this article, we compare the dynamic characteristics of an all-vertical-piled wharf with those of a traditional inshore high-piled wharf through numerical analysis; our research reveals that the vibration period of an all-vertical-piled wharf under cyclic loading is longer than that of an inshore high-piled wharf and is much closer to the period of the loading wave. Therefore, dynamic calculation and analysis should be conducted when designing and calculating the characteristics of an all-vertical-piled wharf. We establish a dynamic finite element model to examine the dynamic response of an all-vertical-piled wharf under wave cyclic loads and compare the results with those under wave equivalent static load; the comparison indicates that dynamic amplification of the structure is evident when the wave dynamic load effect is taken into account. Furthermore, a simplified dynamic numerical method for calculating the dynamic response of an all-vertical-piled wharf is established based on the P-Y curve. Compared with finite element analysis, the simplified method is more convenient to use and applicable to large structural deformation while considering the soil non-linearity. We confirmed that the simplified method has acceptable accuracy and can be used in engineering applications.展开更多
The relationship between capillary pressure and saturation plays a critical role in the characterization of two-phase flow and transport in aquifers and oil reservoirs. This relationship is usually determined under th...The relationship between capillary pressure and saturation plays a critical role in the characterization of two-phase flow and transport in aquifers and oil reservoirs. This relationship is usually determined under the static condition, where capillary pressure is the only function of saturation. However,considerable experiments have suggested that the dependence of capillary pressure on desaturation rate is under the dynamic condition. Thus, a more general description of capillary pressure that includes dynamic capillary effect has been approved widely. A comprehensive understanding of the dynamic capillary effect is needed for the investigation of the two-phase flow in porous media by various methods. In general, dynamic capillary effect in porous media can be studied through the laboratory experiment, pore-to macro-scale modeling, and artificial neural network. Here, main principle and research procedures of each method are reviewed in detail. Then, research progress, disadvantages and advantages are discussed, respectively. In addition, upscaling study from pore-to macro-scale are introduced, which explains the difference between laboratory experiment and pore-scale modeling. At last, several future perspectives and recommendations for optimal solution of dynamic capillary effect are presented.展开更多
The sea surface temperature （SST） in the In- dian Ocean affects the regional climate over the Asian continent mostly through a modulation of the monsoon system. It is still difficult to provide an a priori indicatio...The sea surface temperature （SST） in the In- dian Ocean affects the regional climate over the Asian continent mostly through a modulation of the monsoon system. It is still difficult to provide an a priori indication of the seasonal variability over the Indian Ocean. It is widely recognized that the warm and cold events of SST over the tropical Indian Ocean are strongly linked to those of the equatorial eastern Pacific. In this study, a statistical prediction model has been developed to predict the monthly SST over the tropical Indian Ocean. This model is a linear regression model based on the lag relationship between the SST over the tropical Indian Ocean and the Nino3.4 （5°S-5°N, 170°W-120°W） SST Index. The pre- dictor （i.e., Nino3.4 SST Index） has been operationally predicted by a large size ensemble E1 Nifio and the Southern Oscillation （ENSO） forecast system with cou- pled data assimilation （Leefs_CDA）, which achieves a high predictive skill of up to a 24-month lead time for the equatorial eastern Pacific SST. As a result, the prediction skill of the present statistical model over the tropical In- dian Ocean is better than that of persistence prediction for January 1982 through December 2009.展开更多
Synchronization is one of the most important characteristics of dynamic systems.For this paper,the authors obtained results for the nonlinear systems controller for the custom Synchronization of two 4D systems.The fin...Synchronization is one of the most important characteristics of dynamic systems.For this paper,the authors obtained results for the nonlinear systems controller for the custom Synchronization of two 4D systems.The findings have allowed authors to develop two analytical approaches using the second Lyapunov(Lyp)method and the Gardanomethod.Since the Gardano method does not involve the development of special positive Lyp functions,it is very efficient and convenient to achieve excessive systemSYCR phenomena.Error is overcome by using Gardano and overcoming some problems in Lyp.Thus we get a great investigation into the convergence of error dynamics,the authors in this paper are interested in giving numerical simulations of the proposed model to clarify the results and check them,an important aspect that will be studied is Synchronization Complete hybrid SYCR and anti-Synchronization,by making use of the Lyapunov expansion analysis,a proposed control method is developed to determine the actual.The basic idea in the proposed way is to receive the evolution of between two methods.Finally,the present model has been applied and showing in a new attractor,and the obtained results are compared with other approximate results,and the nearly good coincidence was obtained.展开更多
This paper presents Modified Chebyshev-Picard Iteration(MCPI)methods for long-term integration of the coupled orbit and attitude dynamics.Although most orbit predictions for operational satellites have assumed that th...This paper presents Modified Chebyshev-Picard Iteration(MCPI)methods for long-term integration of the coupled orbit and attitude dynamics.Although most orbit predictions for operational satellites have assumed that the attitude dynamics is decoupled from the orbit dynamics,the fully coupled dynamics is required for the solutions of uncontrolled space debris and space objects with high area-to-mass ratio,for which cross sectional area is constantly changing leading to significant change on the solar radiation pressure and atmospheric drag.MCPI is a set of methods for solution of initial value problems and boundary value problems.The methods refine an orthogonal function approximation of long-time-interval segments of state trajectories iteratively by fusing Chebyshev polynomials with the classical Picard iteration and have been applied to multiple challenging aerospace problems.Through the studies on integrating a torque-free rigid body rotation and a long-term integration of the coupled orbit-attitude dynamics through the effect of solar radiation pressure,MCPI methods are shown to achieve several times speedup over the Runge-Kutta 7(8)methods with several orders of magnitudes of better accuracy.MCPI methods are further optimized by integrating the decoupled dynamics at the beginning of the iteration and coupling the full dynamics when the attitude solutions and orbit solutions are converging during the iteration.The approach of decoupling and then coupling during iterations provides a unique and promising perspective on the way to warm start the solution process for the longterm integration of the coupled orbit-attitude dynamics.Furthermore,an attractive feature of MCPI in maintaining the unity constraint for the integration of quaternions within machine accuracy is illustrated to be very appealing.展开更多
The reduced basis methods (RBM) have been demonstrated as a promising numerical technique for statics problems and are extended to structural dynamic problems in this paper. Direct step-by-step integration and mode su...The reduced basis methods (RBM) have been demonstrated as a promising numerical technique for statics problems and are extended to structural dynamic problems in this paper. Direct step-by-step integration and mode superposition are the most widely used methods in the field of the finite element analysis of structural dynamic response and solid mechanics. Herein these two methods are both transformed into reduced forms according to the proposed reduced basis methods. To generate a reduced surrogate model with small size, a greedy algorithm is suggested to construct sample set and reduced basis space adaptively in a prescribed training parameter space. For mode superposition method, the reduced basis space comprises the truncated eigenvectors from generalized eigenvalue problem associated with selected sample parameters. The reduced generalized eigenvalue problem is obtained by the projection of original generalized eigenvalue problem onto the reduced basis space. In the situation of direct integration, the solutions of the original increment formulation corresponding to the sample set are extracted to construct the reduced basis space. The reduced increment formulation is formed by the same method as mode superposition method. Numerical example is given in Section 5 to validate the efficiency of the presented reduced basis methods for structural dynamic problems.展开更多
This paper estimated the liquefaction potential of a saturated soil deposit subjected to a horizontal seismic excitation at its base using the total stress approach.A comparative analysis between the simplified and th...This paper estimated the liquefaction potential of a saturated soil deposit subjected to a horizontal seismic excitation at its base using the total stress approach.A comparative analysis between the simplified and the nonlinear dynamic methods was used to verify to what extent the simplified method could be reliable.In order to generalise the reliability of the simplified method for any value of the maximum acceleration for the used earthquakes,a correction for the maximum acceleration less than 0.3 g was proposed based on the comparison of safety factor values determined by the dynamic method illustrated by the equivalent linear model with lumped masses and the simplified method for a given profile of soil subjected to 38 earthquakes.The nonlinear behaviour of soil was represented by two hyperbolic models:Hardin and Drnevich,and Masing.To determine the cyclic resistance ratio(CRR),the cone penetration test(CPT) based method,the standard penetration test(SPT) based method,and the shear wave velocity based method were used.The safety factor was calculated as the ratio of CRR/CSR,where CSR represents the cyclic stress ratio.The results of the proposed correction have given smaller values of the safety factor compared to the nonlinear dynamic methods for the maximum acceleration less than 0.3 g.In other words,by considering this correction,the most unfavourable case is always given by the modified simplified method.展开更多
An efficient active-set approach is presented for both nonnegative and general linear programming by adding varying numbers of constraints at each iteration. Computational experiments demonstrate that the proposed app...An efficient active-set approach is presented for both nonnegative and general linear programming by adding varying numbers of constraints at each iteration. Computational experiments demonstrate that the proposed approach is significantly faster than previous active-set and standard linear programming algorithms.展开更多
基金supported by the National Key Research and Develop-ment Program(No.2022YFC3701103)the National Natural Science Foundation of China(Nos.42130714 and 41931287).
文摘The application of nitrogen fertilizers in agricultural fields can lead to the release of nitrogen-containing gases(NCGs),such as NO_(x),NH_(3) and N_(2)O,which can significantly impact regional atmospheric environment and con-tribute to global climate change.However,there remain considerable research gaps in the accurate measurement of NCGs emissions from agricultural fields,hindering the development of effective emission reduction strategies.We improved an open-top dynamic chambers(OTDCs)system and evaluated the performance by comparing the measured and given fluxes of the NCGs.The results showed that the measured fluxes of NO,N_(2)O and NH_(3)were 1%,2%and 7%lower than the given fluxes,respectively.For the determination of NH_(3) concentration,we employed a stripping coil-ion chromatograph(SC-IC)analytical technique,which demonstrated an absorption efficiency for atmospheric NH_(3) exceeding 96.1%across sampling durations of 6 to 60 min.In the summer maize season,we utilized the OTDCs system to measure the exchange fluxes of NO,NH_(3),and N_(2)O from the soil in the North China Plain.Substantial emissions of NO,NH_(3) and N_(2)O were recorded following fertilization,with peaks of 107,309,1239 ng N/(m^(2)·s),respectively.Notably,significant NCGs emissions were observed following sus-tained heavy rainfall one month after fertilization,particularly with NH_(3) peak being 4.5 times higher than that observed immediately after fertilization.Our results demonstrate that the OTDCs system accurately reflects the emission characteristics of soil NCGs and meets the requirements for long-term and continuous flux observation.
文摘The feasibility of using a problem-dependent method to solve systems of second order ODEs is corroborated by an eigen-based theory and a methodology to develop such a numerical method is constructed.The key steps of this methodology are to decouple a system of ODEs of second order into a set of uncoupled ODEs of second order;next,an eigen-dependent method is proposed to approximate the solution of each uncoupled ODE of second order.It is vital to transform all eigen-dependent methods to a problem-dependent method to bypass an Eigen analysis.The development of an eigen-dependent method plays a key role in this methodology so that slow eigenmodes can be accurately integrated while there is no instability or excessive amplitude growth in fast eigenmodes.This can explain why a problem-dependent method can simultaneously combine the explicitness of each step and A-stability.Consequently,huge computational efforts can be saved for solving nonlinear stiff problems.A new family of problem-dependent methods is developed in this work so that the feasibility of the proposed methodology can be affirmed.It has almost the same performance as that of the HHT-αmethod.However,it can save more than 99.5%of CPU demand in approximating a solution for a system of 1000 nonlinear second order ODEs.
基金supported by the National Natural Science Foundation of China(52174003,52374008).
文摘The complex vibration directly affects the dynamic safety of drill string in ultra-deep wells and extra-deep wells.It is important to understand the dynamic characteristics of drill string to ensure the safety of drill string.Due to the super slenderness ratio of drill string,strong nonlinearity implied in dynamic analysis and the complex load environment,dynamic simulation of drill string faces great challenges.At present,many simulation methods have been developed to analyze drill string dynamics,and node iteration method is one of them.The node iteration method has a unique advantage in dealing with the contact characteristics between drill string and borehole wall,but its drawback is that the calculation consumes a considerable amount of time.This paper presents a dynamic simulation method of drilling string in extra-deep well based on successive over-relaxation node iterative method(SOR node iteration method).Through theoretical analysis and numerical examples,the correctness and validity of this method were verified,and the dynamics characteristics of drill string in extra-deep wells were calculated and analyzed.The results demonstrate that,in contrast to the conventional node iteration method,the SOR node iteration method can increase the computational efficiency by 48.2%while achieving comparable results.And the whirl trajectory of the extra-deep well drill string is extremely complicated,the maximum rotational speed downhole is approximately twice the rotational speed on the ground.The dynamic torque increases rapidly at the position of the bottom stabilizer,and the lateral vibration in the middle and lower parts of drill string is relatively intense.
基金National Key R&D Program of China(No.2017YFB1304000)Fundamental Research Funds for the Central Universities,China(No.2232023G-05-1)。
文摘The high-speed winding spindle employs a flexible support system incorporating rubber O-rings.By precisely configuring the structural parameters and the number of the O-rings,the spindle can stably surpass its critical speed points and maintain operational stability across the entire working speed range.However,the support stiffness and damping of rubber O-rings exhibit significant nonlinear frequency dependence.Conventional experimental methods for deriving equivalent stiffness and damping,based on the principle of the forced non-resonance method,require fabricating custom setups for each O-ring specification and conducting vibration tests at varying frequencies,resulting in low efficiency and high costs.This study proposes a hybrid simulation-experimental method for dynamic parameter identification.Firstly,the frequency-dependent dynamic parameters of a specific O-ring support system are experimentally obtained.Subsequently,a corresponding parametric finite element model is established to simulate and solve the equivalent elastic modulus and equivalent stiffness-damping coefficient of this O-ring support system.Ultimately,after iterative simulation,the simulated and experimental results achieve a 99.7%agreement.The parametric finite element model developed herein can directly simulate and inversely estimate frequency-dependent dynamic parameters for O-rings of different specifications but identical elastic modulus.
基金supported by the National Natural Science Foundation of China(Grant No.12272142)Fundamental Research Funds for the Central Universities(Grant No.2172021XXJS048)。
文摘This paper proposes a non-intrusive computational method for mechanical dynamic systems involving a large-scale of interval uncertain parameters,aiming to reduce the computational costs and improve accuracy in determining bounds of system response.The screening method is firstly used to reduce the scale of active uncertain parameters.The sequential high-order polynomials surrogate models are then used to approximate the dynamic system’s response at each time step.To reduce the sampling cost of constructing surrogate model,the interaction effect among uncertain parameters is gradually added to the surrogate model by sequentially incorporating samples from a candidate set,which is composed of vertices and inner grid points.Finally,the points that may produce the bounds of the system response at each time step are searched using the surrogate models.The optimization algorithm is used to locate extreme points,which contribute to determining the inner points producing system response bounds.Additionally,all vertices are also checked using the surrogate models.A vehicle nonlinear dynamic model with 72 uncertain parameters is presented to demonstrate the accuracy and efficiency of the proposed uncertain computational method.
基金supported by the National Natural Science Foundation of China(Nos.52401342 and 12572025)the Fundamental Research Funds for the Central Universities of China(Nos.D5000240076 and G2025KY05171)+1 种基金the Natural Science Basic Research Program of Shaanxi Province(No.2025JCYBMS-026)the Basic Research Programs of Taicang(No.TC2024JC36)。
文摘In recent years,scholars around the world have shown increasing interest in elastic support structures,leading to significant progress in dynamic modeling techniques for pipeline systems.Although multiple analytical approaches exist,engineers increasingly prioritize computationally efficient,precise low-order models for practical implementation.In order to address this need,this study develops an innovative nonlinear dynamic formulation for pipelines accounting for both foundation and boundary nonlinearities.The proposed solution methodology initiates with global mode extraction using the global mode technique,followed by a detailed implementation procedure.Model validation is conducted through a cantilever pipeline case study featuring nonlinear support conditions,where strong agreement between the proposed model's predictions and finiteelement benchmark solutions demonstrates its reliability.Subsequently,a comprehensive parametric study investigates the combined effects of foundation stiffness,boundary constraints,excitation intensity,and nonlinear interaction terms on the vibrational response of the cantilever pipe.This systematic approach yields critical insights for practical engineering designs and applications.
基金Project(41630642)supported by the Key Project of National Natural Science Foundation of ChinaProject(51974360)supported by the National Natural Science Foundation of ChinaProject(2018JJ3656)supported by the Natural Science Foundation of Hunan Province,China。
文摘In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step for the dynamic response simulation of rock mass in a high in-situ stress field.In this paper,stress initialization methods,including their principles and operating procedures for reproducing steady in-situ stress state in LS-DYNA,are first introduced.Then the most popular four methods,i.e.,explicit dynamic relaxation(DR)method,implicit-explicit sequence method,Dynain file method and quasi-static method,are exemplified through a case analysis by using the RHT and plastic hardening rock material models to simulate rock blasting under in-situ stress condition.Based on the simulations,it is concluded that the stress initialization results obtained by implicit-explicit sequence method and dynain file method are closely related to the rock material model,and the explicit DR method has an obvious advantage in solution time when compared to other methods.Besides that,it is recommended to adopt two separate analyses for the whole numerical simulation of rock mass under the combined action of in-situ stress and dynamic disturbance.
基金co-supported by the Foundation of the State Key Laboratory of Aerodynamics(No.RAL202203)the National Key Laboratory Foundation(No.6142202202)the China Postdoctoral Science Foundation(No.2024M754133)。
文摘An advanced Actuator Surface Method(ASM)coupled with Computational Fluid Dynamics(CFD)is developed and applied to the complex unsteady aerodynamic simulation of helicopter.By introducing an improved three-dimensional anisotropic Gaussian kernel,this method effectively addresses the severe aerodynamic load fluctuations commonly associated with traditional Virtual Blade Method(VBM)due to turbulent flow around blade elements.To manage the issues of regional shape and grid cell quantity variations caused by virtual blade sweeping,a universal hybrid grid generation strategy is established without body-fitted and disk interpolation grids,which enhances the computational stability at both blade elements and blade edges.Aerodynamic numerical simulations of helicopter are performed using this method,focusing on rotor/fuselage interaction dominated by rotor wake motion and fuselage blockage effects,Blade-Vortex Interaction(BVI)induced by tip vortices,and maneuvering flights involving collective pitch ramp increases.The results indicate that the advanced ASM demonstrates reliability and robustness in the simulation of complex unsteady flow fields around helicopter.Under similar computational accuracy,the advanced ASM improves computational efficiency by nearly 40 times compared to the oversetgrid-based full Blade-Resolved(B-R)method,and by 6 times compared to the VBM.It shows significant advantages when applied to complex full-aircraft interaction and maneuvering flight conditions that require substantial computational resources.
文摘In the leakage detection of reservoir dam bodies,traditional geophysical methods can only achieve one-time detection.Meanwhile,due to the non-uniqueness of geophysicalin version,how to improve the fidelity of geophysical detection inversion profiles has become a key challenge in the industry.This study aims to construct a long-term real-time monitoring system using the time-lapse resistivity method,reveal the spatiotemporal evolution law of dam leakage,and provide technical support for accurate treatment.By integrating the Internet of Things(IoT),5G technology,and AI technology,real-time data acquisition,real-time transmission,and automatic inversionare realized.Through dynamic imaging analysis of the electrical anomaly characteristics of the leakage area and comparison between corresponding rainfall events,the leakage range and resistivity changes,the reliability and efficiency of dam leakagedete ction are significantly improved.This achieves long-termdynamic monitoring of dam leakage and provides a new perspective for the safe operation and maintenance of reservoirs.
基金Project(52108361)supported by the National Natural Science Foundation of ChinaProjects(BK20231217,BK20220265)supported by the Basic Research Program of Jiangsu Province,China+5 种基金Project(sklhse-KF-2025-D-02)supported by the Open Research Fund Program of the State Key Laboratory of Hydroscience and Engineering,ChinaProject(2023ZB15)supported by the Independent Research Project of the State Key Laboratory of Subtropical Building and Urban Science,ChinaProject(SKLGME023001)supported by the Key Laboratory of Geomechanics and Geotechnical Engineering Safety,the Chinese Academy of SciencesProject(2025A04J3992)supported by the Basic and Applied Basic Research Project of the Guangzhou Science and Technology Bureau,ChinaProject(SKLGP2022Z015)supported by the State Key Laboratory of Geohazard Prevention and Geoenvironment Protection Independent Research Project,ChinaProjects(2023YFS0436,2024NSFSC1715)supported by the Science and Technology Department of Sichuan Province,China。
文摘Steep bedding slopes are widely distributed in Southwestern China’s mountainous regions and have complex seismic responses and instability risks,causing casualties and property losses.Considering the high-seismic-intensity environment,the dynamic failure evolution and instability mechanism of high-steep bedding slopes are simulated via the discrete element method and shaking table test.The dynamic response characteristics and cumulative failure effects of slopes subjected to continuous ground motion are investigated.The results show that the dynamic response characteristics of slopes under continuous earthquakes are influenced by geological and topographic conditions.Elevation has a distinct impact on both the slope interior and surface,with amplification effects more pronounced on the surface.The weak interlayers have different influences on the dynamic amplification effect of slopes.Weak interlayers have dynamic magnification effects on the slope surface at relative elevations of 0-0.33 and 0.82-1.0 but have weakening effects between 0.33 and 0.82.Moreover,the weak interlayers also have controlling effects on the dynamic instability mode of slopes.The characteristics of intergranular contact failure,fracture propagation,and displacement distribution are analyzed to reveal the dynamic failure evolution and instability mechanism through the discrete-element model.The dynamic instability process of slopes includes three stages:fracture initiation(0-0.2g),fracture expansion(0.2g-0.3g),and sliding instability(0.3g-0.6g).This work can provide a valuable reference for the seismic stability and reinforcement of complex slopes.
基金support from the National Natural Science Foundation of China(Grant Nos.52174123&52274222).
文摘This paper presents a framework for constructing surrogate models for sensitivity analysis of structural dynamics behavior.Physical models involving deformation,such as collisions,vibrations,and penetration,are devel-oped using the material point method.To reduce the computational cost of Monte Carlo simulations,response surface models are created as surrogate models for the material point system to approximate its dynamic behavior.An adaptive randomized greedy algorithm is employed to construct a sparse polynomial chaos expansion model with a fixed order,effectively balancing the accuracy and computational efficiency of the surrogate model.Based on the sparse polynomial chaos expansion,sensitivity analysis is conducted using the global finite difference and Sobol methods.Several examples of structural dynamics are provided to demonstrate the effectiveness of the proposed method in addressing structural dynamics problems.
基金supported by the Open Research Fund of Key Laboratory of Geological Hazards on Three Gorges Reservoir Area(China Three Gorges University),Ministry of Education(No.2022KDZ03)the Science and Technology Projects of Yunnan Provincial Science and Technology Department(No.202401AT070328)+1 种基金the Young talents project of“Xingdian Talent Support Program”in Yunnan Province(No.YNWR-QNBJ-2020-019)the Fund Project of China Academy of Railway Sciences Co.,Ltd.(No.2021YJ178)。
文摘0 INTRODUCTION In recent years,modern railways have been actively under construction in the complex mountainous area of Southwest China.However,rockfall poses a significant threat to both construction and operation phases of railway projects(Yan et al.,2023;Chen et al.,2022;Fanos and Pradhan,2018).
基金financially supported by the Major Science and Technology Project of MOT,China(Grant Nos.2013 328 224 070 and 2014 328 224 040)the National Natural Science Foundation of China(Grant No.51409134)
文摘There has been a growing trend in the development of offshore deep-water ports in China. For such deep sea projects, all-vertical-piled wharves are suitable structures and generally located in open waters, greatly affected by wave action. Currently, no systematic studies or simplified numerical methods are available for deriving the dynamic characteristics and dynamic responses of all-vertical-piled wharves under wave cyclic loads. In this article, we compare the dynamic characteristics of an all-vertical-piled wharf with those of a traditional inshore high-piled wharf through numerical analysis; our research reveals that the vibration period of an all-vertical-piled wharf under cyclic loading is longer than that of an inshore high-piled wharf and is much closer to the period of the loading wave. Therefore, dynamic calculation and analysis should be conducted when designing and calculating the characteristics of an all-vertical-piled wharf. We establish a dynamic finite element model to examine the dynamic response of an all-vertical-piled wharf under wave cyclic loads and compare the results with those under wave equivalent static load; the comparison indicates that dynamic amplification of the structure is evident when the wave dynamic load effect is taken into account. Furthermore, a simplified dynamic numerical method for calculating the dynamic response of an all-vertical-piled wharf is established based on the P-Y curve. Compared with finite element analysis, the simplified method is more convenient to use and applicable to large structural deformation while considering the soil non-linearity. We confirmed that the simplified method has acceptable accuracy and can be used in engineering applications.
基金financially supported by the National Natural Science Foundation of China (No. 42102149)the Fundamental Research Funds for the Central Universities (No. 2462021YXZZ005)。
文摘The relationship between capillary pressure and saturation plays a critical role in the characterization of two-phase flow and transport in aquifers and oil reservoirs. This relationship is usually determined under the static condition, where capillary pressure is the only function of saturation. However,considerable experiments have suggested that the dependence of capillary pressure on desaturation rate is under the dynamic condition. Thus, a more general description of capillary pressure that includes dynamic capillary effect has been approved widely. A comprehensive understanding of the dynamic capillary effect is needed for the investigation of the two-phase flow in porous media by various methods. In general, dynamic capillary effect in porous media can be studied through the laboratory experiment, pore-to macro-scale modeling, and artificial neural network. Here, main principle and research procedures of each method are reviewed in detail. Then, research progress, disadvantages and advantages are discussed, respectively. In addition, upscaling study from pore-to macro-scale are introduced, which explains the difference between laboratory experiment and pore-scale modeling. At last, several future perspectives and recommendations for optimal solution of dynamic capillary effect are presented.
基金supported by the National Basic Research Program of China (Grant No. 2012CB417404)the National Natural Science Foundation of China (Grant Nos.41075064 and 41176014)
文摘The sea surface temperature （SST） in the In- dian Ocean affects the regional climate over the Asian continent mostly through a modulation of the monsoon system. It is still difficult to provide an a priori indication of the seasonal variability over the Indian Ocean. It is widely recognized that the warm and cold events of SST over the tropical Indian Ocean are strongly linked to those of the equatorial eastern Pacific. In this study, a statistical prediction model has been developed to predict the monthly SST over the tropical Indian Ocean. This model is a linear regression model based on the lag relationship between the SST over the tropical Indian Ocean and the Nino3.4 （5°S-5°N, 170°W-120°W） SST Index. The pre- dictor （i.e., Nino3.4 SST Index） has been operationally predicted by a large size ensemble E1 Nifio and the Southern Oscillation （ENSO） forecast system with cou- pled data assimilation （Leefs_CDA）, which achieves a high predictive skill of up to a 24-month lead time for the equatorial eastern Pacific SST. As a result, the prediction skill of the present statistical model over the tropical In- dian Ocean is better than that of persistence prediction for January 1982 through December 2009.
文摘Synchronization is one of the most important characteristics of dynamic systems.For this paper,the authors obtained results for the nonlinear systems controller for the custom Synchronization of two 4D systems.The findings have allowed authors to develop two analytical approaches using the second Lyapunov(Lyp)method and the Gardanomethod.Since the Gardano method does not involve the development of special positive Lyp functions,it is very efficient and convenient to achieve excessive systemSYCR phenomena.Error is overcome by using Gardano and overcoming some problems in Lyp.Thus we get a great investigation into the convergence of error dynamics,the authors in this paper are interested in giving numerical simulations of the proposed model to clarify the results and check them,an important aspect that will be studied is Synchronization Complete hybrid SYCR and anti-Synchronization,by making use of the Lyapunov expansion analysis,a proposed control method is developed to determine the actual.The basic idea in the proposed way is to receive the evolution of between two methods.Finally,the present model has been applied and showing in a new attractor,and the obtained results are compared with other approximate results,and the nearly good coincidence was obtained.
文摘This paper presents Modified Chebyshev-Picard Iteration(MCPI)methods for long-term integration of the coupled orbit and attitude dynamics.Although most orbit predictions for operational satellites have assumed that the attitude dynamics is decoupled from the orbit dynamics,the fully coupled dynamics is required for the solutions of uncontrolled space debris and space objects with high area-to-mass ratio,for which cross sectional area is constantly changing leading to significant change on the solar radiation pressure and atmospheric drag.MCPI is a set of methods for solution of initial value problems and boundary value problems.The methods refine an orthogonal function approximation of long-time-interval segments of state trajectories iteratively by fusing Chebyshev polynomials with the classical Picard iteration and have been applied to multiple challenging aerospace problems.Through the studies on integrating a torque-free rigid body rotation and a long-term integration of the coupled orbit-attitude dynamics through the effect of solar radiation pressure,MCPI methods are shown to achieve several times speedup over the Runge-Kutta 7(8)methods with several orders of magnitudes of better accuracy.MCPI methods are further optimized by integrating the decoupled dynamics at the beginning of the iteration and coupling the full dynamics when the attitude solutions and orbit solutions are converging during the iteration.The approach of decoupling and then coupling during iterations provides a unique and promising perspective on the way to warm start the solution process for the longterm integration of the coupled orbit-attitude dynamics.Furthermore,an attractive feature of MCPI in maintaining the unity constraint for the integration of quaternions within machine accuracy is illustrated to be very appealing.
文摘The reduced basis methods (RBM) have been demonstrated as a promising numerical technique for statics problems and are extended to structural dynamic problems in this paper. Direct step-by-step integration and mode superposition are the most widely used methods in the field of the finite element analysis of structural dynamic response and solid mechanics. Herein these two methods are both transformed into reduced forms according to the proposed reduced basis methods. To generate a reduced surrogate model with small size, a greedy algorithm is suggested to construct sample set and reduced basis space adaptively in a prescribed training parameter space. For mode superposition method, the reduced basis space comprises the truncated eigenvectors from generalized eigenvalue problem associated with selected sample parameters. The reduced generalized eigenvalue problem is obtained by the projection of original generalized eigenvalue problem onto the reduced basis space. In the situation of direct integration, the solutions of the original increment formulation corresponding to the sample set are extracted to construct the reduced basis space. The reduced increment formulation is formed by the same method as mode superposition method. Numerical example is given in Section 5 to validate the efficiency of the presented reduced basis methods for structural dynamic problems.
文摘This paper estimated the liquefaction potential of a saturated soil deposit subjected to a horizontal seismic excitation at its base using the total stress approach.A comparative analysis between the simplified and the nonlinear dynamic methods was used to verify to what extent the simplified method could be reliable.In order to generalise the reliability of the simplified method for any value of the maximum acceleration for the used earthquakes,a correction for the maximum acceleration less than 0.3 g was proposed based on the comparison of safety factor values determined by the dynamic method illustrated by the equivalent linear model with lumped masses and the simplified method for a given profile of soil subjected to 38 earthquakes.The nonlinear behaviour of soil was represented by two hyperbolic models:Hardin and Drnevich,and Masing.To determine the cyclic resistance ratio(CRR),the cone penetration test(CPT) based method,the standard penetration test(SPT) based method,and the shear wave velocity based method were used.The safety factor was calculated as the ratio of CRR/CSR,where CSR represents the cyclic stress ratio.The results of the proposed correction have given smaller values of the safety factor compared to the nonlinear dynamic methods for the maximum acceleration less than 0.3 g.In other words,by considering this correction,the most unfavourable case is always given by the modified simplified method.
文摘An efficient active-set approach is presented for both nonnegative and general linear programming by adding varying numbers of constraints at each iteration. Computational experiments demonstrate that the proposed approach is significantly faster than previous active-set and standard linear programming algorithms.
