三种分形和分数阶导数阻尼振动模型的比较研究

标准的整数阶导数方程不能准确描述粘弹性材料的记忆性参考文献[1]和阻尼的分数次幂频率依赖[2],因此分形导数、分数阶导数及正定分数阶导数被用于描述粘弹性介质中的阻尼振动.该文通过分析模型和数值模拟,比较了三种模型描述的振动过程.结果显示,当p小于约O.75或大于约1.9时(p为非整数阶导数的阶数),分形导数模型衰减最快;当P大于约0.75且小于约1.9时,正定分数阶导数模型衰减最快,衰减最慢的分别为分数阶导数模型(p<1)和分形导数模型(p>1).且正定分数阶导数模型衰减快于分数阶导数模型,当p接近2时,两种模型较为相近.     

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.     

Role of Prehistories in the Initial Value Problems of Fractional Viscoelastic Equations

The fractional viscoelastic equation (FVE), which is a second-order differential equation with fractional derivatives describing the dynamical behavior of a single-degree-of-freedom viscoelastic oscillator, is considered. Some viscoelastic damped mechanical systems may be described by FVEs. However, FVEs with conventional nonzero initial values cannot generally be solved. In this paper, the prehistories of the unknown functions before the initial times, referred to as the initial functions, are taken into account to solve FVEs. Mathematically, appropriate initial functions are essential for unique solutions of FVEs. Physically, the initial functions reflect the processes of giving the initial values. FVEs are solved for some initial functions both by analytical and numerical methods. The initial functions affect the solutions of FVEs. It is discussed how the solutions depend on the initial functions. Implication of the solutions to viscoelastic materials will be discussed.     

Role of Prehistories in the Initial Value Problems of Fractional Viscoelastic Equations

The fractional viscoelastic equation (FVE), which is a second-order differential equation with fractional derivatives describing the dynamical behavior of a single-degree-of-freedom viscoelastic oscillator, is considered. Some viscoelastic damped mechanical systems may be described by FVEs. However, FVEs with conventional nonzero initial values cannot generally be solved. In this paper, the prehistories of the unknown functions before the initial times, referred to as the initial functions, are taken into account to solve FVEs. Mathematically, appropriate initial functions are essential for unique solutions of FVEs. Physically, the initial functions reflect the processes of giving the initial values. FVEs are solved for some initial functions both by analytical and numerical methods. The initial functions affect the solutions of FVEs. It is discussed how the solutions depend on the initial functions. Implication of the solutions to viscoelastic materials will be discussed.     

Excluded volume effects and fractional viscoelasticity in polymers

Meccanica - The excluded volume effect is added to a fractional viscoelastic model for modeling fractal polymers. This reveals a physical connection between the fractional time derivative, fractal...     

Characterizing the creep of viscoelastic materials by fractal derivative models

In this paper, we make the first attempt to apply the fractal derivative to modeling viscoelastic behavior. The methodology of scaling transformation is utilized to obtain the creep modulus and relaxation compliance for the proposed fractal Maxwell and Kelvin models. Comparing with the fractional derivatives reported in the literature, the fractal derivative as a local operator has lower calculation costs and memory storage requirements. Moreover, numerical results show that the proposed fractal models require fewer parameters, have simpler mathematical expression and result in higher accuracy than the classical integer-order derivative models. Results further confirm that the proposed fractal models can characterize the creep behavior of viscoelastic 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.     

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.     

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

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

Non-linear viscoelastic behavior of polymer melts interpreted by fractional viscoelastic model

Very recently, researchers dealing with constitutive law pertinent viscoelastic materials put forward the successful idea to introduce viscoelastic laws embedded with fractional calculus, relating the stress function to a real order derivative of the strain function. The latter consideration leads to represent both, relaxation and creep functions, through a power law function. In literature there are many papers in which the best fitting of the peculiar viscoelastic functions using a fractional model is performed. However there are not present studies about best fitting of relaxation function and/or creep function of materials that exhibit a non-linear viscoelastic behavior, as polymer melts, using a fractional model. In this paper the authors propose an advanced model for capturing the non-linear trend of the shear viscosity of polymer melts as function of the shear rate. Results obtained with the fractional model are compared with those obtained using a classical model which involves classical Maxwell elements. The comparison between experimental data and the theoretical model shows a good agreement, emphasizing that fractional model is proper for studying viscoelasticity, even if the material exhibits a non-linear behavior.     

层状分数阶黏弹性饱和地基与梁共同作用的时效研究

饱和地基与梁共同作用问题的研究在力学领域及工程界都具有重要意义.采用分数阶Merchant模型研究饱和地基的流变固结,该模型比常用整数阶黏弹性模型更能精确反映地基的时变特征.基于层状正交各向异性黏弹性饱和地基的固结解答,采用有限元法与边界元法耦合的方法,研究梁与分数阶黏弹性饱和地基的共同作用问题.依据Timoshenk...     

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.     

Wave propagation in a fractional viscoelastic Andrade medium: Diffusive approximation and numerical modeling

This study focuses on the numerical modeling of wave propagation in fractionally-dissipative media. These viscoelastic models are such that the attenuation is frequency-dependent and follows a power law with non-integer exponent within certain frequency regimes. As a prototypical example, the Andrade model is chosen for its simplicity and its satisfactory fits of experimental flow laws in rocks and metals. The corresponding constitutive equation features a fractional derivative in time, a non-local-in-time term that can be expressed as a convolution product whose direct implementation bears substantial memory cost. To circumvent this limitation, a diffusive representation approach is deployed, replacing the convolution product by an integral of a function satisfying a local time-domain ordinary differential equation. An associated quadrature formula yields a local-in-time system of partial differential equations, which is then proven to be well-posed. The properties of the resulting model are also compared to those of the Andrade model. The quadrature scheme associated with the diffusive approximation, and constructed either from a classical polynomial approach or from a constrained optimization method, is investigated. Finally, the benefits of using the latter approach are highlighted as it allows to minimize the discrepancy with the original model. Wave propagation simulations in homogeneous domains are performed within a split formulation framework that yields an optimal stability condition and which features a joint fourth-order time-marching scheme coupled with an exact integration step. A set of numerical experiments is presented to assess the overall approach. Therefore, in this study, the diffusive approximation is demonstrated to provide an efficient framework for the theoretical and numerical investigations of the wave propagation problem associated with the fractional viscoelastic medium considered.     

Constitutive equation with fractional derivatives for the generalized UCM model

A fundamental problem on the constitutive equation with fractional derivatives for the generalized upper convected Maxwell model (UCM) is studied. The existing investigations on the constitutive equation are reviewed and their limitations or deficiencies are highlighted. By utilizing the convected coordinates approach, a mathematically rigorous constitutive equation with fractional derivatives for the generalized UCM model is proposed, which has an explicit expression for the stress tensor. This model can be reduced to the linear generalized Maxwell model with fractional derivatives, the UCM model and some other existing models. In addition, the rheological properties of this proposed model in the start-up of simple shear and elongation flows are investigated. It is shown that this generalized UCM model can describe the various stress evolution processes and the strain hardening effect of the viscoelastic fluids.     

Understanding thixotropic and antithixotropic behavior of viscoelastic micellar solutions and liquid crystalline dispersions. I. The model

A simple model consisting of the Upper Convected Maxwell constitutive equation and a kinetic equation for destruction and construction of structure, first proposed by Fredrickson in 1970, is used here to reproduce the complex rheological behavior of viscoelastic systems that also exhibit thixotropy and rheopexy under shear flow. The model requires five parameters that have physical significance and that can be estimated from rheological measurements. Several steady and unsteady flow situations were analyzed with the model. The model predicts creep behavior, stress relaxation and the presence of thixotropic loops when the sample is subjected to transient stress cycles. Such behavior has been observed with surfactant-based solutions and dispersions. The role of the characteristic time for structure built up, λ, in the extent and shape of the thixotropic loops is demonstrated.     

广义Maxwell黏弹性流体在两平板间的非定常流动

将分数阶微积分运算引入Maxwell黏弹性流体的本构方程，研究了黏弹性流体在两平板问的非定常流动．对于广义Maxwell黏弹性流体的分数阶导数模型，导出了对时间具有分数阶导数的特殊运动方程，利用分数阶微积分的Laplace变换理论，得到了流动的解析解．     

半封闭分数导数粘弹性衬砌隧道的频域响应

混凝土衬砌具有粘弹性性质,以往的经典Kelvin模型、弹性理论和壳体理论都不能刻画其蠕变的全过程。本文基于饱和多孔介质理论,在频率域研究了轴对称荷载和流体压力作用下饱和粘弹性土中半封闭分数导数型衬砌隧洞的稳态动力响应。在引入隧洞部分透水边界条件的基础上,通过分数阶导数粘弹性模型描述衬砌的应力—位移本构关系,并利用衬砌内边界以及接触面的连续性条件,得到了饱和土和衬砌的应力、位移和孔压解答。考察了分数导数阶数、材料参数以及衬砌和土体相对渗透系数的影响。研究表明：分数导数阶数对系统响应影响较大,且依赖于衬砌的材料参数。另外,相对渗透系数对系统响应的影响很大。     

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…