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

For the same configuration, the stiffness of 6-node two dimensional isoparametric is stiffer than that of 8-node two dimensional isoparametric element. This phenomenon may be called 'Relative Stiffness Stiffening Phenomenon.' In this paper, the relative stiffness stiffening phenomenon was studied, and could be corrected by modifying the position of Gauss integral points used in the numerical integration of the stiffness matrix. For the same deformation (bending) energy of 6-node and 8-node two dimensional isoparametric elements, Gauss integral points of 6-node element have to move closer, in comparison with those of 8-node element, in the case of numerical integration along the thickness direction.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.12 no.1
• /
• pp.67-74
• /
• 1999

본 문에서는 선체 중앙부의 유한요소 모델링과 진동해석이 수행되었다. 횡부재와 종통부재가 만나 3차원적으로 연결되어 있는 선체구조는 복잡한 구조적 특성 때문에 모델링에 많은 노력이 필요하다. 선수, 선미부에 비해 비교적 부재간의 접속이 간단한 중앙평행부의 진동해석과 같은 경우에는 모델링 기법을 개발해 사용할 수도 있다. 중앙부 횡부재와 종통부재가 만나는 부분의 접속성과 형상표현을 위해 keypoint, super element(SE) 개념을 도입하였고 형성된 SE 들을 isoparametric mapping 기법을 접속된 3차원 부재용으로 개선하여 유한요소로 분할하였다. 진동해석용으로 형성된 선체중앙부 요소망을 ANSYS로 가시화하였고 자유진동해석을 수행하였다.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.2
• /
• pp.271-277
• /
• 1994

The fatigue crack opening behavior is analyzed using finite element method. Because extremely fine mesh subdivision is required when using constant stress constant strain triangular element, this study uses conventional two dimensional eight node isoparametric elements. Since plasitc zone size is similar to crack propagating length per each load cycle because of relatively large element size, a new analysis model that a crack propagates every two load cycle is suggested. the opening load and crack opening displacement can be obtained accurately by this method.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.22 no.1
• /
• pp.206-214
• /
• 1998

When a crack meets bimaterial interface stress singularity depends on the elastic constants of the adjacent materials. In the present study we are going to describe the finite element formulation for problems with a crack to be embedded in the stiffer material$({\mu}_2/{\mu}_1)$. The cubic isoparametric singular element, represented by adequately shifting the mid-side nodes adjacent to the crack tip is constructed to enclose the crack tip. An alternative method to obtain the optimal position of the mid-side nodes of cubic isoparametric elements is presented. In addition, a proper definition for the stress intensity factors of a crack normal to bimaterial interface is provided. It is based upon near a tip displacement solutions. Models for numerical analysis are two dimensional elastic bodies with a through crack under plain strain. The results obtained are compared with the previous solutions.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.21 no.12
• /
• pp.1973-1983
• /
• 1997

Higher order isoparametric elements are usually used in the finite element analysis because they can represent easily the geometric shape of a complex structure ad can improve the solution quality. When these elements are used, the position of internal nodes affects greatly on the solution accuracy. Decreasing of the accuracy related to the position of internal nodes is due to non-conformal mapping often results in an unstable Jacobian value. This paper, in order to remove this difficulty, suggests a modified shape function which can establish conformal mapping between two coordinate systems. Numerical experiments with the proposed shape function show that a stable solution can be obtained without respect to the position of internal nodes, and offer convenience that one can take arbitrarily the position of internal nodes considering only the geometric shape of element boundaries.

• Structural Engineering and Mechanics
• /
• v.57 no.3
• /
• pp.473-506
• /
• 2016

The aim of this work is the development of a 2D quadrilateral isoparametric finite element model, based on a layerwise approach, for the bending analysis of sandwich plates. The face sheets and the core are modeled individually using, respectively, the first order shear deformation theory and the third-order plate theory. The displacement continuity condition at the interfaces 'face sheets-core' is satisfied. The assumed natural strains method is introduced to avoid an eventual shear locking phenomenon. The developed element is a four-nodded isoparametric element with fifty two degrees-of-freedom (52 DOF). Each face sheet has only two rotational DOF per node and the core has nine DOF per node: six rotational degrees and three translation components which are common for the all sandwich layers. The performance of the proposed element model is assessed by six examples, considering symmetric/unsymmetric composite sandwich plates with different aspect ratios, loadings and boundary conditions. The numerical results obtained are compared with the analytical solutions and the numerical results obtained by other authors. The results indicate that the proposed element model is promising in terms of the accuracy and the convergence speed for both thin and thick plates.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.10
• /
• pp.2398-2406
• /
• 1993

Abdi's numerical method(ref.13) for representing a stress singularity by shifting the mid-side nodes of isoparametric elements is reviewed. A simple technique to obtain the optimal position of the mid-side nodes in quadratic isoparametric finite element is presented. From this technique we can directly obtain the position of the side-nodes adjacent to the crack tip. It is also observed that the present technique provides good accuracy for the expression of the opening displacement and the determination of the mid-side nodes for more wide range of material properties than that obtained by Abdicant the finite element method is applied to determine stress intensity factors for pressurized crack perpendicular to and terminating at the interface of two bonded dissimilar materials. A proper definition for stress intensity factors of a crack perpendicular to bimaterial interface is provided. It is based upon a near-tip displacement solutions on the crack surface for interface crack between two dissimilar materials. Numerical testing is carried out with the eight-node and six-node elements. The results obtained are compared with the previous solutions.

• 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.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.2
• /
• pp.41-49
• /
• 1993

This paper presents a geometrically non-linear and elastoplastic F.E. formulation using a total Lagrangian approach for the two dimensional isoparametric curved beam elements. The beam element is derived by using plane stress elements. The basic element geometry is constructed using the coordinates of the nodes on the element center line and the nodal point normals. The element displacement field is described using two translations of the node on the center line and a rotation about the axes normal to the plane containing the center line of the element. The layered approach is used for the elastoplastic analysis of the plane framed structure with the arbitrary cross section. The iterative load or displacement incremental method for non-linear finite element analysis of the frame structure is used. Numerical examples are presented to demonstrate the behavior and the accuracy of the proposed beam element for geometric and elastoplastic non-linear applications. Comparisons made with present theory and other published data show that tilt' beam element products accurate results with good convergence characteristics.

• Structural Engineering and Mechanics
• /
• v.14 no.3
• /
• pp.245-262
• /
• 2002

The free vibration analysis of stiffened laminated composite plates has been performed using the layered (zigzag) finite element method based on the first order shear deformation theory. The layers of the laminated plate is modeled using nine-node isoparametric degenerated flat shell element. The stiffeners are modeled as three-node isoparametric beam elements based on Timoshenko beam theory. Bilinear in-plane displacement constraints are used to maintain the inter-layer continuity. A special lumping technique is used in deriving the lumped mass matrices. The natural frequencies are extracted using the subspace iteration method. Numerical results are presented for unstiffened laminated plates, stiffened isotropic plates, stiffened symmetric angle-ply laminates, stiffened skew-symmetric angle-ply laminates and stiffened skew-symmetric cross-ply laminates. The effects of fiber orientations (ply angles), number of layers, stiffener depths and degrees of orthotropy are examined.