• Computational Structural Engineering
• /
• v.10 no.4
• /
• pp.217-228
• /
• 1997

The high precision analysis by the p-version of the finite element method are fairly well established as highly efficient method for linear elastic problems, especially in the presence of stress singularity. It has been noted that the merits of the p-version are accuracy, modeling simplicity, robustness, and savings in user's and CPU time. However, little has been done to exploit their benefits in elasto-plastic analysis. In this paper, the p-version finite element model is proposed for the materially nonlinear analysis that is based on the incremental theory of plasticity using the constitutive equation for work-hardening materials, and the associated flow rule. To obtain the solution of nonlinear equation, the Newton-Raphson method and initial stiffness method, etc are used. Several numerical examples are tested with the help of the square plates with cutout, the thick-walled cylinder under internal pressure, and the circular plate with uniformly distributed load. Those results are compared with the theoretical solutions and the numerical solutions of ADINA

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.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.

• Tunnel and Underground Space
• /
• v.19 no.6
• /
• pp.498-506
• /
• 2009

The stability assessment of a concrete liner of a compressed air storage tunnel should be performed by an approach which is different from that commonly used for the liners of road tunnels, since the liner is exposed to high air pressure. In this study, the stability analysis method for the liner of compressed air storage tunnel is proposed based on the elastic and elasto-plastic solutions of the thick-walled cylinder problem. In case of elastic analysis, the yield initiation condition at the inner boundary is considered as the failure condition of the liner, while the condition which results in the extension of yielding zone to a certain depth is taken as a failure indicator of the liner in the elasto-plastic analysis taking Mohr-Coulomb criterion. The application of the proposed method revealed that the influence of the relative magnitude of boundary loads on the stability of liner is considerable. In particular, noting that the estimation of the outer boundary load may be relatively difficult, it is thought that the precise prediction of outer boundary load is very important in the analysis. Accordingly, the emphasis is put on the selection of the liner installation time, which may govern the magnitude of outer boundary load.