• Coupled systems mechanics
• /
• v.12 no.6
• /
• pp.519-530
• /
• 2023

This paper presents a mathematical model suitable for the calculation of laminated glass, i.e. glass plates combined with an interlayer material. The model is based on a beam differential equation for each glass plate and a separate differential equation for the slip in the interlayer. In addition to slip, the model takes into account prestressing force in the interlayer. It is possible to combine the two contributions arbitrarily, which is important because the glass sheet fabrication process changes the stiffness of the interlayer in ways that are not easily predictable and could introduce prestressing of varying magnitude. The model is suitable for reformulation into an inverse procedure for calculation of the relevant parameters. Model consisting of a system of differential-algebraic equations, proved too stiff for cases with the thin interlayer. This novel approach covers the full range of possible stiffnesses of layered glass sheets, i.e., from zero to infinite stiffness of the interlayer. The comparison of numerical and experimental results contributes to the validation of the model.

• 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 KIEE Conference
• /
• 1995.07b
• /
• pp.667-669
• /
• 1995

This paper presents a method to design Kalman filter on continuous stochastic dynamical systems via BPFT(block pulse functions transformation). When we design Kalman filter, minimum error valiance matrix is appeared as a form of nonlinear matrix differential equations. Such equations are very difficult to obtain the solutions. Therefore, in this paper, we simply obtain the solutions of nonlinear matrix differential equations from recursive algebraic equations using BPFT. We believe that the presented method is very attractive and proper for the evaluation of Kalman gain on continuous stochastic dynamical systems.

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

• 제어로봇시스템학회:학술대회논문집
• /
• 1992.10a
• /
• pp.680-685
• /
• 1992

Solar arrays and antennas of the satellite are usually stowed within the dimensions of the launch-vehicle fairing and deployed in the orbit. To solve such multibody dynamic problems, differential equations and algebraic equations are simultaneously solved, and special solution techniques are required. In this paper, Lagrange multipliers associated with the constraints are iteratively computed by monotonically reducing an appropriately defined constraint error vector, and the resulting equation of motion is solved by a well-established ODE technique. Defomable bodies as well as rigid bodies are treated, and applications to satellite solar arrays are explained.

• Transactions of the Korean Society of Machine Tool Engineers
• /
• v.10 no.1
• /
• pp.112-119
• /
• 2001

Design sensitivity is an important is an important device in improving a mechanical system design. A continuum design consists of the shape and orientation design. This research develops the shape and orientation design sensitivity method. The configura-tion design variables of multibody systems define the shape and orientation changes. The equations of motion are directly differentiated to obtain the governing equations for the design sensitivity. The governing equation of the design sensitivity is formulated as an over determined differential algebraic equation and treated as ordinary differential equations on mani-folds. The material derivative of a domain functional is performed to obtain the sensitivity due to shape and orientation changes. The configuration design sensitivities of a fly-ball governor system and a spatial four bar mechanism are obtained using the proposed method and are validated against those obtained from the finite difference method.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1996.10a
• /
• pp.189-196
• /
• 1996

Vibration of a beam translating over supports in vertical motion is investigated in this paper. Equations of motion are formulated using the virtual work principle by regarding the supports as kinematical constraints imposed on an unrestrained beam and by discretizing the beam via the assumed mode method. Differential-algebraic equations of motion are derived and reduced to differential equations in independent generalized coordinates by the generalized coordinate partitioning method. Geometric stiffness of the beam due to translating motion is considered and how the geometric stiffness of beam affects dynamic stability is also investigated. Instability of the beam. in various conditions is also investigated using Floquet theory and then the results are verified through the dynamic response analysis. Results of numerical simulation are presented for various prescribed motions of the beam.

• Journal of the Korean Society for Precision Engineering
• /
• v.17 no.11
• /
• pp.165-175
• /
• 2000

In recent years, it becomes a very important issue to consider the mechanical systems such as high-speed vehicles and railway trains moving on elastic beam structures. In this paper, a general approach, which can predict the dynamic behavior of constrained mechanical system and elastic beam structure, is proposed. Also, various supporting conditions of a foundation support are considered for the elastic beam structures. The elastic structure is assumed to be a nonuniform and linear Bernoulli-Euler beam with proportional damping effect. Combined Differential-Algebraic Equations of motion are derived using multibody dynamics theory and Finite Element Method. The proposed equations of motion can be solved numerically using generalizd coordinate partitioning method and Predictor-Corrector algorithm, which is an implicit multi-step integration method.

• Steel and Composite Structures
• /
• v.26 no.4
• /
• pp.473-484
• /
• 2018

Free vibration of cross-ply laminated plates using a higher-order shear deformation theory is studied. The arbitrary number of layers is oriented in symmetric and anti-symmetric manners. The plate kinematics are based on higher-order shear deformation theory (HSDT) and the vibrational behaviour of multi-layered plates are analysed under simply supported boundary conditions. The differential equations are obtained in terms of displacement and rotational functions by substituting the stress-strain relations and strain-displacement relations in the governing equations and separable method is adopted for these functions to get a set of ordinary differential equations in term of single variable, which are coupled. These displacement and rotational functions are approximated using cubic and quantic splines which results in to the system of algebraic equations with unknown spline coefficients. Incurring the boundary conditions with the algebraic equations, a generalized eigen value problem is obtained. This eigen value problem is solved numerically to find the eigen frequency parameter and associated eigenvectors which are the spline coefficients.The material properties of Kevlar-49/epoxy, Graphite/Epoxy and E-glass epoxy are used to show the parametric effects of the plates aspect ratio, side-to-thickness ratio, stacking sequence, number of lamina and ply orientations on the frequency parameter of the plate. The current results are verified with those results obtained in the previous work and the new results are presented in tables and graphs.

• Coupled systems mechanics
• /
• v.1 no.2
• /
• pp.155-163
• /
• 2012

This paper presents analytical approximate solutions for the initial post-buckling deformation of the pipes in oil and gas wells. The governing differential equation with sinusoidal nonlinearity can be reduced to form a third-order-polynomial nonlinear equation, by coupling of the well-known Maclaurin series expansion and orthogonal Chebyshev polynomials. Analytical approximations to the resulting boundary condition problem are established by combining the Newton's method with the method of harmonic balance. The linearization is performed prior to proceeding with harmonic balancing thus resulting in a set of linear algebraic equations instead of one of non-linear algebraic equations, unlike the classical method of harmonic balance. We are hence able to establish analytical approximate solutions. The approximate formulae for load along axis, and periodic solution are established for derivative of the helix angle at the end of the pipe. Illustrative examples are selected and compared to "reference" solution obtained by the shooting method to substantiate the accuracy and correctness of the approximate analytical approach.