An adaptation of finite linear viscoelasticity theory for rubber-like viscoelasticity by use of a generalized strain measure

In this paper a simplified three-dimensional constitutive equation for viscoelastic rubber-like solids is derived by employing a generalized strain measure and an asymptotic expansion similar to that used by Coleman and Noll (1961) in their derivation of finite linear viscoelasticity (FLV) theory. The first term of the expansion represents exactly the time and strain separability relaxation behavior exhibited by certain soft polymers in the rubbery state and in the transition zone between the glassy and rubbery states. The relaxation spectra of such polymers are said to be deformation independent. Retention of higher order terms of the asymptotic expansion is recommended for treating deformation dependent spectra.Certain assumptions for the solid theory are relaxed in order to obtain a constitutive equation for uncross-linked liquid materials which exhibit large elastic recovery properties.Apart from the strain energyW(I₁,I ₂), which alternatively characterizes the long-time elastic response of solids or the instantaneous elastic response of elastic liquids, only the linear viscoelastic relaxation modulus is required for the first-order theory. Both types of material functions can be obtained, in theory, from simple laboratory testing procedures. The constitutive equations for solids proposed by Chang, Bloch and Tschoegl (1976) and a special form of K-BKZ theory for elastic liquids are shown to be particular cases of the first-order theory.Previously published experimental data on a cross-linked styrene-butadiene rubber (SBR) and an uncross-linked polyisobutylene (PIB) rubber is used to corroborate the theory.     

A thermodynamic approach to non-linear viscoelasticity

We show that Extended Irreversible Thermodynamics (EIT) may be regarded as a unifying thermodynamic scheme which generates a variety of non-linear viscoelastic equations of state. In particular, we show how the so-called eight-constants Oldroyd model, with all its invariant properties, may be derived from the postulates of EIT. It is also shown that the structure of the thermodynamic based rheological equations of state is in complete agreement with that of the representation theorems formulation of constitutive relationships originally proposed by Rivlin and Ericksen. We discuss how more complex rheological models can be derived from higher order approximations of the theory. Finally, the connections between our approach and other related works are briefly discussed.     

Finite rise time step strain modeling of nearly monodisperse polymer melts and solutions

The step shear strain experiment is one of the fundamental transient tests used to characterize the rheology of viscoelastic polymer melts and solutions. Many melts and solutions exhibit homogeneous deformation and stress relaxation; in these cases the transient dynamics can be modeled by completely ignoring momentum effects and imposing singular kinematics. Recently, however, it has been observed that there are certain classes of nearly monodisperse melts and solutions that exhibit anomalous nonhomogeneous deformation and stress relaxation (Morrison and Larson (1990), Larson, Khan, and Raju (1988), Vrentas and Graessley (1982), and Osaki and Kurata (1980)). We demonstrate that, for these classes, a finite rise time must be incorporated, some source of inhomogeneity must be present, and a small amount of added Newtonian viscosity is necessary. We examine five nonlinear and quasilinear models; the Johnson-Segalman, Phan Thien Tanner, Giesekus, White-Metzner, and Larson models. We determine which mathematical features of the models are necessary and/or sufficient to describe the observed experimental behavior.     

A differential constitutive equation for polymer melts

A modification of the Giesekus constitutive equation is derived by incorporating (approximately, via the Peterlin approximation) the finite extensibility of polymer molecules into dumbbell kinetic theory along with the anisotropic hydrodynamic drag suggested by Giesekus. The constitutive equation that is obtained retains much of the simplicity of Giesekus' constitutive equation, but it involves terms that are cubic in the stress as well as those that are quadratic. It is shown that the constitutive equation quantitatively describes the steady elongational viscosity of the IUPAC polymer melt A (including the strain softening of the melt), but it cannot describe the elongational and shear viscosities simultaneously. It is also shown that the constitutive equation satisfies the Lodge-Meissner relation for shear strains less than unity.     

Viscometric flows of polymeric melts treated according to the rational theory of Eckart's continuum. I. General theory

The purely rational theory of Eckart continua (i.e. elastic bodies with a variable relaxed state) is applied to viscometric flows of polymeric melts. The main assumptions are thermodynamic non-interaction of inelastic behaviour and of non-elastic stress, as well as elastic isotropy. After establishing the time-dependent differential equations of viscometric flow, these equations are simplified to a set of algebraic equations describing steady-state flow. From this we deduce two general equations connecting the three elastic steady-state viscometric functions which do not depend upon the elastic behaviour. The law of rubber elasticity used in this paper is described in the Appendix.     

Stress optical measurement of the third normal stress difference in polymer melts under oscillatory shear

Results are reported for the dynamic moduli,G andG, measured mechanically, and the dynamic third normal stress difference, measured optically, of a series bidisperse linear polymer melts under oscillatory shear. Nearly monodisperse hydrogenated polyisoprenes of molecular weights 53000 and 370000 were used to prepare blends with a volume fraction of long polymer, _L, of 0.10, 0.20, 0.30, 0.50, and 0.75. The results demonstrate the applicability of birefringence measurements to solve the longstanding problem of measuring the third normal stress difference in oscillatory flow. The relationship between the third normal stress difference and the shear stress observed for these entangled polymer melts is in agreement with a widely predicted constitutive relationship: the relationship between the first normal stress difference and the shear stress is that of a simple fluid, and the second normal stress difference is proportional to the first. These results demonstrate the potential use of 1,3-birefringence to measure the third normal stress difference in oscillatory flow. Further, the general constitutive equation supported by the present results may be used to determine the dynamic moduli from the measured third normal stress difference in small amplitude oscillatory shear. Directions for future research, including the use of birefringence measurements to determineN ₂/N ₁ in oscillatory shear, are described.     

An experimental study of exit pressures for polymer melts

A slit die apparatus is used to measure exit pressures for five different polymer melts. Viscosity data obtained from the same apparatus agree well with values obtained from a cone-and-plate rheometer or a capillary rheometer. Except for a PVC sample where thermal degradation was found to occur, the exit pressures obtained by linear extrapolation of the measured pressure profiles are all positive, and increase with increasing shear stress. The values of the first normal stress difference calculated according to the exit pressure theory are of the right order of magnitude and in some cases correlate satisfactorily with values measured in a cone-and-plate rheometer. However, the high sensitivity of the exit pressure values to the method of extrapolation and the wild scatter of exit pressure data for some materials make it difficult to use the exit pressure method as a routine procedure for accurate determination of the first normal stress difference.     

Orthogonal superposition of small and large amplitude oscillations upon steady shear flow of polymer fluids

The orthogonal superposition of small and large amplitude oscillations upon steady shear flow of elastic fluids has been considered. Theoretical results, obtained by numerical methods, are based on the Leonov viscoelastic constitutive equation. Steady-state components, amplitudes and phase angles of the oscillatory components of the shear stress, the first and second normal stress differences as functions of shear rate, deformation amplitude and frequency have been calculated. These oscillatory components include the first and third harmonic of the shear stresses and the second harmonic of the normal stresses. In the case of small amplitude superposition, the effect of the steady shear flow upon the frequency-dependent storage modulus and dynamic viscosity has been determined and compared with experimental data available in literature for polymeric solutions. The predicted results have been found to be in fair agreement with the experimental data at low shear rates and only in qualitative agreement at high shear rates and low frequencies. A comparison of the present theoretical results has also been made with the predictions of other theories.In the case of large amplitude superposition, the effect of oscillations upon the steady shear flow characteristics has been determined, indicating that the orthogonal superposition has less influence on the steady state shear stresses and the first difference of normal stresses than the parallel superposition. However, in the orthogonal superposition a more pronounced influence has been observed for the second difference of normal stresses.     

Elongational rheology of polyethylene melts — temporary network constitutive laws

The uniaxial elongational properties of various polyethylenes have been evaluated using an elongational rheometer and a melt-strength apparatus. It is possible to derive the data obtained in elongation from the distribution of relaxation times obtained from oscillatory shearing measurements (linear viscoelasticity), using a Wagner constitutive equation. The effects of the molecular parameters of the samples have been studied, in particular the effect of polydispersity on the shape of the damping function.     

Velocity profiles in rectangular channel flow of liquid crystalline polymer solutions

The motion of aqueous solutions of hydroxypropylcellulose in a shallow channel was visualized by means of an electrochemical technique. Isotropic and anisotropic solutions of approximately the same viscosity were used. The velocity profiles of the liquid crystalline solution were found to be qualitatively very different from those of the isotropic one. In particular, anomalous maxima in the velocity near the side walls were observed.     

Comparison between uniaxial and biaxial elongational flow behavior of viscoelastic fluids as predicted by differential constitutive equations

Behavior of polymer melts in biaxial as well as uniaxial elongational flow is studied based on the predictions of three constitutive models (Leonov, Giesekus, and Larson) with single relaxation mode. Transient elongational viscosities in both flows are calculated for three constitutive models, and steady-state elongational viscosities are obtained as functions of strain rates for the Giesekus and the Larson models.Change of elongational flow behavior with adjustable parameter is investigated in each model. Steady-state viscosities _E and _B are obtained for the Leonov model only when the strain-hardening parameter is smaller than the critical value _cr determined in each flow. In this model, uniaxial elongational viscosity _E increases with increasing strain rate , while biaxial elongational viscosity _B decreases with increasing biaxial strain rate _B . The Giesekus model predictions depend on the anisotropy parameter . _E and _B increase with strain rates for small _B while they decrease for large . When is 0.5, _E in increasing, but _B is decreasing. The Larson model predicts strain-softening behavior for both flows when the chain-contraction parameter > 0.5. On the other hand, when is small, the steady-state viscosities of this model show distinct maximum around = _B = 1.0 with relaxation time . The maximum is more prominent in _E than in _B .     

Measurement of rheological and ultrasonic properties of food and synthetic polymer melts

A slit die viscometer has been used in conjunction with a co-rotating twin screw extruder to study the rheological behaviour of maize grits, potato powder and low density polyethylene, as a function of feed rate, screw speed and temperature. The shear viscosity of both maize and potato decreased with increasing feed rate. Increasing the temperature or screw speed at any given feed rate also reduced the viscosity. The ultrasonic velocity through the material has also been shown to be sensitive to the extruder operating conditions. Overall, the ultrasonic velocity decreased as screw speed and temperature increased and feed rate decreased.     

Pressure driven flow of dilute polymer solutions in a channel: Analytic results for very small and for very large channels

Kinetic theory of dilute macromolecular solutions is applied to pressure driven flow in a small channel where wall- (and interfacial) layers have to be reckoned with. The complete rheology is studied. It turns out that for very small channels both the shear stress and the normal stress are an order of magnitude larger than corresponding quantities in simple shear. On the other hand, when the channel is so wide that the wall layers are very thin in comparison, agreement with results appropriate for simple shear is found. The volume flow rate-pressure difference relation is derived and compared to the prediction which utilizes the slip velocity concept. For very small channels, this concept is five orders of magnitude off, but reproduces asymptotically correct results for very large channels.     

Rheological characteristics of sewage sludge: A granuloviscous material

Problems have been experienced in the screw conveying of a centrifuge dewatered sewage sludge, at about 66% by weight water content. The rheology of the sludge was therefore tested in a viscometer in which the normal stress between the platens could be controlled. Rough platens were used to measure internal shear strengths while smooth steel platens were used to study wall slip.The sludge displayed a number of the characteristics of a granuloviscous material. It was found to possess high initial shear strength, but rapidly softened with shear to lower steady-state strengths. The initial shear strength increased with increasing normal stress, indicating that the sludge was frictional. The initial strength was lower against the smooth platen.The steady-state shear strengths were both rate and stress dependent using the rough platen, and this internal shear strength could be described by a frictional pseudoplastic equation. Using the smooth platen the material was nearly rate independent, and was best described by a Coulomb friction equation. At the higher rates the wall shear strengths were lower than the internal shear strengths.The strengths (both initial and steady-state) were characterised by a high degree of variability from sample to sample, and also showed marked hysteresis with the rough platen. They did not show any consistent variation with moisture content.In terms of screw conveying, the major conclusion is that smooth walls to the flights of the screw and barrel of the conveyor are desirable.     

A “zero-parameter” constitutive relation for simple shear viscoelasticity

Based on the Cox-Merz rule and Eyring's expression for the nonlinear shear viscosity, a Wagner-type constitutive relation with no nontrivial adjustable parameters is proposed for simple shear viscoelasticity. The predictions for a number of non-steady shear flows are worked out analytically. It is shown that most features of shear viscoelasticity are reproduced by the model.     

Comparison of linear viscoelastic constitutive equations for non-uniform polymer melts

Linear-viscoelastic properties of polydisperse and randomly-branched polymer melts were fit with several proposed relaxation functions by non-linear regression. Three polymer systems were investigated, including 1) crosslinked polyethylenes, 2) polydisperse linear poly(dimethylsiloxane)s, and 3) Marlex polyethylenes, which are polydisperse and probably contain long-chain branching. Four relaxation functions were evaluated, including the Rouse, reptation, stretched-exponential, and stretched-exponential-power-law (SEPL) relaxation functions. The SEPL best described each series of polymers, and therefore may be a general relaxation function for non-uniform polymer melts. The flow activation energy for crosslinked polyethylene may be coupled to a breadth-of-relaxation index, indicating that a coupling between a characteristic short relaxation time and longest relaxation time, as suggested by Ngai and Plazek (J. Polym. Sci. Polym. Phys. Ed. 1985, 23:2159–2180), may hold for some non-uniform polymers.     

The viscoelastic constitutive differential operator law in terms of the relaxation and retardation spectra

The case of a linear viscoelastic medium is considered. The Carson transforms of the relaxation and retardation functions are expressed in two different ways, taking on the one hand the differential operator form of the constitutive equation, and on the other hand the generalized mechanical models. By identification we deduce general explicit expressions for the constant coefficients of the differential operator law, in terms of the discrete relaxation and retardation spectra.     

Shear viscosity of polybutyleneterephthalate/liquid crystal polymer blends

Viscosity measurements have been carried out on blends of polybuty-leneterephthalate (PBT) and a liquid crystalline copolyesteramide (LCP). The flow curves of the blends with LCP content larger than 20%, show a behavior similar to that of the pure LCP, with a rapid rise of the viscosity at low shear rates. The viscosity-composition curves exhibit a deep minimum at low LCP content which may be mainly attributed to the lack of interactions between the two phases.     

Driven torsion pendulum for measuring the complex shear modulus in a steady shear flow

To investigate the viscoelastic behavior of fluid dispersions under steady shear flow conditions, an apparatus for parallel superimposed oscillations has been constructed which consists of a rotating cup containing the liquid under investigation in which a torsional pendulum is immersed. By measuring the resonance frequency and bandwidth of the resonator in both liquid and in air, the frequency and steady-shear-rate-dependent complex shear modulus can be obtained. By exchange of the resonator lumps it is possible to use the instrument at four different frequencies: 85, 284, 740, and 2440 Hz while the steady shear rate can be varied from 1 to 55 s^–1. After treatment of the theoretical background, design, and measuring procedure, the calibration with a number of Newtonian liquids is described and the accuracy of the instrument is discussed.Notation a radius of the lump - A geometrical constant - b inner radius of the sample holder - c constant - C ₁, C ₂ apparatus constants - D damping of the pendulum - e _x , e _y , e _z Cartesian basis - e _r , e , e _z orthonormal cylindrical basis - E geometrical constant - E _t , ⁰ E _t , _t relative strain tensor - f function of shear rate - F _t relative deformation tensor - G (t) memory function - G ^* complex shear modulus - G Re(G ^* ) - G Im(G ^* ) - h distance between plates - H ^* transfer function - , functional - i imaginary unit: i ²= – 1 - I moment of inertia - J _exc excitation current - J ₀ amplitude of J _exc - k ^* = k – ik complex wave number - K torsional constant - K fourth order tensor - l length of the lump - L mutual inductance - M _dr driving torque - M _liq torque exerted by the liquid - ⁰ M _liq, _liq steady state and dynamic part of Mliq - n power of the shear rate - p isotropic pressure - Q quality factor - r radial position - R,R ₀, R _c Re(Z ^*, Z ₀ ^* , Z _c ^* ) - s time - t, t time - T temperature - T, ⁰ T, stress tensor - u velocity - U lock-in output - ₀ velocity - V _det detector output voltage - V _sig, V _cr signal and cross-talk part of V _det - x Cartesian coordinate - X , X ₀, X _c Im(Z ^*, Z ₀ ^* , Z _c ^* ) - y Cartesian coordinate - z Cartesian coordinate, axial position     

The role of short chain molecules for the rheology of polystyrene melts.

Atactic polystyrenes of narrow molar mass distribution with average molar masses larger than the critical molar massM _c were mixed with similar polystyrenes of molecular masses lower thanM _c . Linear viscoelastic melt properties of these binary blends were measured with a dynamic viscometer of the concentric cylinder type. One of the experimental findings is that the time-temperature shift factorsa _T are dependent on the composition of the samples. This can be understood, if free volume due to chain-ends is taken into account. A computer-fitted WLF-equation being modified in a proper way leads to the following results: At the glass-transition-temperature the fraction of free volume in polystyrene of infinite molar mass is only 0.015. At a temperature of 180 °C the mean value of the free volume at a chain end is 0.029 nm³ for the polystyrene investigated.