Duality in constitutive formulation of finite-strain elastoplasticity based on F=FeFp and F=FpFe decompositions

A constitutive theory for large elastic–plastic deformations is presented by employing F=F^pF^e decomposition of the total deformation gradient. A duality in constitutive formulation based on this and the well-known Lee's decomposition F=F_eF_p is established for isotropic polycrystalline and single crystal plasticity.     

A large deformation theory for rate-dependent elastic–plastic materials with combined isotropic and kinematic hardening

We have developed a large deformation viscoplasticity theory with combined isotropic and kinematic hardening based on the dual decompositions F=F^eF^p [Kröner, E., 1960. Allgemeine kontinuumstheorie der versetzungen und eigenspannungen. Archive for Rational Mechanics and Analysis 4, 273–334] and [Lion, A., 2000. Constitutive modelling in finite thermoviscoplasticity: a physical approach based on nonlinear rheological models. International Journal of Plasticity 16, 469–494]. The elastic distortion F^e contributes to a standard elastic free-energy ψ^(e), while , the energetic part of F^p, contributes to a defect energy ψ^(p) – these two additive contributions to the total free energy in turn lead to the standard Cauchy stress and a back-stress. Since F^e=FF^p-1 and , the evolution of the Cauchy stress and the back-stress in a deformation-driven problem is governed by evolution equations for F^p and – the two flow rules of the theory.We have also developed a simple, stable, semi-implicit time-integration procedure for the constitutive theory for implementation in displacement-based finite element programs. The procedure that we develop is “simple” in the sense that it only involves the solution of one non-linear equation, rather than a system of non-linear equations. We show that our time-integration procedure is stable for relatively large time steps, is first-order accurate, and is objective.     

Explicit expressions of S(v), H(v) and L(v) for anisotropic elastic materials

The three matrices L(v), S(v) and H(v), appearing frequently in the investigations of the two-dimensional steady state motions of elastic solids, are expressed explicitly in terms of the elastic stiffness for general anisotropic materials. The special cases of monoclinic materials with a plane of symmetry at x₃ = 0, x₁ = 0, and x₂ = 0 are all deduced. Results for orthotropic materials appearing in the literature may be recovered from the present explicit expressions.     

ON THE CRITICAL POINTS OF THE MAP Fp:X→‖AXB-C‖pp

Suppose A,B and C are the bounded linear operators on a Hilbert space H, when A has a generalized inverse A^- such that (AA^-)^*=AA^- and B has a generalized inverse B^- such that (B^-B)^*=B^-B,the general characteristic forms for the critical points of the map F_p:X^→‖AXB-C‖_p^p(1p=2. Similarly, the same question has been discussed for several operators.     

Polycrystalline kinematics: An extension of single crystal kinematics that incorporates initial microstructure

Single crystal F^eF^P kinematics are widely used as the basis for many crystal plasticity models. Within this kinematic framework, geometrically necessary dislocations (GNDs) initially do not exist and then they evolve as needed in the material. A shortcoming of this kinematic model is that there is no rigorous way to define the initial and evolving GND state in the same manner. By augmenting the single crystal F^eF^P kinematics with a geometric argument, a consistent methodology for determining the initial and evolving GND state has been derived. The augmented kinematics describe GND related microstructural features in the undeformed material like low angle sub-grain boundaries and high angle grain boundaries. Therefore these kinematics are particularly applicable to polycrystalline materials.     

On a formulation for anisotropic elastoplasticity at finite strains invariant with respect to the intermediate configuration

The paper is concerned with a formulation of anisotropic finite strain inelasticity based on the multiplicative decomposition of the deformation gradient F=F_eF_p. A major feature of the theory is its invariance with respect to rotations superimposed on the inelastic part of the deformation gradient. The paper motivates and shows how such an invariance can be achieved. At the heart of the formulation is the mixed-variant transformation of the structural tensor, defined as the tensor product of the privileged directions of the material as given in a reference configuration, under the action of F_p. Issues related to the plastic material spin are discussed in detail. It is shown that, in contrast to the isotropic case, any flow function formulated purely in terms of stress quantities, necessarily exhibits a non-vanishing plastic material spin. The possible construction of spin-free rates is discussed as well, where it is shown that the flow rule must then depend not only on the stress but on the strain as well.     

Dynamic Crack Growth Past a Stiff Inclusion: Optical Investigation of Inclusion Eccentricity and Inclusion-matrix Adhesion Strength

Interactions between a dynamically growing matrix crack and a stationary stiff cylindrical inclusion are studied optically. Test specimens with two different bond strengths (weak and strong) and three crack-inclusion eccentricities (e = 0, d/2 and 3d/4, d being inclusion diameter) are studied using reflection mode Coherent Gradient Sensing (CGS) and high-speed photography. These variants produce distinct dynamic crack trajectories and failure behaviors. A weaker inclusion-matrix interface attracts a propagating crack while a stronger one deflects the crack away. The former results in a propagating crack lodging (‘key-hole’) into the inclusion-matrix interface whereas in the latter the crack tends to circumvent the inclusion. When the inclusion is in the prospective crack path, the maximum attained crack speed is much higher in the weakly bonded inclusion cases relative to the strongly bonded counterparts. For a crack propagating towards a weakly bonded inclusion, the effective stress intensity factor (K _e) value remains constant for each inclusion eccentricity considered. But these constant K _e values increase with increasing eccentricity. A distinct drop in K _e occurs when the crack is near the inclusion. In strongly bonded inclusion cases, on the other hand, monotonically increasing K _e before the crack reaches the inclusion is observed. A drop in K _e is seen just before the crack reaches the inclusion. The mode-mixity estimates are of opposite signs for weakly and strongly bonded inclusions in case of the largest eccentricity studied, confirming the observed crack attraction and deflection mechanisms.

The Burgers vector and the flow of screw and edge dislocations in finite-deformation single-crystal plasticity

We consider finite plasticity based on the decomposition F=F^eF^p of the deformation gradient F into elastic and plastic distortions F^e and F^p. Within this framework the macroscopic Burgers vector may be characterized by the tensor field . We derive a natural convected rate for G associated with evolution of F^p and as our main result show that, for a single-crystal,

•

temporal changes in G—as characterized by its convected time derivative—may be decomposed into temporal changes in distributions of screw and edge dislocations on the individual slip systems.

We discuss defect energies dependent on the densities of these distributions and show that corresponding thermodynamic forces are macroscopic counterparts of classical Peach-Koehler forces.     

The hole-pressure problem

There is a need to unify present hypotheses of the nature and role of the hole-pressure,p _e , and thus provide consolidation on which to base future research and understanding. This paper is intended to meet this need. Attention is directed towards the calculation ofp _e from the velocity and stress fields for viscoelastic fluids flowingacross rectangular holes. The constitutive models used are the Newtonian, Second-order and Maxwell models, for values of Reynolds number up to 10 and Weissenberg number up to 0.1.The numerical complications involved are studied through an investigation of the constituent parts ofp _e . Verification of present theory is then sought, from which justification may be derived for the estimation of elasticity fromp _e measurements. Attention is directed towards the predictions of Higashitani and Pritchard and the extension to the Tanner and Pipkin theory for Second-order fluids. The effects of variation of geometric dimensions and flow type uponp _e are also discussed.     

A weakened form of Ilyushin's postulate and the structure of self-consistent Eulerian finite elastoplasticity

From a general standpoint in terms of internal variables, we formulate a general theory of self-consistent Eulerian finite elastoplasticity based on the additive decomposition of the Eulerian strain rate, i.e., D=D^e+D^p, as well as two consistency criteria. In this theory, the elastic behaviour is characterized by an exactly integrable elastic rate equation for D^e with a general form of complementary elastic potential. It is assumed that the yield function depends in a general manner on the Kirchhoff stress and the internal variables. Moreover, the plastic rate equation for D^p and the evolution equation for each internal variable are allowed to assume general forms relying on the just-mentioned variables and the stress rate. It is indicated that two consistency criteria, i.e., the self-consistency for the elastic rate equation and Prager's yielding stationarity, lead to the unique choice of objective rates, i.e., the logarithmic rate.The structure of the above theory is further studied and examined by virtue of a weakened form of Ilyushin's postulate. In a spinning frame defining the logarithmic rate, we introduce the notion of standard elastoplastic strain cycle, which starts at a point not on but inside a yield surface and incorporates only one infinitesimal plastic subpath. We show that this type of strain cycle is always possible. Then, by ruling out strain cycles starting at points on yield surfaces we propose a weakened form of Ilyushin's postulate, which says that the changing rate of the stress work done along every standard strain cycle should be non-negative, whenever the incorporated plastic subpath tends to vanish. By virtue of simple, rigorous procedures, we demonstrate that this weakened form of Ilyushin's postulate is adequate to ensure direct results concerning the normality rule and the convexity of the yield surface in the context of the foregoing Eulerian finite elastoplasticity theory. Specifically, with an exactly integrable elastic rate equation defining D^e, we prove that, in the space of the Kirchhoff stresses, the difference (D−D^e) is just the gradient of the yield function multiplied by a plastic multiplier, and thus bears the very kinematical and physical feature of plastic strain rate. Furthermore, we prove that, in the space of the Kirchhoff stresses, the elastic domain bounded by each yield surface should be convex. The main results are derived in a self-contained manner within the context of an Eulerian theory of finite elastoplasticity, without involving issues concerning how to define intermediate stress-free states and plastic strains, etc.     

Homogenization of an elasticity problem in domains with a net of slender bars near surface

We consider an elasticity problem in a domain Ω⁽⁾=ΩF⁽⁾, where Ω is an open bounded domain in R³, F⁽⁾ is a connected nonperiodic set in Ω like a net of slender bars, and is a parameter characterizing the microstructure of the domain. We consider the case of a surface distribution of the set F⁽⁾, i.e., for sufficiently small , the set F⁽⁾ is concentrated in arbitrary small neighbourhood of a surface Γ. Under a hypothesis on the asymptotic behaviour of the energy functional, we obtain the macroscopic (homogenized) model. To cite this article: M. Goncharenko, L. Pankratov, C. R. Mecanique 331 (2003).     

The decomposition F = FeFp,material symmetry,and plastic irrotationality for solids that are isotropic-viscoplastic or amorphous

This study develops a general framework for discussing both isotropic-viscoplastic materials and amorphous materials. The framework, which allows for large deformations, is based on the Kröner–Lee decomposition of the deformation gradient into elastic and inelastic parts, a system of microforces consistent with its own balance, and a mechanical version of the second law that includes, via the microforces, work performed during inelastic flow. The constitutive theory allows for dependences on the elastic and inelastic parts of the deformation gradient and on the inelastic stretch-rate, but dependences on the inelastic spin are not included. The constitutive equation for the microstress T^p conjugate to inelastic flow – suitably restricted by the second law – and the microforce balance are shown to be together equivalent to a flow rule that includes a back stress due to the variation in the free energy with inelastic deformation. The introduction of a concept of material microstability reduces this flow rule to one of classical Mises-type.In a theory based on the Kröner–Lee decomposition, there are two classes of symmetry transformations available: transformations of the reference configuration and transformations of the relaxed spaces. We discuss the notion of material symmetry for a general class of materials that includes, as special cases, isotropic-viscoplastic solids, and amorphous solids. Essential to this discussion of symmetry is a general constitutive relation for the microstress T^p.The symmetry-based framework allows us to show that for typical boundary-value problems involving isotropic, viscoplastic solids or amorphous solids, if a problem has a solution, then every time- and space-dependent rotation of the relaxed spaces also yields a solution, and it is possible to choose this rotation such that the transformed solution is inelastically spin-free: W^p ≡ 0. Thus, when discussing such materials, we may, without loss in generality, restrict attention to flow rules that are inelastically irrotational.     

Linear Stability Analysis of Resonant Periodic Motions in the Restricted Three-Body Problem

The equations of the restricted three-body problem describe the motion of a massless particle under the influence of two primaries of masses 1 −μ and μ, 0≤μ≤ 1/2, that circle each other with period equal to 2π. When μ=0, the problem admits orbits for the massless particle that are ellipses of eccentricity e with the primary of mass 1 located at one of the focii. If the period is a rational multiple of 2π, denoted 2π p/q, some of these orbits perturb to periodic motions for μ > 0. For typical values of e and p/q, two resonant periodic motions are obtained for μ > 0. We show that the characteristic multipliers of both these motions are given by expressions of the form in the limit μ→ 0. The coefficient C(e,p,q) is analytic in e at e=0 and C(e,p,q)=O(e^|p-q|). The coefficients in front of e^|p-q|, obtained when C(e,p,q) is expanded in powers of e for the two resonant periodic motions, sum to zero. Typically, if one of the two resonant periodic motions is of elliptic type the other is of hyperbolic type. We give similar results for retrograde periodic motions and discuss periodic motions that nearly collide with the primary of mass 1 −μ.     

Effect of thermal losses on the microscopic hyperbolic heat conduction model

The effects of radiative losses on the thermal behavior of thin metal films, as described by the microscopic two-step hyperbolic heat conduction model, are investigated. Different criteria, which determine the ranges within which thermal radiative losses are significant, are derived. It is found that radiative losses from the electron gas are significant in thin films having [(C_R e_e^4/₃ T_￥ ⁴ )/(k_e^1/₃ L^2/₃ G)] 3 4.6 ×10⁷{{C_R \epsilon _e^{{4 \over 3}} T_\infty ^4 } \over {k_e^{{1 \over 3}} L^{{2 \over 3}} G}}\geq 4.6 \times 10^7 for /_o > 4 and F_F < 1 and [(C_R e_e^3/₂ T_￥ ^9/₂)/(k_e^1/₂ L^1/₂ G)] 3 7.4 ×10¹⁰{{C_R \epsilon _e^{{3 \over 2}} T_\infty ^{{9 \over 2}}} \over {k_e^{{1 \over 2}} L^{{1 \over 2}} G}}\geq 7.4 \times 10^{10} for /_o < 4 and F_F > 1.     

An alternative derivation of the polar decomposition

This paper has showed that there is an alternative in deriving polar decomposition, which is quite different from Malvern's method[1]. Derivation procedure presented in this study has been started from the deformation gradient tensor F itself. Meanwhile, derivation procedure by Malvern [I] was started from the difference of the final squared length (ds)² of the element in the deformed configuration and the initial squared length (dS)² of the element in the undeformed configuration.     

Effects of Fuel Lewis Number on Localised Forced Ignition of Globally Stoichiometric Stratified Mixtures: a Numerical Investigation

The influences of fuel Lewis number Le_F on localised forced ignition of globally stoichiometric stratified mixtures have been analysed using three-dimensional compressible Direct Numerical Simulations (DNS) for cases with Le_F ranging from 0.8 to 1.2. The globally stoichiometric stratified mixtures with different values of root-mean-square (rms) equivalence ratio fluctuation (i.e. ?^′= 0.2, 0.4 and 0.6) and the Taylor micro-scale l_? of equivalence ratio ? variation (i.e. l_?/l_f= 2.1, 5.5 and 8.3 with l_f being the Zel’dovich flame thickness of the stoichiometric laminar premixed flame) have been considered for different initial rms values of turbulent velocity u^′. A pseudo-spectral method is used to initialise the equivalence ratio variation following a presumed bi-modal distribution for prescribed values of ?^′ and l_?/l_f for global mean equivalence ratio 〈?〉=1.0. The localised ignition is accounted for by a source term in the energy transport equation that deposits energy for a stipulated time interval. It has been observed that the maximum values of temperature and the fuel reaction rate magnitude increase with decreasing Le_F during the period of external energy deposition. The initial values of Le_F, u^′/S_b(?=1), ?^′ and l_?/l_f have been found to have significant effects on the extent of burning of the stratified mixtures following localised ignition. For a given value of u^′/S_b(?=1), the extent of burning decreases with increasing Le_F. An increase in u^′ leads to a monotonic reduction in the burned gas mass for all values of Le_F in all stratified mixture cases but an opposite trend is observed for the Le_F=0.8 homogeneous mixture. It has been found that an increase in ?^′ has adverse effects on the burned gas mass, whereas the effects of l_?/l_f on the extent of burning are non-monotonic and dependent on ?^′ and Le_F. Detailed physical explanations have been provided for the observed Le_F, u^′/S_b(?=1), ?^′ and l_?/l_f dependences.     

On the orders of the best approximations of integrals of functions by integrals of rank σ

We study the values e _σ(f) of the best approximation of integrals of functions from the spaces L _p (A, dμ) by integrals of rank σ. We determine the orders of the least upper bounds of these values as σ → ∞ in the case where the function ƒ is the product of two nonnegative functions one of which is fixed and the other varies on the unit ball U _p (A) of the space L _p (A, dμ). We consider applications of the obtained results to approximation problems in the spaces S _p ^ϕ. __________ Translated from Neliniini Kolyvannya, Vol. 10, No. 4, pp. 528–559, October–December, 2007.     

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.     

On the Euler Equations in the Critical Triebel-Lizorkin Spaces

In this paper we study the initial value problem of the incompressible Euler equations in ⁿ for initial data belonging to the critical Triebel-Lizorkin spaces, i.e., v ₀ F ⁿ⁺¹ _1,q , q[1, ]. We prove the blow-up criterion of solutions in F ⁿ⁺¹ _1,q for n=2,3. For n=2, in particular, we prove global well-posedness of the Euler equations in F ³ _1,q , q[1, ]. For the proof of these results we establish a sharp Moser-type inequality as well as a commutator-type estimate in these spaces. The key methods are the Littlewood-Paley decomposition and the paradifferential calculus by J. M. Bony.