• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.12 no.4
• /
• pp.537-549
• /
• 1999

선형탄성 파괴해석은 균열을 갖는 변형률 경화재료의 파괴거동을 예측하는데 불충분하기 때문에 최근에는 균열 선단 부에서 대규모 소성 역을 갖는 균열 체에 적용할 수 있는 많은 파괴역학개념이 제안되고 있다. 따라서, 본 연구에서는 대규모항복 조건하의 연성파괴를 보이는 평판을 정확하게 해석할 수 있는 새로운 유한요소모델을 제시하고자 한다. 균열 선단 부의 응력 장을 정의하는데 가장 지배적인 파괴매개변수인 J-적분 값과 소성 역의 크기 및 형상을 J-적분법과 등가영역적분법을 통해 파괴거동을 설명할 수 있도록 증분소성이론에 기초를 둔 p-version 유한요소해석이 채택되었다. 제안된 유한요소모델에 의한 수치해석결과는 이론 해와 h-version 유한요소해석과 비교되었다.

• Journal of Korean Society of Steel Construction
• /
• v.9 no.1 s.30
• /
• pp.81-94
• /
• 1997

본 연구에서는 유한 회전에 의한 효과를 고려한 곡선 보요소를 개발하고, 이 요소를 이용하여 공간뼈대 구조물의 기하학적 비선형 해석을 수행하였다. 이 곡선 보요소는 증분 변위장에 Rodriguez의 2차 유한 회전항을 포함시킴으로써, 유한 회전에 의한 기하학적 평형을 유지하도록 하였다. 대변형 해석을 위하여 Total Lagrangian 방법이 적용되었으며, 비선형 해석을 수행하기 위한 알고리즘으로는, 여러개의 임계점을 갖는 비선형 거동가지도 추적할 수 있도록 하중 및 변위 증분의 조합법이 사용되었다. 공간 뼈대 구조물의 해석 예제를 통하여, 기하학적 비선형 해석에서 발생하는 유한 회전에 의한 효과를 확인하고, 본 연구에서 제안한 유한요소의 효율성 및 비선형 알고리즘으로 선택한 하중 및 변위 증분의 조합법의 적용성을 입증하였다.

• Proceedings of the KIEE Conference
• /
• 2005.04a
• /
• pp.64-66
• /
• 2005

본 논문에서는 압전체의 타원 운동을 이용한 선형 Micro-Positioning Actuator(MPA)에 대한 설계 및 해석 기법을 제안하였다. 또한 MPA의 해석을 위해 3차원 압전체의 유한요소의 정식화를 통한 유한요소 해석기법의 검증을 수행하였다. 검증된 3차원 유한요소 해석을 이용하여 MPA의 공진 특성 즉 impedance, 모드, 그리고 타원운동 해석을 하였다. 마지막으로 제안된 해석기법을 이용하여 기계 시스템을 고려한 MPA의 최종 설계를 하였다.

• Journal of the Korean Geosynthetics Society
• /
• v.20 no.3
• /
• pp.11-20
• /
• 2021

In this paper, the approximately coupled method of finite element method and boundary element method to obtain efficient and accurate analysis results is proposed for a two-dimensional elasto-static problem with a geometrically abruptly changing part. As the finite element of a two-dimensional problem, three-node and four-node plane stress element is applied, and as the boundary element of a two-dimensional problem, three-node boundary element is applied. In the modeling stage, firstly, an entire analysis target object is modeled as finite elements, and then a geometrically abruptly changing part is modeled as boundary elements. The boundary element is defined using the nodes defined for modeling finite elements. In the analysis stage, finite element analysis is firstly performed on a entire analysis target object, and boundary element analysis is automatically performed afterwards. As for the boundary conditions at boundary element analysis, displacement conditions and stress conditions, which are the results of finite element analysis, are applied. As a numerical example, the analysis results for a two-dimensional elasto-static problem, a plate with a crack, are presented and investigated.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2008.05a
• /
• pp.1634-1638
• /
• 2008

하천의 2차원 흐름 해석, 유사이동 해석, 오염확산 해석을 위한 유체의 수치해석법에는 유한요소법, 유한차분법, 유한차분법의 변형인 유한체적법, 경계적분법 등이 있다. 유체에 대한 수치해석 기법으로 전통적으로 가장 많이 사용되고 있는 방법은 유한차분법이지만, 비구조적 요소망(unstructured mesh)을 이용하여 복잡한 형상을 표현하기가 상대적으로 용이한 유한요소법이 다양한 형태의 하천 해석에는 더욱 적합할 것이다. 본 연구에서는 비구조적 요소망을 advanced front method를 이용하여 생성해 보았다. Advanced front method는 해석하고자 하는 영역에 적절한 절점들을 생성한 후 삼각 요소망을 구성하는 grid based advanced front method와 절점들을 생성하지 않고 해석 영역에 삼각 요소를 바로 구성하는 element based advanced front method로 나눌 수 있다. Grid based advanced front method에서 해석 영역에 적절한 절점을 생성하는 방법으로는 일반적인 격자 구조의 절점 생성 방법을 적용하였으며 경계와의 거리가 가까운 절점은 생성되지 않으며, 삼각 요소를 구성할 때에는 두 개의 인접 절점을 비교하여 최적의 삼각 요소를 구성하게 된다. 단 두 개의 인접 절점만을 비교함으로서 비교적 빠른 시간 안에 최적의 삼각 요소망을 구성할 수 있다. 삼각 요소망을 생성한 후에는 Laplacian smoothing을 이용하여 삼각 요소망의 형질을 개선하였다. Element based advanced front method는 외부 경계에서부터 시작된 Front가 내부 영역으로 확대되어지며 각 Front에서 적절한 절점을 직접 생성하여 바로 삼각 요소를 구성하게 된다. Front에서 생성된 절점은 인접 절점들이 있는지 검색하여 인접 절점이 있다면 생성된 절점은 삭제되어지며 인접 절점이 삼각 요소를 위한 나머지 한 점으로 채택되어진다. Front는 외부 경계와 교차되어지지 않아야 하며 또한 연속된 Front를 효율적으로 관리하기 위해 list 자료 구조를 활용하였다.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1995.10a
• /
• pp.221-226
• /
• 1995

터빈 익렬, 펌프 익차, 원형 냉각탑, 치차 등과 같이 동일한 형상이 원주 방향으로 반복되어 있는 순환 대칭 구조물의 진동특성을 유한 요소법을 사용하여 해석하는 경우에 전체구조를 모델링하는 대신에 구조물을 동일한 형상의 부분구조로 분할하여 부분구조 한개만을 모델링하고 분할된 경계에서 적절한 경계조건을 부과하여 진동해석을 수행함으로서 컴퓨터 기억용량을 절감시키고 계산시간을 단축할 수 있는 방법이 널리 사용되고 있다. Orris and Petyt[1]는 부분구조의 양쪽 분할 경계면, 즉 연결 경계상에 있는 절점변위의 상관관계를 복소파동전파식을 이용해서 구하여 부분구조의 감소된 복소강성행렬 및 질량행렬을 만들고 실수부와 허수부를 분리하여 유한요소해석을 수행하는 방법을 제안하였다. 유한요소 프로그램 ANSYS[2]에서는 이와 같은 방법을 사용하고 있다. Thomas[3]는 순회 정규모드를 이용하였고, 참고문헌[4]에서는 순회행렬을 이용하였다. 또한 유한요소 프로그램 MSC/NASTRAN[5]에서는 푸리에 급수를 이용하고 유한요소 절점의 위치 및 변위를 원통 좌표계를 표현하여 순환대칭구조물의 유한요소해석을 수행할 수 있도록 되어있다. 본 논문에서는 순환 대칭구조물의 형상의 주기성과 순환성을 고려하여 이산퓨리에 변환을 이용함으로써 순환대칭구조물의 유한요소진동해석을 체계적으로 저용량의 컴퓨터에서 신속하고 정확하게 수행할 수 있는 방법을 제안하고자 한다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.1
• /
• pp.1-10
• /
• 2004

Adaptive finite element analysis, in which its solution error meets with the user defined allowable error, is recently used to improve the reliability of finite element analysis results. This adaptive analysis is composed of two procedures; one is the error estimation of an analysis result and the other is the reconstruction of finite elements. In the (p-method, an element size is controlled by relocating of nodal positions (r-method) and the order of an element shape function is determined by the hierarchical polynomial (p-method) corresponding to the clement solution error by the enhanced SPR. In order to show the effectiveness and the accuracy of the suggested rp-method, various numerical examples were analyzed and these analysis results were examined by comparing with those obtained by the existed methods.

• Computational Structural Engineering
• /
• v.4 no.2
• /
• pp.9-12
• /
• 1991

공학문제에 있어서, 해석적으로 접근할 수 없었던 많은 경우의 문제들이 유한요소법(Finite Element Methods)의 정형화된 모형화 및 해석과정을 통하여 쉽게 접근되어질 수 있었다. 최근 보다 효율적인 요소개발과 컴퓨터 기술의 발달로 유한요소법은 더욱 효과적인 해석 수단이 되어가고 있다. 그러나 지반공학 문제와 같은 무한영역 문제를 유한요소법으로 해석할 경우, 매우 큰 영역을 모형화하기 위하여 많은 수의 요소가 요구되며 이에 따른 자유도(Degree of Freedom) 수의 증가로 많은 계산시간을 요구하게 된다. 본 고는 무한영역 문제를 효과적으로 모형화하기 위하여 연구, 개발되어진 무한요소(Infinite Element)에 대하여 소개하려 한다. 무한요소의 기본개념과 강성행렬의 형성방법을 보인 후, 기초공학 문제를 예로 하여 이의 적용방법을 간략하게 설명하였다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.4
• /
• pp.409-415
• /
• 2010

A finite element model for punching shear of flat plate structures is presented. A parametric study also has been conducted to verification of influence of several parameters in terms of the flexural reinforcement ratio, slab thickness. Reisnner-Mindlin assumptions are adopted to consider of shear deformation. Layered shell element is considered for the material non-linearities. The finite element model of this study was verified comparing with existing experimental results. The model is able to predict the capacity of the flat plate structures. The punching shear of flat plate structures varied depending on the flexural reinforcement ratio, slab thickness.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.12
• /
• pp.294-302
• /
• 2018

To analyze the behavior of a soil beam under pore water pressure, the results of analytical solutions and finite element analysis (FEM) were compared quantitatively. In contrast to the results of the analytical solution, the horizontal stress obtained from the FEM did not show a symmetrical distribution. On the other hand, the horizontal stress became closer to symmetrical distribution as the number of elements of the soil beam were increased. A comparison of the horizontal stresses from the analytic solution with those obtained from Gaussian points of FEM showed that the magnitude of the tensile stress from the FEM using 3 elements was 6% of the maximum value of the analytical solution and the compressive stress from the FEM using the same elements was 37% of the maximum value of the analytical solution. The magnitude of the tensile stress from the FEM using 6 elements was 61% of the maximum value of the analytical solution and the magnitude of the compressive stress from the FEM using the elements was 83% of the maximum value of the analytical solution. Vertical stresses, which were obtained from the analytical solution, showed a continuous distribution with the depth of the soil beam, whereas the vertical stresses from the FEM showed a discrete distribution corresponding to each element. The results also showed that the average value of the vertical stresses of each element was close to that of the pore water pressure. A comparison of the vertical displacements computed at the near vertical center line of the soil beam from the FEM with those of the analytical solution showed that the magnitude of the vertical displacement from FEM using 3 elements was 35% of the value of the analytical solution and the magnitude of the vertical displacement from FEM using 6 elements was 57% of the value of the analytical solution.