• Interaction and multiscale mechanics
• /
• v.1 no.1
• /
• pp.1-31
• /
• 2008

The present paper presents multiscale modelling via coupling of the discrete and finite element methods. Theoretical formulation of the discrete element method using spherical or cylindrical particles has been briefly reviewed. Basic equations of the finite element method using the explicit time integration have been given. The micr-macro transition for the discrete element method has been discussed. Theoretical formulations for macroscopic stress and strain tensors have been given. Determination of macroscopic constitutive properties using dimensionless micro-macro relationships has been proposed. The formulation of the multiscale DEM/FEM model employing the DEM and FEM in different subdomains of the same body has been presented. The coupling allows the use of partially overlapping DEM and FEM subdomains. The overlap zone in the two coupling algorithms is introduced in order to provide a smooth transition from one discretization method to the other. Coupling between the DEM and FEM subdomains is provided by additional kinematic constraints imposed by means of either the Lagrange multipliers or penalty function method. The coupled DEM/FEM formulation has been implemented in the authors' own numerical program. Good performance of the numerical algorithms has been demonstrated in a number of examples.

• Geomechanics and Engineering
• /
• v.7 no.6
• /
• pp.613-631
• /
• 2014

To complement the method of field-scale seismic ground motion simulations by buried blast techniques, the application and evaluation of the capability of a numerical modeling platform to simulate buried explosion-induced ground motion at a real soil site is presented in this paper. Upon a layout of the experimental setup at a level site wherein multiple charges that were buried over a large-diameter circle and detonated in a planned sequence, the formulation of a numerical model of the soil and the explosives using the finite element code LS-DYNA is developed for the evaluation of the resulting ground motion and surface subsidence. With a compact elastoplastic cap model calibrated for the loess soils on the basis of the site and laboratory test program, numerical solutions are obtained by explicit time integration for various dynamic aspects and their relation with the field blast experiment. Quantitative comparison of the computed ground acceleration time histories at different locations and induced spatial subsidence on the surface afterwards is given for further engineering insights in regard to the capabilities and limitations of both the numerical and experimental approaches.

• Interaction and multiscale mechanics
• /
• v.2 no.2
• /
• pp.129-151
• /
• 2009

This paper presents a meshfree procedure using a convex generalized meshfree (GMF) approximation for the large deformation analysis of particle-reinforced rubber compounds on microscopic level. The convex GMF approximation possesses the weak-Kronecker-delta property that guarantees the continuity of displacement across the material interface in the rubber compounds. The convex approximation also ensures the positive mass in the discrete system and is less sensitive to the meshfree nodal support size and integration order effects. In this study, the convex approximation is generated in the GMF method by choosing the positive and monotonic increasing basis function. In order to impose the periodic boundary condition in the unit cell method for the microscopic analysis, a singular kernel is introduced on the periodic boundary nodes in the construction of GMF approximation. The periodic boundary condition is solved by the transformation method in both explicit and implicit analyses. To simulate the interface de-bonding phenomena in the rubber compound, the cohesive interface element method is employed in corporation with meshfree method in this study. Several numerical examples are presented to demonstrate the effectiveness of the proposed numerical procedure in the large deformation analysis.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2000.04a
• /
• pp.195-204
• /
• 2000

In this paper a complicated structural behavior in collision and its effect of energy translation to the collision bulkhead was examined through a methodology of the numerical simulation to obtain a ideal bow construction and a location of collision bulkhead against head on collision. In the present the bow structure is normally designed in consideration of its specific structural arrangements and internal and external loads in these area such as hydrostatic and dynamic pressure, wave impact and bottom slamming in accordance with the Classification rules, and the specific location of collision bulkhead by SOLAS requirement. By these studies the behavior of the bow collapse due to collision was synthetically evaluated for the different size of tankers and its operational speed limits, and by the result of these simulation it provides the optimal design concept for the low construction to prevent the subsequent plastic deformation onto or near to the collision bulkhead boundary and to determine the rational location of collision bulkhead.

• 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 A
• /
• v.21 no.10
• /
• pp.1589-1597
• /
• 1997

The high-velocity impact forging process with eccentric loading condition is analyzed using the explicit time integration finite element method. In order to consider the strain hardening, strain rate hardening and thermal softening effects, which are frequently observed in high-velocity deformation phenomena, the Johnson-Cook constitutive model is applied to model the workpiece. It is assumed that the material response of the dies is elastic in the study. As a result of the eccentric loading simulation, it is found that the increase of the eccentric ratio and the allowable tilting angle cause the decrease of the maximum forging load and the blow efficiency, and it is also found that the forging load and the blow efficiency generated in the high-velocity impact forging process with three-dimensional geometry can be obtained efficiently.

• Structural Engineering and Mechanics
• /
• v.16 no.2
• /
• pp.153-176
• /
• 2003

From the equation of motion of a "bare" non-uniform beam (without any concentrated elements), an eigenfunction in term of four unknown integration constants can be obtained. When the last eigenfunction is substituted into the three compatible equations, one force-equilibrium equation, one governing equation for each attaching point of the concentrated element, and the boundary equations for the two ends of the beam, a matrix equation of the form [B]{C} = {0} is obtained. The solution of |B| = 0 (where ${\mid}{\cdot}{\mid}$ denotes a determinant) will give the "exact" natural frequencies of the "constrained" beam (carrying any number of point masses or/and concentrated springs) and the substitution of each corresponding values of {C} into the associated eigenfunction for each attaching point will determine the corresponding mode shapes. Since the order of [B] is 4n + 4, where n is the total number of point masses and concentrated springs, the "explicit" mathematical expression for the existing approach becomes lengthily intractable if n > 2. The "numerical assembly method"(NAM) introduced in this paper aims at improving the last drawback of the existing approach. The "exact"solutions in this paper refer to the numerical results obtained from the "continuum" models for the classical analytical approaches rather than from the "discretized" ones for the conventional finite element methods.

• Structural Engineering and Mechanics
• /
• v.36 no.3
• /
• pp.301-319
• /
• 2010

The paper seeks to explore some aspects of the current state of knowledge on progressive collapse in the technical literature covering blast loads and structural analysis procedure applicable to reinforced concrete (RC) buildings. The paper describes the progressive collapse analysis of a commercial RC building located in the city of Riyadh and subjected to different blast scenarios. A 3-D finite element model of the structure was created using LS-DYNA, which uses explicit time integration algorithms for solution. Blast loads were treated as dynamic pressure-time history curves applied to the exterior elements. The inherent shortcomings of notional member removal have been taken care of in the present paper by simulating the damage of structural elements through the use of solid elements with the provision of element erosion. Effects of erosion and cratering are studied for different scenarios of the blast.

• English Language & Literature Teaching
• /
• v.13 no.4
• /
• pp.87-108
• /
• 2007

Given that various classroom contextual factors influence the nature of writing instructional practices, it would be worthwhile to explore these factors to generate better environment for learning to write. Among many factors, this study examined teachers' beliefs and knowledge, which would operate as a very influential contextual factor in that changes in principles and methods of teaching writing would be the results of their underlying beliefs and knowledge related to teaching writing. Three professional teachers who teach second- and third-grade English language learners (ELLs) were interviewed, and the analysis of teacher interviews was conducted. The research findings indicated that basically all of the teachers perceived the role of writing in second language learning as very important, sharing the belief that the ultimate goal of teaching writing is to have their students gain fluency in writing and that some of instructional methods such as integration of writing and other language aspects, content-based writing, and providing scaffolding are important. In addition, some beliefs that two ESL teachers shared included the importance of ample and continuous opportunities to write, vocabulary knowledge, and explicit instruction about writing. Other beliefs, including the importance of creating a comfortable writing environment and opportunities for writing for varied purposes and genres were represented.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.14 no.5
• /
• pp.424-431
• /
• 2004

An implicit direct-time integration method for obtaining transient responses of general dynamic systems is described. The conventional Newmark method cannot be directly applied to state-space first-order differential equations, which contain no explicit acceleration terms. The method proposed here is the state-space Newmark method that incorporates the average velocity concept, and can be applied to an analysis of general dynamic systems that are expressed by state-space first-order differential equations. It is also readily coded into a program. Stability and accuracy analyses indicate that the method is numerically unconditionally stable like the conventional Newmark method, and has a period error of 2nd-order accuracy for small damping and 4th-order for large damping and an amplitude error of 2nd-order, regardless of damping. In addition, its utility and validity are confirmed by two application examples. The results suggest that the proposed state-space Newmark method based on average velocity be generally applied to the analysis of transient responses of general dynamic systems with a high degree of reliability with respect to stability and accuracy.