Comparative study of some constitutive equations characterising non-Newtonian fluids

This paper compares, in a general way, the predictions of the constitutive equations given by Rivlin and Ericksen, Oldroyd, and Walters. Whether we consider the rotational problems in cylindrical co-ordinates or in spherical polar co-ordinates, the effect of the non-Newtonicity on the secondary flows is collected in a single parameterα which can be explicitly expressed in terms of the non-Newtonian parameters that occur in each of the above-mentioned constitutive equations. Thus, for a given value ofα, all the three fluids will have identical secondary flows. It is only through the study of appropriate normal stresses that a Rivlin-Ericksen fluid can be distinguished from the other two fluids which are indistinguishable as long as this non-Newtonian parameter has the same value.

The application of B-P constitutive equations in finite element analysis of high velocity impact

We present in this paper the application of B-P constitutive equations in finite element analysis of high velocity impact. The impact process carries out in so quick time that the heat-conducting can be neglected and meanwhile, the functions of temperature in equations need to be replaced by functions of plastic work. The material constants in the revised equations can be determined by comparison of the one-dimensional calculations with the experiments of Hopkinson bar. It can be seen from the comparison of the calculation with the experiment of a tungsten alloy projectile impacting a three-layer plate that the B-P constitutive equations in that the functions of temperature were replaced by the functions of plastic work can be used to analysis of high velocity impact.

On numerical implementation of an isotropic elastic-viscoplastic constitutive model for bulk metallic glasses

An efficient algorithm within the finite deformation framework is developed for finite element implementation of a recently proposed isotropic, Mohr-Coulomb type material model, which captures the elastic-viscoplastic, pressure sensitive and plastically dilatant response of bulk metallic glasses. The constitutive equations are first reformulated and implemented using an implicit numerical integration procedure based on the backward Euler method. The resulting system of nonlinear algebraic equations is solved by the Newton-Raphson procedure. This is achieved by developing the principal space return mapping technique for the present model which involves simultaneous shearing and dilatation on multiple potential slip systems. The complete stress update algorithm is presented and the expressions for viscoplastic consistent tangent moduli are derived. The stress update scheme and the viscoplastic consistent tangent are implemented in the commercial finite element code ABAQUS/Standard. The accuracy and performance of the numerical implementation are verified by considering several benchmark examples, which includes a simulation of multiple shear bands in a 3D prismatic bar under uniaxial compression.

Dynamic constitutive equation of GFRP obtained by Lagrange experiment

The note presents a method of constructing dynamic constitutive equations of material by means of Lagrange experiment and analysis. Tests were carried out by a light gas gun and the stress history profiles were recorded on multiple Lagrange positions. The dynamic constitutive equations were deduced from the regression of a series of data which was obtained by Lagrange analysis based upon recorded multiple stress histories. Here constitutive equations of glass fibre reinforced phenolic resin composite(GFRP) in uniaxil strain state under dynamic loading are given. The proposed equations of the material agree well with experimental results.

Energy Functions And Basic Equations For Ferroelectrics

Energy functions (or characteristic functions) and basic equations for ferroelectrics in use today are given by those for ordinary dielectrics in the physical and mechanical communications. Based on these basic equations and energy functions, the finite element computation of the nonlinear behavior of the ferroelectrics has been carried out by several research groups. However, it is difficult to process the finite element computation further after domain switching, and the computation results are remarkably deviating from the experimental results. For the crack problem, the iterative solution of the finite element calculation could not converge and the solutions for fields near the crack tip oscillate. In order to finish the calculation smoothly, the finite element formulation should be modified to neglect the equivalent nodal load produced by spontaneous polarization gradient. Meanwhile, certain energy functions for ferroelectrics in use today are not compatible with the constitutive equations of ferroelectrics and need to be modified. This paper proposes a set of new formulae of the energy functions for ferroelectrics. With regard to the new formulae of the energy functions, the new basic equations for ferroelectrics are derived and can reasonably explain the question in the current finite element analysis for ferroelectrics.

A numerical method for solving the dynamic three-dimensional Ericksen-Leslie equations for nematic liquid crystals subject to a strong magnetic field

A finite difference technique, based on a projection method, is developed for solving the dynamic three-dimensional Ericksen-Leslie equations for nematic liquid crystals subject to a strong magnetic field. The governing equations in this situation are derived using primitive variables and are solved using the ideas behind the GENSMAC methodology (Tome and McKee [32]; Tome et al. [34]). The resulting numerical technique is then validated by comparing the numerical solution against an analytic solution for steady three-dimensional flow between two-parallel plates subject to a strong magnetic field. The validated code is then employed to solve channel flow for which there is no analytic solution. (C) 2009 Elsevier B.V. All rights reserved.

Interpretation of the pressuremeter test using numerical models based on deformation tensor equations

The pressuremeter test in boreholes has proven itself as a useful tool in geotechnical explorations, especially comparing its results with those obtained from a mathematical model ruled by a soil representative constitutive equation. The numerical model shown in this paper is aimed to be the reference framework for the interpretation of this test. The model analyses variables such as: the type of response, the initial state, the drainage regime and the constitutive equations. It is a model of finite elements able to work with a mesh without deformation or one adapted to it.

A coupled elastoplastic-damage constitutive model with Lode angle dependent failure criterion

A coupled elastoplastic-damage constitutive model with Lode angle dependent failure criterion for high strain and ballistic applications is presented. A Lode angle dependent function is added to the equivalent plastic strain to failure definition of the Johnson–Cook failure criterion. The weakening in the elastic law and in the Johnson–Cook-like constitutive relation implicitly introduces the Lode angle dependency in the elastoplastic behaviour. The material model is calibrated for precipitation hardened Inconel 718 nickel-base superalloy. The combination of a Lode angle dependent failure criterion with weakened constitutive equations is proven to predict fracture patterns of the mechanical tests performed and provide reliable results. Additionally, the mesh size dependency on the prediction of the fracture patterns was studied, showing that was crucial to predict such patterns

Development of a biaxial compression device for biological samples : preliminary experimental results for a closed cell foam

Biological tissues are subjected to complex loading states in vivo and in order to define constitutive equations that effectively simulate their mechanical behaviour under these loads, it is necessary to obtain data on the tissue's response to multiaxial loading. Single axis and shear testing of biological tissues is often carried out, but biaxial testing is less common. We sought to design and commission a biaxial compression testing device, capable of obtaining repeatable data for biological samples. The apparatus comprised a sealed stainless steel pressure vessel specifically designed such that a state of hydrostatic compression could be created on the test specimen while simultaneously unloading the sample along one axis with an equilibrating tensile pressure. Thus a state of equibiaxial compression was created perpendicular to the long axis of a rectangular sample. For the purpose of calibration and commissioning of the vessel, rectangular samples of closed cell ethylene vinyl acetate (EVA) foam were tested. Each sample was subjected to repeated loading, and nine separate biaxial experiments were carried out to a maximum pressure of 204 kPa (30 psi), with a relaxation time of two hours between them. Calibration testing demonstrated the force applied to the samples had a maximum error of 0.026 N (0.423% of maximum applied force). Under repeated loading, the foam sample demonstrated lower stiffness during the first load cycle. Following this cycle, an increased stiffness, repeatable response was observed with successive loading. While the experimental protocol was developed for EVA foam, preliminary results on this material suggest that this device may be capable of providing test data for biological tissue samples. The load response of the foam was characteristic of closed cell foams, with consolidation during the early loading cycles, then a repeatable load-displacement response upon repeated loading. The repeatability of the test results demonstrated the ability of the test device to provide reproducible test data and the low experimental error in the force demonstrated the reliability of the test data.

Polymer melt dynamics: Microscopic roots of fractional viscoelasticity

The rheological properties of polymer melts and other complex macromolecular fluids are often successfully modeled by phenomenological constitutive equations containing fractional differential operators. We suggest a molecular basis for such fractional equations in terms of the generalized Langevin equation (GLE) that underlies the renormalized Rouse model developed by Schweizer [J. Chem. Phys. 91, 5802 (1989)]. The GLE describes the dynamics of the segments of a tagged chain under the action of random forces originating in the fast fluctuations of the surrounding polymer matrix. By representing these random forces as fractional Gaussian noise, and transforming the GLE into an equivalent diffusion equation for the density of the tagged chain segments, we obtain an analytical expression for the dynamic shear relaxation modulus G(t), which we then show decays as a power law in time. This power-law relaxation is the root of fractional viscoelastic behavior.

The oscillation of a spheroid along its axis in a non-Newtonian fluid

The present paper investigates the nature of the fluid flow when a spheroid is suspended in an infinitely extending elastico-viscous fluid defined by the constitutive equations given by Oldroyd or Rivlin and Ericksen, and is made to perform small amplitude oscillations along its axis. The solution of the vector wave equation is expressed in terms of the solution of the corresponding scalar wave equation, without the use of Heine's function or spheroidal wave functions. Two special cases (i) a sphere and (ii) a spheroid of small ellipticity, are studied in detail.

Prediction of solid block masonry prism compressive strength using FE model

Masonry strength is dependent upon characteristics of the masonry unit,the mortar and the bond between them. Empirical formulae as well as analytical and finite element (FE) models have been developed to predict structural behaviour of masonry. This paper is focused on developing a three dimensional non-linear FE model based on micro-modelling approach to predict masonry prism compressive strength and crack pattern. The proposed FE model uses multi-linear stress-strain relationships to model the non-linear behaviour of solid masonry unit and the mortar. Willam-Warnke's five parameter failure theory developed for modelling the tri-axial behaviour of concrete has been adopted to model the failure of masonry materials. The post failure regime has been modelled by applying orthotropic constitutive equations based on the smeared crack approach. Compressive strength of the masonry prism predicted by the proposed FE model has been compared with experimental values as well as the values predicted by other failure theories and Eurocode formula. The crack pattern predicted by the FE model shows vertical splitting cracks in the prism. The FE model predicts the ultimate failure compressive stress close to 85 of the mean experimental compressive strength value.

A general treatment of two dimensional plane flows, for a micropolar fluid

In this paper we have studied the flow of a micropolar fluid, whose constitutive equations were given by Eringen, in two dimensional plane flow. In two notes, we have discussed the validity of the boundary condition v=a ω and its effect on the entire flow field. We have restricted our study to the case when Stokes' approximation is valid, i. e. slow motion for it is difficult to uncouple the equations in the most general case.