Accurate evaluation of the transverse shear modulus in unidirectional fibrous composites

A new relation for the prediction of the transverse shear modulus in unidirectional fiber composites has been derived. The theoretical results of this relationship are in better agreement with the experiments than those of other relations, existing in the literature. The discrepancies, which are observed among the theoretical predictions and the experimental values, are explained by the consideration of the boundary layers existing between the matrix and the fibers of the composite. A new model, which includes the intermediate phase between the matrix and the fiber, called the mesophase, is considered in order to take into account the above-mentioned layers.     

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 elastic modulus of particulate composites using the concept of a mesophase

A theoretical model for the evaluation of the elastic modulus in particulate composites has been developed. The method takes into account the existence of a mesophase between main phases, which constitutes an important parameter influencing the behaviour of a composite material. This layer between the matrix and filler develops different physico-chemical properties from those of the constituent phases and variable ones along its thickness. The effect of the progressive variation of the elastic modulus of the mesophase on the modulus of the composite was estimated by applying various simple laws of variation. Convenient laws of variation were introduced, varying from a simple one, assuming a linear law, to a more refined one using a parabolic law. Experimental results with particulates, based on iron-filled epoxy composites, compared satisfactorily with other models. However, the model based on a parabolic law was superior to all others on physical grounds.List of symbols E elastic modulus - v Poisson's ratio - volume fraction - r radius - V volume - K bulk modulus - r thickness of mesophase - parameter which depends on the sudden changes of the heat capacity for the filled and unfilled polymer C _p ^f and C _p ⁰      

Flow birefringence and rheological measurements on shear induced micellar structures

Aqueous solutions of cationic surfactant systems with strongly binding counterions show the striking phenomenon of shear induced phase transitions. At low shear rates or angular frequencies, the solutions exhibit Newtonian flow. At high rates of shear, however, the rheological properties change dramatically. Above a well defined threshold value of the velocity gradient, a supermolecular structure can be formed from micellar aggregates. This shear induced structure (SIS) behaves like a gel and exhibits strong flow birefringence. The formation of the shear induced structure is very complicated and depends on the specific conditions of the surfactant system. In this paper we discuss new results which have been obtained from rheological measurements and from flow birefringence data. We examine the stability of the shear induced state as a function of temperature, surfactant concentration and salt concentration and we analyse the effect of solubilisation of alcohols and hydrocarbons. The results are interpreted in terms of a kinetic model which accounts for the observed behavior.Dedicated to the 60. birthday of Prof. H. Harnisch, Hoechst AGPartly presented at the 2nd Conference of European Rheologists, Prague, June 17–20, 1986     

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.     

The rotation of a suspended axisymmetric ellipsoid in a magnetic field

Under the influence of a uniform and parallel magnetic field, a ferromagnetic fiber suspended in a Newtonian fluid rotates to align with the field direction. This study examines the field-induced rotation process for an individual non-Brownian axisymmetric ellipsoid suspended in a stagnant Newtonian fluid. Theoretical predictions are derived by a perturbation analysis for the limiting case where the strength of the applied magnetic field far exceeds the saturation magnetization of the ellipsoid. Numerical calculations are performed for the more general problem of an ellipsoid with known isotropic, non-hysteretic magnetic properties, using nickel and a stainless steel as examples. The analysis encompasses materials with field-induced, nonlinear magnetic properties, distinguishing these results from the simpler cases where the particle magnetization is either independent of, or linearly dependent on, the strength of the applied external field. In this study, predictions indicate that when the ellipsoid is magnetically saturated, the particle rotation is governed by the magnetoviscous time constant, _MV = _s/₀ M _s ² . It is found that the rotation rate depends strongly on the aspect ratio,a/b, of the ellipsoid, but only weakly on the dimensionless magnetization,M _s/H ₀. A geometric parameter for an ellipsoid, defined in eq. (2.5) - a, b major, minor semi-axes of an axisymmetric ellipsoid - D demagnetization tensor for an ellipsoid - D ^M magnetometric demagnetization tensor, the volume-average ofD ^P (r) - D ^P (r) position dependent demagnetization tensor, implicitly defined in eq. (2.12) - D _xx,D _yy,D _zz demagnetization factors, the diagonal elements ofD. Values for ellipsoids are defined in eq. (2.15) - F ^(m) magnetic force exerted on a body in a magnetic field - H _i ;H _i magnetic field inside a ferromagnetic body; magnitude ofH _i - H ₀;H ₀ magnetic field applied by external sources; magnitude ofH ₀ - h _i ;h _ix,h _iy Cartesian components of dimensionless internal magnetic field,h _i =H _i /H ₀ - I moment of inertia tensor - k geometric parameter for hydrodynamic resistance of a body rotating in a Newtonian fluid given in eq. (2.3) - L ^(h);L _z ^(h) hydrodynamic torque exerted on a rotating body; thez-component of the hydrodynamic torque - L ^(m);L _z ^(m) magnetic torque exerted on a magnetic body in a magnetic field, eq. (2.10); thez-component of the magnetic torque - M;M the magnetization, or dipole moment density, of a magnetic material; the magnitude ofM - M _s the saturation value ofM, approached by all ferromagnetic materials asH _i becomes large (figure 3) - m _s the dimensionless saturation magnetization,M _s/H ₀ - r position vector of a point within a ferromagnetic body - s dummy integration variable in eq. (2.5) - t time - U magnetoquasistatic potential energy of a magnetic body in a magnetic field, given in eq. (2.8) - u curve-fitting variable in eq. (4.1);u = logH _i - V volume of a magnetic particle; for an axisymmetric ellipsoid,V = (4/3) ab ² - x, y, z rectangular coordinate axes fixed in the ellipsoid (figure 1) - angle of inclination of the major axis of the ellipsoid with respect toH ₀ - _s viscosity of the Newtonian suspending medium - µ ₀ the magnetic permeability of free space,µ ₀ =4 · 10^–7H/m - _MV the magnetoviscous time constant, a characteristic time for a process involving a competition of viscous and magnetic stresses - the magnetic susceptibility of a magnetic material, = M/H _i - ; _z angular velocity of a rotating body; angular velocity about thez-axis of an ellipsoid, _z=–d/dt     

Different behaviors of concentrated polysaccharide systems in large-amplitude oscillating shear fields

Two series of oscillatory flow tests were carried out on two concentrated polysaccharide systems (hydroxyethyl guar gum and scleroglucan) in order to evaluate the possibilities offered by large-deformation techniques for the discrimination between different classes of macromolecular systems (i.e., entanglement networks and weak gels). Frequently and strain sweeps, as well as combined steady and oscillatory shear tests were performed to analyze the influence of strain amplitude and of the superposed parallel shear rate on the dynamic properties of these materials.     

Experimental studies of suspended ferromagnetic fibers in a magnetic field

The analysis of the rotation of a ferromagnetic ellipsoid suspended in a Newtonian fluid and subjected to a uniform magnetic field is extended to include a long, slender cylindrical fiber which is magnetically saturated. Experimental observations of rotating nickel cylinders with aspect ratiosL/D ranging from 5 to 40 agree with the theoretical predictions that: (1) the proper magnetoviscous time constant for the motion is _MV = _s/µ ₀ M _s ² , (2) larger fiber aspect ratios result in considerably longer orientation times; and (3) the strength of the applied external field has only a slight effect on the overall fiber rotation, and has no effect on the maximum angular velocity achieved. Quantitative agreement of theory and experiments is obtained for fibers withL/D 20; for the shorter fibers, the theory tends to overpredict the fiber rotation rate by as much as 30%. D diameter of the cylinder - D ^P (r) position-dependent demagnetization tensor, implicitly defined in eq. (2.5) - D _xx,D _yy,D _zz volume-averaged demagnetizing factors for an ellipsoid equivalent to a uniformly magnetized cylinder, defined in eq. (2.6) - H _i ;H _i magnetic field inside a ferromagnetic body; magnitude ofH _i - H ₀;H ₀ magnetic field applied by external sources; magnitude ofH ₀ - k geometric parameter in the hydrodynamic resistance of a body rotating in a Newtonian fluid, eq. (2.2) - L length of the cylinder - L ^(h);L _z ^(h) hydrodynamic torque exerted on a rotating body; thez-component ofL ^(h) on the cylinder - L ^(m);L _z ^(m) magnetic torque exerted on a magnetic body in a magnetic field, eq. (2.4); thez-component ofL ^(m) on the cylinder - M the magnetization of a magnetic material - M _s the saturation magnitude ofM, approached by all ferromagnetic materials asH _i becomes large - r position vector of a point within a ferromagnetic body - V volume of a magnetic particle - x, y, z rectangular coordinate axes fixed in the cylinder according to figure 1 - angle of inclination of the axis of the cylinder with respect toH ₀ - shear rate - small parameter of slender body theory,=1/ln (2L/D) - _s constant viscosity of the suspending fluid - µ ₀ the magnetic permeability of free space,µ ₀=4 · 10^–7 H/m - _MV the magnetoviscous time constant, a characteristic time for a process involving a competition of viscous and magnetic stresses - ₁ the first normal-stress coefficient - ; _z angular velocity of a rotating body; angular velocity of a cylinder about thez-axis, _z =– d/dt     

On the relaxation of shear and normal stresses of viscoelastic fluids following constant shear rate experiments

By means of a cone and plate rheometer the relaxation of the shear stress and the first normal stress difference in polymer liquids upon cessation of a constant shear rate were examined. The experiments were conducted mostly in a high shear rate region of relevance for the processing of these materials. The relaxation behavior at these shear rates can only be measured accurately under extremely precise specifications of the rheometer. To determine under which conditions the integral normal thrust is a convenient measure for the relaxing local first normal stress difference the radial distribution of the pressure in the shear gap was measured. The shape of relaxation of both the shear stress and the first normal stress difference could be closely approximated for the entire measured shear rate and time range by a two parameter statistical function. In the range of measured shear rates, one of the parameters, the standard deviationS, is equal for the shear and the normal stress, and is independent of the shear rate within the limit of experimental error. The second parameter, the mean relaxation timet _50, of the shear stress andt _50, of the first normal stress difference, can be calculated approximately from the viscosity function and only a single relaxation experiment.     

A boundary element simulation of the unsteady motion of a sphere in a cylindrical tube containing a viscoelastic fluid

A boundary element method is used to simulate the unsteady motion of a sphere falling under gravity along the centreline of a cylindrical tube containing a viscoelastic fluid. The fluid is modelled by the upper-convected Maxwell constitutive equation. Results show that the viscoelasticity of the liquid leads to a damped oscillation in sphere velocity about its terminal value. The maximum sphere velocity, which occurs in the first overshoot, is approximately proportional to the square root of the Weissenberg number when the ratio of the sphere radius to the tube radius is sufficiently small. Particular attention is also paid to the wall effects. It is shown that a closer wall reduces the oscillatory amplitude of the sphere velocity but increases its frequency. The results suggest that the falling-ball technique, which is now widely used for viscosity measurement, might also be used for the determination of a relaxation time for a viscoelastic fluid.     

Bimodal model of slurry viscosity with application to coal-slurries. Part 2. High shear limit behavior

It is shown that in a truly bimodal coal-water slurry the hydrodynamic interactions between the coarse particles impose on the fine fraction a shear rate higher than that applied externally by the viscometer walls. A semi-empirical function of the coarse volume fraction is obtained for this correction factor to the applied shear rate. The derivation of this shear correction factor is based on lubrication concepts and introduces the maximum packing fraction,ø _m, at which flow can take place.ø _m is obtainable from a simple dry packing experiment. It is shown that the contribution of the coarse particles to the viscosity rise can be successfully described by a viscosity model employing the same concepts used to derive the shear correction factor. The bimodal model is applied in the high shear limit to polymodal coal slurries with a continuous particle size distribution. In the model, the contribution of the coarse particles to the viscosity rise is taken from separate viscosity measurements for the coarse coal particles, while the contribution to the viscosity of the fine coal particles is taken to be that given by the measured viscosity of colloidal suspensions of monomodal rigid spheres. It is shown that there is a ratio of coarse to fine fraction volumes in the continuous size distribution, corresponding to a specific separating particle size, for which the measured viscosities of the polymodal slurries match almost perfectly over the whole solids volume fraction range with the viscosity values obtained using the bimodal approach. The match is found to be relatively insensitive to the precise value of the separating particle size.     

The squeezing flow of a model suspension fluid

The paper is concerned with the squeezing flow of a model suspension fluid. The numerical solution obtained by a time-dependent Boundary Element Method is compared to an asymptotic solution at large radius. It is found that the kinematics are Newtonian in character, and the fibres quickly align themselves radially. Consequently, the squeezing force is only weakly dependent on the initial orientations of the fibres and the device can be used for measuring the effective viscosity of the suspension. The effective viscosity found from the squeezing flow agrees surprisingly well with experimental data and numerical data derived from the falling sphere geometry at low volume fractions ( < 0.1).     

Shrinkage effect on polyethylene based blends

An experimental study of the shrinkage effect in low density-high density polyethylene blends, high density-high density polyethylene blends and low density polyethylene-ethylene vinyl acetate copolymer blends is presented. Viscosity measurements are also included.The strain recovery is analysed as the addition of the interfacial tension effect and the elastic recoverable strain. The results are compared with those obtained for other polyolefin blends observing that in the case of compatible systems as low density-high density polyethylene and high density-high density polyethylene blends the strain recovery and the viscosity follow approximately the additive rule. In the case of more incompatible systems like polyethylene-ethylene vinyl acetate copolymer, polyethylene-polypropylene and others very large recoveries were observed. The formation of minifibres of the dispersed phase during the preparation of the samples is suggested as responsible for the obtained results.     

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.     

Analysis of incompressible three-dimensional flows using the concept of stream tubes in relation with a transformation of the physical domain

This paper presents an analysis of three-dimensional flows based on the concept of stream tubes and streamlines. In the case of incompressible liquids, the physical domain can be transformed, under certain assumptions, into a cylinder where the streamlines are parallel straight lines. In contrast to classical methods, the unknown to be determined is the transformation between the two domains. This analysis generalizes the formulation already proposed for plane and axisymmetric flows.     

Comments on the paper “On the derivation of the constitutive equation of the simple fluid from that of a simple material” by R.R. Huilgol

Huilgol's paper [1] is discussed critically. It is shown that there are significant errors in each section of the paper. His argument depends critically on a refusal to recognize that when a statement is made that a functional of a tensor function is invariant under a group of transformations, these transformations must be constants and that otherwise the statement is meaningless. Furthermore, he consistently ignores the contradictions to which this refusal leads.     

Stokes' first problem for linear viscoelastic fluids with finite memory

We consider Stokes' first problem for a viscoelastic fluid. The memory of the fluid is truncated to a finite time interval and discontinuities in the stress relaxation modulus or its derivatives are allowed at the point of truncation. We investigate secondary waves which are generated by the interaction of these singularities in the memory with earlier waves.Dedicated to Prof. Hanswalter Giesekus on the occasion of his retirement as Editor of Rheologica Acta.     

On the steady simple shear flows of the one-mode Giesekus fluid

Analytical solutions for the plane Couette flow and the plane Poiseuille flow of the one-mode Giesekus fluid without any retardation time have been obtained by considering the domain of definition for each of the two branch solutions which arise due to the presence of the quadratic stress terms in the constitutive equations. For each fixed value of the mobility parametera, the limiting value of the Weissenberg number for the upper branch solution, i.e., the physically realistic solution is determined in terms of the corresponding dimensionless shear stress for the plane Couette flow and in terms of the corresponding dimensionless pressure gradient for the plane Poiseuille flow. In the case of the plane Couette flow, it is shown that fora falling in the range 0a1/2 only the physically realistic solution exists while for 1/2<a 1 a nonphysical solution coexists with the realistic one. In the case of the plane Poiseuille flow, it is shown that the non-physical solution cannot even exist around the center plane of the channel, and the effects of the mobility parameter and the dimensionless pressure gradient on the flow variables are investigated. Possible extensions of the present approach to other steady simple shear flows with and without the introduction of the retardation time are also discussed.     

The concept of a rotational rheometer with helical screw impeller

The paper is concerned with the development of an on-line rheometer for suspensions with settling particles. A new concept of a rotational rheometer with helical screw impeller is discussed. On the basis of model considerations a procedure for the flow curve determination of the investigated fluid has been proposed. Contrary to previous attempts, in the developed procedure no calibration of the instrument is necessary. The preliminary experimental verification of the developed approach proved its applicability.