Discrete dislocation plasticity analysis of the wedge indentation of films

The plane strain indentation of single crystal films on a rigid substrate by a rigid wedge indenter is analyzed using discrete dislocation plasticity. The crystals have three slip systems at ±35.3^° and 90^° with respect to the indentation direction. The analyses are carried out for three values of the film thickness, 2, 10 and , and with the dislocations all of edge character modeled as line singularities in a linear elastic material. The lattice resistance to dislocation motion, dislocation nucleation, dislocation interaction with obstacles and dislocation annihilation are incorporated through a set of constitutive rules. Over the range of indentation depths considered, the indentation pressure for the 10 and thick films decreases with increasing contact size and attains a contact size-independent value for contact lengths . On the other hand, for the films, the indentation pressure first decreases with increasing contact size and subsequently increases as the plastic zone reaches the rigid substrate. For the 10 and thick films sink-in occurs around the indenter, while pile-up occurs in the film when the plastic zone reaches the substrate. Comparisons are made with predictions obtained from other formulations: (i) the contact size-independent indentation pressure is compared with that given by continuum crystal plasticity; (ii) the scaling of the indentation pressure with indentation depth is compared with the relation proposed by Nix and Gao [1998. Indentation size effects in crystalline materials: a law for strain gradient plasticity. J. Mech. Phys. Solids 43, 411-423]; and (iii) the computed contact area is compared with that obtained from the estimation procedure of Oliver and Pharr [1992. An improved technique for determining hardness and elastic-modulus using load and displacement sensing indentation experiments, J. Mater. Res. 7, 1564-1583].     

Limit analysis of multi-layered plates. Part I: The homogenized Love-Kirchhoff model

The purpose of this paper is to determine , the overall homogenized Love-Kirchhoff strength domain of a rigid perfectly plastic multi-layered plate, and to study the relationship between the 3D and the homogenized Love-Kirchhoff plate limit analysis problems. In the Love-Kirchhoff model, the generalized stresses are the in-plane (membrane) and the out-of-plane (flexural) stress field resultants. The homogenization method proposed by Bourgeois [1997. Modélisation numérique des panneaux structuraux légers. Ph.D. Thesis, University Aix-Marseille] and Sab [2003. Yield design of thin periodic plates by a homogenization technique and an application to masonry wall. C. R. Méc. 331, 641-646] for in-plane periodic rigid perfectly plastic plates is justified using the asymptotic expansion method. For laminated plates, an explicit parametric representation of the yield surface is given thanks to the π-function (the plastic dissipation power density function) that describes the local strength domain at each point of the plate. This representation also provides a localization method for the determination of the 3D stress components corresponding to every generalized stress belonging to . For a laminated plate described with a yield function of the form , where σ^u is a positive even function of the out-of-plane coordinate x₃ and is a convex function of the local stress σ, two effective constants and a normalization procedure are introduced. A symmetric sandwich plate consisting of two Von-Mises materials ( in the skins and in the core) is studied. It is found that, for small enough contrast ratios (), the normalized strength domain is close to the one corresponding to a homogeneous Von-Mises plate [Ilyushin, A.-A., 1956. Plasticité. Eyrolles, Paris].     

Nanoindentation of polymeric thin films with an interfacial force microscope

The mechanical properties of interphase regions at bi-material interfaces can be quite different from the surrounding bulk materials. For composite materials, this interphase region is usually thin but plays an important role in their overall mechanical properties. Nanoindentation has become a commonly used experimental technique for measuring the mechanical properties of materials, especially when one of the dimensions is small. However, the extraction of reduced elastic modulus from the nanoindentation of thin films on substrates can pose challenges due to the influence of the substrate. In this study, the nanoindentation of thin films on substrates has been examined with a view to extracting the reduced modulus of thin polymer films.Thin films of (3-aminopropyl)triethoxysilane (C₉H₂₃NO₃Si, γ-APS) were deposited on silicon. An interfacial force microscope (IFM) was used to indent the γ-APS films. The effect of the substrate was studied by considering two very different thicknesses ( and ). The nanoindentation data were analyzed via contact mechanics theories and a finite element analysis that incorporated surface interactions. The analyses showed that nanoindentation experiments can provide reliable values of film modulus when the film is very different from the substrate. It was found that the commonly used rule of thumb that the indentation depth should be less than 10% of the thickness did not eliminate substrate effects for a wide range of material combinations. Instead, it is proposed that the contact radius should be less than 10% of the thickness so that contact mechanics theories for monolithic materials can be used without considering the presence of the substrate. The modulus of γ-APS polymer films and the surface energy between the tungsten tip of the IFM and γ-APS films were extracted and were related to their cure. A completely cured thick γ-APS film had a reduced modulus of . This value falls in the usual range for polymers due to the amorphous nature of the γ-APS films.     

Nanoindentation and the indentation size effect: Kinetics of deformation and strain gradient plasticity

A study of the indentation size effect (ISE) in aluminum and alpha brass is presented. The study employs rate effects to examine the fundamental mechanisms responsible for the ISE. These rate effects are characterized in terms of the rate sensitivity of the hardness, , where H is the hardness and is an effective strain rate in the plastic volume beneath the indenter. can be measured using indentation creep, load relaxation, or rate change experiments. The activation volume V∗, calculated based on which can traditionally be used to compare rate sensitivity data from a hardness test to conventional uniaxial testing, is calculated. Using materials with different stacking fault energy and specimens with different levels of work hardening, we demonstrate how increasing the dislocation density affects V∗; these effects may be taken as a kinetic signature of dislocation strengthening mechanisms. We noticed both H and exhibit an ISE. The course of V∗ vs. H as a result of the ISE is consistent with the course of testing specimens with different level of work hardening. This result was observed in both materials. This suggests that a dislocation mechanism is responsible for the ISE. When the results are fitted to a strain gradient plasticity model, the data at deep indents (microhardness and large nanoindentation) exhibit a straight-line behavior closely identical to literature data. However, for shallow indents (nanoindentation data), the slope of the line severely changes, decreasing by a factor of 10, resulting in a “bilinear behavior”.     

Some new classes of inverse coefficient problems in non-linear mechanics and computational material science

Three classes of inverse coefficient problems arising in engineering mechanics and computational material science are considered. Mathematical models of all considered problems are proposed within the J₂-deformation theory of plasticity. The first class is related to the determination of unknown elastoplastic properties of a beam from a limited number of torsional experiments. The inverse problem here consists of identifying the unknown coefficient g(ξ²) (plasticity function) in the non-linear differential equation of torsional creep −(g(|∇u|²)u_x1)_x1−(g(|∇u|²)u_x2)_x2=2?, x∈Ω⊂R², from the torque (or torsional rigidity) T(?), given experimentally. The second class of inverse problems is related to the identification of elastoplastic properties of a 3D body from spherical indentation tests. In this case one needs to determine unknown Lame coefficients in the system of PDEs of non-linear elasticity, from the measured spherical indentation loading curve P=P(α), obtained during the quasi-static indentation test. In the third model an inverse problem of identifying the unknown coefficient g(ξ²(u)) in the non-linear bending equation is analyzed. The boundary measured data here is assumed to be the deflections w_i[τ_k]?w(λ_i;τ_k), measured during the quasi-static bending process, given by the parameter τ_k, , at some points , of a plate. An existence of weak solutions of all direct problems are derived in appropriate Sobolev spaces, by using monotone potential operator theory. Then monotone iteration schemes for all the linearized direct problems are proposed. Strong convergence of solutions of the linearized problems, as well as rates of convergence is proved. Based on obtained continuity property of the direct problem solution with respect to coefficients, and compactness of the set of admissible coefficients, an existence of quasi-solutions of all considered inverse problems is proved. Some numerical results, useful from the points of view of engineering mechanics and computational material science, are demonstrated.     

Multiaxial yield behaviour of Al replicated foam

Multiaxial experiments are performed on replicated aluminium foam using a custom-built apparatus. The foam structure is isotropic, and features open monomodal pores in average diameter. Plane stress (σ₁, σ₂, σ₃=0) and axisymmetric (σ₁, σ₂=σ₃) yield envelopes are measured using cubical specimens, supplemented by tests on hollow cylindrical and uniaxial samples. In addition to the three stress components at 0.2% offset strain, the computer-controlled testing apparatus also measures the three instantaneous displacement vectors.Results show that the shape of the yield surface is independent of the relative density of the foam in the explored range (13-28%). Strain increment vectors lie, within error, roughly normal to the line traced through data points in stress space. Replicated foams feature asymmetric yield behaviour between tension and compression. The data additionally show an influence on the yield surface of the third stress tensor invariant (i.e., of the Lode angle). Simple general expressions for the yield surface are fitted to the data, leading to conclude that their behaviour is slightly better captured by parabolic rather than elliptic expressions dependent on all three stress invariants.     

Geometrical properties of the vorticity vector derived using large-eddy simulation

In this paper, the geometrical properties of the resolved vorticity vector derived from large-eddy simulation are investigated using a statistical method. Numerical tests have been performed based on a turbulent Couette channel flow using three different dynamic linear and nonlinear subgrid-scale stress models. The geometrical properties of have a significant impact on various physical quantities and processes of the flow. To demonstrate, we examined helicity and helical structure, the attitude of with respect to the eigenframes of the resolved strain rate tensor and negative subgrid-scale stress tensor -τ_ij, enstrophy generation, and local vortex stretching and compression. It is observed that the presence of the wall has a strong anisotropic influence on the alignment patterns between and the eigenvectors of , and between and the resolved vortex stretching vector. Some interesting wall-limiting geometrical alignment patterns and probability density distributions in the form of Dirac delta functions associated with these alignment patterns are reported. To quantify the subgrid-scale modelling effects, the attitude of with respect to the eigenframe of -τ_ij is studied, and the geometrical alignment between and the Euler axis is also investigated. The Euler axis and angle for describing the relative rotation between the eigenframes of -τ_ij and are natural invariants of the rotation matrix, and are found to be effective for characterizing a subgrid-scale stress model and for quantifying the associated subgrid-scale modelling effects on the geometrical properties of .     

Global existence and uniqueness of solutions for the equations of third grade fluids

We consider the equations governing the motion of third grade fluids in . We show global existence of solutions without any smallness assumption, by assuming only that the initial velocity belongs to the Sobolev space H². The uniqueness of such solutions is also proven in dimension two.     

Boundary conditions in small-deformation, single-crystal plasticity that account for the Burgers vector

This paper discusses boundary conditions appropriate to a theory of single-crystal plasticity (Gurtin, J. Mech. Phys. Solids 50 (2002) 5) that includes an accounting for the Burgers vector through energetic and dissipative dependences on the tensor G=curlH^p, with H^p the plastic part in the additive decomposition of the displacement gradient into elastic and plastic parts. This theory results in a flow rule in the form of N coupled second-order partial differential equations for the slip-rates , and, consequently, requires higher-order boundary conditions. Motivated by the virtual-power principle in which the external power contains a boundary-integral linear in the slip-rates, hard-slip conditions in which

(A)

on a subsurface S_hard of the boundary

for all slip systems α are proposed. In this paper we develop a theory that is consistent with that of (Gurtin, 2002), but that leads to an external power containing a boundary-integral linear in the tensor , a result that motivates replacing (A) with the microhard condition

(B)

on the subsurface S_hard.

We show that, interestingly, (B) may be interpreted as the requirement that there be no flow of the Burgers vector across S_hard.What is most important, we establish uniqueness for the underlying initial/boundary-value problem associated with (B); since the conditions (A) are generally stronger than the conditions (B), this result indicates lack of existence for problems based on (A). For that reason, the hard-slip conditions (A) would seem inappropriate as boundary conditions.Finally, we discuss conditions at a grain boundary based on the flow of the Burgers vector at and across the boundary surface.     

A gradient theory of small-deformation isotropic plasticity that accounts for the Burgers vector and for dissipation due to plastic spin

This study develops a gradient theory of small-deformation viscoplasticity based on: a system of microforces consistent with its peculiar balance; a mechanical version of the second law that includes, via the microforces, work performed during viscoplastic flow; a constitutive theory that accounts for the Burgers vector through a free energy dependent on , with H^p the plastic part of the elastic-plastic decomposition of the displacement gradient. The microforce balance and the constitutive equations, restricted by the second law, are shown to be together equivalent to a nonlocal flow rule in the form of a coupled pair of second-order partial differential equations. The first of these is an equation for the plastic strain-rate in which the stress T plays a basic role; the second, which is independent of T, is an equation for the plastic spin. A consequence of this second equation is that the plastic spin vanishes identically when the free energy is independent of, but not generally otherwise. A formal discussion based on experience with other gradient theories suggests that sufficiently far from boundaries solutions should not differ appreciably from classical solutions, but close to microscopically hard boundaries, boundary layers characterized by a large Burgers vector and large plastic spin should form.Because of the nonlocal nature of the flow rule, the classical macroscopic boundary conditions need be supplemented by nonstandard boundary conditions associated with viscoplastic flow. As an aid to solution, a variational formulation of the flow rule is derived.Finally, we sketch a generalization of the theory that allows for isotropic hardening resulting from dissipative constitutive dependences on .     

On the statistical analysis of local fracture stresses in notched bars

Test results for critical local fracture stresses are analysed statistically for both “as-received” and “degraded” pressure-vessel weld metal. The values were determined from the fracture loads of blunt-notch four-point-bend specimens fractured over a range of low test temperatures, making use of results from a finite-element stress analysis of the stress-strain distributions ahead of the notch root. The “degraded” material tested in this work has been austenitized at a high temperature, followed by both prestraining and temper embrittlement. This has led to a situation in which the fracture stress for the “degraded” material is reduced significantly below that for the “as-received” material. The fracture mechanisms are different in that the “degraded” material shows evidence of intergranular fracture as well as cleavage fracture (in coarse grain size) whereas the “as-received” material shows only cleavage fracture (in fine grain size). The critical stress (σ_F) distributions plotted on normal probability paper show that the experimental cumulative distribution function (CDF) is linear for each condition with different mean values: for “as-received” material and for “degraded” material. The values of standard deviation are small and almost identical (33-). The decrease of the local fracture stress after degradation is related to the local fracture micro-mechanisms. Statistical analysis of the results for the two conditions supports the hypothesis that the values of σ_F are essentially single valued, within random experimental errors. A similar analysis of the data treating both conditions as a single population reveals some interesting points relating to statistical modelling and lower-bound estimation for mechanical properties. Comparisons are made with Weibull analysis of the data. A further conclusion is that it is extremely important to base any statistical model on inferences drawn from micro-mechanical modelling of processes, and that examination of “normal” CDFs can often provide good indications of when it is necessary to subject data to further statistical and physical analysis.     

A theory of strain-gradient plasticity for isotropic, plastically irrotational materials. Part I: Small deformations

This study develops a small-deformation theory of strain-gradient plasticity for isotropic materials in the absence of plastic rotation. The theory is based on a system of microstresses consistent with a microforce balance; a mechanical version of the second law that includes, via microstresses, work performed during viscoplastic flow; a constitutive theory that allows:

•

the microstresses to depend on , the gradient of the plastic strain-rate, and

•

the free energy ψ to depend on the Burgers tensor .

The microforce balance when augmented by constitutive relations for the microstresses results in a nonlocal flow rule in the form of a tensorial second-order partial differential equation for the plastic strain. The microstresses are strictly dissipative when ψ is independent of the Burgers tensor, but when ψ depends on G the gradient microstress is partially energetic, and this, in turn, leads to a back stress and (hence) to Bauschinger-effects in the flow rule. It is further shown that dependencies of the microstresses on lead to strengthening and weakening effects in the flow rule.Typical macroscopic boundary conditions are supplemented by nonstandard microscopic boundary conditions associated with flow, and, as an aid to numerical solutions, a weak (virtual power) formulation of the nonlocal flow rule is derived.     

Load capacity of bodies

For the stress analysis in a plastic body Ω, we prove that there exists a maximal positive number C, the load capacity ratio, such that the body will not collapse under any external traction field t bounded by CY₀, where Y₀ is the yield stress. The load capacity ratio depends only on the geometry of the body and is given by