• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.6
• /
• pp.519-525
• /
• 2009

In this study a cohesive zone model was used to simulate the delamination phenomena which occurs by a successive crack initiation and propagation in composite laminates. The cohesive zone model was incorporated to the classical finite element method via cohesive element formulation and then implemented into the user-subroutine UEL of a commercial finite element program Abaqus. To validate the formulation and implementation of the cohesive element the finite element results were compared with the experimental data of double cantilever beam and end notched flexure tests. The numerical results well agree with the experimental load-displacement curves. Also the effect of the elastic stiffness and the size of the cohesive element on the global load-displacement curves were studied numerically. To minimize the mesh-dependency of the crack propagation path and eliminate the zig-zag patterns in the load-displacement curve, cohesive elements should be refined at the crack-tip.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.9 no.3
• /
• pp.65-73
• /
• 1989

Continuum formulations for the expressions of dynamic energy release rates and computational methods for dynamic stress intensity factors are developed for the analysis of dynamic fracture problems subjected to stress wave loading. Explicit volume integral expressions for instantaneous dynamic energy release rates are derived by modeling virtual crack extensions with the dynamic Eulerian-Lagrangian kinematic description. In the finite element applications a finite region around a crack-tip is modeled by using quarter-point singular isoparametric elements, and the volume integrals are evaluated for each crack-tip element during virtual crack extensions while the singularity is maintained. It is shown that the use of the present method is more reliable and accurate for the dynamic fracture analysis than that of other path-independent integral methods when the effects of stress waves are significant.

• Journal of Ocean Engineering and Technology
• /
• v.1 no.2
• /
• pp.113-122
• /
• 1987

Ships and offshore strrctures are exposed to the corrosive surroundings, and the extablishment of the design criteria and the elucidation on the influence by this environment are requested to maintain the safety and to demonstrate the function of the structure. In this paper, the fatigue-crack-growth behavior on the compact tension specimens of quenched, tempered HT80 grade steels and RA36 high tensile steels having a single edge fatigue cracked notch respectively, were investigated under the repeated tensile stress with constant stroke in sea water for the welded parts by shielded metal arc welding. Main results obtained are summerized as follows; 1. The fatigue-crack-growth rates da/dN in sea water appeared to be greater behavior than those in air environment at the same stress intensisy factor range $\DeltaK$. 2. The correlation data of da/dN$\DeltaK$ of the two kinds of high tensile steels in sea water showed no great difference, however, the correlation data of da/dN$\DeltaK/\sigma_y$($\sigma_y$ stands for yield strength of the material) showed that the fatigue-crack-growth behavior of RA36 plate is affected by active path corrosion(APC) mechanism, while that of HT80 grade plate is mainly affected by hydrogen embrittlement mechanism.

• Tunnel and Underground Space
• /
• v.17 no.3 s.68
• /
• pp.234-243
• /
• 2007

In this paper, it was aimed to enhance core processes required in establishing hydrogeological models constructed using borehole investigation results. Water Conducting Feature(WCF) information was extracted from borehole investigation, and HydroDFN model was constructed based on the WCF information. The HydroDFN model was sub-divided by cubic blocks, and equivalent permeability of each sub-divided block was calculated and compared with the results of hydraulic test at the borehole. Through these analysis processes, suggestion for identifying and prescribing WCF parameters in the construction of HydroDFN model was made.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.6
• /
• pp.1792-1798
• /
• 1991

In this paper the path independent $J_{k}$(k=1, 2) integrals are evaluated in a rectilinear anisotropic body for two dimensional case. The relationship among elastic constants are examined. Using those relationship the expression of $J_{2}$ Integral in terms of $K_{I}$ is found to be very simple.e.e.

• Proceedings of the Korean Society of Machine Tool Engineers Conference
• /
• 1997.04a
• /
• pp.156-163
• /
• 1997

Recently, attempts have been made to be join ceramics to metals in order to make up for the brittleness of ceramics. The difference in the coefficients of linear expansion of the two materials joined at high temperature will cause residual stress, which has a strong influence on the strength of the bonded joints. In this paper, the residual stress distribution and stress intensity factors of the ceramic/metal bonded joints were analyzed by 2-dimensional elastic boundary element method. Fracture toughness tests of ceramic/metal bonded joints with an interface crack were carried out. So the advanced method of quantitative strength evaluation for ceramic/metal bonded joints is to be suggested. Fracture surface and crack propagation path were observed using scanning electron microscope.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.11
• /
• pp.1835-1843
• /
• 2001

The energy release rate for an interface crack in anisotropic dissimilar materials was obtained by the eigenfunction expansion method and also was analyzed numerically by the reciprocal work contour integral method. It was shown that the results for orthotropic dissimilar materials are consistent with the other worker's results.

• Journal of the Korean Society of Manufacturing Technology Engineers
• /
• v.18 no.4
• /
• pp.401-409
• /
• 2009

This paper presents a method of calculating the stress intensity factor (K) and crack propagation direction (${\theta}_0$) at the crack-tip that is associated with delamination in the large scale integration(LSI) package. To establish a reasonable strength evaluation method and life prediction, it is necessary to assess fracture parameters under various fracture conditions. Therefore, we conducted quantitative stress singularity analysis considering thermal stress simulating the changes of crack length (a), (h) and (v) in delamination using the 2-dimensional elastic boundary element method (BEM), and from these results predicted crack propagation direction and path.

• The Journal of Engineering Geology
• /
• v.14 no.1
• /
• pp.29-46
• /
• 2004

It is important to identify geometries of fracture that act as a conduit of fluid flow for characterization of ground water flow in fractured rock. Fracture geometries control hydraulic conductivity and stream lines in a rock mass. However, we have difficulties to acquire whole geometric data of fractures in a field scale because of discontinuous distribution of outcrops and impossibility of continuous collecting of subsurface data. Therefore, it is needed to develop a method to describe whole feature of a target fracture geometry. This study suggests a new approach to develop a method to characterize on the whole feature of a target fracture geometry based on the Fourier transform. After sampling of specimens along a target fracture from borehole cores, effective frequencies among roughness components were selected by the Fourier transform on each specimen. Then, the selected effective frequencies were averaged on each frequency. Because the averaged spectrum includes all the frequency profiles of each specimen, it shows the representative components of the fracture roughness of the target fracture. The inverse Fourier transform is conducted to reconstruct an averaged whole roughness feature after low pass filtering. The reconstructed roughness feature also shows the representative roughness of the target subsurface fracture including the geometrical characteristics of each specimen. It also means that overall roughness feature by scaling up of a fracture. In order to identify the characteristics of permeability coefficients along the target fracture, fracture models were constructed based on the reconstructed roughness feature. The computation of permeability coefficient was performed by the homogenization analysis that can calculate accurate permeability coefficients with full consideration of fracture geometry. The results show a range between $10^{-4}{\;}and{\;}10^{-3}{\;}cm/sec$, indicating reasonable values of permeability coefficient along a large fracture. This approach will be effectively applied to the analysis of permeability characteristics along a large fracture as well as identification of the whole feature of a fracture in a field scale.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.22 no.7
• /
• pp.692-700
• /
• 2012

Vibration intensity has been used to localize vibration source of a vibrating system. Not only that, vibration intensity has also been used for structural diagnostic in identifying crack and mounted stiffeners. To clearly identify the location of vibration source and understand the changes of energy transmission path, clear flow visualization is required. Most of previous works used vectors to indicate the magnitude and direction of emerging vibration energy and transmission paths. However, due to the large surface area of a plate like-structure, clear transmission paths cannot be achieved using vector visualization. This becomes an issue when detail vector flow at all locations of the whole plate surface is required. In this study, streamlines visualization is used to clearly indicate the power flow transmission path at all plate surface. By using streamlines representation, not only clear transmission paths are obtained, but also improves the vector visualization which helps us to understand the changes of the energy flow especially for stiffened plates. In this study, vibration intensity computation is firstly compared to previous work to validate the vibration intensity computation. To clearly show the power flow transmission paths, streamlines representation is shown. This representation overcomes the unclear vector direction especially for stiffened plates. Different pattern of energy transmission path can be observed using streamlines representation for stiffened and unstiffened plate. The complex streamlines pattern can also be observed at high resonance frequencies which is unclear by using vector representation.