• Structural Engineering and Mechanics
• /
• v.57 no.2
• /
• pp.327-355
• /
• 2016

A new analytical derivation of the elastodynamic point load solutions for an isotropic multi-layered half-space is presented by means of the precise integration method (PIM) and the approach of dual vector. The time-harmonic external load is prescribed either on the external boundary or in the interior of the solid medium. Starting with the axisymmetric governing motion equations in a cylindrical coordinate system, a second order ordinary differential matrix equation can be gained by making use of the Hankel integral transform. Employing the technique of dual vector, the second order ordinary differential matrix equation can be simplified into a first-order one. The approach of PIM is implemented to obtain the solutions of the ordinary differential matrix equation in the Hankel integral transform domain. The PIM is a highly accurate algorithm to solve sets of first-order ordinary differential equations and any desired accuracy of the dynamic point load solutions can be achieved. The numerical simulation is based on algebraic matrix operation. As a result, the computational effort is reduced to a great extent and the computation is unconditionally stable. Selected numerical trials are given to validate the accuracy and applicability of the proposed approach. More examples are discussed to portray the dependence of the load-displacement response on the isotropic parameters of the multi-layered media, the depth of external load and the frequency of excitation.

• Architectural research
• /
• v.13 no.1
• /
• pp.17-24
• /
• 2011

Isogeometric solution of Poisson equation is provided. NURBS (NonUniform B-spline Surface) is introduced to express both geometry of structure and unknown field of governing equation. The terms of stiffness matrix and load vector are consistently derived with very accurate geometric definition. The validity of the isogeometric formulation is demonstrated by using two numerical examples such as square plate and L-shape plate. From numerical results, the present solutions have a good agreement with analytical and finite element (FE) solutions with the use of a few cells in isogeometric analysis.

• Kyungpook Mathematical Journal
• /
• v.61 no.2
• /
• pp.383-393
• /
• 2021

In this paper, numerical techniques are presented for solving initial value problems of fractional differential equations with variable coefficients. The method is derived by applying a Taylor vector approximation. Moreover, the operational matrix of fractional integration of a Taylor vector is provided in order to transform the continuous equations into a system of algebraic equations. Furthermore, numerical examples demonstrate that this method is applicable and accurate.

• Structural Engineering and Mechanics
• /
• v.29 no.6
• /
• pp.673-687
• /
• 2008

In this paper, a thermo-viscoelastic problem in an infinite isotropic medium in two dimensions in the presence of a point heat source is considered. The fundamental equations of the problems of generalized thermoelasticity including heat sources in a thermo-viscoelastic media have been derived in the form of a vector matrix differential equation in the Laplace-Fourier transform domain for a two dimensional problem. These equations have been solved by the eigenvalue approach. The results have been compared to those available in the existing literature. The graphs have been drawn for different cases.

• Coupled systems mechanics
• /
• v.9 no.1
• /
• pp.77-89
• /
• 2020

This study is dedicated to the dynamic elasticity problem for a semi-strip. The semi-strip is loaded by the dynamic load at the center of its short edge. The conditions of fixing are given on the lateral sides of the semi-strip. The initial problem is reduced to one-dimensional problem with the help of Laplace's and Fourier's integral transforms. The one-dimensional boundary problem is formulated as the vector boundary problem in the transform's domain. Its solution is constructed as the superposition of the general solution for the homogeneous vector equation and the partial solution for the inhomogeneous vector equation. The matrix differential calculation is used for the deriving of the general solution. The partial solution is constructed with the help of Green's matrix-function, which is searched as the bilinear expansion. The case of steady-state oscillations is considered. The problem is reduced to the solving of the singular integral equation. The orthogonalization method is applied for the calculations. The stress state of the semi-strip is investigated for the different values of the frequency.

• Steel and Composite Structures
• /
• v.50 no.2
• /
• pp.123-133
• /
• 2024

The paper deals with the propagation of Rayleigh waves in a nonlocal thermoelastic isotropic layer which is lying over a nonlocal thermoelastic isotropic half-space under the purview of Green-Lindsay model and Eringen's nonlocal elasticity in the presence of voids. The normal mode analysis is employed to the considered equations to obtain vector matrix differential equation which is then solved by eigenvalue approach. The frequency equation of Rayleigh waves is derived and different particular cases are also deduced. The effects of voids and nonlocality on different characteristics of Rayleigh waves are presented graphically.

• Structural Engineering and Mechanics
• /
• v.42 no.1
• /
• pp.39-53
• /
• 2012

This paper presents an alternative way to derive the exact element stiffness matrix for a beam on Winkler foundation and the fixed-end force vector due to a linearly distributed load. The element flexibility matrix is derived first and forms the core of the exact element stiffness matrix. The governing differential compatibility of the problem is derived using the virtual force principle and solved to obtain the exact moment interpolation functions. The matrix virtual force equation is employed to obtain the exact element flexibility matrix using the exact moment interpolation functions. The so-called "natural" element stiffness matrix is obtained by inverting the exact element flexibility matrix. Two numerical examples are used to verify the accuracy and the efficiency of the natural beam element on Winkler foundation.

• 제어로봇시스템학회:학술대회논문집
• /
• 1986.10a
• /
• pp.304-309
• /
• 1986

This paper is about the Inverse Jacobian for the TOP-1 robot. The robot Jacobian is used for the movement in accordance with differental changes. A Matrix and Homogeneous Transformation Matrix, Differential Motion Vector D are applied to Jacobian equation for the movement of the robot in accordance with the minut changes. The solution of Jacobian equation is acquired and applied for the subtle movement of each arms of the robot. The interface with APPLE-II Micro-computer is searched out too. The Software and the interface resulted from this paper are considered to be higly useful in the accurate control on the robot when they are linked with dynamics of robot.

• Journal of Mechanical Science and Technology
• /
• v.16 no.8
• /
• pp.1072-1078
• /
• 2002

The constrained motion requires the determination of constraint force acting on unconstrained systems for satisfying given constraints. Most of the methods to decide the force depend on numerical approaches such that the Lagrange multiplier method, and the other methods need vector analysis or complicated intermediate process. In 1992, Udwadia and Kalaba presented the generalized inverse method to describe the constrained motion as well as to calculate the constraint force. The generalized inverse method has the advantages which do not require any linearization process for the control of nonlinear systems and can explicitly describe the motion of holonomically and/or nongolonomically constrained systems. In this paper, an explicit equation to describe the constrained motion is derived by minimizing the performance index, which is a function of constraint force vector, with respect to the constraint force. At this time, it is shown that the positive-definite weighting matrix in the performance index must be the inverse of mass matrix on the basis of the Gauss's principle and the derived differential equation coincides with the generalized inverse method. The effectiveness of this method is illustrated by means of two numerical applications.

• Proceedings of the KSME Conference
• /
• 2000.11a
• /
• pp.529-534
• /
• 2000

A finite element analysis for a rotating cantilever beam is presented in this study. Based on a dynamic modelling method using the stretch deformation instead of the conventional axial deformation, three linear partial differential equations are derived from Hamilton's principle. Two of the linear differential equations show the coupling effect between stretch and chordwise deformations. The other equation is an uncoupled one for the flapwise deformation. From these partial differential equations and the associated boundary conditions, are derived two weak forms: one is for the chordwise motion and the other is for the flapwise motion. The weak forms are spatially discretized with newly defined two-node beam elements. With the discretized equations or the matrix-vector equations, the behaviours of the natural frequencies are investigated for the variation of the rotating speed.