• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.5 s.176
• /
• pp.1075-1083
• /
• 2000

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. These equations consist of the definition of inelastic strain and the evolution of the state variables introduced to quantify the irreversible processes occurred in the material. With respect to the definition of the inelastic strain, the inelastic constitutive models can be categorized into elastoplastic model, unified viscoplastic model and separated viscoplastic model and the different integration methods have been applied to each category. In the present investigation, the generalized integration method applicable for various types of constitutive equations is developed and implemented 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 line-search technique in combination with Newton method. The strategy to control the time increment for the improvement of the accuracy of the numerical integration is proposed. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method. 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. These equations consist of the definition of inelastic strain and the evolution of the state variables introduced to quantify the irreversible processes occurred in the material. With respect to the definition of the inelastic strain, the inelastic constitutive models can be categorized into elastoplastic model, unified viscoplastic model and separated viscoplastic model and the different integration methods have been applied to each category. In the present investigation, the generalized integration method applicable for various types of constitutive equations is developed and implemented 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 line-search technique in combination with Newton method. The strategy to control the time increment for the improvement of the accuracy of the numerical integration is proposed. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.6
• /
• pp.627-633
• /
• 2010

This paper presents a method for identifying the parameter set of inelastic constitutive equations, which is based on an Evolutionary Algorithm. The advantage of the method is that appropriate parameters can be identified even when the measured data are subject to considerable errors and the model equations are inaccurate. The design of experiments suited for the parameter identification of a material model by Chaboche under the uniaxial loading and stationary temperature conditions was first considered. Then the parameter set of the model was identified by the proposed method from a set of experimental data. In comparison to those by other methods, the resultant stress-strain curves by the proposed method correlated better to the actual material behaviors.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.21 no.7
• /
• pp.1042-1049
• /
• 1997

LMFBR(Liquid Metal Fast Breeder Reactor) vessel is operated under the high temperatures of 500-550.deg. C. Thus, transient thermal loads were severe enough to cause inelastic deformation due to creep-fatigue and plasticity. For reduction of such inelastic deformations, Y-piece structure in the form of a thermal sleeve is used in LMFBR vessel under repeated start-up, service and shut-down conditions. Therefore, a systematic method for inelastic analysis is needed for design of the Y-piece structure subjected to such loading conditions. In the present investigation, finite element analysis of heat transfer and inelastic thermal stress were carried out for the Y-piece structure in LMFBR vessel under service conditions. For such analysis, ABAQUS program was employed based on the elasto-plastic and Chaboche viscoplastic constitutive equations. Based on numerical data obtained from the analysis, creep-fatigue damage estimation according to ASME Code Case N-47 was made and compared to each other. Finally, it was found out that the numerical predictio of damage level due to creep based on Chaboche unified viscoplastic constitutive equation was relatively better compared to elasto-plastic constitutive formulation.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.11
• /
• pp.3516-3524
• /
• 1996

Service conditions for structures at elevated temperatures in nuclear power plant involve transient thermal and mechanical load levels that are severe enough to caeuse inelastic deformations due to creep and plasticity. Therefore, a systematic mehtod of inelastic analysis is needed for the design of structural components in nuclear poser plants subjected to such loading conditions. In the present investigation, the Chabodhe model, one of the unified viscoplastic constitutive equations, was selected for systematic inelastic analysis. The material response was integrated based on GMR ( generallized mid-point rule) time integral scheme and provided to ABAQUS as a material subroutine, UMAT program. By comparing results obtaned from uniaxial analysis using the developed UMAT program with those from Runge-Kutta solutions and experimentaiton, the validity of the adopted Chaboche model and the numerical stability and accuracy of the developed UMAT program were verified. In addition, the developed material subroutine was applied for uniaxial creep and tension analyses for the plate with a hole in the center. The application further demonstrates usefulness of the developed program.

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

• Journal of the Korean Institute of Intelligent Systems
• /
• v.19 no.1
• /
• pp.96-101
• /
• 2009

This paper presents a method for identifying the parameter set of inelastic constitutive equations, which is based on an Evolutionary Algorithm. The advantage of the method is that appropriate parameters can be identified even when the measured data are subject to considerable errors and the model equations are inaccurate. The design of experiments suited for the parameter identification of a material model by Chaboche under the uniaxial loading and stationary temperature conditions was first considered. Then the parameter set of the model was identified by the proposed method from a set of experimental data. In comparison to those by other methods, the resultant stress-strain curves by the proposed method correlated better to the actual material behaviors.

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

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.29 no.3
• /
• pp.237-244
• /
• 2016

This paper presents a nonlinear Moving Least Squares(MLS) difference method for material nonlinearity problem. The MLS difference method, which employs strong formulation involving the fast derivative approximation, discretizes governing partial differential equation based on a node model. However, the conventional MLS difference method cannot explicitly handle constitutive equation since it solves solid mechanics problems by using the Navier's equation that unifies unknowns into one variable, displacement. In this study, a double derivative approximation is devised to treat the constitutive equation of inelastic material in the framework of strong formulation; in fact, it manipulates the first order derivative approximation two times. The equilibrium equation described by the divergence of stress tensor is directly discretized and is linearized by the Newton method; as a result, an iterative procedure is developed to find convergent solution. Stresses and internal variables are calculated and updated by the return mapping algorithm. Effectiveness and stability of the iterative procedure is improved by using algorithmic tangent modulus. The consistency of the double derivative approximation was shown by the reproducing property test. Also, accuracy and stability of the procedure were verified by analyzing inelastic beam under incremental tensile loading.

• Structural Engineering and Mechanics
• /
• v.13 no.2
• /
• pp.135-154
• /
• 2002

This paper presents the convected material frame approach to study the nonlinear behavior of inelastic frame structures. The convected material frame approach is a modification of the co-rotational approximation by incorporating an adaptive convected material frame in the basic definition of the displacement vector and strain tensor. In the formulation, each discrete element is associated with a local coordinate system that rotates and translates with the element. For each load increment, the corresponding strain-displacement and nodal force-stress relationships are defined in the updated local coordinates, and based on the updated element geometry. The rigid body motion and deformation displacements are decoupled for each increment. This modified approach incorporates the geometrical nonlinearities through the continuous updating of the material frame geometry. A generalized nonlinear function is used to derive the inelastic constitutive relation and the kinematic hardening is considered. The equation of motion is integrated by an explicit procedure and it involves only vector assemblage and vector storage in the analysis by assuming a lumped mass matrix of diagonal form. Several numerical examples are demonstrated in close agreement with the solutions obtained by the ANSYS code. Numerical studies show that the proposed approach is capable of investigating large deflection of inelastic planar structures and providing an excellent numerical performance.

• Proceedings of the Korea Concrete Institute Conference
• /
• 2000.10b
• /
• pp.1145-1150
• /
• 2000

Concrete behaves as ductile material at high temperature. The existing stress-strain relationship is not valid at high temperature condition. Thus, stress-strain curve of concrete at high temperature is re-established by modifying Saenz's suggestion in this study. A constitutive model of concrete subjected to elevated temperature is also suggested. The model consists of three components; free thermal stain, mechanical strain and thermal creep strain. As the temperature increase, the thermal creep becomes more critical to the failure of concrete. The thermal creep strain of concrete is derived from the modified power-law relation for the steady state creep. The proposed equation for thermal creep employs a Dorn's temperature compensated time theorem