A tensor method for the derivation of the equations of rigid body dynamics, based on the concepts of continuum mechanics, is presented. The formula of time derivative of the inertia tensor with zero corotational rate ...A tensor method for the derivation of the equations of rigid body dynamics, based on the concepts of continuum mechanics, is presented. The formula of time derivative of the inertia tensor with zero corotational rate is used to prove the equivalences of five methods, namely, Lagrange's equations, Nielsen's equations, Gibbs-Appell's equations, Kane's equations and the generalized momentum type of Kane's equations. Some differential identities on angular velocity and angular acceleration are given.展开更多
In Fluid Structure Interaction(FSI) problems encountered in marine hydrodynamics, the pressure field and the velocity of the rigid body are tightly coupled. This coupling is traditionally resolved in a partitioned man...In Fluid Structure Interaction(FSI) problems encountered in marine hydrodynamics, the pressure field and the velocity of the rigid body are tightly coupled. This coupling is traditionally resolved in a partitioned manner by solving the rigid body motion equations once per nonlinear correction loop, updating the position of the body and solving the fluid flow equations in the new configuration. The partitioned approach requires a large number of nonlinear iteration loops per time–step. In order to enhance the coupling, a monolithic approach is proposed in Finite Volume(FV) framework,where the pressure equation and the rigid body motion equations are solved in a single linear system. The coupling is resolved by solving the rigid body motion equations once per linear solver iteration of the pressure equation, where updated pressure field is used to calculate new forces acting on the body, and by introducing the updated rigid body boundary velocity in to the pressure equation. In this paper the monolithic coupling is validated on a simple 2D heave decay case. Additionally, the method is compared to the traditional partitioned approach(i.e. "strongly coupled" approach) in terms of computational efficiency and accuracy. The comparison is performed on a seakeeping case in regular head waves, and it shows that the monolithic approach achieves similar accuracy with fewer nonlinear correctors per time–step. Hence, significant savings in computational time can be achieved while retaining the same level of accuracy.展开更多
In this paper,we describe a modified RATTLE(M-RATTLE)method for rigid body dynamics directly in Cartesian coordinates.The M-RATTLE method introduces a new way of resetting the coordinates to satisfy the constraints at...In this paper,we describe a modified RATTLE(M-RATTLE)method for rigid body dynamics directly in Cartesian coordinates.The M-RATTLE method introduces a new way of resetting the coordinates to satisfy the constraints at each step,which is designed for the rigid body dynamics calculations in the Cartesian coordinates.M-RATTLE is algebraically equivalent to the RATTLE method and the cost of performing rigid body dynamics by M-RATTLE is independent of the number of constraints.The interaction forces between atoms belonging to the same rigid molecule do not need to be computed and explicit expressions of the constraints of internal degrees of freedom are unnecessary.The performance and sampling results of the proposed method are compared with those of the symplectic splitting method for an isolated rigid benz molecule and for a cluster of twenty-seven benz molecules.展开更多
For the dynamics of a rigid body with a fixed point based on the quaternion and the corresponding generalized momenta, a displacement-based symplectic integration scheme for differential-algebraic equations is propose...For the dynamics of a rigid body with a fixed point based on the quaternion and the corresponding generalized momenta, a displacement-based symplectic integration scheme for differential-algebraic equations is proposed and applied to the Lagrange's equations based on dependent generalized momenta. Numerical experiments show that the algorithm possesses such characters as high precision and preserving system invariants. More importantly, the generalized momenta based Lagrange's equations show unique advantages over the traditional Lagrange's equations in symplectic integrations.展开更多
The distinct element method(DEM)has been used successfully for the dynamic analysis of rigid block sys- tems.One of many difficulties associated with DEM is modeling of damping.In this paper,new procedures are propose...The distinct element method(DEM)has been used successfully for the dynamic analysis of rigid block sys- tems.One of many difficulties associated with DEM is modeling of damping.In this paper,new procedures are proposed for the damping modeling and its numerical implementation in distinct element analysis of rigid muhi-block systems.The stiff- ness proportional damping is constructed for the prescribed damping ratio,based on the non-zero fundamental frequency ef- fective during the time interval while the boundary conditions remain essentially constant.At this time interval,the funda- mental frequency can be estimated without complete eigenvalue analysis.The damping coefficients will vary while the damp- ing ratio remains the same throughout the entire analysis.A new numerical procedure is developed to prevent unnecessary energy loss that can occur during the separation phases.These procedures were implemented in the development of the dis- tinet element method for the dynamic analyses of piled multi-block systems.The analysis results |or the single-block and two-block systems were in a good agreement with the analytic predictions.Applications to the seismic analyses of piled four- block systems revealed that the new procedures can make a significant difference and may lead to much-improved results.展开更多
In this paper the method of design of kinematical and dynamical equations of mechanical systems, applied to numerical ealization, is proposed. The corresponding difference equations, which are obtained, give a guarant...In this paper the method of design of kinematical and dynamical equations of mechanical systems, applied to numerical ealization, is proposed. The corresponding difference equations, which are obtained, give a guarantee of computations with a given precision. The equations of programmed constraints and those of constraint perturbations are defined. The stability of the programmed manifold for numerical solutions of the kinematical and dynamical equations is obtained by corresponding construction of the constraint perturbation equations. The dynamical equations of system with programmed constraints are set up in the form of Lagrange’s equations in generalized coordinates. Certain inverse problems of rigid body dynamics are examined.展开更多
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 orthogonality of eigenvector is a precondition to compute the dynamic responses of linear multi-rigid-flexible-body system using the classical modal analysis method. For a linear multi-rigid-flexible-body system, ...The orthogonality of eigenvector is a precondition to compute the dynamic responses of linear multi-rigid-flexible-body system using the classical modal analysis method. For a linear multi-rigid-flexible-body system, the eigenfunction does not satisfy the orthogonality under ordinary meaning. A new concept--augmented eigenvector is introduced, which is used to overcome the orthogonality problem of eigenvectors of linear multi-rigid-flexible-body system. The constitution method and the orthogonality of augmented eigenvector are expatiated. After the orthogonality of augmented eigenvector is acquired, the coupling of coordinates in dynamics equations can be released, which makes it possible to analyze exactly the dynamic responses of linear multi-rigid-flexible-body system using the classical modal analysis method.展开更多
Based on the assumption that the slope bodies are rigid, the dynamic model of the landsiding (forward model) was put forward. According to the dynamic model, the system equations of Kalman filter were constituted. The...Based on the assumption that the slope bodies are rigid, the dynamic model of the landsiding (forward model) was put forward. According to the dynamic model, the system equations of Kalman filter were constituted. The mechanical status of a slope was hence combined with the monitoring data by Kalman filter. The model uncertainties or model errors could also be considered through some fictitious observation equations. Different from existed methods, the presented method can make use for not only the statistic information contained in the data but also the information provided by the mechanical and geological aspect of slopes. At last a numerical example was given out to show the feasibility of the method. [展开更多
文摘A tensor method for the derivation of the equations of rigid body dynamics, based on the concepts of continuum mechanics, is presented. The formula of time derivative of the inertia tensor with zero corotational rate is used to prove the equivalences of five methods, namely, Lagrange's equations, Nielsen's equations, Gibbs-Appell's equations, Kane's equations and the generalized momentum type of Kane's equations. Some differential identities on angular velocity and angular acceleration are given.
基金sponsored by Bureau Veritas under the administration of Dr.ime Malenica
文摘In Fluid Structure Interaction(FSI) problems encountered in marine hydrodynamics, the pressure field and the velocity of the rigid body are tightly coupled. This coupling is traditionally resolved in a partitioned manner by solving the rigid body motion equations once per nonlinear correction loop, updating the position of the body and solving the fluid flow equations in the new configuration. The partitioned approach requires a large number of nonlinear iteration loops per time–step. In order to enhance the coupling, a monolithic approach is proposed in Finite Volume(FV) framework,where the pressure equation and the rigid body motion equations are solved in a single linear system. The coupling is resolved by solving the rigid body motion equations once per linear solver iteration of the pressure equation, where updated pressure field is used to calculate new forces acting on the body, and by introducing the updated rigid body boundary velocity in to the pressure equation. In this paper the monolithic coupling is validated on a simple 2D heave decay case. Additionally, the method is compared to the traditional partitioned approach(i.e. "strongly coupled" approach) in terms of computational efficiency and accuracy. The comparison is performed on a seakeeping case in regular head waves, and it shows that the monolithic approach achieves similar accuracy with fewer nonlinear correctors per time–step. Hence, significant savings in computational time can be achieved while retaining the same level of accuracy.
文摘In this paper,we describe a modified RATTLE(M-RATTLE)method for rigid body dynamics directly in Cartesian coordinates.The M-RATTLE method introduces a new way of resetting the coordinates to satisfy the constraints at each step,which is designed for the rigid body dynamics calculations in the Cartesian coordinates.M-RATTLE is algebraically equivalent to the RATTLE method and the cost of performing rigid body dynamics by M-RATTLE is independent of the number of constraints.The interaction forces between atoms belonging to the same rigid molecule do not need to be computed and explicit expressions of the constraints of internal degrees of freedom are unnecessary.The performance and sampling results of the proposed method are compared with those of the symplectic splitting method for an isolated rigid benz molecule and for a cluster of twenty-seven benz molecules.
文摘For the dynamics of a rigid body with a fixed point based on the quaternion and the corresponding generalized momenta, a displacement-based symplectic integration scheme for differential-algebraic equations is proposed and applied to the Lagrange's equations based on dependent generalized momenta. Numerical experiments show that the algorithm possesses such characters as high precision and preserving system invariants. More importantly, the generalized momenta based Lagrange's equations show unique advantages over the traditional Lagrange's equations in symplectic integrations.
文摘The distinct element method(DEM)has been used successfully for the dynamic analysis of rigid block sys- tems.One of many difficulties associated with DEM is modeling of damping.In this paper,new procedures are proposed for the damping modeling and its numerical implementation in distinct element analysis of rigid muhi-block systems.The stiff- ness proportional damping is constructed for the prescribed damping ratio,based on the non-zero fundamental frequency ef- fective during the time interval while the boundary conditions remain essentially constant.At this time interval,the funda- mental frequency can be estimated without complete eigenvalue analysis.The damping coefficients will vary while the damp- ing ratio remains the same throughout the entire analysis.A new numerical procedure is developed to prevent unnecessary energy loss that can occur during the separation phases.These procedures were implemented in the development of the dis- tinet element method for the dynamic analyses of piled multi-block systems.The analysis results |or the single-block and two-block systems were in a good agreement with the analytic predictions.Applications to the seismic analyses of piled four- block systems revealed that the new procedures can make a significant difference and may lead to much-improved results.
基金Supported by Russian Fund of Fund amental Investigations(Pr.990101064)and Russian Minister of Educatin
文摘In this paper the method of design of kinematical and dynamical equations of mechanical systems, applied to numerical ealization, is proposed. The corresponding difference equations, which are obtained, give a guarantee of computations with a given precision. The equations of programmed constraints and those of constraint perturbations are defined. The stability of the programmed manifold for numerical solutions of the kinematical and dynamical equations is obtained by corresponding construction of the constraint perturbation equations. The dynamical equations of system with programmed constraints are set up in the form of Lagrange’s equations in generalized coordinates. Certain inverse problems of rigid body dynamics are examined.
基金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.
文摘The orthogonality of eigenvector is a precondition to compute the dynamic responses of linear multi-rigid-flexible-body system using the classical modal analysis method. For a linear multi-rigid-flexible-body system, the eigenfunction does not satisfy the orthogonality under ordinary meaning. A new concept--augmented eigenvector is introduced, which is used to overcome the orthogonality problem of eigenvectors of linear multi-rigid-flexible-body system. The constitution method and the orthogonality of augmented eigenvector are expatiated. After the orthogonality of augmented eigenvector is acquired, the coupling of coordinates in dynamics equations can be released, which makes it possible to analyze exactly the dynamic responses of linear multi-rigid-flexible-body system using the classical modal analysis method.
文摘Based on the assumption that the slope bodies are rigid, the dynamic model of the landsiding (forward model) was put forward. According to the dynamic model, the system equations of Kalman filter were constituted. The mechanical status of a slope was hence combined with the monitoring data by Kalman filter. The model uncertainties or model errors could also be considered through some fictitious observation equations. Different from existed methods, the presented method can make use for not only the statistic information contained in the data but also the information provided by the mechanical and geological aspect of slopes. At last a numerical example was given out to show the feasibility of the method. [