• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2010년도 정기 학술대회
• /
• pp.56-59
• /
• 2010

포화된 다공질 매체의 수치해석은 일반적으로 혼합유한요소방법(Mixed Finite Element Method)이 쓰인다. 이 혼합유한요소 방법은 고체변형과 유체의 이동을 동시에 고려하게 되는데 고체의 변형이 거의 없이 유체만 이동할 경우나 고체와 유체의 변형이 없이 간극수압만 존재할 경우에는 요소잠김현상(Element Locking)이 발생하여 혼합유한요소방법으로 해석하기에는 수치적으로 불안정해 진다. 본 논문에서는 이러한 수치적 불안정성을 해결한 스태거드 방법(Park and Tak 2010)을 소개하고 수치적 효율성을 위해 MPI(Message-Passing Interface) 라이브러리를 이용한 병렬해석 기법이 적용된다.

• Journal of Mechanical Science and Technology
• /
• 제17권6호
• /
• pp.796-804
• /
• 2003

The use of fracture mechanics has traditionally concentrated on crack growth under an opening mechanism. However, many service failures occur from cracks subjected to mixed-mode loading. Hence, it is necessary to evaluate the fatigue behavior under mixed-mode loading. Under mixed-mode loading, not only the fatigue crack propagation rate is of importance, but also the crack propagation direction. In modified range 0.3$\leq$a/W$\leq$0.5, the stress intensity factors (SIFs) of mode I and mode II for the compact tension shear (CTS) specimen were calculated by using elastic finite element analysis. The propagation behavior of the fatigue cracks of cold rolled stainless steels (STS304) under mixed-mode conditions was evaluated by using K$\_$I/ and $_{4}$ (SIFs of mode I and mode II). The maximum tangential stress (MTS) criterion and stress intensity factor were applied to predict the crack propagation direction and the propagation behavior of fatigue cracks.

• Structural Engineering and Mechanics
• /
• 제5권1호
• /
• pp.69-83
• /
• 1997

A new efficient hybrid/mixed thin~moderately thick plate bending element with 6-node (HM6-14) is formulated based on the Reissner-Mindlin plate bending theory. The convergence of this element is proved by error estimate theories and verified by patch test respectively. Numerical studies on such an element as HM6-14 demonstrate that it has remarkable convergence, invariability to geometric distorted mesh situations, to axial rotations, and to node positions, and no "locking" phenomenon in thin plate limit. The present element is suitable to many kinds of shape and thin~moderately thick plate bending problems. Further, in comparison with original hybrid/mixed plate bending element HP4, the present element yields an improvement of solutions. Therefore, it is an efficient element and suitable for the development of adaptive multi-field finite element method (FEM).

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제7권2호
• /
• pp.95-111
• /
• 2003

BDM mixed methods are obtained for a good approximation of velocity for flow equations. In this paper, we study an implementation issue of solving the algebraic system arising from the BDM mixed finite elements. First we discuss post-processing based on the use of Lagrange multipliers to enforce interelement continuity. Furthermore, we establish an equivalence between given mixed methods and projection finite element methods developed by Chen. Finally, we present the implementation of the first order BDM on rectangular grids and show it is as simple as solving the pressure equation.

• Structural Engineering and Mechanics
• /
• 제48권1호
• /
• pp.41-56
• /
• 2013

Initiation and growth of delamination is a great concern of designers of composite structures. Interface elements with de-cohesive constitutive law in the content of continuum damage mechanics can be used to predict initiation and growth of delamination in single and mixed mode conditions. In this paper, an interface element based on the cohesive zone method has been developed to simulate delaminatoin growth of post-buckled laminate under fatigue loading. The model was programmed as the user element and user material by the "User Programmable Features" in ANSYS finite element software. The interface element is a three-dimensional 20 node brick with small thickness. Because of mixed-mode condition of stress field at the delamination-front of post-buckled laminates, a mixed-mode bilinear constitutive law has been used as user material in this model. The constitutive law of interface element has been verified by modelling of a single element. A composite laminate with initial delamination under quasi-static compressive Loading available from literature has been remodeled with the present approach. Moreover, it will be shown that, the closer the delamination to the free surface of laminate, the slower the delamination growth under compressive fatigue loading. The effects of laminate configuration on delamination growth are also investigated.

• 한국전산구조공학회논문집
• /
• 제12권3호
• /
• pp.363-370
• /
• 1999

본 연구에서는 철근콘크리트기둥과 철골보로 이루어진 혼합구조 접합부의 해석에 대한 유한 유소법을 이용한 해석 모델 방법을 제시하였다. 혼합구조 접합부에서 콘크리트와 강판이 접하는 접촉면은 두 접촉면 사이를 부착과 마찰의 개념으로 표현할 수 있는 주-종속 접촉 알고리즘(master-slave contact algorithm)을 이용하여 모델링하였다. 그리고, 휨응력의 지배를 받는 강관에는 비적합 모드 요소를 사용하였다. 본 연구에서의 혼합구조 특징은 보에서 기둥으로 힘의 전달을 원활히 하기 위하여 다이아프램이 사용되었고, 이러한 혼합구조 접합부 모델링 방법에 대한 타당성을 알아보기 위하여 3차원 비선형 해석을 행하여 실험결과와 비교한 결과 잘 일치하는 결과를 얻었다.

• 정보처리학회논문지D
• /
• 제10D권7호
• /
• pp.1145-1148
• /
• 2003

공학문제 해결을 위한 수치적 프로그램에서 원하는 해와 그 해의 변이 값에 대하여 같은 수준의 오차를 유지할 수 있는 기존의 복합유한 요소방법을 소개하고 이에 대한 효과적인 프로그램 재사용을 이용한 Matrix 생성기법을 소개한다. 또한, 원하는 임의의 차수의 기저에 대한 Matrix의 자동 생성기법을 제안한다. 여기서, 자동 생성된 Matrix는 최소한의 nonzero element를 갖고, 이는 Inverse Matriix 형성에 있어서 최소오차와 효율성을 보장한다. 위에서 제안한 MatriBt 생성기법을 최소표면적 문제에 적용하여 본다.

• 대한조선학회논문집
• /
• 제34권4호
• /
• pp.99-107
• /
• 1997

본 논문에서는 박판보 구조물의 전단변형 효과를 보다 정확히 해석할 수 있는 모델인 불균일전단 뒤틀림 보이론에 대하여 혼합형 유한요소 해석법을 적용하였다. 일반적으로 알려져 있는 바와 같이 혼합형 유한요소는 구조물의 변위 및 용력을 해석함에 있어서 변위형 유한요소에 비해서는 정도면에서 보다 균형 잡힌 해석 결과를 제공하는 특성이 있음을 밝힌다. 계산예로서, 단면의 형상이 불연속적으로 변화하는 변단면 박판보에 양단 비틀림 모멘트가 작용하는 전형적인 보의 굽힘-비틀림 연성 거동 문제에 대하여 개발된 유한요소 해석모델을 적용함으로써 그 유용성을 확인하였다.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제5권1호
• /
• pp.55-78
• /
• 2001

Given the anisotropic Poisson equation $-{\nabla}{\cdot}{\mathcal{K}}{\nabla}p=f$, one can convert it into a system of two first order PDEs: the Darcy law for the flux $u=-{\mathcal{K}{\nabla}p$ and conservation of mass ${\nabla}{\cdot}u=f$. A very natural mixed finite volume method for this system is to seek the pressure in the nonconforming P1 space and the Darcy velocity in the lowest order Raviart-Thomas space. The equations for these variables are obtained by integrating the two first order systems over the triangular volumes. In this paper we show that such a method is really a standard finite element method with local recovery of the flux in disguise. As a consequence, we compare two approaches in analyzing finite volume methods (FVM) and shed light on the proper way of analyzing non co-volume type of FVM. Numerical results for Dirichlet and Neumann problems are included.

• 한국생산제조학회지
• /
• 제9권5호
• /
• pp.73-79
• /
• 2000

The applications of fracture mechanics have traditionally concentrated on cracks loaded by tensile stresses, and growing under an opening or mode I mechanism. However, many cases of failures occur from growth of cracks subjected to mixed mode loading. Several criteria have been proposed regarding the crack growth direction under mixed mode loadings. This paper is aimed at investigation of fatigue crack growth behaviour under mixed mode(I＋II) with variation of angle and pre-crack length in two dimensional branched type precrack. Especially the direction of fatigue crack propagation was predicted and effective stress intensity factor was calculated by finite element analysis(FEA. In this paper, the maximum tangential stress(MTS) criterion was used to predict crack growth direction. Not only experiment but also finite element analysis was carried out and the theoretical predictions were compared with experimental results.