Assessment of nonlinear strain measures for extensional and shearing flows of polymer melts

From stress-strain experiments in extensional and shearing flows, nonlinear strain measures and effective damping functions are derived for a polyisobutylene melt. The strain measures determined in planar extensional flow and in simple shear flow coincide. Experimental results are compared with predictions of two molecular theories, the Doi-Edwards model and the molecular stress function approach of Wagner and Schaeffer. Discrepancies between theories and experiment lead to a reconsideration of the classification of extensional flows. The symmetry of the flow field is identified and quantified as an important parameter influencing the strain measure, and a unifying strain measure for general extensional and shearing flows of polymer melts is presented.     

The strain-hardening behaviour of linear and long-chain-branched polyolefin melts in extensional flows

By generalizing the Doi-Edwards model to the Molecular Stress Function theory of Wagner and Schaeffer, the extensional viscosities of polyolefin melts in uniaxial, equibiaxial and planar constant strain-rate experiments starting from the isotropic state can be described quantitatively. While the strain hardening of four linear polymer melts (two high-density polyethylenes, a polystyrene and a polypropylene) can be accounted for by a tube diameter that decreases affinely with the average stretch, the two long-chain-branched polymer melts considered (a low-density polyethylene and a long-chain branched polypropylene) show enhanced strain hardening in extensional flows due to the presence of long-chain branches. This can be quantified by a molecular stress function, the square of which is quadratic in the average stretch and which follows from the junction fluctuation theory of Flory. The ultimate magnitude of the strain-hardening effect is governed by a maximum value of the molecular stress, which is specific to the polymer melt considered and which is the only free non-linear parameter of the theory. Received: 1 June 1999/Accepted: 24 November 1999     

A modification of the convective constraint release mechanism in the molecular stress function model giving enhanced vortex growth

The molecular stress function model with convective constraint release (MSF with CCR) constitutive model [M.H. Wagner, P. Rubio, H. Bastian, The molecular stress function model for polydisperse polymer melts with dissipative convective constraint release, J. Rheol. 45 (2001) 1387] is capable of fitting all viscometric data for IUPAC LDPE, with only two adjustable parameters (with difference found only on reported “steady-state” elongational viscosities). The full MSF with CCR model is implemented in a backwards particle-tracking implementation, using an adaptive method for the computation of relative stretch that reduces simulation time many-fold, with insignificant loss of accuracy. The model is shown to give improved results over earlier versions of the MSF (without CCR) when compared to well-known experimental data from White and Kondo [J.L. White, A. Kondo, Flow patterns in polyethylene and polystyrene melts during extrusion through a die entry region: measurement and interpretation, J. Non-Newtonian Fluid Mech. 3 (1977) 41]; but still to under-predict contraction flow opening angles. The discrepancy is traced to the interaction between the rotational dissipative function and the large stretch levels caused by the contraction flow. A modified combination of dissipative functions in the constraint release mechanism is proposed, which aims to reduce this interaction to allow greater strain hardening in a mixed flow. The modified constraint release mechanism is shown to fit viscometric rheological data equally well, but to give opening angles in the complex contraction flow that are much closer to the experimental data from White and Kondo. It is shown (we believe for the first time) that a constitutive model demonstrates an accurate fit to all planar elongational, uniaxial elongational and shear viscometric data, with a simultaneous agreement with this well-known experimental opening angle data. The sensitivity of results to inaccuracies caused by representing the components of the deformation gradient tensor to finite precision is examined; results are found to be insensitive to even large reductions in the precision used for the representation of components. It is shown that two models that give identical response in elongational flow, and a very similar fit to available shear data, give significantly different results in flows containing a mix of deformation modes. The implication for constitutive models is that evaluation against mixed deformation mode flow data is desirable in addition to evaluation against viscometric measurements.     

Constitutive modeling of temperature and strain rate dependent elastoplastic hardening materials using a corotational rate associated with the plastic deformation

In this paper, a constitutive model with a temperature and strain rate dependent flow stress (Bergstrom hardening rule) and modified Armstrong-Frederick kinematic evolution equation for elastoplastic hardening materials is introduced. Based on the multiplicative decomposition of the deformation gradient,new kinematic relations for the elastic and plastic left stretch tensors as well as the plastic deformation-dependent spin tensor are proposed. Also, a closed-form solution has been obtained for the elastic and plastic left stretch tensors for the simple shear problem.To evaluate model validity, results are compared with known experimental data for SUS 304 stainless steel, which shows a good agreement with the results of the proposed theoretical model.Finally, the stress-deformation curve, as predicted by the model, is plotted for the simple shear problem at room and elevated temperatures using the same material properties for AA5754-O aluminium alloy.     

Flows of constant stretch history for polymeric materials with power-law distributions of relaxation times

It is herein shown that for separable integral constitutive equations with power-law distributions of relaxation times, the streamlines in creeping flow are independent of flow rate.For planar flows of constant stretch history, the stress tensor is the sum of three terms, one proportional to the rate-of-deformation tensor, one to the square of this tensor, and the other to the Jaumann derivative of the rate-of-deformation tensor. The three tensors are the same as occur in the Criminale-Ericksen-Filbey Equation, but the coefficients of these tensors depend not only on the second invariant of the strain rate, but also on another invariant which is a measure of flow strength. With the power-law distribution of relaxation times, each coefficient is equal to the second invariant of the strain rate tensor raised to a power, times a function that depends only on strength of the flow. Axisymmetric flows of constant stretch history are more complicated than the planar flows, because three instead of two nonzero normal components appear in the velocity gradient tensor. For homogeneous axisymmetric flows of constant stretch history, the stress tensor is given by the sum of the same three terms. The coefficients of these terms again depend on the flow strength parameter, but in general the dependences are not the same as in planar flow.     

The nonlinear strain measure of polyisobutylene melt in general biaxial flow and its comparison to the Doi-Edwards model

The nonlinear strain measure of a polyisobutylene (PIB) melt as determined by analysis of uniaxial, planar, ellipsoidal, and equibiaxial extensions is compared to the predictions of the molecular model of Doi and Edwards. It is found that the universal strain function of the Doi-Edwards model is unable to predict the nonlinear behavior of this polymer melt in general extensional flow. The qualitative agreement between predictions and experimental data for the strain dependence of shear stress and first normal stress difference in shear flow that was considered as powerful evidence for the correctness of the Doi-Edwards model seems to be accidental. The exaggerated strain dependence of the model suggests a need to reconsider the assumptions concerning the chain retraction process.Presented at the Golden Jubilee Conference of the British Society of Rheology and Third European Rheology Conference, Edinburgh, 3–7 September, 1990.Dedicated to Professor F.R. Schwarzl on the occasion of his 65th birthday     

A continuum model for intermittent deformation of single crystal micropillars

Strain bursts are often observed during compression tests of single crystal micropillars. In this work, we formulate a new continuum model that accounts for the strain bursts within the framework of crystal plasticity. The strain bursts are separated from the loading stage (nearly elastic loading) by introducing a dimensionless constant in the continuum model, and are detected by load serrations. The boundary conditions in the context of micropillar compression are studied and they are shown to be changing and unpredictable as plastic deformation proceeds. To evaluate the validity of our model, finite element simulations of the uniaxial compression tests on nickel micropillars are performed. Our simulations produce clearly visible strain bursts during the plastic flow and the produced intermittent flows are comparable with the experimental observations. For the bulk crystal, a series of strain bursts is identified in the course of plastic flow, despite an apparently smooth stress–strain response. We also show that the intermittent flow is intensified in the micrometer-scale due to both increasing numbers of the successive strain bursts and increasing amplitude of the strain burst, when the specimen size decreases. Finally, we show that the occurrences of the strain bursts are always associated with negative values of the second-order work.     

Corotational rates in constitutive modeling of elastic-plastic deformation

The principal axes technique is used to develop a new hypoelastic constitutive model for an isotropic elastic solid in finite deformation. The new model is shown to produce solutions that are independent of the choice of objective stress rate. In addition, the new model is found to be equivalent to the isotropic finite elastic model; this is essential if both models describe the same material.

The new hypoelastic model is combined with an isotropic flow rule to form an elastic-plastic rate constitutive equation. Use of the principal axes technique ensures that the stress tensor is coaxial with the elastic stretch tensor and that solutions do not depend on the choice of objective stress rate. The flow rule of von Mises and a parabolic hardening law are used to provide an example of application of the new theory. A solution is obtained for the prescribed deformation of simple rectilinear shear of an isotropic elastic and isotropic elastic-plastic material.     

A derivation of the Mehrabadi-Cowin equations

An alternative derivation to that given by Mehrabadi and Cowin (1978) is presented here for a pair of kinematic equations governing a certain class of flows in the plastic deformation of dilatant granular materials. This class has been described by Spencer (1981) as double shearing flows. In their derivation Mehrabadi and Cowin (1978), prior to presenting the equations relative to rectangular Cartesian coordinates, obtained an intermediate pair of equations relative to a non-orthogonal network of characteristic coordinates. The essential difference between the original and present derivation is that here, the flow rule, expressed relative to rotating, rectangular Cartesian coordinates, is transformed directly to obtain the kinematic equations relative to fixed rectangular Cartesian coordinate axes, without the need to obtain the characteristic equations.     

A generalized anisotropic hardening rule based on the Mroz multi-yield-surface model for pressure insensitive and sensitive materials

In this paper, a generalized anisotropic hardening rule based on the Mroz multi-yield-surface model for pressure insensitive and sensitive materials is derived. The evolution equation for the active yield surface with reference to the memory yield surface is obtained by considering the continuous expansion of the active yield surface during the unloading/reloading process. The incremental constitutive relation based on the associated flow rule is then derived for a general yield function for pressure insensitive and sensitive materials. Detailed incremental constitutive relations for materials based on the Mises yield function, the Hill quadratic anisotropic yield function and the Drucker–Prager yield function are derived as the special cases. The closed-form solutions for one-dimensional stress–plastic strain curves are also derived and plotted for materials under cyclic loading conditions based on the three yield functions. In addition, the closed-form solutions for one-dimensional stress–plastic strain curves for materials based on the isotropic Cazacu–Barlat yield function under cyclic loading conditions are summarized and presented. For materials based on the Mises and the Hill anisotropic yield functions, the stress–plastic strain curves show closed hysteresis loops under uniaxial cyclic loading conditions and the Masing hypothesis is applicable. For materials based on the Drucker–Prager and Cazacu–Barlat yield functions, the stress–plastic strain curves do not close and show the ratcheting effect under uniaxial cyclic loading conditions. The ratcheting effect is due to different strain ranges for a given stress range for the unloading and reloading processes. With these closed-form solutions, the important effects of the yield surface geometry on the cyclic plastic behavior due to the pressure-sensitive yielding or the unsymmetric behavior in tension and compression can be shown unambiguously. The closed form solutions for the Drucker–Prager and Cazacu–Barlat yield functions with the associated flow rule also suggest that a more general anisotropic hardening theory needs to be developed to address the ratcheting effects for a given stress range.     

预延伸非晶聚合物材料各向异性变形局部化

给出一种可描述预延伸各向异性特性的背应力张量三维表达式，引入大变形弹塑性有限元驱动应力法，结合ＢＰＡ８ 链细观分子网络模型，模拟了预延伸各向异性非晶聚合物材料平面应变拉伸变形局部化力学行为．详细讨论了预延伸比（ＩｎｉｔｉａｌＤｒａｗｉｎｇＲａｔｉｏ；ＩＤＲ）和预延伸方向（ＩｎｉｔｉａｌＤｒａｗｉｎｇＤｉｒｅｃｔｉｏｎ；ＩＤＤ）对变形抗力、颈缩规律、剪切带方向以及试件中心部位链延伸比的影响．     

The \vec\omega-\vec{D} fluid: general theory with special emphasis on stationary two dimensional flows

In this paper a theory is presented in which the extra stress tensor is allowed to depend not only on the rate of strain tensor but also on the relative vorticity of the fluid, i.e. on the vorticity relative to the rate of rotation of the principal straining directions. This theory has its origin in an expansion of in terms of kinematic tensors in the limit of stationarity in a material sense (constant stretch history flows). For two dimensional flows of an incompressible fluid three tensors suffice to completely specify . The three material functions which appear can depend only on two invariants, namely the second invariant of and on . Using the predictions of an Oldroyd 8 constant fluid in a homogeneous planar flow of constant stretch history, the three material functions are studied in detail. For the special case of a quasi-Newtonian fluid shear thinning and extension thickening can directly be accounted for in the “viscosity” function. Received: September 26, 1996     

Incorporation of yield surface distortion in finite deformation constitutive modeling of rigid-plastic hardening materials based on the Hencky logarithmic strain

In this paper a constitutive model for rigid-plastic hardening materials based on the Hencky logarithmic strain tensor and its corotational rates is introduced. The distortional hardening is incorporated in the model using a distortional yield function. The flow rule of this model relates the corotational rate of the logarithmic strain to the difference of the Cauchy stress and the back stress tensors employing deformation-induced anisotropy tensor. Based on the Armstrong–Fredrick evolution equation the kinematic hardening constitutive equation of the proposed model expresses the corotational rate of the back stress tensor in terms of the same corotational rate of the logarithmic strain. Using logarithmic, Green–Naghdi and Jaumann corotational rates in the proposed constitutive model, the Cauchy and back stress tensors as well as subsequent yield surfaces are determined for rigid-plastic kinematic, isotropic and distortional hardening materials in the simple shear deformation. The ability of the model to properly represent the sign and magnitude of the normal stress in the simple shear deformation as well as the flattening of yield surface at the loading point and its orientation towards the loading direction are investigated. It is shown that among the different cases of using corotational rates and plastic deformation parameters in the constitutive equations, the results of the model based on the logarithmic rate and accumulated logarithmic strain are in good agreement with anticipated response of the simple shear deformation.     

A non‐linear k–ε model with realizability for prediction of flows around bluff bodies

The incompressible flow around bluff bodies (a square cylinder and a cube) is investigated numerically using turbulence models. A non‐linear k–ε model, which can take into account the anisotropy of turbulence with less CPU time and computer memory then RSM or LES, is adopted as a turbulence model. In tuning of the model coefficients of the non‐linear terms are adjusted through the examination of previous experimental studies in simple shear flows. For the tuning of the coefficient in the eddy viscosity (=C_μ), the realizability constraints are derived in three types of basic 2D flow patterns, namely, a simple shear flow, flow around a saddle and a focal point. C_μ is then determined as a function of the strain and rotation parameters to satisfy the realizability. The turbulence model is first applied to a 2D flow around a square cylinder and the model performance for unsteady flows is examined focussing on the period and the amplitude of the flow oscillation induced by Karman vortex shedding. The applicability of the model to 3D flows is examined through the computation of the flow around a surface‐mounted cubic obstacle. The numerical results show that the present model performs satisfactorily to reproduce complex turbulent flows around bluff bodies. Copyright © 2003 John Wiley & Sons, Ltd.     

Deformations of an elastic,internally constrained material. Part 1: Homogeneous deformations

The nonlinear elastic response of a class of materials for which the deformation is subject to an internal material constraint described in experiments by James F. Bell on the finite deformation of a variety of metals is investigated. The purely kinematical consequences of the Bell constraint are discussed, and restrictions on the full range of compatible deformations are presented in geometrical terms. Then various forms of the constitutive equation relating the stress and stretch tensors for an isotropic elastic Bell material are presented. Inequalities on the mechanical response functions are introduced. The importance of these in applications is demonstrated in several examples throughout the paper.This paper focuses on homogeneous deformations. In a simple illustration of the theory, a generalized form of Bell's empirical rule for uniaxial loading is derived, and some peculiarities in the response under all-around compressive loading are discussed. General formulae for universal relations possible in an isotropic elastic, Bell constrained material are presented. A simple method for the determination of the left stretch tensor for essentially plane problems is illustrated in the solution of the problem of pure shear of a materially uniform rectangular block. A general formula which includes the empirical rule found in pure shear experiments by Bell is derived as a special case. The whole apparatus is then applied in the solution of the general problem of a homogeneous simple shear superimposed on a uniform triaxial stretch; and the great variety of results possible in an isotropic, elastic Bell material is illustrated. The problem of the finite torsion and extension of a thin-walled cylindrical tube is investigated. The results are shown to be consistent with Bell's data for which the rigid body rotation is found to be quite small compared with the gross deformation of the tube. Several universal formulas relating various kinds of stress components to the deformation independently of the material response functions are derived, including a universal rule relating the axial force to the torque.Constitutive equations for hyperelastic Bell materials are derived. The empirical work function studied by Bell is introduced; and a new constitutive equation is derived, which we name Bell's law. On the basis of this law, we then derive exactly Bell's parabolic laws for uniaxial loading and for pure shear. Also, form Bell's law, a simple constitutive equation relating Bell's deviatoric stress tensor to his finite deviatoric strain tensor is obtained. We thereby derive Bell's invariant parabolic law relating the deviatoric stress intensity to the corresponding strain intensity; and, finally, Bell's fundamental law for the work function expressed in these terms is recovered. This rule is the foundation for all of Bell's own theoretical study of the isotropic materials cataloged in his finite strain experiments on metals, all consistent with the internal material constraint studied here.     

On the derivatives of the stretch and rotation with respect to the deformation gradient

In this paper, a result involving the eigenprojections of the right stretch and its derivative with respect to the deformation gradient is derived, and a related result is found for the rotation. As an application, the form of the constitutive law for an isotropic hyperelastic material in the case when the strain energy function is expressed in terms of the right stretch, is shown to follow at once.     

What is behind the plastic strain rate?

The plastic strain rate plays a central role in macroscopic models on elasto-viscoplasticity. In order to discuss the concept behind this quantity, we propose, first, a kinetic toy model to describe the dynamics of sliding layers representative of plastic deformation of single crystalline metals. The dynamic variable is given by the distribution function of relative strains between adjacent layers, and the plastic strain rate emerges as the average hopping rate between energy wells. We demonstrate the behavior of this model under different deformations and how it captures the elastic-to-plastic transition. Second, the kinetic toy model is reduced to a closed evolution equation for the average of the relative strain, allowing us to make a direct link to macroscopic theories. It is shown that the constitutive relation for the plastic strain rate does not only depend on the stress, but also on the macroscopic applied deformation rate, contrary to common practice.     

A cyclic plasticity model for single crystals

The Armstrong–Frederick type kinematic hardening rule was invoked to capture the Bauschinger effect of the cyclic plastic deformation of a single crystal. The yield criterion and flow rule were built on individual slip systems. Material memory was introduced to describe strain range dependent cyclic hardening. The experimental results of copper single crystals were used to evaluate the cyclic plasticity model. It was found that the model was able to accurately describe the cyclic plastic deformation and properly reflect the dislocation substructure evolution. The well-known three distinctive regimes in the cyclic stress–strain curve of the copper single crystals oriented for single slip can be reproduced by using the model. The model can predict the enhanced hardening for crystals oriented for multislip, showing the model's ability to describe anisotropic cyclic plasticity. For a given loading history, the model was able to capture not only the saturated stress–strain response but also the detailed transient stress–strain evolution. The model was used to predict the cyclic plasticity under a high–low loading sequence. Both the stress–strain responses and the microstructural evolution can be appropriately described through the slip system activation.