Viscoplastic approach for rate-dependent failure analysis of concrete joints and interfaces

In this work, a new rate-dependent interface model for computational analysis of quasi-brittle materials like concrete is presented. The model is formulated on the basis of the inviscid elastoplastic model by [Carol, I., Prat, P.C., López, C.M., 1997. “A normal/shear cracking model. Interface implementation for discrete analysis”. Journal of Engineering Mechanics, ASCE, 123 (8), pp. 765–773.]. The rate-dependent extension follows the continuous form of the classical viscoplastic theory by [Perzyna, P., 1966. “Fundamental problems in viscoplasticity”. Advances in Applied Mechanics, 9, pp. 244–368.]. According to [Ponthot, J.P., 1995. “Radial return extensions for viscoplasticity and lubricated friction”. In: Proceedings of International Conference on Structural Mechanics and Reactor Technology SMIRT-13, Porto Alegre, Brazil, (2), pp. 711–722.] and [Etse, G., Carosio, A., 2002. “Diffuse and localized failure predictions of Perzyna viscoplastic models for cohesive-frictional materials”. Latin American Applied Research (32), pp. 21–31.] it includes a consistency parameter and a generalized yield condition for the viscoplastic range that allows an straightforward extension of the full backward Euler method for viscoplastic materials. This approach improves the accuracy and stability of the numerical solution. The model predictions are tested against experimental results on mortar and concrete specimens that cover different stress paths at different strain rates. The results in this work demonstrate, on one hand, the capabilities of the proposed elasto–viscoplastic interface constitutive formulation to predict the rate-dependency of mortar and concrete failure behavior, and, on the other hand, the efficiency of the numerical algorithms developed for the computational implementation of the model that include the consistent tangent operator to improve the convergence rate at the finite element level.     

Sensitivity of simulations depending on material parameter variations

In order to predict the inelastic deformations of shock wave loaded plates simulations using viscoplastic constitutive equations combined with a structural theory are carried out. The most accurate results are determined by comparing calculated and measured plate deflections to each other.     

An implicit consistent algorithm for the integration of thermoviscoplastic constitutive equations in adiabatic conditions and finite deformations

The so-called viscoplastic consistency model, proposed by Wang, Sluys and de Borst, is extended here to the integration of a thermoviscoplastic constitutive equation for J₂ plasticity and adiabatic conditions. The consistency condition in this case includes not only strain rate but also the effect of temperature on the yield function. Using the backward Euler integration scheme to integrate the constitutive equations, an implicit algorithm is proposed, leading to generalized expressions of the classical return mapping algorithm for J₂ plasticity, both for the iterative calculation of the plastic multiplier increment and for the consistent tangent operator when strain rate and temperature are considered also as state variables of the hardening equation. The model was implemented in a commercial finite element code and its performance is demonstrated with the numerical simulation of four Taylor impact tests.     

A theory of strain-gradient plasticity for isotropic, plastically irrotational materials. Part II: Finite deformations

This paper generalizes to finite deformations our companion paper [Gurtin, M.E., Anand, L., 2004. A theory of strain-gradient plasticity for isotropic, plastically irrotational materials. Part I: Small deformations. Journal of the Mechanics and Physics of Solids, submitted]. Specifically, we develop a gradient theory for finite-deformation isotropic viscoplasticity in the absence of plastic spin. The theory is based on the Kröner–Lee decomposition F = F^eF^p of the deformation gradient into elastic and plastic parts; a system of microstresses consistent with a microforce balance; a mechanical version of the second law that includes, via microstresses, work performed during viscoplastic flow; a constitutive theory that allows:

• the microstresses to depend on D^p, the gradient of the plastic stretching,

• the free energy ψ to depend on the Burgers tensor G = F^pCurlF^p.

The microforce balance when augmented by constitutive relations for the microstresses results in a nonlocal flow rule in the form of a tensorial second-order partial differential equation for F^p. The microstresses are strictly dissipative when ψ is independent of the Burgers tensor, but when ψ depends on G the microstresses are partially energetic, and this, in turn, leads to backstresses and (hence) Bauschinger-effects in the flow rule. The typical macroscopic boundary conditions are supplemented by nonstandard microscopic boundary conditions associated with viscoplastic flow, and, as an aid to numerical solution, a weak (virtual power) formulation of the nonlocal flow rule is derived. Finally, the dependences of the microstresses on D^p are shown, analytically, to result in strengthening and possibly weakening of the body induced by viscoplastic flow.     

耦合蠕变损伤的粘塑性模型用于疲劳-蠕变交互作用的研究

该文对550℃下30CrMoNiV5-11转子钢进行了轴向应力控制的带保持时间的疲劳实验.按照Le-maitre有效应力的概念,采用Kachanov损伤演化方程,推导了耦合损伤的粘塑性流动方程和随动硬化方程,使用返回映射和向后欧拉方法实现应力更新,并在有限元软件ANSYS的接口程序UMAT中嵌入用户程序.对试样疲劳-蠕变交互行为进行模拟,结果表明,耦合蠕变损伤变量后,模拟结果较不含损伤变量模型的模拟结果精度大为提高.     

Micromechanical analysis of the finite elastic–viscoplastic response of multiphase composites

Nonlinear thermoelastic–viscoplastic constitutive equations for large deformations with isotropic and directional hardening, are incorporated into a micromechanical finite strain analysis. As a result of this analysis, which is based on the homogenization technique for periodic microstructures, a global thermoinelastic constitutive law is established that governs the overall response of multiphase materials under finite deformations. This constitutive law is expressed in terms of the instantaneous effective mechanical and thermal stress tangent tensors together with the instantaneous global inelastic stress tensor that represents the viscoplastic effects. Results for a thermoinelastic matrix reinforced by a hyperelastic compressible material are given that illustrate the response of fibrous and particulate composites to various types of loading.     

Three-dimensional modeling and numerical analysis of rate-dependent irrecoverable deformation in shape memory alloys

Shape memory alloys (SMAs) provide an attractive solid-state actuation alternative to engineers in various fields due to their ability to exhibit recoverable deformations while under substantial loads. Many constitutive models describing this repeatable phenomenon have been proposed, where some models also capture the effects of rate-independent irrecoverable deformations (i.e., plasticity) in SMAs. In this work, we consider a topic not addressed to date: the generation and evolution of irrecoverable viscoplastic strains in an SMA material. Such strains appear in metals subjected to sufficiently high temperatures. The need to account for these effects in SMAs arises when considering one of two situations: the exposure of a conventional SMA material (e.g., NiTi) to high temperatures for a non-negligible amount of time, as occurs during shape-setting, or the utilization of new high temperature shape memory alloys (HTSMAs), where the elevated transformation temperatures induce transformation and viscoplastic behaviors simultaneously. A new three-dimensional constitutive model based on established SMA and viscoplastic modeling techniques is derived that accounts for these behaviors. The numerical implementation of the model is described in detail. Several finite element analysis (FEA) examples are provided, demonstrating the utility of the new model and its implementation in assessing the effects of viscoplastic behaviors in shape memory alloys.     

DISSIPATION-BASED CONSISTENT RATE-DEPENDENT MODEL FOR CONCRETE

An energy-dissipation based viscoplastic consistency model is presented to describe the performance of concrete under dynamic loading. The development of plasticity is started with the thermodynamic hypotheses in order that the model may have a sound theoretical background. Independent hardening and softening and the rate dependence of concrete are described separately for tension and compression. A modified implicit backward Euler integration scheme is adopted for the numerical computation. Static and dynamic behavior of the material is illustrated with certain numerical examples at material point level and structural level, and compared with existing experimental data. Results validate the effectiveness of the model.     

On Deformations with the Same Principal Stretches

In the context of the finite strain theory, plane isochoric homogeneous deformations are considered. Inspired by two examples (plane elliptical and plane hyperbolic deformations), it is seen that for any such isochoric deformation the corresponding principal stretches are equal to those of simple shear provided there is a certain relation between the amount of shear of the simple shear and the parameters of the general plane deformation. Then, the link is established between any two homogeneous deformations which have identical principal stretches. It involves two rotations, one in the undeformed state and the other in the deformed state. These rotations are determined explicitly for an arbitrary isochoric homogeneous deformation and the simple shear with the same principal stretches.     

The effect of annealing on the viscoplastic response of semicrystalline polymers at finite strains

A constitutive model is developed for the viscoplastic behavior of a semicrystalline polymer at finite strains. A solid polymer is treated as an equivalent heterogeneous network of chains bridged by permanent junctions (physical cross-links, entanglements and lamellar blocks). The network is thought of as an ensemble of meso-regions linked with each other. In the sub-yield region of deformations, junctions between chains in meso-domains slide with respect to their reference positions (which reflects sliding of nodes in the amorphous phase and fine slip of lamellar blocks). Above the yield point, this sliding process is accompanied by displacements of meso-domains in the ensemble with respect to each other (which reflects coarse slip and disintegration of lamellar blocks). To account for the orientation of lamellar blocks in the direction of maximal stresses and formation of micro-fibrils in the post-yield region of deformations (which is observed as strain-hardening of specimens) elastic moduli are assumed to depend on the principal invariants of the right Cauchy–Green tensor for the viscoplastic flow. Stress–strain relations for a semicrystalline polymer are derived by using the laws of thermodynamics. The constitutive equations are determined by six adjustable parameters that are found by matching observations in uniaxial tensile tests on injection-molded isotactic polypropylene at elongations up to 80%. Prior to testing, the specimens were annealed at various temperatures ranging from 110 to 163 °C. Fair agreement is demonstrated between the experimental data and the results of numerical simulation. The effect of annealing temperature on the material parameters is studied in detail.     

材料的弹粘塑性本构关系和应力波的演化

试验表明，大多数工程材料在冲击载荷作用之下的变形一般都同时包含有可恢复的瞬态性弹性变形和不可恢复的粘滞性塑性变形，即其本构关系可以用弹粘塑性模型来描述。本文从内变量理论出发，探讨了时率相关材料的弹粘塑性本构关系的一般特性，建立了增量型的弹粘塑性本构关系的一般理论框架和普适的表达式，并且对两种最常用的本构模型——Bodner-Partom模型和Johnson-Cook模型给出了在一维应变条件下的具体形式。通过计算和讨论一维应变粘塑性靶板中冲击波的衰减机制和应力波的演化规律，特别是考察各种粘塑性本构模型中的材料参数对冲击波的衰减和应力波的演化的影响，得出了一些可以直接应用或具有一定借鉴价值的结果，为研究应力波的其他衰减机制以及在人防工程中智能防护层设计时新材料的选取奠定了基础。     

A theory of amorphous solids undergoing large deformations,with application to polymeric glasses

This paper develops a continuum theory for the elastic–viscoplastic deformation of amorphous solids such as polymeric and metallic glasses. Introducing an internal-state variable that represents the local free-volume associated with certain metastable states, we are able to capture the highly non-linear stress–strain behavior that precedes the yield-peak and gives rise to post-yield strain softening. Our theory explicitly accounts for the dependence of the Helmholtz free energy on the plastic deformation in a thermodynamically consistent manner. This dependence leads directly to a backstress in the underlying flow rule, and allows us to model the rapid strain-hardening response after the initial yield-drop in monotonic deformations, as well as the Bauschinger-type reverse-yielding phenomena typically observed in amorphous polymeric solids upon unloading after large plastic deformations. We have implemented a special set of constitutive equations resulting from the general theory in a finite-element computer program. Using this finite-element program, we apply the specialized equations to model the large-deformation response of the amorphous polymeric solid polycarbonate, at ambient temperature and pressure. We show numerical results to some representative problems, and compare them against corresponding results from physical experiments.     

Design sensitivity analysis for parameters affecting geometry,elastic–viscoplastic material constant and boundary condition by consistent tangent operator-based boundary element method

This paper presents a design sensitivity analysis method by the consistent tangent operator concept-based boundary element implicit algorithm. The design variables for sensitivity analysis include geometry parameters, elastic–viscoplastic material parameters and boundary condition parameters. Based on small strain theory, Perzyna’s elastic–viscoplastic material constitutive relation with a mixed hardening model and two flow functions is considered in the sensitivity analysis. The related elastic–viscoplastic radial return algorithm and the formula of elastic–viscoplastic consistent tangent operator are derived and discussed. Based on the direct differentiation approach, the incremental boundary integral equations and related algorithms for both geometric and elastic–viscoplastic sensitivity analysis are developed. A 2D boundary element program for geometry sensitivity, elastic–viscoplastic material constant sensitivity and boundary condition sensitivity has been developed. Comparison and discussion with the results of this paper, analytical solution and finite element code ANSYS for four plane strain numerical examples are presented finally.     

On internal dissipation inequalities and finite strain inelastic constitutive laws: Theoretical and numerical comparisons

Internal dissipation always occurs in irreversible inelastic deformation processes of materials. The internal dissipation inequalities (specific mathematical forms of the second law of thermodynamics) determine the evolution direction of inelastic processes. Based on different internal dissipation inequalities several finite strain inelastic constitutive laws have been formulated for instance by Simo [Simo, J.C., 1992. Algorithms for static and dynamic multiplicative plasticity that preserve the classical return mapping schemes of the infinitesimal theory. Computer Methods in Applied Mechanics and Engineering 99, 61–112]; Simo and Miehe [Simo, J.C., Miehe, C., 1992. Associative coupled thermoplasticity at finite strains: formulation, numerical analysis and implementation. Computer Methods in Applied Mechanics and Engineering 98, 41–104]; Lion [Lion, A., 1997. A physically based method to represent the thermo-mechanical behavior of elastomers. Acta Mechanica 123, 1–25]; Reese and Govindjee [Reese, S., Govindjee, S., 1998. A theory of finite viscoelasticity and numerical aspects. International Journal of Solids and Structures 35, 3455–3482]; Lin and Schomburg [Lin, R.C., Schomburg, U., 2003. A finite elastic–viscoelastic–elastoplastic material law with damage: theoretical and numerical aspects. Computer Methods in Applied Mechanics and Engineering 192, 1591–1627]; Lin and Brocks [Lin, R.C., Brocks, W., 2004. On a finite strain viscoplastic theory based on a new internal dissipation inequality. International Journal of Plasticity 20, 1281–1311]; and Lin and Brocks [Lin, R.C., Brocks, W., 2005. An extended Chaboche’s viscoplastic law at finite strains: theoretical and numerical aspects. Journal of Materials Science and Technology 21, 145–147]. These constitutive laws are consistent with the second law of thermodynamics. As the internal dissipation inequalities are described in different configurations or coordinate systems, the related constitutive laws are also formulated in the corresponding configurations or coordinate systems. Mathematically, these constitutive laws have very different formulations. Now, a question is whether the constitutive laws provide identical constitutive responses for the same inelastic constitutive problems. In the present work, four types of finite strain viscoelastic and endochronically plastic laws as well as three types of J₂-plasticity laws are formulated based on four types of dissipation inequalities. Then, they are numerically compared for several problems of homogeneous or complex finite deformations. It is demonstrated that for the same inelastic constitutive problem the stress responses are identical for deformation processes without rotations. In the simple shear deformation process with large rotation, the presented viscoelastic and endochronically plastic laws also show almost identical stress responses up to a shear strain of about 100%. The three laws of J₂-plasticity also produce the same shear stresses up to a shear strain of 100%, while different normal stresses are generated even at infinitesimal shear strains. The three J₂-plasticity laws are also compared at three complex finite deformation processes: billet upsetting, cylinder necking and channel forming. For the first two deformation processes similar constitutive responses are obtained, whereas for the third deformation process (with large global rotations) significant differences of constitutive responses can be observed.     

应变控制的热机械疲劳行为的数值模拟

根据高温合金材料的力学性能，以弹粘塑性本构模型为基础，用数值模拟方法研究材料的热机械疲劳循环特性，模型将应变分为弹性应变、温度应变和粘塑性应变三部分，认为材料在高温循环载荷下呈现明显的弹粘塑性特征，根据虚位移原理建立轴对称体的弹粘塑性计算有限元格式，对于循环机械载荷和循环温度载荷，程序中采用了增量法迭代求解，在非线性项中不仅考虑了机械载荷增量的影响，同时也考虑了温度增量的影响，根据应变控制热机械疲劳的特点，发展了应变增量法的有限元计算方法、通过数值模拟，得到材料在各种循环载荷下的应力—应变响应，数值模拟较好地反映了粘塑性变形过程以及温度变化的效应，所描述的不可逆系统在某一时刻的状态完全由当时的状态参数、内变量、承载时间及塑性应变累积量决定，对带缺口试件的模拟结果显示了程序对复杂轴对称结构进行热机械疲劳计算的有效性。     

Constitutive equations for thermomechanical deformations of glassy polymers

Poly-Carbonate (PC) and Poly-Methyl-Methacrylate (PMMA) are lightweight and mechanically tough transparent glassy polymers. Their mechanical behavior at low to moderate strain rates has been well characterized; however, that at high strain rates needs additional work. We propose two modifications to existing pressure-dependent viscoplastic constitutive equations that enable one to simulate better mechanical deformations of PC and PMMA at high strain rates. First, the elastic moduli are taken to depend upon the current temperature and the current effective strain rate. Second, two internal variables are introduced to better characterize the strain softening of the material at high strain rates. A technique to find values of newly introduced material parameters is described. We compute the local temperature rise due to energy dissipated during plastic deformations. The true axial stress vs. the true axial strain curves in uniaxial compression from numerical simulations of the test configurations at high strain rates using the proposed constitutive equations are found to agree well with the experimental results available in the literature.     

Thermomechanical material modelling based on a hybrid free energy density depending on pressure,isochoric deformation and temperature

In order to represent temperature-dependent mechanical material properties in a thermomechanical consistent manner it is common practice to start with the definition of a model for the specific Helmholtz free energy. Its canonical independent variables are the Green strain tensor and the temperature. But to represent calorimetric material properties under isobaric conditions, for example the exothermal behaviour of a curing process or the dependence of the specific heat on the temperature history, the temperature and the pressure should be taken as independent variables. Thus, in the field of calorimetry the Gibbs free energy is usually used as thermodynamic potential whereas in continuum mechanics the Helmholtz free energy is normally applied. In order to simplify the representation of calorimetric phenomena in continuum mechanics a hybrid free energy density is introduced. Its canonical independent variables are the isochoric Green strain tensor, the pressure and the temperature. It is related to the Helmholtz free energy density by a Legendre transformation. In combination with the additive split of the stress power into the sum of isochoric and volumetric terms this approach leads to thermomechanical consistent constitutive models for large deformations. The article closes with applications of this approach to finite thermoelasticity, curing adhesives and the glass transition.     

On the structure of the theory of viscoplasticity

In various formulations of plasticity, there is evident a structure embracing several features, including inviscidity, a yield condition, and a constitutive inequality. By means of these features the constitutive equations of plasticity are derived. In the present paper we introduce a viscoplastic counterpart of the constitutive inequality of plasticity, and we consider its physical significance. We also present a theory of viscoplasticity having a structure similar to that of plasticity and its relation with the Hohenemser-Prager prototype of viscoplastic constitutive relations is considered.     

Paradoxical Bending Behavior of Shearable Nonlinearly Elastic Rods

In nonlinear elasticity the exact geometry of deformation is combined with general constitutive relations. This allows a very sophisticated interaction of deformations in different material directions. Based on the Cosserat theory for planar deformations of nonlinearly elastic rods we demonstrate some paradoxical bending effects caused by a nontrivial interaction of extension, flexure, and shear. The analytical results are illustrated by numerical examples.