• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2005.04a
• /
• pp.383-389
• /
• 2005

This research is to develop a seismic response analysis method for a spent fuel storage cask. FEM model is built for the test model of 1/8 scale spent fuel dry storage cask using available 3D contact conditions in ABAQUS/Explicit. Input load for this analysis os a seismic wave of El-centro earthquake, and the friction and damping coefficients in the analysis condition we obtained from the test result. Penalty and kinematic contact methods of ABAQUS are used for mechanical contact formulation. The analysis method was verified for rocking angle obtained by seismic response tests. The kinematic contact method with an adequate normal contact stiffness showed a good agreement with tests.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.116-120
• /
• 2009

A two-dimensional unsteady inviscid flow solver has been developed for the simulation of complex geometric configurations on adaptive Cartesian meshes. Embedded condition was used for boundary condition and a predictor-corrector explicit time marching scheme was used for time-accurate numerical simulation. The Cartesian mesh generator, which was previously developed for steady problem, was used grid generation for unsteady flow. The solver was based on ALE formulation for body motion. For diminishing the effects of cut-cells, the cell merging method was used. Using cell merging method, it was eliminated the CFL constraints. The conservation problem, which is caused cell-type variation around region swept by solid boundary, was also solved using cell merging method. The results are presented for 2D circular cylinder and missile launching problem.

• Journal of Mechanical Science and Technology
• /
• v.17 no.4
• /
• pp.545-561
• /
• 2003

We use the first-passage-time formulation by Torquato, Kim and Cule ［J. Appl. Phys., Vol. 85, pp. 1560∼1571 (1999) ］, which makes use of the first-passage region in association with the diffusion tracer's Brownian movement, and develop a new efficient Brownian motion simulation method to compute the effective conductivity of digitized composite media. By using the new method, one can remarkably enhance the speed of the Brownian walkers sampling the medium and thus reduce the computation time. In the new method, we specifically choose the first-passage regions such that they coincide with two, four, or eight digitizing units according to the dimensionality of the composite medium and the local configurations around the Brownian walkers. We first obtain explicit solutions for the relevant first-passage-time equations in two-and three-dimensions. We then apply the new method to solve the illustrative benchmark problem of estimating the effective conductivities of the checkerboard-shaped composite media. for both periodic and random configurations. Simulation results show that the new method can reduce the computation time about by an order of magnitude.

• Proceedings of the KIEE Conference
• /
• 1991.07a
• /
• pp.417-421
• /
• 1991

In electromagnetic transient analysis, complex transmission system should be modelled in detail. But in large system, this full representation of power transmission system has a big burden in many sides such as computation time, modelling efforts, etc. It is very required, therefore, in electromagnetic transients studies to represent parts of a complete system in a reduced or an equivalent form. This paper develops the method from which system equivalents may be derived. The suggested method is of an essentially transient form, and allows travelling wave interaction between the equivalent and explicit transmission network to be modelled.

• Proceedings of the Korea Concrete Institute Conference
• /
• 2002.05a
• /
• pp.303-308
• /
• 2002

A governing differential equation (GDE) of PSC flexural member with constant eccentricity considering the long-term losses including concrete creep, shrinkage, and PS steel relaxation is derived based on the two approaches. The first approach utilizes the force and moment equilibrium equations derived based on the geometry of strains of the uniform and curvature strains while the second one utilizes the principle of minimum total potential energy formulation. The identity of the two GDE's is verified by comparing the coefficients consisting of the GDE's. The boundary conditions resulting from the functional analysis of the variational calculus are investigated. Rayleigh-Ritz method provides a way to get the explicit form of the continuous deflection function in which the total potential energy is minimized with respect to the unknown coefficients consisting of the trial functions. As a closure, the analytically calculated results are compared with the experiments and show good agreements.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 1995.10a
• /
• pp.732-735
• /
• 1995

A meshless method which is the new computational method being developed recently, is applied to the simulation of large deformation problems. Among the many types of meshless methods, the Reproducing Kernel particle method (RKPM) is used and the nearly incompressible hyperelastic materials are employed in simulations. The meshless methods can avoid metsh distortions and mesh entanglements that may frequently happen when the mesh-based methods like finite element method are used for the simulations of largely deformed materials. A general features of meshless methods are reviewed and the formulation of RKPM is presented. Next, the performance of explicit RKPM is demonstrated by examples.

• Journal of Institute of Control, Robotics and Systems
• /
• v.14 no.8
• /
• pp.792-799
• /
• 2008

Motion estimation is crucial in a remote control for its convenience or accuracy. Time delays, however, can occur in the problem because data communication is required through a network. In this paper, state estimation problem with uncertain time delayed measurements is addressed. In dynamic system with noise, after taking measurements, it often requires some time until that is available in the filter algorithm. Standard filters not considering this time delays cannot be used since the current measurement is related with a past state. These delayed measurements are solved with augmented extended Kalman filter, and the uncertainty of delayed time is also resolved based on an explicit formulation. The proposed method is analyzed and verified by simulations.

• Journal of Mechanical Science and Technology
• /
• v.17 no.2
• /
• pp.254-268
• /
• 2003

The explicit formulation of the J-integral and the M-integral is constructed in terms of the stress intensity factor and the higher order stress coefficients for Mode II cracks under small or large scale yielding. Furthermore, the stress intensity factor and the higher order stress coefficients as well are computed with the aid of the two-state J- and the M-integral, which is found to be accurate and efficient. It is found that the contribution from the higher order singularities to the J-integral is closely related to the configuration of the plastic zone.

• Structural Engineering and Mechanics
• /
• v.13 no.4
• /
• pp.387-410
• /
• 2002

Formulation of an 8 nodes assumed strain shell element is presented for the analysis of shells. The stiffness matrix based on the Mindlin-Reissner theory is analytically integrated through the thickness. The element is free of membrane and shear locking behavior by using the assumed strain method such that the element performs very well in modeling of thin shell structures. The material is assumed to be isotropic and laminated composite. The element has six degrees of freedom per node and can model the stiffened plates and shells. A great number of numerical testing carried out for the validation of present 8 node shell element are in good agreement with references.

• Structural Engineering and Mechanics
• /
• v.65 no.1
• /
• pp.91-98
• /
• 2018

Three structure-dependent integration methods with no numerical dissipation have been successfully developed for time integration. Although these three integration methods generally have the same numerical properties, such as unconditional stability, second-order accuracy, explicit formulation, no overshoot and no numerical damping, there still exist some different numerical properties. It is found that TLM can only have unconditional stability for linear elastic and stiffness softening systems for zero viscous damping while for nonzero viscous damping it only has unconditional stability for linear elastic systems. Whereas, both CEM and CRM can have unconditional stability for linear elastic and stiffness softening systems for both zero and nonzero viscous damping. However, the most significantly different property among the three integration methods is a weak instability. In fact, both CRM and TLM have a weak instability, which will lead to an adverse overshoot or even a numerical instability in the high frequency responses to nonzero initial conditions. Whereas, CEM possesses no such an adverse weak instability. As a result, the performance of CEM is much better than for CRM and TLM. Notice that a weak instability property of CRM and TLM might severely limit its practical applications.