• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2006.05a
• /
• pp.349-352
• /
• 2006

The plastic deformation of polycrystalline materials is induced by changes of the microstructure when the loading is beyond the critical state of stress. Constitutive models for the crystal plasticity have the common objective which relates microscopic single crystals in the crystallographic texture to the macroscopic continuum point. In this paper, a new consistent tangent stiffness for the anisotropic elasto-viscoplastic analysis of polycrystalline deformation is developed, which can be used in the finite element analysis for the slip-dominated large deformation of polycrystalline materials. In order to calculate the consistent tangent stiffness, the state function is defined based on the consistency condition between the elastic and plastic stress. The rate of shearing increment($\Delta{\gamma}^{\alpha}$) is calculated with satisfying the consistency condition. The consistency condition becomes zero when the trial resolved shear stress($\tau^{{\alpha}^*}$) becomes resolved shear stress($\tau^{\alpha}$) at every step. Iterative method is utilized to calculate the rate of shearing increment based on the implicit backward Euler method. The consistent tangent stiffness can be formulated by differentiating the rate of shearing increment with total strain increment after the instant rate of shearing increment converges. The proposed tangent stiffness is applied to the ABAQUS/Standard by implementing in the ABAQUS/UMAT.

• Journal of Welding and Joining
• /
• v.6 no.4
• /
• pp.35-43
• /
• 1988

Welding residual stresses were calculated by two dimensional thermal elasto-plastic analysis using element method. Complicated plastic behavior during heat transfer was simulated with time. Fist, temperature distributions. To consider time varying behavior of material properties and loading and unloading processes, iterative calculation based on initial stiffness method was carried out. The method proposed by Yamata was used in time increment control which determined the accuracy of claculation. comparison with other caculated and experimental results shows fairly good agreement.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.37 no.2
• /
• pp.64-72
• /
• 1995

A series of cubical triaxial tests with three independent principal stresses was per- formed on Baekma river sand( # 40~100). It was found that the major principal strain at failure remained approximately constant for b values larger than about 0.3 for both the drained and undrained condition, and thereafter increased as b value decreased. The test results showed that the direction of the strain increment at failure form acute angles with the failure surfaces for both the drained and undrained condition. The results were thus not in agreement with the normality condition from classic plasticity theory. Howev- er, it was found that the projections of the plastic strain increment vectors on the octahe- dral plane were perpendicular to the failure surface in that plane. Failure strength in terms of effective stress anlaysis was greatly influenced by the variation of intermediate principal stress and so was failure criterion. The effective stress failure surfaces for both the drained and undrained condition were estimated quite well by use of Lade's failure criterion.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.4
• /
• pp.621-629
• /
• 1992

This study is concerned with an application of the Boundary Element Method to 2-dimensional elastoplastic stress analysis on the material nonlinearities. The boundary integral formulation adopted an initial stress equation in the inelastic term. In order to determine the initial stress increment, the increment of initial elastic strain energy due to elastic increment in stressstrain curve was used as the convergence criterion during iterative process. For the validity of this procedure, the results of B.E.M. with constant elements and NISA with linear elements where compared on the thin plate with 2 edge v-notches under static tension and the thick cylinder under internal pressure. And this paper compared the results of using unmedical integral with the results of using semi-analytical integral on the plastic domain integral.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.8 no.1
• /
• pp.159-164
• /
• 1988

The directions of the principal strain increment, stress, and stress increment during rotation of the principal stress axes at any stress level was studied for $K_0$-consolidated clay using torsion shear apparatus with individual control of the vertical stress, the confining pressure, and the shear stress on hollow cylinder specimens under undrained and drained condition. The torsion shear tests were performed according to predetermined stress-paths, which were chosen to cover over the full range of rotation of principal stress axes. The test results indicated that the strain increment vectors at failure coincided with the stress vectors. That is, the direction of strain increment coincided with the direction of stress increment at small stress levels and with the direction of stress at higher stress levels, which indicated that the behavior of clay was transfered from elastic to plastic as the stress level was increased. The applicability of the elastoplastic theory for modeling of the behavior of clay during rotation of the principal stress axes was given.

• Journal of the Korean Society of Safety
• /
• v.7 no.1
• /
• pp.20-29
• /
• 1992

Effects of strain hardening exponents on the retardation behavior of fatigue crack propagation are experimentally investigated. The retardation of fatigue crack propagation seems to be induced by the crack closure at crack tip. The phenomenon of crack closure becomes remarkable with the increment of strain hardening exponent and magnitude of percent peak load. The ratio of crack growth increment(a$\_$d//w$\_$d/) is influenced by a single overloading (a$\_$d/) and estimated plastic zone size (W$\_$d/=2r$\_$y/) is increased according with the increasing of strain ha.dening exponents. The number of retarded crack growth cycles were (N$\_$d/) decreased as the baseline stress intensity factor .ange( K$\_$b/) was increased. Within the limitation of these experimental results obtained under the single overload, an empirical relation between crack retardation ratio (Nd/N＊), strain hardening exponent (n) and percent peak load (％PL) has been proposed as; Nd/N＊＝ exp [PL $.$ PL$.$A(n)＋B(n) ] where, A(n)＝${\alpha}$n＋${\beta}$, B(n)＝${\gamma}$n＋$\delta$, PL＝％PL/100 and ${\alpha}$＝0.78, ${\beta}$＝0.54, ${\gamma}$＝0.58 and $\delta$＝-0.01, It is interesting to note that all these constants are identical for materials such as aluminum(A3203), steel(S4SC), steel(SS41) and stainless steel(SUS316) used in this experimental study.

• Geomechanics and Engineering
• /
• v.16 no.6
• /
• pp.619-626
• /
• 2018

A simple prediction procedure was investigated for calculating the stresses and displacements of a circular opening. Unlike existed approaches, the proposed approach starts each step with a radius increment. The stress for each annulus could be obtained analytically, while strain increments for each step can be determinate numerically from the compatility equation by finite difference approximation, flow rule and Hooke's law. In the successive manner, the distributions of stresses and displacements could be found. It should be noted that the finial radial stress and displacement were equal to the internal supporting pressure and deformation at the tunnel wall, respectively. By assuming different plastic radii, GRC and the evolution curve of plastic radii and internal supporting pressures could be obtained conveniently. Then the real plastic radius can be calculated by using linear interpolation in the evolution curve. Some numerical and engineering examples were performed to demonstrate the accuracy and validity for the proposed procedure. The comparisons results show that the proposed procedure was faster than that in Lee and Pietrucszczak (2008). The influence of annulus number and dilation on the accuracy of solutions was also investigated. Results show that the larger the annulus number was, the more accurate the solutions were. Solutions in Park et al. (2008) were significantly influenced by dilation.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.5
• /
• pp.1562-1572
• /
• 1996

The rigid-plastic finite element formulation including the inertia force is derived and then the rigid-plastic finite elemnt program considering the inertia effect is developed. In order to consider the strain hardening, strain rate hardening and thermal softening effects which are frequentrly observed in high-velocity deformation phenomena, the Johnson-Cook constitutive odel is applied. The developed program is used to simulate two high-velocity deformation problemss ; rod impact test and hdigh-velocity compression precess. As a result of rod impact test simulation, it is found that the siulated result has a good agreement with the experimental observation. Through the high-velocity compression process simulation. it is also found that the accuracy of the simulated results is dependent upon the time increment size and mesh size.

• Journal of the Korean Society for Precision Engineering
• /
• v.21 no.2
• /
• pp.116-121
• /
• 2004

In rigid-plastic finite element method, there is a heavy computation time and convergence problem. In this study, static-explicit rigid-plastic finite element method will be introduced. This method is the way that restrict the convergence interval. It is expected that various results from the numerical analysis will give very useful information for the design of tools in sheet metal forming 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.