• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.32 no.6A
• /
• pp.411-418
• /
• 2012

In this study, behavior of unpatched and patched cylindrical shells with through-wall cracks has been estimated using numerical experiments, and patching effect of them has been investigated according to various patching parameters. To show credibility of numerical models considered, two ways such as h- and p-methods have been adopted. Also, domain integral method and virtual crack extension method have been considered to calculate energy release rates based on linear elastic fracture mechanics. For examples, the unpatched cylindrical shells with circumferential cracks under remote tension have firstly been analyzed to show the validity of finite element modeling with h-method or p-method, and then the results have been compared with literature values published. Next, the sensitive analysis of patch repaired problems in terms of thickness of patch and adhesive, shear modulus of adhesive, composite material type of patch, crack length, etc. has been carried out.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.32 no.1A
• /
• pp.11-18
• /
• 2012

This paper deals with the numerical determination of the stress intensity factors of cracked aluminum plates under the mixed mode of $K_I$ and $K_{II}$ in glass-epoxy fiber reinforced composites. For the stress intensity factors, two different models are reviewed such as VCCT and two-step extension method. The p-convergent partial layerwise model is adopted to determine the fracture parameters in terms of energy release rates and stress intensity factors. The p-convergent approach is based on the concept of subparametric element. In assumed displacement field, strain-displacement relations and 3-D constitutive equations of a layer are obtained by combination of 2-D and 1-D higher-order shape functions. In the elements, Lobatto shape functions and Gauss-Lobatto technique are employed to interpolate displacement fields and to implement numerical quadrature. Using the models and techniques considered, effects of composite laminate configuration according to inclined angles and adhesive properties on the performance of bonded composite patch are investigated. In addition to these, the out-of-plane bending effect has been investigated across the thickness of patch repaired laminate plates due to the change of neutral axis. The present model provides accuracy and simplicity in terms of stress intensity factors, stress distribution, number of degrees of freedom, and energy release rates as compared with previous works in literatures.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.2
• /
• pp.221-228
• /
• 1989

The objective of this paper is to develop a numerical technique for analyzing crack driving forces in multi-phase materials. The analysis was based on finite element method coupled with a virtual crack extension technique which is known as the most efficient tool in computational fracture mechanics analysis. The modified J-integral method, proposed by Miyamoto and Kikuchi for the analysis of dual-phase material was carried out by subtracting the J-values for contours surrounding each phase boundary from the J-values for overall contour. It was shown that the proposed numerical procedure, based on the modified J-integral coupled with a virtual crack extension technique, can be used as an effective numerical tool for determining crack driving forces in multi-phase materials.

• Computational Structural Engineering
• /
• v.6 no.4
• /
• pp.57-66
• /
• 1993

The p-version crack model based on integrals of Legendre polynomial and virtual crack extension method is proposed with its potential for application to stress intensity factor computations in linear elastic fracture mechanics. The main advantage of this model is that the data preparation effort is minimal because only a small number of elements are used and high accuracy and the rapid convergence can be achieved in the vicinity of crack tip. There are two important findings from this study. Firstly, the limit value, the strain energy of the exact solution, can be estimated with successive three p-version approximations by ascertaining that the approximations enter the asymptotic range. Secondly, the rate of convergence of p-version model is almost twice that of h-version model on the basis of uniform or quasiuniform mesh refinement for the cracked panel problem subjected to tension.

• Journal of the Korean Professional Engineers Association
• /
• v.20 no.3
• /
• pp.5-14
• /
• 1987

연속체의 해석에 있어서, 특별한 경우를 제외하고는, 구조물의 개략적인 거동을 파악해야 될 경우가 종종 있다. 이러한 요구에 부응하기 위해서 강체요소법(Rigid Element Method)이라 불리우는 새로운 해석법이 개발되었다. 강체요소법은 원래 평정연구실에서 벽식프리캐스트 철근콘크리트 구조물의 탄소성해석을 하기 위해서 개발된 해석법에 착안하여, 내수벽과 같은 연속체에 적용함으로서 시작된 수치해석법이다. 그 후 저자들은 도통쉘, 구형쉘 혹은 이들이 조합된 쉘구조물에 적용할 수 있도록 개발 확장하였다. 강체요소법의 기본개념은 연속체의 분해된 각 요소를 강체(rigid body)라고 가정하고, 각 요소들은 요소의 강성으로 치환된 가상스프링으로 서로 연결되어 있다고 가정하여, 이 가상스프링의 거동을 평가함으로서 전체구조물의 거동을 파악하는 해석법이다. 이때 요소의 주변에 취해진 스프링은 해석을 단순화하기 위해서 축력, 면내전단력 및 면외전단력만을 전달한다고 가정하고, 요소의 강체변위(자유도)는 요소내의 임의의 한 점에서 취하며, 이 점에서의 강체변위(rigid displacements)는 요소의 주변에 취해진 스프링을 통하여 다른 요소로 전달된다. 상기와 같은 강체요소법의 개념을 연속체의 탄성 및 탄소성해석에 적용하면, 해석적 개념이 단순할 뿐만 아니라 구조물 전체의 자유도수를 대폭 줄여 컴퓨터 계산시간을 절약할 수 있는 잇점이 있고, 거시적인 모델(macroscopic modeling)과 미시적인 모델 (microscopic modeling)의 중간적인 성격을 가지기 때문에 구조물의 파괴상황에 대해서도 그 개략을 파악할 수 있다. 본 논문에서는 강체요소법을 보다 일반화된 해석법으로 개발, 확장하기 위해서 종전에 단층스프링시스템(single-layer spring system)으로 해석이 어려웠던 문제점들을 보완한 복층프링시스템(double-layer spring system)을 사용함으로서 휨, 비틀림의 효과를 파악할 수 있는 이론적 개념을 적용한 새로운 구요소, 원통요소 및 평면요소를 개발하고, 이러한 강체요소들의 적합매트릭스의 유도 및 해석저긴 방법을 정식화하였다. 또 휨, 비틀림 및 전단력의 효과를 고려한 사각형원통요소 및 능형원 통요소를 이용하여 원통쉘의 탄성 및 탄소성해석할 수 있는 프로그램을 개발하고, 이 프로그램으로 캔틸레버로된 연속형철근콘크리트 원통쉘의 탄성 및 탄소성해석에 적용하여 구조물의 거동에 관한 수치해석의 결과, 즉 내력의 분포, 균열의 진전, 파괴의 상황 및 변형의 상태 등을 파악해 보았다.

• Journal of the Korea Concrete Institute
• /
• v.20 no.2
• /
• pp.147-156
• /
• 2008

In an uniaxial compressive test of a concrete standard specimen (150$\times$300 mm) the crack initiation and extension with the stress increase are the major reason of the failure, which is similar to the breakage of the particle bonding in the simulation by using particle bonded model, especially particle flow code in 3 dimensions (PFC3D) developed by Itasca Consulting Group Inc. That is the main motive to study the possibility of an uniaxial compressive strength test simulation. It is important to investigate the relationship between the micro-parameters and the macro-properties because the 3-dimensional particle bonded model uses the spherical particles to analyze the physical phenomena. Contact bonded model used herein has eight micro-parameters and there are five macro-properties; Young's modulus, Poisson's ratio, uniaxial compressive strength and the crack initiation stress and the ratio concerning the crack propagation with the stress. To simulate the compressive test we made quantitative relationships between the micro-parameters and the macro-properties by using the fractional factorial design and various sensitivity analyses including regression analysis, which result in the good agreement with the previous studies. Also, the stress-stain curve and the crack distribution over the specimen given by PFC3D showed the mechanical behavior of the concrete standard specimen under the uniaxial compression. It is concluded that the particle bonded model can be a good tool for the analyzing the mechanical behavior of concrete under the uniaxial compressive load.