A Natural Generalization of Hypoelasticity and Eulerian Rate Type Formulation of Hyperelasticity

According to the classical hypoelasticity theory, the hypoelasticity tensor, i.e. the fourth order Eulerian constitutive tensor, characterizing the linear relationship between the stretching and an objective stress rate, is dependent on the current stress and must be isotropic. Although the classical hypoelasticity in this sense includes as a particular case the isotropic elasticity, it fails to incorporate any given type of anisotropic elasticity. This implies that one can formulate the isotropic elasticity as an integrable-exactly classical hypoelastic relation, whereas one can in no way do the same for any given type of anisotropic elasticity. A generalization of classical theory is available, which assumes that the material time derivative of the rotated stress is dependent on the rotated Cauchy stress, the rotated stretching and a Lagrangean spin, linear and of the first degree in the latter two. As compared with the original idea of classical hypoelasticity, perhaps the just-mentioned generalization might be somewhat drastic. In this article, we show that, merely replacing the isotropy property of the aforementioned stress-dependent hypoelasticity tensor with the invariance property of the latter under an R-rotating material symmetry group R⋆ G ₀, one may establish a natural generalization of classical theory, which includes all of elasticity. Here R is the rotation tensor in the polar decomposition of the deformation gradient and G ₀ any given initial material symmetry group. In particular, the classical case is recovered whenever the material symmetry is assumed to be isotropic. With the new generalization it is demonstrated that any two non-integrable hypoelastic relations based on any two objective stress rates predict quite different path-dependent responses in nature and hence can in no sense be equivalent. Thus, the non-integrable hypoelastic relations based on any given objective stress rate constitute an independent constitutive class in its own right which is disjoint with and hence distinguishes itself from any class based on another objective stress rate. Only for elasticity, equivalent hypoelastic formulations based on different stress rates may be established. Moreover, universal integrability conditions are derived for all kinds of objective corotational stress rates and for all types of material symmetry. Explicit, simple, integrable-exactly hypoelastic relations based on the newly discovered logarithmic stress rate are presented to characterize hyperelasticity with any given type of material symmetry. It is shown that, to achieve the latter goal, the logarithmic stress rate is the only choice among all infinitely many objective corotational stress rates. This revised version was published online in August 2006 with corrections to the Cover Date.     

On the types and number of plane waves in hypoelastic materials

General principles are formulated for modeling the elastic deformation of materials and analyzing plane waves in nonlinearly elastic materials such as hyperelastic, hypoelastic, and those governed by the general law of elasticity. The results of studying the propagation of plane waves in hypoelastic materials are further outlined. The influence of initial stresses and initial velocities on the types and number of plane waves is studied. Wave effects characteristic of hypoelastic materials are predicted theoretically. One of such effects is blocking of certain types of plane waves by initial stresses __________ Translated from Prikladnaya Mekhanika, Vol. 41, No. 11, pp. 96–107, November 2005.     

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.     

Characteristic cones of the equations in the nonlinear theory of elasticity

The roots of the equation for the characteristic normals for two systems of differential equations in the nonlinear theory of elasticity are investigated. The first model is constructed using a thermodynamic identity. The second is a very simple hypoelastic model (the deviator of the stress-rate tensor is proportional to the deviator of the strain-rate tensor). It is shown that the roots of the equations for the normals to the characteristics for the second model are the same as the first-order terms in the expansion of the roots of the first model with respect to the strain-tensor deviator.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 126–132, May–June, 1974.The author is grateful to S. K. Godunov for discussions.     

Homogenization problems of a hollow cylinder made of elastic materials with discontinuous properties

In this paper, we present the homogenization of an anisotropic hollow layered tube with discontinuous elastic coefficients. We focus on some aspects of technological importance, such as the effective coefficients of anisotropic materials, the behavior of the homogenized displacements and stresses, the discontinuities of in-plane shear, hoop and longitudinal stresses, the homogenization-induced anisotropy in the isotropic case. We conclude that the problem of cylindrically anisotropic tubes under extension, torsion, shearing and pressuring is stable by homogenization and we define the effective tensor of the material elastic coefficients. Some numerical examples confirm the theoretical results.     

Modeling cylindrical waves in nonlinear elastic composites

A procedure of deriving nonlinear wave equations that describe the propagation and interaction of hyperelastic cylindrical waves in composite materials modeled by a mixture with two elastic constituents is outlined. Nonlinearity is introduced by metric coefficients, Cauchy-Green strain tensor, and Murnaghan potential. It is the quadratic nonlinearity of all governing relations. For a configuration (state) dependent on the radial coordinate and independent of the angular and axial coordinates, quadratically nonlinear wave equations for stresses are derived and a relationship between the components of the stress tensor and partial strain gradient is established. Four combinations of physical and geometrical nonlinearities in systems of wave equations are examined. Nonlinear wave equations are explicitly written for three of the combinations __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 6, pp. 63–72, June 2007.     

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.     

Non-oscillation Criteria for Hypoelastic Models under Simple Shear Deformation

Ten objective rates, spinning or non-spinning, are critically examined from the viewpoint of Sturm's theorems in ordinary differential equations. Upon developing implication relations of oscillatory, non-oscillatory, and disconjugate behavior, we establish oscillation and non-oscillation criteria which pick out the objective stress rates that lead to oscillatory and non-oscillatory responses in simple shear deformation, respectively. Among the hypoelastic equations associated with the spinning objective rates examined, the Jaumann equation is an oscillatory minorant, the homogeneous Xiao–Bruhns–Meyers equation is a non-oscillatory majorant, and the homogeneous Green–Naghdi equation is a disconjugate majorant. If (Sturm comparable) non-spinning objective rates are also taken into consideration, then the Durban–Baruch equation becomes an oscillatory minorant, but the other two equations remain to play the same roles. The Jaumann equation is a Sturm majorant for all the other nine homogeneous hypoelastic equations, and the homogeneous Szabó–Balla-2 equation is a Sturm minorant for all the other nine homogeneous hypoelastic equations. Most of the solutions of the zeroth-grade hypoelastic equations at simple shear have already been published, except for those of Szabó and Balla, to which the closed-form solutions are derived here. Moreover, all solutions are extended to include the effect of initial stresses. This revised version was published online in August 2006 with corrections to the Cover Date.     

New Statement of the elasticity problem for a layer

The problem on the equilibrium of an inhomogeneous anisotropic elastic layer is considered. The classical statement of the problem in displacements consists of three partial differential equations with variable coefficients for the three displacements and of three boundary conditions posed at each point of the boundary surface. Sometimes, instead of the statement in displacements, it is convenient to use the classical statement of the problem in stresses [1] or the new statement of the problem in stresses proposed by B. E. Pobedrya [2]. In the case of the problem in stresses, it is necessary to find six components of the stress tensor, which are functions of three coordinates. The choice of the statement of the problem depends on the researcher and, of course, on the specific problem. The fact that there are several statements of the problem makes for a wider choice of the method for solving the problem. In the present paper, for a layer with plane boundary surfaces, we propose a new statement of the problem, which, in contrast to the other two statements indicated above, can be called a mixed statement. The problem for a layer in the new statement consists of a system of three partial differential equations for the three components of the displacement vector of the midplane points. The system is coupled with three integro-differential equations for the three longitudinal components of the stress tensor. Thus, in the new statement, just as in the other statements in stresses, one should find six functions. In the new statement, three of these functions (the displacements of the midplane points) are functions of two coordinates, and the other three functions (the longitudinal components of the stress tensor) are functions of three coordinates. It is shown that all equations in the new statement are the Euler equations for the Reissner functional with additional constraints. After the problem is solved in the new statement, three components of the displacement vector and three transverse components of the stress tensor are determined at each point of the layer. The new statement of the problem can be used to construct various engineering theories of plates made of composite materials.     

The hydrodynamics of convective stellar envelopes

Large-scale movements in convective stellar envelopes are studied by the methods of hydrodynamics, taking account of the compressibility of the gas, rotation, gravitation, and anisotropic turbulent viscosity. From dimensional analysis and observational data, the order of the nondimensional parameters controlling the motion of gaseous masses in the convective zone can be estimated. Estimates of the dimensionless parameters are used to simplify the equations of motion and develop an approximate method of solution. It is assumed that the motion of the gas in the convective envelope is anisotropic and turbulent and that the Reynolds stress tensor is linear in the rates of deformation, while the coefficients of proportionality in a spherical coordinate system form a diagonal tensor. It is assumed that the boundaries of the convective zone are free of turbulent stresses. Three approximations are found for the stream function for meridional motions, and the angular velocity of differential rotation for the case when the forces of turbulent viscosity are very strong. The components of the tensor of turbulent viscosity are assumed to be known functions of the depth of the convective zone.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 27–38, March–April, 1971.     

On the bifurcation and postbifurcation analysis of elastic-plastic solids under general prebifurcation conditions

In this work we have studied the bifurcation and postbifurcation of elastic-plastic solids whose behavior near the critical point could not be idealized as hypoelastic and thus the “hypoelastic comparison solid” concept of R. Hill's theory is no longer applicable. First a simple continuous model is considered in order to illustrate the different possibilities in the stability behavior of the structures considered here. Next, a general three-dimensional stability analysis for a broad class of rate independent elastic-plastic solids is presented. It is found that for all the constitutive theories considered and for all possible prebifurcation solutions, the bifurcation functional is a simple generalization of Hill's. A completely different postbifurcation analysis is needed, however, in the case where the “hypoelastic comparison solid” concept cannot be used.     

Dislocations in second strain gradient elasticity

A second strain gradient elasticity theory is proposed based on first and second gradients of the strain tensor. Such a theory is an extension of first strain gradient elasticity with double stresses. In particular, the strain energy depends on the strain tensor and on the first and second gradient terms of it. Using a simplified but straightforward version of this gradient theory, we can connect it with a static version of Eringen’s nonlocal elasticity. For the first time, it is used to study a screw dislocation and an edge dislocation in second strain gradient elasticity. By means of this second gradient theory it is possible to eliminate both strain and stress singularities. Another important result is that we obtain nonsingular expressions for the force stresses, double stresses and triple stresses produced by a straight screw dislocation and a straight edge dislocation. The components of the force stresses and of the triple stresses have maximum values near the dislocation line and are zero there. On the other hand, the double stresses have maximum values at the dislocation line. The main feature is that it is possible to eliminate all unphysical singularities of physical fields, e.g., dislocation density tensor and elastic bend-twist tensor which are still singular in the first strain gradient elasticity.     

Near-Wall Modification of an Explicit Algebraic Reynolds Stress Model Using Elliptic Blending

The elliptic blending approach is used in order to modify an Explicit Algebraic Reynolds Stress Model so as to reproduce the correct near wall behaviour of the turbulent stresses. The anisotropy stress tensor is expressed as a linear combination of tensor bases whose coefficients are sensitised to the non-local wall-blocking effect through the elliptic blending parameter γ. This parameter is obtained from a separate elliptic equation. The model does not use the distance from the wall thus it can be easily applied to complex geometries. It is validated against detailed DNS data for mean and turbulence quantities for the case of flow and heat transfer between parallel flat plates at three Reynolds numbers as well as against experimental data for the flow in a backward facing step at Re _H = 28,000. The comparison with DNS results or experiments is quite satisfactory and shows the validity of the approach.     

各向同性率无关材料本构关系的不变性表示

在内变量理论的框架下,针对各向同性率无关材料,使用张量函数表示理论建立了塑性应变全量及增量本构关系的最一般的张量不变性表示. 它们均由3个完备不可约的基张量组合构成,这3个基张量分别是应力的零次幂、一次幂和二次幂. 因此得出,塑性应变、塑性应变增量与应力三者共主轴. 通过对基张量的正交化,给出了本构关系式在主应力空间中的几何解释. 进一步,全量(或增量)本构关系中3个组合因子被表达为应力、塑性应变(或塑性应变增量)的不变量的函数. 当塑性应变(或塑性应变增量)的3个不变量之间满足一定关系时,所给出的本构关系将退化为经典的形变理论(或塑性势理论).最后,还讨论它与奇异屈服面理论的关系,当满足一定条件时,两者是一致的.      

Laws of flow with a limiting gradient in anisotropic porous media

The equations of viscoplastic fluid flow through a porous medium are written for all types of anisotropy. It is shown that in anisotropic media the flows with a limiting gradient are characterized by two material tensors: the tensor of permeability (flow resistance) coefficients and the tensor of limiting gradients. A complex of laboratory measurements for determining the tensors of permeability coefficients and limiting gradients is considered for all types of anisotropic media. It is shown that the tensors of permeability coefficients and limiting gradients are coaxial. Conditions of flow onset and fluid flow laws are formulated for media with monoclinic and triclinic symmetries of flow characteristics.     

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.     

Three-dimensional phenomenological thermodynamic model of pseudoelasticity of shape memory alloys at finite strains

The familiar small strain thermodynamic 3D theory of isotropic pseudoelasticity proposed by Raniecki and Lexcellent is generalized to account for geometrical effects. The Mandel concept of mobile isoclinic, natural reference configurations is used in order to accomplish multiplicative decomposition of total deformation gradient into elastic and phase transformation (p.t.) parts, and resulting from it the additive decomposition of Eulerian strain rate tensor. The hypoelastic rate relations of elasticity involving elastic strain rate are derived consistent with hyperelastic relations resulting from free energy potential. It is shown that use of Jaumann corotational rate of stress tensor in rate constitutive equations formulation proves to be convenient. The formal equation for p.t. strain rate , describing p.t. deformation effects is proposed, based on experimental evidence. Phase transformation kinetics relations are presented in objective form. The field, coupled problem of thermomechanics is specified in rate weak form (rate principle of virtual work, and rate principle of heat transport). It is shown how information on the material behavior and motion inseparably enters the rate virtual work principle through the familiar bridging equation involving Eulerian rate of nominal stress tensor.

---

     

Effect of Allowance for Small‐Scale Perturbations of the Carrier Phase on Properties of the System of Equations of a Two‐Fluid Flow with Incompressible Phases

The effect of the fluctuating components of kinetic energy and stress tensor of the carrier phase, which were previously obtained by the cell technique, on the properties of the system of equations of a gas–liquid flow with incompressible phases is considered. It is shown that the characteristic properties of this system and also the possibility of modeling the Zuber–Findlay empirical relation are determined by the tensor of fluctuating stresses of the carrier phase.