• Journal of Powder Materials
• /
• v.29 no.4
• /
• pp.303-308
• /
• 2022

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 return-mapping 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.

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

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.6
• /
• pp.1455-1464
• /
• 1994

A new and efficient algorithm for the integration of the thermal-elasto-plastic constitutive equation is proposed. While it falls into the category of the return mapping method, the algorithm adopts the three point approximation of plastic corrector within one time increment step. The results of its application to a von Mises-type thermal-elasto-plastic model with combined hardening and temperature-dependent material properties show that the accurate iso-error maps are obtained for both angular and radial errors. The accuracy achieved is because the predicted stress increment in a single step calculation follows the exact value closely not only at the end of the step but also through the whole path. Also, the comparison of the computational time for the new and other algorithms shows that the new one is very efficient.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.5 no.3
• /
• pp.57-64
• /
• 2001

For the dynamic nonlinear analysis of stiffened plate and shell structures, total Lagrangian formulation is presented based upon the degenerated shell element considering finite rotation effects. Assumed strain concept is adopted in order to overcome shear locking phenomena and to eliminate spurious zero energy mode. In the elasto-plastic analysis, the return mapping algorithm based on the consistent elasto-plastic tangent modulus is applied to collapse analysis of shell structures. Newmark integration method is used for dynamic nonlinear analysis of shell structures under dynamic forces.

• Journal of the Korean Geosynthetics Society
• /
• v.5 no.2
• /
• pp.9-15
• /
• 2006

The numerical analysis program using artificial displacement method is developed to analyze the deformation behavior of excavated behind ground of retention wall. The elasto-plastic model suggested by Drucker-Prager was used to represent soil behavior and the model's solution was obtained from the return mapping method. To validate of the program, the predicted results by the numerical analysis and the measured results by a field test are compared. The results of numerical analysis showed good agreement with the measured results in field and theoretical values.

• Proceedings of the IEEK Conference
• /
• 2002.06a
• /
• pp.307-310
• /
• 2002

This paper is proposed due to the fact that resolving domain name with conventional input method to PDA is not so convenient. The substance of this paper is that a user Pronounces tile character which represents tile domain name and that vocalized character is transferred through the GATEWAY, where DNS service can be received in return. PDA receives, compress and send the voice to the GATEWAY, Then, the GATEWAY uncompress, recognizes the voice, converts to characters, search for the mapping entry After mapping to the mapping entry, the GATEWAY sends the DNS request. Combining two entities makes DNS service based on the human voice possible.

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

• Transactions of Materials Processing
• /
• v.7 no.3
• /
• pp.233-245
• /
• 1998

A new approach has been proposed for the incremental analysis of the nonsteady state large deformation of planar anisotropic elastic-plastic sheet forming. A mathematical brief review of a constitutive law for the incremental deformation theory has been presented from flow theory using the minimum plastic work path for elastic-plastic material. Since the material embedded coordinate system(Lagrangian quantity) is used in the proposed theory the stress integration procedure is completely objective. A new return mapping algorithm has been also developed from the general midpoint rule so as to achieve numerically large strain increment by successive control of yield function residuals. Some numerical tests for the return mapping algorithm were performed using Barlat's six component anisotropic stress potential. Performance of the proposed algorithm was shown to be good and stable for a large strain increment, For planar anisotropic sheet forming updating algorithm of planar anisotropic axes has been newly proposed. In order to show the effectiveness and validity of the present formulation earing simulation for a cylindrical cup drawing and front fender stamping analysis are performed. From the results it has been shown that the present formulation can provide a good basis for analysis for analysis of elastic-plastic sheet metal forming processes.

• Journal of the Korean Geotechnical Society
• /
• v.22 no.9
• /
• pp.45-59
• /
• 2006

The successful performance of any numerical geotechnical simulation depends on the accuracy and efficiency of the numerical implementation of constitutive model used to simulate the stress-strain (constitutive) response of the soil. The corner stone of the numerical implementation of constitutive models is the numerical integration of the incremental form of soil-plasticity constitutive equations over a discrete sequence of time steps. In this paper a well known two-surface soil plasticity model is implemented using a generalized implicit return mapping algorithm to arbitrary convex yield surfaces referred to as the Closest-Point-Projection method (CPPM). The two-surface model describes the nonlinear behavior of coarse-grained materials by incorporating a bounding surface concept together with isotropic and kinematic hardening as well as fabric formulation to account for the effect of fabric formation on the unloading response. In the course of investigating the performance of the CPPM integration method, it is proven that the algorithm is an accurate, robust, and efficient integration technique useful in finite element contexts. It is also shown that the algorithm produces a consistent tangent operator $\frac{d\sigma}{d\varepsilon}$ during the iterative process with quadratic convergence rate of the global iteration process.

• Journal of Powder Materials
• /
• v.30 no.1
• /
• pp.29-34
• /
• 2023

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.