Elastoplastic deformation of elements of an isotropic solid along paths of small curvature: Constitutive equations incorporating the stress mode

Constitutive equations relating the components of the stress tensor in a Eulerian coordinate system and the linear components of the finite-strain tensor are derived. These stress and strain measures are energy-consistent. It is assumed that the stress deviator is coaxial with the plastic-strain differential deviator and that the first invariants of the stress and strain tensors are in a nonlinear relationship. In the case of combined elastoplastic deformation of elements of the body, this relationship, as well as the relationship between the second invariants of the stress and strain deviators, is determined from fundamental tests on a tubular specimen subjected to proportional loading at several values of stress mode angle (the third invariant of the stress deviator). Methods to individualize these relationships are proposed. The initial assumptions are experimentally validated. The constitutive equations derived underlie an algorithm for solving boundary-value problems __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 6, pp. 43–55, June 2007.     

The evolution of Hooke's law due to texture development in FCC polycrystals

Initially isotropic aggregates of crystalline grains show a texture-induced anisotropy of both their inelastic and elastic behavior when submitted to large inelastic deformations. The latter, however, is normally neglected, although experiments as well as numerical simulations clearly show a strong alteration of the elastic properties for certain materials. The main purpose of the work is to formulate a phenomenological model for the evolution of the elastic properties of cubic crystal aggregates. The effective elastic properties are determined by orientation averages of the local elasticity tensors. Arithmetic, geometric, and harmonic averages are compared. It can be shown that for cubic crystal aggregates all of these averages depend on the same irreducible fourth-order tensor, which represents the purely anisotropic portion of the effective elasticity tensor. Coupled equations for the flow rule and the evolution of the anisotropic part of the elasticity tensor are formulated. The flow rule is based on an anisotropic norm of the stress deviator defined by means of the elastic anisotropy. In the evolution equation for the anisotropic part of the elasticity tensor the direction of the rate of change depends only on the inelastic rate of deformation. The evolution equation is derived according to the theory of isotropic tensor functions. The transition from an elastically isotropic initial state to a (path-dependent) final anisotropic state is discussed for polycrystalline copper. The predictions of the model are compared with micro–macro simulations based on the Taylor–Lin model and experimental data.     

Tresca-type plastic materials in the theory of hypoelasticity II. Optical constitutive equations and birefringence in simple shear

Behavior of a Tresca type plastic dielectric is investigated theoretically from a continuum mechanical point of view. The optical constitutive equations are defined as special cases of a hypo-elastic dielectric of grade two. The singularity condition of the constitutive equations satisfies the Tresca yield criterion. The index deviator tensor is proportional to the stress deviator tensor and, then, the birefringence and the extinction angle are expressed by the stress deviator. Their numerical variations with the angle of shear in simple shear deformation are shown.     

Theory of 4D-media with stationary dislocations

In earlier studies, the authors showed that an application of classical methods of mechanics of deformable media to the study of properties of 4D-space-time continuum permit stating consistent models of nonholonomic media mechanics consistent with the first and second laws of thermodynamics. In the present paper, we show that the classical methods of continuum mechanics are also promising when modeling physical processes. It is shown that, just as in the three-dimensional theory of stationary dislocations, there exist dislocations of three types for a generalized 4D-medium. They correspond to the decomposition of the free distortion tensor into a spherical tensor, a deviator tensor, and a pseudotensor of rotations. We interpret several particular models, thus showing that the proposed model describes the spectrum of known physical interactions: electromagnetic, strong, weak, and gravitational. We show that the resolving equations include the Maxwell equations of electrodynamics and the Yukawa equations for strong interactions as subsystems.     

Determining the axisymmetric elastoplastic state of thin shells with allowance for the third invariant of the stress deviator

A technique to determine the axisymmetric elastoplastic state of thin shells with allowance for the third invariant of the stress deviator is developed. The technique is based on the theory of thin shells that takes into account transverse shear and torsional strains. Plastic equations that relate the components of the stress tensor in Eulerian coordinates with the linear components of the finite-strain tensor are used as constitutive equations. The nonlinear scalar functions in the constitutive equations are found from base tests on tubular specimens under proportional loading for different stress modes. The boundary-value problem is solved by numerically integrating a system of ordinary differential equations     

Identification of symmetry type of linear elastic stiffness tensor in an arbitrarily orientated coordinate system

We develop a method through the mirror plane (MP) to identify the symmetry type of linear elastic stiffness tensor whose components are given with respect to an arbitrarily oriented coordinate system. The method is based on the irreducible decomposition of high-order tensor into a set of deviators and the multipole representation of a deviator into a scalar and a unit-vector set. Since a unit-vector depends on two Euler angles, we can illustrate the MP normals of the elastic tensor as zeros of a characteristic function on a unit disk and identify its symmetry immediately, which is clearer and simpler than the methods proposed before. Furthermore, by finding the common MPs of three unit-vector sets using Fortran recipes, we can also analytically recognize the symmetry type first and then recover the natural coordinate system associated with the linear elastic tensor. The structures of linear elastic stiffness tensors of real materials with all possible anisotropies are investigated in detail.     

On a representation of the solution of the linear elasticity equations

A new representation of the stress tensor in the linear theory of elasticity is proposed. The representation satisfies the equilibrium equations and the compatibility conditions for strains. In this representation, the stress tensor is expressed in terms of a harmonic vector. The second boundary-value problem for an elastic half-space and elastic layer is considered as an example.Translated from Prikladnaya Mekhanika, Vol. 40, No. 11, pp. 85–91, November 2004.This revised version was published online in April 2005 with a corrected cover date.     

Constitutive equations for describing the elastoplastic deformation of elements of a body along small-curvature paths in view of the stress mode

Equations relating the components of the stress and strain tensors (constitutive equations) are formulated in terms of Euler coordinates. The equations describe the finite elastoplastic deformation of an isotropic body along paths of small curvature. It is assumed that the stress deviator is coaxial with the plastic-strain differential deviator. The relationships between the first and second invariants of the stress and strain tensors in the case of complex elastoplastic deformation of the body’s elements are determined from base tests on tubular specimens loaded along rectilinear paths for several values of the stress mode angle. Methods for specification of these relationships are proposed. The assumptions adopted to derive the constitutive equations are validated experimentally __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 4, pp. 62–72, April 2006.     

Conditional averaging of the equations of flow in random composite porous media

The problem of conditional averaging of the transport equations is solved for a neutral impurity in a composite medium with random porosity and impurity diffusion tensor. An unclosed system of conditionally averaged equations is constructed and closed using the globally averaged equations. The average impurity concentration fields for the individual phases of the composite medium and the phase-continuum interaction characteristics are calculated.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No, 1, pp. 75–81, January–February, 1987.     

Nonlinear analysis of beam–column joints using softened truss model

The aim of this study is to expand the application of the nonlinear softened truss model for membrane elements on beam–column joints. The softened truss model employs three equations for equilibrium, three for compatibility and four equations for the constitutive laws of materials. The constitutive equations for both the concrete and steel are based on the actually observed stress–strain relationships. The model has three important attributes. The first is the nonlinear association of stress and strain. The second, and conceivably more noteworthy, is the softening of concrete in compression due to tensile strains in the perpendicular direction. The third is that the influence of the concrete tensile stresses between cracks on the average stress–strain relationship for reinforcing steel and the influence of orthogonal tensile stresses on the compression stress–strain relationship for concrete can be considered in the model. For beam–column joints, one of the most important factors influencing the behaviour is certainly the bond conditions of the beam bars. In this study, the softened truss model is expanded to take into account the influence of this important factor into account. In the revised version of the model, full strain compatibility does not exist between the steel reinforcement and the surrounding concrete and thus the factors influencing the bond-slip between concrete and reinforcement is adequately considered. The improved softened truss model is applied on 51 exterior beam–column joint tests. It is apparent from the results that the revised model gives very accurate predictions of the shear strength of joints and is an improvement on the existing version of the model proposed by Hsu.     

On the stability of the flow of a viscoelastic fluid down an inclined plane

The paper presents an analysis of laminar flow of a film of viscoelastic fluid flowing under gravity down an infinite inclined plane. It is assumed that the mechanical behavior of the fluid can be represented by a generalized Maxwell model, whose constitutive equation contains a time derivative of the deviator of the stress tensor in the Jaumann sense [1. 2]. The equations of motion of the viscoelastic fluid considered here admit an exact solution for the case of rectilinear laminar flow with a plane free boundary. The stability of this flow with respect to surface waves is investigated by the method of successive approximations described in [3, 4].     

Construction of Dynamically Equivalent Time-Invariant Forms for Time-Periodic Systems

In this study dynamically equivalent time-invariant forms are obtained for linear and non-linear systems with periodically varying coefficients via Lyapunov–Floquet (L–F) transformation. These forms are equivalent in the sense that the local stability and bifurcation characteristics are identical for both systems in the entire parameter space. It is well known that the L–F transformation converts a linear periodic first order system into a time-invariant one. In the first part of this study a set of linear second order periodic equations is converted into an equivalent set of time-independent second order equations through a sequence of linear transformations. Then the transformations are applied to a time-periodic quadratic Hamiltonian to obtain its equivalent time-invariant form. In the second part, time-invariant forms of nonlinear equations are studied. The application of L–F transformation to a quasi-linear periodic equation converts the linear part to a time-invariant form and leaves the non-linear part with time-periodic coefficients. Dynamically equivalent time-invariant forms are obtained via time-periodic center manifold reduction and time-dependent normal form theory. Such forms are constructed for general hyperbolic systems and for some simple critical cases, including that of one zero eigenvalue and a purely imaginary pair. As a physical example of these techniques, a single and a double inverted pendulum subjected to periodic parametric excitation are considered. The results thus obtained are verified by numerical simulation.     

Large eddy simulation of natural convection along a vertical isothermal surface

Large eddy simulations of natural convection along a vertical isothermal surface have been carried out using a parallel CFD code SMAFS (Smoke Movement And Flame Spread) developed by the first author to study the dynamics of the natural convection flow and the associated convective heat transfer, with sub-grid scale turbulence modeled using the Smagorinsky model. In the computation, the filtered governing equations are discretized using finite volume method, with the variables at the cell faces in the finite volume discrete equations approximated by a second order bounded QUICK scheme and the diffusion term computed based on central difference scheme. The computation was time marched explicitly, with momentum equations solved based on a second order fractional-step Adams–Bashford scheme and enthalpy computed using a second order Runge–Kutta scheme. The Poisson equation for pressure from the continuity equation was solved using a multi-grid solver. The results including the temperature and velocity profiles of the boundary layer and the local heat transfer rate are analyzed. Comparison is made with experimental data and good agreement is found.     

Some families of exact solutions of the equations of two-dimensional shallow water theory

Results are presented of a study of the exact solutions of the equations of two-dimensional unsteady and steady shallow water theory, based on the group properties of these equations. The first part presents the group properties of the equations in question; the second part presents the invariant solutions of these equations.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, Vol. 10, No. 6, pp. 62–71, November–December, 1969.The author wishes to thank L. V. Ovsyannikov and N. Kh. Ibragimov for valuable guidance in carrying out this study.     

Dynamic compactness of coupled thermoelastic waves in continuously nonuniform anisotropic media

In the present work, the dynamic problem of coupled thermoelasticity with the most general type of nonuniformity and anisotropy is analyzed. The hyperbolic nature of the system of equations of coupled thermoelasticity is demonstrated, effects of extinction of separate waves by superposition of elastic and thermoelastic wave fronts are investigated, and the interrelationship of different orders of discontinuity of stresses, displacements, and temperature is determined. The case of the uncoupled problem of thermoelasticity is especially analyzed. Sufficient conditions are obtained for the dynamic density for wave processes in thermoelasticity, previously investigated for boundary value problems of hyperbolic systems of second order differential equations [1], andelastic stress waves [2] are obtained. The generally accepted system of tensor notation for the theory of thermoelasticity is used [3].Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 154–163, May–June, 1981.     

Turbulent air flow over a moving curved surface

A numerical model of turbulent air flow over a curved surface is described. The model is based on two-dimensional nonlinear Reynolds equations and continuity equations written in a coordinate system moving with the profile of the curved surface. The Reynolds stresses are represented in the form of the product of the isotropic turbulent viscosity coefficient, which increases linearly with height, and the deformation tensor of the mean velocity field. Flow over a stationary sinusoidal surface and a sinusoidal gravity wave on water is simulated. The structure of the velocity and pressure wave fields is obtained. The differences in flow over stationary and moving surfaces are analyzed.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 20–24, September–October, 1986.     

Equations of laminar boundary layer in a two-phase medium

The motion of a two-phase medium in which the carrier component has low viscosity is considered. The equations obtained in [1], to which the viscous stress tensor in the fluid is added, are used. The boundary layer method [2] makes it possible to obtain asymptotic equations for the wall region. These equations have different forms depending on the characteristic values of the dimensionless determining parameters.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 51–60, January–February, 1979.I thank A. N. Kraiko for discussing the work.     

A new two-point deformation tensor and its relation to the classical kinematical framework and the stress concept

Starting from the issue of what is the correct form for a Legendre transformation of the strain energy in terms of Eulerian and two-point tensor variables we introduce a new two-point deformation tensor, namely H=(F−F^−T)/2, as a possible deformation measure involving points in two distinct configurations. The Lie derivative of H is work conjugate to the first Piola–Kirchhoff stress tensor P. The deformation measure H leads to straightforward manipulations within a two-point setting such as the derivation of the virtual work equation and its linearization required for finite element implementation. The manipulations are analogous to those used for the Lagrangian and Eulerian frameworks. It is also shown that the Legendre transformation in terms of two-point tensors and spatial tensors require Lie derivatives. As an illustrative example we propose a simple Saint Venant–Kirchhoff type of a strain-energy function in terms of H. The constitutive model leads to physically meaningful results also for the large compressive strain domain, which is not the case for the classical Saint Venant–Kirchhoff material.