• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.24 no.3
• /
• pp.243-291
• /
• 2020

In this paper, multi-block generalized backward differentiation methods for numerical solutions of ordinary differential and differential algebraic equations are introduced. This class of linear multi-block methods is implemented as multi-block boundary value methods (MB₂ VMs). The root distribution of the stability polynomial of the new class of methods are determined using the Wiener-Hopf factorization of a matrix polynomial for the purpose of their correct implementation. Numerical tests, showing the potential of such methods for output of multi-block of solutions of the ordinary differential equations in the new approach are also reported herein. The methods which output multi-block of solutions of the ordinary differential equations on application, are unlike the conventional linear multistep methods which output a solution at a point or the conventional boundary value methods and multi-block methods which output only a block of solutions per step. The MB₂ VMs introduced herein is a novel approach at developing very large scale integration methods (VLSIM) in the numerical solution of differential equations.

• Journal of applied mathematics & informatics
• /
• v.9 no.2
• /
• pp.561-577
• /
• 2002

In this paper, a two-dimensional finite volume unstructured mesh method (FVUM) based on a triangular background interpolation mesh is developed for analysing the evolution of the saltwater intrusion into single and multiple coastal aquifer systems. The model formulation consists of a ground-water flow equation and a salt transport equation. These coupled and non-linear partial differential equations are transformed by FVUM into a system of differential/algebraic equations, which is solved using backward differentiation formulas of order one through five. Simulation results are compared with previously published solutions where good agreement is observed.

• Coupled systems mechanics
• /
• v.3 no.3
• /
• pp.291-304
• /
• 2014

Paper discusess a dynamic engineering problem of a mass attached to a pendulum sliding along a cable. In this problem the pendulum mass and the cable are coupled together in a model described by a system of differential algebraic equations (DAE). In the paper we have presented formulation of the system of differential equations that models the problem and determination of the initial conditions. The developed model is general in a sense of free choice of support location, elastic cable properties, pendulum length and inclusion of braking forces. Examples illustrate and validate the model.

• Journal of the Korean Institute of Telematics and Electronics
• /
• v.16 no.5
• /
• pp.18-26
• /
• 1979

This paper extends some well-known system theories on algebraic matrix Lyapunov and Riccati equations. These extended results contain two point boundary conditions in matrix differential equations and include conventional results as special cases. Necessary and sufficient conditions are derived under which linear systems are stabilizable with feedback gains derived from periodic two-point boundary matrix differential equations. An iterative computation method for two-point boundary differential Riccati equations is given with an initial guess method. The results in this paper are related to periodic feedback controls and also to the quadratic cost problem with a discrete state penalty.

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

A formulation for flexible multibody systems (MBS) is investigated, where rigid MBS substructures are coupled with flexible bodies described by a nonlinear finite element (FE) approach. Several aspects that turned out to be crucial for the presented approach are discussed. The system describing equations are given in differential algebraic form (DAE), where many sophisticated solvers exist. In this paper the performance of several solvers is investigated regarding their suitability for the application to the usually highly stiff DAE. The substructures are connected with each other by nonlinear algebraic constraint equations. Further, partial derivatives of the constraints are required, which often leads to extensive algebraic trans-formations. Handcoding of analytically determined derivatives is compared to an approach utilizing algorithmic differentiation.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 1992.10a
• /
• pp.349-353
• /
• 1992

The purpose of this paper is to develop methods for the dynamic analysis of multibody system that consist of interconnected rigid and deformable component. The equations of motion are derived by using the Lagrange's equation and finite element theory for the elastic mechanism systems. The type of equation of motion is the differential algebraic equation included kinematic nonlinear algebraic equation. The generalized coordinate partitioning method is used for solving this equation. To show the validity of this analysis solver, couple of models were canalized and those results were compared with the commercial package(ADAMS).

• International Journal of Precision Engineering and Manufacturing
• /
• v.3 no.4
• /
• pp.54-63
• /
• 2002

In recent years, mechanical systems such as high speed vehicles and railway trains moving on elastic beam structures have become a very important issue to consider. In this paper, a general approach, which can predict the dynamic behavior of a constrained mechanical system moving on a flexible beam structure, is proposed. Various supporting conditions for the foundation support are considered for the elastic beam structure. The elastic structure is assumed to be a non-uniform and linear Bernoulli-Euler beam with a proportional damping effect. Combined differential-algebraic equation of motion is derived using the multi-body dynamics theory and the finite element method. The proposed equations of motion can be solved numerically using the generalized coordinate partitioning method and predictor-corrector algorithm, which is an implicit multi-step integration method.

• International Journal of Precision Engineering and Manufacturing
• /
• v.3 no.4
• /
• pp.64-71
• /
• 2002

Recently, mechanical systems such as a high-speed vehicles and railway trains moving on flexible beam structures have become a very important issue to consider. Using the general approach proposed in the first part of this paper, it is possible to predict motion of the constrained mechanical system and the elastic structure, with various kinds of foundation supporting conditions. Combined differential-algebraic equation of motion derived from both multibody dynamics theory and finite element method can be analyzed numerically using a generalized coordinate partitioning algorithm. To verify the validity of this approach, results from the simply supported elastic beam subjected to a moving load are compared with the exact solution from a reference. Finally, parametric study is conducted for a moving vehicle model on a simply supported 3-span bridge.

• Structural Engineering and Mechanics
• /
• v.5 no.3
• /
• pp.283-295
• /
• 1997

This paper presents a high precision integration method for the dynamic response analysis of structures with holonomic constraints. A detail recursive scheme suitable for algebraic and differential equations (ADEs) which incorporates generalized forces is established. The matrix exponential involved in the scheme is calculated precisely using $2^N$ algorithm. The Taylor expansions of the nonlinear term concerned with state variables of the structure and the generalized constraint forces of the ADEs are derived and consequently, their particular integrals are obtained. The accuracy and effectiveness of the present method is demonstrated by two numerical examples, a plane truss with circular slot at its tip point and a slewing flexible cantilever beam which is currently interesting in optimal control of robot manipulators.

• The Transactions of the Korean Institute of Electrical Engineers D
• /
• v.52 no.4
• /
• pp.198-204
• /
• 2003

In this paper, we presented a new algebraic iterative algorithm for the optimal control of the nonlinear systems. The algorithm is based on two steps. The first step transforms nonlinear optimal control problem into a sequence of linear optimal control problem using the quasilinearization method. In the second step, TPBCP(two point boundary condition problem) is solved by algebraic equations instead of differential equations using the new integral operational matrix of BPF(block pulse functions). The proposed algorithm is simple and efficient in computation for the optimal control of nonlinear systems and is less error value than that by the conventional matrix. In computer simulation, the algorithm was verified through the optimal control design of synchronous machine connected to an infinite bus.