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

In this paper, Chaotic motions of a curved pipe conveying oscillatory flow are theoretically investigated. The nonlinear partial differential equation of motion is derived by Newton's method. The transformed nonlinear ordinary differential equation is a type of Hill's equation, which have the parametric and external excitation. Bifurcation curves of chaotic motion of the piping systems are obtained by applying Melnikov's method. Poincare maps numerically demonstrate theoretical results and show transverse homoclinic orbit of the chaotic motion.

• Structural Engineering and Mechanics
• /
• v.39 no.1
• /
• pp.21-46
• /
• 2011

The natural frequencies of continuous systems depend on the governing partial differential equation and can be numerically estimated using the finite element method. The accuracy and convergence of the finite element method depends on the choice of basis functions. A basis function will generally perform better if it is closely linked to the problem physics. The stiffness matrix is the same for either static or dynamic loading, hence the basis function can be chosen such that it satisfies the static part of the governing differential equation. However, in the case of a rotating beam, an exact closed form solution for the static part of the governing differential equation is not known. In this paper, we try to find an approximate solution for the static part of the governing differential equation for an uniform rotating beam. The error resulting from the approximation is minimized to generate relations between the constants assumed in the solution. This new function is used as a basis function which gives rise to shape functions which depend on position of the element in the beam, material, geometric properties and rotational speed of the beam. The results of finite element analysis with the new basis functions are verified with published literature for uniform and tapered rotating beams under different boundary conditions. Numerical results clearly show the advantage of the current approach at high rotation speeds with a reduction of 10 to 33% in the degrees of freedom required for convergence of the first five modes to four decimal places for an uniform rotating cantilever beam.

• Journal of applied mathematics & informatics
• /
• v.27 no.1_2
• /
• pp.441-452
• /
• 2009

In this paper, a fifth order numerical method is presented for solving singularly perturbed differential-difference equations with negative shift. In recent papers the term negative shift has been using for delay. Similar boundary value problems are associated with expected first exit time problem of the membrane, potential in models for neuron and in variational problems in control theory. In the numerical treatment for such type of boundary value problems, first we use Taylor approximation to tackle terms containing small shifts which converts it to a boundary value problem for singularly perturbed differential equation. The two point boundary value problem is transformed into general first order ordinary differential equation system. A discrete approximation of a fifth order compact difference scheme is presented for the first order system and is solved using the boundary conditions. Several numerical examples are solved and compared with exact solution. It is observed that present method approximates the exact solution very well.

• Journal of applied mathematics & informatics
• /
• v.32 no.1_2
• /
• pp.39-53
• /
• 2014

In this paper, a non-asymptotic method is presented for solving singularly perturbed delay differential equations whose solution exhibits a boundary layer behavior. The second order singularly perturbed delay differential equation is replaced by an asymptotically equivalent first order neutral type delay differential equation. Then, Simpson's integration formula and linear interpolation are employed to get three term recurrence relation which is solved easily by Discrete Invariant Imbedding Algorithm. Some numerical examples are given to validate the computational efficiency of the proposed numerical scheme for various values of the delay and perturbation parameters.

• Journal of applied mathematics & informatics
• /
• v.30 no.5_6
• /
• pp.981-992
• /
• 2012

Integro-differential equations arise in modeling various physical and engineering problems. Several numerical and analytical methods have been developed to solving such equations. We introduce the NHPM for solving nonlinear integro-differential equations. Several examples for solving integro-differential equations are presented to illustrate the efficiency of the proposed NHPM.

• Journal of Korean Society of Steel Construction
• /
• v.15 no.3
• /
• pp.313-320
• /
• 2003

The governing equation of the post-buckling of shear-deformable uniform columns under a combined load consisting of a uniformly distributed axial load and a concentrated load at a free end was derived and the post-buckling analysis was investigated by using differential transformation. The loads were obtained for various end-slopes. The results obtained by the present method agree well with published results. In this paper, the differential transformation method was illustrated through its application to the non-linear differential equation of the post-buckling. It is expected that applications of the method to more challenging problems will are expected follow in future to ensue.

• Nonlinear Functional Analysis and Applications
• /
• v.28 no.2
• /
• pp.449-471
• /
• 2023

In this paper, we study a linear system of homogeneous commensurate /incommensurate fractional-order differential equations by developing a new semi-analytical scheme. In particular, by decoupling the system into two fractional-order differential equations, so that the first equation of order (δ + γ), while the second equation depends on the solution for the first equation, we have solved the under consideration system, where 0 < δ, γ ≤ 1. With the help of using the Adomian decomposition method (ADM), we obtain the general solution. The efficiency of this method is verified by solving several numerical examples.

• Structural Engineering and Mechanics
• /
• v.17 no.1
• /
• pp.1-14
• /
• 2004

Numerical solution to linear bending analysis of circular plates is obtained by the method of harmonic differential quadrature (HDQ). In the method of differential quadrature (DQ), partial space derivatives of a function appearing in a differential equation are approximated by means of a polynomial expressed as the weighted linear sum of the function values at a preselected grid of discrete points. The method of HDQ that was used in the paper proposes a very simple algebraic formula to determine the weighting coefficients required by differential quadrature approximation without restricting the choice of mesh grids. Applying this concept to the governing differential equation of circular plate gives a set of linear simultaneous equations. Bending moments, stresses values in radial and tangential directions and vertical deflections are found for two different types of load. In the present study, the axisymmetric bending behavior is considered. Both the clamped and the simply supported edges are considered as boundary conditions. The obtained results are compared with existing solutions available from analytical and other numerical results such as finite elements and finite differences methods. A comparison between the HDQ results and the finite difference solutions for one example plate problem is also made. The method presented gives accurate results and is computationally efficient.

• Computational Structural Engineering
• /
• v.6 no.4
• /
• pp.83-88
• /
• 1993

A time-discontinuous Galerkin method based upon using a finite element formulation in time has evolved. This method, working from the differential equation viewpoint, is different from those which have been generally used. They admit discontinuities with respect to the time variable at each time step. In particular, the elements can be chosen arbitrarily at each time step with no connection with the elements corresponding to the previous step. Interpolation functions and weighting functions are taken to be discontinuous across inter-element boundaries. These methods lead to a unconditional stable higher-order accurate ordinary differential equation solver.

• Journal of Mechanical Science and Technology
• /
• v.14 no.2
• /
• pp.168-176
• /
• 2000

Related to the error dynamics of an adaptive system, averaging theorems are developed for coupled differential equations which consist of ordinary differential equations and a parabolic partial differential equation. The results are then applied to the convergence analysis of the parameter estimate errors in the model reference adaptive control of a nonautonomous parabolic partial differential equation with lowly time-varying parameters.