• Journal of Korean Tunnelling and Underground Space Association
• /
• v.10 no.4
• /
• pp.349-360
• /
• 2008

In this paper, the modeling procedure of random field with an elasto-plastic finite element algorithm and probability of failure on closely-spaced tunnel were investigated. Local average subdivision (LAS) method which can generate discrete random variables fast and accurately as well as change the resolution in certain region was used. And correlated value allocating and weighted average method were suggested to implement geometrical characteristics of tunnel. After the probability of failure on the test problem was thoroughly investigated using random finite element method, the results were compared with the deterministic strength reduction factor method and single random variable method. Of particular importance in this work, is the conclusion that the probability of failure determined by simplified probabilistic analysis, in which spatial variability is ignored by assuming perfect correlation, can be estimated from the safety factor determined by strength reduction factor method. Also, single random variable method can lead to unconservative estimates of the probability of failure.

• Proceedings of the Korean Society for Rock Mechanics Conference
• /
• 2000.09a
• /
• pp.117-128
• /
• 2000

The development of proper joint model, which can describe real phenomena exactly and still can be used easily, is one of the most important element for the analysis of the mechanical and hydraulic behavior of discontinuous rock mass. In this study, an elasto-plastic constitutive model of joint behavior considering asperity degradation was extended with the concept of first and second order asperities. The proposed model was implemented to numerical code with discrete finite joint element. The parametric study with the various asperity angles and degradation coefficients showed that the model can reproduce the shear behavior of typical rough joints well. Results of laboratory monotonic and cyclic shear tests were compared with those of numerical tests to validate the model. The hydraulic model considering the relations between gouge production and aperture was introduced to the mechanical model. In an attempt to examine the performance of the model, comparative numerical test was conducted. Permeability between joint surfaces increased rapidly at the first stage, but became nearly constant with increasing shear displacement due to gouge production and uniform variation of aperture distribution.

• Tunnel and Underground Space
• /
• v.10 no.3
• /
• pp.366-377
• /
• 2000

The development of proper joint model, which can describe real phenomena exactly and still can be used easily, is one of the most important element for the analysis of the mechanical and hydraulic behavior of discontinuous rock mass. In this study, an elasto-plastic constitutive model of joint behavior considering asperity degradation was extended with the concept of first and second order asperities. The proposed model was implemented to numerical code with discrete finite joint element. The parametric study with the various asperity angles and degradation coefficients showed that the model can reproduce the shear behavior of typical rough joints well. Results of laboratory monotonic and cyclic shear tests were compared with those of numerical tests to validate the model. The hydraulic model considering the relations between gouge production and aperture was introduced to the mechanical mode1. In an attempt to examine the performance of the model, comparative numerical test was conducted. Permeability between joint surfaces increased rapidly at the first stage, but became nearly constant with increasing shear displacement due to gouge production and uniform variation of aperture distribution.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.2
• /
• pp.175-182
• /
• 2012

In this paper we are concerned with the performance of structural stability of torsion in square cross section of a beam with holes. The critical load is defined as the smallest load at which the equilibrium of the structure fails to be stable as the load is slowly increased from zero. The beams subjected to torsion are frequently encountered in general structures and these forces influence to the stability of structure. The boundary element method is found to be very efficient and accurate for the analysis of torsion problems including complex boundary conditions with respect to its simplicity and generality. In this paper, it is required to derive the boundary element formulation for torsion problem and integrate directly on the discrete boundary. To investigate the validity of the developed computer program, three distinctly solid cross-sections which are elliptical, rectangular and triangular one are analyzed, and comparisons are made with analytical approaches where these can also be used.

• Computers and Concrete
• /
• v.8 no.1
• /
• pp.71-96
• /
• 2011

Existing reinforced concrete (RC) beams can be strengthened with externally bolted steel plates to the sides of beams. The effectiveness of this type of bolted side-plate (BSP) beam can however be affected by partial interaction between the steel plates and RC beams due to the mechanical slip of bolts. To avoid over-estimation of the flexural strength and ensure accurate prediction of the load-deformation response of the beams, the effect of partial interaction has to be properly considered. In this paper, a special non-linear macro-finite-element model that takes into account the effects of partial interaction is proposed. The RC beam and the steel plates are modelled as two different elements, interacting through discrete groups of bolts. A layered method is adopted for the formulation of the RC beam and steel plate elements, while a special non-linear model based on a kinematic hardening assumption for the bolts is used to simulate the bolt group effect. The computer program SiBAN was developed based on the proposed approach. Comparison with the available experimental results shows that SiBAN can accurately predict the partial interaction behaviour of the BSP beams. Further numerical simulations show that the interaction between the RC beam and the steel plates is greatly reduced by the formation of plastic hinges and should be considered in analyses of the strengthened beams.

• Transactions of Materials Processing
• /
• v.22 no.3
• /
• pp.150-157
• /
• 2013

Fiber metal laminates (FMLs) are layered materials comprised of thin metal sheets and fiber reinforced plastic (FRP). This paper presents the numerical study of the formability enhancement of FMLs composed of an aluminum alloy and self-reinforced polypropylene (SRPP) composite. In this study, a numerical simulation based on finite element (FE) modeling is proposed to evaluate the formability of FMLs using ABAQUS/Explicit. The FE model, which included a single layer of solid and shell elements to model the blank, used discrete layers of the solid element with a contact model and shell elements with a friction based model for the aluminum alloy-composite interface conditions. This method allowed the description of each layer of FMLs and was able to simulate the interaction between the layers. It is noted through this research that the proposed numerical simulation described properly the formability enhancement of the FMLs and the simulation results showed good agreement with experimental results.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2003.11a
• /
• pp.272-280
• /
• 2003

A definite shut-in pressure in hydraulic fracturing techniques is needed for obtaining the correct information on the in-situ stress regimes in rock masses. The relation between the behaviour of hydraulically induced fractures and the condition of remote stress is considered to be major reasons of an ambiguous shut-in pressure in hydraulic fracturing pressure-time history curves. This paper describes the results of a series of numerical analyses carried out using UDEC(Universal Distinct Element Code, Itasca), which is based on the discrete element method, to compare several methods for determining the shut-in pressure during hydraulic fracturing. The fully coupling of hydraulic and mechanical analysis was applied, and the effects of four different discontinuity geometries in numerical modelling have been investigated for this purpose. The effects of different remote stress regimes and different physical properties on hydraulic fracture propagation have been also analyzed. Several methods for obtaining shut-in pressure from the ambiguous shut-in curves have been applied to all the numerical models. The graphical intersection methods, such as (P vs. t) method, (P vs. log(t)) method, (log(P) vs. log(t)) method, give smaller values of the shut-in pressure than the statistical method, (dP/dt vs. P). Care should be taken in selecting a method for shut-in pressure, because there can be existed a stress anomaly around the wellbore and fracturing from the wellbore by a constant flow rate may have a more complicate mechanism.

• Journal of Ocean Engineering and Technology
• /
• v.32 no.3
• /
• pp.177-183
• /
• 2018

This study simulated ice load and the motion response of a moored semi-submersible rig in pack-ice conditions using a finite element method. Ice flows of random size and shape were modeled, and interactions for ice-sea, ice-structure, ice-ice were simulated using a simplified method. Parameters for the simplified method such as drag force coefficient and the pressure-penetration relation were obtained based on the result of detailed analysis using the coupled Eulerian-Lagrangian method. The mooring lines were modeled by spring elements based on their stiffness. As a result of the simulation over 1,400 seconds, the force and motion response of the rig were obtained and validated using discrete elements and compared with the results found by the Krylov State Research Centre.

• Computers and Concrete
• /
• v.21 no.5
• /
• pp.593-605
• /
• 2018

A nonlinear finite element modeling approach is developed to assess the behavior of a dowel bar embedded on a single concrete block substrate, subjected to monotonic loading. In this approach, a discrete representation of the steel reinforcing bar is considered, using beam finite elements with nonlinear material behavior. The bar is connected to the concrete embedment through nonlinear Winkler spring elements. This modeling approach can only be used if a new constitutive model is developed for the spring elements, to simulate the deformability and strength of the concrete substrate. To define this constitutive model, an extensive literature review was conducted, as well as 3 experimental tests, in order to select the experimental data which can be used in the calibration of the model. Based on this data, an empirical model was established to predict the global dowel response, for a wide range of bar diameters and concrete strengths. This empirical model provided the information needed for calibration of the nonlinear Winkler spring model, valid for dowel displacements up to 4 mm. This new constitutive model is composed by 5 stages, in order to reproduce the concrete substrate response.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.20 no.10
• /
• pp.923-930
• /
• 2010

Finite element models of dynamic systems can be updated in two stages. In the first stage, mass and stiffness matrices are updated neglecting damping, and in the second stage, damping matrices are estimated with the mass and stiffness matrices fixed. Three methods to estimate damping matrices for this purpose are proposed in this paper. The methods include one for proportional damping systems and two for non-proportional damping systems. Method 1 utilizes orthogonality of normal modes and estimates damping matrices using the modal parameters extracted from the measured responses. Method 2 estimates damping matrices from impedance matrices which are the inverse of FRF matrices. Method 3 estimates damping using the equation which relates a damping matrix to the difference between the analytical and measured FRFs. The characteristics of the three methods are investigated by applying them to simulated discrete system data and experimental cantilever beam data.