Deployment of buoy systems is one of the most important procedures for the operation of buoy system. In the present study, a single-point mooring buoy system which contains surface buoy, cable segments with components...Deployment of buoy systems is one of the most important procedures for the operation of buoy system. In the present study, a single-point mooring buoy system which contains surface buoy, cable segments with components, anchor and so on is modeled by applying multi-body dynamics method. The motion equations are developed in discrete node description and fully Cartesian coordinates. Then numerical method is used to solve the ordinary differential equations and dynamics simulations are achieved while anchor is casting from board. The trajectories and velocities of different nodes without current and with current in buoy system are obtained. The transient tension force of each part of the cable is analyzed in the process of deployment. Numerical results indicate that the transient payload increases to a peak value when the anchor is touching the seabed and the maximum tension force will vary with different floating configuration. This work is helpful for design and deployment planning of buoy system.展开更多
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 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.展开更多
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 order to analyze the influence of configuration parameters on dynamic characteristics of machine tools in the working space, the configuration parameters have been suggested based on the orthogonal experiment metho...In order to analyze the influence of configuration parameters on dynamic characteristics of machine tools in the working space, the configuration parameters have been suggested based on the orthogonal experiment method. Dynamic analysis of a milling machine, which is newly designed for producing turbine blades, has been conducted by utilizing the modal synthesis method. The finite element model is verified and updated by experimental modal analysis （EMA） of the machine tool. The result gained by modal synthesis method is compared with whole-model finite element method （FEM） result as well. According to the orthogonal experiment method, four configuration parameters of machine tool are considered as four factors for dynamic characteristics. The influence of configuration parameters on the first three natural frequencies is obtained by range analysis. It is pointed out that configuration parameter is the most important factor affecting the fundamental frequency of machine tools, and configuration parameter has less effect on lower-order modes of the system than others. The combination of configuration parameters which makes the fundamental frequency reach the maximum value is provided. Through demonstration, the conclusion can be drawn that the influence of configuration parameters on the natural frequencies of machine tools can be analyzed explicitly by the orthogonal experiment method, which offers a new method for estimating the dynamic characteristics of machine tools.展开更多
The impact dynamics of a flexible multibody system is investigated. By using a partition method, the system is divided into two parts, the local impact region and the region away from the impact. The two parts are con...The impact dynamics of a flexible multibody system is investigated. By using a partition method, the system is divided into two parts, the local impact region and the region away from the impact. The two parts are connected by specific boundary conditions, and the system after partition is equivalent to the original system. According to the rigid-flexible coupling dynamic theory of multibody system, system＇s rigid-flexible coupling dynamic equations without impact are derived. A local impulse method for establishing the initial impact conditions is proposed. It satisfies the compatibility con- ditions for contact constraints and the actual physical situation of the impact process of flexible bodies. Based on the contact constraint method, system＇s impact dynamic equa- tions are derived in a differential-algebraic form. The contact/separation criterion and the algorithm are given. An impact dynamic simulation is given. The results show that system＇s dynamic behaviors including the energy, the deformations, the displacements, and the impact force during the impact process change dramatically. The impact makes great effects on the global dynamics of the system during and after impact.展开更多
In re-entry, the drilling riser hanging to the holding vessel takes on a free hanging state, waiting to be moved from the initial random position to the wellhead. For the re-entry, dynamics calculation is often done t...In re-entry, the drilling riser hanging to the holding vessel takes on a free hanging state, waiting to be moved from the initial random position to the wellhead. For the re-entry, dynamics calculation is often done to predict the riser motion or evaluate the structural safety. A dynamics calculation method based on Flexible Segment Model （FSM） is proposed for free hanging marine risers. In FSM, a riser is discretized into a series of flexible segments. For each flexible segment, its deflection feature and external forces are analyzed independently. For the whole riser, the nonlinear governing equations are listed according to the moment equilibrium at nodes. For the solution of the nonlinear equations, a linearization iteration scheme is provided in the paper. Owing to its flexibility, each segment can match a long part of the riser body, which enables that good results can be obtained even with a small number of segments. Moreover, the linearization iteration scheme can avoid widely used Newton-Rapson iteration scheme in which the calculation stability is influenced by the initial points. The FSM-based dynamics calculation is timesaving and stable, so suitable for the shape prediction or real-time control of free hanging marine risers.展开更多
There are vast constraint equations in conventional dynamics analysis of deployable structures,which lead to differential-algebraic equations(DAEs)solved hard.To reduce the difficulty of solving and the amount of equa...There are vast constraint equations in conventional dynamics analysis of deployable structures,which lead to differential-algebraic equations(DAEs)solved hard.To reduce the difficulty of solving and the amount of equations,a new flexible multibody dynamics analysis methodology of deployable structures with scissor-like elements(SLEs)is presented.Firstly,a precise model of a flexible bar of SLE is established by the higher order shear deformable beam element based on the absolute nodal coordinate formulation(ANCF),and the master/slave freedom method is used to obtain the dynamics equations of SLEs without constraint equations.Secondly,according to features of deployable structures,the specification matrix method(SMM)is proposed to eliminate the constraint equations among SLEs in the frame of ANCF.With this method,the inner and the boundary nodal coordinates of element characteristic matrices can be separated simply and efficiently,especially on condition that there are vast nodal coordinates.So the element characteristic matrices can be added end to end circularly.Thus,the dynamic model of deployable structure reduces dimension and can be assembled without any constraint equation.Next,a new iteration procedure for the generalized-a algorithm is presented to solve the ordinary differential equations(ODEs)of deployable structure.Finally,the proposed methodology is used to analyze the flexible multi-body dynamics of a planar linear array deployable structure based on three scissor-like elements.The simulation results show that flexibility has a significant influence on the deployment motion of the deployable structure.The proposed methodology indeed reduce the difficulty of solving and the amount of equations by eliminating redundant degrees of freedom and the constraint equations in scissor-like elements and among scissor-like elements.展开更多
A three-dimensional discrete element model of the connective type is presented. Moreover,a three-dimensional numerical analysis code,which can carry out the transitional pro- cess from connective model(for continuum)t...A three-dimensional discrete element model of the connective type is presented. Moreover,a three-dimensional numerical analysis code,which can carry out the transitional pro- cess from connective model(for continuum)to contact model(for non-continuum),is developed for simulating the mechanical process from continuum to non-continuum.The wave propagation process in a concrete block(as continuum)made of cement grout under impact loading is numer- ically simulated with this code.By comparing its numerical results with those by LS-DYNA,the calculation accuracy of the model and algorithm is proved.Furthermore,the failure process of the concrete block under quasi-static loading is demonstrated,showing the basic dynamic tran- sitional process from continuum to non-continuum.The results of calculation can be displayed by animation.The damage modes are similar to the experimental results.The two numerical examples above prove that our model and its code are powerful and efficient in simulating the dynamic failure problems accompanying the transition from continuum to non-continuum.It also shows that the discrete element method(DEM)will have broad prospects for development and application.展开更多
Multibody system dynamics provides a strong tool for the estimation of dynamic performances and the optimization of multisystem robot design. It can be described with differential algebraic equations(DAEs). In this pa...Multibody system dynamics provides a strong tool for the estimation of dynamic performances and the optimization of multisystem robot design. It can be described with differential algebraic equations(DAEs). In this paper, a particle swarm optimization(PSO) method is introduced to solve and control a symplectic multibody system for the first time. It is first combined with the symplectic method to solve problems in uncontrolled and controlled robotic arm systems. It is shown that the results conserve the energy and keep the constraints of the chaotic motion, which demonstrates the efficiency, accuracy, and time-saving ability of the method. To make the system move along the pre-planned path, which is a functional extremum problem, a double-PSO-based instantaneous optimal control is introduced. Examples are performed to test the effectiveness of the double-PSO-based instantaneous optimal control. The results show that the method has high accuracy, a fast convergence speed, and a wide range of applications.All the above verify the immense potential applications of the PSO method in multibody system dynamics.展开更多
As a novel kind of particle method for explicit dynamics,the finite particle method(FPM)does not require the formation or solution of global matrices,and the evaluations of the element equivalent forces and particle d...As a novel kind of particle method for explicit dynamics,the finite particle method(FPM)does not require the formation or solution of global matrices,and the evaluations of the element equivalent forces and particle displacements are decoupled in nature,thus making this method suitable for parallelization.The FPM also requires an acceleration strategy to overcome the heavy computational burden of its explicit framework for time-dependent dynamic analysis.To this end,a GPU-accelerated parallel strategy for the FPM is proposed in this paper.By taking advantage of the independence of each step of the FPM workflow,a generic parallelized computational framework for multiple types of analysis is established.Using the Compute Unified Device Architecture(CUDA),the GPU implementations of the main tasks of the FPM,such as evaluating and assembling the element equivalent forces and solving the kinematic equations for particles,are elaborated through careful thread management and memory optimization.Performance tests show that speedup ratios of 8,25 and 48 are achieved for beams,hexahedral solids and triangular shells,respectively.For examples consisting of explicit dynamic analyses of shells and solids,comparisons with Abaqus using 1 to 8 CPU cores validate the accuracy of the results and demonstrate a maximum speed improvement of a factor of 11.2.展开更多
In this paper, by defining new state vectors and developing new transfer matrices of various elements mov- ing in space, the discrete time transfer matrix method of multi-rigid-flexible-body system is expanded to stud...In this paper, by defining new state vectors and developing new transfer matrices of various elements mov- ing in space, the discrete time transfer matrix method of multi-rigid-flexible-body system is expanded to study the dynamics of multibody system with flexible beams moving in space. Formulations and numerical example of a rigid- flexible-body three pendulums system moving in space are given to validate the method. Using the new method to study the dynamics of multi-rigid-flexible-body system mov- ing in space, the global dynamics equations of system are not needed, the orders of involved matrices of the system are very low and the computational speed is high, irrespec- tive of the size of the system. The new method is simple, straightforward, practical, and provides a powerful tool for multi-rigid-flexible-body system dynamics.展开更多
The hybrid dynamics of multi-rigid-body and multi-flexible-body system becomes the mainstream of multi-body dynamics.Currently there lacks a compact approach to model the hybrid dynamics,especially in modern machine t...The hybrid dynamics of multi-rigid-body and multi-flexible-body system becomes the mainstream of multi-body dynamics.Currently there lacks a compact approach to model the hybrid dynamics,especially in modern machine tool application,due to the difficulty of solving the hybrid equations or the limitation of current software when dealing with the hybrid dynamics.The extended transfer matrix method（E-TMM）,which extends elements in three-dimensional space with higher matrixes,is proposed to simplify the modeling process of the hybrid dynamics.The E-TMM modeling approaches of 3 basic elements including 3D vibrant rigid body,joint and flexible body are studied in details.A parallel mill-turn tool spindle head unit driven by dual-linear motors is chosen as a plant to demonstrate the E-TMM modeling process.By using E-TMM,the spindle head unit is simplified as a topological network consisting of the three types of element,i.e.,3D vibrant rigid body,joint and flexible body,including 11 rigid bodies,14 joints and 1 3D-Timoshenko beam.Then the dynamic model of the system can be easily obtained by deducing the element-network by means of state vector transformation.The dynamic characteristics of the spindle head,such as natural frequencies,dynamic flexibility,etc.can be predicted by solving the obtained model.Experiment verification indicates that the E-TMM is valid with enough accuracy in the dynamic analysis of the parallel mill-turn tool spindle head.The E-TMM is capable of modeling the dynamics of machine tool structure with no requirements of deducing and solving the sophisticated differential equations.Moreover,the E-TMM provides a simple and elegant tool for hybrid dynamic analysis in future dynamic design of machine tools.展开更多
The impact problem of a flexible multibody system is a non-smooth, high-transient, and strong-nonlinear dynamic process with variable boundary. How to model the contact/impact process accurately and efficiently is one...The impact problem of a flexible multibody system is a non-smooth, high-transient, and strong-nonlinear dynamic process with variable boundary. How to model the contact/impact process accurately and efficiently is one of the main difficulties in many engineering applications. The numerical approaches being used widely in impact analysis are mainly from two fields： multibody system dynamics （MBS） and computational solid mechanics （CSM）. Approaches based on MBS provide a more efficient yet less accurate analysis of the contact/impact problems, while approaches based on CSM are well suited for particularly high accuracy needs, yet require very high computational effort. To bridge the gap between accuracy and efficiency in the dynamic simulation of a flexible multibody system with contacts/impacts, a partition method is presented considering that the contact body is divided into two parts, an impact region and a non-impact region. The impact region is modeled using the finite element method to guarantee the local accuracy, while the non-impact region is modeled using the modal reduction approach to raise the global efficiency. A three-dimensional rod-plate impact experiment is designed and performed to validate the numerical results. The principle for how to partition the contact bodies is proposed： the maximum radius of the impact region can be estimated by an analytical method, and the modal truncation orders of the non-impact region can be estimated by the highest frequency of the signal measured. The simulation results using the presented method are in good agreement with the experimental results. It shows that this method is an effec-rive formulation considering both accuracy and efficiency. Moreover, a more complicated multibody impact problem of a crank slider mechanism is investigated to strengthen this conclusion.展开更多
In this paper,Runge-Kutta Discontinuous Galerkin（RKDG） finite element method is presented to solve the onedimensional inviscid compressible gas dynamic equations in a Lagrangian coordinate.The equations are discreti...In this paper,Runge-Kutta Discontinuous Galerkin（RKDG） finite element method is presented to solve the onedimensional inviscid compressible gas dynamic equations in a Lagrangian coordinate.The equations are discretized by the DG method in space and the temporal discretization is accomplished by the total variation diminishing Runge-Kutta method.A limiter based on the characteristic field decomposition is applied to maintain stability and non-oscillatory property of the RKDG method.For multi-medium fluid simulation,the two cells adjacent to the interface are treated differently from other cells.At first,a linear Riemann solver is applied to calculate the numerical ？ux at the interface.Numerical examples show that there is some oscillation in the vicinity of the interface.Then a nonlinear Riemann solver based on the characteristic formulation of the equation and the discontinuity relations is adopted to calculate the numerical ？ux at the interface,which suppresses the oscillation successfully.Several single-medium and multi-medium fluid examples are given to demonstrate the reliability and efficiency of the algorithm.展开更多
The difficulty in solving stochastic dynamics problems lies in the need for a large number of repeated computations of deterministic dynamic equations,which has been a challenge in stochastic dynamics analysis and was...The difficulty in solving stochastic dynamics problems lies in the need for a large number of repeated computations of deterministic dynamic equations,which has been a challenge in stochastic dynamics analysis and was discussed in this study.To efficiently and accurately compute the exponential of the dynamics state matrix and the matrix functions due to external loads,an adaptively filtered precise integration method was proposed,which inherits the high precision of the precise integrationmethod,improves the computational efficiency and saves the memory required.Moreover,the perturbation method was introduced to avoid repeated computations of matrix exponential and terms due to external loads.Based on the filtering and perturbation techniques,an adaptively filtered precise integration method considering perturbation for stochastic dynamics problems was developed.Two numerical experiments,including a model of phononic crystal and a bridge model considering random parameters,were performed to test the performance of the proposed method in terms of accuracy and efficiency.Numerical results show that the accuracy and efficiency of the proposed method are better than those of the existing precise integration method,the Newmark-βmethod and the Wilson-θmethod.展开更多
This paper presents a computationally efficient real-time trajectory planning framework for typical unmanned combat aerial vehicle （UCAV） performing autonomous air-to-surface （A/S） attack. It combines the benefits...This paper presents a computationally efficient real-time trajectory planning framework for typical unmanned combat aerial vehicle （UCAV） performing autonomous air-to-surface （A/S） attack. It combines the benefits of inverse dynamics optimization method and receding horizon optimal control technique. Firstly, the ground attack trajectory planning problem is mathematically formulated as a receding horizon optimal control problem （RHC-OCP）. In particular, an approximate elliptic launch acceptable region （LAR） model is proposed to model the critical weapon delivery constraints. Secondly, a planning algorithm based on inverse dynamics optimization, which has high computational efficiency and good convergence properties, is developed to solve the RHCOCP in real-time. Thirdly, in order to improve robustness and adaptivity in a dynamic and uncer- tain environment, a two-degree-of-freedom （2-DOF） receding horizon control architecture is introduced and a regular real-time update strategy is proposed as well, and the real-time feedback can be achieved and the not-converged situations can be handled. Finally, numerical simulations demon- strate the efficiency of this framework, and the results also show that the presented technique is well suited for real-time implementation in dynamic and uncertain environment.展开更多
Multi-body dynamics,relative coordinates and graph theory are combined to analyze the structure of a vehicle suspension.The dynamic equations of the left front suspension system are derived for modeling.First,The pure...Multi-body dynamics,relative coordinates and graph theory are combined to analyze the structure of a vehicle suspension.The dynamic equations of the left front suspension system are derived for modeling.First,The pure tire theory model is used as the input criteria of the suspension multibody system dynamic model in order to simulate the suspension K&C characteristics test.Then,it is important to verify the accuracy of this model by comparing and analyzing the experimental data and simulation results.The results show that the model has high precision and can predict the performance of the vehicle.It also provides a new solution for the vehicle dynamic modeling.展开更多
The dynamics, stability and control problem of a kind of infinite dimensional system are studied in the functional space with the method of modern Mathematics. First, the dynamical control model of the distributed par...The dynamics, stability and control problem of a kind of infinite dimensional system are studied in the functional space with the method of modern Mathematics. First, the dynamical control model of the distributed parameter system with multi-body flexible and multi-topological structure was established which has damping, gyroscopic parts and constrained damping. Secondly, the necessary and sufficient condition of controllability and observability, the stability theory and asymptotic property of the system were obtained. These results expand the theory of the field about the dynamics and control of the system with multi-body flexible structure, and have important engineering significance.展开更多
The decay dynamic of an excited quantum emitter(QE)is one of the most important contents in quantum optics.It has been widely applied in the field of quantum computing and quantum state manipulation.When the electroma...The decay dynamic of an excited quantum emitter(QE)is one of the most important contents in quantum optics.It has been widely applied in the field of quantum computing and quantum state manipulation.When the electromagnetic environment is described by several pseudomodes,the effective Hamiltonian method based on the multi-mode Jaynes-Cummings model provides a clear physical picture and a simple and convenient way to solve the decay dynamics.However,in previous studies,only the resonant modes are taken into account,while the non-resonant contributions are ignored.In this work,we study the applicability and accuracy of the effective Hamiltonian method for the decay dynamics.We consider different coupling strengths between a two-level QE and a gold nanosphere.The results for dynamics by the resolvent operator technique are used as a reference.Numerical results show that the effective Hamiltonian method provides accurate results when the two-level QE is resonant with the plasmon.However,when the detuning is large,the effective Hamiltonian method is not accurate.In addition,the effective Hamiltonian method cannot be applied when there is a bound state between the QE and the plasmon.These results are of great significance to the study of the decay dynamics in micro-nano structures described by quasi-normal modes.展开更多
基金supported by the National Natural Science Foundation of China (Grant No. 51175484)the Science Foundation of Shandong Province (Grant No. ZR2010EM052)
文摘Deployment of buoy systems is one of the most important procedures for the operation of buoy system. In the present study, a single-point mooring buoy system which contains surface buoy, cable segments with components, anchor and so on is modeled by applying multi-body dynamics method. The motion equations are developed in discrete node description and fully Cartesian coordinates. Then numerical method is used to solve the ordinary differential equations and dynamics simulations are achieved while anchor is casting from board. The trajectories and velocities of different nodes without current and with current in buoy system are obtained. The transient tension force of each part of the cable is analyzed in the process of deployment. Numerical results indicate that the transient payload increases to a peak value when the anchor is touching the seabed and the maximum tension force will vary with different floating configuration. This work is helpful for design and deployment planning of buoy system.
基金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.
文摘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(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.
基金Important National Science & Technology Specific Projects (2009ZX04001-073)National Natural Science Foundation of China (51105025)
文摘In order to analyze the influence of configuration parameters on dynamic characteristics of machine tools in the working space, the configuration parameters have been suggested based on the orthogonal experiment method. Dynamic analysis of a milling machine, which is newly designed for producing turbine blades, has been conducted by utilizing the modal synthesis method. The finite element model is verified and updated by experimental modal analysis （EMA） of the machine tool. The result gained by modal synthesis method is compared with whole-model finite element method （FEM） result as well. According to the orthogonal experiment method, four configuration parameters of machine tool are considered as four factors for dynamic characteristics. The influence of configuration parameters on the first three natural frequencies is obtained by range analysis. It is pointed out that configuration parameter is the most important factor affecting the fundamental frequency of machine tools, and configuration parameter has less effect on lower-order modes of the system than others. The combination of configuration parameters which makes the fundamental frequency reach the maximum value is provided. Through demonstration, the conclusion can be drawn that the influence of configuration parameters on the natural frequencies of machine tools can be analyzed explicitly by the orthogonal experiment method, which offers a new method for estimating the dynamic characteristics of machine tools.
基金supported by the National Natural Science Foundation of China(Nos.11132007,11272155,and 10772085)the Fundamental Research Funds for the Central Universities(No.30920130112009)the 333 Project of Jiangsu Province of China(No.BRA2011172)
文摘The impact dynamics of a flexible multibody system is investigated. By using a partition method, the system is divided into two parts, the local impact region and the region away from the impact. The two parts are connected by specific boundary conditions, and the system after partition is equivalent to the original system. According to the rigid-flexible coupling dynamic theory of multibody system, system＇s rigid-flexible coupling dynamic equations without impact are derived. A local impulse method for establishing the initial impact conditions is proposed. It satisfies the compatibility con- ditions for contact constraints and the actual physical situation of the impact process of flexible bodies. Based on the contact constraint method, system＇s impact dynamic equa- tions are derived in a differential-algebraic form. The contact/separation criterion and the algorithm are given. An impact dynamic simulation is given. The results show that system＇s dynamic behaviors including the energy, the deformations, the displacements, and the impact force during the impact process change dramatically. The impact makes great effects on the global dynamics of the system during and after impact.
基金supported by the National Natural Science Foundation of China (Grant No. 51009092)the Doctoral Foundation of Education Ministry of China (Grant No. 20090073120013)the Scientific Research Foundation of State Education Ministry for the Returned Overseas Chinese Scholars
文摘In re-entry, the drilling riser hanging to the holding vessel takes on a free hanging state, waiting to be moved from the initial random position to the wellhead. For the re-entry, dynamics calculation is often done to predict the riser motion or evaluate the structural safety. A dynamics calculation method based on Flexible Segment Model （FSM） is proposed for free hanging marine risers. In FSM, a riser is discretized into a series of flexible segments. For each flexible segment, its deflection feature and external forces are analyzed independently. For the whole riser, the nonlinear governing equations are listed according to the moment equilibrium at nodes. For the solution of the nonlinear equations, a linearization iteration scheme is provided in the paper. Owing to its flexibility, each segment can match a long part of the riser body, which enables that good results can be obtained even with a small number of segments. Moreover, the linearization iteration scheme can avoid widely used Newton-Rapson iteration scheme in which the calculation stability is influenced by the initial points. The FSM-based dynamics calculation is timesaving and stable, so suitable for the shape prediction or real-time control of free hanging marine risers.
基金Supported by National Natural Science Foundation of China(Grant No.51175422)
文摘There are vast constraint equations in conventional dynamics analysis of deployable structures,which lead to differential-algebraic equations(DAEs)solved hard.To reduce the difficulty of solving and the amount of equations,a new flexible multibody dynamics analysis methodology of deployable structures with scissor-like elements(SLEs)is presented.Firstly,a precise model of a flexible bar of SLE is established by the higher order shear deformable beam element based on the absolute nodal coordinate formulation(ANCF),and the master/slave freedom method is used to obtain the dynamics equations of SLEs without constraint equations.Secondly,according to features of deployable structures,the specification matrix method(SMM)is proposed to eliminate the constraint equations among SLEs in the frame of ANCF.With this method,the inner and the boundary nodal coordinates of element characteristic matrices can be separated simply and efficiently,especially on condition that there are vast nodal coordinates.So the element characteristic matrices can be added end to end circularly.Thus,the dynamic model of deployable structure reduces dimension and can be assembled without any constraint equation.Next,a new iteration procedure for the generalized-a algorithm is presented to solve the ordinary differential equations(ODEs)of deployable structure.Finally,the proposed methodology is used to analyze the flexible multi-body dynamics of a planar linear array deployable structure based on three scissor-like elements.The simulation results show that flexibility has a significant influence on the deployment motion of the deployable structure.The proposed methodology indeed reduce the difficulty of solving and the amount of equations by eliminating redundant degrees of freedom and the constraint equations in scissor-like elements and among scissor-like elements.
基金Project supported by the National Natural Science Foundation of China(Nos.59978005 and 10232024)the National Distinguished Youth Fund of China(No.10025212).
文摘A three-dimensional discrete element model of the connective type is presented. Moreover,a three-dimensional numerical analysis code,which can carry out the transitional pro- cess from connective model(for continuum)to contact model(for non-continuum),is developed for simulating the mechanical process from continuum to non-continuum.The wave propagation process in a concrete block(as continuum)made of cement grout under impact loading is numer- ically simulated with this code.By comparing its numerical results with those by LS-DYNA,the calculation accuracy of the model and algorithm is proved.Furthermore,the failure process of the concrete block under quasi-static loading is demonstrated,showing the basic dynamic tran- sitional process from continuum to non-continuum.The results of calculation can be displayed by animation.The damage modes are similar to the experimental results.The two numerical examples above prove that our model and its code are powerful and efficient in simulating the dynamic failure problems accompanying the transition from continuum to non-continuum.It also shows that the discrete element method(DEM)will have broad prospects for development and application.
基金Project supported by the National Natural Science Foundation of China(Nos.91648101 and11672233)the Northwestern Polytechnical University(NPU)Foundation for Fundamental Research(No.3102017AX008)the National Training Program of Innovation and Entrepreneurship for Undergraduates(No.S201710699033)
文摘Multibody system dynamics provides a strong tool for the estimation of dynamic performances and the optimization of multisystem robot design. It can be described with differential algebraic equations(DAEs). In this paper, a particle swarm optimization(PSO) method is introduced to solve and control a symplectic multibody system for the first time. It is first combined with the symplectic method to solve problems in uncontrolled and controlled robotic arm systems. It is shown that the results conserve the energy and keep the constraints of the chaotic motion, which demonstrates the efficiency, accuracy, and time-saving ability of the method. To make the system move along the pre-planned path, which is a functional extremum problem, a double-PSO-based instantaneous optimal control is introduced. Examples are performed to test the effectiveness of the double-PSO-based instantaneous optimal control. The results show that the method has high accuracy, a fast convergence speed, and a wide range of applications.All the above verify the immense potential applications of the PSO method in multibody system dynamics.
基金the financial support provided by the National Key Research and Development Program of China(Grant No.2016YFC0800200)the National Natural Science Foundation of China(Grant Nos.51578494 and 51778568)the Fundamental Research Funds for the Central Universities(Grant No.2019QNA4043).
文摘As a novel kind of particle method for explicit dynamics,the finite particle method(FPM)does not require the formation or solution of global matrices,and the evaluations of the element equivalent forces and particle displacements are decoupled in nature,thus making this method suitable for parallelization.The FPM also requires an acceleration strategy to overcome the heavy computational burden of its explicit framework for time-dependent dynamic analysis.To this end,a GPU-accelerated parallel strategy for the FPM is proposed in this paper.By taking advantage of the independence of each step of the FPM workflow,a generic parallelized computational framework for multiple types of analysis is established.Using the Compute Unified Device Architecture(CUDA),the GPU implementations of the main tasks of the FPM,such as evaluating and assembling the element equivalent forces and solving the kinematic equations for particles,are elaborated through careful thread management and memory optimization.Performance tests show that speedup ratios of 8,25 and 48 are achieved for beams,hexahedral solids and triangular shells,respectively.For examples consisting of explicit dynamic analyses of shells and solids,comparisons with Abaqus using 1 to 8 CPU cores validate the accuracy of the results and demonstrate a maximum speed improvement of a factor of 11.2.
基金supported by the Natural Science Foundation of China Government (10902051)the Natural Science Foundation of Jiangsu Province (BK2008046)the German Science Foundation
文摘In this paper, by defining new state vectors and developing new transfer matrices of various elements mov- ing in space, the discrete time transfer matrix method of multi-rigid-flexible-body system is expanded to study the dynamics of multibody system with flexible beams moving in space. Formulations and numerical example of a rigid- flexible-body three pendulums system moving in space are given to validate the method. Using the new method to study the dynamics of multi-rigid-flexible-body system mov- ing in space, the global dynamics equations of system are not needed, the orders of involved matrices of the system are very low and the computational speed is high, irrespec- tive of the size of the system. The new method is simple, straightforward, practical, and provides a powerful tool for multi-rigid-flexible-body system dynamics.
基金supported by National Key Technology R&D Program of China (Grant No. 2006BAF01B09)the Research Fund for Doctoral Program of Higher Education of China (Grant No. 200800060010)
文摘The hybrid dynamics of multi-rigid-body and multi-flexible-body system becomes the mainstream of multi-body dynamics.Currently there lacks a compact approach to model the hybrid dynamics,especially in modern machine tool application,due to the difficulty of solving the hybrid equations or the limitation of current software when dealing with the hybrid dynamics.The extended transfer matrix method（E-TMM）,which extends elements in three-dimensional space with higher matrixes,is proposed to simplify the modeling process of the hybrid dynamics.The E-TMM modeling approaches of 3 basic elements including 3D vibrant rigid body,joint and flexible body are studied in details.A parallel mill-turn tool spindle head unit driven by dual-linear motors is chosen as a plant to demonstrate the E-TMM modeling process.By using E-TMM,the spindle head unit is simplified as a topological network consisting of the three types of element,i.e.,3D vibrant rigid body,joint and flexible body,including 11 rigid bodies,14 joints and 1 3D-Timoshenko beam.Then the dynamic model of the system can be easily obtained by deducing the element-network by means of state vector transformation.The dynamic characteristics of the spindle head,such as natural frequencies,dynamic flexibility,etc.can be predicted by solving the obtained model.Experiment verification indicates that the E-TMM is valid with enough accuracy in the dynamic analysis of the parallel mill-turn tool spindle head.The E-TMM is capable of modeling the dynamics of machine tool structure with no requirements of deducing and solving the sophisticated differential equations.Moreover,the E-TMM provides a simple and elegant tool for hybrid dynamic analysis in future dynamic design of machine tools.
基金supported by the National Natural Science Foundation of China (Grants 11772188, 11132007)
文摘The impact problem of a flexible multibody system is a non-smooth, high-transient, and strong-nonlinear dynamic process with variable boundary. How to model the contact/impact process accurately and efficiently is one of the main difficulties in many engineering applications. The numerical approaches being used widely in impact analysis are mainly from two fields： multibody system dynamics （MBS） and computational solid mechanics （CSM）. Approaches based on MBS provide a more efficient yet less accurate analysis of the contact/impact problems, while approaches based on CSM are well suited for particularly high accuracy needs, yet require very high computational effort. To bridge the gap between accuracy and efficiency in the dynamic simulation of a flexible multibody system with contacts/impacts, a partition method is presented considering that the contact body is divided into two parts, an impact region and a non-impact region. The impact region is modeled using the finite element method to guarantee the local accuracy, while the non-impact region is modeled using the modal reduction approach to raise the global efficiency. A three-dimensional rod-plate impact experiment is designed and performed to validate the numerical results. The principle for how to partition the contact bodies is proposed： the maximum radius of the impact region can be estimated by an analytical method, and the modal truncation orders of the non-impact region can be estimated by the highest frequency of the signal measured. The simulation results using the presented method are in good agreement with the experimental results. It shows that this method is an effec-rive formulation considering both accuracy and efficiency. Moreover, a more complicated multibody impact problem of a crank slider mechanism is investigated to strengthen this conclusion.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 11261035,11171038,and 10771019)the Science Reaearch Foundation of Institute of Higher Education of Inner Mongolia Autonomous Region,China (Grant No. NJZZ12198)the Natural Science Foundation of Inner Mongolia Autonomous Region,China (Grant No. 2012MS0102)
文摘In this paper,Runge-Kutta Discontinuous Galerkin（RKDG） finite element method is presented to solve the onedimensional inviscid compressible gas dynamic equations in a Lagrangian coordinate.The equations are discretized by the DG method in space and the temporal discretization is accomplished by the total variation diminishing Runge-Kutta method.A limiter based on the characteristic field decomposition is applied to maintain stability and non-oscillatory property of the RKDG method.For multi-medium fluid simulation,the two cells adjacent to the interface are treated differently from other cells.At first,a linear Riemann solver is applied to calculate the numerical ？ux at the interface.Numerical examples show that there is some oscillation in the vicinity of the interface.Then a nonlinear Riemann solver based on the characteristic formulation of the equation and the discontinuity relations is adopted to calculate the numerical ？ux at the interface,which suppresses the oscillation successfully.Several single-medium and multi-medium fluid examples are given to demonstrate the reliability and efficiency of the algorithm.
基金the support of the National Natural Science Foundation of China(Grant Nos.11472067 and 51609034)the Science Foundation of Liaoning Province of China(No.2021-MS-119)+1 种基金the Dalian Youth Science and Technology Star Project(No.2018RQ06)the Fundamental Research Funds for the Central Universities(Grant No.DUT20GJ216).
文摘The difficulty in solving stochastic dynamics problems lies in the need for a large number of repeated computations of deterministic dynamic equations,which has been a challenge in stochastic dynamics analysis and was discussed in this study.To efficiently and accurately compute the exponential of the dynamics state matrix and the matrix functions due to external loads,an adaptively filtered precise integration method was proposed,which inherits the high precision of the precise integrationmethod,improves the computational efficiency and saves the memory required.Moreover,the perturbation method was introduced to avoid repeated computations of matrix exponential and terms due to external loads.Based on the filtering and perturbation techniques,an adaptively filtered precise integration method considering perturbation for stochastic dynamics problems was developed.Two numerical experiments,including a model of phononic crystal and a bridge model considering random parameters,were performed to test the performance of the proposed method in terms of accuracy and efficiency.Numerical results show that the accuracy and efficiency of the proposed method are better than those of the existing precise integration method,the Newmark-βmethod and the Wilson-θmethod.
基金supported by the National Defense Foundation of China(No.403060103)
文摘This paper presents a computationally efficient real-time trajectory planning framework for typical unmanned combat aerial vehicle （UCAV） performing autonomous air-to-surface （A/S） attack. It combines the benefits of inverse dynamics optimization method and receding horizon optimal control technique. Firstly, the ground attack trajectory planning problem is mathematically formulated as a receding horizon optimal control problem （RHC-OCP）. In particular, an approximate elliptic launch acceptable region （LAR） model is proposed to model the critical weapon delivery constraints. Secondly, a planning algorithm based on inverse dynamics optimization, which has high computational efficiency and good convergence properties, is developed to solve the RHCOCP in real-time. Thirdly, in order to improve robustness and adaptivity in a dynamic and uncer- tain environment, a two-degree-of-freedom （2-DOF） receding horizon control architecture is introduced and a regular real-time update strategy is proposed as well, and the real-time feedback can be achieved and the not-converged situations can be handled. Finally, numerical simulations demon- strate the efficiency of this framework, and the results also show that the presented technique is well suited for real-time implementation in dynamic and uncertain environment.
基金Supported by the National Key Research and Development Program of China(2017YFB0103801)
文摘Multi-body dynamics,relative coordinates and graph theory are combined to analyze the structure of a vehicle suspension.The dynamic equations of the left front suspension system are derived for modeling.First,The pure tire theory model is used as the input criteria of the suspension multibody system dynamic model in order to simulate the suspension K&C characteristics test.Then,it is important to verify the accuracy of this model by comparing and analyzing the experimental data and simulation results.The results show that the model has high precision and can predict the performance of the vehicle.It also provides a new solution for the vehicle dynamic modeling.
文摘The dynamics, stability and control problem of a kind of infinite dimensional system are studied in the functional space with the method of modern Mathematics. First, the dynamical control model of the distributed parameter system with multi-body flexible and multi-topological structure was established which has damping, gyroscopic parts and constrained damping. Secondly, the necessary and sufficient condition of controllability and observability, the stability theory and asymptotic property of the system were obtained. These results expand the theory of the field about the dynamics and control of the system with multi-body flexible structure, and have important engineering significance.
基金Project supported by the National Natural Science Foundation of China(11964010,11564013 and 11464014)the Natural Science Foundation of Hunan Province(2020JJ4495)+1 种基金the Scientific Research Fund of Hunan Provincial Education Department(22A0377 and 21A0333)the Jishou University Innovation Foundation for Postgraduate(Jdy20038)。
文摘The decay dynamic of an excited quantum emitter(QE)is one of the most important contents in quantum optics.It has been widely applied in the field of quantum computing and quantum state manipulation.When the electromagnetic environment is described by several pseudomodes,the effective Hamiltonian method based on the multi-mode Jaynes-Cummings model provides a clear physical picture and a simple and convenient way to solve the decay dynamics.However,in previous studies,only the resonant modes are taken into account,while the non-resonant contributions are ignored.In this work,we study the applicability and accuracy of the effective Hamiltonian method for the decay dynamics.We consider different coupling strengths between a two-level QE and a gold nanosphere.The results for dynamics by the resolvent operator technique are used as a reference.Numerical results show that the effective Hamiltonian method provides accurate results when the two-level QE is resonant with the plasmon.However,when the detuning is large,the effective Hamiltonian method is not accurate.In addition,the effective Hamiltonian method cannot be applied when there is a bound state between the QE and the plasmon.These results are of great significance to the study of the decay dynamics in micro-nano structures described by quasi-normal modes.
