• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.22 no.2
• /
• pp.444-450
• /
• 1998

In the present investigation, it is attempted to derive a yield function and associated flow rules of loose metal powders to predict plastic deformation and density change during compaction. The loose metal powders yield by shear stress as well as hydrostatic stress and the yield strength is much smaller in tension than compression. Therefore, a yield function for the powders is expressed as a shifted ellipse toward the negative direction in the hydrostatic stress axis in the space defined by the two stresses. Each of parameters A, B and .delta. used in the yield function is expressed as a function of relative density and it is determined by uniaxial strain and hydrostatic compressions using Cu powder. Flow rules obtained by imposing the normality rule to the yield function are applied to the analyses of unidirectional, bidirectional and hydrostatic compressions, resulting in an excellent agreement with experiments. The yield function is further examined by checking volume changes in plane stain, uniaxial strain and shear deformations.

• Proceedings of the KSME Conference
• /
• 2001.11a
• /
• pp.746-752
• /
• 2001

Design sensitivity analysis for nonlinear structural problems has been emerged in the last decade as a glowing area of engineering research. As a result, theoretical formulations and computational algorithms have already developed for design sensitivity of nonlinear structural problems. There is not enough research for practical nonlinear problems using multi-element, due to difficulties of implementation into FEA. Therefore, nonlinear response analysis for stiffened shell which consists of Mindlin plate and Timoshenko beam, was considered. Specially, it presents the backward-Euler method which is adopted to describe an exact yield state in the stress computation procedure. Then, design sensitivity analysis of nonlinear structures, particularly elasto-perfectly-plastic structure, is developed using direct differentiation method. The accuracy of the developed sensitivity analysis was compared with the central finite difference method. Finally, on the basis of above results, design improvement for stiffened shell is suggested.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2003.04a
• /
• pp.45-52
• /
• 2003

A finite element procedure to estimate ultimate strength of space frames considering spread of plasticity is presented. The improved displacement field is introduced based on inclusion of second order terms of finite rotations. All the nonlinear terms due to bending and torsional moment as well as axial force are precisely considered. The concept of plastic hinge is introduced and the incremental load/displacement method is applied for the elasto-plastic analysis. The initial yield surface is defined based on the residual stress and the full plastification surface is considered under the combined action of axial force, bending and torsional moments. The elasto-plastic stiffness matrices are derived using the flow rule and the normality condition of the limit function. Finite element solutions for ultimate strength of space frames are compared with available solutions and experimental results.

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

• Geomechanics and Engineering
• /
• v.29 no.6
• /
• pp.615-626
• /
• 2022

The soft clays are widely distributed, and one of the prominent engineering problems is the creep behavior. In order to predict the creep deformation of soft clays in an easier and more acceptable way, a simple creep constitutive model has been proposed in this paper. Firstly, the triaxial creep test data indicated that, the strain-time (𝜀-t) curve showing in the 𝜀-lgt space can be divided into two lines with different slopes, and the time referring to the demarcation point is named as t_EOP. Thereafter, the strain increments occurred after the time t_EOP are totally assumed to be the creep components, and the elastic and plastic strains had occurred before t_EOP. A hyperbolic equation expressing the relationship between creep volumetric strain, stress and time is proposed, with several triaxial creep test data of soft clays verifying the applicability. Additionally, the creep flow law is suggested to be similar with the plastic flow law of the modified Cam-Clay model, and the proposed volumetric strain equation is used to deduced the scaling factor for creep strains. Therefore, a creep constitutive model is thereby established, and verified by successfully predicting the creep principal strains of triaxial specimens.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.46 no.6
• /
• pp.21-28
• /
• 2004

This study aims to investigate the effects of partially distributed loads on the dynamic behaviour of steel parabolic arches by using the elasto-plastic finite element model based on the Von Mises yield criteria and the Prandtl-Reuss How rule. For this purpose, the vertical and the radial load conditions were considered as a distributed loading and the loading range is varied from 40% to 100% of arch span. Normal arch and arch with initial deflection were studied. The initial deflection of arch was assumed by the sinusoidal motile of ${\omega}_i\;=\;{\\omega}_O$ sin ($n{\pi}x/L$). Several numerical examples were tested considering symmetric initial deflection when the maximum initial deflection at the apex is fixed as L/1000. The analysis resluts showed that the maximum deflection at the apex of arch was occurred when 70% of arch span was loaded. The maximum deflection at the quarter point of arch span was occurred when 50% of arch span was loaded. It is known that the optimal rise to span ratio between 0.2 and 0.3 when the vertical or radial distributed load is applied. It is verified that the influence of initial deflection of radial load case is more serious than that of vertical load case.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1997.10a
• /
• pp.76-79
• /
• 1997

In this study, an elasto-plastic finite element code, ESFORM, was developed to analyze sheet stamping processes. A formulation of 4-node degenerated shell element was implemented in the code. Workpiece materials were assumed to have planar anisotropy, and governed by associated flow rule. Explicit time integration method was employed to save computation time and reduce the required computer memory. Penalty method was used to describe interface behavior between workpiece and rigid die. Deep drawing of square cup and front finder stamping processes were simulated by ESFORM>

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

• Computers and Concrete
• /
• v.4 no.5
• /
• pp.377-402
• /
• 2007

The paper presents quasi-static plane strain FE-simulations of strain localization in reinforced concrete beams without stirrups. The material was modeled with two different isotropic continuum crack models: an elasto-plastic and a damage one. In case of elasto-plasticity, linear Drucker-Prager criterion with a non-associated flow rule was defined in the compressive regime and a Rankine criterion with an associated flow rule was adopted in the tensile regime. In the case of a damage model, the degradation of the material due to micro-cracking was described with a single scalar damage parameter. To ensure the mesh-independence and to capture size effects, both criteria were enhanced in a softening regime by nonlocal terms. Thus, a characteristic length of micro-structure was included. The effect of a characteristic length, reinforcement ratio, bond-slip stiffness, fracture energy and beam size on strain localization was investigated. The numerical results with reinforced concrete beams were quantitatively compared with corresponding laboratory tests by Walraven (1978).

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.2
• /
• pp.248-258
• /
• 1992

Three-dimensional rigid-plastic finite element formulation based on the membrane theory was described and a computer program for large deformation analysis was developed. In the formulation, normal and planar anisotropy of sheet material and rotation of the principal axes of anisotropy was taken into consideration. Sheet metal was assumed to be rigid-plastic material obeying Hill's quadratic yield criterion and its associated flow rule. Deep drawing process, as a preliminary test, for normal anisotropic material was analyzed in order to examine the validity of developed finite element program. The results were consistent with the existing finite element solutions or experimental data. The present study was mainly concerned with the influence of planar anisotropy on deformation behaviour. Finite element analysis and experiment were carried out for the whole process of deep drawing of planar anisotropic material. The computational and experimental results on the shape of ear, strain distribution and punch load were in good agreement.