• Proceedings of the Korean Geotechical Society Conference
• /
• 1991.10a
• /
• pp.179-194
• /
• 1991

This paper explains outline of a behavior analysis program for underground structures, and its application to a tunnel problem. The program can deal with elasto-plastic behavior of medium and supporting structures, discontinuous behavior due to existing joint, creation and propagation of cracks. in-situ loading condition, and incremental behavior due to stepwise excavation, etc. The program also has additional capabilities such as graphic output of mesh, displacement pattern, stress condition, and safety factor contour, and automatic mesh generation during the excavation steps.

• Structural Engineering and Mechanics
• /
• v.14 no.4
• /
• pp.401-416
• /
• 2002

This paper presents a domain decomposition method for buckling analysis of rectangular Kirchhoff plates subjected to uniaxial inplane load and with mixed edge support conditions. A plate is decomposed into two rectangular subdomains along the change of the discontinuous support conditions. The automated Ritz method is employed to derive the governing eigenvalue equation for the plate system. Compatibility conditions are imposed for transverse displacement and slope along the interface of the two subdomains by modifying the Ritz trial functions. The resulting Ritz function ensures that the transverse displacement and slope are continuous along the entire interface of the two subdomains. The validity and accuracy of the proposed method are verified with convergence and comparison studies. Buckling results are presented for several selected rectangular plates with various combination of mixed edge support conditions.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2008.03a
• /
• pp.181-188
• /
• 2008

In this study, estimation of behavioral characteristics between twin tunnels was performed through the series of laboratory experiment on the small scale tunnel model. In the model test, the experimental parameters were geological conditions, center to center distance between twin tunnels, application of discontinuous inclination. To estimated behavior of pillar and load-displacement relationship by model tests and numerical analyses.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2007.04a
• /
• pp.487-492
• /
• 2007

The enhancement of the service life of damaged or cracked structures is a major issue for researchers and engineers. The hierarchical void element with the integrals of Legend polynomials is used to characterize the fracture behavior of unpatched crack as well as repaired crack with bonded composite patches by computing the stress intensity factors and stress contours at the crack tip. The numerical approach is based on the v-version degenerate shell element including the theory of anisotropic laminated composites. Since the equivalent single layer approach is adopted in this study, the proposed element is necessary to represent a discontinuous crack part as a continuum body with zero stiffness of materials. Thus the aspect ratio of this element to represent the crack should be extremely slender. The sensitivity of numerical solution with respect to energy release rate, displacement and stress has been tested to show the robustness of hierarchical void element as the aspect ratio is increased up to 2000. The stiffness derivative method and displacement extrapolation method have been applied to calculate the stress intensity factors of Mode I problem.

• Structural Engineering and Mechanics
• /
• v.12 no.1
• /
• pp.17-34
• /
• 2001

The frictionless contact problem for a layered composite which consists of two elastic layers having different elastic constants and heights resting on two simple supports is considered. The external load is applied to the layered composite through a rigid stamp. For values of the resultant compressive force P acting on the stamp vertically which are less than a critical value $P_{cr}$ and for small flexibility of the layered composite, the continuous contact along the layer - the layer and the stamp - the layered composite is maintained. However, if the flexibility of the layered composite increases and if tensile tractions are not allowed on the interface, for P > $P_{cr}$, a separation may be occurred between the stamp and the layered composite or two elastic layers interface along a certain finite region. The problem is formulated and solved for both cases by using Theory of Elasticity and Integral Transform Technique. Numerical results for $P_{cr}$, separation initiation distance, contact stresses, distances determining the separation area, and the vertical displacement in the separation zone between two elastic layers are given.

• Structural Engineering and Mechanics
• /
• v.15 no.3
• /
• pp.331-344
• /
• 2003

The plane contact problem for two infinite elastic layers whose elastic constants and heights are different is considered. The layers lying on a Winkler foundation are acted upon by symmetrical distributed loads whose lengths are 2a applied to the upper layer and uniform vertical body forces due to the effect of gravity in the layers. It is assumed that the contact between two elastic layers is frictionless and that only compressive normal tractions can be transmitted through the interface. The contact along the interface will be continuous if the value of the load factor, ${\lambda}$, is less than a critical value. However, interface separation takes place if it exceeds this critical value. First, the problem of continuous contact is solved and the value of the critical load factor, ${\lambda}_{cr}$, is determined. Then, the discontinuous contact problem is formulated in terms of a singular integral equation. Numerical solutions for contact stress distribution, the size of the separation areas, critical load factor and separation distance, and vertical displacement in the separation zone are given for various dimensionless quantities and distributed loads.

• Nuclear Engineering and Technology
• /
• v.40 no.5
• /
• pp.429-438
• /
• 2008

The objective of the present study is the depth and layout optimizations of a single layer, high level radioactive waste repository in a discontinuous rock mass with special joint set arrangements. A single layer repository model, considering variations in the repository depths, pitches, and tunnel spacings, is used to analyze the thermomechanical interaction behavior. It is assumed that the repository is constructed in saturated granite with joints; the PWR spent fuel in a disposal canister is installed in a deposition drift which is then sealed with compacted bentonite; and the backfill material is filled in the repository tunnel. The decay heat generated by the high level radioactive wastes governs the thermomechanical behavior of the near field rock mass of the repository. The temperature and displacement behavior of the repository is influenced more by the pitch variations than the tunnel spacing and repository depth. However, the stress behavior is influenced more by the repository depth variations than the pitch and tunnel spacing. For the final selection of the tunnel spacing, pitch, and repository depth, other aspects such as the nuclide migration through a groundwater flow path, construction costs, operation costs, and so on should be considered.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2008.03a
• /
• pp.833-841
• /
• 2008

Recently, the upheaval generation in the road which is under service had been reported. Due to the upheaval generation, total 4 lanes were forced to curtail to 3 lanes, and traffic was delayed. In normal situation of cut-slopes in korea, that condition is hard to detect since most cut-slopes contain discontinuous material, that is rock. Common collapses in rock-slopes is wedge failure, plane failure and toppling failure which is all individual mechanism of discontinuous rock mass. In contrast, such upheaval in the road in front of cut-slope can be generated only when circular movement is triggered within the cut-slope. In this sense, rock-slopes barely show any kind of movement in the road locates at the front of them. Numerical analysis is general method in simulation of slope displacement and evaluation of safety. However, numerical analysis programs which are related with rock-slopes are not able to simulate such upheaval movement because that programs are based on discontinuous modeling mechanism. In addition, although numerical analysis programs which are based on FEM/FDM and thus utilize continuous modeling mechanism are able to simulate circular movement and upheaval situation, they have weakness in reflecting discontinuities of rock-slope itself. In this study, detailed in-site investigation and numerical analysis based on in-site condition were performed in order to expect upheaval movement in the road. In this procedure, the FLAC program which uses continuous modeling method was utilized, and new approach reflecting discontinuity developed toward the road with a ubiquitous joint model was tried to derive reliable analysis result.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.6 no.3
• /
• pp.11-21
• /
• 2002

The evaluation of displacement ductility is performed by direct method through tracking the inelastic hysteretic behavior of RC bridge columns subject to cyclic loading using a flexibility-based fiber element mode. To reasonably track the inelastic behavior until the RC bridge column reaches its ultimate state, the average stress-average strain relations and joint elements, which agree well with experiments, are modified and applied considering the tension stiffening behavior and discontinuous displacement between the column and its base. In addition the evaluation of displacement ductility is performed by a direct method easily applicable to numerical analysis. Locations for the integration points, values for the post-crushing concrete strength and low-cycle fatigue failure of longitudinal reinforcement that affect the calculation of yielding and ultimate displacements are proposed for the application to flexibility-based fiber element model. Since less than 10% of error occurs during the displacement ductility analysis, the yielding and ultimate displacements evaluated by the applied analysis method and model appear to be valid.

• Proceedings of the KSME Conference
• /
• 2004.11a
• /
• pp.202-207
• /
• 2004

An elastic-plastic finite element analysis of fatigue crack closure is performed for plane strain conditions. The stabilization behavior of crack opening level and the effect of mesh size on the crack opening stress are investigated. In order to obtain a stabilized crack opening level for plane strain conditions, the crack must be advanced through approximately four times the initial monotonic plastic zone. The crack opening load tends to increase with the decrease of mesh size. The mesh size nearly equal to the theoretical plane strain cyclic plastic zone size may provide reasonable numerical results comparable with experimental crack opening data. The crack opening behavior is influenced by the crack growth increment and discontinuous opening behavior is observed. A procedure to predict the most appropriate mesh size for different stress ratio is suggested. Crack opening loads predicted by the FE analysis based on the procedure suggested resulted in good agreement with experimental ones within the error of 5 %. Effect of the distance behind the crack tip on the crack opening load determined by the ASTM compliance offset method based on the load-displacement relation and by the rotational offset method based on the load-differential displacement relation is investigated. Optimal gage location and method to determine the crack opening load is suggested.