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

This study deals with numerical determination of stress intensity factors of adhesively patch-repaired plates with cracks at V-shaped or semicircular notches. The p-convergent anisoparametric model are considered and then three-dimensional virtual crack closure technique is presented using formulations of anisoparametric elements. In assumed displacement fields of an element, strain-displacement relations and three-dimensional constitutive equations are derived with three-dimensional hierarchical shape functions expanded from one-dimensional Lobatto functions. Transfinite mapping technique is used to represent a circular boundary. The present model provides accuracy and simplicity in terms of stress concentration factor, stress distribution, the number of degrees of freedom, and non-dimensional stress intensity factor as compared with previous works in literatures. Stress intensity factors obtained by the three-dimensional virtual crack closure technique are estimated with respect to the variation of width of finite plate, radius of notch root, angular inclination of V-shaped notch, and crack length.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.30 no.1
• /
• pp.31-37
• /
• 2017

In this paper, computation results of reducible modeling, stress recovery and energy release rate were compared with the results of VABS, Virtual Crack Closure Technique. The result of stress recovery analysis for 1-D model including the stiffness matrix is compared with stress results of three-dimensional 3-D FEM. Energy release rate of composite beam with longitudinal cracks is calculated and compare verifications of numerical analysis results of 3-D FEM and VABS. The procedure of calculating energy release rate through dimensional reduction and stress recovery is intended to be efficient and be utilized in the life-cycle of high-altitude uav's wing, wind blades and tilt rotor blade.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.4
• /
• pp.293-300
• /
• 2012

In this study, the delamination analysis has been implemented to investigate the initiation and propagation of crack in composite laminates composed of orthotropic materials. A simple modeling was achieved by moving nodal technique without re-meshing work when crack propagation occurred. This paper aims at achieving two specific objectives. The first is to suggest a very simple modeling scheme compared with those applied to conventional h-FEM based models. To verify the performance of the proposed model, analysis of double cantilever beams with composite materials was implemented and then the results were compared with reference values in literatures. The second one is to investigate the behavior of interior delamination problems using the proposed model. To complete these objectives, the full-discrete-layer model based on Lobatto shape functions was considered and energy release rates were calculated using three-dimensional VCCT(virtual crack closure technique) based on linear elastic fracture mechanics.

• 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.

• Journal of Ocean Engineering and Technology
• /
• v.16 no.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.

• Journal of the Korean Society for Precision Engineering
• /
• v.27 no.10
• /
• pp.69-76
• /
• 2010

This study examined how the fracture toughness is affected according to the variation of the initial crack length and the fiber arranged angle using FEA method and experimental method. Therefore, the energy release rates were calculated and compared by J-integral method and VCCT(Virtual Crack Closure Technique). The results of fracture toughness test verified these results. At this time, the locus method was used in order to determine the energy release rate. When the results of FEA were compared with those of experiment, all of those decreased with the increase of angle between load and the fiber arranged direction. The decrease was due to reducing maximum load and stiffness, and the reason of reduction has been judged that the inplane shear stress.

• Composites Research
• /
• v.29 no.1
• /
• pp.16-23
• /
• 2016

In this paper, numerical results of sectional analysis, stress recovery and energy release rate were compared with the results of VABS, 3-D FEM through the blade analysis library. The result of stress recovery analysis for one-dimensional model including the stiffness matrix is compared with stress results of three-dimensional FEM. We discuss the configuration of the blade analysis library and compare verifications of numerical analysis results of VABS. Blade analysis library through dimensional reduction and stress recovery is intended to be utilized in conjunction with pre- and post-processing of the analysis program of the composite blade, high-altitude uav's wing, wind blades and tilt rotor blade.