• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.91-98
• /
• 2004

In two dimensional incompressible flows, a preconditioning technique called Hierarchical Iterative Procedure(HIP) has been implemented on a stabilized finite element formulation. The stabilization has been peformed by a modified residual method proposed by Illinca et. al.[3]. The stabilization which is necessary to escape from the LBB constraint renders an equal order formulation. In this paper, we increased the order of interpolation whithin an element up to cubic. The conjugate gradient squared(CGS) method is used for the outer iteration, and the HIP for the preconditioning for the incompressible Navier-Stokes equation. The hierarchical elements has been used to achieve a higher order accuracy in fluid flow analyses, but a proper efficient iterative procedure for higher order finite element formulation has not been available so far. The numerical results by the present HIP for the lid driven cavity flow showed the present procedure to be stable, very efficient and useful in flow analyses in conjunction with hierarchical elements.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.6 no.4
• /
• pp.1000-1013
• /
• 2014

This paper presents some of the efforts by the authors towards numerical prediction of springing of ships. A time-domain Higher Order Boundary Element Method (HOBEM) based on cubic shape function is first presented to solve a complete second-order problem in terms of wave steepness and ship motions in a consistent manner. In order to avoid high order derivatives on the body surfaces, e.g. mj-terms, a new formulation of the Boundary Value Problem in a body-fixed coordinate system has been proposed instead of traditional formulation in inertial coordinate system. The local steady flow effects on the unsteady waves are taken into account. Double-body flow is used as the basis flow which is an appropriate approximation for ships with moderate forward speed. This numerical model was used to estimate the complete second order wave excitation of springing of a displacement ship at constant forward speeds.

• Structural Engineering and Mechanics
• /
• v.50 no.2
• /
• pp.181-199
• /
• 2014

Locking phenomenon is a mesh problem and can be staved off with mesh refinement. If the studier is not preferred going to the solution with increasing mesh size or the computer memory can stack over flow than using higher order plate finite element or using integration techniques is a solution for this problem. The purpose of this paper is to show the shear locking phenomenon can be avoided by increase low order finite element mesh size of the plates and to study shear locking-free analysis of thick plates using Mindlin's theory by using higher order displacement shape function and to determine the effects of various parameters such as the thickness/span ratio, mesh size on the linear responses of thick plates subjected to uniformly distributed loads. A computer program using finite element method is coded in C++ to analyze the plates clamped or simply supported along all four edges. In the analysis, 4-, 8- and 17-noded quadrilateral finite elements are used. It is concluded that 17-noded finite element converges to exact results much faster than 8-noded finite element, and that it is better to use 17-noded finite element for shear-locking free analysis of plates.

• Steel and Composite Structures
• /
• v.28 no.4
• /
• pp.415-426
• /
• 2018

This paper presents a multi-layer finite element for buckling and free vibration analyses of laminated beams based on a higher-order layer-wise theory. An N-layer beam element with (9N + 7) degrees-of-freedom is proposed for analyses. Delamination and slip between the layers are not allowed. Element matrices for the single- and multi-layer beam elements are derived by Lagrange's equations. Buckling loads and natural frequencies are calculated for different end conditions and lamina stacking. Comparisons are made to show the accuracy of proposed element.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.4
• /
• pp.559-566
• /
• 2003

In this study, we present a new highly accurate two-dimensional curved composite beam element. The present element, which is based on the Hellinger-Reissner variational principle and classical lamination theory, employs consistent stress parameters corresponding to cubic displacement polynomials with additional nodeless degrees to resolve the numerical difficulties due to the spurious constraints. The stress parameters are eliminated and the nodeless degrees are condensed out to obtain the (9x9) element stiffness matrix. It should be noted that the stacking sequences without transverse deformation to the load plane makes a two dimensional analysis of curved composite beams practically useful . Several numerical examples confirm the superior locking-free behavior of the present higher-order laminated curved beam element.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.12
• /
• pp.2982-2994
• /
• 1993

Low-velocity impact responses of composite laminates are investigated using the finite element method based on various theories. In two-dimensional nonlinear analysis, a displacement field considering higher order shear deformation and large deflection of the laminate is assumed and a finite element formulation is developed using a C$^{o}$-continuous 9-node plate element. Also, three-dimensional linear analysis based on the infinitesimal strain-displacement assumptions is performed using 8-node brick elements with incompatible modes. A modified Hertzian contact law is incorporated into the finite element program to evaluate the impact force. In the time integration, the Newmark constant acceleration algorithm is used in conjuction with successive iterations within each time step. Numerical results from static analysis as well as the impact response analysis are presented including impact force histories, deflections, strains in the laminate. Impact responses according to two typical low-velocity impact conditions are compared each other.

• Structural Engineering and Mechanics
• /
• v.13 no.6
• /
• pp.695-711
• /
• 2002

Here, the dynamic response characteristics of thick cross-ply laminated composite cylindrical shells are studied using a higher-order displacement model. The formulation accounts for the nonlinear variation of the in-plane and transverse displacements through the thickness, and abrupt discontinuity in slope of the in-plane displacements at any interface. The effect of inplane and rotary inertia terms is included. The analysis is carried out using finite element approach. The influences of various terms in the higher-order displacement field on the free vibrations, and transient dynamic response characteristics of cylindrical composite shells subjected to thermal and mechanical loads are analyzed.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2004.05a
• /
• pp.216-221
• /
• 2004

In two dimensional incompressible flaws, a preconditioning technique called Hierarchical Iterative Procedure(HIP) has been implemented on a SUPG finite element formulation. By using the SUPG formulation, one can escape from the LBB constraint and hence achieve an equal order formulation. In this paper, we increased the order of interpolation up to cubic. The conjugate gradient squared(CGS) method is used for the outer iteration, and the HIP for the preconditioning for the incompressible Navier-Stokes equation. The hierarchical elements has been used to achieve a higher order accuracy in fluid flaw analyses, but a proper efficient iterative procedure for higher order finite element formulation has not been available so far. The numerical results by the present HIP for the lid driven cavity flaw showed the present procedure to be stable, very efficient and useful in flaw analyses in conjunction with hierarchical elements.

• Journal of Ocean Engineering and Technology
• /
• v.17 no.5 s.54
• /
• pp.11-18
• /
• 2003

In two-dimensional incompressible flows, a preconditioning technique called Hierarchical Iterative Procedure (HIP) has been implemented on a SUPG finite element formulation. By using the SUPG formulation, one can escape from the LBB constraint hence, achieving an equal order formulation. In this paper, we increased the order of interpolation up to cubic. The conjugate gradient squared (CGS) method is used for the outer iteration, and the HIP for the preconditioning for the incompressible Navier-Stokes equation. The hierarchical elements have been used to achieve a higher order accuracy in fluid flow analyses, but a proper and efficient iterative procedure for higher order finite element formulation has not been available, thus far. The numerical results by the present HIP for the lid driven cavity flow showed the present procedure to be stable, very efficient, and useful in flow analyses, in conjunction with hierarchical elements.

• Journal of computational fluids engineering
• /
• v.9 no.3
• /
• pp.57-65
• /
• 2004

The incompressible Navier-Stokes equations in two dimensions are stabilized by a modified residual method, and then discretized by hierarchical elements. The stabilization is necessary to escape from the Ladyzhenskaya-Babuska-Brezzi(LBB) constraint and hence to achieve an equal order formulation. To expedite a standard iterative method such as the conjugate gradient squared(CGS) method, a preconditioning technique called the Hierarchical Iterative Procedure(HIP) has been applied. In this paper, we increased the order of interpolation within an element up to cubic. The hierarchical elements have been used to achieve a higher order accuracy in fluid flow analyses, but a proper efficient iterative procedure for higher order finite element formulation has not been available so far The numerical results by the present HIP for the lid driven cavity flow and others showed the present procedure to be stable, very efficient and useful in flow analyses in conjunction with hierarchical elements.