• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.3
• /
• pp.311-319
• /
• 2003

This paper aims at investigating the exact dynamic response of the system undergoing a exponential impact force from the viewpoints of conservations of momentum and energy. The midpoint method applied in the Newmark's family algorithm is found to be identical to the case of the application of the trapezoidal method which provides conservations of momentum and energy. For the linear impact force the mid point, the trapezoidal and the (n+1) point method exactly meet the conservation characteristics independent of the size of integration interval. On the other hand, constants for the dynamic motion resulting from the nonlinear impact are underestimated or overestimated by these method mentioned above. To overcome this indispensible error, the Simpson 1/3 method as one of multi step methods whose advantages is to use longer time interval with the same number of evaluation functions is adopted for the exact conservations of momentum and energy. Moreover, the suggested method is expected to expand the similar algorithm for the general dynamic motion including finite rotations.

• 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.