广义黏弹组合模型的等效性及其基本性质

对广义黏弹组合模型的等效性及其基本性质进行了研究,结果表明:广义Maxwell模型由于并联结构,其蠕变柔量很难得到. 但提出了求解广义Maxwell模型蠕变柔量的方法,给出了其具体表达式,且在此基础上证明了广义Maxwell模型与Kelvin链的等效性,建立了这两个模型物理参数间的转换关系式. 证明了广义黏弹组合模型的一个基本性质:当将模型的松弛时间谱和延迟时间谱由大到小顺序排列时,松弛时间与延迟时间互不相等,且相互交织,两相邻松弛时间中间有且仅有一个延迟时间,同时,两相邻延迟时间中间有且仅有一个松弛时间;当两者同阶相比时,延迟时间总是大于松弛时间. 这一基本性质明确了使用广义黏弹组合模型来描述现实中某种特定材料的黏弹性行为时,该材料必须具备的基本条件,因此,它可作为这类流变模型在工程应用中的一个实用判据. Wiechert模型和广义中村模型、广义Jeffreys模型和广义N-K模型、Maxwell链和广义Kelvin模型之间的等效性可作为特例.最后,实例验证了所提出的求解广义黏弹组合模型蠕变柔量方法及其基本性质.      

非定常微分型黏弹性本构模型

在应力作用下, 材料的力学参数随着微观结构的变化而变化, 需要考虑参数的时间效应. 利用黏滞系数随时间变化的黏性元件, 构造出非定常Maxwell模型、非定常Kelvin模型和非定常Zener模型. 求解非定常模型的微分型本构方程得到它们的松弛模量、蠕变柔量和卸载方程. 结果表明, 可以把常见的经验松弛函数和经验蠕变函数视为非定常微分型本构模型.      

关于黏弹性材料的广义Maxwell模型

采用流变力学分析黏弹性材料的流变特性时，常要用到广义Maxwell模型 表达的应力松弛模量. 而从试验中获得的应力松弛模量，其表达式常为 Kohlrausch-William-Watts function(KWW函数)形式. 通过把KWW函数和广义Maxwell模型的拟合问题转化为两 矩阵相等的求解问题后，又把两矩阵的相等等价于两矩阵差值向量的一阶范数为无穷小的问 题，并通过引入广义逆矩阵，求得两矩阵差值向量的一阶范数的最小值，最后以一阶范数的 最小值为目标函数，松弛时间为约束条件，利用单纯形法对两矩阵差值向量的一阶范数的最 小值优化，从而提出了一种针对黏弹材料的KWW函数与广义Maxwell模型转换的计算方法. 借助于MATLAB软件，实现了对黏弹材料的广义Maxwell模型的拟合.     

考虑应变率的广义压电热弹理论及其应用

工程中大量材料的形变介于弹性与黏性之间, 既具有弹性固体特性, 又具有黏性流体特点, 即为黏弹性. 黏弹性使得材料出现很多力学松弛现象, 如应变松弛、滞后损耗等行为. 在研究受热载荷作用的多场耦合问题的瞬态响应时, 考虑此类问题中的热松弛和应变松弛现象, 对准确描述其瞬态响应尤为重要. 针对广义压电热弹问题的瞬态响应, 尽管已有学者建立了考虑热松弛的广义压电热弹模型, 但迄今, 尚未计入应变松弛. 本文中, 考虑到材料变形时的应变松弛, 通过引入应变率, 在Chandrasekharaiah广义压电热弹理论的基础之上, 经拓展, 建立了考虑应变率的广义压电热弹理论. 借助热力学定律, 给出了理论的建立过程并得到了相应的状态方程及控制方程. 在本构方程中, 引入了应变松弛时间与应变率的乘积项, 同时, 分别在本构方程和能量方程中引入了热松弛时间因子. 其后, 该理论被用于研究受移动热源作用的压电热弹一维问题的动态响应问题. 采用拉普拉斯变换及其数值反变换, 对问题进行了求解, 得到了不同应变松弛时间和热源移动速度下的瞬态响应, 即无量纲温度、位移、应力和电势的分布规律, 并重点考察了应变率对各物理量的影响效应, 将结果以图形形式进行了表示. 结果表明: 应变率对温度、位移、应力和电势的分布规律有显著影响.      

低密度开孔泡沫蠕变松弛特性分析

利用三维Voronoi模型和有限元方法分析了胞壁材料具有粘弹特性的低密度开孔泡沫的蠕变和应力松弛行为.采用了三参量标准线性固体模型来描述胞壁材料的粘弹特性.所得结果表明.低密度开孔泡沫具有与其胞壁材料相同的松弛时间,当相对密度较低时(低于1%)开孔泡沫的松弛模量与胞壁材料的松弛模量和泡沫相对密度平方成正比.此外,计算结果还表明,低密度开孔泡沫在较小的初始应力条件下具有与其胞壁材料相同的延迟时间.其蠕变柔度与胞壁材料的蠕变柔度和泡沫相对密度平方倒数基本成正比.但随着初始应力值的增大,泡沫的延迟时间将会显著增加.     

坡形加热下的二维广义磁热黏弹性问题研究

运用具有一个热松弛时间的广义热黏弹性理论，研究了处于均布磁场中的二维磁热黏弹 性问题. 运用Laplace变换(对时间变量)和Fourier变换(对于一个空间变量)，得到了变 换域内场量的精确表达式，并把结果应用到表面受到坡形加热的半空间问题. 应用 数值逆变换得到了时间-空间域内场量的解，对丙烯酸塑料 给出场量的响应图. 并把运用广义热黏弹性理论所得的结果与传统热黏弹性理 论及热弹性理论下的结果进行了比较.     

复合固体推进剂拉伸蠕变柔量计算的新方法

复合固体推进剂的流变特性可用Wiechert体描述,其拉伸蠕变柔量只需求解一个一元n次方程就可得到.此方程的根被证实全是实根,且两两互异.算例验证表明了文中给出的蠕变柔量具有很高的精度,但计算量小,程序实现简单.此方法还可以用于求解固体推进剂的体积蠕变柔量.     

A GENERALIZED PIEZOELECTRIC-THERMOELASTIC THEORY WITH STRAIN RATE AND ITS APPLICATION 1)

工程中大量材料的形变介于弹性与黏性之间, 既具有弹性固体特性, 又具有黏性流体特点, 即为黏弹性. 黏弹性使得材料出现很多力学松弛现象, 如应变松弛、滞后损耗等行为. 在研究受热载荷作用的多场耦合问题的瞬态响应时, 考虑此类问题中的热松弛和应变松弛现象, 对准确描述其瞬态响应尤为重要. 针对广义压电热弹问题的瞬态响应, 尽管已有学者建立了考虑热松弛的广义压电热弹模型, 但迄今, 尚未计入应变松弛. 本文中, 考虑到材料变形时的应变松弛, 通过引入应变率, 在Chandrasekharaiah广义压电热弹理论的基础之上, 经拓展, 建立了考虑应变率的广义压电热弹理论. 借助热力学定律, 给出了理论的建立过程并得到了相应的状态方程及控制方程. 在本构方程中, 引入了应变松弛时间与应变率的乘积项, 同时, 分别在本构方程和能量方程中引入了热松弛时间因子. 其后, 该理论被用于研究受移动热源作用的压电热弹一维问题的动态响应问题. 采用拉普拉斯变换及其数值反变换, 对问题进行了求解, 得到了不同应变松弛时间和热源移动速度下的瞬态响应, 即无量纲温度、位移、应力和电势的分布规律, 并重点考察了应变率对各物理量的影响效应, 将结果以图形形式进行了表示. 结果表明: 应变率对温度、位移、应力和电势的分布规律有显著影响.     

鲫鱼尾鳍的黏弹性力学性质测量

鱼类在游动过程中表现出高效率和机动灵活的特点。鱼类通过拍动鱼鳍获得动力,并对游动过程进行控制,因此鱼鳍的力学性质影响着鱼类在游动中的表现。在以往的研究中,研究者多将鱼鳍假设为弹性材料。本文通过单轴拉伸后的松弛实验测量了鲫鱼尾鳍的黏弹性力学性质。在松弛实验中,拉力随着时间的增加而逐渐减小。在实验的前100s时间内,拉力衰减至最大拉力的75%。本文采用五参数的线性黏弹性模型对松弛实验的数据进行了拟合。基于拟合得到的模型,发现在快速起动及巡游过程中,鲫鱼尾鳍的黏弹性性质能够增加鲫鱼尾鳍的表观刚度,同时在巡游过程中,由于黏性引起的能量耗散非常小。     

黏弹-Perzyna黏塑性有限元法应力更新隐式算法

黏弹-黏塑性耦合模型的黏弹性部分由弹簧、黏壶和Kelvin链串联而成,黏塑性部分为双曲线型DruckerPrager屈服函数、各向同性硬化和Perzyna黏塑性流动模型。基于黏弹性蠕变柔度,通过定义与弹性问题相对应的与时间增量相关的黏弹性剪切模量和体积模量,导出增量递推形式的本构方程。为保证算法的收敛和稳定性,把Perzyna黏塑性流动方程转化为与弹塑性相似的一致性条件,建立黏塑性增量因子单侧逼近其收敛值的N-R迭代算法。最后,给出应力更新完全隐式算法和最终计算公式。分别采用黏弹性、黏弹-塑性和黏弹-黏塑性本构关系对一地基蠕变模型进行三维有限元分析和比较,结果表明,本文算法具有较高的计算效率和稳定性。     

Investigation of relaxation properties of polymer melts by comparison of relaxation time spectra calculated with different algorithms

In this paper a recently introduced algorithm (Brabec and Schausberger, 1995) for the calculation of relaxation time spectra is compared with two standard methods, i.e., Weese's regularization, and Baumgaertel's and Winter's regression algorithm. A reasonable agreement between those three algorithms is found for the relaxation properties of mono-, polydisperse, bi-, and multimodal polystyrene samples. All three numerical methods reproduce the relaxation properties for long and medium times correctly, but they show some disagreement at short times because of sparse experimental data. The high numerical accuracy opens the possibility to test and improve the physical models which underlie the calculations. The good agreement of the different algorithms suggests that small inconsistencies to physical models are not due to a failure of the numerical methods, but due to an insufficiency of the generalized Maxwell model.     

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.     

网纹红土分数阶应力松弛模型

网纹红土松弛特性具有非线性。基于分数阶微积分理论,探讨了能更准确描述网纹红土松弛非线性特性和全过程的FVMS(Fractional Voigt and Maxwell model in series)松弛模型和FVMP(Fractional Voigt and Maxwell model in parallel)松弛模型及其理论解,进而应用提出的模型对三轴松弛试验实测数据进行反演,讨论了分数阶阶数的敏感性,并与西原模型和Burgers模型进行对比分析。研究结果表明,建立的四元件分数阶松弛本构模型应用于网纹红土应力松弛特性分析是有效可行的,模型灵活且精度更高,参数确定简便,发现分数阶阶数对应力松弛量的影响较大,但其对FVMS模型和FVMP模型松弛速率的影响不同,为实际工程长期稳定性分析提供了参考。     

Flow of a non-Newtonian fluid between intersecting planes of which one is moving

We study the flow of an Oldroyd-B fluid between two intersecting plates, one of which is fixed and the other moving along its plane. This problem was first considered by Strauss (1975) for the Maxwell fluid using a similarity transformation. We find that even in the case of a Maxwell fluid, which can be obtained by setting a specific parameter, say , in the Oldroyd-B model to zero, our results disagree with those of Strauss (1975). We find that circulating cells are present, adjacent to the stationary plate while Strauss (1975) finds them adjacent to the moving plate. We also delineate the effect of the coefficient , which is a measure of the elasticity of the flow, on the flow pattern. We find that an increase in the elastic parameter reduces the cellular structure.     

Rheological representation of fractional order viscoelastic material models

We develop rheological representations, i.e., discrete spectrum models, for the fractional derivative viscoelastic element (fractional dashpot or springpot). Our representations are generalized Maxwell models or series of Kelvin-Voigt units, which, however, maintain the number of parameters of the corresponding fractional order model. Accordingly, the number of parameters of the rheological representation is independent of the number of rheological units. We prove that the representations converge to the corresponding fractional model in the limit as the number of units tends to infinity. The representations extend to compound fractional derivative models such as the fractional Maxwell model, fractional Kelvin-Voigt model, and fractional standard linear solid. Computational experiments show that the rheological representations are accurate approximations of the fractional order models even for a small number of units.     

Analysis of thermal responses in a two-dimensional porous medium caused by pulse heat flux

In this article, the generalized model for thermoelastic waves with two relaxation times is utilized to compute the increment of temperature, the displacement components, the stress components, and the changes in the volume fraction field in a two-dimensional porous medium. By using the Fourier-Laplace transform and the eigenvalue method, the considered variables are obtained analytically. The derived approach is estimated with numerical outcomes which are applied to the porous media with a geometrical simplification. The numerical results for the considered variables are performed and presented graphically. Finally, the outcomes are represented graphically to display the difference among the classical dynamical(CD) coupled, the Lord-Shulman(LS), and the Green-Lindsay(GL) models.     

Structure-energy analysis of models of nonlinearly viscoelastic materials with several aging and viscosity functions

We consider generalized one-dimensional Maxwell and Kelvin-Voigt models of viscoelastic materials in which the properties of elastic and viscous elements are determined by the corresponding secant moduli and viscosity coefficients, which are functions of the parameters determined by the deformation process. In contrast to the nonlinear endochronic theory of aging viscoelastic materials (NETAVEM), in which one and the same aging function is used to describe the properties of all elastic elements and one and the same viscosity function is used to describe the properties of all viscous elements [1, 2], it is assumed that the type of these functions is distinct for each elementary model. For the generalized Maxwell and Kelvin-Voigt models under study, we obtain representations of the specific work of internal forces as the sum of four terms of different physical meaning. There representations are similar to those given in [1, 2] for NETAVEM. An example of construction of viscoelasticity constitutive relations containing two aging functions and one viscosity function is given for a material whose properties are sensitive to the strain rate. The simultaneous use of several aging and viscosity functions to describe the properties of structure elements of the model and the use of several components of specific work as arguments of these functions allows us to extend the scope of the models under study.     

Relaxation and retardation functions of the Maxwell model with fractional derivatives

A four-parameter Maxwell model is formulated with fractional derivatives of different orders of the stress and strain using the Riemann-Liouville definition. This model is used to determine the relaxation and retardation functions. The relaxation function was found in the time domain with the help of a power law series; a direct solution was used in the Laplace domain. The solution can be presented as a product of a power law term and the Mittag-Leffler function. The retardation function is determined via Laplace transformation and is solely a power law type.The investigation of the relaxation function shows that it is strongly monotonic. This explains why the model with fractional derivatives is consistent with thermodynamic principles.This type of rheological constitutive equation shows fluid behavior only in the case of a fractional derivative of the stress and a first order derivative of the strain. In all other cases the viscosity does not reach a stationary value.In a comparison with other relaxation functions like the exponential function or the Kohlrausch-Williams-Watts function, the investigated model has no terminal relaxation time. The time parameter of the fractional Maxwell model is determined by the intersection point of the short- and long-rime asymptotes of the relaxation function.