• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.11
• /
• pp.1907-1916
• /
• 2003

During decades, there has been much progress in understanding of the inelastic behavior of the materials and numerous inelastic constitutive equations have been developed. The complexity of these constitutive equations generally requires a stable and accurate numerical method. To obtain the increment of state variable, its evolution laws are linearized by several approximation methods, such as general midpoint rule(GMR) or general trapezoidal rule(GTR). In this investigation, semi-implicit integration schemes using GTR and GMR were developed and implemented into ABAQUS by means of UMAT subroutine. The comparison of integration schemes was conducted on the simple tension case, and simple shear case and nonproportional loading case. The fully implicit integration(FI) was the most stable but amplified the truncation error when the nonlinearity of state variable is strong. The semi-implicit integration using GTR gave the most accurate results at tension and shear problem. The numerical solutions with refined time increment were always placed between results of GTR and those of FI. GTR integration with adjusting midpoint parameter can be recommended as the best integration method for viscoplastic equation considering nonlinear kinematic hardening.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.36 no.4
• /
• pp.437-442
• /
• 2012

To analyze a multibody system, this paper proposes an implicit numerical integration method for joint coordinates subsystem synthesis method. To verify the proposed method, a multibody model for an unmanned robot vehicle, which consists of six identical independent suspension systems, is developed. The symbolic method is applied to compute the system Jacobian matrix for the implicit integration method. The proposed method is also verified by performing rough terrain run-over simulation in comparison with the conventional implicit integration method. In addition, to evaluate the efficiency of the proposed method, the CPU time obtained by using this method is compared with that obtained by using the conventional implicit method.

• Structural Engineering and Mechanics
• /
• v.75 no.5
• /
• pp.541-557
• /
• 2020

A novel family of controllable, dissipative structure-dependent integration methods is derived from an eigen-based theory, where the concept of the eigenmode can give a solid theoretical basis for the feasibility of this type of integration methods. In fact, the concepts of eigen-decomposition and modal superposition are involved in solving a multiple degree of freedom system. The total solution of a coupled equation of motion consists of each modal solution of the uncoupled equation of motion. Hence, an eigen-dependent integration method is proposed to solve each modal equation of motion and an approximate solution can be yielded via modal superposition with only the first few modes of interest for inertial problems. All the eigen-dependent integration methods combine to form a structure-dependent integration method. Some key assumptions and new techniques are combined to successfully develop this family of integration methods. In addition, this family of integration methods can be either explicitly or implicitly implemented. Except for stability property, both explicit and implicit implementations have almost the same numerical properties. An explicit implementation is more computationally efficient than for an implicit implementation since it can combine unconditional stability and explicit formulation simultaneously. As a result, an explicit implementation is preferred over an implicit implementation. This family of integration methods can have the same numerical properties as those of the WBZ-α method for linear elastic systems. Besides, its stability and accuracy performance for solving nonlinear systems is also almost the same as those of the WBZ-α method. It is evident from numerical experiments that an explicit implementation of this family of integration methods can save many computational efforts when compared to conventional implicit methods, such as the WBZ-α method.

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

In this paper, a program for real time simulation of a vehicle is developed. The program uses relative coordinates and BEF(Backward Difference Formula) numerical integration method. Numerical tests showed that the proposed implicit method is more stable in carring out the numerical integration for vehicl dynamics than the explicit method. Hardware requirements for real time simulation are suggested. Algorithms of parallel processing is developed with DSP (digital signal processor).

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2000.04a
• /
• pp.151-158
• /
• 2000

The use of unconditionally stable implicit time integration techniques for pseudo-dynamic tests has been recently proposed and advanced by several researchers such as Thewalt and Mahin Nakashima and Shing. The developed implicit algorithms are based on a-Method of Hugest et al. In this paper a concise summary and explanation of implicit method for Pseudo dynamic test is presented. Especially The a-C method developed by shing at al. has been in-depth evaluated for this study. Important parameters of the a-C method have been analyzed by the simulation test.

• Journal of computational fluids engineering
• /
• v.16 no.4
• /
• pp.72-83
• /
• 2011

A high-order discontinuous Galerkin method for the two-dimensional compressible Navier-Stokes equations was developed on unstructured triangular meshes. For this purpose, the BR2 methd(the second Bassi and Rebay discretization) was adopted for space discretization and an implicit Euler backward method was used for time integration. Numerical tests were conducted to estimate the convergence order of the numerical solutions of the Poiseuille flow for which analytic solutions are available for comparison. Also, the flows around a flat plate, a 2-D circular cylinder, and an NACA0012 airfoil were numerically simulated. The numerical results showed that the present implicit discontinuous Galerkin method is an efficient method to obtain very accurate numerical solutions of the compressible Navier-Stokes equations on unstructured meshes.

• Structural Engineering and Mechanics
• /
• v.28 no.5
• /
• pp.549-570
• /
• 2008

Numerical integration is an efficient approach for nonlinear dynamic analysis. In this paper, general category of the implicit integration errors will be discussed. In order to decrease the errors, Dynamic Relaxation method with modified time step (MFT) will be used. This procedure leads to an alternative algorithm which is very general and can be utilized with any implicit integration scheme. For numerical verification of the proposed technique, some single and multi degrees of freedom nonlinear dynamic systems will be analyzed. Moreover, results are compared with both exact and other available solutions. Suitable accuracy, high efficiency, simplicity, vector operations and automatic procedures are the main merits of the new algorithm in solving nonlinear dynamic problems.

• Proceedings of the Korea Concrete Institute Conference
• /
• 2000.04a
• /
• pp.617-622
• /
• 2000

The use of unconditionally stable implicit time integration techniques for pseudo-dynamic test has been recently proposed and advanced by several researchers such as Thewalt and Mahin, Nakashima and Shing, etc. The developed implicit algorithms are based on the $\alpha$-Method of Huges et al. In this paper, a concise summary and explanation of implicit method for Pseudo dynamic tese is presented. Especially, The $\alpha$-C method developed by shing et al. has been in-depth evaluated for this study. Important parameters of the $\alpha$-C method have been analyzed by the simulation test.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.7
• /
• pp.1-8
• /
• 2003

Numerical characteristics of air chemistry associated with hypersonic flows are described and are compared with those of hydrogen oxygen combustion, applying the partially implicit time integration method to air chemistry. This paper reveals that the time integration of air chemistry needs a chemical Jacobian for stable calculations. However the positive real eigenvalues in air chemistry are relatively smaller than those of hydrogen combustion, and the numerical integration is less sensitive than that with combustion. lt is also found that the application of the partia1ly irnplicit method reduces the computing time without numerical instabilities.

• Wind and Structures
• /
• v.34 no.1
• /
• pp.59-72
• /
• 2022

In this work the numerical results of the flow around a 5:1 rectangular cylinder at Reynolds numbers 3 000 and 40 000, zero angle of attack and smooth incoming flow condition are presented. Implicit Large Eddy Simulations (ILES) have been performed with a high-order accurate spatial scheme and an implicit high-order accurate time integration method. The spatial approximation is based on a discontinuous Galerkin (dG) method, while the time integration exploits a linearly-implicit Rosenbrock-type Runge-Kutta scheme. The aim of this work is to show the feasibility of high-fidelity flow simulations with a moderate number of DOFs and large time step sizes. Moreover, the effect of different parameters, i.e., dimension of the computational domain, mesh type, grid resolution, boundary conditions, time step size and polynomial approximation, on the results accuracy is investigated. Our best dG result at Re=3 000 perfectly agrees with a reference DNS obtained using Nek5000 and about 40 times more degrees of freedom. The Re=40 000 computations, which are strongly under-resolved, show a reasonable correspondence with the experimental data of Mannini et al. (2017) and the LES of Zhang and Xu (2020).