• 한국지반공학회지:지반
• /
• 제12권6호
• /
• pp.87-100
• /
• 1996

본 연구에서는 동반논문(오세붕 & 이승래, 1996)에서 정식화한 비등방경화 구성모델에 대한 내재적인 응력적분 알고리즘의 정확성과 효율성을 검증하였다. 비배수 삼축압축시험경로에 대한 오차평가를 통하여 정확도해석을 수행하였고 수치적인 굴착해석예제를 수행함으로써 정확도 및 수렴도를 분석하였다. 그 결과 제안된 알고리즘이 비등방경화 구성관계에 대하여 음력을 정화하게 적분하고 Newton 법을 이용한 비선형 해석시에 점근적인 2차 수렴도를 확보함을 알 수 있었다. 그리고 이러한 검증을 토대로 지반의 초기조건 및 시공단계를 고려한 유한요소법을 이용하여 실제 굴착문제를 해석하였다. 비등방경화 구성관계를 이용함으로써, 추정된 벽체의 변위는 Cainflay모델에 비하여 실제와 더 유사한 해석결과를 얻을 수 있었다.

• Structural Engineering and Mechanics
• /
• 제38권3호
• /
• pp.283-300
• /
• 2011

A new conception of fundamental tasks in dynamics of the bridge-track-train systems (BTT), with the aim to evaluate moving load's models adequacy, has been developed. The 2D physical models of BTT systems, corresponding to the fundamental tasks, have been worked out taking into account one-way constraints between the moving unsprung masses and the track. A method for deriving the implicit equations of motion, governing vibrations of BTT systems' models, as well as algorithms for numerical integration of these equations, leading to the solutions of high accuracy and relatively short times of simulations, have been also developed. The derived equations and formulated algorithms constitute the basis for numerical simulation of vibrations of the considered systems.

• 한국전산유체공학회지
• /
• 제19권1호
• /
• pp.27-33
• /
• 2014

A new and efficient implicit scheme is proposed to obtain a steady-state solution in time integration and the comparison of characteristics with the approximation ways for the implicit method to solve the incompressible Navier-Stokes equations is provided. The conservative, finite-volume cell-vertex upwind scheme and artificial compressibility method using dual time stepping for time accuracy is applied in this paper. The numerical results obtained indicate that the direct application of Jacobian matrix to the Lower and upper sweeps of implicit LU-SGS leads to better performance as well as convergence regardless of CFL number and true time step than explicit scheme and approximation of Jacobian matrix. The flow simulation around box in uniform flow with unstructured meshes is demonstrated to check the validity of the current formulation.

• 대한기계학회논문집A
• /
• 제27권9호
• /
• pp.1562-1570
• /
• 2003

The prediction of the inelastic behavior of the structure is an essential part of reliability assessment procedure, because most of the failures are induced by the inelastic deformation, such as creep and plastic deformation. During decades, there has been much progress in understanding of the inelastic behavior of the materials and a lot of inelastic constitutive equations have been developed. The complexity of these constitutive equations generally requires a stable and accurate numerical method. The radial return mapping is one of the most robust integration scheme currently used. Nonlinear kinematic hardening model of Armstrong-Fredrick type has recovery term and the direction of kinematic hardening increment is not parallel to that of plastic strain increment. In this case, The conventional radial return mapping method cannot be applied directly. In this investigation, we expanded the radial return mapping method to consider the nonlinear kinematic hardening model and implemented this integration scheme into ABAQUS by means of UMAT subroutine. The solution of the non-linear system of algebraic equations arising from time discretization with the generalized midpoint rule is determined using Newton method and bisection method. Using dynamic yield condition derived from linearization of flow rule, the integration scheme for elastoplastic and viscoplastic constitutive model was unified. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1998년도 추계 학술대회논문집
• /
• pp.48-54
• /
• 1998

Three-dimensional incompressible Navier-Stokes equations have been solved by the node-centered finite volume method with unstructured hybrid grids. The pressure-velocity coupling is handled by the artificial compressibility algorithm and convective fluxes are obtained by Roe's flux difference splitting scheme with linear reconstruction of the solutions. Euler implicit method is used for time-integration. The viscous terms are discretised in a manner to handle any kind of grids such as tetrahedra, prisms, pyramids, hexahedra, or mixed-element grid. The numerical efficiency and accuracy of the present method is critically evaluated for several example problems.

• 천문학회지
• /
• 제40권4호
• /
• pp.91-97
• /
• 2007

The Fokker-Planck (FP) model is one of the commonly used methods for studies of the dynamical evolution of dense spherical stellar systems such as globular clusters and galactic nuclei. The FP model is numerically stable in most cases, but we find that it encounters numerical difficulties rather often when the effects of tidal shocks are included in two-dimensional (energy and angular momentum space) version of the FP model or when the initial condition is extreme (e.g., a very large cluster mass and a small cluster radius). To avoid such a problem, we have developed a new integration scheme for a two-dimensional FP equation by adopting an Alternating Direction Implicit (ADI) method given in the Douglas-Rachford split form. We find that our ADI method reduces the computing time by a factor of ${\sim}2$ compared to the fully implicit method, and resolves problems of numerical instability.

• Journal of applied mathematics & informatics
• /
• 제8권1호
• /
• pp.165-180
• /
• 2001

This paper presents a new difference scheme for numerical solution of stiff system of ODE’s. The present study is mainly motivated to develop an absolutely stable numerical method with a high order of approximation. In this work a double implicit A-stable difference scheme with the sixth order of approximation is suggested. Another purpose of this study is to introduce automatic choice of the integration step size of the difference scheme which is derived from the proposed scheme and the one step scheme of the fourth order of approximation. The algorithm was tested by means of solving the Kreiss problem and a chemical kinetics problem. The behavior of the gas explosive mixture (H₂+ O₂) in a closed space with a mobile piston is considered in test problem 2. It is our conclusion that a hydrogen-operated engine will permit to decrease the emitted levels of hazardous atmospheric pollutants.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제23권2호
• /
• pp.93-114
• /
• 2019

In this article, we introduce a finite difference method for solving the Navier-Stokes equations in rectangular domains. The method is proved to be energy stable and shown to be second-order accurate in several benchmark problems. Due to the guaranteed stability and the second order accuracy, the method can be a reliable tool in real-time simulations and physics-based animations with very dynamic fluid motion. We first discuss a simple convection equation, on which many standard explicit methods fail to be energy stable. Our method is an implicit Runge-Kutta method that preserves the energy for inviscid fluid and does not increase the energy for viscous fluid. Integration-by-parts in space is essential to achieve the energy stability, and we could achieve the integration-by-parts in discrete level by using the Marker-And-Cell configuration and central finite differences. The method, which is implicit and second-order accurate, extends our previous method [1] that was explicit and first-order accurate. It satisfies the energy stability and assumes rectangular domains. We acknowledge that the assumption on domains is restrictive, but the method is one of the few methods that are fully stable and second-order accurate.

• 한국생산제조학회지
• /
• 제23권4호
• /
• pp.337-344
• /
• 2014

While implicit integration methods such as Newton's method have excellent stability for the analysis of stiff and constrained mechanical systems, they have the drawback that the evaluation and LU-factorization of the system Jacobian matrix required at every time step are time-consuming. This paper proposes a Jacobian update-free Newton's method in order to overcome these defects. Because the motions of all bodies in a vehicle model are limited with respect to the chassis body, the equations are formulated with respect to the moving chassis-body reference frame instead of the fixed inertial reference frame. This makes the system Jacobian remain nearly constant, and thus allows the Newton's method to be free from the Jacobian update. Consequently, the proposed method significantly decreases the computational cost of the vehicle dynamic simulation. This paper provides detailed generalized formulation procedures for the equations of motion, constraint equations, and generalized forces of the proposed method.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2006년도 추계 학술대회논문집
• /
• pp.31-35
• /
• 2006

The delta-formulation of the Navier-Stokes equations has been popularly used in the aerodynamics area. Implicit algorithm can be easily implemented in that by using Taylor series expansion. This formulation is extended for an unsteady analysis by using a dual-time integration. In the meanwhile, the incompressible flows with heat transfers which occur in the area of thermo-hydraulics have been solved by a segregated algorithm such as the SIMPLE method, where each equation is discretised by using an under-relaxed deferred correction method and solved sequentially. In this study, the dual-time delta formulation is implemented in the segregated Navier-Stokes solver which is based on the collocated cell-centerd scheme with un unstructured mesh FVM. The pressure correction equation is derived by the SIMPLE method. From this study, it was found that the Euler dual-time method in the delta formulation can be combined with the SIMPLE method.