On instabilities of single-integral constitutive equations for viscoelastic liquids

Various types of instabilities are exposed in this paper for time-strain separable single-integral viscoelastic constitutive equations (CE's). They were distinguished into two groups and defined as Hadamard and dissipative type of instabilities. As for the Hadamard-type, previously obtained criteria are found to be necessary only. They are necessary and sufficient only for thermodynamic stability. Improved, stricter Hadamard stability criteria are described briefly in this paper, and then applied to study of stability of several CE's. It is shown that the Currie potential with the K-BKZ equation and the model proposed by Papanastasiou et al. are Hadamard unstable. In the case of dissipative stability, the necessary and sufficient condition for stress boundedness in any regular flow with a given history, is proved. Then, this criterion was applied to the neoHookean, Mooney, and Yen and McIntire specifications of the general K-BKZ model, to exhibit unbounded solutions. In addition, Larson-Monroe potential which is later proved to be Hadamard unstable but satisfies the above criterion of boundedness, is shown to have unstable decreasing branch in steady simple shear flow. At present, to the authors' knowledge, there is no viscoelastic single-integral CE of factorable type proposed in the literature which can satisfy all the Hadamard and dissipative stability criteria.     

黏塑性本构计算的稳定性分析

提出本构方程计算方法的稳定性问题,针对黏塑性本构计算的显式精确算法的稳定性进行分析,发现该算法并非无条件稳定,使用小扰动方法给出了其计算稳定的必要条件,稳定性条件对数值计算中的时间步长提出限制要求。通过有限元算例验证了分析的正确性,计算结果也表明理论推导得到的稳定性公式能够准确预测满足计算稳定性条件要求的最大时间步长与各参数之间关系。     

Mathematical characteristics of the pom-pom model

 Recently, in order to describe the complex rheological behavior of polymer melts with long side branches like low density polyethylene, new constitutive equations called the pom-pom equations have been derived in the integral/differential form and also in the simplified differential type by McLeish and Larson on the basis of the reptation dynamics with simplified branch structure taken into account. In this study, mathematical stability analysis under short and high frequency wave disturbances has been performed for these constitutive equations. It is proved that the differential model is globally Hadamard stable, as long as the orientation tensor remains positive definite or the smooth strain history in the flow is previously given. However both versions of the model are Hadamard unstable if we neglect the arm withdrawal in the case of maximum backbone stretch. It is also dissipatively unstable, since the steady shear flow curves exhibit non-monotonic dependence on shear rate. Additionally, in the flow regime of creep shear flow where the applied constant shear stress exceeds the maximum achievable value in the steady flow curves, the constitutive equations exhibit severe instability that the solution possesses strong discontinuity at the moment of change of chain dynamics mechanisms. Received: 14 August 2001 Accepted: 18 October 2001     

Torsional flow stability of highly dilute polymer solutions

We examine the torsional flow stability, to axisymmetric disturbances, of a variety of multimode and non-linear constitutive models in a bounded parallel plate geometry. The analysis is facilitated by the construction of a regular perturbation scheme in the ratio of polymer to total viscosity. As a model for Boger fluids this corresponds to the assumption that the Boger fluid is highly dilute. The consequent simplification provided by the perturbation scheme allows us to examine the effects of a discrete spectrum of relaxation times, shear thinning, second normal stress difference, and finite extensional viscosity on the torsional instability.     

On non-linear stability in unconstrained non-linear elasticity

A method of obtaining a full three-dimensional non-linear Hadamard stability analysis of inhomogeneous deformations of arbitrary, unconstrained, hyperelastic materials is presented. The analysis is an extension of that given by Chen and Haughton (Proc. Roy. Soc. London A 459 (2003) 137) for two-dimensional incompressible problems. The process that we present replaces the second variation condition expressed as an integral involving a quadratic in three arbitrary perturbations, with an equivalent sixth-order system of ordinary differential equations. The positive definiteness condition is thereby reduced to the simple numerical evaluation of zeros of a well-behaved function. The general theory is illustrated by applying it to the problem of the inflation of a thick-walled spherical shell. The present analysis provides a simpler alternative approach to bifurcation problems approached by using the incremental equations of non-linear elasticity.     

Analysis of diffuse plastic stability in tubes and sheets

The diffuse plastic instability in tubes and sheets under biaxial stress conditions is examined by the use of perturbation methods. Very general constitutive relationships for material properties are used. This requires the inclusion of first order changes in the strain directions inside the patch and also treatment of the material anisotropy and strain rate sensitivity in addition to strain hardening. The inclusion of variations in strain direction is found to alter the form of the characteristic equation for stability from first order to second order but both roots are real for all cases investigated. The value of strain hardening at which the largest root becomes significantly positive is almost the same as that reached when changes in strain directors are ignored. However the strain hardening at which this root becomes formally zero can be very different. The former condition is considered to be of more practical importance than the latter. By this test the stability increases rapidly above a strain rate sensitivity of about 0.1.     

A variational approach for materially stable anisotropic hyperelasticity

In this work we propose an anisotropic stored energy function which satisfies a priori the Legendre–Hadamard condition, which is strongly related to the material stability of the constitutive equations. In the linearized case this condition implies positive wave speeds. The Legendre–Hadamard condition plays also an important role for the (local) existence of solutions in the neighborhood of stationary points. We apply the proposed hyperelastic energies to soft tissues and compare the formulation with existing models which have been used for the calculation of medial collateral ligament and arterial walls. In our numerical and analytical investigations we discuss the distribution of wave speeds for a sequence of deformation states containing some essential stress–strain characteristics of the compared models.     

On secondary loops in LAOS via self-intersection of Lissajous–Bowditch curves

When the shear stress measured in large amplitude oscillatory shear (LAOS) deformation is represented as a 2-D Lissajous–Bowditch curve, the corresponding trajectory can appear to self-intersect and form secondary loops. This self-intersection is a general consequence of a strongly nonlinear material response to the imposed oscillatory forcing and can be observed for various material systems and constitutive models. We derive the mathematical criteria for the formation of secondary loops, quantify the location of the apparent intersection, and furthermore suggest a qualitative physical understanding for the associated nonlinear material behavior. We show that when secondary loops appear in the viscous projection of the stress response (the 2-D plot of stress vs. strain rate), they are best interpreted by understanding the corresponding elastic response (the 2-D projection of stress vs. strain). The analysis shows clearly that sufficiently strong elastic nonlinearity is required to observe secondary loops on the conjugate viscous projection. Such a strong elastic nonlinearity physically corresponds to a nonlinear viscoelastic shear stress overshoot in which existing stress is unloaded more quickly than new deformation is accumulated. This general understanding of secondary loops in LAOS flows can be applied to various molecular configurations and microstructures such as polymer solutions, polymer melts, soft glassy materials, and other structured fluids.     

Variational formulation and stability analysis of a three dimensional superelastic model for shape memory alloys

This paper presents a variational framework for the three-dimensional macroscopic modelling of superelastic shape memory alloys in an isothermal setting. Phase transformation is accounted through a unique second order tensorial internal variable, acting as the transformation strain. Postulating the total strain energy density as the sum of a free energy and a dissipated energy, the model depends on two material scalar functions of the norm of the transformation strain and a material scalar constant. Appropriate calibration of these material functions allows to render a wide range of constitutive behaviours including stress-softening and stress-hardening. The quasi-static evolution problem of a domain is formulated in terms of two physical principles based on the total energy of the system: a stability criterion, which selects the local minima of the total energy, and an energy balance condition, which ensures the consistency of the evolution of the total energy with respect to the external loadings. The local phase transformation laws in terms of Kuhn–Tucker relations are deduced from the first-order stability condition and the energy balance condition.The response of the model is illustrated with a numerical traction–torsion test performed on a thin-walled cylinder. Evolutions of homogeneous states are given for proportional and non-proportional loadings. Influence of the stress-hardening/softening properties on the evolution of the transformation domain is emphasized. Finally, in view of an identification process, the issue of stability of homogeneous states in a multi-dimensional setting is answered based on the study of second-order derivative of the total energy. Explicit necessary and sufficient conditions of stability are provided.     

Analysis of material instabilities in inelastic solids by incremental energy minimization and relaxation methods: evolving deformation microstructures in finite plasticity

We propose an approach to the definition and analysis of material instabilities in rate-independent standard dissipative solids at finite strains based on finite-step-sized incremental energy minimization principles. The point of departure is a recently developed constitutive minimization principle for standard dissipative materials that optimizes a generalized incremental work function with respect to the internal variables. In an incremental setting at finite time steps this variational problem defines a quasi-hyperelastic stress potential. The existence of this potential allows to be recast a typical incremental boundary-value problem of quasi-static inelasticity into a principle of minimum incremental energy for standard dissipative solids. Mathematical existence theorems for sufficiently regular minimizers then induce a definition of the material stability of the inelastic material response in terms of the sequentially weakly lower semicontinuity of the incremental variational functional. As a consequence, the incremental material stability of standard dissipative solids may be defined in terms of the quasi-convexity or the rank-one convexity of the incremental stress potential. This global definition includes the classical local Hadamard condition but is more general. Furthermore, the variational setting opens up the possibility to analyze the post-critical development of deformation microstructures in non-stable inelastic materials based on energy relaxation methods. We outline minimization principles of quasi- and rank-one convexifications of incremental non-convex stress potentials for standard dissipative solids. The general concepts are applied to the analysis of evolving deformation microstructures in single-slip plasticity. For this canonical model problem, we outline details of the constitutive variational formulation and develop numerical and semi-analytical solution methods for a first-level rank-one convexification. A set of representative numerical investigations analyze the development of deformation microstructures in the form of rank-one laminates in single slip plasticity for homogeneous macro-deformation modes as well as inhomogeneous macroscopic boundary-value problems. The well-posedness of the relaxed variational formulation is indicated by an independence of typical finite element solutions on the mesh-size.     

On the criterion for the absolute stability of the control system

In this paper, firstly we give the criterion for the absolute stability of the second canonical form for the control system, including the equation of the longitudinal motions of a plane as a particular example. The corresponding result in [8], [9] is a particular example given in this paper. Secondly, we give the criteria for the absolute stability of the first canonical form in the usual case and in the critical case. Finally, we give some criteria for the absolute stability of the general form for the direct control system.All the results in this paper merely depend upon the relations between the parameters of the system itself to give an explicit algebraic discriminant.     

Overspeed burst of elastoviscoplastic rotating disks – Part I: Analytical and numerical stability analyses

Burst of rotating disks in case of overspeed is investigated. The certification of turbo-engines requires to demonstrate the integrity of disks at higher rotation speeds than the maximum rotation speed reachable in service. The determination of the burst speed by analysis can help to reduce the number of tests required for the certification. This prediction can be established by non-linear stability analyses of finite element simulations. Non-linearities originate from (i) the material behaviour described by elastoviscoplastic constitutive equations, (ii) geometric changes accounted for by the finite strain formulation. In this work, loss of uniqueness and loss of stability criteria from [Hill, R., 1958. A general theory of uniqueness and stability in elastic-plastic solids. J. Mech. Phys. Solids 6, 236–249] are applied. The loss of stability criterion is restricted to the case of rotating disks and compared to several simple widely used material based failure criteria. 3D simulations of rotating metal disks are performed for a given elastoviscoplastic behaviour and the stability criteria are evaluated. The sensitivity of disk stability to material parameters, such as yield criterion, hardening and viscosity is evaluated in the case of a nickel based superalloy.     

Scattered fracture of porous materials with brittle skeleton

A model of damage accumulation in a porous medium with a brittle skeleton saturated with a compressible fluid is formulated in the isothermal approximation. The model takes account of the skeleton elastic energy transformation into the surface energy of microcracks. In the case of arbitrary deformations of an anisotropic material, constitutive equations are obtained in a general form that is necessary and sufficient for the objectivity and thermodynamic consistency principles to be satisfied. We also formulate the kinetics equation ensuring that the scattered fracture dissipation is nonnegative for any loading history. For small deviations from the initial state, we propose an elastic potential which permits describing the principal characteristics of the behavior of a saturated porous medium with a brittle skeleton. We study the acoustic properties of the material under study and find their relationship with the strength criterion depending on the accumulated damage and the material current deformation. We consider the problem of scattered fracture of a saturated porous material in a neighborhood of a spherical cavity. We show that the cavity failure occurs if the Hadamard condition is violated.     

Practical comparison of differential viscoelastic constitutive equations in finite element analysis of planar 4:1 contraction flow

Finite element modeling of planar 4:1 contraction flow (isothermal incompressible and creeping) around a sharp entrance corner is performed for favored differential constitutive equations such as the Maxwell, Leonov, Giesekus, FENE-P, Larson, White-Metzner models and the Phan Thien-Tanner model of exponential and linear types. We have implemented the discrete elastic viscous stress splitting and streamline upwinding algorithms in the basic computational scheme in order to augment stability at high flow rate. For each constitutive model, we have obtained the upper limit of the Deborah number under which numerical convergence is guaranteed. All the computational results are analyzed according to consequences of mathematical analyses for constitutive equations from the viewpoint of stability. It is verified that in general the constitutive equations proven globally stable yield convergent numerical solutions for higher Deborah number flows. Therefore one can get solutions for relatively high Deborah number flows when the Leonov, the Phan Thien-Tanner, or the Giesekus constitutive equation is employed as the viscoelastic field equation. The close relationship of numerical convergence with mathematical stability of the model equations is also clearly demonstrated.     

Three-dimensional disturbances of planar viscometric flow

This paper examines three-dimensional disturbances of a plane steady shear flow of simple fluids with short memory. Under the assumption of nearly-viscometric flow, constitutive equations are derived and then a general form of the Reynolds-Orr energy equation is obtained. With the aid of this derived energy formula, sufficient conditions are generated for the stability of three-dimensional disturbances of the planar viscometric flow. These conditions are analysed and a comparison is made with the corresponding two-dimensional stability problem. There is a strong indication that the basic flow is less stable against three-dimensional disturbances than against two-dimensional ones.     

On the surface instability of a highly elastic half-space

The phenomenon of surface instability of an isotropic half-space under biaxial plane stress is studied for compressible elastic materials in finite strain. Euler's method is used to derive the general form of the stability criterion, and analytical details are exhibited by special application to the class of hyperelastic Hadamard materials in two complementary cases: (i) the full solution is derived for the compressible, neo-Hookean members, and (ii) the plane deformation solution is provided for every isotropic, elastic material and specific results are presented for the full Hadamard class. Results appropriate to incompressible Mooney-Rivlin materials are herein obtained as special limit cases. Several theorems are established and some of the conclusions are illustrated graphically.     

GLOBAL ASYMPTOTIC STABILITY CONDITIONS OF DELAYED NEURAL NETWORKS

1　IntroductionandProblemEductionRecently,thestudiesofthestabilityforcellularneuralnetworks (CNNs)anddelayedcellularneuralnetworks (DCNNs)haveattractedattentionsofresearchersandseveralimportantresultshavebeenobtained .MostpapersdealtwithcompletelystableCNNsandDCNNsthataresuitableforimageprocessingapplications.CNNshavebeenwidelyappliedtoimageprocessing ,toprocessmovingimages,onemustintroducedelaysinthesignalstransmittedamongthecells.Buttimedelaysmayleadtoanoscillationphenomenonand ,furt…     

Thermal softening induced plastic instability in rate-dependent materials

A linear perturbation analysis is performed for a class of rate-dependent materials, such as the Johnson-Cook model, in which the rate contribution to the stress can be separated from that of the plastic strain and temperature and in which the temperature rises adiabatically. The analysis is facilitated by perturbing both the rate of momentum equation and the momentum equation. An identical material stability/instability criterion is deduced from the characteristic spectral equations for one-dimensional deformation, one-dimensional shearing, and general three-dimensional field equations, and thus shows that the instability derived here is a material constitutive instability.The criteria indicate that the materials become unstable once the thermal softening overcomes the strain hardening, regardless of the strain rate. The strain rate enters the criteria through its effects on the accumulated temperature and the current stress. Based on the criterion, the three-dimensional instability surface is established in the space of plastic strain, plastic strain rate, and temperature. Instability surface is shown as a material property and independent of deformation histories or modes. Both necking and shear banding are simulated to validate the excellent predictive capability of the criterion.     

Constraints, Constitutive Limits, and Instability in Finite Thermoelasticity

This paper examines all the possible types of thermomechanical constraints in finite-deformational elasticity. By a thermomechanical constraint we mean a functional relationship between a mechanical variable, either the deformation gradient or the stress, and a thermal variable, temperature, entropy or one of the energy potentials; internal energy, Helmholtz free energy, Gibbs free energy or enthalpy. It is shown that for the temperature-deformation, entropy-stress, enthalpy-deformation, and Helmholtz free energy-stress constraints equilibrium states are unstable, in the sense that certain perturbations of the equilibrium state grow exponentially. By considering the constrained materials as constitutive limits of unconstrained materials, it is shown that the instability is associated with the violation of the Legendre–Hadamard condition on the internal energy. The entropy-deformation, temperature-stress, internal energy-stress, and Gibbs free energy-deformation constraints do not exhibit this instability. It is proposed that stability of the rest state (or equivalently convexity of internal energy) is a necessary requirement for a model to be physically valid, and hence entropy-deformation, temperature-stress, internal energy-stress, and Gibbs free energy-deformation constraints are physical, whereas temperature-deformation constraints (including the customary notion of thermal expansion that density is a function of temperature only), entropy-stress constraints, enthalpy-deformation constraints, and Helmholtz free energy-stress constraints are not.