A fixed-point iteration is proposed to integrate the stress and state variables in the incremental analysis of plastic deformation. The Conventional Newton–Raphson method requires a second-order derivative of the yield function to generate a complicated code, and the convergence cannot be guaranteed beforehand. The proposed fixed-point iteration does not require a second-order derivative of the yield function, and convergence is ensured for a given strain increment. The fixed-point iteration is easier to implement, and the computational time is shortened compared with the Newton–Raphson method. The plane-stress condition is considered for the biaxial loading conditions to confirm the convergence of the fixed-point iteration. 3-dimensional tensile specimen is considered to compare the computational times in the ABAQUS/explicit finite element analysis.

The plastic deformation behavior of additively manufactured anisotropic structures are analyzed using the finite element method (FEM). Hill’s quadratic anisotropic yield function is used, and a modified return-mapping method based on dual potential is presented. The plane stress biaxial loading condition is considered to investigate the number of iterations required for the convergence of the Newton-Raphson method during plastic deformation analysis. In this study, incompressible plastic deformation is considered, and the associated flow rule is assumed. The modified returnmapping method is implemented using the ABAQUS UMAT subroutine and effective in reducing the number of iterations in the Newton-Raphson method. The anisotropic tensile behavior is computed using the 3-dimensional FEM for two tensile specimens manufactured along orthogonal additive directions.

A progressive failure analysis procedure for composite laminates is developed in here and in the companion paper. An anisotropic plastic constitutive model for fiber-reinforced composite material, is developed, which is simple and efficient to be implemented into computer program for a predictive analysis procedure of composites. In current development of the constitutive model, an incremental elastic-plastic constitutive model is adopted to represent progressively the nonlinear material behavior of composite materials until a material failure is predicted. An anisotropic initial yield criterion is established that includes the effects of different yield strengths in each material direction, and between tension and compression. Anisotropic work-hardening model and subsequent yield surface are developed to describe material behavior beyond the initial yield under the general loading condition. The current model is implemented into a computer code, which is Predictive Analysis for Composite Structures (PACS), and is presented in the companion paper. The accuracy and efficiency of the anisotropic plastic constitutive model are verified by solving a number of various fiber-reinforced composite laminates with and without geometric discontinuity. The comparisons of the numerical results to the experimental and other numerical results available in the literature indicate the validity and efficiency of the developed model.

An orthotropic plastic constitutive model for fiber-reinforced composite material, is developed, which is simple and efficient to be implemented into computer program for a predictive analysis procedure of composite laminates. An orthotropic initial yield criterion, as well as work-hardening model and subsequent yield surface are established that includes the effects of different yield strengths in each material direction, and between tension and compression. The current model is implemented into a computer code, which is Predictive Analysis for Composite Structures (PACS). The accuracy and efficiency of the anisotropic plastic constitutive model and the computer program PACS are verified by solving a number of various fiber-reinforced composite laminates. The comparisons of the numerical results to the experimental and other numerical results available in the literature indicate the validity and efficiency of the developed model.

콘크리트의 미세공극 혹은 미세균열의 발생과 성장은 콘크리트의 점차적인 물성 저하를 야기한다. 이와같은 손상은 이방성을 가지며 소성과 함께 콘크리트의 비선형거동을 일으키는 주요원인이 된다. 본 논문은 콘크리트의 탄소성 변형 및 손상을 고려하여 콘크리트의 이방성 손상거동을 해석할 수 있는 콘크리트 연속체 손상모델의 개발에 관한 연구이다. 등가 탄성 에너지원리를 이용하여 이방 손상텐서로 표현된 유효탄성텐서를 구하고, 이를 포함하고 있는 열역학 법칙의 자유에너지함수와 소산포텐셜로부터 손상의 전개법칙을 유도한 후, 손상에너지해방률의 함수로 표현한 손상면을 적용하므로써 콘크리트의 이상성손상을 효율적으로 해석 할 수 있는 구성방정식을 유도하였다. 또한 이방성 손상모델에 콘크리트의 소성모델을 도입시켜 탄소성 변형 및 손상을 함께 고려할 수 있는 콘크리트의 연속체 손상모델을 개발하였다. 개발된 손상모델을 유한요소해석 프로그램에 적용하여 1축 및 2축의 여러 조합응력을 받는 콘크리트 모형을 유한요소해석하였으며, 실험결과 또는 타 모델과의 비교로부터 손상모델의 타당성을 검증하였다.