On the plastic wave speeds in rate-independent elastic–plastic materials with anisotropic elasticity

This paper presents a theoretical study of the speeds of plastic waves in rate-independent elastic–plastic materials with anisotropic elasticity. It is shown that for a given propagation direction the plastic wave speeds are equal to or lower than the corresponding elastic speeds, and a simple expression is provided for the bound on the difference between the elastic and the plastic wave speeds. The bound is given as a function of the plastic modulus and the magnitude of a vector defined by the current stress state and the propagation direction. For elastic–plastic materials with cubic symmetry and with tetragonal symmetry, the upper and lower bounds on the plastic wave speeds are obtained without numerically solving an eigenvalue problem. Numerical examples of materials with cubic symmetry (copper) and with tetragonal symmetry (tin) are presented as a validation of the proposed bounds. The lower bound proposed here on the minimum plastic wave speed may also be used as an efficient alternative to the bifurcation analysis at early stages of plastic deformation for the determination of the loss of ellipticity.     

An anisotropic elastic–viscoplastic model for soft clays

Experimental evidences have shown deficiencies of the existing overstress and creep models for viscous behaviour of natural soft clay. The purpose of this paper is to develop a modelling method for viscous behaviour of soft clays without these deficiencies. A new anisotropic elastic–viscoplastic model is extended from overstress theory of Perzyna. A scaling function based on the experimental results of constant strain-rate oedometer tests is adopted, which allows viscoplastic strain-rate occurring whether the stress state is inside or outside of the yielding surface. The inherent and induced anisotropy is modelled using the formulations of yield surface with kinematic hardening and rotation (S-CLAY1). The parameter determination is straightforward and no additional experimental test is needed, compared to the Modified Cam Clay model. Parameters determined from two types of tests (i.e., the constant strain-rate oedometer test and the 24 h standard oedometer test) are examined. Experimental verifications are carried out using the constant strain-rate and creep tests on St. Herblain clay. All comparisons between predicted and measured results demonstrate that the proposed model can successfully reproduce the anisotropic and viscous behaviours of natural soft clays under different loading conditions.     

Orthotropic elastic–plastic material model for paper materials

Paper and paperboard generally exhibit anisotropic and non-linear mechanical material behaviour. In this work, the development of an orthotropic elastic–plastic constitutive model, suitable for modelling of the material behaviour of paper is presented. The anisotropic material behaviour is introduced into the model by orthotropic elasticity and an isotropic plasticity equivalent transformation tensor. A parabolic stress–strain relation is adopted to describe the hardening of the material. The experimental and numerical procedures for evaluation of the required material parameters for the model are described. Uniaxial tensile testing in three different inplane material directions provides the calibration of the material parameters under plane stress conditions. The numerical implementation of the material model is presented and the model is shown to perform well in agreement with experimentally observed mechanical behaviour of paper.     

Derivatives of the three-dimensional Green’s functions for anisotropic materials

A concise formulation is presented for the derivatives of Green’s functions of three-dimensional generally anisotropic elastic materials. Direct calculation for derivatives of the Green’s function on the Cartesian coordinate system is a common practice, which, however, usually leads to a complicated course. In this paper the Green’s function derived by Ting and Lee [Ting, T.C.T., Lee, V.G., 1997. The three-dimensional elastostatic Green’s function for general anisotropic linear elastic solids. The Quarterly Journal of Mechanics and Applied Mathematics 50 (3) 407–426] is extended to obtain the derivatives. Using a spherical coordinate system, the Green’s function can be shown as the composition of two independent functions, one depends only on the radial distance of the field point to the origin and the other is in spherical angles. The method of derivation is based on the total differential scheme and then takes its partial differentiation accordingly. With the application of Cauchy residue theorem, the contour integral can be evaluated in terms of the Stroh eigenvalues of a sextic equation. For the degenerate case, evaluation of residues at multiple poles is also given. Applications of the present result are made to examine the Green’s functions and stress components for isotropic and transversely isotropic materials. The results are in exact agreement with existing solutions.     

Structural shakedown for elastic–plastic materials with hardening saturation surface

An elastic–plastic material model with internal variables and thermodynamic potential, not admitting hardening states out of a saturation surface, is assumed as a basis to formulate a statical Melan-type shakedown theorem. Grounding on the optimality conditions relative to the shakedown load multiplier problem for a structure subjected to cyclic loads, the impending inadaptation collapse mechanism at the shakedown limit state is analyzed and discussed. It is shown that the adopted model is able to catch ratchetting collapse mode at a structural level. Numerical results for a simple structure are finally reported.     

A new class of micro–macro models for elastic–viscoplastic heterogeneous materials

The determination of the effective behavior of heterogeneous materials from the properties of the components and the microstructure constitutes a major task in the design of new materials and the modeling of their mechanical behavior. In real heterogeneous materials, the simultaneous presence of instantaneous mechanisms (elasticity) and time dependent ones (non-linear viscoplasticity) leads to a complex space–time coupling between the mechanical fields, difficult to represent in a simple and efficient way. In this work, a new self-consistent model is proposed, starting from the integral equation for a translated strain rate field. The chosen translated field is the (compatible) viscoplastic strain rate of the (fictitious) viscoplastic heterogeneous medium submitted to a uniform (unknown) boundary condition. The self-consistency condition allows to define these boundary conditions so that a relative simple and compact strain rate concentration equation is obtained. This equation is explained in terms of interactions between an inclusion and a matrix, which lead to interesting conclusions. The model is first applied to the case of two-phase composites with isotropic, linear and incompressible viscoelastic properties. In that case, an exact self-consistent solution using the Laplace–Carson transform is available. The agreement between both approaches appears quite good. Results for elastic–viscoplastic BCC polycrystals are also presented and compared with results obtained from Kröner–Weng's and Paquin et al. (Arch. Appl. Mech. 69 (1999) 14)'s model.     

Homogenized elastic–viscoplastic behavior of anisotropic open-porous bodies with pore pressure

Constitutive modeling is studied for the homogenized elastic–viscoplastic behavior of pore-pressurized anisotropic open-porous bodies made of metallic base solids at small strains and rotations. For this purpose, by describing micro–macro relations relevant to periodic unit cells of anisotropic open-porous bodies subjected to pore pressure, constitutive features are discussed for the viscoplastic macrostrain rate in steady states. On the basis of the constitutive features found, the viscoplastic macrostrain rate is represented as an anisotropic function of Terzaghi’s effective stress, which is shown using Hill’s macrohomogeneity condition. The resulting viscoplastic equation is used to simulate the homogenized elastic–viscoplastic behavior of an ultrafine plate-fin structure subjected to uniaxial/biaxial loading in addition to pore pressure. The corresponding finite element homogenization analysis is also performed for comparison. It is demonstrated that the developed viscoplastic equation simulates well the anisotropic effect of pore pressure in the viscoplastic range in spite of there being no anisotropic factor and no fitting parameter in Terzaghi’s effective stress itself.     

A new formulation of the effective elastic–plastic response of two-phase particulate composite materials

In this paper the double-inclusion model, originally developed to determine effective linear elastic properties of composite materials, is reformulated and extended to predict the effective nonlinear elastic–plastic response of two-phase particulate composites reinforced with spherical particles. The resulting problem of elastic–plastic deformation of a double-inclusion embedded in an infinite reference medium subjected to an incrementally applied far-field strain is solved by the finite element method. The proposed double-inclusion model is evaluated by comparison of the model predictions to the available exact results obtained by the direct approach using representative volume elements containing many particles. It is found that the double-inclusion formulation is capable of providing accurate prediction of the effective elastic–plastic response of two-phase particulate composites at moderate particle volume fractions.     

On the forward and backward in time problems in the Kelvin–Voigt thermoviscoelastic materials

This paper studies the uniqueness of solutions to the forward and backward in time boundary value problems associated with the Kelvin–Voigt viscoelastic model of the thermoelastic materials. For thermoviscoelastic materials with a center of symmetry, it is shown the uniqueness of solutions to the forward in time boundary value problems without any assumptions upon the thermoviscoelastic constitutive coefficients other than the symmetry properties and those induced by the dissipation inequality. While for the final boundary value problems two uniqueness theorems are presented: the first one is essentially based on the assumption that the specific heat is of negative definite sign, while the second is established in the class of displacement–temperature variation fields whose dissipation energy has a temporal behavior lower than an appropriate growing exponential.     

On the impact of a rigid–plastic missile into rigid or elastic target

Here we carry out a systematic parametric study of a uniform cylindrical missile impacting rigid or elastic structures. We give an analytical result for the impact force in case of rigid target. A new parameter, the damage potential is introduced and it is shown that this single dimensionless combination of the parameters describes the course of the impact in this simplest case. For elastic target structures, we also show numerically that the course of the reaction force, the maximum target displacement and the duration of the impact depend primarily on the same dimensionless parameter with a secondary effect of the missile to target mass ratio and the relative stiffness of the target. The rigid target assumption is not always conservative with regard to the reaction force due to target vibration. We find a resonant effect in the maximum target displacement as the function of the missile to target mass ratio. The motivation of our work is rooted in the investigation of aircraft fuselage impact into robust structures like the containment of a nuclear power plant.     

On the stress–strain states of cellular materials under high loading rates

A virtual Taylor impact of cellular materials is analyzed with a wave propagation technique, i.e. the Lagrangian analysis method, of which the main advantage is that no pre-assumed constitutive relationship is required. Time histories of particle velocity, local strain, and stress profiles are calculated to present the local stress–strain history curves, from which the dynamic stress–strain states are obtained.The present results reveal that the dynamic-rigid-plastic hardening(D-R-PH) material model introduced in a previous study of our group is in good agreement with the dynamic stress–strain states under high loading rates obtained by the Lagrangian analysis method. It directly reflects the effectiveness and feasibility of the D-R-PH material model for the cellular materials under high loading rates.     

Fluid–structure interaction in water-filled thin pipes of anisotropic composite materials

The effects of elastic anisotropy in piping materials on fluid–structure interaction are studied for water-filled carbon-fiber reinforced thin plastic pipes. When an impact is introduced to water in a pipe, there are two waves traveling at different speeds. A primary wave corresponding to a breathing mode of pipe travels slowly and a precursor wave corresponding to a longitudinal mode of pipe travels fast. An anisotropic stress–strain relationship of piping materials has been taken into account to describe the propagation of primary and precursor waves in the carbon-fiber reinforced thin plastic pipes. The wave speeds and strains in the axial and hoop directions are calculated as a function of carbon-fiber winding angles and compared with the experimental data. As the winding angle increases, the primary wave speed increases due to the increased stiffness in the hoop direction, while the precursor wave speed decreases. The magnitudes of precursor waves are much smaller than those of primary waves so that the effect of precursor waves on the deformation of pipe is not significant. The primary wave generates the hoop strain accompanying the opposite-signed axial strain through the coupling compliance of pipe. The magnitude of hoop strain induced by the primary waves decreases with increasing the winding angle due to the increased hoop stiffness of pipe. The magnitude of axial strain is small at low and high winding angles where the coupling compliance is small.     

On Green’s function for a bimaterial elastic half-plane

The problem of a point force acting in a composite, two-dimensional, isotropic elastic half-plane is considered. An exact solution is obtained, using Mellin transforms and the Melan solution for a point force in a homogeneous half-plane.     

New integral formulation and self-consistent modeling of elastic–viscoplastic heterogeneous materials

Predicting the overall behavior of heterogeneous materials, from their local properties at the scale of heterogeneities, represents a critical step in the design and modeling of new materials. Within this framework, an internal variables approach for scale transition problem in elastic–viscoplastic case is introduced. The proposed micromechanical model is based on establishing a new system of field equations from which two Navier’s equations are obtained. Combining these equations leads to a single integral equation which contains, on the one hand, modified Green operators associated with elastic and viscoplastic reference homogeneous media, and secondly, elastic and viscoplastic fluctuations. This new integral equation is thus adapted to self-consistent scale transition methods. By using the self-consistent approximation we obtain the concentration law and the overall elastic–viscoplastic behavior of the material. The model is first applied to the case of two-phase materials with isotropic, linear and compressible viscoelastic properties. Results for elastic–viscoplastic two-phase materials are also presented and compared with exact results and variational methods.     

Homogenization of elastic–viscoplastic heterogeneous materials: Self-consistent and Mori-Tanaka schemes

This paper deals with the prediction of the macroscopic behavior of a multiphase elastic–viscoplastic material. The proposed homogenization schemes are based on an interaction law postulated by Molinari et al. [Molinari, A., Ahzi, S., Kouddane, R. 1997. On the self-consistent modelling of elastic–plastic behavior of polycrystals. Mech. Mater., 26, 43–62]. Self-consistent schemes are developed to describe the behavior of disordered aggregates. The Mori-Tanaka approach is used to capture the behavior of composite materials, where one phase can be clearly identified as the matrix. The proposed schemes are developed within a general framework where compressible elasticity and anisotropy of the materials are taken into account. Inclusions can have various shapes and orientations. Illustrations of the homogenization procedure are given for a two-phase composite materials. Comparisons between results of the literature and predictions based on the interaction law are performed and have demonstrated the efficiency of the proposed homogenization schemes.     

A note on divergence and flutter instabilities in elastic–plastic materials

Dynamic stability of uniform straining of a nonlinear elastic solid is known to require that all eigenvalues of the acoustic tensor associated with the tangent elastic moduli be real and nonnegative. The focus of this note is to what extent this conclusion applies to time-independent, elastoplastic materials. Nonlinearity of the elastic–plastic constitutive law imposes limits on validity of a solution to the linear problem for which the acoustic tensor is determined. The effect of those limits on the conclusions about instability is examined.     

Lamé problem for a weakly anisotropic body with cubic symmetry of elastic properties

The Lamé problem is solved for a body with cubic symmetry of elastic properties and the elastic anisotropy parameter is determined. In the case of plane deformation, the stresses in a ring are found to terms of the first order in the small anisotropic parameter. Stresses in a ring of KCl under internal pressure are calculated.     

Shakedown theorems for elastic–plastic solids in the framework of gradient plasticity

Static and kinematic shakedown theorems are given for a class of generalized standard materials endowed with a hardening saturation surface in the framework of strain gradient plasticity. The so-called residual-based gradient plasticity theory is employed. The hardening law admits a hardening potential, which is a C¹-continuous function of a set of kinematic internal variables and of their spatial gradients, and is required to satisfy a global sign restriction (but not to be necessarily convex). The totally produced, the accumulated and the freely moving dislocations per unit volume, distinguished as statistically stored and geometrically necessary ones, are in this way accounted for. Like for a generalized standard material, the shakedown safety factor is found to depend on the (generally size dependent) yield and saturation limits, but not on the particular differential-type hardening law of the material.     

Overall longitudinal shear elastic modulus of a 1–3 composite with anisotropic constituents

The overall properties of a binary elastic periodic fiber-reinforced composite, with transversely isotropic constituents in an anti-plane-strain deformation state, are studied here for a cell periodicity of square type. This analysis considers four different orientations of the axis of transverse isotropy of constituents with respect to the direction of fibers. Each case is characterized by very simple closed-form expressions for the effective coefficients, which were obtained using the asymptotic homogenization method. Local problems defined on a periodic square unit cell are solved using Weierstrassian and Natanzon’s functions and perturbation theory relative to small anisotropy. In the isotropic limit, comparison with rigorous bounds and some well-known mixing rules are made. Also, comparisons with finite element calculations show that the derived closed-form formulae provide excellent results even for large anisotropy.     

On the Regularity Set and Angular Integrability for the Navier–Stokes Equation

We investigate the size of the regular set for suitable weak solutions of the Navier–Stokes equation, in the sense of Caffarelli–Kohn–Nirenberg (Commun Pure Appl Math 35:771–831, 1982). We consider initial data in weighted Lebesgue spaces with mixed radial-angular integrability, and we prove that the regular set increases if the data have higher angular integrability, invading the whole half space \({\{t > 0\}}\) in an appropriate limit. In particular, we obtain that if the \({L^{2}}\) norm with weight \({|x|^{-\frac12}}\) of the data tends to 0, the regular set invades \({\{t > 0\}}\); this result improves Theorem D of Caffarelli et al. (Commun Pure Appl Math 35:771–831, 1982).