• Computers and Concrete
• /
• v.8 no.4
• /
• pp.371-388
• /
• 2011

In this paper, time-continuous constitutive equations for strain rate-dependent materials are presented first, among which those for the overstress and the consistency viscoplastic models are considered. By allowing the stress states to be outside the yield surface, the overstress viscoplastic model directly defines the flow rule for viscoplastic strain rate. In comparison, a rate-dependent yield surface is defined in the consistency viscoplastic model, so that the standard Kuhn-Tucker loading/unloading condition still remains true for rate-dependent plasticity. Based on the formulation of the consistency viscoplasticity, a computational elasto-viscoplastic constitutive model is proposed for the short fiber-reinforced fresh cementitious paste for extrusion purpose. The proposed constitutive model adopts the von-Mises yield criterion, the associated flow rule and nonlinear strain rate-hardening law. It is found that the predicted flow stresses of the extrudable fresh cementitious paste agree well with experimental results. The rate-form constitutive equations are then integrated into an incremental formulation, which is implemented into a numerical framework based on ANSYS/LS-DYNA finite element code. Then, a series of upsetting and ram extrusion processes are simulated. It is found that the predicted forming load-time data are in good agreement with experimental results, suggesting that the proposed constitutive model could describe the elasto-viscoplastic behavior of the short fiber-reinforced extrudable fresh cementitious paste.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.53 no.3
• /
• pp.19-26
• /
• 2011

In order to estimate the viscous effects of clay over a wide range of strain levels, we confirmed the performance of a viscoelastic-viscoplastic constitutive model by simulating cyclic undrained triaxial tests to determine the cyclic strength and deformation characteristics of natural marine clay. The viscoelastic-viscoplastic constitutive model was then incorporated into an effective stress-based seismic response analysis to estimate the effects of an intermediate clay layer on the behavior of sand layers. Seismic response was simulated by the cyclic viscoelastic-viscoplastic constitutive model created with data recorded at Rokko Island, Kobe, Japan. The results show that a cyclic viscoelastic-viscoplastic constitutive model can provide a good description of dynamic behavior including viscoelastic effects, within a small strain range.

• Journal of the Korean Geotechnical Society
• /
• v.18 no.4
• /
• pp.55-64
• /
• 2002

In the present study, in order to apply a cyclic viscoelastic-viscoplastic constitutive model to multi-layered ground conditions during large earthquake, the numerical simulations of the 1995 Hyogoken Nanbu Earthquake at Port Island, Kobe, Japan, were performed by the seismic response analysis. From the seismic response analysis, it was found that the acceleration calculated from the cyclic elasto-viscoplastic model and cyclic viscoelastic-viscoplastic models for clay was in close agreement with the recorded accelerations at the Port Island down-hole array, and the cyclic elastic-viscoplastic and viscoelastic-viscoplastic constitutive models showed little different behavior characteristics near clay layer. Thus, the propriety of viscoplastic model for clay was convinced. Therefore, it can be concluded that a cyclic viscoelastic-viscoplastic constitutive model can give a good description of the amplification and also it showed accurate damping characteristics of clay during large event which induces plastic deformation in large strain range.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.48 no.7
• /
• pp.45-52
• /
• 2006

A cyclic viscoelastic-viscoplastic constitutive model for clay is incorporated into an effective stress based seismic response analysis to describe viscous effect of clay layer to sand layer during earthquake. The seismic response against main shock of 1995 Hyogoken Nambu Earthquake is analyzed in the present study. Acceleration responses in both clay layer and just upper liquefiable sand layer are damped due to viscous effect of clay. A cyclic viscoelastic-viscoplastic constitutive model for clay was implemented into a FEM code, and $Newmark{\beta}$ method was employed for the time discretization in the finite element formulation. Seismic responses were simulated by numerical method with recorded data at Port Island, Kobe, Japan. As results of this study, it was found that a cyclic viscoelastic-viscoplastic constitutive model can give good description of dynamic behavior characteristics including viscoelastic effect.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.29 no.2 s.233
• /
• pp.253-261
• /
• 2005

In decades, a substantial body of work on a unified viscoplastic model which considers the mechanism of plastic deformation and creep deformation has developed. The systematic scheme for numerical analysis of unified model is necessary because the dominant failure mechanism is the defect growth and coalescence in materials. In the present study, the unified viscoplastic model for materials with defects suggested by Suquet and Michel was employed for numerical analysis. The constitutive equations are integrated based on the generalized mid-point rule and implemented into a finite element program (ABAQUS) by means of user-defined subroutine (UMAT). To evaluate the validity of the developed UMAT code and the assessment of the adopted viscoplastic model, the results obtained from the UMAT code was compared with the numerical reference solution and experimental data. The unit cell analysis also has been investigated to study the effect of strain rate, temperature, stress triaxiality and initial defect volume fraction on the growth and coalescence of the defect.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.9
• /
• pp.1555-1561
• /
• 2003

Finite element analysis using modern constitutive equation is one of the most general tools to simulate the deformation behavior and to predict the life of the structure. Constitutive equation becomes complicated so as to predict the material behavior more accurately than the classical models. Because of the complexity of constitutive model, numerical treatment becomes so difficult that the calculation should be verified carefully. One-element tests, simple tension or simple shear, are usually used to verify the accuracy of finite element analysis using complicated constitutive model. Since this test is mainly focused on the time integration scheme, it is also necessary to verify the equilibrium iteration using material stiffness matrix and to compare FE results with solution of structures. In this investigation, viscoplastic solution of thick walled cylinder was derived considering axial constraints and was compared with the finite element analysis. All the numerical solutions showed a good coincidence with FE results. This numerical solution can be used as a verification tool for newly developed FE code with complicated constitutive model.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.44 no.4
• /
• pp.114-122
• /
• 2002

Up to now, many constitutive models for clay have been proposed and studied based on the elasto-plastic or elasto-viscoplastic theory and it has been recognized that the effect of time on the loading process is a salient feature. In the present study, cyclic behavior characteristics of clay was studied with a viscoelastic-viscoplastic constitutive model for clay based on the non-linear kinematic hardening rule. In order to examine the behavior of clay several cyclic untrained triaxial tests and also their numerical simulations were performed. As results of that, it was found that the proposed model can well describe cyclic behaviors of clay such as frequency dependent characteristics, and have the high feasibility of numerical simulation for dynamic analysis.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1995.10a
• /
• pp.178-184
• /
• 1995

A new combined isotropic/kinematic and orthotropic hardening viscoplastic model is proposed which can account for not only differential orientations but also preferred orientations of grains in n metal at finite plastic deformations with an introduction of multiple spin (rate of rotation) concept within the general framework of the model, the effects of anisotropy and constitutive spins will be discussed in conjunction with a closed-form solution for simple shear in n rigid-plastic material, which will be used to simulate experimental data of Montheillet, et al. (1984) for fixed-end tortion tests at finite plastic deformations.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.5 s.176
• /
• 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.

• Journal of Theoretical and Applied Mechanics
• /
• v.3 no.1
• /
• pp.45-65
• /
• 2002

A constitutive model on oorthotropic thermo-elasto-viscoplasticity for fiber-reinforced composite materials Is illustrated, and their thermomechanical responses are predicted with the fully-coupled finite element formulation. The unmixing-mixing scheme can be adopted with the multipartite matrix method as the constitutive model. Basic assumptions based upon the composite micromechanics are postulated, and the strain components of thermal expansion due to temperature change are included In the formulation. Also. more than two sets of mechanical variables, which represent the deformation states of multipartite matrix can be introduced arbitrarily. In particular, the unmixing-mixing scheme can be used with any well-known isotropic viscoplastic theory of the matrix material. The scheme unnecessitates the complex processes for developing an orthotropic viscoplastic theory. The governing equations based on fully-coupled thermomechanics are derived with constitutive arrangement by the unmixing-mixing concept. By considering some auxiliary conditions, the Initial-boundary value problem Is completely set up. As a tool of numerical analyses, the finite element method Is used with isoparametric Interpolation fer the displacement and the temperature fields. The equation of mutton and the energy conservation equation are spatially discretized, and then the time marching techniques such as the Newmark method and the Crank-Nicolson technique are applied. To solve the ultimate nonlinear simultaneous equations, a successive iteration algorithm is constructed with subincrementing technique. As a numerical study, a series of analyses are performed with the main focus on the thermomechanical coupling effect in composite materials. The progress of viscoplastic deformation, the stress-strain relation, and the temperature History are careful1y examined when composite laminates are subjected to repeated cyclic loading.