• 대한기계학회논문집A
• /
• 제24권5호
• /
• 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.

• 한국전산구조공학회논문집
• /
• 제23권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.

• 대한기계학회논문집A
• /
• 제21권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.

• 대한기계학회논문집A
• /
• 제20권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.

• 대한기계학회논문집A
• /
• 제27권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.

• 한국지능시스템학회논문지
• /
• 제19권1호
• /
• pp.96-101
• /
• 2009

본 논문에서는 제안된 진화적 알고리즘을 바탕으로 한 비탄성 구성방정식의 파라미터를 결정하기 위한 방법을 제시한다. 이 방법의 장점은 오차를 갖고 있는 측정된 데이터들이나 모델 방정식들이 부정확하더라도 적절한 파라미터들이 결정되어진다는 것이다. 실험설계는 단축하중과 일정 온도조건하의 샤보쉬 재료모델의 파라미터 결정에 적합하였다. 동시에 모델의 파라미터들은 실험데이터들과 제안한 방법에 의한 값들과 일치하였다. 다른 방법들에 의한 값들과 비교해 본 결과, 제안한 방법에 의한 응력-변형률 선도는 실제적인 재료거동에 비해 좋게 나타났다.

• 대한기계학회논문집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.

• 한국전산구조공학회논문집
• /
• 제29권3호
• /
• pp.237-244
• /
• 2016

본 연구는 재료비선형 문제를 다루기 위한 비선형 MLS 차분법의 정식화 과정을 제시한다. MLS 차분법은 절점모델을 기반으로 고속 미분근사식을 활용하여 지배 미분방정식을 직접 이산화 하는데, 변수를 변위로 일원화한 Navier 방정식을 사용하여 탄성재료 문제를 다룬 기존의 MLS 차분법은 재료의 구성방정식을 별도로 고려할 수 없다. 본 연구에서는 비선형 재료의 구성방정식을 반영할 수 있는 강정식화를 위해 1차 미분근사를 반복 사용하는 겹미분근사를 고안했다. 응력의 발산으로 표현되는 평형방정식을 그대로 이산화하고 Newton 방법을 적용하여 반복계산을 통해 수렴해를 찾는 비선형 알고리즘을 제시했다. 응력 계산과 내부변수의 갱신은 return mapping 알고리즘을 활용하였고, 알고리즘 접선계수(algorithmic tangent modulus)의 적용을 통해 빠르고 안정적인 반복계산이 가능하도록 하였다. 재생성 시험을 통해 겹미분근사의 정당성을 검증했고, 비선형재료에 대한 인장문제의 해석을 통해 개발된 비선형 MLS 차분 알고리즘의 정확성과 안정성을 확인하였다.

• Structural Engineering and Mechanics
• /
• 제13권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.

• 한국콘크리트학회:학술대회논문집
• /
• 한국콘크리트학회 2000년도 가을 학술발표회 논문집(II)
• /
• 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