• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.5B
• /
• pp.591-594
• /
• 2008

Present study aims at the development of a reasonable boundary condition for a structure over inclined seabed in case of the numerical model with quadrilateral mesh system. The technique for the inclined impermeable/permeable boundary in the quadrilateral mesh is newly proposed. The new technique and LES-WASS-3D model (Hur and Lee, 2007) have been used for the investigation of the dynamics of fluid field, and validated through the comparison with a typical stair-type boundary condition. 3-Dimensional numerical model with Large Eddy Simulation is called LES-WASS-3D, and is able to simulate directly interaction of WAve Structure Sea bed/Sandy beach.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.11
• /
• pp.205-212
• /
• 2018

In order to overcome shortcomings of commercial analysis program for solving certain geotechnical problems, finite element method adopting isoparametric quadrilateral element was selected as a tool for analyzing soil behavior and calculating process was programmed. Two examples were considered in order to verify reliability of the developed program. One of the two examples is the case of acting isotropic confining pressure on finite element and the other is the case of acting shear stress on the sides of the finite element. Isoparametric quadrilateral element was considered as the finite element and displacements in the element can be expressed by node displacements and shape functions in the considered element. Calculating process for determining strain which is defined by derivatives using global coordinates was coded using the Jacobian and the natural coordinates. Four point Gauss rule was adopted to convert double integral which defines stiffness of the element into numerical integration. As a result of executing analysis of the finite element under isotropic confining pressure, calculated stress corresponding to four Gauss points and center of the element were equal to the confining pressure. In addition, according to the analyzed results for the element under shear stress, horizontal stresses and vertical stresses were varied with positions in the element and the magnitudes and distribution pattern of the stresses were thought to be rational.

• Structural Engineering and Mechanics
• /
• v.41 no.2
• /
• pp.205-229
• /
• 2012

This paper presents a new nine-node Lagrangian quadrilateral plate bending element (MQP9) using the Integrated Force Method (IFM) for the analysis of thin and moderately thick plate bending problems. Three degrees of freedom: transverse displacement w and two rotations ${\theta}_x$ and ${\theta}_y$ are considered at each node of the element. The Mindlin-Reissner theory has been employed in the formulation which accounts the effect of shear deformation. Many standard plate bending benchmark problems have been analyzed using the new element MQP9 for various grid sizes via Integrated Force Method to estimate defections and bending moments. These results of the new element MQP9 are compared with those of similar displacement-based plate bending elements available in the literature. The results are also compared with exact solutions. It is observed that the presented new element MQP9 is free from shear locking and produced, in general, excellent results in all plate bending benchmark problems considered.

• Computational Structural Engineering
• /
• v.3 no.4
• /
• pp.91-104
• /
• 1990

Static performance of quadrilateral plate elements was compared through numerical experiments. Sixteen plate elements were selected for comparison from the literature, which were displacement elements, equilibrium elements, mixed elements or hybrid elements based on the Kirchhoff theory or the Mindlin theory. Thin plate bending problems, such as square plate problems, rhombic plate problems, circular plate problems and cantilevered plate problems, were modeled by various meshes and solved under various kinds of boundary conditions. Kirchhoff elements were not so good as Mindlin elements in view of efficiency and convergence. Hinton's elements resulted in the best results for the problems considered with respect to efficiency, convergence and reliability but in some problems they also resulted in more or less inaccurate solutions.

• Structural Engineering and Mechanics
• /
• v.17 no.3_4
• /
• pp.379-391
• /
• 2004

This paper deals with the development of h-version adaptive mesh refinement and recovery strategy using variable-node elements and its application to various engineering field problems with 2D quadrilateral and 3D hexahedral models. The variable-node elements which have variable mid-side nodes on edges or faces are effectively used in overcoming some problems in connecting the different layer patterns of the transition zone between the refined and coarse mesh. A modified recovery technique of gradients adequate for variable-node elements and proper selection of error norms for each engineering field problems are proposed. In the region in which the error is greater than the permissible refinement error, the mesh is locally refined by subdivision. Reversely, in some parts of the domain having the error smaller than the permissible recovery error, the mesh is locally recovered (coarsened) by combination. Hierarchical structures (e.g. quadtrees and octrees) and element-based storage structures are composed to perform this adaptive process of refinement and recovery. Some numerical examples of a 3D heat conduction analysis of the concrete with hydration heat and a 2D flow analysis of vortex shedding show effectiveness and validity of the proposed scheme.

• Computational Structural Engineering
• /
• v.8 no.4
• /
• pp.189-198
• /
• 1995

The conventional finite element formulations for incompressible materials show pressure oscillations or pressure modes in four-node quadrilateral elements of commonly used displacement and pressure interpolations. The criterion for the stability in the pressure solution is the so-called Babugka-Brezzi stability condition, and the above elements do not satisfy this condition. In this study, a pressure continuity residual based on the pressure discontinuity at element interfaces is used to study the stabilization of pressure solutions in bilinear displacement-constant pressure four-node quadrilateral elements. This pressure residual is implemented in Q1P0 element derived from the conventional incompressible elasticity. The pressure solutions can be stable with the pressure residual though they exhibit sensitivity to the stabilization parameters. Parametric study for the solution stabilization is also discussed.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.5
• /
• pp.29-37
• /
• 1993

The extension of the weighted integral method in the area of stochastic finite element analysis is presented. The use of weighted integral method in numerical analysis was extended to CST(constant strain triangle) element by Deodatis to calculate the response variability of 2D stochastic systems. In this paper, the extension of the weighted integral method for general plane-elements is represented. It has been shown that the same mesh used in the deterministic FE analysis can be used in the stochastic FE analysis. Furthermore, because the CST element is a special case which has constant strain-displacement matrix the mingling of CST elements with the other quadrilateral elements in the analysis may also be possible.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.23 no.7 s.166
• /
• pp.1205-1214
• /
• 1999

A new quadrilateral plane stress element is developed for efficient and accurate analysis of thermal stress problems. It is convenient to use the same mesh and the same shape functions for thermal analysis and stress analysis. But, because of the inconsistency between deformation related strain field and thermal strain field, oscillatory responses and considerable errors in stresses are resulted in. To avoid undesired oscillations, strain approximation is enhanced by supplementing several assumed strain terms based on the variational principle. Thermal deformation is incorporated into the generalized mixed variational principle for displacement, strain and stress fields, and basic equations for the modified enhanced assumed strain method are derived. For the stress approximation of bilinear elements, the $5{\beta}$ version of Pian and Sumihara is adopted. The numerical results for several problems show that the present element behaves well and reduces oscillatory responses. it also results in almost the same magnitude of error as compared with the quadratic element.

• Korean Journal of Computational Design and Engineering
• /
• v.10 no.1
• /
• pp.61-69
• /
• 2005

In the global ship structure and vibration analysis, the FE(finite element) analysis model is required in the early design stage before the 3D CAD model is defined. And the analysis model generation process is a time-consuming job and takes much more time than the engineering work itself. In particular, ship structure has too many associated structural members such as stringers, stiffness and girders etc. These structural members should be satisfied as the constraints in analysis modeling. Therefore it is necessary to support generation of analysis model with satisfying these constraints as an automatic manner. For the effective support of the global ship analysis modeling, a method to generate analysis model using initial design information within ship design process, that hull form offset data and compartment data, is developed. In order to easily handle initial design information and FE model information, flexible data structure is proposed. An automatic quadrilateral mesh generation algorithm using initial design information to satisfy the constraints imposed on the ship structure is also proposed. The proposed data structure and mesh generation algorithm are applied for the various type of vessels for the usability test. Through this test, we have verified the stability and usefulness of this system including mesh generation algorithm.

• Proceedings of the KSME Conference
• /
• 2000.04a
• /
• pp.351-356
• /
• 2000

Various transition elements are generally used for the effective analysis of a complicated mechanical structure. In this paper, a solid-to-beam transition finite element which connects a continuum element and a $c^1-continuity$ beam element each other is proposed. The shape functions of the transition finite elements, which a 8-noded hexahedral solid element fur 3D analysis and a 4-noded quadrilateral plane element fur 2D analysis are connected to a Euler's beam element, are explicitely formulated. In order to show the effectiveness and convergence characteristics of the proposed transition elements. numerical tests are performed for various examples and their results are compared with those obtained by other methods. As the result of this study. following conclusions are obtained: (1)The proposed transition finite elements show the monotonic convergence characteristics because of having used the compatible displacement folds. (2)As being used the transition element in the finite element analysis, the finite element modelings are more convenient and the analysis results are more accurate because of the formulation characteristies of the Euler's beam element.