• Tunnel and Underground Space
• /
• v.32 no.6
• /
• pp.598-610
• /
• 2022

Avoiding a fault zone would be a best practice for safety in underground construction, which is only sometimes possible because of many restrictions and other field conditions. For instance, there is an ongoing conception of Korea-Japan subsea tunnels that inevitably cross a massive fault system in the Korea Strait. Therefore it was deemed necessary to find an efficient way of predicting the likely behaviour of underground structures under fault slip. This paper presents the findings from simple numerical analysis for investigating the stress induced at a normal fault with a dip of 45 degrees. We used a boundary element software that assumed constant displacement discontinuity, which allowed the displacement to be estimated separately at both the fault's hangingwall and footwall sides. The results suggested that a principal stress rotation of 45 degrees occurred at the edges of the fault during the slip, which was in agreement with the phenomenon for fault plane suggested in the body of literature. A simple numerical procedure presented in this paper could be adopted to investigate other fault-related issues associated with underground structure construction.

• Computers and Concrete
• /
• v.23 no.3
• /
• pp.161-170
• /
• 2019

This paper aims at presenting a numerical method for estimating the elastic modulus of recycled concrete with crushed aggregates. In the method, polygonal aggregates following a given sieve curve are generated, and placed into a square simulation element with the aid of the periodic boundary condition and the overlap criterion of two polygonal aggregates. The mesostructure of recycled concrete is reconstructed by embedding an old interfacial transition zone (ITZ) layer inside each recycled aggregate and by coating all the aggregates with a new ITZ layer. The square simulation element is discretized into a regular grid and a representative point is selected from each sub-element. The iterative method is combined with the fast Fourier transform to evaluate the elastic modulus of recycled concrete. After the validity of the numerical method is verified with experimental results, a sensitivity analysis is conducted to evaluate the effects of key factors on the elastic modulus of recycled concrete. Numerical results show that the elastic modulus of recycled concrete increases with the increase of the total aggregate content and the elastic moduli of old and new ITZ but decreases with increasing the replacement ratio of recycled aggregate and the thicknesses of old and new ITZ. It is also shown that, for a replacement ratio of recycled aggregate smaller than 0.3, the elastic modulus of recycled concrete is reduced by no more than 10%.

• Journal of Biomedical Engineering Research
• /
• v.21 no.4
• /
• pp.373-384
• /
• 2000

The goal of this paper is to develop a computational method to predict cancellous bone density distributions based upon continuum levels of volumetric strain. Volumetric strain is defined as the summation of normal strains, excluding shear strains, within an elastic range of loadings. Volumetric strain at a particular location in a cancellous structure changes with changes of the boundary conditions (prescribed displacements, tractions, and pressure). This change in the volumetric strain is postulated to predict the adaptive change in the bone apparent density. This bone remodeling theory based on volumetric strain is then used with the finite element method to compute the apparent density distribution for cancellous bone in both lumbar spine and proximal femur using an iterative algorithm, considering the dead zone of strain stimuli. The apparent density distribution of cancellous bone predicted by this method has the same pattern as experimental data reported in the literature (Wolff 1892, Keller et al. 1989, Cody et al. 1992). The resulting bone apparent density distributions predict Young's modulus and strength distributions throughout cancellous bone in agreement with the literature (Keller et al. 1989, Carter and Hayes 1977). The method was convergent and sensitive to changes in boundary conditions. Therefore, the computational algorithm of the present study appears to be a useful approach to predict the apparent density distribution of cancellous bone (i.e. a numerical approximation for Wolff's Law)

• Korean Journal of Remote Sensing
• /
• v.17 no.1
• /
• pp.1-14
• /
• 2001

The enhancement of the refractive index structure parameter $C_n^2$ often occurs where vertical gradients of virtual potential temperature ${\theta}_v$ and mixing ratio q have their maximum values. The $C_n^2$ can be a very useful parameter for estimating the convective boundary layer(CBL) height. The behavior of $C_n^2$ peaks, often used to locate the height of mixed layer, was investigated in the present study. In addition, a new method to determine the CBL height objectively using both $C_n^2$ and vertical air velocity variance ${\sigma}_w$ data of UHF radar was also suggested. The present analysis showed that the $C_n^2$ peaks in the backscatter intensity profiles often occurred not only at the top of the CBL but also at the top of a residual layer or at a cloud layer. The $C_n^2$ peaks corresponding to the CBL heights were slightly lower than the CBL heights derived from rawinsonde sounding data when vertical mixing owing to weak solar heating was not significant and the height of strong vertical ${\theta}_v$ gradients were not consistent with that of strong vertical q gradients. However, the $C_n^2$ peaks corresponding to the CBL heights were in good agreement with the rawinsonde-estimated CBL hegiths when vertical mixing owing to solar heating was significant and the vertical gradient of both ${\theta}_v$ and q in the entrainment zone was very strong. The maximum backscatter intensity method, which determines the height of $C_n^2$ peak as the CBL height, correctly estimated the CBL height when the $C_n^2$ profile had single peak, but this method erroneously estimated the CBL height when there was a residual layer or a cloud layer over the top of the CBL. The new method distinguished when there the CBL height from the peak due a cloud layer or a residual layer using both $C_n^2$ and ${\sigma}_w$ data, and correctly estimated the CBL height. As for estimation of diurnal variation of the CBL height, the new method backscatter intensity method even if the vertical profile of backscatter intensity had two peaks from the CBL height and a residual layer or a cloud layer.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.21 no.4
• /
• pp.551-561
• /
• 2018

Footprint is defined as ground area that is projected from the outer edges of the battle space protected by a defence system. This concept can be effectively used for making decisions on site selection of anti missile systems to defend against enemy's ballistic missiles. In this paper, simulations of ballistic missile trajectories based on various launch conditions are performed first and then the footprint is derived with engagement zone set as a boundary condition. Results of the simulation with various relative positions between the defense system and defended asset are also presented. The proposed method, in which the trajectories are generated based on launch point of the ballistic missile, has an advantage of approximating the defended area close to reality. Two applications are introduced in the present paper to describe how the derivation of defended area could be utilized in deployment decision of defense systems.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.7
• /
• pp.997-1004
• /
• 2001

Burner of Flat Flame type expects the uniform flame distribution and NOx reduction. The characteristics of Flat Flame Burner become different according to swirl number in the burner throat. Experiments were focused on swirl effect by four types of swirler with different swirl numbers (0, 0.26, 0.6 and 1.24). It shows many different flow patterns according to swirl number using PIV(Particle Image Velocimetry) method. The flow of burner with swirler is recirculated by pressure difference between its center and outside. Recirculated air makes stable in flame, and reduced pollutant gas. In case of swirl number 0, main flow passes through axial direction. As swirl number increased, The backward flow develops in the center part of burner and Flow gas recirculates. This is caused by radial flow momentum becomes larger than axial flow by swirled air and the pressure at center drops against surrounding. As swirl number increases, the radial and axial velocity was confirmed to be larger than low swirl numbers. And turbulence intensity have similar pattern. The CTRZ(Central Toroidal Recirculation Zone) is shown evidently when y/D=1 and S=1.24. The boundary-layer between main flow and recirculated flow is shown that the width is seen to be decreased as swirl number increased.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.7 s.178
• /
• pp.1650-1660
• /
• 2000

Continuous casting process has been adopted increasingly in recent years to save both energy and labor. It has experienced a rapid development in the production of semi-finished steel products, rep lacing the conventional route of ingot casting and rolling. In order to achieve this merit, however, more studies about the mechanism between roll and slab are needed. In this paper, a dynamic bulging in steel cast slabs was simulated by considering the solidification and heat transfer. This study is to prevent internal cracks of a slab in a bending and unbending zone. The value of moving strand shell bulging between two supporting rollers under the ferrostatic pressure and slab-self weight has been calculated in terms of creep and elasto-plasticity. The strain and strain rate distributions in solidified shell undergoing a series of bulging are calculated with working boundary conditions.

• Transactions of Materials Processing
• /
• v.33 no.3
• /
• pp.193-199
• /
• 2024

The friction welding characteristics of stainless steels, mainly used in energy and chemical plant industries due to its excellent corrosion resistance and high strength, was evaluated in this study. Friction welding was introduced and conducted at a rotation speed of 2,000 RPM, friction pressure of 30 MPa, burn-off length of 5 mm and upset pressure of 110 ~ 200 MPa on rod typed specimens. The grain boundary characteristics distributions such a grain size, shape, misorientation angle and kernel average misorientation of the welds were clarified by electron backscattering diffraction method. The application of friction welding on SUS304 alloy resulted in a significant refinement of the grain size in the weld zone (5.11 mm) compared to that of the base material (48.09 mm). The mechanical properties of the welds, on the other hand, appeared to be relatively low or similar to those of the base material, which were mainly caused by dislocation density in the initial material and grain refinement in the welds.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2005.06a
• /
• pp.965-968
• /
• 2005

Roll drafting operation causes variations in the linear density of bundles because the bundle flow cannot be controlled completely by roll pairs. Defects occurring in this operation bring about many problems successively in the next processes. In this paper, we attempt to analyze the draft dynamics and the linear density irregularity based on the governing equation of a bundle motion that has been suggested in our previous studies. For analyzing the dynamic characteristics of the roll drafting operation, it is indispensable to investigate a transient state in time domain before the bundle flux reaches a steady state. However, since governing equations of bundle flow consisting of continuity and motion equations turn out to be nonlinear, and coupled between variables, the solutions for a transient state cannot be obtained by an analytical method. Therefore, we use the Finite Difference Method(FDM), particularly, the FTBS(Forward-Time Backward-Space) difference method. Then, the total equations system yields to an algebraic equations system and is solved under given initial and boundary conditions in an iterative fashion. From the simulation results, we confirm that state variables show different behavior in the transient state; e.g., the velocity distribution in the flow field changes more quickly the linear density distribution. During a transient flow in a drafting zone, the output irregularity is influenced differently by the disturbances, e.g., the variation in input bundle thickness, the drafting speed, and the draft ratio.

• Tunnel and Underground Space
• /
• v.32 no.5
• /
• pp.312-326
• /
• 2022

The generalized Hoek-Brown (GHB) criterion, which is recognized as one of the standard failure conditions for rock mass, is specialized for rock engineering applications and covers a wide range of rock mass conditions. Accordingly, many research efforts have been devoted to the incorporation of this criterion into the stability analysis of rock structures. In this study, the slip-line analysis method, which is a kind of elastoplastic analysis method, is combined with the GHB failure criterion to derive analytical equations that can easily calculate the plastic radius and stress distribution in the vicinity of the circular tunnel. In the process of derivation of related formulas, it is assumed that the behavior of rock mass after failure is perfectly plastic and the in-situ stress condition is hydrostatic. In the formulation, it is revealed that the plastic radius can be calculated analytically using the two respective tangential friction angles corresponding to the stress conditions at tunnel wall and elastic-plastic boundary. It is also shown that the plastic radius and stress distribution calculated using the derived analytical equations coincide with the results of Lee & Pietruszczak's numerical method published in 2008. In the latter part of this paper, the influence of the quality of the rock mass on the size of the plastic zone, the stress distribution, and the change of the tangential friction angle was investigated using the derived analytical equations.