材料生长与变形的连续介质模型:弹性本构方程

本文着重探讨了生长变形体连续介质理论中的本构模型。首先列出了描述生长变形体能量平衡的微分方程以及熵不等式；以此为基础，通过将密度演化的历史作为独立的本构变量扩展了理性力学的因果性公理与决定性公理，具体而详细地推导了简单材料的生长弹性本构方程，给出了这些本构方程中的相关本构变量之间的约束不等式，得到了“生长变形体的自由能与其密度成反比”的结论，并从热力学上对这一结果进行了定性的解释。最后，文中对几个尚待解决的问题进行了说明，并指出了今后的发展方向。     

A pragmatic approach for deriving constitutive equations endowed with pom–pom attributes

The most important rheological and mathematical features of the pom–pom model are presently used to compare and improve other constitutive models such as the Giesekus and Phan-Thien–Tanner models. A pragmatic methodology is selected that allows derivation of simple constitutive equations, which are suited to possible software implementation. Alterations to the double convected pom–pom, Phan-Thien–Tanner and Giesekus models are proposed and assessed in rheometric flows by comparing model predictions to experimental data.

Mathematical modeling of magnetorheological fluids

Magnetorheological (MR) fluids are a class of smart materials whose rheological properties may be rapidly modified by the application of a magnetic field. These materials typically consist of micron-sized ferrous particles dispersed in a fluid. In the present paper, we consider the full system of equations as well as the Clausius-Duhem inequality for moving isotropic MR fluids in an electro-magnetic field. We present the material constitutive relations for a non-Newtonian incompressible MR fluid. To illustrate the validity of the constitutive relations, the flow of a MR fluid between two parallel fixed plates under the influence of a constant magnetic field perpendicular to the flow direction is considered.Received: 14 July 2003, Accepted: 18 May 2004, Published online: 22 February 2005 Correspondence to: A. Dorfmann     

Linear Equations of Motion in Finite Elasticity

In this note we show that it may be possible and useful to construct valid strain-energy functions that lead directly to linear equilibrium equations for problems in isotropic homogeneous unconstrained nonlinear elasticity. While it is possible to make some general progress the final outcome will depend on the geometry and kinematics of the problem under consideration. Specific examples are given to show how exact solutions, via the linear equations of motion, can be found to non-trivial problems for physically meaningful constitutive models.      

On the predictive/fitting capabilities of the advanced differential constitutive equations for linear polyethylene melts

This contribution examines the capabilities of three differential constitutive models (XPP, PTT-XPP, and modified Leonov) in predicting rheological properties of two virtually linear polyethylene materials (HDPE Tipelin FS 450-26, mLLDPE Exact 0201) with specific attention to both steady as well as transient shear and uniaxial elongational flow situations. For each situation the (dis)advantages of the individual models are discussed and both, qualitative and quantitative model efficiency evaluation has been carried out.     

CONSTITUTIVE RELATION OF UNSATURATED SOIL BY USE OF THE MIXTURE THEORY(Ⅰ)-NONLINEAR CONSTITUTIVE EQUATIONS AND FIELD EQUATIONS

The nonlinear constitutive equations and field equations of unsaturated soils were constructed on the basis of mixture theory. The soils were treated as the mixture composed of three constituents. First, from the researches of soil mechanics, some basic assumptions about the unsaturated soil mixture were made, and the entropy inequality of unsaturated soil mixture was derived. Then, with the common method usually used to deal with the constitutive problems in mixture theory, the nonlinear constitutive equations were obtained. Finally, putting the constitutive equations of constituents into the balance equations of momentum, the nonlinear field equations of constituents were set up. The balance equation of energy of unsaturated soil was also given, and thus the complete equations for solving the thermodynamic process of unsaturated soil was formed. Foundation items: the National Natural Science Foundation of China (59678003); Special Research Plan of the Education Department of Shaanxi Province (01JK178) Biographies: HUANG Yi (1936-) ZHANG Yin-ke (1964-)     

Constitutive equations for fiber suspensions in viscoelastic media

The kinetic theory of elastic dumbells with a friction factor that depends on the fiber orientation is used to obtain constitutive equations for fiber suspensions in a polymer matrix. We followed the approach of Fan (X.J. Fan, in P. Moldenaers and R. Keunings (Eds.), Theoretical and Applied Rheology, Proceedings XIth International Congress on Rheology, Brussels, Elsevier, Amsterdam, 1992, pp. 850–852), and derived equations for polymer solutions based on the FENE-P, FENE-CR, and Giesekus models. Start-up and steady-state free shear flows are studied to explore the effects of the fiber-polymer coupling as well as the fiber volume fraction. Predictions based on different types of closure approximations for the fourth-order fiber orientation tensor are also discussed.     

The recoil of rigid PVC

Polyvinyl chloride (PVC) is widely used in industry, but it is difficult to find any extensive discussion of theological relations that describe it. The present paper discusses the behaviour of rigid PVC in extensional deformations at various temperatures. In the tests, a step elongation was applied at an initial time, then after a delay the specimen was cut, permitting recoil. Both linear and nonlinear strain regimes were studied; the linear relaxation properties (relaxation function, storage modulus) were cross-checked with eccentric-disk measurements. In the non-linear strain regime a single-integral constitutive equation of the KBKZ Wagner type was used. Separability of time and strain effects was demonstrated in our tests and so a damping function could be found which was only a function of strain. Video recordings of recoil were made, and detailed predictions of the strain-time behaviour were checked against experiments.Various constitutive equation proposals were used in the comparison, and the separated kernel integral irreversible model performed reasonably well; some other models of a differential type, and a Doi-Edwards model were not as accurate.Dedicated to Professor Arthur S. Lodge on the occasion of his 70th birthday and his retirement from the University of Wisconsin.     

Aeroelastic instability of composite wings by the consideration of different structural constitutive assumptions

The aeroelastic instability of composite wings modeled as Circumferentially Asymmetric Stiffness (CAS) thin-walled composite beams with closed cross-section is carried out. The objective has been to investigate the effects of different assumptions of constitutive equations on the aeroelastic instability behavior. Non-classical effects such as restrained warping and transverse shear are included in the beam model. The unsteady incompressible airloads are presented using Wagner׳s function. A comparison of the results based on different constitutive equations for a number of configurations including three types of stacking sequence for a box cross-section and two types of stacking sequence for a biconvex cross-section, is performed. The effects of the values of twist as well as twist-bending stiffness coefficients have been studied carefully on the results. As an outcome of this investigation it is revealed that the different choices of structural constitutive equations which result in different values of stiffness quantities; namely, twist and twist-bending stiffness, significantly affect the predicted results. For example, a difference of up to 45% in the aeroelastic critical speed has been observed between different sets of constitutive equations in some cases.     

Multi-scale constitutive modeling of the small deformations of semi-crystalline polymers

A multi-scale constitutive model for the small deformations of semi-crystalline polymers such as high density Polyethylene is presented. Each macroscopic material point is supposed to be the center of a representative volume element which is an aggregate of randomly oriented composite inclusions. Each inclusion consists of a stack of parallel crystalline lamellae with their adjacent amorphous layers.Micro-mechanically based constitutive equations are developed for each phase. A viscoplastic model is used for the crystalline lamellae. A new nonlinear viscoelastic model for the amorphous phase behavior is proposed. The model takes into account the fact that the presence of crystallites confines the amorphous phase in extremely thin layers where the concentration of chain entanglements is very high. This gives rise to a stress contribution due to elastic distortion of the chains. It is shown that the introduction of chains’ elastic distortion can explain the viscoelastic behavior of crystalline polymers. The stress contribution from elastic stretching of the tie molecules linking the neighboring lamellae is also taken into account.Next, a constitutive model for a single inclusion considered as a laminated composite is proposed. The macroscopic stress-strain behavior for the whole RVE is found via a Sachs homogenization scheme (uniform stress throughout the material is assumed).Computational algorithms are developed based on fully implicit time-discretization schemes.     

On the fully developed tube flow of a class of non-linear viscoelastic fluids

The fully developed pipe flow of a class of non-linear viscoelastic fluids is investigated. Analytical expressions are derived for the stress components, the friction factor and the velocity field. The friction factor which depends on the Deborah and Reynolds numbers is substantially smaller than the corresponding value for the Newtonian flow field with implications concerning the volume flow rate. We show that non-affine models in the class of constitutive equations considered such as Johnson-Segalman and some versions of the Phan-Thien-Tanner models are not representative of physically realistic flow fields for all Deborah numbers. For a fixed value of the slippage factor they predict physically admissible flow fields only for a limited range of Deborah numbers smaller than a critical Deborah number. The latter is a function of the slippage.     

Constitutive modeling of discontinuities by means of discrete and continuum approximations and damage models

In this paper the mathematical modeling of discontinuities using the discrete approximation and the continuum approximation with weak discontinuities is presented. First, the kinematics of discontinuities is discussed, then two constitutive models based on the continuum damage mechanics theory are developed. The first model is an isotropic damage model and is used in the discrete approximation. The second model is an anisotropic damage model and is used in the continuum approximation. These models are characterized for weighing the mode of failure in the failure criterion. An energy analysis is proposed to establish the equations that relate the parameters of both constitutive models; the fulfillment of the involved equations guarantee that both models are energetically equivalent. It is concluded that the proposed models are suitable to reproduce the constitutive behavior of discontinuities.     

An efficient algorithm for strain history tracking in finite element computations of non-Newtonian fluids with integral constitutive equations

A new finite element technique has been developed for employing integral-type constitutive equations in non-Newtonian flow simulations. The present method uses conventional quadrilateral elements for the interpolation of velocity components, so that it can conveniently handle viscoelastic flows with both open and closed streamlines (recirculating regions). A Picard iteration scheme with either flow rate or elasticity increment is used to treat the non-Newtonian stresses as pseudo-body forces, and an efficient and consistent predictor-corrector scheme is adopted for both the particle-tracking and strain tensor calculations. The new method has been used to simulate entry flows of polymer melts in circular abrupt contractions using the K-BKZ integral constitutive model. Results are in very good agreement with existing numerical data. The important question of mesh refinement and convergence for integral models in complex flow at high flow rate has also been addressed, and satisfactory convergence and mesh-independent results are obtained. In addition, the present method is relatively inexpensive and in the meantime can reach higher elasticity levels without numerical instability, compared with the best available similar calculations in the literature.     

Numerical study of flows of complex fluids between eccentric cylinders using transformation functions

In this paper, we investigate fluid flows between eccentric cylinders by means of two stream‐tube analyses. The first method considers a one‐to‐one global transformation function that allows the physical domain to be transformed into a mapped domain, used as computational domain, that involves concentric streamlines. The second approach uses local transformations and domain decomposition techniques to deal with mixed flow regimes. Both formulations are particularly adapted for handling time‐dependent constitutive equations, since particle‐tracking problems are avoided. Mass conservation is verified in both formulations and the relevant numerical procedure can be carried out using simple meshes built on the mapped streamlines. Fluids obeying anelastic and viscoelastic constitutive equations are considered in the calculations. The numerical results are consistent with those in the literature for the flow rates tested. Application of the method to the K‐BKZ memory‐integral constitutive equation highlights significant differences between the model predictions and those provided by more simple rheological models. Copyright © 2002 John Wiley & Sons, Ltd.     

Nonisothermal two‐ and three‐dimensional flow simulations of inelastic and viscoelastic fluids by a finite‐volume method

This paper applies the finite‐volume method to computations of steady flows of viscous and viscoelastic incompressible fluids in complex two and three‐dimensional geometries. The materials adopted in the study obey different constitutive laws: Newtonian, purely viscous Carreau–Yasuda as also Upper‐Convected Maxwell and Phan‐Thien/Tanner differential models, with a Williams–Landel–Ferry (WLF) equation for temperature dependence. Specific analyses are made depending on the rheological model. A staggered grid is used for discretizing the equations and unknowns. Stockage possibilities allow us to solve problems involving a great number of degrees of freedom, up to 1 500 000 unknowns with a desk computer. In relation to the fluid properties, our numerical simulations provide flow characteristics for various 2D and 3D configurations and demonstrate the possibilities of the code to solve problems involving complex nonlinear constitutive equations with thermal effects. Copyright © 2009 John Wiley & Sons, Ltd.     

Plane strain bending of strain-stiffening rubber-like rectangular beams

This paper is concerned with investigation of the effects of strain-stiffening for the classical problem of plane strain bending by an end moment of a rectangular beam composed of an incompressible isotropic nonlinearly elastic material. For a variety of specific strain-energy densities that give rise to strain-stiffening in the stress–stretch response, the stresses and resultant moments are obtained explicitly. While such results are well known for classical constitutive models such as the Mooney-Rivlin and neo-Hookean models, our primary focus is on materials that undergo severe strain-stiffening in the stress–stretch response. In particular, we consider in detail two phenomenological constitutive models that reflect limiting chain extensibility at the molecular level and involve constraints on the deformation. The amount of bending that beams composed of such materials can sustain is limited by the constraint. Potential applications of the results to the biomechanics of soft tissues are indicated.     

CONSTITUTIVE RELATION OF UNSATURATED SOIL BY USE OF THE MIXTURE THEORY(Ⅰ)—NONLINEAR CONSTITUTIVE EQUATIONS AND FIELD EQUATIONS

IntroductionSoilisthemostcommonlyusedconstructionmaterialincivilengineeringandhydraulicengineering .Thecharacteristicsofsoilhavebeeninvestigatingfornearlyonehundredyears.Butbecauseofitscomplexstructure,changeableenvironmentandbeingsensitivetotheoutsideconditions,thesoiloftenshowsvariedproperties[1,2 ].Themaindifficultytothedevelopmentofgeotechnicalmechanicsishowtosetupconstitutiveequationswhichcouldsatisfactorilyaccountforengineeringpropertiesofsoil[3].Manyconstitutivemodelshavebeenformedinth…     

Transonic Euler computation in streamfunction co-ordinates

Numerical simulation by a finite element method is used to examine the problem of the rotating flow of a viscoelastic fluid in a cylindrical vessel agitated with a paddle impeller. The mathematical model consists of a viscoelastic constitutive equation of Oldroyd B type coupled to the hydrodynamic equations expressed in a rotating frame. This system is solved by using an unsteady approach for velocity, pressure and stress fields. For Reynolds numbers in the range 0.1–10, viscoelastic effects are taken into account up to a Deborah number De of 1.33 and viscoelasticity and inertia cross-effects are studied. Examining the velocity and stress fields as well as the power consumption, it is found that their evolutions are significantly different for low and moderate inertia. These results confirm the trends of experimental studies and show the specific contribution of elasticity without interference of the pseudoplastic character found in actual fluids.