麦秆发泡包装衬垫非线性粘弹性模型及参数识别

由非线性粘弹性材料的一维Volterra-Frecher积分型本构方程推导得出了非线性粘弹性材料的一维微分型本构方程,建立了麦秆发泡包装衬垫的非线性粘弹性模型。根据在一种跌落高度、一种衬垫厚度和静态应力条件下的麦秆发泡包装衬垫的动态压缩试验数据,识别出了其缓冲性能模型参数。     

Higher-order fractional constitutive equations of viscoelastic materials involving three different parameters and their relaxation and creep functions

Two higher-order fractional viscoelastic material models consisting of the fractional Voigt model (FVM) and the fractional Maxwell model (FMM) are considered. Their higher-order fractional constitutive equations are derived due to the models’ constructions. We call them the higher-order fractional constitutive equations because they contain three different fractional parameters and the maximum order of equations is more than one. The relaxation and creep functions of the higher-order fractional constitutive equations are obtained by Laplace transform method. As particular cases, the analytical solutions of standard (integer-order) quadratic constitutive equations are contained. The generalized Mittag–Leffler function and H-Fox function play an important role in the solutions of the higher-order fractional constitutive equations. Finally, experimental data of human cranial bone are used to fit with the models given by this paper. The fitting plots show that the models given in the paper are efficient in describing the property of viscoelastic materials.     

A systematic methodology for creep master curve construction using the stepped isostress method (SSM): a numerical assessment

Long-term creep of viscoelastic materials is experimentally inferred through accelerating techniques based on the time–temperature superposition principle (TTSP) or on the time–stress superposition principle (TSSP). According to these principles, a given property measured for short times at a higher temperature or higher stress level remains the same as that obtained for longer times at a lower temperature or lower stress level, except that the curves are shifted parallel to the horizontal axis, matching a master curve. These procedures enable the construction of creep master curves with short-term experimental tests.The Stepped Isostress Method (SSM) is an evolution of the classical TSSP method. Higher reduction of the required number of test specimens to obtain the master curve is achieved by the SSM technique, since only one specimen is necessary. The classical approach, using creep tests, demands at least one specimen per each stress level to produce a set of creep curves upon which TSSP is applied to obtain the master curve.This work proposes an analytical method to process the SSM raw data. The method is validated using numerical simulations to reproduce the SSM tests based on two different viscoelastic models. One model represents the viscoelastic behavior of a graphite/epoxy laminate and the other represents an adhesive based on epoxy resin.     

The effect of annealing on the elastoplastic and viscoelastic responses of isotactic polypropylene

Observations are reported on isotactic polypropylene (i) in a series of tensile tests with a constant strain rate on specimens annealed for 24 h at various temperatures in the range from 110 to 150 °C, (ii) in two series of creep tests in the subyield region of deformations on samples not subjected to thermal treatment and on specimens annealed at 140 °C, and (iii) in a series of tensile relaxation tests on non-annealed specimens. Constitutive equations are derived for the elastoplastic and non-linear viscoelastic responses of semicrystalline polymers. A polymer is treated as an equivalent transient network of macro-molecules bridged by junctions (physical cross-links, entanglements and lamellar blocks). The network is assumed to be highly heterogeneous, and it is thought of as an ensemble of meso-regions with different activation energies for separation of strands from temporary nodes. The elastoplastic behavior is modelled as sliding of junctions in meso-domains with respect to their reference positions driven by macro-deformation. The viscoelastic response is attributed to detachment of active strands from temporary junctions and attachment of dangling chains to the network. Constitutive equations for isothermal deformations with small strains are derived by using the laws of thermodynamics. Adjustable parameters in the stress–strain relations are found by fitting the experimental data.     

粘弹性对称铺设层合板的非线性动力响应

基于Karman板理论和线粘弹性Boltzmann叠加原理,建立了粘弹性对称铺设层合板的非线性积分—偏微分动力学方程。针对材料具积分型本构关系以及松弛模量为Prony级数的形式,应用Galerkin技术、Newmark方法和Newton-Cotes方法,给出了求解粘弹性层合结构非线性动力学问题的一种有效的数值算法。具体地求解了若干算例,且与相关文献进行了比较。     

A relationship between non-exponential stress relaxation and delayed elasticity in the viscoelastic process in amorphous solids: Illustration on a chalcogenide glass

Inorganic glasses are viscoelastic materials since they exhibit, below as well as above their glass transition temperature, a viscoelastic deformation under stress, which can be decomposed into a sum of an elastic part, an inelastic (or viscous) part and a delayed elastic part. The delayed elastic part is responsible for the non-linear primary creep stage observed during creep tests. During a stress relaxation test, the strain, imposed, is initially fully elastic, but is transformed, as the stress relaxes, into an inelastic and a delayed elastic strains. For linear viscoelastic materials, if the stress relaxation function can be fitted by a stretched exponential function, the evolution of each part of the strain can be predicted using the Boltzmann superposition principle. We develop here the equations of these evolutions, and we illustrate their accuracy by comparing them with experimental evolutions measured on GeSe₉ glass fibers. We illustrate also, by simple equations, the relationship between any kind of relaxation function based on additive contribution of different relaxation processes and the delayed elastic contribution to stress relaxation: the delayed elasticity is directly correlated to the dispersion of relaxations times of the processes involved during relaxation.     

橡胶材料的非线性黏弹性本构方程

橡胶是一种非线性黏弹性材料,准确描述其非线性黏弹性力学响应的本构方程是橡胶材料及制品设计优化的关键。文中基于超弹性模型和并行流变模型(PRF)描述了橡胶材料的非线性黏弹性响应特征,重点探讨了由实验数据确定PRF本构方程材料参数的方法。首先通过单轴拉伸实验数据拟合得到超弹性模型,将应力松弛实验数据拟合得到表征材料线性黏弹性的模型-Prony级数,再将Prony级数转化为初始的PRF模型,进而不断优化得到PRF模型的准确材料参数,最后进行实验验证。结果表明,PRF模型计算的不同应变下的应力松弛数据与实验数据之间的误差仅为0.067%,PRF能准确地描述橡胶材料非线性响应的应力松弛行为。     

Application of Time-Based Fractional Calculus Methods to Viscoelastic Creep and Stress Relaxation of Materials

The quasi-static viscoelastic response of polymeric materialsis investigated utilizing constitutive models based on fractionalcalculus. Time-based fractional calculus analysis techniques areemphasized. Analytic solutions to quasi-static boundary value problemsin which the viscoelastic behavior is characterized by thefour-parameter fractional calculus-based solid model are given. Varioussets of data from the literature are fit with existing and newfractional calculus-based constitutive equations.     

A procedure for the computational investigation of stress-relaxation phenomena

The analysis of stress-relaxation tests in viscoelastic materials is discussed. Stress-relaxation tests are performed by means of a two-stage process, as a stretching stage up to an assumed strain level, according to a given strain rate, followed by a constant strain stage. Major attention is devoted to the influence of experimental variables, namely the strain level and the strain rate during stretching stage, on the subsequent relaxation phenomena. A constitutive analysis to interpret the behaviour shown by experimental tests is performed, accounting for the influence of the above mentioned experimental variables, to get an actual definition of time dependent behaviour of the material. This procedure is applied to investigate stress-relaxation phenomena in polymeric materials and soft biological tissues, as ligaments. The results confirm the necessity to analyse data from stress-relaxation tests taking into account the experimental variables that characterize the different stages, and evaluating the influence of the specific variables outlined.     

Numerical modeling of micro- and macro-behavior of viscoelastic porous materials

This paper presents a numerical method for solving the two-dimensional problem of a polygonal linear viscoelastic domain containing an arbitrary number of non-overlapping circular holes of arbitrary sizes. The solution of the problem is based on the use of the correspondence principle. The governing equation for the problem in the Laplace domain is a complex hypersingular boundary integral equation written in terms of the unknown transformed displacements on the boundaries of the holes and the exterior boundaries of the finite body. No specific physical model is involved in the governing equation, which means that the method is capable of handling a variety of viscoelastic models. A truncated complex Fourier series with coefficients dependent on the transform parameter is used to approximate the unknown transformed displacements on the boundaries of the holes. A truncated complex series of Chebyshev polynomials with coefficients dependent on the transform parameter is used to approximate the unknown transformed displacements on the straight boundaries of the finite body. A system of linear algebraic equations is formed using the overspecification method. The viscoelastic stresses and displacements are calculated through the viscoelastic analogs of the Kolosov–Muskhelishvili potentials, and an analytical inverse Laplace transform is used to provide the time domain solution. Using the concept of representative volume, the effective viscoelastic properties of an equivalent homogeneous material are then found directly from the corresponding constitutive equations for the average field values. Several examples are given to demonstrate the accuracy of the method. The results for the stresses and displacements are compared with the numerical solutions obtained by commercial finite element software (ANSYS). The results for the effective properties are compared with those obtained with the self-consistent and Mori–Tanaka schemes.     

A dynamic variational multiscale method for viscoelasticity using linear tetrahedral elements

In this article, we develop a dynamic version of the variational multiscale (D‐VMS) stabilization for nearly/fully incompressible solid dynamics simulations of viscoelastic materials. The constitutive models considered here are based on Prony series expansions, which are rather common in the practice of finite element simulations, especially in industrial/commercial applications. Our method is based on a mixed formulation, in which the momentum equation is complemented by a pressure equation in rate form. The unknown pressure, displacement, and velocity are approximated with piecewise linear, continuous finite element functions. To prevent spurious oscillations, the pressure equation is augmented with a stabilization operator specifically designed for viscoelastic problems, in that it depends on the viscoelastic dissipation. We demonstrate the robustness, stability, and accuracy properties of the proposed method with extensive numerical tests in the case of linear and finite deformations.     

Modeling of complex modulus of bituminous mixtures measured in tension/compression to estimate secant modulus in indirect tensile test

For mechanistic based pavement design, it is necessary to know the stiffness of materials used for road construction. In the European standard EN 12697-26, several experimental tests are proposed to measure the modulus of bituminous mixtures. However, bituminous materials exhibit strong viscoelastic behavior. Hence, the stiffness of this specific material depends not only on the sample geometry but also on the loading law (strain or stress curve versus time). As a consequence, comparison between these different tests, performed in time or frequency domain, is not straight forward. The present paper focuses on measurement of secant modulus using the Indirect tensile (IT) test and comparison with complex modulus. In the IT test, the loading law is not systematically controlled. So, it is important to investigate the influence of the loading waveform on the test result. For this study, the case of a High Modulus Bituminous Mixture for base course (French EME) has been considered. Its viscoelastic behavior has been firstly determined using complex modulus measurements and using Prony series model. An intermediate step, based on an original master curve construction method is used. Doing so, the viscoelastic model is validated on the time domain by simulating direct tensile tests. The case of IT test is then considered. The Prony series model is conformed against IT test results for which the force waveform is known. For this last test, the dependance of the Poisson ratio with temperature and loading time is highlighted. Assuming a purely elastic behavior in isotropic compression, a formula, able to derive the viscoelastic Poisson ratio from the complex modulus is presented. Finally, a theoretical parametric study considering the IT loading waveform is undertaken. It appears that the correction factor given in the standard EN 12697-26 cannot be applied at all temperatures. Actually, this correction factor should be material dependent. This study has been carried out in the framework of collaboration between LCPC and USIRF (Union des Syndicats de l’Industrie Routière française).     

泡沫塑料包装衬垫缓冲性能建模

本文从非线性粘弹性物质的多重积分型本构方程出发,引入塑性应变,推导了粘弹塑性物质的微分型本构方程.进一步考虑材料的损伤特性,建立了泡沫塑料包装衬垫的缓冲性能模型.根据在一种跌落高度、衬垫厚度和静态应力组合条件下的国产可发聚苯乙烯沫泡塑料的连续多次跌落冲击试验数据,识别了其缓冲性能模型参数.动力学计算与试验结果的比较表明,该模型能较好地反映可发聚苯乙烯泡沫塑料包装衬垫在连续多次跌落冲击过程中的加速度、速度、位移变化,冲击能量的吸收性和塑性变形等基本特性.并且仅根据在一种组合条件下的试验数据所识别的这一模型,在其它组合条件下其缓冲性能与试验结果也吻合较好.     

Modeling short- and long-term time-dependent nonlinear behavior of polyethylene

A new methodology for developing macromechanical constitutive formulations for time-dependent materials is presented in this article. In particular, two phenomenological constitutive models for polymer materials are illustrated, describing time-dependent and nonlinear mechanical behavior. In this new approach, short-term creep test data are used for modeling both short-term and long-term responses. The differential form of a model is used to simulate typical nonlinear viscoelastic polymeric behavior using a combination of springs and dashpots. Unified plasticity theory is then used to develop the second model, which is a nonlinear viscoplastic one. Least squares fitting is applied for the determination of material parameters for both models, based on experimental results. Due to practical constraints, experimental data are usually available for short-term time-frames. In the presented proposed formulation, the material parameters determined from short-term testing are used to obtain material parameter relationships for predicting the long-term material response. This is done by extending short-term information for longer time frames. Finally, theoretical and experimental results of tensile tests on polyethylene subjected to various load levels and test times are compared and discussed. Very good agreement of the modeling results with experimental data shows that the developed formulation provides a flexible and reliable framework for predicting load responses of polymers.     

Dynamic compression testing of particle-reinforced polymer roll cover materials

Two polymer composite roll cover materials were studied using dynamic compression tests. The main focus was determination of viscoelastic properties of the materials and development of mathematical models to describe the behavior of these materials in the contact area between two rolls. The compression tests were conducted using short rise time pulses with different durations in servohydraulic materials testing machines. In the modeling, combinations of standard elastic and viscoelastic elements were used together with the Boltzmann superposition principle. A simple spring-dashpot model was found to fit sufficiently to the experimental data with relaxation (retardation) times ranging from a few milliseconds in transient loading tests to tens of hours in static compression tests     

Validity of constitutive equations used for calculation of stresses in cooling castings

Abstract

The fundamentals of the constitutive equations used for modelling the mechanical behaviour of metallic materials over wide temperature ranges are briefly reviewed. Examples of stress-calculation problems in cooling castings are given. Four constitutive relationships used for practical calculation of the stresses are analysed. An experimental assessment of their validity, based on a comparison of stress values predicted by the models with those obtained for ring-shaped, aluminium alloy castings with metallic cores, as well as with some published data concerning beam-shaped, low-carbon steel castings, is presented. It is concluded that the essential simplifications assumed in deriving the models lead to significant disagreements with experimental results in most of the investigated cases.

MST/58     

A variational framework for fiber‐reinforced viscoelastic soft tissues

The mechanical properties of soft biological tissues vary depending on how the internal structure is organized. Classical examples of tissues are ligaments, tendons, skin, arteries, and annulus fibrous. The main element of such tissues is the fibers which are responsible for the tissue resistance and the main mechanical characteristic is their viscoelastic anisotropic behavior. The objective of this paper is to extend an existing model for isotropic viscoelastic materials in order to include anisotropy provided by fiber reinforcement. The incorporation of the fiber allows the mechanical behavior of these tissues to be simulated. The model is based on a variational framework in which its mechanical behavior is described by a free energy incremental potential whose local minimization provides the constraints for the internal variable updates for each load increment. The main advantage of this variational approach is the ability to represent different material models depending on the choice of suitable potential functions. Finally, the model is implemented in a finite‐element code in order to perform numerical tests to show the ability of the proposed model to represent fiber‐reinforced materials. The material parameters used in the tests were obtained through parameter identification using experimental data available in the literature. Copyright © 2011 John Wiley & Sons, Ltd.     

Viscoelasticity and viscoplasticity of semicrystalline polymers: Structure–property relations for high-density polyethylene

Observations are reported on two commercial grades of high-density polyethylene (HDPE) in uniaxial tensile tests, relaxation tests, creep tests and cyclic tests with a strain-controlled deformation program. Constitutive equations are derived for the viscoelastic and viscoplastic responses of semicrystalline polymers at three-dimensional deformation with small strains. A polymer is modeled as a two-phase continuum consisting of a crystalline skeleton and an amorphous phase treated as a transient network of chains. Its viscoelastic response is associated with thermally activated rearrangement of strands in the temporary network. The viscoplastic behavior reflects fine and coarse slip of lamellar stacks and sliding of junctions between chains in the network. Adjustable parameters in the stress–strain relations are found by fitting the experimental data. The study focuses on the effect of molecular weight of HDPE on its mechanical properties.     

Viscoelastic Characterisation of the Temporomandibular Joint Disc in Bovines

Abstract: The aim of this paper was to develop a biomechanical model to reproduce the viscoelastic behaviour of the material of the bovine temporomandibular joint (TMJ) disc, as a substitute material to be applied in subsequent experimental studies of possible human TMJ disorders. As the mechanical properties of the disc change due to different factors, this study is focussed on evaluating its viscoelastic response under compression in stress–relaxation, creep and cycling loading tests, using an equipment that allows us to simulate the real function conditions of the material in a liquid medium at 37 °C. In order to fit the data obtained from the relaxation and creep tests, generalised Maxwell and Kelvin models of eight terms are proposed, leading to Prony’s series of practical interest for subsequent numerical calculations.