• International Journal of Naval Architecture and Ocean Engineering
• /
• v.8 no.4
• /
• pp.320-329
• /
• 2016

An implicit method of solving the six degree-of-freedom rigid body motion equations based on the second order Adams-Bashforth-Moulten method was utilised as an improvement over the leapfrog scheme by making modifications to the rigid body motion solver libraries directly. The implementation will depend on predictor-corrector steps still residing within the hybrid Pressure Implicit with Splitting of Operators - Semi-Implicit Method for Pressure Linked Equations (PIMPLE) outer corrector loops to ensure strong coupling between fluid and motion. Aitken's under-relaxation is also introduced in this study to optimise the convergence rate and stability of the coupled solver. The resulting coupled solver ran on a free floating object tutorial test case when converged matches the original solver. It further allows a varying 70%-80% reduction in simulation times compared using a fixed under-relaxation to achieve the required stability.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.33 no.12
• /
• pp.1433-1441
• /
• 2009

In multibody dynamic analysis, one of the most important problems is to reduce computation times for real-time simulation. This paper presents the derivation procedure of equations of motion of a 6${\times}$6 autonomous vehicle in terms of chassis local coordinates which do not require coordinates transformation matrix to enhance efficiency for real-time dynamic analysis. Also, equations of motion are derived using the VT(velocity transformation) technique and symbolic computation method coded by MATLAB. The Jacobian matrix of the equations of motion of a system is derived from symbolic operations to apply the implicit integration method. The analysis results were compared with ADAMS results to verify the accuracy and approve the feasibility of real time analysis.

• Structural Engineering and Mechanics
• /
• v.38 no.3
• /
• pp.283-300
• /
• 2011

A new conception of fundamental tasks in dynamics of the bridge-track-train systems (BTT), with the aim to evaluate moving load's models adequacy, has been developed. The 2D physical models of BTT systems, corresponding to the fundamental tasks, have been worked out taking into account one-way constraints between the moving unsprung masses and the track. A method for deriving the implicit equations of motion, governing vibrations of BTT systems' models, as well as algorithms for numerical integration of these equations, leading to the solutions of high accuracy and relatively short times of simulations, have been also developed. The derived equations and formulated algorithms constitute the basis for numerical simulation of vibrations of the considered systems.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2006.04a
• /
• pp.389-394
• /
• 2006

Fully flexible cell preserves Hamiltonian in structure, so the symplectic time integrator is applied to the equations of motion. Primarily, generalized leapfrog time integration (GLF) is applicable, but the equations of motion by GLF have some of implicit formulas. The implicit formulas give rise to a complicate calculation for coding and need an iteration process. In this paper, the time integration formulas are obtained for the fully flexible cell molecular dynamics simulation by using the splitting time integration. It separates flexible cell Hamiltonian into terms corresponding to each of Hamiltonian term, so the simple and completely explicit recursion formula was obtained. The explicit formulas are easy to implementation for coding and may be reduced the integration time because they are not need iteration process. We are going to compare the resulting splitting time integration with the implicit generalized leapfrog time integration.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.3
• /
• pp.8-16
• /
• 2002

The implicit formulation of rotor dynamics for helicopter flight simulation has been derived and and presented. The generalized vector kinematics regarding the relative motion between coordinates were expressed as a unified matrix operation and applied to get the inertial velocities and accelerations at arbitaty rotor blade span position. Based on these results the rotor aeromechanic equations for flapping dynamics, lead-lag dynamics and torque dynamics were formulated as an implicit form. Spatial integration methods of rotor dynamic equations along blade span and the expanded applicability of the present implicit formulations for arbitrary hings geometry and hinge sequences have been investigated. Time integration methods for present DAE(Differential Algebraic Equation) to calculate dynamic response calculation are recommenaded as future works.

• Proceedings of the KIEE Conference
• /
• 1988.07a
• /
• pp.950-953
• /
• 1988

The equations of motion for linearly elastic bodies undergoing large displacement motion are derived. This produces a set of equations which are efficient to numerically integrate. The equations for the elastic bodies are formulated and simplified to provide as much efficiency as possible in their numerical solution. A futher efficiency is obtained through the use of floating reference frame. The equation are presented in two forms for numerical integration. 1) Explicit numerical integration 2) Implicit numerical integration. In this paper, there was used the numerical integration. The implicit numerical integration is extended to solved second order equation, futher reducing the numerical effort required. The formulation given is seen to be occulate and is expected to be efficient for many types of problems.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.12
• /
• pp.1966-1973
• /
• 2001

The Jacobian matrix of the equations of motion of a system using velocity transformation technique is derived via variation methods to apply the implicit integration algorithm, DASSL. The concept of generalized coordinate partitioning is used to parameterize the constraint set with independent generalized coordinates. DASSL is applied to determine independent generalized coordinates and velocities. Dependent generalized coordinates, velocities, accelerations and Lagrange multipliers are explicitly retained in the formulation to satisfy all of the governing kinematic and dynamic equations. The derived Jacobian matrix of a system is proved to be valid and accurate both analytically and through solution of numerical examples.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.31 no.8
• /
• pp.826-832
• /
• 2007

Fully flexible cell preserves Hamiltonian in structure so that the symplectic time integrator is applicable to the equations of motion. In the direct formulation of fully flexible cell N-Sigma-T ensemble, a generalized leapfrog time integration (GLF) is applicable for fully flexible cell simulation, but the equations of motion by GLF has structure of implicit algorithm. In this paper, the time integration formula is derived for the fully flexible cell molecular dynamics simulation by using the splitting time integration. It separates flexible cell Hamiltonian into terms corresponding to each of Hamiltonian term. Thus the simple and completely explicit recursion formula was obtained. We compare the performance and the result of present splitting time integration with those of the implicit generalized leapfrog time integration.

• Transactions of the Korean Society of Machine Tool Engineers
• /
• v.10 no.2
• /
• pp.18-28
• /
• 2001

Super computers has been utilized to carry out vehicle dynamics in real time. This research propose an implicit integra-tion method for vehicle state variables. Newton chord method is empolyed to solve the equations of motion and con-straints. The equations of motion and constraints are formulated such that the Jacobian matrix for Newton chord method is needed to be computed only once for a dynamic analysis. Numerical experiments showed that the Jacobian matrix generat-ed at the initial time could have been utilized for the Newton chord iterations throughout simulations under various driving conditions. Convergence analysis of Newton chord method with the proposed Jacobian updating method is carried out. The proposed algorithm yielded accurate solutions for a prototype vehicle multibody model in realtime on a 400 MHz PC compatible.

• Structural Engineering and Mechanics
• /
• v.75 no.3
• /
• pp.357-367
• /
• 2020

The objective of this article is investigation of dynamic response of thick multilayer functionally graded (FG) beam under generalized dynamic forces. The plane stress problem is exploited to describe the constitutive equation of thick FG beam to get realistic and accurate response. Applied dynamic forces are assumed to be sinusoidal harmonic, sinusoidal pulse or triangle in time domain and point load. Equations of motion of deep FG beam are derived based on the Hamilton principle from kinematic relations and constitutive equations of plane stress problem. The numerical finite element procedure is adopted to discretize the space domain of structure and transform partial differential equations of motion to ordinary differential equations in time domain. Numerical time integration method is used to solve the system of equations in time domain and find the time responses. Numerical parametric studies are performed to illustrate effects of force type, graduation parameter, geometrical and stacking sequence of layers on the time response of deep multilayer FG beams.