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

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

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

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

A phenomenological constitutive law for the nonlinear viscoelastic behaviour of epoxy resins in the glassy state

We propose a constitutive model to describe the nonlinear viscoelasticity of epoxy resins in the glassy state. Experimental tests have shown that nonlinear viscoelasticity rules the cyclic behaviour of epoxy resins before the material strength is reached. The proposed model allows the simulation of this cyclic behaviour and, in particular, of the flex which characterises the stress–strain curve upon unloading. The consistent integration algorithm, based on the Central Difference Scheme, of the proposed constitutive law is given. As an example, the model is successfully used to numerically homogenize the cyclic behaviour of hollow sphere-filled epoxy resins (i.e., syntactic foams) by means of unit cell models.     

Plane surface suddenly set in motion in a viscoelastic fluid with fractional Maxwell model

The fractional calculus approach in the constitutive relationship model of viscoelastic fluid is introduced. The flow near a wall suddenly set in motion is studied for a non-Newtonian viscoelastic fluid with the fractional Maxwell model. Exact solutions of velocity and stress are obtained by using the discrete inverse Laplace transform of the sequential fractional derivatives. It is found that the effect of the fractional orders in the constitutive relationship on the flow field is significant. The results show that for small times there are appreciable viscoelastic effects on the shear stress at the plate, for large times the viscoelastic effects become weak. The project supported by the National Natural Science Foundation of China (10002003), Foundation for University Key Teacher by the Ministry of Education, Research Fund for the Doctoral Program of Higher Education     

On constitutive equations for fiber-reinforced nonlinearly viscoelastic solids

In this paper, we are interested in developing constitutive equations for fiber-reinforced nonlinearly viscoelastic bodies, in particular for transversely isotropic nonlinearly viscoelastic solids. It follows from results in the theory of algebraic invariants that constitutive equations for such materials can be expressed in terms of functions of 18 independent invariants associated with deformation and fiber orientation. These invariants are analyzed, and we obtain restrictions such as positivity of some of them.     

Quasi-static and dynamical analysis for viscoelastic Timoshenko beam with fractional derivative constitutive relation

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｆｒａｃｔｉｏｎａｌｄｅｒｉｖａｔｉｖｅｃｏｎｓｔｉｔｕｔｉｖｅｍｏｄｅｌｓｏｆａｖｉｓｃｏｅｌａｓｔｉｃｍａｔｅｒｉａｌｗｅｒｅｐｒｏｐｏｓｅｄｂｙＧｅｍｅｎｔａｔｆｉｒｓｔｉｎ 1 93 0’ｓ[1].Ｓｉｎｃｅ 1 980’ｓ,ｔｈｅｍｏｄｅｌｓｈａｖｅｒｅｃｅｉｖｅｄｉｎｃｒｅａｓｉｎｇａｔｔｅｎｔｉｏｎ[2 ,3].Ｏｎｌｙａｆｅｗｐａｒａｍｅｔｅｒｓａｒｅｃｏｎｔａｉｎｅｄｉｎｔｈｅｍｏｄｅｌｓａｎｄｔｈｅｍｏｄｅｌｓｃａｎｄｅｓｃｒｉｂｅｔｈｅｍｅｃｈａｎｉｃａｌｃｈａｒａｃ…     

A micromechanical analysis of the nonlinear elastic and viscoelastic constitutive relation of a polymer filled with rigid particles

Micromechanical theory is applied to study the nonlinear elastic and viscoelastic constitutive relations of polymeric matrix filled with high rigidity solid particles. It is shown that Eshelby's method can be extended to the case of nonlinear matrix and Eshelby's tensor still exists provided that Poisson's ratio of the nonlinear matrix assumes constant value in deforming process and the rigidity of elastic filling particles is much higher than that of the matrix. A new method for averaging process is proposed to overcome the difficulty that occured in applying the ordinary equivalent inclusion method or the self-consistant method to nonlinear matrices. A rather simple constitutive equation is obtained finally and the strengthening effect of solid particles to composites is investigated. The work supported by the LNM, Institute of Mechanics, Chinese Academy of Sciences and by the National Natural Science Foundation of China     

Study on the constitutive equation with fractional derivative for the viscoelastic fluids – Modified Jeffreys model and its application

Based on the classic linear viscoelastic Jeffreys model, a modified Jeffreys model is suggested. The corresponding five-parameter equation with fractional derivatives of different orders of the stress and rate of strain is stated and the characteristic material functions of the linear viscoelasticity theory, such as the dynamic moduli, are derived. The comparison between the measured dynamic moduli of Sesbania gel and xanthan gum and the theoretical predictions of the proposed phenomenological model shows an excellent agreement. Received: 26 August 1997 Accepted: 26 May 1998     

基于分数导数本构的粘弹性土层—隧洞衬砌系统的稳态动力响应

将混凝土衬砌材料视为具有分数导数本构的粘弹性体,在频率域研究了深埋圆形隧洞粘弹性土的稳态动力响应。利用衬砌的内边界条件和混凝土衬砌与土体界面处的连续性条件,得到了粘弹性土体和衬砌稳态振动的应力和位移解析解。进行了隧洞衬砌厚度、土体阻尼比、本构阶数、材料参数比等对土体的位移和应力幅值影响的算例分析,结果表明:材料参数比对系统动力特性有较大影响,随着材料参数比的增加,位移和应力幅值减小;随着土体的阻尼比、衬砌厚度的增加,位移和应力幅值减小;当材料参数比Tσ/Tε=3时,随着分数导数阶数增大,响应幅值减小。     

分数阶黏弹性土层中分数阶三维轴对称桩的竖向振动

将桩基和土体视为三维连续介质,桩基和土体的应力一应变关系采用分数阶黏弹性模型描述。在三维轴对称情况下,利用三维弹性理论和连续介质力学理论,运用分离变量法和分数阶导数的性质,得到了分数阶黏弹性土层中分数阶黏弹性桩基的三维轴对称解;并分析了相关参数对桩顸动态刚度和等效阻尼的影响。研究结果表明：与土体相比,桩基的相关参量对桩顶复刚度的影响较大;桩基和土体的密度比、模量比对桩顶复刚度都有较大的影响。     

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…     

Theoretical analysis of microscopic oil displacement mechanism by viscoelastic polymer solution

The microscopic oil displacement mechanism in viscoelastic polymer flooding is theoretically analyzed with mechanical method. The effects of viscoelasticity of polymer solution on such three kinds of residual oil as in pore throat, in sudden expansion pore path, and in dead end are analyzed. Results show that the critical radius of mobile residual oil for viscoelastic polymer solution is larger than that for viscous polymer solution, which makes the oil that is immobile in viscous polymer flooding displaced under the condition of viscoelastic polymer solution. The viscous polymer solution hardly displaces the oil in dead ends. However, when the effect of viscoelasticity is considered, the residual oil in sudden expansion pore paths and dead ends can be partly displaced. A dimensionless parameter is suggested to denote the relative dominance of gravity and capillary pressure. The larger the dimensionless parameter, the more accurate the increment expressions.     

On thermodynamically consistent constitutive equations for fiber-reinforced nonlinearly viscoelastic solids with application to biomechanics

In this paper, we are interested in developing thermodynamically consistent constitutive equations for fiber-reinforced nonlinearly viscoelastic bodies, in particular for transversely isotropic nonlinearly viscoelastic solids, in isothermal processes. It follows from results in the theory of algebraic invariants that constitutive equations for such materials can be expressed in terms of functions of 18 independent invariants associated with deformation and fiber orientation: 10 of them are isotropic invariants and 8 of them are associated with the deformation and the orientation of the fiber. Among the 8 anisotropic invariants just 6 are related to the viscoelastic response. The terms in the Cauchy stress tensor associated to these 6 invariants are analyzed with respect to thermodynamical consistency, and we obtain restrictions for the corresponding constitutive coefficients. This framework is applied to viscoelastic potentials within the context of biomaterials.     

On the application of fractional calculus for the formulation of viscoelastic Reddy beam

Meccanica - The focus of the current work is to present the bending analysis of visco-elastic beams based on Reddy’s third-order shear deformation theory. Fractional calculus is taken into...     

Limitations of the Phan-Thien non-linear network viscoelastic model

The failure of the Phan-Thien non-linear network viscoelastic model to predict realistic stress-strain relationships under certain conditions is illustrated in two examples. Care in the application of the model when applied in situations where high shear rates are expected is indicated.     

An unsymmetric constitutive equation for anisotropic viscoelastic fluid

A continuum constitutive theory of corotational derivative type is developed for the anisotropic viscoelastic fluid–liquid crystalline (LC) polymers. A concept of anisotropic viscoelastic simple fluid is introduced. The stress tensor instead of the velocity gradient tensor D in the classic Leslie–Ericksen theory is described by the first Rivlin–Ericksen tensor A and a spin tensor W measured with respect to a co-rotational coordinate system. A model LCP-H on this theory is proposed and the characteristic unsymmetric behaviour of the shear stress is predicted for LC polymer liquids. Two shear stresses thereby in shear flow of LC polymer liquids lead to internal vortex flow and rotational flow. The conclusion could be of theoretical meaning for the modern liquid crystalline display technology. By using the equation, extrusion–extensional flows of the fluid are studied for fiber spinning of LC polymer melts, the elongational viscosity vs. extension rate with variation of shear rate is given in figures. A considerable increase of elongational viscosity and bifurcation behaviour are observed when the orientational motion of the director vector is considered. The contraction of extrudate of LC polymer melts is caused by the high elongational viscosity. For anisotropic viscoelastic fluids, an important advance has been made in the investigation on the constitutive equation on the basis of which a series of new anisotropic non-Newtonian fluid problems can be addressed. The project supported by the National Natural Science Foundation of China (10372100, 19832050) (Key project). The English text was polished by Yunming Chen.     

A “quasi-elastic” affine formulation for the homogenised behaviour of nonlinear viscoelastic polycrystals and composites

The derivation of the overall behaviour of nonlinear viscoelastic (or rate-dependent elastoplastic) heterogeneous materials requires a linearisation of the constitutive equations around uniform per phase stress (or strain) histories. The resulting Linear Comparison Material (LCM) has to be linear thermoviscoelastic to fully retain the viscoelastic nature of phase interactions. Instead of the exact treatment of this LCM (i.e., correspondence principle and inverse Laplace transforms) as proposed by the “classical” affine formulation, an approximate treatment is proposed here. First considering Maxwellian behaviour, comparisons for a single phase as well as for two-phase materials (with “parallel” and disordered morphologies) show that the “direct inversion method” of Laplace transforms, initially proposed by Schapery (1962), has to be adapted to fit correctly exact responses to creep loading while a more general method is proposed for other loading paths. When applied to nonlinear viscoelastic heterogeneous materials, this approximate inversion method gives rise to a new formulation which is consistent with the classical affine one for the steady-state regimes. In the transient regime, it leads to a significantly more efficient numerical resolution, the LCM associated to the step by step procedure being no more thermoviscoelastic but thermoelastic. Various comparisons for nonlinear viscoelastic polycrystals responses to creep as well as relaxation loadings show that this “quasi-elastic” formulation yields results very close to classical affine ones, even for high contrasts.     

An empirical constitutive law for concentrated colloidal suspensions in the approach of the glass transition

Concentrated, non-crystallizing colloidal suspensions in their approach of the glass state exhibit distinct dynamics patterns. These patterns suggest a powerlaw rheological constitutive model for near-glass viscoelasticity, as presented here. The rheological parameters used for this model originate in the mode-coupling theory. The proposed constitutive model provides explicit expressions for the steady shear viscosity, the steady normal stress coefficient, the modulus-compliance relation, and the α peak of G″. The relaxation pattern distinctly differs from gelation.     

Accurate and straightforward symplectic approach for fracture analysis of fractional viscoelastic media

An accurate and straightforward symplectic method is presented for the fracture analysis of fractional two-dimensional(2D) viscoelastic media. The fractional Kelvin-Zener constitutive model is used to describe the time-dependent behavior of viscoelastic materials. Within the framework of symplectic elasticity, the governing equations in the Hamiltonian form for the frequency domain(s-domain) can be directly and rigorously calculated. In the s-domain, the analytical solutions of the displacement ...