A recent review publication  presented an extensive and comprehensive asse- ssment of the phenomenological relations of Poissons ratios (PRs)   to the behavior and responses of contemporary materials under specific loading con- ditions. The present review and analysis paper is intended as a theoretical mechanics complement covering mathematical and physical modeling of a single original elastic and of six time and process (i.e. path and stress) dependent viscoelastic PR definitions as well as a seventh special path independent one.
Viscoelasticity had its origins with the seminal studies of Kelvin, Maxwell, Voigt, Boltzmann and Volterra during the period of 1865 to 1913  -  . However, it did not start to flourish until the mid 1940s when general constitutive relations and elastic/viscoelastic correspondence principles (EVCP) were starting to be developed  -  . These works were followed by the works in the partial listings of Refs.  -  . There remain two major contentious issues regarding viscoelastic PRs  -  and both, while summarized in this section, are analyzed and discussed in detail subsequently.
1) PR dependence on time and on stress and/or strain history
2) Contributions of the loading phase in determining material properties
The importance of starting loading phases on subsequent viscoelastic responses and consequently on material characterization (moduli/com- pliances) has been demonstrated through analytical simulations in  ,  and through experimental investigations in  -  . Less restrictive conclusions based on the “en-times-rule” have been reached in     based on comparisons of strain responses. While such approxi- mations may be satisfactory for long time stress analyses, the neglect of loading phases produces significant errors in viscoelastic characterizations   . In  it is stated that measurements of the loading times are subject to errors with larger errors corresponding to shorter times. The analysis is based on constant stress loadings without inclusion of loading cycles in the data analysis.
Finally, it is important to note that experimental determinations of viscoe-
Table 1. Viscoelastic Poisson’s ratios and the EVCP.
NOTES: All PRs depend on time, stresses and stress-time histories, and these results apply equally to homogeneous and nonhomogeneous media. None of the PRs represent universal material properties as they are path dependent, i.e. INot in proper form in FT space for EVCP. See Table 3. IISee the Section “The physics of time invariant stresses, etc.” for the physical difficulties associated with achieving constant strains or stresses or constant rates for. IIIInverse FT has double convolution integral in space. See Eqs. (32) for some physically unattainable forms, which do not lead to EVCP. IVBased on natural strains and, therefore, nonlinear. VSimilar EVCP difficulties as Case I. VIIf and only if. See Section “The physics of time invariant...” VIIThe Class VII PRs are path and strain independent and only functions of time.
lastic PR time functions are based on specific loading or strain histories and, consequently, are limited to those conditions and are non-exportable to other loadings, since viscoelastic PRs are path dependent and hence non-unique as seen in Figure 1. In this figure, distinct PRs are displayed that were obtained from experimental strain data under different loading conditions on identical material specimen. Similar multiple PR time functions were reported in  .
The Class VII PR definition avoids loading path dependence problems, but still leads to double convolution integrals in the time space that are extremely CPU intensive.
A recent paper  , while not directly related to elastic PRs, analyzes the effects of geometry on changes in stress directions and demonstrates how usual
Figure 1. Case I Poisson’s ratios based on creep and relaxation experiments  .
tension strains can reverse direction, a phenomenon which the authors denote by “negative stiffness.”
In  an analysis is presented on how to interpret loading over-and under- shoots when extracting viscoelastic material properties from experimental data. An ancillary benefit is provided by demonstrating the importance of including loading phases in any analysis to determine viscoelastic material properties. Thus reinforcing the theoretical quasi-static and dynamic analysis results of   that detail analytically the important effects of loading start-up phases on material characterizations.
However, as it is shown in  -  the entire issue of viscoelastic PRs can be by-passed by formulating 3-D analytical and experimental material chara- cterizations in terms of 1-D loadings, and shear and bulk relaxation moduli or creep compliances. These characterizations are accomplished in either real time or Fourier or Laplace integral transform spaces and specifically include the associated loading cycles and totally steer clear of viscoelastic PR use. They are based on measurements of time, the 1-D stress and the 2-D strains and lead to universal viscoelastic material characterizations.
Comment 1. Aside from Poisson’s great native genius, what made his namesake ratio universally accepted and most useful is that for all but seven decades of its two century existence the application were limited to elastic materials where for each of them it is a distinct constant. Unfortunately, nature was significantly less kind and far less magnanimous when it comes to viscoelasticity.
Detailed analyses and discussions of the above topics are presented in subsequent sections.
2. Linear Analysis
2.1. General Considerations
Elastic and viscoelastic constitutive relations based on moduli and/or compli- ances are derived from fundamental principles, i.e. the thermodynamic laws pertaining to energy and entropy  . Thus, moduli and compliance parameters are thermodynamic derivatives in the same manner as for instance specific heats, etc.
In linear elastic materials moduli and compliances are time, path and stress independent and stress-strain relations involving PRs must be and are derived from these basic equations. Furthermore, in elasticity where PRs are constants and constitutive relations are algebraic, PRs are unique material properties and derivations are a straight forward matter.
In viscoelastic media where PRs are time, stress and stress history (path) de- pendent and where stress-strain relations involve time integrals multiple com- plications arise. For one, numerous PR definitions are possible and derivations may be carried out in either time space or in the integral transform space yielding in some instances mutually exclusive results and PR definitions. Secondly, since viscoelastic media dissipate energy loading histories take on prime importance and must be considered  -  .
Figure 2 is a schematic representation of a loading path and moduli at various temperatures. The indicated times represent the following: is the time when relaxation begins, the time time when loading reaches the desired steady- state condition and the beginning time when the fully relaxed modulus starts. Clearly, the time at which the same loading path intercepts the modulus varies with temperature thereby eliciting distinct responses. Secondly, and equally importantly, if the buildup loading path from is altered or if, then the viscoelastic responses will change accordingly  -  . In  it is stated that transient rise times do not influence characterization if data is collected after ten times the rise time. However, such an approximation ignores the relation of the rise time to the time when relaxation and/or creep begin-see Figure 2 and is not born out by either the subsequent exact analytical simulations of  and  or the experimental results presented in  -  . The rise time phenomenon is analyzed in detail in a later section.
The PR analyses and formulations are generalized to non-homogeneous linear viscoelastic materials, however the main focus is on PR time and process dependence. The results apply equally to homogeneous viscoelastic media and PRs are path dependent for either homogeneous or non-homogeneous viscoelastic media.
Table 1 and Table 2 depict the various PR definitions of the six listed categories (classes) in real time space or in integral transform space. All are implicitly path (stress history) dependent, while Classes II and VI are implicitly so being respectively functions of constant strain and of constant stress loadings. Class IV involves a nonlinear definition involving natural strains and thus inapplicable to linear viscoelasticity.
Thus, while some of the viscoelastic PR models are defined in terms of con- sistent mathematical transient loading models, which are physically unattainable, they fail to translate into realistic and observable physical entities.
It will be shown that PRs, although extremely important in elastic material characterizations because there they are constants for each material, are not fundamental viscoelastic quantities in the same sense as relaxation moduli and creep compliances. Instead, they are multiply defined process dependent quan- tities as detailed in Table 1 and Table 2 and not interchangeable between divers loading conditions. Consequently, viscoelastic material characterizations are best expressed in terms of moduli or compliances     .
2.2. Elastic PRs and Constitutive Relations
Since elastic conditions form the ICs of any viscoelastic formulation and to properly establish EVCPs, it is necessary to first review pertinent elastic PR developments. Consider a Cartesian system with coordinates, with and where. The Einstein tensor notation is applied with pairs of repeated indices indicating summation and underlined ones signifying no summations. Isotropic and anisotropic elastic material behavior may be characterized in terms Young's moduli, moduli, shear and bulk moduli, compliances and PRs      . In linear homogeneous elasticity each of these parameters are constants uniquely associated with elastic material behavior. The PRs in question are the original legacy Class I ones
Figure 2. Schematic of typical viscoelastic moduli and loading pattern.
Table 2. Viscoelastic Poisson’s ratio definitions in time and Fourier transform (FT)  spaces.
defined1 in  and  as the ratio of two isothermal orthogonal strains2 or
The anisotropic homogenous elastic constitutive relations read
with where is a reference temperature when the thermal trains are zero. In the isothermal cases.
In particular, for a 1-D loading and all other
with the equal to constants that are specific to each elastic medium. Gen- erally except when.
However, for non-homogeneous elastic materials the usually constant 1-D PRs become spatial functions with values depending still on distinct material properties but also additionally on the position within the elas- tic body, such that Eqs. (3) turn into 
For linear homogeneous elastic media, it can be readily proven that
see  to  . However, in non homogeneous elastic bodies no such restrictions apply  . Nor is the PR an absolute constant for each iso- tropic isothermal non-homogeneous elastic medium, but rather varies spatially in a manner similar to the corresponding moduli and compliances  .
Table 3 and Figure 3 display the various moduli for the isotropic isothermal homogenous elastic media based on Hooke’s law. These are seen to vary from an incompressible material at to a rigid one in shear with only tension or compressive strains at a PR of 1. Parenthetically, one must add that these limits do not apply in linear viscoelasticity, except at for those viscoelastic media that have elastic initial conditions.
It can be further demonstrated that for linear elastic media, the following hold true
and similar relations for the Lamé parameters (constants)
Alternate elastic constitutive relations are obtained by inverting the set (2) which yields
Table 3. Relations between isotropic linear elastic moduli and Poisson’s ratios.
Figure 3. Elastic moduli and Poisson’s ratio.
2.3. 3-D Hooke’s Law and Poisson Ratio Induced Nonlinearities
In linear 1-D elasticity, the isotropic isothermal Hooke’s law simplifies to
where is the elastic compliance.
Furthermore, in conservative linear systems Maxwell's reciprocal theorem    applies and the loading order is immaterial for linear elastic media. Therefore, in 2 or 3-D linear isotropic elasticity the strain is deter- mined by assuming or prescribing that is applied first and then followed by the other two normal stresses. As derived in  , geometric considerations reveal the inherent PR nonlinearities based on a linear elastic Hooke’s law
Taking the FT of these linearized isotropic elastic constitutive relations yields
Comment 2. It should be noted that in most current viscoelastic materials, such as high polymers, rubbers, etc., the inequality (linearization term) in (11) and (12) is often violated because of low and values comparable to and sometimes smaller than. The inclusion of this approximation is crucial to linear PR elasticity and the admissibility of the EVCP based on PRs. On the other hand, constitutive relations in terms of moduli and/or compliances are formulated from fundamental deviatoric and volumetric laws deformation without the need to invoke the linearization assumption.
The net effect of this linearization term is that linear elastic materials possess linear modulus/compliance constitutive relations, but potentially nonlinear ones when PR formulations are included.
However, the above linearization must be enforced since the EVCP depends on integral transforms, which can only be applied to linear systems. Additionally and separately, in order to formulate the EVCP in terms of PRs it is necessary to find a term by term and parameter by parameter viscoelastic match in the FT space to the elastic transformed relations (12), with a corresponding valid and plausible viscoelastic set of constitutive relations in the time space.
Alternately, the troublesome nonlinearity can be totally avoided if Hooke’s law is derived from fundamental deviatoric (shape change) and volumetric (volume changes) constitutive relations and by applying the three normal stresses simultaneously, to whit
Substitution of (14) into (13) yields
and one equation for each of the other two normal strains. Eq. (15) and its other two cousins are unambiguously linear and do not suffer from the nonlinearities of the PR relations (11). Furthermore in the latter formulation all three normal loads are applied simultaneously which is of prime importance in dissipate systems such as viscoelastic or plastic media.
2.4. Viscoelastic PRs and Constitutive Relations
Whereas in isotropic isothermal homogeneous linear elasticity the PR values are constants and material specific, the circumstances are far different for their vis- coelastic cousins under the same specifications.
Consider the general linear anisotropic isothermal viscoelastic relations  -  , which are the counterparts to Eqs. (2)
Consequently, the Class I viscoelastic PR for an isothermal process when becomes
Note that even in the simplest 1-D isotropic problem when only one stress derivative is present and if it is equal to a constant value3 for and if the material is homogeneous, then the viscoelastic PR is still a time function, i.e.. This is due to the fact that for an isotropic homogeneous ma- terial
Comment 3. In other instances when the stress or strain derivatives are time functions, they are “trapped” in the integrals and the PRs become time, stress and stress history dependent, i.e. process dependent     . Therefore, the viscoelastic PRs are no longer pure unique material property characterizations in the sense of moduli and compliances, and are, therefore, unsuitable as material descriptors. This applies to Classes I throu- gh VI. ClassVI (1-D constant stress) and Class II (1-D constant strain) are special limited cases, and in reality none of the PRs are exportable from one loading condition or from one Class to another including time independent loadings.
In the very special case when, the shear modulus is related to the bulk modulus through a constant multiplier with the requirement that their relaxation times are equal  . Real materials gene- rally exhibit except for the auxetic ones where the inequality reversed. In either instance Eq. (18a) is satisfied.
The relaxation moduli and creep compliances are as a matter of convenience generally represented by Prony series  of the type
with similar expressions for and.
Fourier transform (FT)  are introduced in preference to Laplace trans- forms (LT) as they are more inclusive and extends from, where all state variables are at rest, to the current time at which the solution is sought. Con- sequently, in general analyses such as those presented here ICs at need not be included and need not be specified until solutions of particular problems are undertaken. FTs were first introduced to viscoelasticity in order to derive the integral EVCP in Ref.  and are defined here as
The FTs are equivalent to the two sided LTs as stated by the Rodrigues formula and provided all integrals exist these transforms are related to each other by  
Application of FTs converts the convolution integral isothermal constitutive relations (16) to
Class VI PR expressions are derived in  , while all other Categories have defined expressions as seen in  . The Class VII PRs have their genesis in a definition starting in the FT space, based on a form of the constitutive relations that can be interpreted as supporting an EVCP. However, they are devoid of physical meaning and the stress-strain relations do not transform into a form in time space that can be derived independently from thermodynamic first prin- ciples. They represent a group of material property time functions independent of stress and strain histories. Nevertheless, their utility is far less than the relaxation moduli and creep compliances that define them because they are specific to a limited 1-D loading condition from which one cannot generalize the PR functions.
The Elastic-Viscoelastic Correspondence Principle or Analogy (EVCP)4
The EVCP, when applicable, is a powerful protocol for solving linear visco- elastic problems based on similarities between integral transform (IT) elastic and viscoelastic formulations and solutions. EVCP comes in two varieties: separation of variables  and the considerably less restrictive one based on integral transforms  . The IT EVCP consists of substituting complex viscoelastic moduli/compliances for their elastic equivalences in Laplace (LT) or Fourier (FT) transforms elastic solutions. This means that any analytical and some numerical elastic solutions are easily convertible to their equivalent viscoelastic solutions. In those instances where material properties do not conform to con- volution integral constitutive relations, approximate EVCPs have been derived that are based on piecewise convolution integral representations   .
The EVCP can be formally stated as:
1) Given an elastic body with the FT of its general solution
2) Then for viscoelastic convolution5 integral constitutive relations given by (16) and with identical boundary conditions, the viscoelastic solution is
are universally valid provided the transforms and their inverses exist.
Comment 4. The EVCP protocol  for rigorously establishing Eqs. (23) to (25) is to first derive the pertinent viscoelastic governing relations and their generalized constitutive relations and solutions in the time space and then take their FT and compare these transforms to the FTs of the corres- ponding elastic solutions.
In order to properly highlight the difficulties of generating an EVCP involving viscoelastic PRs, it is necessary to first review the process by which Eqs. (2), the 3-D general Hooke’s law with PRs is derived. Elasticity is a conservative system and its constitutive relations are algebraic. Whereas the viscoelastic system is a dissipative one with memory leading to integral differential constitutive relations in time, such as (16).
However, the linearization described in Section 3.3 must be enforced since the EVCP depends on integral transforms, which can only be applied to linear systems. Additionally and separately, in order to formulate the EVCP in terms of PRs it is necessary to find a term by term and parameter by parameter viscoelastic match in the FT space to the elastic transformed relations (12), with a corresponding valid and plausible viscoelastic set of constitutive relations in the time space.
In a non-conservative system such as a linear viscoelastic body, the luxury of choice of free sequencing of loads is no longer available, as different loading schedules result in distinct responses, i.e. Maxwel’s reciprocal theorem is inapp- licable. This is due to the material's memory, and consequent energy dissipation, which is expressed by integral type constitutive relations rather than the elastic algebraic ones. The viscoelastic ones also do not allow for the legitimate types of elastic algebraic approximations and manipulations of Eqs. (10) and (11), thus dictating the need for different approaches  .
Since the correspondence principle has been proven inapplicable when PRs are involved    , the usual elastic relations for shear and bulk moduli have no viscoelastic counterparts in terms of PRs, i.e.
but without PRs, the expressions
and Eqs. (25) are valid.
However, if is substituted for in (26) then the resulting relations
are valid subject to the limitations of the non-uniqueness of the viscoelastic PRs due to their stress and path (stress history) dependence. Difficulties with gen- erating a proper EVCP in terms of PRs stem from the following concepts, which will be established below in detail, see Eqs. (32).
1) The time and separate stress history dependence of the viscoelastic PRs, i.e. experimental conditions do not necessarily reflect in-service situations.
2) The inherent assemblage of individual material contributors in the elastic constitutive relations when PRs are involved calling for a triple product of ma- terial property parameters and stress tensors.
3) The viscoelastic PR triple product while obligatory in the Fourier transform space for establishing an EVCP has limited physical counterparts in the real time space nor any relation to the Boltzmann superposition principle  , except for one model (32b).
In terms of moduli/compliances the 3-D Hooke’s law and the corresponding linear viscoelastic constitutive relations are both derivable from first principles  to read for elastic media
However, when the isotropic Hooke’s law is assembled with elastic PRs, it is in an almost ad hoc manner of one direction at a time to yield6 
Specifically, the pitfalls in producing a coherent EVCP arise from the triple product elastic combinations and their FTs7
that requires a viscoelastic triple FT product counterpart
, which is generated in the integral transform space (FT or LT). However, in order to achieve any physical and mathematical credence this triple product needs to be generated from real time combinations to conform with its prescribed FT inversions. Therefore, starting with the need- ed triple PR FT product, the required isotropic isothermal constitutive relations in the time space must be the inverse FT
where represent any class PR for which the FT triple pro- duct can be generated8. The above four expressions indicate how the constitutive relations must be formulated in the real time space so that the desired FT triple products can be realized. If the isotropic 1-D elastic protocol of Eqs. (10) is used in assembling equivalent viscoelastic PR constitutive relations, then the constru- ction becomes
It then follows from this inverse heuristic approach that
There are three Classes of PRs that can produce the form (34), namely Class II, III and VI. However, neither Class II or VI produce a physically valid expre- ssion in the real time space. Class III, as pointed out in earlier, has only one convolution form out three possible ones that produce a physical counterpart (32b).
The definition of the Class VI viscoelastic PR leads to
Note that for the Class VI PRs the loading path is defined by prescribing it as. However, it must be remembered the this Class is based on the absolute neglect of the loading cycle and that a time independent stress is not realistically achievable for all times.
Some of the ambiguities arising from the Class III PRs can be resolved, i.e. removed, by defining a Class VII PR as
It is the only PR class that is entirely free from strain, stress and path restri- ctions since its definition is not based on the legacy Class I   , formulation or on any of its derivative definitions of Classes II to VI.
Comment 5. The Class VII is an artificially defined isothermal viscoelastic PR devoid of strain or stress associations and hence path independent while still remaining a time dependent universal martial property. It obeys the construction mandated by Eqs. (32) necessary for the existence of EVCPs. This Class PR should, therefore, be the one used if it is desired to express constitutive relations in terms of PRs. However, since this PR form nece- ssitates a priori knowledge of two compliances that completely characte- rize the isotropic viscoelastic medium, it is fair to ask why bother with the PR. Furthermore, since Class VII PRs do not involve strains they are not experimentally measurable.
At least three more than rhetorical questions remain:
1) What, if any, is the physical meaning of the Class VI PRs,?
2) What is their relation to the strains?
3) What is the physical meaning, if any, of the double convolution integrals of Eqs. (32a) and (32c)?
The answers to all three questions are “none” indicating that they cannot be obtained directly through experimentally determined strains and stresses. Fur- thermore, due to the presence of a triple product in the above constitutive rela- tions these PRs are not a portal to EVCP formulations. Even though the Class VI PRs appear as pure universal material property descriptors their utility function is less than that of its parts, i.e. the compliances, since these PRs are first and foremost process dependent and, therefore, restricted and limited in scope in this case to time independent 1-D stresses. (See the next Section on starting loading sequences to achieve time independent stresses, strains, their time rates, etc., in a dissipative material with memory effects.)
There is a certain irony accompanying the creation of the Class VI PR based on a time independent 1-D stress field
The intent is to produce a viscoelastic look alike to the structure of the elastic constitutive relations containing a mixture of moduli, stresses and the original Class I elastic PR  . What emerges as a Class VI PR is a function independent of strains and their history and one that is only dependent on certain moduli and compliances. Of course, in the elastic case the same path following Eqs. (31) leads to for and all other9
For a viscoelastic medium, can be similarly interpreted as
with the real time convolution integral property
It should be noted that the Class VI PR is a special case of the general Class III PR with and. Neither Class III nor Class VI PRs have a physical definition whose values could be measured experimentally in the time domain. While it is certainly possible to measure the strains in real time, the time functions associated with this PR class are complicated and the evaluation of the integral in (41) is difficult. The second integral of (32) clearly represents the proper form of the triple FT product involving PRs in the time domain. The character of the time function defining can be deduced from the inversion of the FTs, Eq. (36). Since the Prony series for and are known expressions determined from generalized Kelvin or Maxwell models. However, since Eq. (36) requires the a priori determination of the compliances, one can justifiable inquire why bother with when the characterization is complete with the knowledge of either set of the or. Particularly, since detailed knowledge of these non-universal PRs still does not produce a useful EVCP and either the moduli or the compliances do so handily.
In the viscoelastic constitutive relations there appears to be no clear protocol to derive directly in real time space either of the first two versions of double convolution integrals displayed in (32) that can be based on thermodynamic first principles or any mechanical or electrical simulation model generated under the three fundamental thermodynamic laws or on the Boltzmann superposition principle  . Only through mimicking Hooke’s law in the FT space can Eqs. (32) be so postulated in an inverse fashion. In a larger sense then, the intro- duction of PRs into viscoelastic stress-strain relations can be interpreted as an artificial act of faith not based on first principles, i.e. with PRs related to thermodynamic derivatives, or even to achievable physical characterizations and experiments, as will be seen subsequently.
Comment 6. None of the five PR categories defined in    fit these needed real time prescriptions, nor does the sixth class fit the cons- traints imposed by (32). Consequently, it would appear that no universal EVCP based on PRs can be rationally deduced or justified.
The expressions (26) are valid in a very limited sense for Class II PRs, but are only applicable there and cannot be exported to the general Class I expressions (24) with and, therefore, are of no use beyond the one and only case of a time independent strain component, say (see next Section). Additionally, any expressions for and derived from (30) based on Class II PRs apply only to conditions of time independent strains or and cannot be generalized to any other strain field. Their complete lack of universality makes such modulus expressions useless from a practical point of view. Similar universal compliance (and modulus) relations may be found in Refs  -  .
2.5. Comments on PR’s Dependence on Time
The following is a partial list of the cardinal difficulties that are directly attribu- table to their time dependence and that are associated with the use of viscoelastic PRs:
1) At least seven mutually independent PR definitions are available.
2) Difficulties stemming from the improper assembly of constitutive relations containing PRs in the time space resulting in combinations
for which no integral transform can obtained that leads to a proper EVCP expression. See Eqs. (33) for proper time space forms required to realize the needed double integral transform products and triple time function products (32).
3) Conversely independent assemblies in the integral transform space (FT and LT) lead to some combinations displayed in Eqs. (32) such that the triple pro- ducts have no physical counterpart in the real time space.
4) Additionally and separately viscoelastic PRs are also stress and stress history (i.e. path) dependent.
2.6. The Physics of Time Invariant Stresses, Strains and of Their Derivatives
The Influence of Constant Stress/Strain Models on PRs
The Class II and VI PR formulations are predicated on the mathematically correct prescription of constant 1-D strain (II) or constant 1-D stress (VI), such that
where is the Heaviside step function with defined properties 
The special function has a discontinuity at, but possesses a derivative i.e. another special function-the Dirac delta function or impulse fun- ction 
Consequently the Class II and VI mathematical models are well defined but loading (path) dependent. Unfortunately the actual physics of this loading problem are considerably less forgiving, as shown analytically in   and experimentally in  -  . An instantaneous displacement requires an infinite acceleration and an associated infinite force in zero time with similar physically unattainable conditions for constant stress and strain rates. Further as analytically proved in  , even the analytical modeled application of Dirac delta function loads lead to singularities in time derivatives of displacements for linear and nonlinear elastic and viscoelastic media.
In addition, testing machine crossheads have inertia and require finite times to develop displacements and loads as demonstrated in  . In Section 4 data is displayed for experimental loading cycle results of high polymer tension coupons obtained from four different INSTRON machines. The detailed analyses show that the rise times to constant strain rates actually occur considerably after creep and relaxation have begun at 0.001 to 0.01 sec, consequently rendering the mathematical model of Class II and VI physically unrealizable and, therefore, making the associated PR functions physically unachievable and, therefore, unacceptable. Furthermore, in   analytical studies clearly show the heavy influence of various loading protocols on the accuracies of viscoelastic material characterization.
2.7. Penalties Associated with Instantaneous Loading Models
2.7.1. Dynamic Loadings
The instantaneous loading concept can be modeled mathematically  through the use of Dirac delta and Heaviside functions  and  while the gradual loading process to achieve a steady or time dependent loads requires more than one functional representation For instance, consider the case of dynamic loadings leading to wave propagations. Figure 4 illustrates the experimentally measured 1-D effects of an impact load on a bar at (0,t). In particular note the green curve which is the impact force. The other two curves represent accelerometer readings at two distinct points at. The impact loading, then, can be represented schematically as in Figure 5 and by Eqs. (46) to (48)
Figure 4. Impulse force and accelerometer responses  .
Figure 5. Initial instantaneous and gradual loading phase schematic.
with four conditions on
The inherent difficulties associated with the and functions are due to the latter's double values and hence discontinuities at. Consider their troublesome properties
plus those of (45) and the more beneficial finite one
This effectively means that for the mathematical model, if the IC properly consists of an instantaneous impulse force
then it will propagate into the linear and nonlinear elastic and viscoelastic displacement solutions, such that making the temporal deriva- tives that determine velocities and accelerations unbounded at when physically both are finite  .
In view of consistency requirements one should expect no less than an infinite acceleration to produce an instantaneous impulse force.
2.7.2. Quasi-Static Loadings: Creep, Relaxation, “Constant” Strain Rate, SHM
By contrast the quasi-static models differ from that dynamic ones essentially in the fact that after the completion of the time dependent loading cycle they re- main at a constant value. Let represent the set of state variables, such that. Then the model representation become
With the mathematical model, one is again faced with a discontinuity due to the double valued Heaviside function at the temporal origin resulting in
Similar conclusions to those enunciated in Section 3.7.1 apply here as well.
Comment 7. In the final analysis the impulsive loading and sudden displa- cement initial problems are reminiscent of the non-existing sharp edged gusts, instantaneous penetration of atmospheric disturbances and zero time lift and drag buildups encountered in aeroelasticity. All such aerodynamically improbable physical models have been effectively corrected by the introduction of the three distinct delay functions of Theodorsen, Küssner and Wagner  . In equal measures, the solid mechanics physical aspects command that no less be done for the preceding ICs and possible subsequent changes in loading functions.
2.8. Analysis of Time Dependent Experimental Ramp Loadings
Refs.  and  contain analytical simulations that show the various regions of viscoelastic moduli, and conversely of compliances. The simulation studies proved that the ratio and the extent as well as the slope of Region are the most important contributors to optimum viscoelastic material damping properties as seen in Figure 6. Of course, the shape and regions of the relaxation modulus curves are directly dependent on the sensitivity of the experiments that are used to determine such modulus properties. Furthermore, since PRs are related to moduli as seen in Eqs. (28) and (32) they exhibit similar sensitivities.
In Ref.  a systematic analysis of starting transients influence on constitutive relations is undertaken by analytically investigating simulation models for simple tension experiments. Since numerical values of material parameters are
Figure 6. Viscoelastic modulus regions  .
sought, meaningful experiments need to be devised which can be solved analy- tically with symbolic values for the as yet unknown material parameters. Consi- der a “simple” 1-D tension or compression (without buckling) creep experi- ment (Figure 4). In a Cartesian coordinate system with and the loaded direction. Whether or not shear is present at in the direction has no bearing on the loading function formulation, but, of courses, influences internal stress distributions.
The loading function is defined in three time domains as
where is a constant amplitude, i.e. or of (43), and is one of the following functions
with the Heaviside unit step function. While mathematically all three functions are admissible, the only physically “reasonable” function is the third one, although other similar forms are equally acceptable, as it properly defines a gradual load rise from 0 to unity with physically required vanishing slopes at and.
The time necessary to achieve the constant load is dictated by labo- ratory equipment used to induce loadings. While has no predetermined relation to, the time when relaxation begins, i.e. the relaxation modulus (the instantaneous elastic Young’s modulus) for and for, the determination of moduli values are strongly influenced by the relative position of and.
The next two figure from  present creep functions determined by solving the governing relations including the loading phase for a given load case and different values of the loading time. For illustration purposes the strain cur- ves were established prescribing that the material characteristics are known. The results obtained following the proposed procedure were compared to this refe- rence curve. Times were chosen here by evaluating the increase in compared to and are indicated as percentages in the plot legends. Since the values are not significant on their own, they will be identified on the plots by the relative increase in from.
Figure 7 and Figure 8 from Ref.  represent creep curves obtained through analytical simulations by assuming that the loading conforms to an ideal situation, i.e. by replacing by other loading functions as detailed in (55) and (56). If the strain values used when solving the system of equations correspond indeed to Loading Case 1, the creep curve is recovered with good accuracy (less than 0.5% error for terms in Prony series). However, if a more realistic situation is considered (Loading Cases 2 and 3), the material characteristics are determined with some error depending on. In Figure 7 the strains were obtained from Loading Case 2 and all the points were left out of the computation. As expected, the largest deviation appears in the values of. Almost identical results were obtained using strains corresponding to Loading Case 3 and displayed in Figure 8. A number of additional examples may be found in  . The influence of testing machine cross head speeds is discussed in the next Section.
2.9. The Influence of the Loading Phase on Constitutive Relations and PRs
Starting with the linear anisotropic viscoelastic constitutive relations (16) and specializing them to 1-D loading isothermal isotropic conditions yields
For a hypothetical constant stress as defined by (43), Eq. (57) can be integrated to yield
Figure 7. Loading Case 2 creep function determined assuming loading Case 1  .
Figure 8. Loading Case 3 Creep Function determined assuming loading Case 1  .
Eq. (58) is correct if and only if for and and all other for all. Therefore, it does not apply when the loading phase (55) is taken into account since then Eq. (57) is no longer correct and must be properly interpreted.
If the loading cycle (55) is taken into account then Eq. (57) transforms into
Therefore, the error in strain responses generated by disregarding the starting stress transients is given by
and is shown in Figure 9 for experimental data obtained in  . These results show that the effects of the starting transients decay with time but that they also represent important contributions to inaccuracies over significant time periods in excess of, the time to reach steady-state loading conditions.
However, a much more compelling case can be made for the linear relaxa- tion/creep functions or moduli/compliances, which are more pervasive and universal properties than specific strain responses   .
2.10. The Case III PR Caveat
The Class III PR integral transform definition, Table 3, while lending itself more expeditiously to the formation of an EVCP has troublesome idiosyncrasies and a dichotomy of meanings in the time space. Note that for Case III PRs, for 1-D loadings and time independent temperatures the convolution integrals prevail and then and only then
Figure 9. Error in viscoelastic strains by disregarding starting load transients for data from  .
which makes the expression independent of but not of the 1-D loading process (path) that specifically creates it. Eq. (63) does not apply to mul- tiple loadings where it becomes a function of all the pertinent normal stresses, i.e.
for isotropic materials.
The normal stresses can be described by continuous time functions, provided the loading functions are not discontinuous. Then for multi-D loadings and general anisotropic linear viscoelastic media, the Class III PR is given by
Consequently, in all these instances the Class III PRs are also dependent on all non-zero stresses and on their histories as well as explicitly on time. Thus, different experimental determinations of PRs involving diverse stress and strain histories, lead to distinct PR functions for the same viscoelastic medium.
The neglect of the loading cycle in the period leads to serious ina- ccuracies and discrepancies in response patterns particularly for slower loading rates as demonstrated analytically in   and confirmed experimentally in  -  . See Eqs. (43)-(56) and Figure 2, Figure 7, Figure 8 and Figure 10, which are pertinent to the above discussion.
Matching analyses for constant and for can also be per- formed and similar conclusions can be reached for each of these loading cases.
Additionally, Category III PRs and their similar cousins carry with them an implicit contradiction. Consider a 1-D loading, where this function is integrable for. See Table 1. Upon FT inversion one is lead to the limited relation
when in fact it should read for, cf. Eq. (36).
Figure 10. Time to reach constant strain rate of 50 mm/min  .
Eqs. (67) to (69) are identical if and only if the loading is 1-D and the stress is limited to, i.e. a time independent stress component, which leads to significant characterization errors   . In either exact or approximate representation Class III PRs are, therefore, also stress history (path, process) dependent.
Furthermore, if the constitutive relations are of the linear non-convolution type then for the 1-D loading
and this PR is clearly implicitly stress and stress history (path) dependent, since generally. The non-convolution phenomenon takes place under the frequent time dependent temperature conditions when and then
In summary, it is seen that
1) Eq. (63) is needed for the constructions of an EVCP in the transform space and its 1-D form is stress independent
2) Eq. (69) is the proper form of the constitutive relations in the real time space and is stress and stress history dependent
3) Proper care must be exercised in using and interpreting Eqs. (63) and (69) since they represent contradictory representations
Comment 8. While the preceding analyses establish Class III PR protocols that lead to the EVCP, they provide no clues as to what the proper accom- panying PR functions are. The PR analysis is built on the superposition of responses to 1-D loading(s) without specification as what it/they should be, i.e. constant stress, strain, their time derivatives, etc. It has been established that viscoelastic PRs, unlike moduli and compliances, are time, stress and stress history (path) dependent and hence non-unique. This ambiguity renders viscoelastic PRs unsuited for general material characterization. Consequently, constitutive relations and EVCPs based on universal path independent moduli and/or compliances remain the material descriptors of choice.
In Figure 1, the Class I PRs for the same material are presented for four distinct experimental loading conditions (2 relaxation and 2 creep experiments) but with identical specimen, types of experimental loadings and boundary conditions  . Calculated from Eq. (1) and based on actual measured strains, the resulting PRs for these four different 1-D loadings exhibit time variations that are significantly distinct to testify to the viscoelastic PR inability to be classified as a unique material property due to its strong dependence on loading histories even in relatively small strain and short time ranges. In the two relaxation experiments and for over 95% of the operating time intervals, the viscoelastic PR values exceed 0.5, which is the positive upper limit for linear elasticity. In the creep experiment for the same material the PR values are all less than a half. The only difference is the loading rate and the PRs were computed directly from the measured strains, without regard to any a priori material property determina- tions. Similarly, the same experimental data from the four distinct loadings will yield four distinct curves of Class III PRs.
Similar differences in PR time function obtained in creep and relaxation were reported in  .
Attempts at determining the time time when steady-state loadings are reached based on examination of strains as seen in Figure 9 have been disappointing and misleading  . This protocol produced physically indefensible results with erroneous trends  . On the other hand a new protocol based on strain accelerations calculated from measured strains produce physical results Figure 11 and Figure 12. However, the determination of moduli and/or compliances without PRs avoided any problems by num- erically integrating the governing over the entire time range while specifically including the actual loading phases and measured strains    .
Comment 9. The great universal utility of elastic Poisson’s ratios due to their characteristically constant values for each homogeneous elastic mate- rial at a given temperature is lost in viscoelasticity. Conclusive experi- mental and analytical evidence undisputedly indicates that viscoelastic PRs are time, stress and stress history dependent and, therefore, not unique material descriptors. Consequently, they are not exportable or interchangeable from one loading condition to another.
Furthermore, the analysis demonstrates that the elastic-viscoelastic corres- pondence principle can only be expressed in terms of relaxation moduli or creep compliances and does not involve PRs. The equivalent triple product of Fourier transforms of PRs, moduli/compliances and strains/stresses necessary for the establishment of the EVCP, while mimicking Hooke's law, does not invert to a proper set of constitutive relations in the time space and, therefore, no EVCP involving PRs is possible.
Additionally and separately, analyses and experiments show that starting loading build ups, including achievable testing machine cross head speeds, are of significant importance to mandate their inclusion in material characterization protocols. Their neglect in the definitions of Class II and VI PRs, while mathematically defensible, are physically unattainable and unrealistic. Consequently, Class II and VI definitions only exist as idealized mathematical models but
Figure 11. Dynamic experimental impulse force/acceleration  .
Figure 12. Time to reach steady-state strain rates based on acceleration measurements  .
cannot be duplicated physically for as evidenced by the typical experimental results shown in Figure 10. On the other hand, if the cross head accelerations are used then reasonable times to reach steady-state strain rates can be calculated, as seen in Figure 12. Consequently, the latter should be the protocol of choice for the determination of.
Even if these two mathematical models were physically reproducible, the PRs so determined would only be applicable to the particular loading conditions under which they were experimentally determined. Consequently, the inherently specialized Class II and VI behaviors cannot be generalized and/or extended to any other loading processes since they produce process specific PRs as evidenced by Figure 1. For that matter, nor can any of the viscoelastic PR Classes I through VI be generalized as they are process specific and not universal material descri- ptors such as relaxation moduli and creep compliances or relaxation and creep functions.
The pervasive advantages and disadvantages of viscoelastic PRs vis-à-vis the EVCP and their lack of universality are summarized in Table 4.
For linear elastic and viscoelastic media:
1) Elastic PRs are of fundamental importance as linear material descriptors because they are constants for each material without regard as to what experi- mental protocol is used to determine their value. However, such is not the case of their linear viscoelastic counterparts as they are time, stress and stress history (path, process) dependent. Hence they have distinct values and functionality for each and every process. Consequently, experimentally obtained viscoelastic PRs do not match more complicated service loading conditions.
2) Elastic and viscoelastic PRs are inherently nonlinear, cf. (11) and (12). Their linear use is based on the explicit approximation of their linearity as governed by the ratio.
3) All viscoelastic PRs are time dependent. Additionally, by definition PRs of all Classes, except VII, are time, stress and stress history (path, process) depen- dent and, therefore, non-unique, i.e. dependent on loading history. All consti- tutive relations involving PRs written in the time space are very computationally intensive compared to expressions involving only moduli or compliances.
4) The principal difference between PRs and moduli/compliances is that the former in all Classes except VII are an algebraic ratio of mutually perpendicular tension and/or compression strains, whereas the latter are ratios of stresses to their responding strains. This renders the PRs nonlinear functions of the two strains.
5) Unlike elastic PRs and whether in a homogeneous or non-homogeneous medium, all of the six Classes of viscoelastic Poisson’s ratios are represented by functions of time, stresses and stress histories. Due to their process dependency they are non-unique material descriptors and as such unsuitable for general material characterizations. Therefore, characterizations in terms of relaxation moduli and/or creep compliances are the characterization vehicles of choice as they represent universal functions with parameters specific to each material and while
Table 4. Summary of viscoelastic PRs and moduli advantages and disadvantages.
Note: All PRs are time dependent. Additionally, by definition PRs of all Classes, except VII, are time, stress and stress history (path, process) dependent and, therefore, non-unique. All constitutive relations involving PRs written in the time space are very computationally intensive compared to expressions involving only moduli or compliances.
time dependent they are process independent.
6) The indisputable proof of the explicit time and path dependence of visco- elastic PRs can be found in experimental results, such as those of Figure 1, where PRs due to different 1-D “constant” loading conditions are displayed. The various PR functions obtained from creep, relaxation and prescribed strain rate experiments bear no relations to each other and are specific and distinct expre- ssions for each loading event.
7) Instantaneous non-zero loads, although mathematically rigorously defin- able by Dirac delta and/or Heaviside step functions, are physically unattainable. Consequently, the mathematically acceptable Class II and VI PR models are physically nonexistent and experimentally non-producible since they are solely based on one time independent stress or strain tensor.
8) Except in wave propagation problems, elastic loading histories are generally unimportant, since these are conservative materials without dissipation and memory. Viscoelastic materials, on the other hand, continuously dissipate ener- gy and loading histories are of prime significance to their responses.
9) When the proper loading cycles that are physically necessary to achieve constant values of one or more state variables are included in the analyses, sub- stantially different viscoelastic responses are achieved when compared to instan- taneous loading cycles.
10) Analyses establish that Class III PRs are capable to produce EVCP pro- tocols, on condition that physically acceptable double convolution integral constitutive relations are established (32b). However, these analyses fail to provide directions for establishing the unique, proper and necessary viscoelastic Class III PR functions of time, stress and stress history. Consequently, what process is to define the PRs “constant” strain, stress, their time derivatives, etc.? Hence, this inherent ambiguity renders viscoelastic path dependent non-unique PRs unsuited for general viscoelastic material characterization.
11) If viscoelastic PRs are to be used, then they must be obtained experi- mentally by loading and time processes that are exact detailed duplicates of the desired field conditions.
12) The Class VII is an artificially defined isothermal viscoelastic PR devoid of strain or stress associations and hence path independent while still remaining a time dependent universal material property descriptor. It obeys the construction mandated by Eqs. (32) necessary for the existence of EVCPs and does not suffer from the loading cycle problems associated with the other Classes. (See the next Item below.) However, its use in the time space is subject to intensive CPU use involving triple function products and double convolution integrals.
13) Since the viscoelastic PRs are time, stress and stress history dependent, unless the experimentally determined history exactly matches the application loading path no correspondence in properties can be established.
14) Some of the triple Fourier transform products, Eq. (32) necessary for establishing any PR based elastic viscoelastic correspondence principle have no physical real time counterpart constitutive relations, thus proving their nonexi- stence of the EVCP based on these PRs.
15) The evaluation of all the PR double convolution time integral constitutive relations stemming from the triple FT products of (32) are highly CPU intensive, thus impractical for use in analytical and/or finite element analyses. By contrast, constitutive relations without PRs but in terms of moduli or compliances lead to more efficient single convolution double product time integrals.
16) The validity of EVCPs without PRs based solely on moduli/compliances or on relaxation/creep functions remains unaffected and should be considered the universal unambiguous material characterization of choice.
17) If viscoelastic PRs are to be used in routine stress/strain analyses then their experimental determinations must exactly duplicate each and every actual service time and loading history.
18) In any viscoelastic formulation there is no need for the inclusion of PRs as all constitutive relation material properties (isotropic and anisotropic, homo- geneous and nonhomogeneous, linear and nonlinear) can be unambiguously characterized by their moduli and/or compliances thus removing any contro- versial descriptions due to the viscoelastic PR’s time and loading history dependence.
19) Finally, elastic solutions that are heavily dependent on PRs, such as thin plates, shells  , Timoshenko beams  , are best recast in their viscoelastic reincarnations solely in terms of relaxation moduli or creep compliances.
*Professor Emeritus of Aerospace Engineering and Senior Academic Lead for Computational Structural/Solid, Mechanics at NCSA. Fellow of the American Institute of Aeronautics and Astronautics (AIAA) and Fellow of the American Society for Composites (ASC).
2For example in a 1-D stress field with, the strain is the one in the loaded direction and is the transverse strain component.
3See the Section “The physics of time invariant stresses, strains and their derivatives” on the physical aspects and implications of time invariant stresses, strains or their time derivatives, particularly during actual loading phases.
4Also see the section “The physics of time invariant stresses, etc.”
5The EVCP is inapplicable if the material property are considered temperature dependent and if the temperature is a time function, i.e..
6cf. Eq. (11).
7This is but one of the two triple products in the constitutive relations. The other one involves and and it experiences an identical formulation in the constitutive relations for volume chan- ges.
8Eqs. (36) apply equally to linear homogeneous viscoelastic media.
9For non-auxetic elastic materials, and while. For auxetic materials the reverse inequalities hold and the elastic PR’s are in the range. However, viscoelastic PRs have no such restrictions except at, since the ICs are elastic conditions.