• Proceedings of the Korean Geotechical Society Conference
• /
• 2000.11a
• /
• pp.705-712
• /
• 2000

The fluid generally flows through fractures in crystalline rocks where most of underground storage facilities are constructed because of their low hydraulic conductivities. The fractured rock is better to be conceptualized with a discrete fracture concept, rather continuum approach. In the aspect of fluid flow in underground, the simultaneous flow of groundwater and gas should be considered in the cases of generation and leakage of gas in nuclear waste disposal facilities, air sparging process and soil vapor extraction for eliminating contaminants in soil or rock pore, and pneumatic fracturing for the improvement of permeability of rock mass. For the purpose of appropriate analysis of groundwater-gas flow, this study presents an unsteady-state multi-phase FEM fracture network simulator. Numerical simulation has been also conducted to investigate the hydraulic head distribution and air tightness around Ulsan LPG storage cavern. The recorded hydraulic head at the observation well Y was -5 to -10 m. From the results obtained by the developed model, it shows that the discrete fracture model yielded hydraulic head of -10 m, whereas great discrepancy with the field data was observed in the case of equivalent continuum modeling. The air tightness of individual fractures around cavern was examined according to two different operating pressures and as a result, only several numbers of fractures neighboring the cavern did not satisfy the criteria of air tightness at 882 kPa of cavern pressure. In the meantime, when operating pressure is 710.5 kPa, the most areas did not satisfy air tightness criteria. Finally, in the case of gas leaking from cavern to the surrounding rocks, the resulted hydraulic head and flowing pattern was changed and, therefore, gas was leaked out from the cavern ceiling and groundwater was flowed into the cavern through the walls.

• Steel and Composite Structures
• /
• v.20 no.3
• /
• pp.623-649
• /
• 2016

Most of the early studies on plates vibration are focused on two-dimensional theories, these theories reduce the dimensions of problems from three to two by introducing some assumptions in mathematical modeling leading to simpler expressions and derivation of solutions. However, these simplifications inherently bring errors and therefore may lead to unreliable results for relatively thick plates. The main objective of this research paper is to present 3-D elasticity solution for free vibration analysis of continuously graded carbon nanotube-reinforced (CGCNTR) rectangular plates resting on two-parameter elastic foundations. The volume fractions of oriented, straight single-walled carbon nanotubes (SWCNTs) are assumed to be graded in the thickness direction. In this study, an equivalent continuum model based on the Eshelby-Mori-Tanaka approach is employed to estimate the effective constitutive law of the elastic isotropic medium (matrix) with oriented, straight carbon nanotubes (CNTs). The proposed rectangular plates have two opposite edges simply supported, while all possible combinations of free, simply supported and clamped boundary conditions are applied to the other two edges. The formulations are based on the three-dimensional elasticity theory. A semi-analytical approach composed of differential quadrature method (DQM) and series solution is adopted to solve the equations of motion. The fast rate of convergence of the method is demonstrated and comparison studies are carried out to establish its very high accuracy and versatility. The 2-D differential quadrature method as an efficient and accurate numerical tool is used to discretize the governing equations and to implement the boundary conditions. The convergence of the method is demonstrated and to validate the results, comparisons are made between the present results and results reported by well-known references for special cases treated before, have confirmed accuracy and efficiency of the present approach. The novelty of the present work is to exploit Eshelby-Mori-Tanaka approach in order to reveal the impacts of the volume fractions of oriented CNTs, different CNTs distributions, various coefficients of foundation and different combinations of free, simply supported and clamped boundary conditions on the vibrational characteristics of CGCNTR rectangular plates. The new results can be used as benchmark solutions for future researches.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2005.05a
• /
• pp.122-125
• /
• 2005

The AutoForm previously used the membrane element and it accomplished sheet metal forming analysis. The membrane analysis has been widely applied to various sheet metal forming processes because of its saving time effectiveness. However, it's well known that the membrane analysis can not provides correct information for the processes which considerable bending effects. From this time research it tried to compare the formation analysis result which uses the shell element which is applied newly in the AutoForm and actual products. The shell element is compromise method between continuum analysis and membrane analysis. The Finite element method by using shell element is the most economical numerical method. From analysis results, FEA by using shell element can estimate accurately the problems happened in actual auto-body panel.

• The Journal of Engineering Geology
• /
• v.15 no.3
• /
• pp.257-268
• /
• 2005

In situ rock mass is generally heterogeneous and discontinuous, with varying degrees of strength along the planes of weakness. The planes of weakness such as joints, faults, cracks and bedding planes, control the strength and deformation characteristics of the rock mass. Subsequently, the stability of underground opening depends upon the spatial distribution of discontinuities and their mechanical properties in relation with geometrical shape of openins as well as the mechanical properties of intact rock materials. Understanding the behaviour of a discontinuous rock mass remains a key issue for improving excavation design in hiかy stressed environments. Although recent advances in rock mechanics have provided guidelines for the design of underground opening in isotropic rock mass, prediction and control of deformation in discontinuous rock masses are still unclear. In this study, parametric study was performed to investigate the plastic zone size, stress distribution and deformation behavior around underground opening in a discontinuous rock mass using a continuum joint model. The solutions were obtained by an elasto-plastic finite difference analysis, employing the Mohr-Coulomb failure criteria. Non-associated flow rule and perfectly plastic material behavior are also assumed.

• Journal of the Korean Geotechnical Society
• /
• v.19 no.4
• /
• pp.167-178
• /
• 2003

Damage and overbreak of the remaining rock induced by blasting can not be avoided during tunnel construction which may result in either short-term or long-term tunnel instability. Therefore, in this paper, a methodology to take into account the effect of blast-induced damage in tunnel stability assessment is proposed. Dynamic numerical analysis was executed to evaluate damage and overbreak of the remaining rock for the most common blasting pattern in road tunnel. Rock damage was quantified by utilizing the damage variable factor which is adopted proposed in continuum damage mechanics. The damaged rock stiffness and the damaged failure criteria are used to consider the effect of rock damage in tunnel stability analysis. The damaged geological strength index of the damaged rock was newly proposed from the relationship between deformation modulus and geological strength index. Also the Hoek-Brown failure criteria of the damaged rock was obtained using the damaged geological strength index. Analysing the tunnel stability with the consideration of the blast-induced damage of remaining rock, it was found that the extend of plastic zone and deformation increased compared to the undamaged rock. Therefore the short-term or long-term tunnel stability will be threatened when the rock damage from blasting is ignored in the tunnel stability analysis.

• Earthquakes and Structures
• /
• v.24 no.6
• /
• pp.393-404
• /
• 2023

This paper presents the application of hybrid-simulation-based adapter elements for the non-linear two-scale analysis of reinforced concrete frames with masonry infills under seismic-like demands. The approach provides communication and distribution of the computations carried out by two or more remote or locally distributed numerical models connected through the OpenFresco Framework. The modeling consists of a global analysis formed by macro-elements to represent frames and walls, and to reduce global degrees of freedom, portions of the structure that require advanced analysis are substituted by experimental elements and dimensional couplings acting as interfaces with their respective sub-assemblies. The local sub-assemblies are modeled by solid finite elements where the non-linear behavior of concrete matrix and masonry infill adopt a continuum damage representation and the reinforcement steel a discrete one, the conditions at interfaces between concrete and masonry are considered through a contact model. The methodology is illustrated through the analysis of a frame-wall system subjected to lateral loads comparing the results of using macro-elements, finite element model and experimental observations. Finally, to further assess and validate the methodology proposed, the paper presents the pushover analysis of two more complex structures applying both modeling scales to obtain their corresponding capacity curves.

• Geomechanics and Engineering
• /
• v.2 no.2
• /
• pp.125-140
• /
• 2010

Slope stability analysis of the right abutment of a railway bridge proposed at about 350 m above the ground level, crossing a river and connecting two huge hillocks in the Himalayas, India is presented in this paper. The site is located in a highly active seismic zone. The rock slopes are intensely jointed and the joint spacing and orientation are varying at different locations. Static slope stability of the rock slope is studied using equivalent continuum approach through the most commonly used commercial numerical tools like FLAC and SLOPE/W of GEOSTUDIO. The factor of safety for the slope under static conditions was 1.88 and it was reduced by 46% with the application of earthquake loads in pseudo-static analysis. The results obtained from the slope stability analyses confirmed the global stability of the slope. However, it is very likely that there could be possibility of wedge failures at some of the pier locations. This paper also presents the results from kinematics of right abutment slope for the wedge failure analysis based on stereographic projections. Based on the kinematics, it is recommended to flatten the slope from 50o to 43o to avoid wedge failures at all pier locations.

• Advanced Composite Materials
• /
• v.18 no.1
• /
• pp.77-93
• /
• 2009

This paper investigated the damage transition mechanism between the fiber-breaking mode and the fiber-avoiding crack mode when the fiber-length is reduced in the unidirectional discontinuous carbon fiber-reinforced-plastics (CFRP) composites. The critical fiber-length for the transition is a key parameter for the manufacturing of flexible and high-strength CFRP composites with thermoset resin, because below this limit, we cannot take full advantage of the superior strength properties of fibers. For this discussion, we presented a numerical model for the microscopic damage and fracture of unidirectional discontinuous fiber-reinforced plastics. The model addressed the microscopic damage generated in these composites; the matrix crack with continuum damage mechanics model and the fiber breakage with the Weibull model for fiber strengths. With this numerical model, the damage transition behavior was discussed when the fiber length was varied. The comparison revealed that the length of discontinuous fibers in composites influences the formation and growth of the cluster of fiber-end damage, which causes the damage mode transition. Since the composite strength is significantly reduced below the critical fiber-length for the transition to fiber-avoiding crack mode, we should understand the damage mode transition appropriately with the analysis on the cluster growth of fiber-end damage.

• Journal of Mechanical Science and Technology
• /
• v.14 no.2
• /
• pp.236-250
• /
• 2000

Compression ignition direct injection diesel engines employed a high pressure injection system have been developed as a measure to improve a fuel efficiency and reduce harmful emissions. In order to understand the effects of the pressure variation, many experimental works have been done, however there are many difficulties to get data in engine condition. This work gives numerical results for the high pressure effects on spray characteristics in wide or limited space with near walls. The gas phase is modelled by Eulerian continuum conservation equations of mass, momentum, energy and fuel vapour fraction. The liquid phase is modelled using the discrete droplet model approach in Lagrangian form and the drop behavior on a wall is calculated with a new droplet-wall interaction model based on the experiments observing individual drops. The droplet distributions, vapour fractions and gas flows are shown in various injection pressure cases. In free spray case which the injection spray has no wall impaction, the spray dispersion and vapour fraction increase and drop sizes decrease with increasing injection pressure. The same phenomena appears more clearly in wall impaction cases.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.625-630
• /
• 2000

Numerical and experimental investigations are peformed for the rarefied gas flows in pumping channels of a helical-type drag pump. Modern turbomolecular pumps include a drag stage in the discharge side, operating roughly in $10^{-2}{\sim}10Torr$. The flow occurring in the pumping channel develops from the molecular transition to slip flow traveling downstream. Two different numerical methods are used in this analysis: the first one is a continuum approach in solving the Navier-Stokes equations with slip boundary conditions, and the second one is a stochastic particle approach through the use of the direct simulation Monte Carlo(DSMC) method. The flow in a pumping channel is three-dimensional(3D), and the main difficulty in modeling a 3D case comes from the rotating frame of reference. Thus, trajectories of particles are no longer straight lines. In the Present DSMC method, trajectories of particles are calculated by integrating a system of differential equations including the Coriolis and centrifugal forces. Our study is the first instance to analyze the rarefied gas flows in rotating frame in the presence of noninertial effects.