• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.893-898
• /
• 2004

This research proposes an implementation method of linearized equations of motion for multibody systems with closed loops. The null space of the constraint Jacobian is first pre multiplied to the equations of motion to eliminate the Lagrange multiplier and the equations of motion are reduced down to a minimum set of ordinary differential equations. The resulting differential equations are functions of all relative coordinates, velocities, and accelerations. Since the coordinates, velocities, and accelerations are tightly coupled by the position, velocity, and acceleration level constraints, direct substitution of the relationships among these variables yields very complicated equations to be implemented. As a consequence, the reduced equations of motion are perturbed with respect to the variations of all coordinates, velocities, and accelerations, which are coupled by the constraints. The position, velocity and acceleration level constraints are also perturbed to obtain the relationships between the variations of all relative coordinates, velocities, and accelerations and variations of the independent ones. The perturbed constraint equations are then simultaneously solved for variations of all coordinates, velocities, and accelerations only in terms of the variations of the independent coordinates, velocities, and accelerations. Finally, the relationships between the variations of all coordinates, velocities, accelerations and these of the independent ones are substituted into the variational equations of motion to obtain the linearized equations of motion only in terms of the independent coordinate, velocity, and acceleration variations.

• 제어로봇시스템학회:학술대회논문집
• /
• 2000.10a
• /
• pp.265-265
• /
• 2000

In this paper the systematic modeling method of general wheeled mobile robot is proposed. First we show how to describe kinematics properties of wheeled mobile robot in the method formulating constraint equations using Basic Homogeneous Transform(BHT) which is used mainly the kinematics modeling of manipulator, and, under assumption it's provided part of nullvector in given constraint equations, find kinematics model of mobile robot related to actuators in real robot.

• Journal of the Korean Society for Precision Engineering
• /
• v.16 no.4 s.97
• /
• pp.26-36
• /
• 1999

A general procedure for determining the optimum location of suspension hard points with respect to kinematic design parametes is presented. Suspensions are modeled as connection of rigid bodies by ideal kinematic joints. Constraint equations of the kinematic joints are expressed in terms of the generalized coordinates and hard points. By directly differentiating the constraint equations with respect to the hard points, kinematic sencitivity equations are obtained. In order to cope with algebraic complexity associated with the differentiation process, a symbolic computation technique is used. A performance index is defined in terms of static design parameters such as camber, caster, toe, ect.. Gradient of the performance index can be analytically computed from the kinematic sensitivity equations. Optimization results show the effectiveness and validity of the procedure, which is applicable to any type of suspension if its kinematic configurations are given.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.5 no.1
• /
• pp.38-48
• /
• 1997

A method for performing kinematic design sensitivity analysis of vehicle suspension systems is presented. For modeling of vehicle suspensions, the multibody dynamic formulation is adopted, where suspensions are assumed as combination of rigid bodies and ideal frictionless joints. In a relative joint coordinate setting, kinematic constraint equations are obtained by imposing cut-joints that transform closed-loop shape suspension systems into open-loop systems. By directly differentiating the constraint equations with respect to kinematic design variables, such as length of bodies, notion axis, etc., sensitivity equations are derived. By solving the sensitivity equations, sensitivity of static design factors that can be used for design improvement, can be obtained. The validity and usefulness of the method are demonstrated through an example where kinematic sensitivity analysis of a MacPherson strut suspension of performed.

• 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.

• Proceedings of the KSME Conference
• /
• 2001.06b
• /
• pp.465-470
• /
• 2001

This paper presents an algorithm which seeks steady-state equilibrium positions of constrained multibody systems driven by constant generalized speeds. Since the relative coordinates are employed, the constraint equations at cut joints are incorporated into the formulation. The proposed algorithm leads to nonlinear equations that need to be solved iteratively. This algorithm should satisfy both types of conditions: the force equilibrium equations and the kinematic constraint equations. To verify the effectiveness of the proposed algorithm, two numerical examples are solved and the results are compared with those of a commercial program. This method, compared to the conventional method of using dynamic analysis, has the advantage of computational efficiency and stability.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.5 no.1
• /
• pp.154-161
• /
• 1997

This paper presents the method to eliminate the constraint reaction in the Lagrange multiplier form equation of motion by using a generalized coordinate driveder from the velocity constraint equation. This method introduces a matrix method by considering the m dimensional space spanned by the rows of the constraint jacobian matrix. The orthogonal vectors defining the constraint manifold are projected to null vectors by the tangential vectors defined on the constraint manifold. Therefore the orthogonal projection matrix is defined by the tangential vectors. For correcting the generalized position coordinate, the optimization problem is formulated. And this correction process is analyzed by the quasi Newton method. Finally this method is verified through 3 dimensional vehicle model.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 1996.04a
• /
• pp.661-665
• /
• 1996

An inverse kinemetics problem of a reclaimer which digs and transports ironstones or coals in the raw yard is investigated. Because of the special features of the reclaimer of which scooping buckets are attached around the rotating drum at the end of boom, kinematic redundancy occurs in determining the joint varialbes For a given reclaiming point in space the forward kinematics yields 3 equations, however the number of involved variables in the equations are four. A plane equation approximating the surface near a reclaiming point is obtained by considering 8 adjacent points surrounding the reclaiming point. One extra equation to overcome redunduncyis further obtained from the condition that the normal vector at a reclaiming point is perpendicular to the plane. An approximate solution for a simplified problem is first discussed, Numerical solution for the oritinal nonlinear porblem with a constraint equation is also investigated. Finally a closed form solution which is not exact but sufficiently close enough is proposed by exploiting geometric constraint.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.2 no.4
• /
• pp.80-90
• /
• 1994

Basic constraint equations derived from orthogonality conditions between a pair of body-fixed vectors and a body-fixed vector or a vector between two bodies are reformulated by using relative coordinate kinematics between two adjacent reference frames. Arithmetic numbers of operations required to compute derivatives of the constraint equations are drastically reduced. A mixed formulation of relative and cartesian coordinates is developed to further simplify derivatives of the constraints. Advantages and disadvantages of the new formulation are discussed. Possible singularity problem of para llelism constraints is resolved by introducing an extra generalized coordinate. Kinematic analysis of a McPherson strut suspension system are carried out to illustrate use and efficiency of the new formulation.

• Proceedings of the KSR Conference
• /
• 2005.11a
• /
• pp.191-196
• /
• 2005

In this study, to describe vehicle-track-bridge dynamic interaction phenomena with 1/4 vehicle model, nonlinear Hertzian contact spring and nonlinear contact damper are introduced. In this approach external loads acting on 1/4 vehicle model are self weight of vehicle and geometry information of running surface. The constraint equation on contact surface is implemented by Penalty method. Also, to improve the numerical stability and to maintain accuracy of solution, the artificial damper and the reaction from constraint violation are introduced. A nonlinear time integration method, in this study, Newmark method is adopted for both equations of vehicles and structure. And to reduce the error caused by inadequate time step size, adaptive time-stepping technique is partially introduced. As the nonlinear Hertzian contact spring has no resistance to tensile force, the bouncing phenomena of wheelset can be described. 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.