• Journal of the Korea Society for Simulation
• /
• v.18 no.2
• /
• pp.57-64
• /
• 2009

Stable and effective constraint enforcement system is one of the crucial components for physically-based dynamic simulations. This paper presents analysis and evaluation for traditional constraint enforcement systems(Lagrange Multiplier method, Baumgarte stabilization method, Post-stabilization method, Implicit constraint enforcement method, Fast projection method) to provide a guideline to users who need to integrate a suitable constraint enforcement system into their dynamic simulations. The mathematical formulations for traditional constraint enforcement systems are presented in this paper. This paper describes a summary of evaluation which consists of constraint error comparison, computational cost, and dynamic behavior analysis to verify the efficiency of each traditional constraint enforcement system.

• 제어로봇시스템학회:학술대회논문집
• /
• 1991.10a
• /
• pp.376-382
• /
• 1991

In the analysis of constrained holonomic systems, the Lagange multiplier method yields a system of second-order ordinary differential equations of motion and algebraic constraint equations. Conventional holonomic or nonholonomic constraints are defined as geometric constraints in this paper. Previous works concentrate on the geometric constraints. However, if the total energy of a dynamic system can be computed from the initial energy plus the time integral of the energy input rate due to external or internal forces, then the total energy can be artificially treated as a constraint. The violation of the total energy constraint due to numerical errors can be used as information to control these errors. It is a necessary condition for accurate simulation that both geometric and energy constraints be satisfied. When geometric constraint control is combined with energy constraint control, numerical simulation of a constrained dynamic system becomes more accurate. A new convenient and effective method to implement energy constraint control in numerical simulation is developed based on the geometric interpretation of the relation between constraints in the phase space. Several combinations of energy constraint control with either Baumgarte's Constraint Violation Stabilization Method (CVSM) are also addressed.

• Journal of Korea Multimedia Society
• /
• v.11 no.4
• /
• pp.516-524
• /
• 2008

This paper presents a new modeling technique for the simulation of cloth specific characteristics with a set of hard constraints using an implicit constraint enforcement scheme. A conventional explicit Baumgarte constraint stabilization method has several defects. It requires users to pick problem-dependent coefficients to achieve fast convergence and has inherent stabilization limits. The proposed implicit constraint enforcement method is stable with large time steps, does not require problem dependent feed-back parameters, and guarantees the natural physics-based motion of an object. In addition, its computational complexity is the same as the explicit Baumgarte method. This paper describes a formulation of implicit constraint enforcement and provides a constraint error analysis. The modeling technique for complex components of cloth such as seams, buttons, sharp creases, wrinkles, and prevention of excessive elongation are explained. Combined with an adaptive constraint activation scheme, the results using the proposed method show the substantial enhancement of the realism of cloth simulations with a corresponding savings in computational cost.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.4D
• /
• pp.671-678
• /
• 2006

In this study, to describe vehicle-structure dynamic interaction phenomena with 1/4 vehicle model, nonlinear Hertzian contact spring and nonlinear contact damper are adopted. The external loads acting on 1/4 vehicle model are selfweight of vehicle and geometry information of running surface. The constraint equation on contact surface is implemented by the Penalty method with stabilization and the reaction from constraint violation. To describe pitching motion of various vehicles two types of the displacement constraint equations are exerted to connect between car bodies and between bogie frames, i.e., the rigid body connection and the rigid body connection with pin, respectively. For the time integration of dynamic equations of vehicles and structure Newmark time integration scheme is adopted. To reduce the error caused by inadequate time step size, adaptive time-stepping technique is also adopted. Thus, it is expected that more versatile dynamic interaction phenomena can be described by this approach and it can be applied to various railway dynamic problems with low computational cost.

• Journal of Mechanical Science and Technology
• /
• v.19 no.spc1
• /
• pp.320-327
• /
• 2005

A new mechanism for hyper redundant manipulator (HRM) is presented, which comprises of serially assembled compound three-bar linkages (CTL). The CTL mechanism has some unique properties. This paper presents the forward and inverse kinematics of this mechanism and shows the simulation of the HRM havig 9 CTL units. The recursive algorithm of the inverse kinematics that the author originally developed is employed. It is fast and stable ; moreover, it enables us to obtain a solution in which the end-point of the HRM is controlled by a portion of joints. It also presents the method of the dynamical analysis. There exist kinematical constraints in the proposed closed linkage mechanism. In the dynamic analysis constraints are sufficiently sustained by the constraint stabilization method that the author developed. The mechanical structure of the HRM having some CTL units that is under construction is shown.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.91-98
• /
• 2004

In two dimensional incompressible flows, a preconditioning technique called Hierarchical Iterative Procedure(HIP) has been implemented on a stabilized finite element formulation. The stabilization has been peformed by a modified residual method proposed by Illinca et. al.[3]. The stabilization which is necessary to escape from the LBB constraint renders an equal order formulation. In this paper, we increased the order of interpolation whithin an element up to cubic. The conjugate gradient squared(CGS) method is used for the outer iteration, and the HIP for the preconditioning for the incompressible Navier-Stokes equation. The hierarchical elements has been used to achieve a higher order accuracy in fluid flow analyses, but a proper efficient iterative procedure for higher order finite element formulation has not been available so far. The numerical results by the present HIP for the lid driven cavity flow showed the present procedure to be stable, very efficient and useful in flow analyses in conjunction with hierarchical elements.

• Journal of the Korean Society of Industry Convergence
• /
• v.4 no.1
• /
• pp.87-94
• /
• 2001

In multibody dynamics, differential and algebraic equations which can satisfy both equation of motion and kinematic constraint equation should be solved. To solve this equation, coordinate partitioning method and constraint stabilization method are commonly used. The coordinate partitioning method divides the coordinate into independent and dependent coordinates. The most typical coordinate partitioning method arc LU decomposition, QR decomposition, projection method and SVD(sigular value decomposition).The objective of this research is to find a efficient coordinate partitioning method in flexible multibody systems and a hybrid decomposition algorithm which employs both LU and projection methods is proposed. The accuracy of the solution algorithm is checked with a slider-crank mechanism.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.21 no.8
• /
• pp.1311-1321
• /
• 1997

In multibody dynamics, differential and algebraic equations which can satisfy both equation of motion and kinematic constraint equation should be solved. To solve these equations, coordinate partitioning method and constraint stabilization method are commonly used. In the coordinate partitioning method, the coordinates are divided into independent and dependent and coordinates. The most typical coordinate partitioning method are LU decomposition, QR decomposition, and SVD (singular value decomposition). The objective of this research is to find an efficient coordinate partitioning method in the dynamic analysis of flexible multibody systems. Comparing two coordinate partitioning methods, i.e. LU and QR decomposition in the flexible multibody systems, a new hybrid coordinate partitioning method is suggested for the flexible multibody analysis.

• Journal of the Korean Society for Railway
• /
• v.5 no.3
• /
• pp.174-180
• /
• 2002

In multibody dynamics, DAE(Differential Algebraic Equations) that combine differential equations of motion and kinematic constraint equations should be solved. To solve these equations, either coordinate partitioning method or constraint stabilization method is commonly used. The most typical coordinate partitioning methods are LU decomposition, QR decomposition, and SVD(singular value decomposition). The objective of this research is to suggest a hybrid coordinate partitioning method in the dynamic analysis of multibody systems containing singular configurations. Two coordinate partitioning methods, i.e. LU decomposition and QR decomposition for constrained multibody systems, are combined for a new hybrid coordinate partitioning method. The proposed hybrid method reduces the simulation time while keeping accuracy of the solution.

• Journal of computational fluids engineering
• /
• v.9 no.3
• /
• pp.57-65
• /
• 2004

The incompressible Navier-Stokes equations in two dimensions are stabilized by a modified residual method, and then discretized by hierarchical elements. The stabilization is necessary to escape from the Ladyzhenskaya-Babuska-Brezzi(LBB) constraint and hence to achieve an equal order formulation. To expedite a standard iterative method such as the conjugate gradient squared(CGS) method, a preconditioning technique called the Hierarchical Iterative Procedure(HIP) has been applied. In this paper, we increased the order of interpolation within an element up to cubic. The hierarchical elements have been used to achieve a higher order accuracy in fluid flow analyses, but a proper efficient iterative procedure for higher order finite element formulation has not been available so far The numerical results by the present HIP for the lid driven cavity flow and others showed the present procedure to be stable, very efficient and useful in flow analyses in conjunction with hierarchical elements.