• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.9
• /
• pp.1209-1217
• /
• 2002

The use of a two dimensional hierarchical elements are investigated for the incompressible flow computation. The construction of hierarchical elements are explained by both a geometric configuration and a determination of degrees of freedom. Also a systematic treatment of essential boundary values has been developed for the degrees of freedom corresponding to higher order terms. The numerical study for the poisson problem shows that the computation with hierarchical higher order elements can increase the convergence rate and accuracy of finite element solutions in more efficient manner than the use of standard first order element. for Stokes and Cavity flow cases, a mixed version of penalty function approach has been introduced in connection with the hierarchical elements. Solutions from hierarchical elements showed better resolutions with consistent trends in both mesh shapes and the order of elements.

• Structural Engineering and Mechanics
• /
• v.58 no.6
• /
• pp.949-966
• /
• 2016

The second order analysis taking place due to non-linear behavior of the structures under the mechanical and geometric factors through implementing exact and approximate methods is an indispensible issue in the analysis of such structures. Among the exact methods is the slope-deflection method that due to its simplicity and efficiency of its relationships has always been in consideration. By solving the differential equations of the modified slope-deflection method in which the effect of axial compressive force is considered, the stiffness matrix including trigonometric entries would be obtained. The complexity of computations with trigonometric functions causes replacement with their Maclaurin expansion. In most cases only the first two terms of this expansion are used but to obtain more accurate results, more elements are needed. In this paper, the effect of utilizing higher order terms of Maclaurin expansion on reducing the number of required elements and attaining more rapid convergence with less error is investigated for the Bernoulli beam with various boundary conditions. The results indicate that when using only one element along the beam length, utilizing higher order terms in Maclaurin expansion would reduce the relative error in determining the critical buckling load and kinematic parameters in the second order analysis.

• International Journal of Ocean System Engineering
• /
• v.1 no.3
• /
• pp.155-164
• /
• 2011

The second-order motion characteristics of a semi-submersible are investigated in regular waves. A higher-order boundary element method in a frequency domain and a finite element method in a time-domain were applied to the numerical analysis of the nonlinear hydrodynamic force and motion characteristics of semi-submersibles in view point of potential flow. Various aspects of nonlinear effects on the heave and roll of a semi-submersible were numerically investigated and some selected cases were compared with the model test data.

• Earthquakes and Structures
• /
• v.16 no.3
• /
• pp.359-368
• /
• 2019

This paper focuses on the study of dynamic analysis of thick plates resting on Winkler foundation. The governing equation is derived from Mindlin's theory. This study is a parametric analysis of the reflections of the thickness / span ratio, the aspect ratio and the boundary conditions on the earthquake excitations are studied. In the analysis, finite element method is used for spatial integration and the Newmark-${\beta}$ method is used for the time integration. While using finite element method, a new element is used. This element is 17-noded and it's formulation is derived from using higher order displacement shape functions. C++ program is used for the analyses. Graphs are presented to help engineers in the design of thick plates subjected to earthquake excitations. It is concluded that the 17-noded finite element is used in the earthquake analysis of thick plates. It is shown that the changes in the aspect ratio are more effective than the changes in the aspect ratio. The center displacements of the reinforced concrete thick clamped plates for b/a=1, and t/a=0.2, and for b/a=2, and t/a=0.2, reached their absolute maximum values of 0.00244 mm at 3.48 s, and of 0.00444 mm at 3.48 s, respectively.

• International Journal of Aerospace System Engineering
• /
• v.3 no.1
• /
• pp.10-16
• /
• 2016

A trigonometric Layerwise higher order shear deformation theory (TLHSDT) is developed and implemented for free vibration and buckling analysis of laminated composite and sandwich plates by analytical and finite element formulation. The present model assumes parabolic variation of out-plane stresses through the depth of the plate and also accomplish the zero transverse shear stresses over the surface of the plate. Thus a need of shear correction factor is obviated. The present zigzag model able to meet the transverse shear stress continuity and zigzag form of in-plane displacement continuity at the plate interfaces. Hence, botheration of shear correction coefficient is neglected. In the case of analytical method, the governing differential equation and boundary conditions are obtained from the principle of virtual work. For the finite element formulation, an efficient eight noded $C^0$ continuous isoparametric serendipity element is established and employed to examine the dynamic analysis. Like FSDT, the considered mathematical model possesses similar number of variables and which decides the present models computationally more effective. Several numerical predictions are carried out and results are compared with those of other existing numerical approaches.

• Steel and Composite Structures
• /
• v.45 no.3
• /
• pp.331-348
• /
• 2022

The main goal of this article is to develop the finite element formulation based on the nonlocal strain gradient and the refined higher-order deformation theory employing a new function f(z) to investigate the static bending and free vibration of functionally graded porous (FGP) nanobeams. The proposed model considers the simultaneous effects of two parameters: nonlocal and strain gradient coefficients. The nanobeam is made by FGP material that exists in un-even and logarithmic-uneven distribution. The governing equation of the nanobeam is established based on Hamilton's principle. The authors use a 2-node beam element, each node with 8 degrees of freedom (DOFs) approximated by the C¹ and C² continuous Hermit functions to obtain the elemental stiffness matrix and mass matrix. The accuracy of the proposed model is tested by comparison with the results of reputable published works. From here, the influences of the parameters: nonlocal elasticity, strain gradient, porosity, and boundary conditions are studied.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.3 no.1
• /
• pp.89-96
• /
• 1995

The application of the 4-pole parameter method with 1 - D theory is acceptable for intake system analysis. However, the limitaion appears during the analysis of complicated intake system since this method is developed based on the plane wave thoery. For the intake system analysis, the usage of BEM(Boundary Element Method) is introduced describing its disadvantage. To combine benefits of both method. a hybrid method is introduced. This hybrid method consists of the 4-pole parameter with I-D theory and BEM. The developed method is applied to an automobile intake system analysis to obtain the transmission loss.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.18 no.4
• /
• pp.337-350
• /
• 2014

In this paper, we consider discontinuous Galerkin finite element methods with interior penalty term to approximate the solution of nonlinear parabolic problems with mixed boundary conditions. We construct the finite element spaces of the piecewise polynomials on which we define fully discrete discontinuous Galerkin approximations using the Crank-Nicolson method. To analyze the error estimates, we construct an appropriate projection which allows us to obtain the optimal order of a priori ${\ell}^{\infty}(L^2)$ error estimates of discontinuous Galerkin approximations in both spatial and temporal directions.

• Structural Engineering and Mechanics
• /
• v.31 no.5
• /
• pp.481-506
• /
• 2009

The present study is concerned with the stochastic linear free vibration study of laminated composite plate embedded with piezoelectric layers with random material properties. The system equations are derived using higher order shear deformation theory. The lamina material properties of the laminate are modeled as basic random variables for accurate prediction of the system behavior. A $C^0$ finite element is used for spatial descretization of the laminate. First order Taylor series based mean centered perturbation technique in conjunction with finite element method is outlined for the problem. The outlined probabilistic approach is used to obtain typical numerical results, i.e., the mean and standard deviation of natural frequency. Different combinations of simply supported, clamped and free boundary conditions are considered. The effect of side to thickness ratio, aspect ratio, lamination scheme on scattering of natural frequency is studied. The results are compared with those available in literature and an independent Monte Carlo simulation.

• Steel and Composite Structures
• /
• v.24 no.2
• /
• pp.151-176
• /
• 2017

This paper deals with general equations of motion for free vibration analysis response of thick three-layer doubly curved sandwich panels (DCSP) under simply supported boundary conditions (BCs) using higher order shear deformation theory. In this model, the face sheets are orthotropic laminated composite that follow the first order shear deformation theory (FSDT) based on Rissners-Mindlin (RM) kinematics field. The core is made of orthotropic material and its in-plane transverse displacements are modeled using the third order of the Taylor's series extension. It provides the potentiality for considering both compressible and incompressible cores. To find these equations and boundary conditions, Hamilton's principle is used. Also, the effect of trapezoidal shape factor for cross-section of curved panel element ($1{\pm}z/R$) is considered. The natural frequency parameters of DCSP are obtained using Galerkin Method. Convergence studies are performed with the appropriate formulas in general form for three-layer sandwich plate, cylindrical and spherical shells (both deep and shallow). The influences of core stiffness, ratio of core to face sheets thickness and radii of curvatures are investigated. Finally, for the first time, an optimum range for the core to face sheet stiffness ratio by considering the existence of in-plane stress which significantly affects the natural frequencies of DCSP are presented.