• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.25 no.4
• /
• pp.262-295
• /
• 2021

This paper introduces mono-implicit general linear methods, a special class of general linear methods, which are implicit in the output solution for the numerical integration of differential algebraic equations. We show how L-stable inherent Runge-Kutta members can be derived. The procedures for implementation have been discussed. The numerical test on the problem considered shows that the methods have improved accuracy when compared to RADAU IIA and the results from MATLAB ode15s, which have been taken as our reference solution.

• Structural Engineering and Mechanics
• /
• v.28 no.2
• /
• pp.221-238
• /
• 2008

Numerical solution to buckling analysis of beams and columns are obtained by the method of differential quadrature (DQ) and harmonic differential quadrature (HDQ) for various support conditions considering the variation of flexural rigidity. The solution technique is applied to find the buckling load of fully or partially embedded columns such as piles. A simple semi- inverse method of DQ or HDQ is proposed for determining the flexural rigidities at various sections of non-prismatic column ( pile) partially and fully embedded given the buckling load, buckled shape and sub-grade reaction of the soil. The obtained results are compared with the existing solutions available from other numerical methods and analytical results. In addition, this paper also uses a recently developed technique, known as the differential transformation (DT) to determine the critical buckling load of fully or partially supported heavy prismatic piles as well as fully supported non-prismatic piles. In solving the problem, governing differential equation is converted to algebraic equations using differential transformation methods (DT) which must be solved together with applied boundary conditions. The symbolic programming package, Mathematica is ideally suitable to solve such recursive equations by considering fairly large number of terms.

• 제어로봇시스템학회:학술대회논문집
• /
• 1993.10a
• /
• pp.926-930
• /
• 1993

An efficient method of formulating the equations of motion of multibody systems is presented. The equations of motion for each body are formulated by using Newton-Eulerian approach in their generic form. And then a transformation matrix which relates the global coordinates and relative coordinates is introduced to rewrite the equations of motion in terms of relative coordinates. When appropriate set of kinematic constraints equations in terms of relative coordinates is provided, the resulting differential and algebraic equations are obtained in a suitable form for computer implementation. The system geometry or topology is effectively described by using the path matrix and reference body operator.

• Steel and Composite Structures
• /
• v.35 no.5
• /
• pp.671-685
• /
• 2020

This manuscript presents impacts of gradation of material functions and axial load functions on critical buckling loads and mode shapes of functionally graded (FG) thin and thick beams by using higher order shear deformation theory, for the first time. Volume fractions of metal and ceramic materials are assumed to be distributed through a beam thickness by both sigmoid law and symmetric power functions. Ceramic-metal-ceramic (CMC) and metal-ceramic-metal (MCM) symmetric distributions are proposed relative to mid-plane of the beam structure. The axial compressive load is depicted by constant, linear, and parabolic continuous functions through the axial direction. The equilibrium governing equations are derived by using Hamilton's principles. Numerical differential quadrature method (DQM) is developed to discretize the spatial domain and covert the governing variable coefficients differential equations and boundary conditions to system of algebraic equations. Algebraic equations are formed as a generalized matrix eigenvalue problem, that will be solved to get eigenvalues (buckling loads) and eigenvectors (mode shapes). The proposed model is verified with respectable published work. Numerical results depict influences of gradation function, gradation parameter, axial load function, slenderness ratio and boundary conditions on critical buckling loads and mode-shapes of FG beam structure. It is found that gradation types have different effects on the critical buckling. The proposed model can be effective in analysis and design of structure beam element subject to distributed axial compressive load, such as, spacecraft, nuclear structure, and naval structure.

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

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

• Structural Engineering and Mechanics
• /
• v.43 no.1
• /
• pp.105-126
• /
• 2012

Based on the first-order shear deformation theory (FSDT), and the virtual work principle, an elastic analysis for axisymmetric clamped-clamped Pressurized thick truncated conical shells made of functionally graded materials have been performed. The governing equations are a system of nonhomogeneous ordinary differential equations with variable coefficients. Using the matched asymptotic method (MAM) of the perturbation theory, these equations could be converted into a system of algebraic equations with variable coefficients and two systems of differential equations with constant coefficients. For different FGM conical angles, displacements and stresses along the radius and length have been calculated and plotted.

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

• The Transactions of The Korean Institute of Electrical Engineers
• /
• v.56 no.6
• /
• pp.1000-1006
• /
• 2007

This paper presents a methodology to calculate an optimal solution of equilibrium to differential algebraic equations for power systems. It employs a nonlinear interior point method to solve the optimization formulation which includes dynamic equations representing the two-axis synchronous generator model with AVR and speed governing controls, algebraic equations, and steady-state nonlinear loads. This paper also adopts two algorithms for the improvement of solution convergence. In power system analysis and control, equilibrium optimization (EOPT) is applicable for diverse purposes that need the consideration of dynamic model characteristics at a steady-state condition.

• Journal of Institute of Control, Robotics and Systems
• /
• v.14 no.12
• /
• pp.1270-1277
• /
• 2008

A Haar wavelet based numerical method for solving singularly perturbed linear time invariant system is presented in this paper. The reduced pure slow and pure fast subsystems are obtained by decoupling the singularly perturbed system and differential matrix equations are converted into algebraic Sylvester matrix equations via Haar wavelet technique. The operational matrix of integration and its inverse matrix are utilized to reduce the computational time to the solution of algebraic matrix equations. Finally a numerical example is given to demonstrate the validity and applicability of the proposed method.