• 대한수학회지
• /
• 제53권6호
• /
• pp.1411-1429
• /
• 2016

We consider the nite volume element method for elliptic problems using isoparametric bilinear elements on quadrilateral meshes. A gradient recovery method is presented by using the patch interpolation technique. Based on some superclose estimates, we prove that the recovered gradient $R({\nabla}u_h)$ possesses the superconvergence: ${\parallel}{\nabla}u-R({\nabla}u_h){\parallel}=O(h^2){\parallel}u{\parallel}_3$. Finally, some numerical examples are provided to illustrate our theoretical analysis.

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

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

• 한국공간정보시스템학회:학술대회논문집
• /
• 한국공간정보시스템학회 2004년도 춘계 학술발표회 논문집 제1권1호(통권1호)
• /
• pp.222-229
• /
• 2004

This study purports to analyze equally stressed surfaces in tension-membrane structures through a geometrically nonlinear approach. It adopts the formulation of a 4-node quadrilateral isoparametric plane stress element considering the orthotropic characteristic of membrane textures. Tension structures, which include cables and tension membranes, such as a cable dome initially exhibit unstable conditions because no initial internal stiffness such as bending stiffness is present. Such a structural system requires prestressing to the tension members to attain a stable state. A tension-membrane structure retains a stable three dimensional curved surface as a structural shape. This analytical process for finding the geometry is referred to as Shape Finding Analysis. In order to assess the validity of this study, we examine equally stressed surfaces of saddle and catenary shape shell structures and carry out pertinent stress analyses

• 대한토목학회논문집
• /
• 제8권2호
• /
• pp.25-32
• /
• 1988

본 연구에서는 8절점 28자유도를 갖는 사각형 고차 판 유한 요소를 면내고차 변위를 고려하여 3차원 연속체로부터 유도하였다. 요소매트릭스들은 판의 운동방정식으로부터 Galerkin 가중잔차법으로 유도하고 감차적분을 수행하여 구하였다. 고차 판 유한요소를 이용하여 판의 처짐해석과 자유진동해석을 수행한 결과 판의 두께와 경게조건에 관계없이 매우 정확한 결과를 나타내었다.

• Structural Engineering and Mechanics
• /
• 제37권6호
• /
• pp.633-648
• /
• 2011

This paper investigates the aspect-ratio locking of the isoparametric 8-node quadrilateral (QUAD8) element. An important finding is that, if finite element solution is carried out with in exact arithmetic (i.e., with no truncation and round off errors), the locking tendency of the element is completely avoided even for aspect-ratios as high as 100000. The current finite element codes mostly use floating point arithmetic. Thus, they can only avoid this locking for aspect-ratios up to 100 or 1000. A novel method is proposed in the paper to avoid aspect-ratio locking in floating point computations. In this method, the offending terms of the strain-displacement matrix (i.e., $\mathbf{B}$-matrix) are multiplied by suitable scaling factors to avoid ill-conditioning of stiffness matrix. Numerical examples are presented to demonstrate the efficacy of the method. The examples reveal that aspect-ratio locking is avoided even for aspect-ratios as high as 100000.

• Structural Engineering and Mechanics
• /
• 제57권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.

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

• Journal of Mechanical Science and Technology
• /
• 제15권11호
• /
• pp.1499-1506
• /
• 2001

Various transition elements are used in general for the effective finite element analysis of complicated mechanical structures. In this paper, a solid-to-beam transition finite element, which can b e used for connecting a C1-continuity beam element to a continuum solid element, is proposed. The shape functions of the transition finite element are derived to meet the compatibility condition, and a transition element equation is formulated by the conventional finite element procedure. In order to show the effectiveness and convergence characteristics of the proposed transition element, numerical tests are performed for various examples. As a result of this study, following conclusions are obtained. (1) The proposed transition element, which meets the compatibility of the primary variables, exhibits excellent accuracy. (2) In case of using the proposed transition element, the number of nodes in the finite element model may be considerably reduced and the model construction becomes more convenient. (3) This formulation method can be applied to the usage of higher order elements.

• 대한기계학회논문집A
• /
• 제26권5호
• /
• pp.904-913
• /
• 2002

In this paper, the dynamic stress intensity factors of fracture mechanics are numerically computed in time domain using the FEM. For which the finite element formulations are derived applying the Galerkin method to the equations of motion in convolution integral as has been presented in the previous paper. To assure the strain fields of r$^{-1}$ 2/ singularity near the crack tip, the triangular quarter-point singular elements are imbedded in the finite element mesh discretized by the isoparametric quadratic quadrilateral elements. Two-dimensional problems of the elastodynamic fracture mechanics under the impact load are solved and compared with the existing numerical and analytical solutions, being shown that numerical results of good accuracy are obtained by the presented method.

• Structural Engineering and Mechanics
• /
• 제5권4호
• /
• pp.471-490
• /
• 1997

In this paper, the classical Irons' patch tests which have been generally accepted for the convergence proof of a finite element are performed for Mindlin plate bending elements with a special emphasis on the nonconforming elements. The elements considered are 4-node and 8-node quadrilateral isoparametric elements which have been dominantly used for the analyses of plate bending problems. It was recognized from the patch tests that some nonconforming Mindlin plate elements pass all the cases of patch tests even though nonconforming elements do not preserve conformity. Then, the clues for the Mindlin plate element to pass the Irons' patch tests are investigated. Also, the convergent characteristics of some nonconforming Mindlin plate elements that do not pass the Irons' patch tests are examined by weak patch tests. The convergence tests are performed on the benchmark numerical problems for both nonconforming elements which pass the patch tests and which do not. Some conclusions on the relationship between the patch test and convergence of nonconforming Mindlin plate elements are drawn.