Constitutive equations of non-isothermal polymer melt

Constitutive equations of non-isothermal polymer melt are presented by the analysis of entropic free energy contribution of the macromolecular chains, which are treated as elastic dumbbell models. With describing non-isothermal dumbbell spring, as the function of temperature, the non-linear elastic coefficient expression causes the appearance of temperature gradient in stress constitutive equations. Following the constitutive equation of Hookean dumbbell model, non-isothermal stress constitutive equations of FENE and FENE-P models are derived. In deriving process of constitutive equations, the second moment approximation is used to closure FENE model. Using the non-isothermal constitutive equations, numerical simulations of polymer flow through shear cavity and planar contraction cavity are presented. And the distributions of correlative stress functions and the effects of different temperatures on stress functions are discussed. The present results are shown to explore the non-isothermal constitutive equations of elastic dumbbell models, and to search more accurately describing way of non-isothermal polymer melt.     

Recent advances in the numerical modeling of constitutive relations

In this paper we present an overview of the recent developments in the area of numerical and finite element modeling of nonlinear constitutive relations. The paper discusses elastic, hyperelastic, elastoplastic and anisotropic plastic material models. In the hyperelastic model an emphasis is given to the method by which the incompressibility constraint is applied. A systematic and general procedure for the numerical treatment of hyperelastic model is presented. In the elastoplastic model both infinitesimal and large strain cases are discussed. Various concerns and implications in extending infinitesimal theories into large strain case are pointed out. In the anisotropic elastoplastic case, emphasis is given to the practicality of proposed theories and its feasible and economical use in the finite element environment.     

Elasto-viscoplastic fluid flow in tubes of arbitrary cross-section

Elasto-viscoplastic fluid flow in tubes of arbitrary cross-section is explored. A constitutive equation superposing elasticity on the viscoplastic behavior through a linear combination of the non-linear simplified Phan–Thien–Tanner constitutive model of viscoelasticity and the non-linear Bingham constitutive equation of viscoplasticity is developed. The equation of motion is solved analytically for the longitudinal field for steady flow in tubes of non-circular cross-section. The evolution of the plug and stagnant zones, hallmarks of viscoplastic behavior, is studied when elasticity is present, and the rate of flow is determined in terms of the Weissenberg and Bingham numbers. Elasticity tends to enhance the rate of flow for given viscoplastic conditions, and stagnant and plug zone configurations are substantially altered.     

A criterion for the strong ellipticity of the equilibrium equations of an isotropic non-linearly elastic material

A new criterion for satisfying the condition of strong ellipticity of the equilibrium equations of a homogeneous isotropic elastic medium is obtained, consisting of twelve unconditional and the same number of conditional inequalities. The criterion differs from previously known criteria in the fact that it is more effective for the analytical verification of the condition of strong ellipticity. A number of sufficient criteria of the strong ellipticity of the equations of the non-linear statics of an elastic body is presented. It is shown that the method used in this paper is also effective in the case of an incompressible material. A new constitutive inequality in the non-linear theory of elasticity is proposed that, in particular, ensures the strong ellipticity of the equilibrium equations and the stability of an elastic body in a state of hydrostatic compression. A mechanical interpretation of some implications of the above mentioned inequality is given. Using the criterion obtained, the condition of strong ellipticity is verified for a number of specific models of compressible elastic materials and examples of the violation of this condition are also considered.     

On the use of characteristic variables in viscoelastic flow problems

The system of partial differential equations governing the flow of an upper converted Maxwell fluid is known to be of mixed elliptic–hyperbolic type. The hyperbolic nature of the constitutive equation requires that, where appropriate, inflow conditions are prescribed in order to obtain a well-posed problem. Although there are three convective derivatives in the constitutive equation there are only two characteristic quantities whichare transported along the streamlines. These characteristicquantities are identified. A spectral element method is describedin which continuity of the characteristic variables is usedto couple the extra stress components between contiguous elements.The continuity of the characteristic variables is treated asa constraint on the constitutive equation. These conditionsdo not necessarily impose continuity on the extra-stress components.The velocity and pressure follow from the doubly constrainedweak formulation which enforces a divergence-free velocityfield and irrotational polymeric stress forces. This meansthat both the pressure and the extra-stress tensor are discontinuous.Numerical results are presented to demonstrate this procedure.The theory is applied to the upper convected Maxwell modelwith vanishing Reynolds number. No regularization techniques such as streamline upwind Petrov Galerkin (SUPG), elastic viscous split stress (EVSS) or explicitly elliptic momentum equation(EEME) are used.     

An efficient numerical method to solve 2-D interval bi-modular problems via orthogonal polynomial expansion

An interval analysis of uncertain bi-modular problems is presented by addressing the efficiency of deterministic solution and reduction of computational cost on the non-linear FE iteration. Firstly, the singularity of 2-D conventional bi-modular constitutive matrix is pointed out via a concise mathematical illustration, and is removed via a complement of shear modulus consistent with the coaxial condition. A new FE model with a full rank constitutive matrix is developed to solve deterministic bi-modular problems, which is well performed in the numerical tests, particularly in term of convergence. Secondly, an orthogonal polynomial expansion based surrogate is constructed to alleviate the heavy computational burden caused by repeated non-linear FE solution in the optimization process for bounds estimation. Numerical examples are given to illustrate the accuracy and efficiency of proposed approach, and a good accordance can be observed between the results obtained by the proposed approach and reference solutions.     

The sector constants of continuous state branching processes with immigration

Continuous state branching processes with immigration are studied. We are particularly concerned with the associated (non-symmetric) Dirichlet form. After observing that gamma distributions are only reversible distributions for this class of models, we prove that every generalized gamma convolution is a stationary distribution of the process with suitably chosen branching mechanism and with continuous immigration. For such non-reversible processes, the strong sector condition is discussed in terms of a characteristic called the Thorin measure. In addition, some connections with notion from non-commutative probability theory will be pointed out through calculations involving the Stieltjes transform.     

Instability prevention at the modeling of large elasto-plastic strains under cyclic loading

The cyclic instability phenomenon is investigated at the modeling of large elasto-plastic strains. The possible causes of the cyclic instability and conditions ensuring cyclical stability of elasto-plastic models are analyzed for the case of large strains. Among the possible causes of the cyclic instability the following are considered: the method of strain decomposition on elastic and plastic parts; the constitutive law for the elastic deformation (hypo- and hyper-elasticity); the constitutive equation for the plastic deformation; the constitutive relation for the plastic spin; kinematic hardening law, in particular, the type of the objective rate in the generalized Prager's law. Predictions of 50 various models of the elasto-plastic material have been compared in order to find the causes of the cyclic instability. Two test problems are considered: cyclic simple shear, combined cyclic simple shear and tension-compression. Results of numerical experiments for the various material models are presented and discussed. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Rate‐type Model for electro‐rheological Material Behaviour consistent with Extended Thermodynamics: Application to a steady viscometric flow

Based on extended thermodynamics, a class of rate‐type constitutive equations for electrorheological materials is presented. It is shown unlike many recent works using standard methods of continuum mechanics, that this approach, which consists in postulating expressions which relate the time derivative of the Cauchy stress to the motion, takes well account of the viscoelastic‐plastic effects as well as the coupling between the motion and the electric field. By considering special cases it turns out that the derived nonlinear algebraic model embeds a variety of existing differential models in the literature. The accuracy of the model is illustrated by considering a steady shear flow in which the electric field is applied perpendicularly to the flow direction. Analytical solutions with various parameter variations are obtained and discussed. The model predicts well the behaviour associated with material instabilities observed in viscometric flows for which the existence of the yield stress is required as confirmed by many experiments.     

On the existence of similar solutions for a three dimensional flow problem

Using a topological approach, the existence of similar solutions of a system of non-linear, third order, coupled ordinary differential equations which characterise the steady, three dimensional flow of an incompressible, viscous fluid due to a stretching flat surface is discussed.     

无力磁场和Beltrami流

场矢量处处平行其旋度的管式矢量场具有复杂的拓扑结构．而且通常表现出混饨行为．本文提供了具有常数比例因子的这类矢量场在三种基本坐标系内的解析解,并指出完全导电流体内的一个Beltrami流能维持一个定常无力磁场,如果磁场矢量处处平行速度．     

聚合物材料损伤理论

本文研究了外界环境如载荷、温度、湿度、辐射、腐蚀介质、非相容介质等对聚合物材料变形的影响,给出了微分型和积分型损伤演化方程;讨论了非线性弹性体和非线性粘弹性损伤变形以及腐蚀损伤.给出了粘弹性损伤变分原理.     

Torsional buckling and post-buckling of columns made of aluminium alloy

The paper concerns torsional buckling and the initial post-buckling of axially compressed thin-walled aluminium alloy columns with bisymmetrical cross-section. It is assumed that the column material behaviour is described by the Ramberg–Osgood constitutive equation in non-linear elastic range. The stationary total energy principle is used to derive the governing non-linear differential equation. An approximate solution of the equation determined by means of the perturbation approach allows to determine the buckling loads and the initial post-buckling behaviour. Numerical examples dealing with simply supported I-column are presented and the effect of material elastic non-linearity on the critical loads and initial post-buckling behaviour are compared to the linear solution.     

Orientational Design of Nonorthotropic Two-Dimensional Materials

A linear elastic nonorthotropic two-dimensional material with six independent elastic parameters is considered. Different criteria for comparison of constitutive matrices, which may be mutually rotated, are discussed. Optimal material orientations of a nonorthotropic material, which correspond to extreme values of the energy density, are determined. The results obtained are extended to some nonlinear elastic material models.     

Towards variational constitutive updates for finite strain plasticity models based on non-associative evolution equations

In this contribution, first steps towards variational constitutive updates for finite strain plasticity theory based on non-associative evolution equations are presented. These schemes allow to compute the unknown state variables such as the plastic part of the deformation gradient, together with the deformation mapping, by means of a fully variational minimization principle. Therefore, standard optimization algorithms can be applied to the numerical implementation leading to a very robust and efficient numerical implementation. Particularly, for highly non-linear, singular or nearly ill-posed physical models like that corresponding to crystal plasticity showing a large number of possible active slip planes, this is a significant advantage compared to standard constitutive updates such as the by now classical return-mapping algorithm. While variational constitutive updates have been successfully derived for associative plasticity models, their extension to more complex constitutive laws, particularly to those featuring non-associative evolution equations, is highly challenging. In the present contribution, a certain class of non-associative finite strain plasticity models is discussed and recast into a variationally consistent format. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical solution of viscoplastic constitutive equations with internal state variables Part I: algorithms and implementation

This article treats the initial-boundary-value problem of viscoplasticity using unified constitutive models without a yield surface. Semi-discretization with the finite element method (FEM) leads to a system of differential-algebraic equations (DAE) with strongly non-linear evolution equations for the internal state variables. A special family of partitioned Runge–Kutta methods is introduced which allows an efficient time integration of the semidiscrete system. Coefficients for methods of order one, two, and three are given. Finally, numerical results for some two- and three-dimensional examples using the model of Hart are presented. In a second part we will give the theoretical background and state a proof of convergence for the algorithm presented in this paper. © 1997 by B. G. Teubner Stuttgart–John Wiley & Sons Ltd.     

Effects of partial slip,viscous dissipation and Joule heating on Von Kármán flow and heat transfer of an electrically conducting non-Newtonian fluid

The steady Von Kármán flow and heat transfer of an electrically conducting non-Newtonian fluid is extended to the case where the disk surface admits partial slip. The fluid is subjected to an external uniform magnetic field perpendicular to the plane of the disk. The constitutive equation of the non-Newtonian fluid is modeled by that for a Reiner–Rivlin fluid. The momentum equations give rise to highly non-linear boundary value problem. Numerical solutions for the governing non-linear equations are obtained over the entire range of the physical parameters. The effects of slip, magnetic parameter and non-Newtonian fluid characteristics on the velocity and temperature fields are discussed in detail and shown graphically. Emphasis has been laid to study the effects of viscous dissipation and Joule heating on the thermal boundary layer. It is interesting to find that the non-Newtonian cross-viscous parameter has an opposite effect to that of the slip and the magnetic parameter on the velocity and the temperature fields.     

研究金属中激波构造与衰减的一个物理模型

本文给出了研究金属中激波构造与衰减的一个物理模型.为了建立高速形变下材料的本构方程和研究激波过渡带的构造,需要考虑二个独立的理论方面.首先,将比内能分解成弹性压缩能和弹性形变能,而将形变能作为弹性应变和熵的函数展开到三阶项,其中考虑了热与机械能的耦合效应.其次,从位错动力学角度建议了一个塑性松弛函数以便描述高温、高压下塑性流动的特性.另外,本文给出了一个常微分方程组用以计算定态激波过渡带中各状态变量的分布以及激波的厚度.倘若假定在激波上熵的跳跃可以忽略,并用Hugoniot压缩模量代替等熵压缩摸量,可以获得一个分析解.最后,本文还提出了求解平板对称碰撞中激波波头衰减的一个近似方法。     

Couette and Poiseuille flows of an Oldroyd 6-constant fluid with magnetic field

In this paper, we develop a set of differential equations describing the steady flow of an Oldroyd 6-constant magnetohydrodynamic fluid. The fluid is electrically conducting in the presence of a uniform transverse magnetic field. The developed non-linear differential equation takes into account the effect of the material constants and the applied magnetic field. We presented the solution for three types of steady flows, namely,

(i)

Couette flow

(ii)

Poiseuille flow and

(iii)

generalized Couette flow.

Homotopy analysis method (HAM) is used to solve the non-linear differential equation analytically. It is found from the present analysis that for steady flow the obtained solutions are strongly dependent on the material constants (non-Newtonian parameters) which is different from the model of Oldroyd 3-constant fluid. Numerical solutions are also given and compared with the solutions by HAM.     

Numerical Studies of Viscoelastic Flow Using the Software OpenFOAM

Viscoelastic fluids are a special class of non-Newtonian fluids. There are several types of viscoelastic fluid models, and all of them have a complex rheological response in comparison to Newtonian fluids. This response can be viewed as a combination of viscous and elastic effects and non-linear phenomena. This complex physics makes a numerical simulation a rather challenging task, even in simple test-cases. Studies presented in this paper are numerical studies of the viscoelastic fluid flow in several test cases. These studies have been done in OpenFOAM, an open-source CFD package. Implementation of viscoelastic models and a solver is only available in a community driven version of software (OpenFOAM-ext). One of the goals of research in this paper was to test the solver and models on some simple test cases. We considered start-up and pulsating flows of viscoelastic fluid in a channel and a circular pipe. The important thing is that an analytical solution can be found in these cases, making in possible to test all aspects of numerical simulation in OpenFOAM. Obtained results showed an excellent agreement with the analytical solution for both velocity and stress components. These results encouraged authors' motivation and a choice to use OpenFOAM for simulation of viscoelastic flows. We hope that our research will make a contribution to the OpenFOAM community. Our plan for the further research is a simulation of blood flow in arteries with the viscoelastic solver. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)