Formulation and integration of the standard linear viscoelastic solid with fractional order rate laws

A physically sound three-dimensional anisotropic formulation of the standard linear viscoelastic solid with integer or fractional order rate laws for a finite set of the pertinent internal variables is presented. It is shown that the internal variables can be expressed in terms of the strain as convolution integrals with kernels of Mittag–Leffler function type. A time integration scheme, based on the Generalized Midpoint rule together with the Grünwald algorithm for numerical fractional differentiation, for integration of the constitutive response is developed. The predictive capability of the viscoelastic model for describing creep, relaxation and damped dynamic responses is investigated both analytically and numerically. The algorithm and the present general linear viscoelastic model are implemented into the general purpose finite element code Abaqus. The algorithm is then used together with an explicit difference scheme for integration of structural responses. In numerical examples, the quasi-static and damped responses of a viscoelastic ballast material that is subjected to loads simulating the overrolling of a train are investigated.     

Nonlinear Fractional Order Viscoelasticity at Large Strains

In this paper, we formulate a fractional order viscoelastic model for large deformations and develop an algorithm for the integration of the constitutive response. The model is based on the multiplicative split of the deformation gradient into elastic and viscous parts. Further, the stress response is considered to be composed of a nonequilibrium part and an equilibrium part. The viscous part of the deformation gradient (here regarded as an internal variable) is governed by a nonlinear rate equation of fractional order. To solve the rate equation the finite element method in time is used in combination with Newton iterations. The method can handle nonuniform time meshes and uses sparse quadrature for the calculations of the fractional order integral. Moreover, the proposed model is compared to another large deformation viscoelastic model with a linear rate equation of fractional order. This is done by computing constitutive responses as well as structural dynamic responses of fictitious rubber materials.     

Nonlinear Fractional Order Viscoelasticity at Large Strains

In this paper, we formulate a fractional order viscoelastic model for large deformations and develop an algorithm for the integration of the constitutive response. The model is based on the multiplicative split of the deformation gradient into elastic and viscous parts. Further, the stress response is considered to be composed of a nonequilibrium part and an equilibrium part. The viscous part of the deformation gradient (here regarded as an internal variable) is governed by a nonlinear rate equation of fractional order. To solve the rate equation the finite element method in time is used in combination with Newton iterations. The method can handle nonuniform time meshes and uses sparse quadrature for the calculations of the fractional order integral. Moreover, the proposed model is compared to another large deformation viscoelastic model with a linear rate equation of fractional order. This is done by computing constitutive responses as well as structural dynamic responses of fictitious rubber materials.     

分数阶导数粘弹性模型的有限元法

在分析分数阶导数三元件模型理论的基础上,把分数阶导数三元件模型引入有限元模型中,推导出具有分数阶导数三元件本构关系的粘弹性结构动力学有限元格式。同时,应用分数阶导数型粘弹性结构动力学方程的数值算法求解了该有限元格式的数值解。并以二维沥青路面结构为例进行了路面动态粘弹性响应分析。算例分析表明,该方法能够正确有效地进行路面动态粘弹性分析。     

A fractional differential constitutive model for dynamic stress intensity factors of an anti-plane crack in viscoelastic materials

Fractional differential constitutive relationships are introduced to depict the history of dynamic stress inten- sity factors （DSIFs） for a semi-infinite crack in infinite viscoelastic material subjected to anti-plane shear impact load. The basic equations which govern the anti-plane deformation behavior are converted to a fractional wave-like equation. By utilizing Laplace and Fourier integral transforms, the fractional wave-like equation is cast into an ordinary differential equation （ODE）. The unknown function in the solution of ODE is obtained by applying Fourier transform directly to the boundary conditions of fractional wave-like equation in Laplace domain instead of solving dual integral equations. Analytical solutions of DSIFs in Laplace domain are derived by Wiener-Hopf technique and the numerical solutions of DSIFs in time domain are obtained by Talbot algorithm. The effects of four parameters α, β, b1, b2 of the fractional dif- ferential constitutive model on DSIFs are discussed. The numerical results show that the present fractional differential constitutive model can well describe the behavior of DSIFs of anti-plane fracture in viscoelastic materials, and the model is also compatible with solutions of DSIFs of anti-plane fracture in elastic materials.     

A new fractional derivative model for linearly viscoelastic materials and parameter identification via genetic algorithms

In this study, a new nested model consisting of springs and “spring pots” is proposed to better simulate the viscoelastic behavior of polymeric damping materials in the frequency domain. First, the one-dimensional constitutive equation that consists of ten parameters is derived. The dynamical mechanical properties, which are the storage modulus, the loss modulus, and the loss factor, are obtained from this equation. Then, the low- and high-frequency behavior of this model is investigated. Moreover, a new methodology to identify the unknown parameters that appear in the fractional derivative model, from the experimental data for the Wicket plot, by using genetic algorithms (GAs) is presented. This approach does not require shifting of the experimental data; therefore, possible errors that may arise are eliminated. The new model is fitted to experimental data for several polymeric damping materials that exist in the literature, in order to verify its success. The results are presented in a graphical form with comparison to the already existing models.     

A numerical method for fractional integral with applications

IntroductionThefractionalcalculushasalonghistoryandthereareamassofworkstodiscussthefractionalderivativesandfractionalintegralswitharbitrary (realorcomplex)order[1- 3 ].Thefractionalcalculushasawideapplicationbackground ,especiallyinthefieldsofchemistry ,electromagnetics,materialscienceandmechanics.Forexample,Gement[4 ]proposedthefractionalderivativeconstitutivemodelsofaviscoelasticmaterialatfirst.Themodelshavereceivedincreasingattention[5 - 7].Onlyafewparametersarecontainedinthemodelsandthemo…     

On modeling the micro-indentation response of an amorphous polymer

Interest in instrumented indentation experiments as a means to estimate mechanical properties has grown rapidly in recent years. Although numerous nano/micro-indentation experimental studies on polymeric materials have been reported in the literature, a corresponding methodology for extracting material property information from the experimental data does not exist. This situation for polymeric materials exists primarily because baseline numerical analyses of sharp indentation using appropriate large deformation constitutive models for the nonlinear viscoelastic–plastic response of these materials appear not to have been previously reported in the literature. An existing, widely used theory for amorphous polymers (e.g. [Boyce, M., Parks, D., Argon, A.S., 1988. Large inelastic deformation of glassy polymers. Part 1: Rate dependent constitutive model. Mechanics of Materials 7, 15–33; Arruda, E.M., Boyce, M.C., 1993. Evolution of plastic anisotropy in amorphous polymers during finite straining. International Journal of Plasticity 9, 697–720]) has been recently found to lack sufficient richness to enable one to quantitatively reproduce the major features of the indentation load-versus-depth curves for some common amorphous polymers [Gearing, B.P., 2002. Constitutive equations and failure criteria for amorphous polymeric solids. Ph.D. thesis, Massachusetts Institute of Technology].This study develops a new continuum model for the viscoelastic–plastic deformation of amorphous polymeric solids. We have applied the constitutive model to capture salient features of the mechanical response of the amorphous polymeric solid poly(methyl methacrylate) (PMMA) at ambient temperature and stress states under which this material does not exhibit crazing. We have conducted compression-tension strain-controlled experiments, as well as stress-controlled compression-creep experiments, and these experiments are used to calibrate the material parameters in the constitutive model for PMMA.We have implemented our constitutive model in a finite-element computer program, and using this finite-element program we have simulated micro-indentation experiments on PMMA. We show that our constitutive model and finite element simulations reproduce the experimentally-measured indentation load-versus-depth response with reasonable accuracy.     

分数算子描述的黏弹性体力学问题数值方法

讨论由黎曼-刘维尔 (Riemann-Liouville)分数导数描述的黏弹性体力学问题的数 值方法. 该方法利用黎曼-刘维尔分数 导数定义中核函数的特性，并结合被积函数在单步中的逼近以及Newmark型数值法，建立 了分数导数计算公式. 算例表明，该算法具有收敛快、精度高、稳定性好和易于 应用和改进的优点. 在对动态系统的瞬态响应分析和有限元分析格式中，算法都获得了 满意的结果.     

FE formulation for the viscoelastic body modeled by fractional constitutive law

This paper presents finite element (FE) formulation of the viscoelastic materials described by fractional constitutive law. The time-domain three-dimensional constitutive equation is constructed. The FE equations are set up by equations are solved by numerical integration method. The numerical algorithm developed by the authors for Liouville-Riemann's fractional derivative was adopted to formulate FE procedures and extended to solve the more general case of the hereditary integration. The numerical examples were given to show the correctness and effectiveness of the integration algorithm. The project supported by the Ministry of Education of China for the returned overseas Chinese scholars     

Time-dependent behavior of passive skeletal muscle

An isotropic three-dimensional nonlinear viscoelastic model is developed to simulate the time-dependent behavior of passive skeletal muscle. The development of the model is stimulated by experimental data that characterize the response during simple uniaxial stress cyclic loading and unloading. Of particular interest is the rate-dependent response, the recovery of muscle properties from the preconditioned to the unconditioned state and stress relaxation at constant stretch during loading and unloading. The model considers the material to be a composite of a nonlinear hyperelastic component in parallel with a nonlinear dissipative component. The strain energy and the corresponding stress measures are separated additively into hyperelastic and dissipative parts. In contrast to standard nonlinear inelastic models, here the dissipative component is modeled using an evolution equation that combines rate-independent and rate-dependent responses smoothly with no finite elastic range. Large deformation evolution equations for the distortional deformations in the elastic and in the dissipative component are presented. A robust, strongly objective numerical integration algorithm is used to model rate-dependent and rate-independent inelastic responses. The constitutive formulation is specialized to simulate the experimental data. The nonlinear viscoelastic model accurately represents the time-dependent passive response of skeletal muscle.     

Finite deformation of a polymer: experiments and modeling

The response of a polymer (polytetrafluoroethylene) to quasi-static and dynamic loading is determined and modeled. The polytetrafluoroethylene is extremely ductile and highly nonlinear in elastic as well as plastic behaviors including elastic unloading. Constitutive model developed earlier by Khan, Huang and Liang (KHL) is extended to include the responses of polymeric materials. The strain rate hardening, creep, and relaxation behaviors of polytetrafluoroethylene were determined through extensive experimental study. Based on the observation that both viscoelastic and viscoplastic deformation of polytetrafluoroethylene are time dependent and nonlinear, a phenomenalogical viscoelasto–plastic constitutive model is presented by a series connection of a viscoelastic deformation module (represented by three elements standard solid spring dashpot model), and a viscoplastic deformation module represented by KHL model. The KHL module is affected only when the stress exceeds the initial yield stress. The comparison between the predictions from the extended model and experimental data for uniaxial static and dynamic compression, creep and relaxation demonstrate that the proposed constitutive model is able to represent the observed time dependent mechanical behavior of polytetrafluoroethylene polytetrafluoroethylene qualitatively and quantitatively.     

分数导数型本构关系描述粘弹性梁的振动分析

本文研究粘弹性梁在周期激励作用下的受迫振动问题。梁的材料满足Kelvin-Volgt分数导数型本构关系。基于动力学方程、本构关系和应变－位移关系建立了小变形粘弹性梁的振动方程。采用分离变量法分析粘弹性梁的自由振动,导出模态坐标满足的常微分－积分方程和模态函数满足的常微分方程,对于两端简支的截面梁给出了固有频率和模态函数。对于简谐激励作用下粘弹性梁的受迫振动,利用模态叠加得到了稳态响应。最后给出数值算例说明本文方法的应用。     

Modeling and simulation of three‐dimensional extrusion swelling of viscoelastic fluids with PTT,Giesekus and FENE‐P constitutive models

The investigation of the extrusion swelling mechanism of viscoelastic fluids has both scientific and industrial interest. However, it has been traditionally difficult to afford theoretical and experimental researches to this problem. The numerical methodology based on the penalty finite element method with a decoupled algorithm is presented in the study to simulate three‐dimensional extrusion swelling of viscoelastic fluids flowing through out of a circular die. The rheological responses of viscoelastic fluids are described by using three kinds of differential constitutive models including the Phan‐Thien Tanner model, the Giesekus model, and the finite extensible nonlinear elastic dumbbell with a Peterlin closure approximation model. A streamface‐streamline method is introduced to adjust the swelling free surface. The calculation stability is improved by using the discrete elastic‐viscous split stress algorithm with the inconsistent streamline‐upwind scheme. The essential flow characteristics of viscoelastic fluids are predicted by using the proposed numerical method, and the mechanism of swelling phenomenon is further discussed.Copyright © 2013 John Wiley & Sons, Ltd.     

一种适合橡胶类材料的非线性粘弹性本构模型

借助非线性流变模型建立大变形情况非线性粘弹材料的本构关系,考虑到大多数橡胶类材料具有的几乎不可压缩性,以及体积响应和剪切响应的流变性能不同,将变形梯度乘法分解为等容部分和体积变形部分,给出了一种适合橡胶类材料的非线性粘弹性本构模型,并模拟了粘滞效应。对于极快或极慢的过程,该模型退化为橡胶弹性理论;在小变形情况下退化为经典的广义Maxwell粘弹性材料。模型与热力学第二定律相容,适合于大规模数值分析。     

On the large deflections of linear viscoelastic beams

The quasi-static response and the stored and dissipated energies due to large deflections of a slender inextensible beam made of a linear viscoelastic material and subjected to a time-dependent inclined concentrated load at the free end are investigated. The beam cross-section is considered prismatic, the self-weight is disregarded and the material is initially stress free. The set of four first-order non-linear partial integro-differential equations obtained from geometrical compatibility, equilibrium of forces and moments, and linear viscoelastic constitutive relation is numerically solved using a one-parameter shooting method combined with a fourth-order Runge-Kutta algorithm. An analytical expression is derived to divide the energy supplied by the external load into conserved and dissipated parts. For the case study presented, a three-parameter solid linear viscoelastic constitutive model is employed and a step load is applied. The variables are made non-dimensional, so four parameters govern the problem: the ratio between the final and initial relaxation moduli, the load magnitude, the angle of inclination and the unloading time. A finite-element model is also performed to compare and validate the analytical and numerical formulations. Results are presented for encastré curvature and tip displacement versus time, geometrical configuration, load versus tip displacement, total work done by the external force, stored and dissipated energies versus time, energy per unit length along arc length for three times and versus time for two beam sections.     

On the behavior of a three-dimensional fractional viscoelastic constitutive model

In this paper a three-dimensional isotropic fractional viscoelastic model is examined. It is shown that if different time scales for the volumetric and deviatoric components are assumed, the Poisson ratio is time varying function; in particular viscoelastic Poisson ratio may be obtained both increasing and decreasing with time. Moreover, it is shown that, from a theoretical point of view, one-dimensional fractional constitutive laws for normal stress and strain components are not correct to fit uniaxial experimental test, unless the time scale of deviatoric and volumetric are equal. Finally, the model is proved to satisfy correspondence principles also for the viscoelastic Poisson’s ratio and some issues about thermodynamic consistency of the model are addressed.     

Finite Element Formulation of Viscoelastic Constitutive Equations Using Fractional Time Derivatives

Fractional time derivatives are used to deduce a generalization ofviscoelastic constitutive equations of differential operator type. Theseso-called fractional constitutive equations result in improvedcurve-fitting properties, especially when experimental data from longtime intervals or spanning several frequency decades need to be fitted.Compared to integer-order time derivative concepts less parameters arerequired. In addition, fractional constitutive equations lead to causalbehavior and the concept of fractional derivatives can be physicallyjustified providing a foundation of fractional constitutive equations.First, three-dimensional fractional constitutive equations based onthe Grünwaldian formulation are derived and their implementationinto an elastic FE code is demonstrated. Then, parameter identificationsfor the fractional 3-parameter model in the time domain as well as inthe frequency domain are carried out and compared to integer-orderderivative constitutive equations. As a result the improved performanceof fractional constitutive equations becomes obvious. Finally, theidentified material model is used to perform an FE time steppinganalysis of a viscoelastic structure.     

Time domain modeling of damping using anelastic displacement fields and fractional calculus

A fractional derivative model of linear viscoelasticity based on the decomposition of the displacement field into an anelastic part and elastic part is developed. The evolution equation for the anelastic part is then a differential equation of fractional order in time. By using a fractional order evolution equation for the anelastic strain the present model becomes very flexible for describing the weak frequency dependence of damping characteristics. To illustrate the modeling capability, the model parameters are fit to available frequency domain data for a high damping polymer. By studying the relaxation modulus and the relaxation spectrum the material parameters of the present viscoelastic model are given physical meaning. The use of this viscoelastic model in structural modeling is discussed and the corresponding finite element equations are outlined, including the treatment of boundary conditions. The anelastic displacement field is mathematically coupled to the total displacement field through a convolution integral with a kernel of Mittag–Leffler function type. Finally a time step algorithm for solving the finite element equations are developed and some numerical examples are presented.