• Journal of KSNVE
• /
• v.9 no.1
• /
• pp.196-205
• /
• 1999

The existence. bifurcation. and the orbital stability of periodic motions, which is called nonlinear normal mode, in a nonlinear dual mass Hamiltonian system. which has 6th order homogeneous polynomial as a nonlinear term. are studied in this paper. By direct integration of the equations of motion. Poincare Map. which is a mapping of a phase trajectory onto 2 dimensional surface in 4 dimensional phase space. is obtained. And via the Birkhoff-Gustavson canonical transformation, the analytic expression of the invariant curves in the Poincare Map is derived for small value of energy. It is found that the nonlinear system. which is considered in this paper. has 2 or 4 nonlinear normal modes depending on the value of nonlinear parameter. The Poincare Map clearly shows that the bifurcation modes are stable while the mode from which they bifurcated out changes from stable to unstable.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.11 no.3
• /
• pp.419-428
• /
• 2000

Multiconductor transmission line(MTL) equations are solved by FDTD(Finite-Difference Time-Domain) method to predict crosstalk and fields to transmission line coupling on lossy multiconductor transmission lines terminated in arbitrary linear loads. Skin effect losses as well as dc losses are included in the analysis. In order to increase computational efficiency, the convolution integral of internal impedance of conductors and the line currents is computed by using Prony method. For boundary conditions of MTLs terminated in linear loads, state-variable formulation is adopted. The simulated results by FDTD method are compared with the measured ones obtained by using TEM cell. The predictions are in good agreement with the measurements. In addition, it has been found that skin effect losses as well as dc losses of the conductors should be included for accurate predictions on relatively high loss transmission lines such as PCB. It has also been found that dc losses and skin-effect losses affect late-time responses and early-time responses, respectively.

• Steel and Composite Structures
• /
• v.26 no.6
• /
• pp.673-691
• /
• 2018

The main goal of this research is to examine the in-plane and out-of-plane forced vibration of a curved nanocomposite microbeam. The in-plane and out-of-plane displacements of the structure are considered based on the first order shear deformation theory (FSDT). The curved microbeam is reinforced by functionally graded carbon nanotubes (FG-CNTs) and thus the extended rule of mixture is employed to estimate the effective material properties of the structure. Also, the small scale effect is captured using the strain gradient theory. The structure is rested on a nonlinear orthotropic viscoelastic foundation and is subjected to concentrated transverse harmonic external force, thermal and magnetic loads. The derivation of the governing equations is performed using energy method and Hamilton's principle. Differential quadrature (DQ) method along with integral quadrature (IQ) and Newmark methods are employed to solve the problem. The effect of various parameters such as volume fraction and distribution type of CNTs, boundary conditions, elastic foundation, temperature changes, material length scale parameters, magnetic field, central angle and width to thickness ratio are studied on the frequency and force responses of the structure. The results indicate that the highest frequency and lowest vibration amplitude belongs to FGX distribution type while the inverse condition is observed for FGO distribution type. In addition, the hardening-type response of the structure with FGX distribution type is more intense with respect to the other distribution types.

• Structural Engineering and Mechanics
• /
• v.64 no.4
• /
• pp.391-402
• /
• 2017

In this paper, a new nonlocal trigonometric shear deformation theory is proposed for thermal buckling response of nanosize functionally graded (FG) nano-plates resting on two-parameter elastic foundation under various types of thermal environments. This theory uses for the first time, undetermined integral variables and it contains only four unknowns, that is even less than the first shear deformation theory (FSDT). It is considered that the FG nano-plate is exposed to uniform, linear and sinusoidal temperature rises. Mori-Tanaka model is utilized to define the gradually variation of material properties along the plate thickness. Nonlocal elasticity theory of Eringen is employed to capture the size influences. Through the stationary potential energy the governing equations are derived for a refined nonlocal four-variable shear deformation plate theory and then solved analytically. A variety of examples is proposed to demonstrate the importance of elastic foundation parameters, various temperature fields, nonlocality, material composition, aspect and side-to-thickness ratios on critical stability temperatures of FG nano-plate.

• Earthquakes and Structures
• /
• v.15 no.4
• /
• pp.369-382
• /
• 2018

A free vibration analysis and wave propagation of functionally graded porous beams has been presented in this work using a high order hyperbolic shear deformation theory. Unlike other conventional shear deformation theories, a new displacement field that introduces indeterminate integral variables has been used to minimize the number of unknowns. The constituent materials of the beam are assumed gradually variable along the direction of height according to a simple power law distribution in terms of the volume fractions of the constituents. The variation of the pores in the direction of the thickness influences the mechanical properties. It is therefore necessary to predict the effect of porosity on vibratory behavior and wave velocity of FG beams in this study. A new function of the porosity factor has been developed. Hamilton's principle is used for the development of wave propagation equations in the functionally graded beam. The analytical dispersion relationship of the FG beam is obtained by solving an eigenvalue problem. Illustrative numerical examples are given to show the effects of volume fraction distributions, beam height, wave number, and porosity on free vibration and wave propagation in a functionally graded beam.

• Journal of the Korean Society for Advanced Composite Structures
• /
• v.4 no.1
• /
• pp.1-8
• /
• 2013

The finite element analysis (FEA) is a numerical technique to find solutions of field problems. A field problem is approximated by differential equations or integral expressions. In a finite element, the field quantity is allowed to have a simple spatial variation in terms of linear or polynomial functions. This paper represents a review and an accuracy-study of the finite element method comparing the FEA results with the exact solution. The exact solutions were calculated by solid mechanics and FEA using matrix stiffness method. For this study, simple bar and cantilever models were considered to evaluate four types of basic elements - constant strain triangle (CST), linear strain triangle (LST), bi-linear-rectangle(Q4),and quadratic-rectangle(Q8). The bar model was subjected to uniaxial loading whereas in case of the cantilever model moment loading was used. In the uniaxial loading case, all basic element results of the displacement and stress in x-direction agreed well with the exact solutions. In the moment loading case, the displacement in y-direction using LST and Q8 elements were acceptable compared to the exact solution, but CST and Q4 elements had to be improved by the mesh refinement.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.38 no.1
• /
• pp.17-23
• /
• 2014

Graphics processing units (GPUs) are ideal for solving problems involving parallel data computations. In this study, the GPU is used for effectively carrying out a multi-body dynamic simulation with particle dynamics. The Hilber-Hushes-Taylor (HHT) implicit integration algorithm is used to solve the integral equations. For detecting collisions among particles, the spatial subdivision algorithm and discrete-element methods (DEM) are employed. The developed program is verified by comparing its results with those of ADAMS. The numerical efficiencies of the serial program using the CPU and the parallel program using the GPU are compared in terms of the number of particles, and it is observed that when the number of particles is greater, more computing time is saved by using the GPU. In the present example, when the number of particles is 1,300, the computational speed of the parallel analysis program is about 5 times faster than that of the serial analysis program.

• Journal of the korean Society of Automotive Engineers
• /
• v.2 no.1
• /
• pp.21-31
• /
• 1980

Two phase boundary layer equations of laminar filmwise condensation are solved by an approximate integral method under the following condition; saturated vapour flows vertically downward over a cooled surface of uniform temperature, the condensate film is so thin that the inertia and convection terms are neglected. The following conclusions are drawn under the above assumptions. 1. free convection In case of the linear temperature profile in a liquid film, numerical results for the average coefficients of heat transfer may be expressed as N $u_{m}$=4/3,(G $r_{l}$ /4.H)$^{1}$4/ and in case of the quadratic profile, numerical results may be expressed as N $u_{m}$=2/1.682,(G $r_{l}$ /H)$^{1}$4/. 2. Forced convection When the temperature profile is assumed to be linear in a liquid film, numerical results fir the average heat transfer coefficients may be expressed as N $u_{m}$=(A, R $e_{l}$ /H)$^{1}$2/. This expression is compared with the experimental results hitherto reported; For theoretical Nusselt number (N $u_{m}$)$_{th}$<2*10$^{4}$, the experimental Nusselt number (N $u_{m}$)$_{exp}$ is on the average larger than theoretical Nusselt number (N $u_{m}$)$_{th}$ by 30%. For (N $u_{m}$)$_{th}$>2*10$^{4}$, experimental Nusselt number (N $u_{m}$)$_{exp}$ is about 1.6 times as large as theoretical Nusselt number (N $u_{m}$)$_{th}$. These large deviation may be caused by the presence of turbulence in the liquid film. In case of the quadratic temperature profile in a liquid film, numerical results for the average coefficients of heat transfer may be expressed as N $u_{m}$'=(2,A,Re/H)$^{1}$2/. This formular shows that theoretical Nusselt number (N $u_{m}$)$_{th}$ is larger than experimental Nusselt number (N $u_{m}$)$_{exp}$ by 60%. It is speculated that when the temperature difference between cooled surface and saturated vapour is small, temperature profile in a liquid film is quadratic.quadratic.. quadratic.quadratic..atic..

• Proceedings of the Korean Geotechical Society Conference
• /
• 1998.05a
• /
• pp.153-204
• /
• 1998

Prediction of lining loads due to tunnelling is one of the major issues to be addressed in the design of a tunnel. The objective of this study is to investigate rational and realistic design loads on tunnel linings. factors influencing the lining load are summarized and discussed. The instruments for measuring the lining loads are reviewed and discussed because field measurements are often necessary to verify the design methods. Tunnel construction in the City of Edmonton has been very active for storm and sanitary purposes. Since the early 1970's, the city has also been developing an underground Light Rail Transit system. The load measurements obtained from these tunnels are compared with the results from the existing design methods. However, none of the existing methods are totally satisfactory, Therefore, there is some room for improvement in the prediction of lining loads. The convergence-confinement method is reviewed and applied to a case history of a tunnel in Edmonton. The convergence curves are obtained from 2-D finite element analyses using three different material models and theoretical equations. The limitation of the convergence-confinement method is discussed by comparing these curves with the field measurements. Three-dimensional finite element analyses are performed to gain a better understanding of stress and displacement behaviour near the tunnel face. An improved design method is proposed based on the review of existing design methods and the performance of numerical analyses. A specific method or combination of two different methods is suggested for the estimation of lining loads for different conditions of tunnelling. A method to determine the stress reduction factor is described. Typical values of dimensionless load factors nD/H for tunnels in Edmonton are obtained from parametric analyses. Finally, the loads calculated using the proposed method are compared with field measurements collected from various tunnels in terms of soil types and construction methods to verify the method. The proposed method gives a reasonable approximation of the lining loads. The proposed method is recommended as an approximate guideline for the design of tunnels, but the results should be confirmed by field measurements due to the uncertainties of the ground and lining properties and the construction procedures, This is the reason that in-situ monitoring should be an integral part of the design procedure.

• Journal of Korea Water Resources Association
• /
• v.48 no.4
• /
• pp.299-308
• /
• 2015

The aim of this study is that a theoretical formula for estimating the one-dimensional longitudinal dispersion coefficient is derived based on a transverse distribution equation for the depth averaged stream-wise velocity in open channel. In "Part I. Theoretical equation for stream-wise velocity" which is the former volume of this article, the velocity distribution equation is derived analytically based on the Shiono-Knight Method (SKM). And then incorporating the velocity distribution equation into a triple integral formula which was proposed by Fischer (1968), the one-dimensional longitudinal dispersion coefficient can be derived theoretically in "Part II. Longitudinal dispersion coefficient" which is the latter volume of this article. The proposed equations for the velocity distribution and the longitudinal dispersion coefficient are verified by using observed data set. As a result, the non-dimensional longitudinal dispersion coefficient is inversely proportional to square of the Manning's roughness coefficient and the non-dimensional transverse dispersion coefficient, and is directly proportional to square of the aspect ratio (channel width to depth).