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

• Computers and Concrete
• /
• 제9권1호
• /
• pp.21-33
• /
• 2012

This paper presents the numerical simulation of the rigid 12.6 mm diameter kinetic energy ogive-nosed projectile impact on plain and fiber reinforced concrete (FRC) targets with compressive strengths from 45 to 235 MPa, using a three-dimensional finite element code LS-DYNA. A combined dynamic constitutive model, describing the compressive and tensile damage of concrete, is implemented. A modified Johnson_Holmquist_Cook (MJHC) constitutive relationship and damage model are incorporated to simulate the concrete behavior under compression. A tensile damage model is added to the MJHC model to analyze the dynamic fracture behavior of concrete in tension, due to blast loading. As a consequence, the impact damage in targets made of plain and fiber reinforced concrete with same matrix material under same impact velocities (650 m/s) are obtained. Moreover, the damage distribution of concrete after penetration is procured to compare with the experimental results. Numerical simulations provide a reasonable prediction on concrete damage in both compression and tension.

• 한국소성가공학회:학술대회논문집
• /
• 한국소성가공학회 2003년도 추계학술대회논문집
• /
• pp.191-194
• /
• 2003

The results from a systematic study of the response of a Ti-6Al-4V alloy under quasi-static and dynamic loading at different strain rates and temperatures are presented. It has been shown that the work-hardening rate decreased as the strain rate and the strain increased. The correlations and predictions using modified KHL (Khan-Huang-Liang) viscoplastic constitutive model are compared with those from JC (Johnson-Cook) model and experimental observations. Overall, KHL model correlations and predictions compared much more favorably than the corresponding JC model predictions and correlations.

• 소성∙가공
• /
• 제21권4호
• /
• pp.221-227
• /
• 2012

The automotive and electronic industries have seriously considered the use of magnesium alloys because of their excellent properties such as strength to weight ratio, EMI shielding capability, etc. However, it is difficult to form magnesium alloys at room temperature because of the mechanical deformation related to twinning. Hence, magnesium alloys are normally formed at elevated temperatures. In this study, a temperature dependent constitutive model, the C-H/V model, for the magnesium alloy AZ31B sheet is proposed. A hardening law based on nonlinear kinematic and H/V(Hollomon/Voce) hardening model is used to properly characterize the Bauschinger effect and the stabilization of the flow stress. Material parameters were determined from a series of uni-axial cyclic experiments(C-T-C) with the temperature ranging between 150 and $250^{\circ}C$. The developed models are fit to experimental data and a comparison is made.

• Computers and Concrete
• /
• 제25권2호
• /
• pp.181-192
• /
• 2020

Nowadays, steel fiber reinforced concretes (SFRCs) are widely used in practical applications. Significant experimental research has thus been carried out to determine the constitutive equations that represent the behavior of SFRCs under multiaxial loadings. However, numerical modelling of SFRCs via FEM has been challenging due to the complexities of the implementation of these constitutive equations. In this study, following the literature, a plasticity model is constructed for the behavior of SFRCs that involves the Willam-Warnke failure surface with the relevant evolution laws and a non-associated flow rule for determining the plastic deformations. For the precise (yet rapid) integration of the constitutive equations, an explicit substepping scheme consisting of yield intersection and drift correction algorithms is employed and thus implemented in ABAQUS via UMAT. The FEM model includes various material parameters that are determined from the experimental data. Three sets of parameters are used in the numerical simulations. While the first set is from the experiments that are conducted in this study on SFRC specimens with various contents of steel fibers, the other two sets are from the experiments reported in the literature. The response of SFRCs under multiaxial compression obtained from various numerical simulations are compared with the experimental data. The good agreement between numerical results and the experimental data indicates that not only the adopted plasticity model represents the behavior of SFRCs very well but also the implemented integration scheme can be employed in practical applications of SFRCs.

• Coupled systems mechanics
• /
• 제12권4호
• /
• pp.297-313
• /
• 2023

This paper deals with delamination analysis of non-linear viscoelastic multilayered beam subjected to bending in the plane of the layers. For this purpose, first, a non-linear viscoelastic model is presented. In order to take into account the non-linear viscoelastic behaviour, a non-linear spring and a non-linear dashpot are assembled in series with a linear spring connected in parallel to a linear dashpot. The behaviours of the non-linear spring and dashpot are described by applying non-linear stress-strain and stress-rate of strain relationships, respectively. The constitutive law of the model is derived. Due to the non-linear spring and dashpot, the constitutive law is non-linear. This law is used for describing the time-dependent mechanical behaviour of the beam under consideration. The material properties involved in the constitutive law vary along the beam length due to the continuous material inhomogeneity of the layers. Solution of the strain energy release rate for the delamination is obtained by analyzing the balance of the energy with considering of the non-linear viscoelastic behaviour. The strain energy release rate is found also by using the complementary strain energy for verification. A parametric study is carried-out by using the solution obtained. The solutions derived and the results obtained help to understand the time-dependent delamination of non-linear viscoelastic beams under loading in the plane of layers.

• KCI Concrete Journal
• /
• 제13권1호
• /
• pp.19-25
• /
• 2001

Based on the orthotropic hypoelasticity formulation, a constitutive material model of concrete taking account of triaxial stress state is presented. In this model, the ultimate strength surface of concrete in triaxial stress space is described by the Hsieh's four-parameter surface. On the other hand, the different ultimate strength surface of concrete in strain space is proposed in order to account for increasing ductility in high confinement pressure. Compressive ascending and descending behavior of concrete is considered. Concrete cracking behavior is considered as a smeared crack model, and after cracking, the tensile strain-softening behavior and the shear mechanism of cracked concrete are considered. The proposed constitutive model of concrete is compared with some results obtained from tests under the states of uniaxial, biaxial, and triaxial stresses. In triaxial compressive tests, the peak compressive stress from the predicted results agrees well with the experimental results, and ductility response under high confining pressure matches well the experimental result. The reinforcing bars embedded in concrete are considered as an isoparametric line element which could be easily incorporated into the isoparametric solid element of concrete, and the average stress - average strain relationship of the bar embedded in concrete is considered. From numerical examples for a reinforced concrete simple beam and a structural beam type member, the stress state of concrete in the vicinity of talc critical region is investigated.

• 복합신소재구조학회 논문집
• /
• 제5권4호
• /
• pp.11-17
• /
• 2014

A progressive failure analysis procedure for composite laminates is completed in here. An anisotropic plastic constitutive model for fiber-reinforced composite material is implemented into computer program for a predictive analysis procedure of composite laminates. Also, in order to describe material behavior beyond the initial yield, the anisotropic work-hardening model and subsequent yield surface are 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 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.

• 제어로봇시스템학회:학술대회논문집
• /
• 제어로봇시스템학회 1988년도 한국자동제어학술회의논문집(국제학술편); 한국전력공사연수원, 서울; 21-22 Oct. 1988
• /
• pp.1018-1021
• /
• 1988

This study is concerned with modeling an elastomer constitutive relation by utilizing the truncated Volterra series. Actual experimental data from the Instron Tester are obtained for combined input, i.e. constant strain rate followed by a constant strain input. These data are then estimated for step inputs and utilized for the truncated Volterra series models. One second order and one third order truncated Volterra series models have been employed to estimated the force-displacement relation which is one of the prominent properities to characterize the viscoelastic material. The third order Volterra series model has better results, compared with those of the second order Volterra series model.

• Journal of Mechanical Science and Technology
• /
• 제20권5호
• /
• pp.621-633
• /
• 2006

To reflect the size effect of material $(1\sim15{\mu}m)$ during plastic deformation of polycrystalline copper, a constitutive equation which includes the strain gradient plasticity theory and intrinsic material length model is coupled with the finite element analysis and applied to plane strain deformation problem. The method of least square has been used to calculate the strain gradient at each element during deformation and the effect of distributed force on the strain gradient is investigated as well. It shows when material size is less than the intrinsic material length $(1.54{\mu}m)$, its deformation behavior is quite different compared with that computed from the conventional plasticity. The generation of strain gradient is greatly suppressed, but it appears again as the material size increases. Results also reveal that the strain gradient leads to deformation hardening. The distributed force plays a role to amplify the strain gradient distribution.