• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2002년도 가을 학술발표회 논문집
• /
• pp.120-128
• /
• 2002

In order to analyze shell structures more accurately and effectively, a hybrid 4-node quadrilateral shell element is formulated. The element includes the frilling degrees of freedom and the independent parameter terms of the stress resultants are appropriately selected to overcome some of the shortcomings of the standard 4-node quadrilateral elements. In order to show the accuracy and convergent characteristics of the proposed shell element, three numerical examples are analyzed and the results are compared with the existed. As a result of this study, following conclusions are obtained. (1)Analysis results by the proposed element are less sensitive to the element geometric distortion. (2)The proposed element does not produce any spurious zero-energy modes

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1995년도 가을 학술발표회 논문집
• /
• pp.153-160
• /
• 1995

Mixed finite element formulations for incompressible materials show pressure oscillations or pressure modes in four-node quadrilateral elements. The criterion for the stability in the pressure solution is the so-called Babufka-Brezzi stability condition, and the four-node elements based on mixed variational principles do not appear to satisfy this condition. In this study, a pressure continuity residual based on the pressure discontinuity at element edges is used to study the stabilization of pressure solutions in bilinear displacement-constant pressure four-node quadrilateral elements. It is shown that the pressure solutions, although stable, exhibit sensitivity to the stabilization parameters.

• Structural Engineering and Mechanics
• /
• 제29권5호
• /
• pp.599-602
• /
• 2008

• Structural Engineering and Mechanics
• /
• 제33권4호
• /
• pp.485-502
• /
• 2009

A new 8-node serendipity quadrilateral plate bending element (MQP8) based on the Mindlin-Reissner theory for the analysis of thin and moderately thick plate bending problems using Integrated Force Method is presented in this paper. The performance of this new element (MQP8) is studied for accuracy and convergence by analyzing many standard benchmark plate bending problems. This new element MQP8 performs excellent in both thin and moderately thick plate bending situations. And also this element is free from spurious/zero energy modes and free from shear locking problem.

• Nuclear Engineering and Technology
• /
• 제49권1호
• /
• pp.29-42
• /
• 2017

The purpose of the present study is the presentation of the appropriate element and shape function in the solution of the neutron diffusion equation in two-dimensional (2D) geometries. To this end, the multigroup neutron diffusion equation is solved using the Galerkin finite element method in both rectangular and hexagonal reactor cores. The spatial discretization of the equation is performed using unstructured triangular and quadrilateral finite elements. Calculations are performed using both linear and quadratic approximations of shape function in the Galerkin finite element method, based on which results are compared. Using the power iteration method, the neutron flux distributions with the corresponding eigenvalue are obtained. The results are then validated against the valid results for IAEA-2D and BIBLIS-2D benchmark problems. To investigate the dependency of the results to the type and number of the elements, and shape function order, a sensitivity analysis of the calculations to the mentioned parameters is performed. It is shown that the triangular elements and second order of the shape function in each element give the best results in comparison to the other states.

• 대한기계학회논문집
• /
• 제17권12호
• /
• pp.2995-3006
• /
• 1993

An automatic mesh generation scheme has been developed for finite element analysis with two-dimensional, quadrilateral elements. The basic strategies of the method are to transform the analysis domain into loops with key nodes and the loops are recursively subdivided into subloops with the use of best split lines. Finally by using the basic loop operators, the meshes are completed. In this algorithm an eight-node loop operator is proposed, which is useful in the area where the change of element size is large and the splitting criteria for subdividing the loops have also been modified to the existing algorithms. Lines, arcs, and cubic spline curves are used to define the boundaries of analysis domain. Sample meshes for several geometries are presented to demonstrate the robustness of the algorithm.

• Structural Engineering and Mechanics
• /
• 제6권6호
• /
• pp.711-725
• /
• 1998

Mixed finite element formulations for incompressible materials show pressure oscillations or pressure modes in four-node quadrilateral elements. The criterion for the stability in the pressure solution is the so-called Babu$\check{s}$ka-Brezzi stability condition, and the four-node elements based on mixed variational principles do not appear to satisfy this condition. In this study, a pressure continuity residual based on the pressure discontinuity at element edges proposed by Hughes and Franca is used to study the stabilization of pressure solutions in bilinear displacement-constant pressure four-node quadrilateral elements. Also, a solid mechanics problem is presented by which the stability of mixed elements can be studied. It is shown that the pressure solutions, although stable, are shown to exhibit sensitivity to the stabilization parameters.

• 한국산학기술학회논문지
• /
• 제19권11호
• /
• pp.205-212
• /
• 2018

본 연구에서는 상용 해석프로그램에서 구현하기에 어려움이 있는 지반공학적 문제를 해결하기 위한 쉽고 직관적인 해석 프로그램 개발의 일환으로 계산의 정확도가 상대적으로 높은 요소를 채택한 유한요소법을 정식화 하고 해석과정을 프로그램화 하였다. 개발된 프로그램의 계산과정에 있어서의 신뢰성 확인을 위해 두 가지 예에 대한 해석을 수행하고 결과분석을 해 보았는데 첫 번째 예는 등방구속압이 요소에 작용하는 경우이고 나머지 예는 전단응력이 요소의 측면에 작용하는 경우이다. 유한요소를 구성하는 요소로는 등매개변수 사변형 요소를 채택하였는데 요소내의 변위는 요소의 절점변위와 형상함수로 표현된다. 전체좌표(global coordinate)에 의한 미분계수로 표현되는 변형률을 얻기 위해 자코비언과 자연좌표(natural coordinate)를 이용하는 계산과정을 코딩하였다. 요소의 강성행렬을 정의하는 이중적분식을 수치적분으로 변환시키기 위해 4점 가우스 구적법을 적용하였다. 개발된 프로그램의 계산과정 검증을 위해 등방구속압이 작용하는 요소에 대한 해석을 수행한 결과 요소내의 네 개의 가우스점과 요소 중앙에 대해 계산된 응력값이 등방구속압과 일치됨을 알 수 있었다. 개발된 프로그램의 계산과정 검증을 위해 전단응력이 작용하는 요소에 대한 해석을 수행한 결과 요소내에 발생되는 횡방향응력 및 연직응력이 위치에 따라 변화됨을 알 수 있었으며 외력에 대한 발생응력의 크기 및 분포양상이 합리적임을 알 수 있었다.

• Structural Engineering and Mechanics
• /
• 제14권6호
• /
• pp.699-712
• /
• 2002

In this paper, a new quadrilateral 5-node non-conforming membrane element with drilling degrees of freedom is presented. The main advantage of these elements is the relatively small number of integration points to evaluate a stiffness matrix comparing to the existing transition membrane elements (CLM elements). Moreover, the presented elements pass the patch test by virtue of the Direct Modification Method incorporated into the element formulation. The presented 5-node elements are proved to be very efficient when used in the local mesh refinement for the in-plane structures which have stress concentrations. And some numerical studies also show the good performance of the new element developed in this study.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2000년도 춘계학술대회논문집A
• /
• 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.