• 대한기계학회논문집A
• /
• 제22권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.

• 대한기계학회논문집
• /
• 제17권10호
• /
• pp.2398-2406
• /
• 1993

Abdi's numerical method(ref.13) for representing a stress singularity by shifting the mid-side nodes of isoparametric elements is reviewed. A simple technique to obtain the optimal position of the mid-side nodes in quadratic isoparametric finite element is presented. From this technique we can directly obtain the position of the side-nodes adjacent to the crack tip. It is also observed that the present technique provides good accuracy for the expression of the opening displacement and the determination of the mid-side nodes for more wide range of material properties than that obtained by Abdicant the finite element method is applied to determine stress intensity factors for pressurized crack perpendicular to and terminating at the interface of two bonded dissimilar materials. A proper definition for stress intensity factors of a crack perpendicular to bimaterial interface is provided. It is based upon a near-tip displacement solutions on the crack surface for interface crack between two dissimilar materials. Numerical testing is carried out with the eight-node and six-node elements. The results obtained are compared with the previous solutions.

• Structural Engineering and Mechanics
• /
• 제43권3호
• /
• pp.321-336
• /
• 2012

As an improvement on the isoparametric element method, the derivation presented in this paper is close to that done by Wang (1990) for the 2-D finite element. We extend this idea to solve 3-D crack problems in this paper. A new displacement modelling is constructed with local solutions of three-dimensional cracks and a quasi-compatible isoparametric element for three-dimensional fracture mechanics analysis is presented. The stress intensity factors can be solved directly by means of the present method without any post-processing. A new method for calculating the stress intensity factors of three-dimensional cracks with complex geometries and loads is obtained. Numerical examples are given to demonstrate the validity of the present method. The accuracy of the results obtained by the proposed element is demonstrated by solving several crack problems. The results illustrate that this method not only saves much calculating time but also increases the accuracy of solutions. Because this quasi-compatible finite element of 3-D cracks contains any singularities and easily meets the requirement of compatibility, it can be easily implemented and incorporated into existing finite element codes.

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

• 한국전산구조공학회논문집
• /
• 제25권5호
• /
• pp.413-420
• /
• 2012

본 논문은 2차원 선형탄성 직접 경계요소법에서 저매개변수 요소를 사용할 때 Kernel의 적분방법에 대하여 논의하였다. 일반적으로 등매개변수 요소의 경우 형상함수로 통칭되는 해의 기저함수와 요소의 적분을 위해 사용되는 사상함수를 동일하게 사용한다. 그러나 본 논문에서는 사상함수의 차수를 낮게 취하여 순수기저절점을 도입하고 그때 직접 경계요소의 Kernel을 적분하기 위한 방법이 모색되었다. 일반적으로 경계요소법의 적분 Kernel의 경우 Log수치적분과 코쉬주치(Cauchy principal value) 등을 통해 해결하는데, 본 논문에서는 대수적 조작을 통해 적분값의 정확도를 높일 수 있도록 새로운 수식을 유도하였다. 본 연구에서 저매개변수 기반의 직접 경계요소에 대한 강건성과 정확도를 검증하기 위해 2차원 타원형 편미분방정식으로 표현되는 평면응력과 평면변형문제에 대해 적용하였다. 적용 예제로는 단순연결영역(simple connected region)의 대표적 문제인 캔틸레버보와 다중연결영역(multiple connected region)의 대표적인 문제인 개구부가 있는 사각평면에 대해 각각 수치해석을 수행한 결과 대폭적인 자유도의 감소에 비해 정확도 측면에는 기존의 방법과 차이가 없음을 볼 수 있었다. 본 논문에서 제시된 방법은 기저함수 고차화 저매개변수 직접 경계요소법(subparametric high order boundary element)과 이에 기초를 둔 저매개변수 고차 이중경계요소법(subparametric high order dual boundary element)의 초석이 될 수 있을 것이다.

• 한국전산구조공학회논문집
• /
• 제17권1호
• /
• pp.59-65
• /
• 2004

응력확대계수는 균열진전경로의 수치해석적 연구에서 널리 사용되고 있다. 그러나 많은 경우에서 균열선단주위 응력의 급수전개식 중 이어지는 항은 정량적으로 중요하다. 따라서 본 연구에서는 이의 항을 계산하기 위하여 등매개 2차특이요소를 이용한 유한요소해석을 수행하였다. 일례로 단축하중을 받는 무한 직방성체 내 경사균열에 대하여 균열요소크기와 균열경사각을 달리 주어 가며 해석을 수행하였으며, 수치해석결과는 이론해와 비교하여 잘 일치하고 있다.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 추계학술대회논문집A
• /
• pp.234-239
• /
• 2001

The stress intensity factors have been widely used in numerical studies of crack growth direction. However in many cases, omissive terms of the series expansion are quantitatively significant, so we consider the computation of such terms. For this purpose, we used the finite element method with isometric quadratic quarter-point elements. For examples, infinite square plate with a slant crack subjected to a uniaxial load is analyzed. The numerical analysis were performed for the wide range of crack tip element lengths and inclined angles. The numerical results obtained are compared with the theoretical solutions. Also they were accurate and efficient.

• Structural Engineering and Mechanics
• /
• 제2권3호
• /
• pp.227-244
• /
• 1994

The purpose of this investigation is to propose a numerical simulation method of the crack propagation behavior in 3-dimensionl elastic body. The simulation method is based on the displacement-type finite element method, and the linear fracture theory is introduced. The results from the proposed method are compared with those from the structural experiments, and the good coincidences between them are shown in this paper. At the same time, 2-dimensional analysis is also done, and the results are compared with those obtained from 3-dimensional analysis and the structural experiments.

• 대한기계학회논문집
• /
• 제10권5호
• /
• pp.787-797
• /
• 1986

본 논문에서는 축대칭 선형 문제의 경계적분법에 대한 일반화한 정식화 과정 및 수치적 접근방법이 제시되었으며 정식화 과정 중 Navier 방정식의 기본해로부터 도 출되는 변위 및 표면적 Kernel을 구하는 Hankel 변환법을 이용한 $\ulcorner$직접축대칭접근법 $\lrcorner$과 3차원 Kevin 해로부터 원주경로 따라 적분한 $\ulcorner$3차원 접근법$\lrcorner$이 비교 검토되었 다.

• 한국콘크리트학회:학술대회논문집
• /
• 한국콘크리트학회 1998년도 가을 학술발표논문집(II)
• /
• pp.529-534
• /
• 1998

To find out an adequate failure criterion in two-dimensional linear elastic crack problems, finite element programs, SED, which determine stress intensity factors $K_I, K_{II}$, crack angle and peak load by the minimum strain energy density failure criterion were developed. In this program, the conventional quadratic isoparametric elements were used in all regions except the crack tip zone where triangular singular elements with 6 nodes were used. The results of SED were compared with the results of those which followed by the maximum circumferential tensile stress criteria and those by the maximum energy release rate criteria and those by Jenq and Shah`s experiments of the same geometry and material properties. The maximum energy release rate criteria were better close to those of the Jenq and Shah`s experiments than the maximum circumferential tensile stress criteria and the minimum strain energy density criteria.