• Economic and Environmental Geology
• /
• v.28 no.5
• /
• pp.503-512
• /
• 1995

Groundwater flow in fracture networks is simulated using a discrete fracture (DF) model which assume that groundwater flows only through the fracture network. This assumption is available if the permeability of rock matrix is very low. It is almost impossible to describe fracture networks perfectly, so a stochastic approach is used. The stochastic approach assumes that the characteristic parameters in fracture network have special distribution patterns. The stochastic model generates fracture networks with some characteristic parameters. The finite element method is used to compute fracture flows. One-dimensional line element is the element type of the finite elements. The simulation results are shown by dominant flow paths in the fracture network. The dominant flow path can be found from the simulated groundwater flow field. The model developed in this study provides the tool to estimate the influences of characteristic parameters on groundwater flow in fracture networks. The influences of some characteristic parameters on the frcture flow are estimated by the Monte Carlo simulation based on 30 realizations.

• Journal of the Korean Society for Precision Engineering
• /
• v.13 no.11
• /
• pp.114-123
• /
• 1996

The major objective of the present paper is to estabilish analytical technique in order to closely understand and analyze the actual shearing process. First of all, isothermal and non-isothermal FE-simulation of the shearing process are carried out using finite element software DEFORM. Based on preliminary simulation using DEFORM, the finite element program to analyze two dimensional shearing process is developed. The ductile fracture criterion and the element kill method are also used to estimate if and where a fracture will occur and to investigate the features of the sheared surface in shearing process. It can be seen that the developed program combined with the ductile fracture criterion and element kill method has enabled the achievement of FE-simulation from initial stage to final stage of shearing process. The effects of punch-die clearance on shearing process are also investigated. In order to verify the effectiveness of the proposed technique the simulation results are compared with the known expermental data. It is found that the results of the present work are in close agreement with the published experimental results.

• Structural Engineering and Mechanics
• /
• v.89 no.2
• /
• pp.199-211
• /
• 2024

The loading capacity of engineering structures/components reduces after the initiation and propagation of crack eventually leads to the final failure. Hence, it becomes essential to deal with the crack and its effects at the design and simulation stages itself, by detecting the prone area of the fracture. The phase-field (PF) method has been accepted widely in simulating fracture problems in complex geometries. However, most of the PF methods are formulated with second order continuity theoryinvolving C⁰ continuity. In the present study, PF method based on fourth-order (i.e., higher order) theory, maintaining C¹ continuity has been proposed for ductile fracture simulation. The formulation includes fourth-order derivative terms of phase field variable, varying between 0 and 1. Applications of fourth-order PF theory to ductile fracture simulation resulted in novelty in this area. The proposed formulation is numerically solved using a two-dimensional finite element (FE) framework in 3-layered manner system. The solutions thus obtained from the proposed fourth order theory for different benchmark problems portray the improvement in the accuracy of the numerical results and are well matched with experimental results available in the literature. These results are also compared with second-order PF theory and a comparison study demonstrated the robustness of the proposed model in capturing ductile behaviour close to experimental observations.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.38 no.7
• /
• pp.727-734
• /
• 2014

This paper describes multi-scale based ductile fracture simulation using finite element (FE) damage analysis. The maximum and crack initiation loads of cracked components were predicted using proposed virtual testing method. To apply the local approach criteria for ductile fracture, stress-modified fracture strain model was adopted as the damage criteria with modified calibration technique that only requires tensile and fracture toughness test data. Element-size-dependent critical damage model is also introduced to apply the proposed ductile fracture simulation to large-scale components. The results of the simulation were compared with those of the tests on SA333 Gr. 6 full-scale pipes at $288^{\circ}C$, performed by the Battelle Memorial Institute.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1997.03a
• /
• pp.59-62
• /
• 1997

Recently, quasimolecular dynamics has been successfully used to simulate the deformation characteristic of actual size material. In quasimolecular dynamics, which is an attempt to bridge the gab between atomistic and continuum simulations, molecules are aggregated into large units, called quasimolecules, to simulate the large scale material behavior. In this paper, a numerical simulation using quasimolecular dynamics has been performed to investigate the laminar composite material fracture and crack propagation behaviors in bending process of laminar composite material which is made of fictitious materials. The simulation of the bending of laminar composite material has clarified the effects of strength of material at outer surface upon the fracture behviors of the specimen.

• Journal of International Society for Simulation Surgery
• /
• v.3 no.2
• /
• pp.80-83
• /
• 2016

A 37-year-old male was assaulted and complained of severe periorbital swelling. Physical examination revealed that there were limitation of eyeball movement on upper gaze, diplopia, and hypoesthesia on the infraorbital nerve innervating region. Three-dimensional (3D) computed tomography (CT) of facial bone exhibited the fracture of orbital floor accompanying the significant amount of orbital contents' herniation extending to the far posterior part. To recover the orbital volume and restore orbital floor without threatening the optic nerve, preoperative simplified simulation was applied. The posterior margin of the fractured orbit was delineated with simulation technique using cross-linkage between the coronal and sagittal sections based on the referential axial view of the CT scans. Dissection, reduction of orbital contents, and insertion of the absorbable mesh plate molded after the prefabricated template by the simulation technique was performed. Extensive orbital floor defect was successfully reconstructed and there were no serious complications. The purpose of this report is to emphasize the necessity of preoperative simulation in case of restoring the extensive orbital floor defect.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 1994.10a
• /
• pp.229-234
• /
• 1994

The major objective of the present paper is to estabilish finite element simulation technique in order to further analyze the shearing process. For this the ductile fracture criterion and element kill method are used in the present work. It is postulated that a fracture initiation is based on the magnitude of local effective strain. The features of sheared surface are easily observed by the element kill method. The simulation results are compard with existing experimental results. It is found that the results of the present work are close agreement with the existing results.

• Nuclear Engineering and Technology
• /
• v.53 no.3
• /
• pp.932-941
• /
• 2021

This paper suggests a fracture simulation method for SFR metallic fuel pin under accident condition. Two major failure mechanisms - creep damage and eutectic penetration - are implemented in the suggested method. To simulate damaged element, stress-reduction concept to reduce stiffness of the damaged element is applied. Using the proposed method, the failure size of cladding can be predicted in addition to the failure time and failure site. To verify the suggested method, Whole-pin furnace (WPF) test and TREAT-M test conducted at Argonne National Laboratory (ANL) are simulated. In all cases, predicted results and experimental results are overall in good agreement. Based on the simulation result, the effect of eutectic-penetration depth representing failure behavior on failure size is studied.

• Computers and Concrete
• /
• v.9 no.1
• /
• pp.51-61
• /
• 2012

According to the results of three-point bending tests of rubberized concrete and plain concrete, the parameters such as total fracture energy ($G_F$), initial fracture energy ($G_f$), and tensile strength ($f_t$) are obtained for concrete material. Using ABAQUS software and a bilinear softening fictitious crack model, the crack propagation process was simulated and compared to the experimental results. It is found that the increase of AE hit count has a similar trend with the increase of energy dissipation in FEM simulation. For two types of concretes, both experimental results and numerical simulation indicate that the rubberized concrete has a better fracture resistance.

• Proceedings of the KSME Conference
• /
• 2001.11a
• /
• pp.155-160
• /
• 2001

This paper describes a probabilistic fracture mechanics (PFM) analysis based on Monte Carlo (MC) simulation. In the analysis of CANDU pressure tube, it is necessary to perform the PFM analyses based on statistical consideration of flaw generation time. A depth and an aspect ratio of initial semi-elliptical surface crack, a fracture toughness value, delayed hydride cracking (DHC) velocity, and flaw generation time are assumed to be probabilistic variables. In all the analyses, degradation of fracture toughness due to neutron irradiation is considered. Also, the failure criteria considered are plastic collapse, unstable fracture and crack penetration. For the crack growth by DHC, the failure probability was evaluated in due consideration of flaw generation time.