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

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.12 no.2
• /
• pp.185-199
• /
• 1999

본 연구에서는 새로운 spline 유한대판 요소를 제안하였다. 제안된 정식화는 등매개 개념에 의해 기하학적 형상과 변위장을 가정함에 있어 길이방향은 3차의 B-spline 곡선으로, 횡방향에 대해서는 Lagrange 다항식에 의해 표현된다. 이 논문은 평판과 쉘해석에 있어서의 등매개 spline 유한대판 요소의 개선에 목적을 두고 있다. 이 새로운 요소는 스트립의 내부 절점에서 6개의 자유도를 갖는 합-응력 감절점 쉘 요소로부터 유도하였다. 스트립의 기하학적 형상은 강체 회전에 대한 정의에 위배되지 않고도 두께 방향을 따라 Jacobian이 일정하다는 가정을 따랐으며 고체역학에서 정의되는 면내 회전을 penalty 함수에 의한 구속조건으로 간주하여 면내 회전에 관계된 자유도를 생성하였다. 제안된 요소에 대하여 쉘의 전형적인 문제에 대한 수치예제를 보였으며 이 스트립 요소의 성능을 평가하였다.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.11 no.4
• /
• pp.231-237
• /
• 1999

A geometric non-linear finite element formulation of spatial ocean cable under currents is presented using multiple noded curved cable elements. Tangent stiffness and mass matrices for the isoparametric cable ele¬ment are derived and the initial equilibrium state of ocean cable subjected to self-weights, buoyancy, and current as well as support motions is determined using the load incremental method. Free vibration analysis of ocean cables is performed based on the initial equilibrium configuration. Numerical examples are presented and discussed in order to demonstrate the feasibility of the present finite element method and investigate dynamic characteristics of ocean cables.

• Journal of KSNVE
• /
• v.10 no.1
• /
• pp.123-128
• /
• 2000

This paper describes the optimization of the hexagonal array transducer using finite element method. The transducer consists of the disc type sensors. Three dimensional beam patterns of each element and the array transducer are analysed using the finite element code ATILA. Beam patterns were analyzed for the disc type transducer. To optimize beam patterns of the array transducer, Chebyshev polynomial weight is applied to each element. In case of applying optimized weight, a 30 degree width beam pattern is presented at 10kHz. This paper also includes the effect of rubber filling material instead of using the water inside the transducer array.

• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.2
• /
• pp.73-78
• /
• 1992

Numerical method by Abdi et al. for obtaining a stress singularity near a crack perpendicular to the interface between two elastic materials is reviewed. More efficient and simple method to obtain crack singularity by shifting a mid-side node of 8-node isoparametric element is presented. It is observed that the present analysis provides increased accuracy for the expression of the opening displacement and the determination of the optimal position of the mid-side node for a wide range of material properties.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.11 no.1
• /
• pp.69-78
• /
• 1991

This paper concentrates on the analysis of deflection of the reinforced concrete flexural members under monotonic loading. Concrete is treated as an orthotropic nonlinear material. The concept of equivalent strain and crack strain are used to establish independent stress-strain relationships in the directions of orthotropy. Steel is modeled as an elstoplastic material, and von Mises failure criterion is applied. The finite element computer program for the nonlinear analysis of the deflection of RC flexural members under monotonic loading is developed. The accuracy and reliability of the numerical procedure is demonstrated by the FEM analysis and experiments of the under reinforced concrete beams over the entire loading range up to failure.

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

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.7 no.2
• /
• pp.21-28
• /
• 1987

The paper describes the analysis of large displacement problems including instability phenomena. The element used in this is a degenerated isoparametric shell element with eight nodes. Total Lagrangian formulation has been adopted in this study using Newton-Raphson iteration method with incremental load. The linear stability analyses performed usually for the initial position can be repeated at several advanced fundamental states on the non-linear buckling path. Thus a current estimate of the failure load is given. The numerical examples of a cylindrical panel under uniform load, simply supported plate under axial load, and clamped plate under uniform load are carried out. The examples applying degenerated isoparametric elements to bifurcation buckling and nonlinear collapse problems are also performed.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.8 no.2
• /
• pp.147-158
• /
• 2004

For stiffened plates composed of composite materials, many researchers have used a finite element method which connected isoparametric plate elements and beam elements. However, the finite element method is difficult to reflect local behavior of stiffener because beam elements are transferred stiffness for nodal point of plate elements, especially the application is limited in case of laminated composite structures. In this paper, for analysis of laminated composite stiffened plates, 3D shell elements for stiffener and plate are employed. Reissner-Mindlin's first order shear deformation theory is considered in this study. But when thickness will be thin, isoparamatric plate bending element based on the theory of Reissner-Mindlin is generated by transverse shear locking. To eliminate the shear locking and virtual zero energy mode, the substitute shear strain field is used. A deflection distribution is investigated for simple supported rectangular and skew stiffened laminated composite plates with arbitrary orientation stiffener as not only variation of slenderness and aspect ratio of the plate but also variation of skew angle of skew stiffened plates.