Limit state of structural elements during inelastic deformation

It is shown that, in the case of an axisymmetric stress state, the solution of the statically definable boundary-value problem for an ideal rigid-plastic body using the Mises-Schleicher strength criterion is extended to the rigid-creep model with any specified creep-rupture strength and corresponds to the limit state of a real creeping body.     

Further criteria for elastic waves in anisotropic media

Two additional criteria for the existence of cusp points on elastic wave surfaces are developed.A previously published method [1] is extended to give a simple necessary and sufficient condition for cusps about (1, 1, 0) axes in cubic and tetragonal media. This criterion is plausibly adapted to provide a simple inequality applicable to any section of slowness surface represented by separable quadratic and quartic equations.Two tables of numerical examples are presented.     

On the modelling of anisotropic elastic and inelastic material behaviour at large deformation

The purpose of this work is the formulation and discussion of an approach to the modelling of anisotropic elastic and inelastic material behaviour at large deformation. This is done in the framework of a thermodynamic, internal-variable-based formulation for such a behaviour. In particular, the formulation pursued here is based on a model for plastic or inelastic deformation as a transformation of local reference configuration for each material element. This represents a slight generalization of its modelling as an elastic material isomorphism pursued in earlier work, allowing one in particular to incorporate the effects of isotropic continuum damage directly into the formulation. As for the remaining deformation- and stress-like internal variables of the formulation, these are modelled in a fashion formally analogous to so-called structure tensors. On this basis, it is shown in particular that, while neither the Mandel nor back stress is generally so, the stress measure thermodynamically conjugate to the plastic “velocity gradient”, containing the difference of these two stress measures, is always symmetric with respect to the Euclidean metric, i.e., even in the case of classical or induced anisotropic elastic or inelastic material behaviour. Further, in the context of the assumption that the intermediate configuration is materially uniform, it is shown that the stress measure thermodynamically conjugate to the plastic velocity gradient is directly related to the Eshelby stress. Finally, the approach is applied to the formulation of metal plasticity with isotropic kinematic hardening.     

Rapoport-leas two-phase flow model for anisotropic porous media

The Rapoport-Leas mathematical model of two-phase flow is generalized to include the case of anisotropic porous media. The formula for the capillary pressure, which specifies the relationship between the phase pressures, contains a scalar function of a vector argument. In order to determine the scalar function, the capillary pressure tensor and the tensor inverse to the tensor of characteristic linear dimensions are introduced. The capillary pressure is determined by the contraction of the second-rank tensors with a unit vector collinear to the phase pressure gradients, also assumed to be collinear. It is shown that the saturation function introduced for isotropic porous media (Leverett function) can be generalized to include anisotropic media and is now determined by a fourth-rank tensor. Generalized expressions for the Leverett and relative phase permeability functions are given for orthotropic and transversely isotropic media with account for the hysteresis of the phase permeabilities and capillary pressure.     

A general interface model for a three-dimensional curved thin anisotropic interphase between two anisotropic media

An arbitrarily curved three-dimensional anisotropic thin interphase between two anisotropic solids is considered. The purpose of this study is to model this interphase as a surface between its two neighbouring media by means of appropriately devised interface conditions on it. The analysis is carried out in the setting of unsteady heat conduction and dynamic elasticity, and makes use of the simple idea of a Taylor expansion of the relevant fields in thin regions. It consists of a generalization of a previous study by Bövik [1994. On the modelling of thin interface layers in elastic and acoustic scattering problems. Q. J. Mech. Appl. Math. 47, 17-42] which was confined to the isotropic setting. The remarkable feature of the presently derived anisotropic interface model is that formally it has a more compact form than that of Bövik's isotropic version. This is achieved by a judicious choice of surface differential operators which have been used in the derivation, and makes possible to show that several previously known classical interface models are recovered as special cases of the one obtained in this study, once suitable assumptions are made on the magnitude of the conductivity and elasticity tensors of the interphase.     

Random-field model for the elasticity tensor of anisotropic random media

This Note deals with the construction of a non-Gaussian positive definite matrix-valued random field whose mathematical properties allow the fourth-order elasticity tensor of random non homogeneous anisotropic three dimensional elastic media to be modelled. If the usual parametric probabilistic approach was used, then 21 mutually dependent random fields should be modelled and identified by using experimental data. Such an approach would be very difficult because the systems of the marginal probability distributions of these random fields have to be identified due to the fact that, for a boundary value problem, the displacement field of the random medium is a non-linear mapping of the random elasticity tensor. The theory presented in this paper allows such a probabilistic model of the fourth-order elasticity tensor field to be constructed and depends only of four scalar parameters: three spatial correlation lengths and one parameter allowing the level of the random fluctuations to be controlled. To cite this article: C. Soize, C. R. Mecanique 332 (2004).

Résumé

On présente la construction d'un champ aléatoire à valeurs dans les matrices définies positives dont les propriétés mathématiques permettent de modéliser le tenseur d'élasticité du quatrième ordre des mileux élastiques anisotropes tridimensionnels aléatoires. Si l'approche probabiliste paramétrique usuelle était utilisée, alors il serait nécessaire de modéliser et d'identifier à l'aide de données expérimentales 21 champs aléatoires mutuellement dépendants. Une telle approche serait très difficile de part le fait que le système de lois marginales de ces champs aléatoires doit être identifié parce que, pour un problème aux limites, le champ de déplacement est une transformation non linéaire du tenseur d'élasticité. La théorie présentée dans ce papier permet de construire une modélisation probabiliste du champ de tenseur d'élasticité qui ne dépend que de quatre paramètres scalaires : trois échelles de corrélation spatiale et un paramètre permettant de contrôler le niveau des fluctuations aléatoires. Pour citer cet article : C. Soize, C. R. Mecanique 332 (2004).     

Formulation of a finite deformation model for the dynamic response of open cell biphasic media

The paper illustrates a biphasic formulation which addresses the dynamic response of fluid saturated porous biphasic media at finite deformations with no restriction on the compressibility of the fluid and of the solid skeleton. The proposed model exploits four state fields of purely kinematic nature: the displacements of the solid phase, the velocity of the fluid, the density of the fluid and an additional macroscopic scalar field, termed effective Jacobian, associated with the effective volumetric deformation of the solid phase.The governing equations are characterized by the property of being all expressed in the reference configuration of the solid phase and by the property of employing only work-conjugate variables, thus avoiding the use of a total Cauchy stress tensor.In particular, the set of governing equations includes a momentum balance equation associated with the effective Jacobian field. This equation, differently from the closure-equations proposed by other authors which express a saturation constraint or a porosity balance, is derived as a stationarity condition on account of a least-action variational principle.     

Criteria for elastic waves in anisotropic media — a consolidation

The geometry relating to the tangent plane at a stationary point on a surface has been used to re-examine various criteria for the existence of parabolic points (inflexions in 2D-sections) on the outermost sheet of the slowness surface for elastic waves in anisotropic media.Previous results obtained by the authors, exact [4,5], and approximate [6], are related in detail. Further approximations, based on geometrical properties, are derived; one of these proves equivalent to a sufficient condition first applied by McCurdy [8].Numerical investigation shows that, over a wide range of anisotropy, the simply applied approximate criterion of [6] is sufficient and within the accuracy of observation of the elastic stiffnesses.     

A physically-based constitutive model for anisotropic damage in rubber-toughened glassy polymers during finite deformation

The present work focuses on the development of a physically-based model for large deformation stress-strain response and anisotropic damage in rubber-toughened glassy polymers. The main features leading to a microstructural evolution (regarding cavitation, void aspect ratio, matrix plastic anisotropy and rubbery phase deformation) in rubber-toughened glassy polymers are introduced in the proposed constitutive model. The constitutive response of the glassy polymer matrix is modelled using the hyperelastic-viscoplastic model of [Boyce et al., 1988] and [Boyce et al., 2000]. The deformation mechanisms of the matrix material are accounted for by two resistances: an elastic-viscoplastic isotropic intermolecular resistance acting in parallel with a visco-hyperelastic anisotropic network resistance, each resistance being modified to account for damage effects by void growth with a variation of the void aspect ratio. The effective contribution of the hyperelastic particles to the overall composite behaviour is taken into account by treating the overall system in a composite scheme framework. The capabilities of the proposed constitutive model are checked by comparing experimental data with numerical simulations. The deformation behaviour of rubber-toughened poly(methyl methacrylate) was investigated experimentally in tension at a temperature of 80 °C and for different constant true strain rates monitored by a video-controlled technique. The reinforcing phase is of the soft core-hard shell type and its diameter is of the order of one hundred nanometers. The particle volume fraction was adjusted from 15% to 45% by increments of 5%. The stress-strain response and the inelastic volumetric strain are found to depend markedly on particle volume fraction. For a wide range of rubber volume fractions, the model simulations are in good agreement with the experimental results. Finally, a parametric analysis demonstrates the importance of accounting for void shape, matrix plastic anisotropy and rubber content.     

Generalized rheological model of shear deformation of structured condensed media

A relaxation-type macrorheological equation is constructed, which describes shear deformation of elastoplastoviscous media containing microvoids. The equation corresponds to regimes of steady-state creep and elastic and plastic strains, which pass to the regime of viscous flow with limited growth of voids at high shear strain rates. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 5, pp. 87–95, September–October, 2007.     

Incremental deformation model for a rod

A nonlinear deformation model for a rod with rigid cross sections is proposed. A complete system of local incremental equations, a variational equation equivalent to this system, and an equation of virtual work are formulated. Numerical analysis of the deformation of a ring transmission is performed. Institute of Computer Modeling, Siberian Division, Russian Academy of Sciences, Krasnoyarsk 660036. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 4, pp. 229–235, July–August, 1999.     

Incremental deformation model for a shell

A nonlinear deformation model for a shell with rigid transverse fibers is proposed. A complete system of incremental equations, a variational equation equivalent to this system, and a particular equation of virtual work are formulated. Numerical analysis of the nonlinear deformation of a spherical dome is performed using the complete equation. Institute of Computer Modeling, Siberian Division, Russian Academy of Sciences, Krasnoyarsk 660036. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 5, pp. 202–207, September–October, 1999.     

Experimental stress analysis for process of nonsteady inelastic deformation by photoviscoplasticity

The proposed photoviscoplasticity is a method for determining the stress distribution in nonsteady inelastic deformation during creep. The governing fundamental relations are derived by considering the effect of time as involved in the viscoplastic strain rate and the stress rate, and by considering the distinctive deformation properties of celluloid as model material. The validity of the fundamental relations are evaluated according to calibration test by using celluloid. The proposed method is applied for two practical applications of the compression of blocks by elastic punches and the compression of strip having a circular hole or semicircular notches. The time-dependent variation of stress states can be analyzed with the same accuracy as in the photoelasticity.     

Non-quadratic Kelvin modes based plasticity criteria for anisotropic materials

Novel (non-quadratic) plasticity criteria based on Kelvin modes are formulated here for anisotropic materials. As an example, such a macroscopic criterion is applied with success to the case of FCC nickel-base single crystals. Indeed, relying on the cubic symmetry of the material, the Kelvin decomposition of elasticity tensor easily allows for the definition of an objective and loading independent criterion. The criterion identification is performed from different loading cases for CMSX2 single crystal superalloy. Tension-torsion yield surfaces at room temperature and yield stress dependence on crystal orientation are modeled. The Kelvin modes based criterion is compared to experimental data, to Hill and Barlat and coworkers macroscopic criteria and to Schmid law predictions. The results show that a simple three-parameter yield function built thanks to von Mises equivalent Kelvin stresses accounts for a satisfying plasticity criterion for such alloys.Non-quadratic norm ∥·∥_a plasticity framework is addressed. Intrinsic generalizations of Hershey-Hosford criterion are proposed for cubic material symmetry.     

Strain localization analysis using a large deformation anisotropic elastic–plastic model coupled with damage

Sheet metal forming processes generally involve large deformations together with complex loading sequences. In order to improve numerical simulation predictions of sheet part forming, physically-based constitutive models are often required. The main objective of this paper is to analyze the strain localization phenomenon during the plastic deformation of sheet metals in the context of such advanced constitutive models. Most often, an accurate prediction of localization requires damage to be considered in the finite element simulation. For this purpose, an advanced, anisotropic elastic–plastic model, formulated within the large strain framework and taking strain-path changes into account, has been coupled with an isotropic damage model. This coupling is carried out within the framework of continuum damage mechanics. In order to detect the strain localization during sheet metal forming, Rice’s localization criterion has been considered, thus predicting the limit strains at the occurrence of shear bands as well as their orientation. The coupled elastic–plastic-damage model has been implemented in Abaqus/implicit. The application of the model to the prediction of Forming Limit Diagrams (FLDs) provided results that are consistent with the literature and emphasized the impact of the hardening model on the strain-path dependency of the FLD. The fully three-dimensional formulation adopted in the numerical development allowed for some new results – e.g. the out-of-plane orientation of the normal to the localization band, as well as more realistic values for its in-plane orientation.     

Distribution of random slip-line directions under inelastic deformation

An expression for the probability density of random slip-line directions on longitudinal thin sections of inelastic specimens under tension is obtained. Various hypotheses on random angles of inclination of slip planes are used. The numerical results are compared with known experimental data.