On the yield functions in Lagrangian stress space of elastic–plastic materials with internal constraints

In the present paper it is shown that the elastic range in the second Piola–Kirchhoff stress space can be chosen in a hyperplane which is through the origin of Lagrangian stress space and perpendicular to the normal of the constraint manifold at the plastic configuration, if the determinate stress response of the elastic–plastic material with simple internal constraints with some condition is correctly chosen, otherwise, it is in general in a hypersurface and the normal flow rule by Il yushin’s postulate will have an indeterminate part. The choice of determinate stress response is probable because of its indeterminacy. Therefore the yield function should be a function of the second Piola–Kirchhoff stress lying in the hyperplane so that it is more simple and the back stress as the geometric center of the elastic range in general is inside the elastic range. Finally some examples are concerned. The project supported by the National Natural Science Foundation of China (10272055).     

Well-Posedness of Surface Wave Equations Above a Viscoelastic Fluid

In this article, we prove the local well-posedness, for arbitrary initial data with certain regularity assumptions, of the equations of a Viscoelastic Fluid of Johnson–Segalman type in a domain with a free surface. Managing more general constitutive laws is also briefly depicted. The 2D geometry is defined by a solid fixed bottom and an upper free boundary submitted to surface tension. The proof relies on a Lagrangian formulation. First we solve two intermediate problems through a fixed point using mainly (Allain in Appl Math Optim 16:37–50, 1987) for the Navier–Stokes part. Then we solve the whole Lagrangian problem on [0, T ₀] for T ₀ small enough through a contraction mapping. Since the Lagrangian solution is regular enough and the change of coordinates invertible, we can come back to an Eulerian one.     

Objective Tensor Rates and Applications in Formulation of Hyperelastic Relations

Following Ogden, a class of objective (Lagrangian and Eulerian) tensors is identified among the second-rank tensors characterizing continuum deformation, but a more general definition of objectivity than that used by Ogden is introduced. Time rates of tensors are determined using convective rates. Sufficient conditions of objectivity are obtained for convective rates of objective tensors. Objective convective rates of strain tensors are used to introduce pairs of symmetric stress and strain tensors conjugate in a generalized sense. The classical definitions of conjugate Lagrangian (after Hill) and Eulerian (after Xiao et al.) stress and strain tensors are particular cases of the definition of conjugacy of stress and strain tensors in the generalized sense used in the present paper. Pairs of objective stress and strain tensors conjugate in the generalized sense are used to formulate constitutive relations for a hyperelastic medium. A family of objective generalized strain tensors is introduced, which is broader than Hill’s family of strain tensors. The basic forms of the hyperelastic constitutive relations are obtained with the aid of pairs of Lagrangian stress and strain tensors conjugate after Hill (the strain tensors in these pairs belong to the family of generalized strain tensors). A method is presented for generating reduced forms of the constitutive relations with the aid of pairs of Lagrangian and Eulerian stress and strain tensors conjugate in the generalized sense which are obtained from pairs of Lagrangian tensors conjugate after Hill by mapping tensor fields on one configuration of a deformable body to tensor fields on another configuration.      

Geodesics in the Space of Measure-Preserving Maps and Plans

We study Brenier’s variational models for incompressible Euler equations. These models give rise to a relaxation of the Arnold distance in the space of measure-preserving maps and, more generally, measure-preserving plans. We analyze the properties of the relaxed distance, we show a close link between the Lagrangian and the Eulerian model, and we derive necessary and sufficient optimality conditions for minimizers. These conditions take into account a modified Lagrangian induced by the pressure field. Moreover, adapting some ideas of Shnirelman, we show that, even for non-deterministic final conditions, generalized flows can be approximated in energy by flows associated to measure-preserving maps.     

Forms of null Lagrangians in field theories of continuum mechanics

The divergence representation of a null Lagrangian that is regular in a star-shaped domain is used to obtain its general expression containing field gradients of order ≤ 1 in the case of spacetime of arbitrary dimension. It is shown that for a static three-component field in the three-dimensional space, a null Lagrangian can contain up to 15 independent elements in total. The general form of a null Lagrangian in the four-dimensional Minkowski spacetime is obtained (the number of physical field variables is assumed arbitrary). A complete theory of the null Lagrangian for the n-dimensional spacetime manifold (including the four-dimensional Minkowski spacetime as a special case) is given. Null Lagrangians are then used as a basis for solving an important variational problem of an integrating factor. This problem involves searching for factors that depend on the spacetime variables, field variables, and their gradients and, for a given system of partial differential equations, ensure the equality between the scalar product of a vector multiplier by the system vector and some divergence expression for arbitrary field variables and, hence, allow one to formulate a divergence conservation law on solutions to the system.     

Influence of elastic material compressibility on parameters of an expanding spherical stress wave

We investigated the influence of elastic material compressibility on parameters of an expanding spherical stress wave. The material compressibility is represented by Poisson’s ratio, ν, in this paper. The stress wave is generated by a pressure produced inside a spherical cavity surrounded by the isotropic elastic material. The analytical closed form formulae determining the dynamic state of the mechanical parameters (displacement, particle velocity, strains, stresses, and material density) in the material have been derived. These formulae were obtained for surge pressure p(t) = p ₀ = const inside the cavity. From analysis of these formulae, it is shown that the Poisson’s ratio substantially influences the course of material parameters in space and time. All parameters intensively decrease in space together with an increase of the Lagrangian coordinate, r. On the contrary, these parameters oscillate versus time around their static values. These oscillations decay in the course of time. We can mark out two ranges of parameter ν values in which vibrations of the parameters are “damped” at a different rate. Thus, Poisson’s ratio in the range below about 0.4 causes intense decay of parameter oscillations. On the other hand in the range 0.4 < ν < 0.5, i.e. in quasi-incompressible materials, the “damping” of parameter vibrations is very low. In the limiting case when ν = 0.5, i.e. in the incompressible material, “damping” vanishes, and the parameters harmonically oscillate around their static values. The abnormal behaviour of the material occurs in the range 0.4 < ν < 0.5. In this case, an insignificant increase of Poisson’s ratio causes a considerable increase of the parameter vibration amplitude and decrease of vibration “damping”.      

Evaluation of a Modified Reynolds Stress Model for Turbulent Dispersed Two-Phase Flows Including Two-Way Coupling

A modified Reynolds stress turbulence model for the pressure rate of strain can be derived for dispersed two-phase flows taking into account gas-particle interaction. The transport equations for the Reynolds stresses as well as the equation for the fluctuating pressure can be derived starting from the multiphase Navier–Stokes equations. The unknown pressure rate of strain correlation in the Reynolds stress equations is then modelled by considering the multiphase equation for the fluctuating pressure. This leads to a multiphase pressure rate of strain model. The extra particle interaction source terms occurring in the model for the pressure rate of strain can be constructed easily, with no noticeable extra computational cost. Eulerian–Lagrangian simulation results of a turbulent dispersed two-phase jet are presented to show the differences in results with and without the new two-way coupling terms.     

Acoustic Lagrangian velocity measurement in a turbulent air jet

Measuring Lagrangian velocities in a turbulent flow is of a great interest for turbulence modeling. We report measurements made in an axisymmetric turbulent air jet at Reynolds number R _λ ≃ 320, using acoustical Doppler scattering. Helium-filled soap bubbles are used as Lagrangian tracers. We describe an experimental setup which allows the simultaneous measurement of the full three-component Lagrangian velocity and the longitudinal Eulerian one. Lagrangian velocity probability density functions (PDF) are found Gaussian, close to Eulerian ones. Velocity correlations are analysed as well as the statistical dependence between components.     

A discrete element model for simulating saturated granular soil

A numerical model is developed to simulate saturated granular soil, based on the discrete element method. Soil particles are represented by Lagrangian discrete elements, and pore fluid, by appropriate discrete elements which represent alternately Lagrangian mass of water and Eulerian volume of space. Macro-scale behavior of the model is verified by simulating undrained biaxial compression tests. Micro-scale behavior is compared to previous literature through pore pressure pattern visualization during shear tests. It is demonstrated that dynamic pore pressure patterns are generated by superposed stress waves. These pore-pressure patterns travel much faster than average drainage rate of the pore fluid and may initiate soil fabric change, ultimately leading to liquefaction in loose sands. Thus, this work demonstrates a tool to roughly link dynamic stress wave patterns to initiation of liquefaction phenomena.     

Possibilities and Limitations of Computer Simulations of Industrial Turbulent Dispersed Multiphase Flows

We give an overview on the usage of computer simulations in industrial turbulent dispersed multiphase flows. We present a few examples of industrial flows: bubble columns and bubbly pipe flows, stirred tanks, cyclones, and a fluid catalytic cracking unit. The fluid catalytic cracking unit is used to illustrate the complexity of the physical phenomena involved, and the possibilities and limitations of the different approaches used: Eulerian–Lagrangian (particle-tracking) and Eulerian–Eulerian (two-fluid). In the first approach, the continuous phase is solved using either RANS simulations (Reynolds-Averaged Navier–Stokes simulations) or DNS/LES (Direct Numerical Simulations/Large-Eddy Simulations), and the individual particles are tracked. In the second approach, the dispersed phase is averaged, leading to two sets equations, which are quite similar to the RANS equations of single-phase flows. The Eulerian–Eulerian approach is the most commonly used in industrial applications, however, it requires a significant amount of modelling. Eulerian–Lagrangian RANS can be simpler to use; in particular in situations involving complex boundary conditions, polydisperse flows and agglomeration/breakup. The key issue for the success of the simulations is to have good models for the complex physics involved. A major weakness is the lack of good models for: the turbulence modification promoted by the particles, the inter-particle interactions, and the near-wall effects. Eulerian–Lagrangian DNS/LES can play an important role as a research tool, in order to get a better physical understanding, and to improve the models used in the RANS simulations (either Eulerian–Eulerian or Eulerian–Lagrangian).     

Sediment flow in inclined vessels calculated using a multiphase particle-in-cell model for dense particle flows

Sedimentation of particles in an inclined vessel is predicted using a two-dimensional, incompressible, multiphase particle-in-cell (MP-PIC) method. The numerical technique solves the governing equations of the fluid phase using a continuum model and those of the particle phase using a Lagrangian model. Mapping particle properties to an Eulerian grid and then mapping back computed stress tensors to particle positions allows a complete solution of sedimentation from a dilute mixture to close-pack. The solution scheme allows for distributions of types, sizes and density of particles, with no numerical diffusion from the Lagrangian particle calculations. The MP-PIC solution method captures the physics of inclined sedimentation which includes the clarified fluid layer under the upper wall, a dense mixture layer above the bottom wall, and instabilities which produce waves at the clarified fluid and suspension interface. Measured and calculated sedimentation rates are in good agreement.     

Sediment flow in inclined vessels calculated using a multiphase particle-in-cell model for dense particle flows

Sedimentation of particles in an inclined vessel is predicted using a two-dimensional, incompressible, multiphase particle-in-cell (MP-PIC) method. The numerical technique solves the governing equations of the fluid phase using a continuum model and those of the particle phase using a Lagrangian model. Mapping particle properties to an Eulerian grid and then mapping back computed stress tensors to particle positions allows a complete solution of sedimentation from a dilute mixture to close-pack. The solution scheme allows for distributions of types, sizes and density of particles, with no numerical diffusion from the Lagrangian particle calculations. The MP-PIC solution method captures the physics of inclined sedimentation which includes the clarified fluid layer under the upper wall, a dense mixture layer above the bottom wall, and instabilities which produce waves at the clarified fluid and suspension interface. Measured and calculated sedimentation rates are in good agreement.     

The Anisotropic Lagrangian Averaged Euler and Navier-Stokes Equations

 The purpose of this paper is twofold. First, we give a derivation of the Lagrangian averaged Euler (LAE-α) and Navier-Stokes (LANS-α) equations. This theory involves a spatial scale α and the equations are designed to accurately capture the dynamics of the Euler and Navier-Stokes equations at length scales larger than α, while averaging the motion at scales smaller than α. The derivation involves an averaging procedure that combines ideas from both the material (Lagrangian) and spatial (Eulerian) viewpoints. This framework allows the use of a variant of G. I. Taylor's ``frozen turbulence' hypothesis as the foundation for the model equations; more precisely, the derivation is based on the strong physical assumption that fluctutations are frozen into the mean flow. In this article, we use this hypothesis to derive the averaged Lagrangian for the theory, and all the terms up to and including order α² are accounted for. The equations come in both an isotropic and anisotropic version. The anisotropic equations are a coupled system of PDEs (partial differential equations) for the mean velocity field and the Lagrangian covariance tensor. In earlier works by Foias, Holm & Titi [10], and ourselves [16], an analysis of the isotropic equations has been given. In the second part of this paper, we establish local in time well-posedness of the anisotropic LANS-α equations using quasilinear PDE type methods. (Accepted September 2, 2002) Published online November 26, 2002 Dedicated to Stuart Antman on the occasion of his 60th birthday Communicated by S. MüLLER     

Static equilibrium of an elastic orthotropic medium with an arbitrarily oriented triaxial ellipsoidal inclusion

The stress state of an elastic orthotropic medium with an arbitrarily oriented triaxial ellipsoidal inclusion is analyzed. A solution is obtained using the triple Fourier transform and the Fourier-transformed Green’s function for an infinite anisotropic medium. The high efficiency of the approach is demonstrated by solving the problem for a transversely isotropic material with a spheroidal cavity for which the exact solution is known. A numerical analysis is conducted to study the stress distribution over the surface of the inclusion with different orientations in the orthotropic space. It is revealed that in some cases the orientation of the inclusion has a strong effect on the stress concentration __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 4, pp. 55–61, April 2007.     

A measurement of three-dimensional strains around a creep-crack tip

A new method for measuring the strain around the high-temperature creep-crack tip is proposed in the present paper. A grid pattern was described in a space of about 30 μm with a diamond stylus in the area ahead of the precrack of the specimen. Then, the distortion of the grids and the change in specimen thickness which were induced by the creep deformation were measured by means of a photo-microscope and a roughness-measuring system, respectively. The three-dimensional components of the strain were calculated using the Lagrangian equation, into which the above measurement was introduced. The route of the creep-crack extension was examined in association with the local strain measured by the proposed method.     

Experimentally Evaluating Equilibrium Stress in Uniaxial Tests

Many models use the equilibrium stress, also sometimes known as the back stress, in characterizing the response of both polymeric and non-polymeric materials. We study the characteristics of the equilibrium and show that the tangent modulus and local Poisson’s ratio at equilibrium both are rate independent for common modeling assumptions. This fact is used to propose a method based on uniaxial tension or compression to measure the equilibrium stress, and the associated point’s tangent modulus and local Poisson’s ratio. The method is based on cyclic loading and identification of similar states with vastly different loading rates. The method is used to characterize the equilibrium stress in glassy polycarbonate, and the results are studied in regard to the possible error for such a measurement. The method is faster than most other proposed methods for calculating the equilibrium stress, and provides additional measurements of parameters at equilibrium that are normally not obtained.     

钢丝绳捻制成形接触问题的有限元分析

基于虚功原理建立了钢丝绳捻制成形过程中考虑接触摩擦时的刚度方程,将小球分裂算法用于接触搜索,Augmented Lagrange法用于计算接触力．以钢丝绳一次捻制成形过程为例,分析了摩擦系数和自扭转系数对接触应力和加工应力域变的影响．研究结果表明。摩擦系数对剪应力的影响比对等效应力和等效塑性应变的影响大;考虑摩擦时钢线表面的轴向残余应力计算结果比无摩擦条件下的计算结果更接近实验结果．