Nonlinear Constitutive Equations

Herrman, W. and Nunziato, J.W., Nonlinear Constitutive Equations, in Dynamic Response of Materials to Intense Impulsive Loading (edited by Chou, P.C. and Hopkins, A.K.), Air Force Materials Laboratory, Wright-Paterson AFB, 1973, Chap. 5. 

It has been shown that the thermodynamic foundations of plasticity may be considered within the framework of the continuum mechanics of materials with memory. A nonlinear material with memory is defined by a system of constitutive equations in which some state functions such as the stress tension or the internal energy, the heat flux, etc., are determined as functionals of a function which represents the time history of the local configuration of a material particle. 

Control of nutrient transport dictates significant coupling between transported components in G1 epithelia. This complicates solute transport analysis by requiring a multicomponent description. Flux equations written for each component constitute a nonlinear system in which the coupling nonlinearities are embodied in the coefficients modifying individual transport contributions to flux. 

Often times concentrated polymeric solutions cannot be treated as Newtonian fluids, however, and this tends to offset the simplifications which result from the creeping flow approximation and the fact that the boundaries are well defined. The complex rheological behavior of polymeric solutions and melts requires that nonlinear constitutive equations, such as Eqs. (l)-(5), be used (White and Metzner, 1963) ... 

We recall that our wave equation represents a long wave approximation to the behavior of a structured media (atomic lattice, periodically layered composite, bar of finite thickness), and does not contain information about the processes at small scales which are effectively homogenized out. When the model at the microlevel is nonlinear, one expects essential interaction between different scales which in turn complicates any universal homogenization procedure. In this case, the macro model is often formulated on the basis of some phenomenological constitutive hypotheses nonlinear elasticity with nonconvex energy is a theory of this type. 

The analysis of polymer processing is reduced to the balance equations, mass or continuity, energy, momentum and species and to some constitutive equations such as viscosity models, thermal conductivity models, etc. Our main interest is to solve this coupled nonlinear system of equations as accurately as possible with the least amount of computational effort. In order to do this, we simplify the geometry, we apply boundary and initial conditions, we make some physical simplifications and finally we chose an appropriate constitutive equations for the problem. At the end, we will arrive at a mathematical formulation for the problem represented by a certain function, say / (x, T, p, u,...), valid for a domain V. Due to the fact that it is impossible to obtain an exact solution over the entire domain, we must introduce discretization, for example, a grid. The grid is just a domain partition, such as points for finite difference methods, or elements for finite elements. Independent of whether the domain is divided into elements or points, the solution of the problem is always reduced to a discreet solution of the problem variables at the points or nodal pointsinxxnodes. The choice of grid, i.e., type of element, number of points or nodes, directly affects the solution of the problem. 

Fiber segments are nonlinear, exhibiting an exponential constitutive equation that makes them much stiffer in tension than in compression. 

The nonlinear viscoelastic models (VE), which utilize continuum mechanics arguments to cast constitutive equations in coordinate frame-invariant form, thus enabling them to describe all flows steady and dynamic shear as well as extensional. The objective of the polymer scientists researching these nonlinear VE empirical models is to develop constitutive equations that predict all the observed rheological phenomena. 

Together with Eq. 3.3-17, Eq. 3.3-16 is the White-Metzner constitutive equation, which has been used frequently as a nonlinear viscoelastic model. Of course, for small deformations, X(i) = dx/dt, and the single Maxwell fluid equation (Eq. 3.3-9) is obtained. 

Molecular theories, utilizing physically reasonable but approximate molecular models, can be used to specify the stress tensor expressions in nonlinear viscoelastic constitutive equations for polymer melts. These theories, called kinetic theories of polymers, are, of course, much more complex than, say, the kinetic theory of gases. Nevertheless, like the latter, they simplify the complicated physical realities of the substances involved, and we use approximate cartoon representations of macromolecular dynamics to describe the real response of these substances. Because of the relative simplicity of the models, a number of response parameters have to be chosen by trial and error to represent the real response. Unfortunately, such parameters are material specific, and we are unable to predict or specify from them the specific values of the corresponding parameters of other... 

The constitutive equations discussed previously contain both linear and nonlinear response parameters. Both have to be evaluated experimentally. The first five to ten terms... 

Such contours are then compared with numerically derived ones obtained in the following fashion a nonlinear response constitutive equation is selected (60,61,64,69) to be used with the equation of motion for a given PIF. The numerical solution computational... 

In the example given, the constitutive equation used is a multimode Phan Tien Tanner (PTT). It requires the evaluation of both linear and nonlinear material-response parameters. The linear parameters are a sufficient number of the discrete relaxation spectrum 2, and r]i pairs, which are evaluated from small-strain dynamic experiments. The values of the two nonlinear material-response parameters are evaluated as follows. Three semiarbitrary initial values of the two nonlinear parameters are chosen and the principal normal stress difference field is calculated for each of them using the equation of motion and the multimode PTT. They are compared at each field point (i, j) to the experimentally obtained normal stress difference and used in the following cost function F... 

Thus, adequate determination of nonlinear rheological parameters can be obtained, using industrial polymer processing-relevant flows, albeit with very substantial computational efforts, virtually assuring the relevance of the use of the constitutive equation for solving other complex processing flows. 

Petrie and Ito (84) used numerical methods to analyze the dynamic deformation of axisymmetric cylindrical HDPE parisons and estimate final thickness. One of the early and important contributions to parison inflation simulation came from DeLorenzi et al. (85-89), who studied thermoforming and isothermal and nonisothermal parison inflation with both two- and three-dimensional formulation, using FEM with a hyperelastic, solidlike constitutive model. Hyperelastic constitutive models (i.e., models that account for the strains that go beyond the linear elastic into the nonlinear elastic region) were also used, among others, by Charrier (90) and by Marckmann et al. (91), who developed a three-dimensional dynamic FEM procedure using a nonlinear hyperelastic Mooney-Rivlin membrane, and who also used a viscoelastic model (92). However, as was pointed out by Laroche et al. (93), hyperelastic constitutive equations do not allow for time dependence and strain-rate dependence. Thus, their assumption of quasi-static equilibrium during parison inflation, and overpredicts stresses because they cannot account for stress relaxation furthermore, the solutions are prone to numerical instabilities. Hyperelastic models like viscoplastic models do allow for strain hardening, however, which is a very important element of the actual inflation process. 

The three constant Oldroyd model is a nonlinear constitutive equation of the differential corrotational type, such as the Zaremba-Fromm-Dewitt (ZFD) fluid (Eq. 3.3-11). [For details, see R. B. Bird, R. C. Armstrong, and O. Hassager, Dynamics of Polymeric Liquids, Second Edition, Vol. 1, Wiley, New York, 1987, Table 7.3-2.]... 

In the first place, the averaged model equations are highly nonlinear and require sophisticated numerical analysis for solution. For example, the attempt to obtain numerical solutions for motions of polymeric liquids, based upon simple continuum, constitutive equations, is still not entirely successful after more than 10 years of intensive effort by a number of research groups worldwide [27]. It is possible, and one may certainly hope, that model equations derived from a sound description of the underlying microscale physics will behave better mathematically and be easier to solve, but one should not underestimate the difficulty of obtaining numerical solutions in the absence of a clear qualitative understanding of the behavior of the materials. 

This introduces SI, a minor variant of SI, which is conventionally (and confusingly) called MKS and is used in nonlinear optics. Equation (2.7.12), linking D to E, is the "first constitutive equation." The magnetic case is similar The magnetic induction B is the appropriately scaled sum of the magnetic field El and the magnetization M ... 

The bond graph method defines the structure and constitutive equations of the system. Standard bond graph elements are used to build a model of the structure of the system. Suitable computer programs are available to generate the governing equations, and alternative methods have also been developed for deriving equivalent block diagrams, which can represent nonlinear systems. 

There are two general types of constitutive equations for fluids Newtonian and non-Newtonian. For Newtonian fluids, the relation between the stress tensor, t, and the rate of deformation tensor or the shear stress is linear. For non-Newtonian fluids the relation between the stress tensor and the rate of deformation tensor is nonlinear. The various Newtonian and non-Newtonian rheologies of fluids are shown in Figure 12.2. There are four types of behavior (1) Newtonian, (2) pseudo-plastic, (3) Bingham plastic, and (4) dilatent. The reasons for these different rheological behaviors will also be discussed in subsequent sections of this chapter. But first it is necessary to relate the stress tensor to the rate of deformation tensor. 

One may use the linear viscoelastic data as a pure rheological characterization, and relate the viscoelastic parameters to some processing or final properties of the material inder study. Furthermore, linear viscoelasticity and nonlinear viscoelasticity are not different fields that would be disconnected in most cases, a linear viscoelastic function (relaxation fimction, memory function or distribution of relaxation times) is used as the kernel of non linear constitutive equations, either of the differential or integral form. That means that if we could define a general nonlinear constitutive equation that would work for all flexible chains, the knowledge of a single linear viscoelastic function would lead to all rheological properties. 

From a numerical viev point, rapid progress has been made in the last few years in studies generally devoted to the entry flow problem, together with the use of more and more realistic constitutive equations for the fluids. Consequently, more complexity was involved for the munerical problem, in relation to the nonlinearity induced by the rheological model in the governing equations. The use of nonlinear constitutive models required approximate methods for solving the equations, such as finite element techniques, even for isothermal and steady-state conditions related to a simple flow geometry. 

Equations (3-32)-(3-34) are equivalent to the so-called Oldroyd-B equation. The Oldroyd-B equation is a simple, but qualitatively useful, constitutive equation for dilute solutions of macromolecules (see Section 3.6.2). Refinements to the simple elastic dumbbell model, such as the effects of the nonlinearity of the force-extension relationship at high extensions, are discussed in Section 3.6.2.2.I. 

As remarked earlier, the nonlinear viscoelastic behavior of entangled wormy micellar solutions is similar to that of entangled flexible polymer molecules. Cates and coworkers (Cates 1990 Spenley et al. 1993, 1996) derived a full constitutive equation for entangled wormy micellar solutions, based on suitably modified reptation ideas. The stress tensor obtained from this theory is (Spenley et al. 1993)... 

A mathematical expression relating forces and deformation motions in a material is known as a constitutive equation. However, the establishment of constitutive equations can be a rather difficult task in most cases. For example, the dependence of both the viscosity and the memory effects of polymer melts and concentrated solutions on the shear rate renders it difficult to establish constitute equations, even in the cases of simple geometries. A rigorous treatment of the flow of these materials requires the use of fluid mechanics theories related to the nonlinear behavior of complex materials. However, in this chapter we aim only to emphasize important qualitative aspects of the flow of polymer melts and solutions that, conventionally interpreted, may explain the nonlinear behavior of polymers for some types of flows. Numerous books are available in which the reader will find rigorous approaches, and the corresponding references, to the subject matter discussed here (1-16). 

Many constitutive equations have been proposed in addition to those indicated above, which are special cases of fluids with memory. Most of these expressions arise from the generalization of linear viscoelasticity equations to nonlinear processes whenever they obey the material indifference principle. However, these generalizations are not unique, because there are many equations that reduce to the same linear equation. It should be noted that a determined choice among the possible generalizations may be suitable for certain types of fluids or special kinds of deformations. In any case, the use of relatively simple expressions is justified by the fact that they can predict, at least qualitatively, the behavior of complex fluids. 

---