Modeling of curing processes based on a multi-field potential. Single- and multiscale aspects

This paper provides a continuum mechanical model for the curing of polymers, including the incompressibility effects arising at the late stages of the process. For this purpose, the free energy density functional is split into a deviatoric and a volumetric part, and a multifield formulation is inserted. An integral formulation of the functional is used to depict the time-dependent material behavior. The model is also coupled with the multiscale finite element method, a numerical approach serving for the modeling of heterogeneous materials with a highly oscillatory microstructure. The effects of the proposed extensions are illustrated on the basis of several numerical examples concerned with the study of the influence of Poisson’s ratio on the curing process and the behavior of the microheterogeneous polymers.     

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

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

A new constitutive equation for solid propellant with the effects of aging and viscoelastic Poisson’s ratio

A new constitutive equation for solid propellant with the effects of aging and viscoelastic Poisson’s ratio is proposed. Effects of thermo-oxidative aging and viscoelastic Poisson’s ratio are considered in this comprehensive constitutive equation with two sets of reduced time system coping with the time and temperature dependence. In order to simulate the single and combined effects of aging and viscoelastic Poisson’s ratio, constitutive equation is rewritten into an incremental form and implemented in the user subroutine UMAT at the platform of finite element code ABAQUS. Detailed procedure for acquiring the parameters in constitutive equation is introduced and conducted for the subsequent applied analysis. Two typical loading cases during the service life of solid rocket motor and four sets of combined constitutive models are simulated. Von Mises strain and stress distribution and their changes versus time are utilized as the main analysis index. The results show that the effects of aging and viscoelastic Poisson’s ratio or their combinations will improve or decrease the level and change the distribution of Von Mises strain and stress in varying degrees.     

Dynamic stiffness of chemically and physically ageing rubber vibration isolators in the audible frequency range

The constitutive equations of chemically and physically ageing rubber in the audible frequency range are modelled as a function of ageing temperature, ageing time, actual temperature, time and frequency. The constitutive equations are derived by assuming nearly incompressible material with elastic spherical response and viscoelastic deviatoric response, using Mittag-Leffler relaxation function of fractional derivative type, the main advantage being the minimum material parameters needed to successfully fit experimental data over a broad frequency range. The material is furthermore assumed essentially entropic and thermo-mechanically simple while using a modified William–Landel–Ferry shift function to take into account temperature dependence and physical ageing, with fractional free volume evolution modelled by a nonlinear, fractional differential equation with relaxation time identical to that of the stress response and related to the fractional free volume by Doolittle equation. Physical ageing is a reversible ageing process, including trapping and freeing of polymer chain ends, polymer chain reorganizations and free volume changes. In contrast, chemical ageing is an irreversible process, mainly attributed to oxygen reaction with polymer network either damaging the network by scission or reformation of new polymer links. The chemical ageing is modelled by inner variables that are determined by inner fractional evolution equations. Finally, the model parameters are fitted to measurements results of natural rubber over a broad audible frequency range, and various parameter studies are performed including comparison with results obtained by ordinary, non-fractional ageing evolution differential equations.     

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.     

A Numerical Study of Oscillating Peristaltic Flow of Generalized Maxwell Viscoelastic Fluids Through a Porous Medium

This article presents a numerical study on oscillating peristaltic flow of generalized Maxwell fluids through a porous medium. A sinusoidal model is employed for the oscillating flow regime. A modified Darcy-Brinkman model is utilized to simulate the flow of a generalized Maxwell fluid in a homogenous, isotropic porous medium. The governing equations are simplified by assuming long wavelength and low Reynolds number approximations. The numerical and approximate analytical solutions of the problem are obtained by a semi-numerical technique, namely the homotopy perturbation method. The influence of the dominating physical parameters such as fractional Maxwell parameter, relaxation time, amplitude ratio, and permeability parameter on the flow characteristics are depicted graphically. The size of the trapped bolus is slightly enhanced by increasing the magnitude of permeability parameter whereas it is decreased with increasing amplitude ratio. Furthermore, it is shown that in the entire pumping region and the free pumping region, both volumetric flow rate and pressure decrease with increasing relaxation time, whereas in the co-pumping region, the volumetric flow rate is elevated with rising magnitude of relaxation time.     

非线性粘弹性平面问题的数学模型及其失记效应

本文以位移为基本未知量,利用非线性粘弹性力学中的Leaderman本构关系和线性几何假设,建立了非线性粘弹性平面问题的数学模型;在粘弹性泊松比为常数的情况下,利用Titchmarsh定理和Laplace变换法证明了非线性粘弹性平面问题与非线性弹性平面问题之间存在着某些对应关系,对应关系为粘弹性问题的求解提供了一种新的思路,利用这些关系可直接从相应弹性问题获得粘弹性平面问题的部分响应,与传统的时域有限差分法相比,计算时间明显缩短,另外,对应关系也揭示了粘弹性结构的失记效应,即结构的部分响应仅与外部输入的现时值有关,而与其历史无关。     

Physical Principles of Methods for Measuring Viscoelastic Properties

The main methods used to measure viscoelastic properties of materials in a wide range of frequencies from 10^?4 to 10⁶ Hz are reviewed. It is demonstrated that the accuracy of many experimental methods can be increased by taking into account the form factors, which depend on the specimen type. An example of the form factor for a cylindrical specimen is provided, which is determined numerically on the basis of a two-dimensional deformation model taking into account the specimen geometry and Poisson’s ratio. The importance of the precise determination of Poisson’s ratio for rubber-like and complex-structured materials is demonstrated. Requirements to such measurements and a setup satisfying these requirements are described. Two methods for measuring viscoelastic properties of living tissues (compliance and disturbance propagation velocity) are considered. Based on the developed method of measuring these parameters for materials with a fixed thickness, methods for the creation of a unified standard of measurements of viscoelastic characteristics of living tissues are proposed.     

Fractional Derivative Viscoelasticity at Large Deformations

A time domain viscoelastic model for large three-dimensional responses underisothermal conditions is presented. Internal variables with fractional orderevolution equations are used to model the time dependent part of the response. By using fractional order rate laws, the characteristics of the timedependency of many polymeric materials can be described using relatively fewparameters. Moreover, here we take into account that polymeric materials are often used in applications where the small deformations approximation does nothold (e.g., suspensions, vibration isolators and rubber bushings). A numerical algorithm for the constitutive response is developed and implemented into a finite element code forstructural dynamics. The algorithm calculates the fractional derivatives by means of the Grünwald–Lubich approach.Analytical and numerical calculations of the constitutive response in the nonlinearregime are presented and compared. The dynamicstructural response of a viscoelastic bar as well as the quasi-static response of athick walled tube are computed, including both geometrically and materiallynonlinear effects. Moreover, it isshown that by applying relatively small load magnitudes, the responses ofthe linear viscoelastic model are recovered.     

Discussion of the published manuscript lazopoulos K.A., karaoulanis D., lazopoulos Α.Κ. (2016) «on fractional modelling of viscoelastic mechanical systems», mechanics research communications, 78, pp. 1–5

The paper [20] attempts to solve a one dimensional model of the Zener type viscoelastic model (see Fig. 1 of the paper) replacing the classical time derivatives with a new fractional time derivative that was presented earlier by two of the authors (see ref. [1], [9] and [10] cited in the work). In addition, the authors present the classical solution of what is called in the paper “conventional viscoelastic model” and the solution due to Caputo fractional derivative, which is the solution of what the authors call “existing fractional order derivative viscoelastic models”. The paper seems to have several problems, typos and mistakes that will be discussed as follows     

Waterhammer modelling of viscoelastic pipes with a time-dependent Poisson's ratio

Poisson's ratio in viscoelastic materials is a function of time. However, recently developed waterhammer models of viscoelastic pipes consider it constant. This simplifying assumption avoids cumbersome calculations of double convolution integrals which appear if Poisson's ratio is time-dependent. The present research develops a mathematical model taking the time dependency of Poisson's ratio into account for linear viscoelastic pipes. Poisson's ratio is written in terms of relaxation function and bulk modulus which is assumed to be constant. The relaxation function is obtained from creep function given as the viscoelastic property data of pipe material. The results obtained from the present waterhammer model are compared with the experimental data for two different flow rates. The comparison reveals that with the application of the time-dependent Poisson's ratio and unsteady friction, the viscoelastic data of mechanical tests can directly be used for waterhammer analysis with less need for the calibration of the flow configuration. It was also shown that the creep curve calibrated based on the present model is closer to the actual creep curve than that calibrated based on previous models.     

Stability analysis of pipes conveying fluid with fractional viscoelastic model

Divergence and flutter instabilities of pipes conveying fluid with fractional viscoelastic model has been investigated in the present work. Attention is concentrated on the boundaries of the stability. Based on the Euler–Bernoulli beam theory for structural dynamics, viscoelastic fractional model for damping and, plug flow model for fluid flow, equation of motion has been derived. The effects of gravity, and distributed follower forces are also considered. By transferring the equation of motion to the Laplace domain and using the Galerkin method, the characteristic equations are obtained. By solving the eigenvalue problem, frequencies and dampings of the system have been obtained versus flow velocity. Some numerical test cases have been studied with viscoelastic fractional model and the effect of the fractional derivative order and the retardation time is investigated for various boundary conditions.

     

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

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

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

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

Electroosmosis modulated peristaltic biorheological flow through an asymmetric microchannel: mathematical model

A theoretical study is presented of peristaltic hydrodynamics of an aqueous electrolytic non-Newtonian Jeffrey bio-rheological fluid through an asymmetric microchannel under an applied axial electric field. An analytical approach is adopted to obtain the closed form solution for velocity, volumetric flow, pressure difference and stream function. The analysis is also restricted under the low Reynolds number assumption (Stokes flow) and lubrication theory approximations (large wavelength). Small ionic Peclét number and Debye–Hückel linearization (i.e. wall zeta potential ≤ 25 mV) are also considered to simplify the Nernst–Planck and Poisson–Boltzmann equations. Streamline plots are also presented for the different electro-osmotic parameter, varying magnitudes of the electric field (both aiding and opposing cases) and for different values of the ratio of relaxation to retardation time parameter. Comparisons are also included between the Newtonian and general non-Newtonian Jeffrey fluid cases. The results presented here may be of fundamental interest towards designing lab-on-a-chip devices for flow mixing, cell manipulation, micro-scale pumps etc. Trapping is shown to be more sensitive to an electric field (aiding, opposing and neutral) rather than the electro-osmotic parameter and viscoelastic relaxation to retardation ratio parameter. The results may also help towards the design of organ-on-a-chip like devices for better drug design.     

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.     

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