• Journal of the Korean Society of Safety
• /
• v.11 no.4
• /
• pp.143-150
• /
• 1996

This is an explicit-Implicit, finite element analysis for linear as well as nonlinear hygrothermal stress problems. Additional features, such as moisture diffusion equation, crack element and virtual crack extension(VCE ) method for evaluating J-integral are implemented in this program. The Linear Elastic Fracture Mechanics(LEFM) Theory is employed to estimate the crack driving force under the transient condition for and existing crack. Pores in materials are assumed to be saturated with moisture in the liquid form at the room temperature, which may vaporize as the temperature increases. The vaporization effects on the crack driving force are also studied. The Ideal gas equation is employed to estimate the thermodynamic pressure due to vaporization at each time step after solving basic nodal values. A set of field equations governing the time dependent response of porous media are derived from balance laws based on the mixture theory Darcy's law Is assumed for the fluid flow through the porous media. Perzyna's viscoplastic model incorporating the Von-Mises yield criterion are implemented. The Green-Naghdi stress rate is used for the invariant of stress tensor under superposed rigid body motion. Isotropic elements are used for the spatial discretization and an iterative scheme based on the full newton-Raphson method is used for solving the nonlinear governing equations.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.5A
• /
• pp.861-872
• /
• 2006

The element free Galerkin method is extended by the local partition of unity method to model the cohesive cracks in two dimensional continuum. The shape function of a particle whose domain of influence is completely cut by a crack is enriched by the step enrichment function. If the domain of influence contains a crack tip inside, it is enriched by a branch enrichment function which does not have the LEFM stress singularity. The discrete equations are obtained directly from the standard Galerkin method since the enrichment is only for the displacement field, which satisfies the local partition of unity. Because only particles whose domains of influence are influenced by a crack are enriched, the system matrix is still sparse so that the increase of the computational cost is minimized. The condition for crack growth in dynamic problems is obtained from the material instability; when the acoustic tensor loses the positive definiteness, a cohesive crack is inserted to the point so as to change the continuum to a discontiuum. The crack speed is naturally obtained from the criterion. It is found that this method is more accurate and converges faster than the classical meshless methods which are based on the visibility concept. In this paper, several well-known static and dynamic problems were solved to verify the method.

• Economic and Environmental Geology
• /
• v.39 no.2 s.177
• /
• pp.111-128
• /
• 2006

Three dimensional crustal structure and source features of earthquake hypocenters on the Korean peninsula were investigated using P and S-wave travel time tomography. The main goal of this research was to find Vp/Vs anomalies at earthquake hypocenters as well as those of crustal structure of basins and deep tectonic settings. This allowed fer the extrapolation of more detailed seismotectonic force from the Korean peninsula. The earthquake hypocenters were found to have high Vp/Vs ratio discrepancies (VRD) at the vertical sections. High V/p/Vs ratios were also found in the sedimentary basins and beneath the Chugaryong Rift Zone (CRZ), which was due to mantle plume that subsequently solidified with many fractures and faults which were saturated with connate water. The hypocenters of most earthquakes were found in the upper crust for Youngwol (YE), Kyongju (KE), Hongsung (HE), Kaesong (KSE), Daekwan (DKE), and Daehung (DHE) earthquakes, but near the subcrust or the Moho Discontinuity for Mt. Songni (SE), Sariwon (SRE) and Mt. Jiri (JE) earthquakes. Especially, we found hot springs of the Daekwan, Daehung and Unsan regions coincide with high VRD. Also, this cannot rule out the possibility that there are some partial meltings in the subcrust of this region. High VRD might indicate that many faults and fractures with connate water were dehydrated when earthquakes took place, reducing shear modulus in the hypocenter areas. This is can be explained by due to the fact that a point source which is represented by the moment tensor that may involve changes in volume, shear fracture, and rigidity. High Vp/Vs ratio discrepancies (VRD) were also found beneath Mt. Backdu beneath 40 km, indicating that magma chamber existed beneath Mt. Backdu is reducing shear modulus of S-wave velocity.

• Tunnel and Underground Space
• /
• v.23 no.5
• /
• pp.383-391
• /
• 2013

A number of true triaxial tests on rock samples have been conducted since the late 1960 and their results strongly suggest that the intermediate principal stress has a considerable effect on rock strength. Based on these experimental evidence, various 3-D rock failure criteria accounting for the effect of the intermediate principal stress have been proposed. Most of the 3-D failure criteria, however, are focused on the phenomenological description of the rock strength from the true triaxial tests, so that the associated strength parameters have little physical meaning. In order to confirm the likelihood that the intermediate principal stress dependency of rock strength is related to the presence of weak planes and their distribution to the preferred orientation, true triaxial tests are simulated with the transversely isotropic rock model. The conventional Mohr-Coulomb criterion is extended to its anisotropic version by incorporating the concept of microstructure tensor. With the anisotropic Mohr-Coulomb criterion, the critical plane approach is applied to calculate the strength of the transversely isotropic rock model and the orientation of the fracture plane. This investigation hints that the spatial distribution of microstructural planes with respect to the principal stress triad is closely related to the intermediate principal stress dependency of rock strength.

• Tunnel and Underground Space
• /
• v.13 no.1
• /
• pp.52-63
• /
• 2003

Discontinuities such as faults, fractures and joints in rock mass play the dominant role in the mechanical and hydraulic properties of the rock mass. The key factors that influence on the flow of groundwater are hydraulic and geometric characteristics of discontinuities and their connectivity. In this study, a program that analyzes groundwater flow in the 3D discontinuity network was developed on the assumption that the discontinuity characteristics such as density, trace length, orientation and aperture have particular distribution functions. This program generates discontinuities in a three-dimensional space and analyzes their connectivity and groundwater flow. Due to the limited computing capacity In this study, REV was not exactly determined, but it was inferred to be greater than 25$\times$25$\times$25 ㎥. By calculating the extent of aperture that influences on the groundwater flow, it was found that the discontinuities with the aperture smaller than 30% of the mean aperture had little influence on the groundwater flow. In addition, there was little difference in the equivalent hydraulic conductivity for the the two cases when considering and not considering the boundary effect. It was because the groundwater flow was mostly influenced by the discontinuities with large aperture. Among the parameters considered in this study, the length, aperture, and orientation of discontinuities had the greatest influence on the equivalent hydraulic conductivity of rock mass in their order. In case of existence of a fault in rock mass, elements of the equivalent hydraulic conductivity tensor parallel to the fault fairly increased in their magnitude but those perpendicular to the fault were increased in a very small amount at the first stage and then converged.