• Journal of Mechanical Science and Technology
• /
• 제16권2호
• /
• pp.182-191
• /
• 2002

The reliable stress intensity factor analysis is required for fracture mechanics design or safety evaluation of mechanical joints at which cracks often initiate and grow. It has been reported that cracks in mechanical joints usually nucleate as corner cracks at the faying surface of joints and grow as elliptical arc through cracks. In this paper, three dimensional finite element analyses are performed for elliptical arc through cracks in mechanical joints. Thereafter stress intensity factors along elliptical crack front including two surface points are determined by the virtual crack closure technique. Virtual crack closure technique is a method to calculate stress intensity factor using the finite element analysis and can be applied to non-orthogonal mesh. As a result, the effects of clearance on the stress intensity factor are investigated and crack shape are then predicted.

• Journal of Mechanical Science and Technology
• /
• 제18권11호
• /
• pp.1996-2008
• /
• 2004

A seamless analysis of material behavior incorporating complex geometry and crack- tip modeling is one of greatly interesting topics in engineering and computational fracture mechanics fields. However, there are still large gaps between the industrial applications and fundamental academic studies due to a time consuming detailed modeling. In order to resolve this problem, a numerical method to calculate an energy release rate by virtual crack closure technique was proposed in this paper. Both free mesh method and finite element method have been utilized and, thereafter, robust local and global elements for various geometries and boundary conditions were generated. A validity of the proposed method has been demonstrated through a series of fracture mechanics analyses without tedious crack-tip meshing.

• 콘크리트학회논문집
• /
• 제11권5호
• /
• pp.79-88
• /
• 1999

The purpose of this paper is to develop a mathematical expression for computing crack angles based on reinforcement volumes in the longitudinal and transverse directions, member end-fixity and length-to-width aspect ratio. For this a reinforced concrete beam-column element is assumed to possess a series of potential crack planes represented by a number of differential truss elements. Depending on the boundary condition, a constant angle truss or a variable angle truss is employed to model the cracked structural concrete member. The truss models are then analyzed using the virtual work method of analysis to relate forces and deformations. Rigorous and simplified solution schemes are presented. An equation to estimate the theoretical crack angle is derived by considering the energy minimization on the virtual work done over both the shear and flexural components the energy minimization on the virtual work done over both the shear and flexural components of truss models. The crack angle in this study is defined as the steepest one among fan-shaped angles measured from the longitudinal axis of the member to the diagonal crack. The theoretical crack angle predictions are validated against experimentally observed crack angle reported by previous researchers in the literature. Good agreement between theory and experiment is obtained.

• 한국해양공학회지
• /
• 제16권6호
• /
• pp.55-59
• /
• 2002

Delamination means that cracking occurs on the interface layer between composite laminates. In this paper, to predict the delamination growth in composite laminates subjected to low-velocity impact, the unit load method was introduced, and an eighteen-node 3-D finite element analysis, based on assumed strain mixed formulation, was conducted. Strain energy release rate, necessary to determine the delamination growth, was calculated by using the virtual crack closure technique. The unit load method saves the computation time more than the re-meshing method. The virtual crack closure technique enables the strain energy release rate to be easily calculated, because information of the singular stress field near the crack tip is not required. Hence, the delamination growth in composite laminates that are subjected to low-velocity impact can be efficiently predicted using the above-mentioned methods.

• 한국정밀공학회지
• /
• 제27권3호
• /
• pp.80-88
• /
• 2010

A simple computational model for modeling of subsurface crack growth under cyclic contact loading is presented. In this model, it is assumed that the initial fatigue crack will initiate in the region of the maximum equivalent stress at certain depth under the contacting surface. The position and magnitude of the maximum equivalent stress are determined by using the equivalent contact model, which is based on the Hertzian contact conditions with frictional forces. The virtual crack extension method is used for simulation of the fatigue crack growth from the initial crack up to the formation of the surface pit due to contact fatigue. The relationships between the stress intensity factor and crack length are then determined for various combinations of equivalent contact radii and loadings.

• 대한기계학회논문집
• /
• 제13권2호
• /
• pp.221-228
• /
• 1989

본 연구에서는 가상균열진전법과 수정 J-적분식을 연계하여 사용할 수 있는 유한요소 프로그램을 개발하고, 이 프로그램을 기존 결과가 있는 2상재료에 적응하여 타당성을 검토한 후, 다상재료의 경우에도 확장하였고, 실제 상변태가 발생하는 중수(heavy water)형 원자로 압력관(pressure tube)의 균열해석에 적용하여 보았다.

• Structural Engineering and Mechanics
• /
• 제2권3호
• /
• pp.245-256
• /
• 1994

A new axisymmetric crack model is proposed on the basis of p-version of the finite element method limited to theory of small scale yielding. To this end, axisymmetric stress element is formulated by integrals of Legendre polynomial which has hierarchical nature and orthogonality relationship. The virtual crack extension method has been adopted to calculate the stress intensity factors for 3-D axisymmetric cracked bodies where the potential energy change as a function of position along the crack front is calculated. The sensitivity with respect to the aspect ratio and Poisson locking has been tested to ascertain the robustness of p-version axisymmetric element. Also, the limit value that is an exact solution obtained by FEM when degree of freedom is infinite can be estimated using the extrapolation equation based on error prediction in energy norm. Numerical examples of thick-walled cylinder, axisymmetric crack in a round bar and internal part-thorough cracked pipes are tested with high precision.

• 대한토목학회논문집
• /
• 제32권1A호
• /
• pp.39-48
• /
• 2012

본 논문에서는 V형 노치 및 반원형 노치 균열을 갖는 패치보강 적층판의 응력확대계수 산정을 위하여 수치해석적 방법을 사용한다. p-수렴 비등매개변수 모델이 고려되고, 이와 같은 비등매개변수 모델의 결과를 활용한 3차원 가상균열닫힘법에 대한 식이 표현된다. 1차원 로바토 함수로부터 확장된 3차원 계층적 형상함수를 가지고서, 임의의 요소에서의 변위장의 변위-변형률 관계와 3차원 구성방정식이 표현된다. 원형경계의 기하형상을 나태내기 위해 초유한사상기법을 사용한다. 응력집중계수, 응력분포, 자유도, 그리고 무차원 응력확대계수 등의 항목에 대해서, 제안된 모델의 정확도와 단순성이 기존의 결과들과의 비교를 통해 설명된다. 균열 적층판의 폭, 높이, 노치근입부의 반경, V형 노치의 경사각, 균열길이 등의 변화에 따른 응력확대계수가 산정된다.

• Structural Engineering and Mechanics
• /
• 제50권5호
• /
• pp.665-677
• /
• 2014

A numerical method, using a special mixed finite element associated with the virtual crack extension technique, has been developed to evaluate the energy release rate for kinking cracks. The element is two dimensional 7-node mixed finite element with 5 displacement nodes and 2 stress nodes. The mixed finite element ensures the continuity of stress and displacement vectors on the coherent part and the free edge effect. This element has been formulated starting from a parent element in a natural plane with the aim to model different types of cracks with various orientations. Example problems with kinking cracks in a homogeneous material and bimaterial are presented to assess the computational accuracies.

• 대한기계학회논문집
• /
• 제17권9호
• /
• pp.2172-2180
• /
• 1993

Under general loadings, including body forces, crack-face tractions and thermal loading, the energy release rate equation for a two-dimensional cracked body is presented. Defining the virtual crack extension as the variation of the geometry, the equation is directly derived by a shape design sensitivity of the potential energy. Although the form of the derived energy release rate equation is different from other researchers's results, the three example show that the former is exactly the same as the latter. However, the final integral equation do not involve the derivative of the displacement on the crack surface and crack tip region, thereby improving the numerical accuracy in the computation of the energy relase rate. Moreover, as it was derived from the governing equation including non-linear elasticity without special assumptions, the energy release rate of a elasto-plastic fracture can be obtained and any numerical stress analysis method can be applied.