A Frame-Independent Solution to Saint-Venant’s Flexure Problem

The paper illustrates a solution approach for the Saint-Venant flexure problem which preserves a pure objective tensor form, thus yielding, for sections of arbitrary geometry, representations of stress and displacement fields that exploit exclusively frame-independent quantities. The implications of the availability of an objective solution to the shear warpage problem are discussed and supplemented by several analytical and numerical solutions. The derivation of tensor expressions for the shear center and the shear flexibility tensor is also illustrated. Furthermore, a Cesaro-like integration procedure is provided whereby the derivation of a frame-independent representation of the displacements field for the shear loading case is systematically carried out via the use of Gibbs’ algebra. The objective framework presented in this paper is further exploited in a companion article (Serpieri, in J. Elast. (2013)) to prove the coincidence of energetic and kinematic definitions of the shear flexibility tensor and of the shear principal axes.     

The symmetry approximation for nonsymmetric permeability tensors and its consequences for mass transport

It is well known that the permeability has a tensor character. In practical applications, this is accounted for by the introduction of three principal permeabilities — three scalars — and three mutually orthogonal principal axes. In this paper, it is investigated whether this is always the exact way of describing anisotropy and, if not, what the consequences of the principal axes approximation are for flow and transport. First, it is shown that spatial upscaling may result in nonsymmetric large-scale permeability tensors, for which principal axes do not exist. However, it is possible to define generalized principal axes: three principal axes for the flux and three for the pressure gradient, with only three principal permeabilities. Since nonsymmetric permeability tensors are undesirable in practical applications, an approximation method making the nonsymmetric permeability symmetric is introduced. The important conclusion is then that the exact large-scale flux and large-scale pressure gradient do not have the same directions as the approximate flux and approximate pressure gradient. A practical consequence is that the principal axes approximation results in a difference between flux and transport direction. When considering miscible displacement or transport of mass dissolved in groundwater, the velocity component normal to the flux direction may be considered as a contribution to the transverse macro dispersion.     

Generalized Darcy's Law and the Structure of the Phase and Relative Phase Permeabilities for Two-Phase Flows through Anisotropic Porous Media

For two-phase immiscible fluid flows a generalized Darcy's law is written in invariant tensor form for crystallographic point symmetry groups and anisotropic textures. The representation of the phase permeability coefficient tensors and the structure of the expressions for the relative phase permeabilities are analyzed for all symmetry groups. The relation between the phase and absolute permeability coefficient tensors is specified by a fourth-rank tensor with the external symmetry coinciding with external symmetry of the phase permeability tensors. It is shown that the external symmetry of the phase permeability coefficient tensors can differ from the external symmetry of the absolute permeability tensor. For triclinic and monoclinic symmetry groups it is shown that the phase permeability coefficient tensors may not be coaxial with each other and with the absolute permeability tensor; moreover, the directions of the principal axes of the phase permeability coefficient tensors can depend on the saturation.     

Anisotropy effects in two-phase flows through porous media

Effects manifested in two-phase flows through anisotropic porous reservoirs with monoclinic and triclinic characteristics are analyzed. It is shown that in two-phase flows through media with monoclinic and triclinic symmetries of flow characteristics the position of the principal axes of the phase permeability tensors depends on the saturation and does not coincide with the position of the principal axes of the absolute permeability tensor in single-phase flows and that going over from single-to two-phase flow may lead to a change in the symmetry group of the flow characteristics. A general representation of the phase permeability tensor components is presented and formulas are given for the diagonal and nondiagonal components of the relative phase permeabilities, which are universal and can be used for anisotropic media with any type of anisotropy (symmetry) of flow characteristics. A complex of laboratory tests for finding the nondiagonal components of the phase and relative phase permeability tensors is discussed.     

Generalization of strain-gradient theory to finite elastic deformation for isotropic materials

This paper concerns finite deformation in the strain-gradient continuum. In order to take account of the geometric nonlinearity, the original strain-gradient theory which is based on the infinitesimal strain tensor is rewritten given the Green–Lagrange strain tensor. Following introducing the generalized isotropic Saint Venant–Kirchhoff material model for the strain-gradient elasticity, the boundary value problem is investigated in not only the material configuration but also the spatial configuration building upon the principle of virtual work for a three-dimensional solid. By presenting one example, the convergence of the strain-gradient and classical theories is studied.     

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.     

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     

Potentiality of isotropic nonlinear tensor functions relating two deviators

In the theory of constitutive relations, isotropic quadratic nonlinear tensor functions modeling media with second-order effects, in particular, with misalignment of the force and kinematic tensors, are considered. It is very interesting to consider tensor functions with a scalar potential relating two symmetric deviators of rank two. In this case, the potentiality conditions are integrated, and it is shown that the first integral contains two arbitrary functions of the quadratic invariant of the tensor argument and one arbitrary function of the cubic invariant. A tensorially nonlinear generalization of the rigid-viscoplasticmodel (a two-contact Binghamsolid) is carried out.     

Application of corotational rates of the logarithmic strain in constitutive modeling of hardening materials at finite deformations

In this paper a finite deformation constitutive model for rigid plastic hardening materials based on the logarithmic strain tensor is introduced. The flow rule of this constitutive model relates the corotational rate of the logarithmic strain tensor to the difference of the deviatoric Cauchy stress and the back stress tensors. The evolution equation for the kinematic hardening of this model relates the corotational rate of the back stress tensor to the corotational rate of the logarithmic strain tensor. Using Jaumann, Green–Naghdi, Eulerian and logarithmic corotational rates in the proposed constitutive model, stress–strain responses and subsequent yield surfaces are determined for rigid plastic kinematic and isotropic hardening materials in the simple shear problem at finite deformations.     

使用主值空间表示的各向同性塑性本构方程

针对各向同性材料,在内变量为标量的假定下,应用张量函数表示定理给出了其塑性应变增量的不变性表示.它的3个不可约基张量取决于应力张量、相互正交且共主轴.建立3个基张量构成的张量子空间与三维主值空间的对应关系,将共主轴的张量采用笛卡尔坐标系中的矢量描述,矢量在不同坐标系下的分量均为张量的一组不可约不变量.定义塑性应变增量对应的矢量为内变量增量,使用张量函数表示理论得到,内变量演化方程除取决于应力对应的矢量和内变量本身外,还取决于应力增量在张量子空间中的投影,该投影就是应力对应矢量的增量,因此,本构方程归结为确定主值空间中矢量之间的关系.最后表明,三维主值空间与张量子空间中的流动法则是等价的.     

The role of a realistic volume constraint in modelling a two dimensional granular assembly

This paper presents a deceptively simple mathematical model for the deformation of granular materials composed of rigid particles. The model captures many of the diverse features of the behaviour of such a material and emphasises the importance of volume constraints in situations where the deformation is mainly by particle rearrangement. It is constructed using a simple dissipation function and a rather more complicated dilatancy rule containing an updateable reference strain. This allows the solid-like and fluid-like properties of granular materials to be reconciled in a single model.The model has been used to simulate experiments that use an analogue of an ideal granular material [Joer, H.A., Lanier, J., Fahey, M., 1998. Deformation of granular materials due to rotation of principal axes. Geotechnique 48 (5), 605-619] consisting of a two dimensional assembly of thin PVC rods. These experiments clearly illustrate: partially reversible dilatancy in direct shear tests; cyclic shearing leading to liquefaction in constant volume shear tests; and non-coaxiality of the principal axes of stress and strain increment in circular loading tests. These radically different modes of deformation provide a challenging data set that allows the model's potential to be clearly demonstrated.The authors believe that the comparison of these experimental results and our simulations give strong support to the assertion that volume changes associated with shear deformation are responsible for the rotational kinematic hardening seen in granular materials, and hence, the non-coaxiality of the stress and strain-rate tensors.     

Basis free relations for the conjugate stresses of the strains based on the right stretch tensor

In this paper, general relations between two different stress tensors T^f and T^g, respectively conjugate to strain measure tensors f(U) and g(U) are found. The strain class f(U) is based on the right stretch tensor U which includes the Seth–Hill strain tensors. The method is based on the definition of energy conjugacy and Hill’s principal axis method. The relations are derived for the cases of distinct as well as coalescent principal stretches. As a special case, conjugate stresses of the Seth–Hill strain measures are then more investigated in their general form. The relations are first obtained in the principal axes of the tensor U. Then they are used to obtain basis free tensorial equations between different conjugate stresses. These basis free equations between two conjugate stresses are obtained through the comparison of the relations between their components in the principal axes, with a possible tensor expansion relation between the stresses with unknown coefficients, the unknown coefficients to be obtained. In this regard, some relations are also obtained for T⁽⁰⁾ which is the stress conjugate to the logarithmic strain tensor lnU.     

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.     

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.     

NONLINEAR MICRO－MECHANICAL MODEL FOR PLAIN WOVEN FABRIC

The warp yarns and weft yarns of plain woven fabric which, being the principal axes of material of fabric, are orthogonal in the original configuration, but are obliquely crossed in the deformed configuration in general. The orthotropic constitutive model is unsuitable for fabric. In the oblique principal axes system the relations between loaded stress vectors and stress tensor are investigated, the stress fields of micro-weaving structures of fabric due to pure shear are carefully studied and, finally, a nonlinear micro-mechanical model for plain woven fabric is proposed. This model can accurately describe the nonlinear mechanical behavior of fabric observed in experiments. Under the assumption of small deformation and linearity of mechanical properties of fabric the model will degenerate into the existing linear model.     

Dynamic rotation and stretch tensors from a dynamic polar decomposition

The local rigid-body component of continuum deformation is typically characterized by the rotation tensor, obtained from the polar decomposition of the deformation gradient. Beyond its well-known merits, the polar rotation tensor also has a lesser known dynamical inconsistency: it does not satisfy the fundamental superposition principle of rigid-body rotations over adjacent time intervals. As a consequence, the polar rotation diverts from the observed mean material rotation of fibers in fluids, and introduces a purely kinematic memory effect into computed material rotation. Here we derive a generalized polar decomposition for linear processes that yields a unique, dynamically consistent rotation component, the dynamic rotation tensor, for the deformation gradient. The left dynamic stretch tensor is objective, and shares the principal strain values and axes with its classic polar counterpart. Unlike its classic polar counterpart, however, the dynamic stretch tensor evolves in time without spin. The dynamic rotation tensor further decomposes into a spatially constant mean rotation tensor and a dynamically consistent relative rotation tensor that is objective for planar deformations. We also obtain simple expressions for dynamic analogues of Cauchy's mean rotation angle that characterize a deforming body objectively.     

A kinetic interpretation of the section‐averaged Saint‐Venant system for natural river hydraulics

The classical Saint‐Venant system is well suited for the modeling of dam breaks, hydraulic jumps, reservoir emptying, flooding etc. For many applications, the extension of the Saint‐Venant system to the case of non‐rectangular channels is necessary and this section‐averaged Saint‐Venant system exhibits additional source terms. The main difficulty of these equations consists of the discretization of these source terms. In this paper we propose a kinetic interpretation for the section averaged Saint‐Venant system and derive an associated numerical scheme. The numerical scheme—2nd order in space and time—preserves the positivity of the water height, and is well‐balanced. Numerical results including comparisons with analytic and experimental test problems illustrate the accuracy and the robustness of the numerical algorithm. Copyright © 2010 John Wiley & Sons, Ltd.     

INCREMENTAL MICRO-MECHANICAL MODEL OF PLAIN WOVEN FABRIC

I. INTRODUCTION Previous research on woven fabric ignored the micro-weaving structures in fabric and modeled thewoven fabric as a ?exible orthotropic plate[1??3]. However, such a model is unable to predict certainbuckling phenomena of fabric commonly ob…