The effect of shear deformation on the post-buckling behavior of imperfect nonlinear viscoelastic columns

Summary The post-buckling behavior of imperfect columns made of nonlinear viscoelastic materials is investigated, taking into account the effect of shear deformation. The material is modeled according to the Leaderman representation of nonlinear viscoelasticity. Solutions are developed, within the elastica and the shear deformation theories, in order to calculate the growth in time of the total deflection. The numerical results establish the importance of the shear and the nonlinear viscoelasticity effects, and of the h/ℓ ratio in the column post-buckling behavior. Accepted for publication 11 November 1996     

Stretch-induced wrinkling of polyethylene thin sheets: Experiments and modeling

This paper presents a study on stretch-induced wrinkling of thin polyethylene sheets when subjected to uniaxial stretch with two clamped ends. Three-dimensional digital image correlation was used to measure the wrinkling deformation. It was observed that the wrinkle amplitude increased as the nominal strain increased up to around 10%, but then decreased at larger strain levels. This behavior is consistent with results of finite element simulations for a hyperelastic thin sheet reported previously (Nayyar et al., 2011). However, wrinkles in the polyethylene sheet were not fully flattened out at large strains (>30%) as predicted for the hyperelastic sheet, but exhibited a residual wrinkle whose amplitude depended on the loading rate. This is attributed to the viscoelastic response of the material. Two different viscoelastic models were adopted in finite element simulations to study the effects of viscoelasticity on wrinkling and to improve the agreement with the experiments, including residual wrinkles and rate dependence. It is found that a parallel network model of nonlinear viscoelasticity is suitable for simulating the constitutive behavior and stretch-induced wrinkling of the polyethylene sheets.     

On uniaxial extension of an elasto-viscous cylinder

Summary This paper deals with slow processes of extension of an elasto-visoous cylinder under conditions where the theory of linear viscoelasticity is applicable. Kinematic dependences are given and four nonstationary problems are solved concerning the extension of a cylinder when one of the following parameters is constant: strain rate, extension rate, stress, and tensile force.Experiments were performed on a low-molecular polyisobutylene at 25° and constant extension rates. The changes in total deformation and stresses, and after unloading, in highly elastic and irreversible deformation, were determined. The effect of the forces of surface tension on the elastic recovery of the samples was taken into account. The formulas of the linear theory of viscoelasticity are shown to apply within the range of extension rates studied. The viscosity and rubbery elasticity modulus were determined and found to agree well with the results obtained during steady shear flows.     

Molecular simulation guided constitutive modeling on finite strain viscoelasticity of elastomers

Viscoelasticity characterizes the most important mechanical behavior of elastomers. Understanding the viscoelasticity, especially finite strain viscoelasticity, of elastomers is the key for continuation of their dedicated use in industrial applications. In this work, we present a mechanistic and physics-based constitutive model to describe and design the finite strain viscoelastic behavior of elastomers. Mathematically, the viscoelasticity of elastomers has been decomposed into hyperelastic and viscous parts, which are attributed to the nonlinear deformation of the cross-linked polymer network and the diffusion of free chains, respectively. The hyperelastic deformation of a cross-linked polymer network is governed by the cross-linking density, the molecular weight of the polymer strands between cross-linkages, and the amount of entanglements between different chains, which we observe through large scale molecular dynamics (MD) simulations. Moreover, a recently developed non-affine network model (Davidson and Goulbourne, 2013) is confirmed in the current work to be able to capture these key physical mechanisms using MD simulation. The energy dissipation during a loading and unloading process of elastomers is governed by the diffusion of free chains, which can be understood through their reptation dynamics. The viscous stress can be formulated using the classical tube model (Doi and Edwards, 1986); however, it cannot be used to capture the energy dissipation during finite deformation. By considering the tube deformation during this process, as observed from the MD simulations, we propose a modified tube model to account for the finite deformation behavior of free chains. Combing the non-affine network model for hyperelasticity and modified tube model for viscosity, both understood by molecular simulations, we develop a mechanism-based constitutive model for finite strain viscoelasticity of elastomers. All the parameters in the proposed constitutive model have physical meanings, which are signatures of polymer chemistry, physics or dynamics. Therefore, parametric materials design concepts can be easily gleaned from the model, which is also demonstrated in this study. The finite strain viscoelasticity obtained from our simulations agrees qualitatively with experimental data on both un-vulcanized and vulcanized rubbers, which captures the effects of cross-linking density, the molecular weight of the polymer chain and the strain rate.     

Approximate description of the inelastic deformation of an isotropic material with allowance for the stress mode

The elastoplastic deformation of an isotropic material is described using constitutive equations and allowing for the stress mode. The equations include two nonlinear functions that relate the first and second invariants of the stress and linear-strain tensors to the stress mode angle. It is proposed to use a linear rather than nonlinear relationship between the first invariants of the tensors. This simplification is validated by comparing calculated and experimental strains under loading with constant and variable stress mode angle     

Calculating the Linear Creep Strains of Viscoelastic Fibers under Tension

The linearity domain for the viscoelastic properties of high-molecular organic fibers is determined. The linearity criteria are coincidence of experimental compliance curves and linearity of isochronic creep curves. Statistical criteria are used to establish linearity. The influence function in the constitutive equation of linear viscoelasticity is an Abel-type power kernel. The calculated and experimental creep strains are in good agreement both at the initial stage of deformation and after long-term loading__________Translated from Prikladnaya Mekhanika, Vol. 41, No. 5, pp. 97–106, May 2005.     

Adequacy of a nonlinear theory of viscoelasticity

Some features of the behavior of viscoelastic materials whose existence leads to the choice of nonlinear constitutive relations are discussed. A classification of such constitutive relations is given and a number of requirements imposed by practice on their adequacy are formulated. A nonlinear theory of viscoelasticity is proposed; this theory offers the advantages over the theory in which stresses are expressed in terms of strains by integral operators of increasing multiplicity. By a one-dimensional example, it is shown that the constitutive operator relations are reciprocal.     

Compressible rubber materials: experiments and simulations

Porous rubber materials are often used in automotive industries. In this paper, a carbon black-filled one is investigated, which is used, for example, as sealing. Such materials are distinguished by viscoelastic behaviour and by a structural compressibility induced by the porous structure. To identify the material behaviour, uniaxial tension tests and hydrostatic compression tests are performed. Therein the main focus of attention lies on the basic elasticity and on the viscoelasticity in the whole loading range. An important observation of these tests is the viscoelastic behaviour under hydrostatic compression, which has to be included in the material model. Because of the two-phase character of cellular rubber, the theory of porous media is taken into account. To model the structural compressibility, a volumetric–isochore split of the deformation gradient is used. Therein the volumetric part includes the aspect of the point of compaction. Finally, the concept of finite viscoelasticity is applied introducing an intermediate configuration. Because of the viscoelastic behaviour under hydrostatic compression, the volumetric–isochore split is taken into account for the nonequilibrium parts, too. Nonlinear relaxation functions are used to model the process-dependent relaxation times and the highly nonlinear behaviour with respect to the deformation and feedrate. The material parameters of the model are estimated using a stochastic identification algorithm.     

Volume changes in glassy polymers

Summary Constitutive equations are derived for compressible glassy polymers at non-isothermal loading with finite strains. The model is based on the theory of temporary networks in its version of adaptive links concept. The stress–strain relations are applied to the analysis of uniaxial extension of a viscoelastic bar. Explicit formulas are developed for time-dependent Young's modulus and Poisson's ratio of the bar at small strains. Results of numerical simulation are compared with experimental data for polycarbonate, polyethylene, and poly(methyl methacrylate). It is demonstrated that (i) longitudinal stresses do not affect the specific free volume in the region of linear viscoelasticity at strains up to about 0.2%, and cause substantial changes in the free volume in the region of nonlinear viscoelasticity at strains about 1.0%; (ii) in the latter case, the increment of the free-volume fraction is proportional to the increase in the specific volume. Received 3 April 1998; accepted for publication 22 June 1998     

Plane nonaxisymmetric deformation of a viscoelastic cylinder with consideration of its physical and geometric nonlinearity

The nonaxisymmetric plane problem of the nonlinear theory of viscoelasticity is solved for a cylinder reinforced by an elastic circular shell. The cylinder has an internal cut resembling a Maltese cross in shape. The identification of the nonlinear endochronous theory of aging viscoelastic materials is conducted by a genetic algorithm method on the basis a nonmonotonic experimental stress-strain dependence. Some numerical results obtained for the stress-strain state of this cylinder under the action of internal pressure are discussed with consideration of the above physical nonlinearity and the finite logarithmic strains.     

Temperature effects on the dynamic response of viscoelastic dielectric elastomer

Viscoelasticity and temperature can significantly affect the performance of a dielectric elastomer. In the current study, we use a thermodynamic model to describe the effect of temperature and viscoelasticity on the electromechanical response undergoing a cyclic electric load by taking into account of the temperature dependent dielectric constant. Because of the significant viscoelasticity in the dielectric elastomer, the deformation and the nominal electric displacement can not keep in phase with the electric field at low frequencies. The results show that the magnitude of the cyclic electromechanical actuation strain increases with the decrease of the temperature and decreases with the increasing frequency, and viscoelasticity can result in significant hysteresis for dielectric elastomers under a relative low temperature and a low frequency.     

Influence of viscoelasticity on drop deformation and orientation in shear flow: Part 1. Stationary states

The influence of matrix and droplet viscoelasticity on the steady deformation and orientation of a single droplet subjected to simple shear is investigated microscopically. Experimental data are obtained in the velocity–vorticity and velocity–velocity gradient plane. A constant viscosity Boger fluid is used, as well as a shear-thinning viscoelastic fluid. These materials are described by means of an Oldroyd-B, Giesekus, Ellis, or multi-mode Giesekus constitutive equation. The drop-to-matrix viscosity ratio is 1.5. The numerical simulations in 3D are performed with a volume-of-fluid algorithm and focus on capillary numbers 0.15 and 0.35. In the case of a viscoelastic matrix, viscoelastic stress fields, computed at varying Deborah numbers, show maxima slightly above the drop tip at the back and below the tip at the front. At both capillary numbers, the simulations with the Oldroyd-B constitutive equation predict the experimentally observed phenomena that matrix viscoelasticity significantly suppresses droplet deformation and promotes droplet orientation. These two effects saturate experimentally at high Deborah numbers. Experimentally, the high Deborah numbers are achieved by decreasing the droplet radius with other parameters unchanged. At the higher capillary and Deborah numbers, the use of the Giesekus model with a small amount of shear-thinning dampens the stationary state deformation slightly and increases the angle of orientation. Droplet viscoelasticity on the other hand hardly affects the steady droplet deformation and orientation, both experimentally and numerically, even at moderate to high capillary and Deborah numbers.     

On Scalar Functionals of Nonlinear tensorial Constitutive Equations in the Theory of Plasticity for Nonisothermal Rectilinear Processes of Deformation

Based on base experiments formulated, methods are proposed to specify the scalar functional in the nonlinear equations that relate generalized stresses and finite strains in the theory of plasticity. The base experiments are conducted and the functionals are specified. It is shown that the nonlinear tensorial constitutive equations can be used to describe a nonisothermal process of deformation along a rectilinear path, which is distinct from the base ones, at high temperatures that cause creep strains     

Law of damage accumulation and fracture criteria in highly filled polymer materials

We present the results of a large series of experiments aimed at the study of laws of damage accumulation and fracture in highly filled polymer materials under loading conditions of various types: monotone, repeated, low- and high-cycle, with varying type of stress state, dynamic (in general, more than 50 programs implemented on specimens from one lot of material). The data obtained in these test allow one to make conclusions about the constitutive role of the attained maximum of strain intensity when estimating the accumulated damage in the process of uniaxial tension by various programs (in particular, an additional cyclic deformation below the preliminary attained strain maximum does not affect the limit values of strain and stress in the subsequent active extension), about the strong influence of the stress state on the deformation and fracture, about the specific features of the nonlinear behavior of the material under the shock loading conditions and its influence on the repeated deformation. All tests are described (with an accuracy acceptable in practical calculations, both with respect to stresses and strains in the process of loading and at the moment of fracture) in the framework of the same model of nonlinear viscoelasticity with the same set of constants. The constants of the proposed model are calculated according to a relatively simple algorithm by using the results of standard uniaxial tension tests with constant values of the strain rate and hydrostatic pressure (each test for 2–3 levels of these parameters chosen from the ranges proposed in applications, each loading lasts until the fracture occurs, and one of the tests contains an intermediate interval of total loading and repeated loading) and one axial shock compression test if there are dynamic problems in the applications. The model is based on the use of the criterion fracture parameter which, in the class of proportional loading processes, is the sum of partial increments of the strain intensity on active segments of the process (where the strain intensity is at its historical maximum) with the form of the stress state and the intensity of strain rates taken into account.     

Nonlinear stress and deformation analysis of thin current-carrying strip shells

The paper analyzes the nonlinear deformation of a current-carrying thin shell in coupled electromagnetic and mechanical fields. The nonlinear magnetoelastic kinetic equations, physical equations, geometric equations, electrodynamic equations, expressions for the Lorentz force of a current-carrying thin shell in a coupled field are given. The normal Cauchy form nonlinear differential equations that include ten basic unknown functions are obtained by the variable replacement method. The difference and quasi-linearization methods are used to reduce the nonlinear magnetoelastic equations to a sequence of quasilinear differential equations that can be solved by discrete orthogonalization. Numerical solutions for the stresses and strains in a current-carrying thin strip shell with two edges simply supported are obtained as an example. The dependence of the stresses and strains in the current-carrying thin strip shell on the electromagnetic parameters is discussed. In a special case, it is shown that the deformation of the shell can be controlled by changing the electromagnetic parameters     

Predicting the nonlinear hereditary viscoelasticity of polymers

A mathematical model for the nonlinear hereditary viscoelasticity of polymer materials is proposed to predict deformation processes of various complexity — from simple relaxation and simple creep to complex deformation-relaxation and reverse relaxation processes with alternative loading and unloading. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 6, pp. 147–157, November–December, 2007.     

含随机参数的非线性黏弹性有限元计算

进行了粗粒土与结构接触面单调和循环加载试验，基于宏细观测量结果, 扩展了 损伤概念以 描述该类接触面在受载过程中的物态演化, 及由于物态演化导致的力学特性从初始状态到最终 稳定状态的连续变化过程. 揭示了接触面损伤的细观物理基础主要是接触面内土的颗粒破碎 和剪切压密这两种物态演化；指出接触面的剪胀体应变可以划分为可逆性和不可逆性剪胀体 应变两部分，其中不可逆性剪胀体应变可作为接触面损伤发展的宏观量度，因此其归一化 形式可作为一种损伤因子的定义；提出了建立粗粒土与结构接触面一种损伤本构关系的基本思路.     

切应力对牛骨蠕变变形的影响研究

骨尤其是湿骨,在恒定载荷作用下会发生蠕变变形。为了确定切应力是否影响骨的蠕变变形,采用对骨薄板试样分别施加集中载荷和均布载荷的方式,测量试样挠度实时的变化曲线。结果显示,在载荷恒定时,骨试样的挠度随时间不断增加,体现了典型的蠕变特性。集中荷载下骨的蠕变变形远大于均布荷载下骨的蠕变变形,湿骨的蠕变位移比干骨高近7倍。分析认为,对试样粘弹性性质的影响不仅有正应力的作用,也有切应力的作用;切应力产生的蠕变变形约为正应力所产生蠕变变形的0.85倍。     

A thermo-viscoelastic–viscoplastic–viscodamage constitutive model for asphaltic materials

A temperature-dependent viscodamage model is proposed and coupled to the temperature-dependent Schapery’s nonlinear viscoelasticity and the temperature-dependent Perzyna’s viscoplasticity constitutive model presented in Abu Al-Rub et al., 2009, Huang et al., in press in order to model the nonlinear constitutive behavior of asphalt mixes. The thermo-viscodamage model is formulated to be a function of temperature, total effective strain, and the damage driving force which is expressed in terms of the stress invariants of the effective stress in the undamaged configuration. This expression for the damage force allows for the distinction between the influence of compression and extension loading conditions on damage nucleation and growth. A systematic procedure for obtaining the thermo-viscodamage model parameters using creep test data at different stress levels and different temperatures is presented. The recursive-iterative and radial return algorithms are used for the numerical implementation of the nonlinear viscoelasticity and viscoplasticity models, respectively, whereas the viscodamage model is implemented using the effective (undamaged) configuration concept. Numerical algorithms are implemented in the well-known finite element code Abaqus via the user material subroutine UMAT. The model is then calibrated and verified by comparing the model predictions with experimental data that include creep-recovery, creep, and uniaxial constant strain rate tests over a range of temperatures, stress levels, and strain rates. It is shown that the presented constitutive model is capable of predicting the nonlinear behavior of asphaltic mixes under different loading conditions.