• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.3 no.4
• /
• pp.109-122
• /
• 1983

An accurate thick beam element (TB4) which includes the effects of the shear deformation and rotatory inertia has been degenerated from the three dimensional continuum by employing the Timoshenko beam assumptions. The proposed TB4 element has four nodes and two degrees of freedom at each node, totally eight degrees of freedom. The transverse deflection W and plane rotation ${\theta}$ with the cubic interpolation functions are selected as nodal variables. The element characteristics are formulated by discretizing the beam equations of motion, using the Galerkin weighted residual method, and are numerically integrated by the reduced shear integration technique, using the three-point Gauss quadrature with the various shear coefficients. Several numerical examples are analyzed to demonstrate the accuracy and the monotonic convergence behavior of the proposed TB4 beam element. The result indicates that the TB4 element shows the more excellent performance and the monotonic convergence behavior than the other existing Timoshenko beam type elements for the whole range of the beam aspect ratios, in both static and free vibration analyses.

• Structural Engineering and Mechanics
• /
• v.61 no.2
• /
• pp.193-207
• /
• 2017

In this study, the vibration of an electrostatically actuated micro cantilever beam is analyzed in which a viscoelastic layer covers a portion of the micro beam length. This proposed model is considered as the main element of mass and pollutant micro sensors. The nonlinear motion equation is extracted by means of Hamilton principle, considering nonlinear shortening effect for Euler-Bernoulli beam. The non-linear effects of electrostatic excitation, geometry and inertia have been taken into account. The viscoelastic model is assumed as Kelvin-Voigt model. The motion equation is discretized by Galerkin approach. The linear free vibration mode shapes of non-uniform micro beam i.e. the linear mode shape of the system by considering the geometric and inertia effects of viscoelastic layer, have been employed as comparison function in the process of the motion equation discretization. The discretized equation of motion is solved by the use of multiple scale method of perturbation theory and the results are compared with the results of numerical Runge-Kutta approach. The frequency response variations for different lengths and thicknesses of the viscoelastic layer have been founded. The results indicate that if a constant volume of viscoelastic layer is to be deposited on the micro beam for mass or gas sensor applications, then a modified configuration may be found by using the analysis of this paper.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.209-212
• /
• 2002

Since Berkhoff proposed the mild-slope equation in 1972, it has widely been used for calculation of shallow water wave transformation. Recently, it was extended to give an extended mild-slope equation, which includes the bottom slope squared term and bottom curvature term so as to be capable of modeling wave transformation on rapidly varying topography. These equations were derived by integrating the Laplace equation vertically. In the present study, we develop a finite element model to solve the Laplace equation directly while keeping the same computational efficiency as the mild-slope equation. This model assumes the vertical variation of wave potential as a cosine hyperbolic function as done in the derivation of the mild-slope equation, and the Galerkin method is used to discretize . The computational domain was discretized with proper finite elements, while the radiation condition at infinity was treated by introducing the concept of an infinite element. The upper boundary condition can be either free surface or a solid structure. The applicability of the developed model was verified through example analyses of two-dimensional wave reflection and transmission. .

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.6 no.1
• /
• pp.61-68
• /
• 1986

A set of the shape functions which perfectly satisfy the homogeneous Euler Equations has been proposed for deep beam problems. A finite element beam model using the proposed shape functions has been derived by the Galerkin weighted residual method and used to analyze the numerical examples without reduced shear integration, to show the accuracy and efficiency of the proposed shape functions. The result shows that the finite element model using the proposed shape functions gives very accurate solutions for both static and free vibration analyses. The concept of the proposed shape functions is thought to be applied for the finite element analysis of the elasto-static problems.

• Ocean Systems Engineering
• /
• v.11 no.1
• /
• pp.59-81
• /
• 2021

The nonlinear formulation using the principle of virtual work-energy for free vibration of a large-sag extensible catenary riser in two dimensions is presented in this paper. A support at one end is hinged and the other is a free-sliding roller in the horizontal direction. The catenary riser has a large-sag configuration in the static equilibrium state and is assumed to displace with large amplitude to the motion state. The total virtual work of the catenary riser system involves the virtual strain energy due to bending, the virtual strain energy due to axial deformation, the virtual work done by the effective weight, and the inertia forces. The nonlinear equations of motion for two-dimensional free vibration in the Cartesian coordinate system is developed based on the difference between the Euler's equations in the static state and the displaced state. The linear and nonlinear stiffness matrices of the catenary riser are obtained and the eigenvalue problem is solved using the Galerkin finite element procedure. The natural frequencies and mode shapes are obtained. The results are validated with regard to the reference research addressing the accuracy and efficiency of the proposed nonlinear formulation. The numerical results for free vibration and the effect of the nonlinear behavior for catenary riser are presented.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.3 s.174
• /
• pp.803-809
• /
• 2000

For the effective analysis of an engineering problem, meshless methods which require only positioning finite points without the element meshing recently have been proposed and being studied extensively. Meshless methods have difficulty in imposing essential boundary conditions directly, because non-interpolate shape functions originated from an approximation process are used. So some techniques, which are Lagrange multiplier method, modified variational principles and coupling with finite elements and so on, were introduced in order to impose essential boundary conditions. In spite of these methods, imposition of essential boundary conditions have still many problems like as non-positive definiteness, inaccuracy and negation of meshless characteristics. In this paper, we propose a new method which modifies shape function. Through numerical tests, convergence, accuracy and validity of this method are compared with the standard EFGM which uses Lagrange multiplier method or modified variational principles. According to this study, the proposed method shows the comparable accuracy and efficiency.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.11
• /
• pp.1538-1551
• /
• 1997

Kinetic energy conservation for fixed and moving grids is examined in time-accurate finite element computation of fully unsteady inviscid flows. As numerical algorithms, fractional step method (FSM) and modified SIMPLE are used. To simulate the flow in moving grid system, arbitrary Lagrangian-Eulerian (ALE) method is adopted. In the present study, the energy conserving time integration rule for finite element algorithm is proposed and discussed schematically. It is shown that the discretization by Crank-Nicolson in time and Galerkin (central difference) in space must be used to ensure energy conservation. The developed code has been tested for a standing vortex in fixed or moving grid system, sloshing in a tank and propagation of a solitary wave, and has been shown to be a completely energy conserving algorithm.

• Structural Engineering and Mechanics
• /
• v.7 no.2
• /
• pp.181-202
• /
• 1999

Dynamic response of axisymmetric arbitrary laminated composite cylindrical shell of finite length, using three-dimensional elasticity equations are studied. The shell is simply supported at both ends. The highly coupled partial differential equations are reduced to ordinary differential equations (ODE) with variable coefficients by means of trigonometric function expansion in axial direction. For cylindrical shell under dynamic load, the resulting differential equations are solved by Galerkin finite element method, In this solution, the continuity conditions between any two layer is satisfied. It is found that the difference between elasticity solution (ES) and higher order shear deformation theory (HSD) become higher for a symmetric laminations than their unsymmetric counterpart. That is due to the effect of bending-streching coupling. It is also found that due to the discontinuity of inplane stresses at the interface of the laminate, the slope of transverse normal and shear stresses aren't continuous across the interface. For free vibration analysis, through dividing each layer into thin laminas, the variable coefficients in ODE become constants and the resulting equations can be solved exactly. It is shown that the natural frequency of symmetric angle-ply are generally higher than their antisymmetric counterpart. Also the results are in good agreement with similar results found in literatures.

• Structural Engineering and Mechanics
• /
• v.47 no.5
• /
• pp.599-622
• /
• 2013

Seismic response of two dimensional liquid tanks is numerically simulated using fully nonlinear velocity potential theory. Galerkin-weighted-residual based finite element method is used for solving the governing Laplace equation with fully nonlinear free surface boundary conditions and also for velocity recovery. Based on mixed Eulerian-Lagrangian (MEL) method, fourth order explicit Runge-Kutta scheme is used for time integration of free surface boundary conditions. A cubic-spline fitted regridding technique is used at every time step to eliminate possible numerical instabilities on account of Lagrangian node induced mesh distortion. An artificial surface damping term is used which mimics the viscosity induced damping and brings in numerical stability. Four earthquake motions have been suitably selected to study the effect of frequency content on the dynamic response of tank-liquid system. The nonlinear seismic response vis-a-vis linear response of rectangular liquid tank has been studied. The impulsive and convective components of hydrodynamic forces, e.g., base shear, overturning base moment and pressure distribution on tank-wall are quantified. It is observed that the convective response of tank-liquid system is very much sensitive to the frequency content of the ground motion. Such sensitivity is more pronounced in shallow tanks.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1999.04a
• /
• pp.321-326
• /
• 1999

In this study, fatigue crack propagation problem of plate under multi-axial loading is mainly considered To analyze this special problem, recently developed technique called EFGM(Element-Free Galerkin Method), one of the Meshfree Methods, and general fatigue crack growth raw herein Paris law are used Using the Implemented scheme, paths of fatigue cracks by constant-amplitude load fluctuation and multiple-crack growth behavior are examined. The failure mechanism of steel plate due to crack propagation is studied. As a result, an algorithm that treats multiple fatigue crack problems is proposed. A numerical example shows that the prediction of growing paths can be achieved successfully and efficiently by proposed algorithm.