• Structural Engineering and Mechanics
• /
• v.60 no.5
• /
• pp.809-828
• /
• 2016

In deep open pit mines, slope stability is very important. Particularly, increasing the depths increase the risks in mines having weak rock mass. Blasting operations in this type of open pits may have a negative impact on slope stability. Several or combination of methods can be used in order to enable better analysis in this type of deep open-pit mines. Numerical modeling is one of these options. Many complex problems can be integrated into numerical methods at the same time and analysis, solutions can be performed on a single model. Rock failure criterions and rock models are used in numerical modeling. Hoek-Brown and Mohr-Coulomb terms are the two most commonly used rock failure conditions. In this study, mine planning and discontinuity conditions of a lignite mine facing two big landslides previously, has been investigated. Moreover, the presence of some damage before starting the study was identified in surrounding structures. The primary research of this study is on slope study. In slope stability analysis, numerical modeling methods with Hoek-Brown and Mohr-Coulomb failure criterions were used separately. Preparing the input data to the numerical model, the outcomes of patented-blast vibration minimization method, developed by co-author was used. The analysis showed that, the model prepared by applying Hoek-Brown failure criterion, failed in the stage of 10. However, the model prepared by using Mohr-Coulomb failure criterion did not fail even in the stage 17. Examining the full research field, there has been ongoing production in this mine without any failure and damage to surface structures.

• The Journal of Engineering Geology
• /
• v.19 no.2
• /
• pp.205-215
• /
• 2009

The conventional numerical models to analyze flow in subsurface porous media under the transient state usually generate numerical oscillation and unstability due to local flux domain for critical cases such as infiltration into initially dry soil during rainfall period. In this case, it is required refined mesh and small time step, but it decrease efficiency of computation. In this study, numerical unstability in discontinuity domain is removed by applying particle tracking algorithm to simulate unsteady subsurface flow with inflow boundary condition. Finally the hybrid LE FEM improving numerical stability is proposed. The hypothetical domains with unsteady uniform and nonuniform flow field were used to demonstrated algorithm verification. In comparison with analytic solution, we obtained reasonable results and conducted simulation of hypothetical 3-D recharge/pumping area. The proposed algorithm can simulate saturated/unsaturated porous media with more practical problems and will greatly contribute to accuracy and stability of numerical computation.

• Korean Journal of Computational Design and Engineering
• /
• v.13 no.4
• /
• pp.296-304
• /
• 2008

Recently, various attempts are being done to apply off-line programming system to field of paint robot. But most commercial simulation softwares have problems that are slow simulation speed and not support various painting paramenters on simulation. This paper proposes enhanced paint simulation method for off-line programming system. For these, this method used the mathematical model of flux field from a previous research. The flux field has the flux distribution function, which reflects on the feature of paint spray. A previous research derived this flux distribution function for an integral function and calculated paint thickness function for an integral function. But if flux distribution function is defined as an integral function, it is inadequate to use for real-time simulation because a number of calculation is needed for estimation of paint thickness distribution. Therefore, we defined the flux distribution function by numerical method for reducing a mount of calculation for estimation of paint thickness. We derived the equation of paint thickness function analytically for reducing a mount of calculation from the paint distribution function defined by numerical method. In order to prove proposed paint simulation method this paper compares the simulated and measured thickness. From this comparison this paper show that paint thickness distribution is predicted precisely by proposed spray paint simulation process.

• Structural Engineering and Mechanics
• /
• v.41 no.2
• /
• pp.157-181
• /
• 2012

The paper presents a review of the application of the newly proposed mixed finite element model for seismic simulation of different types of composite frame structures. To evaluate the performance of the element, a comparison with displacement-based and force-based models is conducted. The study revealed that the mixed model is superior to the others in terms of both speed of convergence and numerical stability, and is therefore considered the most practical approach for modeling of composite structures. In this model, the element is derived using independent force and displacement shape functions. The nonlinear response of the frame element is based on the section discretization into fibers with uniaxial material models. The interfacial behavior is modeled using an inelastic interface element. Numerical examples to clarify the advantages of the model are presented for the following structural applications: anchored reinforcing bar problems, composite steel-concrete girders with deformable shear connectors, beam on elastic foundation elements, R/C girders strengthened with FRP sheets, R/C beam-columns with bond-slip, and prestressed concrete girders. These studies confirmed that the model represents a major advancement over existing elements in simulating the inelastic behavior of composite structures.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.12 no.3
• /
• pp.465-474
• /
• 1999

A recently developed numerical method based on a volume integral formulation is applied to calculate the accurate stress intensity factors at the crack tips in unbounded isotropic solids in the presence of multiple anisotropic inclusions and cracks subject to external loads. In this paper, a detailed analysis of the stress intensity factors are carried out for an unbounded isotropic matrix containing an orthotropic cylindrical inclusion and a crack. The accuracy and effectiveness of the new method are examined through comparison with results obtained from analytical method and finite element method using ANSYS. It is demonstrated that this new method is very accurate and effective for solving plane elastostatic problems in unbounded solids containing anisotropic inclusions and cracks.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2006.11a
• /
• pp.934-937
• /
• 2006

This investigation presents a topology formulation to design optimal poroelastic acoustic foams to maximize absorbing ability. For successful formulation, a single set of equations based on Biot's theory is adopted and an appropriate material interpolation strategy is newly developed. Because there was no earlier attempt to solve poroelastic acoustic foam design problems in topology optimization setting, many challenging issues including modeling and interpolation must be addressed. First, the simulation accuracy by a proposed unified model encompassing acoustic air and poroelastic material was checked against analytical and numerical results. Then a material interpolation scheme yielding a distinct acoustic air-poroelastic material distribution was developed. Using the proposed model and interpolation scheme, the topology optimization of a two-dimensional poroelastic acoustic foam for maximizing its absorption coefficient was carried out. Numerical results show that the absorption capacity of an optimized foam layout considerably increases in comparison with a nominal foam layout.

• Proceedings of the KSME Conference
• /
• 2007.05a
• /
• pp.1663-1666
• /
• 2007

The finite element method is applied to analyze the deformation mechanisms in the closed-cell Al foam under the compression. The modeling of the real cellular structure proceeds with the concept of the reverse engineering. First of all, the small, $10{\times}\;10{\times}\;10mm^3$ sized specimens of the closed-cell Al foam are prepared. The micro focus X-ray CTsystem of SHIMADZU Corp. is used to scan the full structures of the specimens. The scanned structures are converted to the geometric surfaces and solids through the software for 3-D scan data processing, RapidFormTMof INUS Tech. Inc. Then the solid meshes are directly generated on the converted geometric solids for the finite element analysis. The large elastic-plastic deformation and 3-D contact problems for the Al cellular material are considered. The clear and successful analysis for the deformation mechanisms in the closed-cell Al foam is carried out through the comparison of the numerical results in this research with the referred experimental ones.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.7
• /
• pp.1426-1437
• /
• 2002

Using the Lagrangian explicit time-integration finite element code NET3D which can treat three-dimensional high-velocity impact problems, oblique penetration processes of long rod projectile with yaw against thin plate are simulated. Through the comparison of simulation result with experimental result and other code's computational result, the adaptability and accuracy of NET3D is evaluated under the complex situation in which yaw angle and oblique angle exist simultaneously. Main research contents to be handled in this paper include the followings. First, the accuracy and efficiency estimation of NET3D code result obtained from the oblique penetration simulations of long rod projectile with yaw against thin plate. Second, the effect of increasing impact velocity. Third, the effect of initial yaw for the spaced-plate target. Residual velocities, residual lengths, angular velocities, and final deformed configurations obtained from the NET3D computations are compared with the experimental results and other code's computational results such as Eulerian code MESA and Lagrangian code EPIC. As a result of comparisons, it has been found that NET3D code is superior to EPIC code and MESA code in the prediction capability of residual velocity and residual length of penetrator. The key features obtained from the experiment can be successfully reproduced through NET3D simulations. Throughout the study, the applicability and accuracy of NET3D as a metallic armor system design tool is verified.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.5 s.236
• /
• pp.609-616
• /
• 2005

While the SIMPLE algorithm is most widely used for the simulations of flow phenomena that take place in the industrial equipments or the manufacturing processes, it is less adopted for the simulations of the free surface flow. Though the SIMPLE algorithm is free from the limitation of time step, the free surface behavior imposes the restriction on the time step. As a result, the explicit schemes are faster than the implicit scheme in terms of computation time when the same time step is applied to, since the implicit scheme includes the numerical method to solve the simultaneous equations in its procedure. If the computation time of SIMPLE algorithm can be reduced when it is applied to the unsteady free surface flow problems, the calculation can be carried out in the more stable way and, in the design process, the process variables can be controlled based on the more accurate data base. In this study, a modified SIMPLE algorithm is presented fur the free surface flow. The broken water column problem is adopted for the validation of the modified algorithm (MoSIMPLE) and for comparison to the conventional SIMPLE algorithm.

• Journal of computational fluids engineering
• /
• v.5 no.1
• /
• pp.1-7
• /
• 2000

In general, FDM(finite difference method) and FVM(finite volume method) are used for analyzing the fluid flow numerically. However it is difficult to apply them to problems involving complex geometries, multi-connected domains, and complex boundary conditions. On the contrary, FEM(finite element method) with coordinates transformation for the unstructured grid is effective for the complex geometries. Most of previous studies have used commercial codes such as KIVA or STAR-CD for the flow analyses in the engine cylinder, and these codes are mostly based on the FVM. In the present study, using the FEM for three-dimensional, unsteady, and incompressible Navier-Stokes equation, the velocity and pressure fields in the engine cylinder have been numerically analyzed. As a numerical algorithm, 4-step time-splitting method is used and ALE(arbitrary Lagrangian Eulerian) method is adopted for moving grids. In the Piston-Cylinder, the calculated results show good agreement in comparison with those by the FVM and the experimental results by the LDA.