On criticism of the nature of objectivity in classical continuum physics

Murdoch (J. Elasticity 60, 233-242, 2000) showed that restrictions imposed upon response functions by material frame-indifference are the consequences of five distinct aspects of observer agreement (that is, of objectivity) and involve only proper orthogonal tensors. Accordingly it is unnecessary to invoke the principle of invariance under superposed rigid motions (in the sense of one observer, two motions), which imposes a restriction upon nature. Liu (Continuum Mech. Thermodyn. 16, 177-183, 2003, and Continuum Mech. Thermodyn. 17, 125-133, 2005) has challenged, misinterpreted and misrepresented the content of both Murdochs work and this work. Here all criticisms of Liu are answered, his counter-examples are used to amplify the tenets of Murdochs work, and a key modelling issue in the controversy is indicated. Further, the response function restrictions for a given observer, derived on the basis of considering other observers, are shown to be independent of possible differences in the scales of mass, length, and time employed by other observers.Received: 2 March 2004, Accepted: 20 August 2004, Published online: 24 March 2005PACS: 83.10Ff     

Effect of material nonhomogeneities on the HRR dominance

This work studies the asymptotic stress and displacement fields near the tip of a stationary crack in an elastic–plastic nonhomogeneous material with the emphasis on the effect of material nonhomogeneities on the dominance of the crack tip field. While the HRR singular field still prevails near the crack tip if the material properties are continuous and piecewise continuously differentiable, a simple asymptotic analysis shows that the size of the HRR dominance zone decreases with increasing magnitude of material property gradients. The HRR field dominates at points that satisfy |α⁻¹ ∂α/∂x_δ|1/r, |α⁻¹ ∂²α/(∂x_δ ∂x_γ)|1/r², |n⁻¹ ∂n/∂x_δ|1/[r|ln(r/A)|] and |n⁻¹ ∂²n/(∂x_δ ∂x_γ)|1/[r²|ln(r/A)|], in addition to other general requirements for asymptotic solutions, where α is a material property in the Ramberg–Osgood model, n is the strain hardening exponent, r is the distance from the crack tip, x_δ are Cartesian coordinates, and A is a length parameter. For linear hardening materials, the crack tip field dominates at points that satisfy |E_tan⁻¹ ∂E_tan/∂x_δ|1/r, |E_tan⁻¹ ∂²E_tan/(∂x_δ ∂x_γ)|1/r², |E⁻¹ ∂E/∂x_δ|1/r, and |E⁻¹ ∂²E/(∂x_δ ∂x_γ)|1/r², where E_tan is the tangent modulus and E is Young’s modulus.     

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.     

Higher-order stress and grain size effects due to self-energy of geometrically necessary dislocations

The higher-order stress work-conjugate to slip gradient in single crystals at small strains is derived based on the self-energy of geometrically necessary dislocations (GNDs). It is shown that this higher-order stress changes stepwise as a function of in-plane slip gradient and therefore significantly influences the onset of initial yielding in polycrystals. The higher-order stress based on the self-energy of GNDs is then incorporated into the strain gradient plasticity theory of Gurtin [2002. A gradient theory of single-crystal viscoplasticity that accounts for geometrically necessary dislocations. J. Mech. Phys. Solids 50, 5-32] and applied to single-slip-oriented 2D and 3D model crystal grains of size D. It is thus found that the self-energy of GNDs gives a D^-1-dependent term for the averaged resolved shear stress in such a model grain under yielding. Using published experimental data for several polycrystalline metals, it is demonstrated that the D^-1-dependent term successfully explains the grain size dependence of initial yield stress and the dislocation cell size dependence of flow stress in the submicron to several-micron range of grain and cell sizes.     

Experimental investigation of initial and subsequent yield surfaces for laminated metal matrix composites

The aim of this work is to construct yield surfaces to describe initial yielding and characterize hardening behavior of a highly anisotropic material. A methodology for constructing yield surfaces for isotropic materials using axial–torsion loading is extended to highly anisotropic materials. The technique uses a sensitive definition of yielding based on permanent strain rather than offset strain, and enables multiple yield points and multiple yield surfaces to be conducted on a single specimen. A target value of 20 × 10⁻⁶ is used for Al₂O₃ fiber reinforced aluminum laminates having a fiber volume fraction of 0.55. Sixteen radial probes are used to define the yield locus in the axial–shear stress plane. Initial yield surfaces for [0₄], [90₄], and [0/90]₂ fibrous aluminum laminates are well described by ellipses in the axial–shear stress plane having aspect ratios of 10, 2.5, and 3.3, respectively. For reference, the aspect ratio of the Mises ellipse for an isotropic material is 1.73. Initial yield surfaces do not have a tension–compression asymmetry. Four overload profiles (plus, ex, hourglass, and zee) are applied to characterize hardening of a [0/90]₂ laminate by constructing 30 subsequent yield surfaces. Parameters to describe the center and axes of an ellipse are regressed to the yield points. The results clearly indicate that kinematic hardening dominates so that material state evolution can be described by tracking the center of the yield locus. For a nonproportional overload of (σ, τ) = (500, 70) MPa, the center of the yield locus translated to (σ, τ) = (430, 37) MPa and the ellipse major axis was only 110 MPa.     

Elasticity and strength of partially sintered ceramics

A discrete element model for the elastic and fracture behavior of partially sintered ceramics is presented. It accounts for the granular character of the material when a large amount of porosity (typically >0.2-0.4) is left after sintering. The model uses elastic force-displacement laws to represent the bond formed between particles during sintering. Bond fracture in tension and shearing is accounted for in the model. Realistic numerical microstructures are generated using a sintering model on random particle packings. In particular, packings with fugitive pore formers are used to create partially sintered microstructures with large pores. The effective elastic response and the strength of these microstructures are calculated in tension and compression. The link between important microstructural features such as bond size or coordination number and macroscopic behavior is investigated. In particular, it is shown that porosity alone is not sufficient to account for the mechanical properties of a partially sintered material.     

On the crack-tip stress singularity of interfacial cracks in transversely isotropic piezoelectric bimaterials

In this paper, characteristics of the interface crack-tip stress and electric displacement fields in transversely isotropic piezoelectric bimaterials are studied. The authors have proven, within the framework of the generalized Stroh formalism for piezoelectric bimaterials, that there is no coexistence of the parameters (oscillating) and κ (non-oscillating) in the interface crack-tip generalized stress field for all transversely isotropic piezoelectric bimaterials. This leads to the classification of piezoelectric bimaterials into one group that exhibits the oscillating property in the interface crack-tip generalized stress field and the other that does not. Fifteen (15) pair-combinations of six (6) piezoelectric materials PZT-4, PZT-5H, PZT-6B, PZT-7A, P-7, and BaTiO₃, which are commonly used in practice, are numerically analyzed in this study, and the results backup the above theoretical conclusions. Moreover, the associated eigenvectors for such material systems (with either =0 or κ=0) are also obtained numerically, and the result show that there still exist four linear independent associate eigenvectors for each bimaterial.     

On the r-value of textured sheet metals

The r-value of a sheet metal is a measure of plastic anisotropy frequently used for prediction of performance in deep-drawing. It has also figured prominently in the literature for validation of theories where the predicted angular dependence of r is compared with the measured dependence. As plastic anisotropy in sheet metals is caused mainly by the preferred orientations of grains within the polycrystalline metal, it is natural to ask how r would depend on the orientation distribution function (ODF) w which defines the crystallographic texture of the polycrystal. In this paper a general formula relating r to w is derived for textured sheet metals whose plastic flow behavior is governed by a plastic potential f(σ, w), the anisotropic part of which depends linearly on the texture coefficients; here σ denotes the deviator of the Cauchy stress. Specific forms of this formula for orthorhombic sheets of cubic and of hexagonal metals are explicitly given.     

Homogenization and Global Responses of Inhomogeneous Spherical Nonlinear Elastic Shells

Homogenization of radially inhomogeneous spherical nonlinear elastic shells subject to internal pressure is studied. The equivalent homogeneous material is defined in such a way that it gives rise to exactly the same global response to the pressure load as that of the inhomogeneous shell. For a shell with general strain–energy function and inhomogeniety, the strain–energy function of the equivalent homogeneous material is determined explicitly. The resulting formula is used to study layered composite shells. The equivalent homogeneous material for an infinitely fine layered composite shell is examined, and is found to give not only the same global response, but also the same average stress field as the composite shell does.     

Bounds on the effective conductivity of statistically isotropic multicomponent materials and random cell polycrystals

Upper and lower bounds on the effective conductivity of statistically isotropic multicomponent materials in d dimensions (d=2 or 3) are constructed from the minimum energy principles and appropriate trial fields. The trial fields, involving harmonic potentials and free parameters to be optimized, lead to the bounds containing up to three-point correlation information about the microgeometry of a composite. The bounds are applied to give estimates for the symmetric cell materials, which are optimal over some ranges of parameters, and asymmetric multicoated spheres, which yield the exact effective conductivity in certain cases. The results also agree with many known ones. New bounds for random cell polycrystals are obtained and illustrated on a number of polycrystalline aggregates.     

On the effective stress in unsaturated porous continua with double porosity

Using mixture theory we formulate the balance laws for unsaturated porous media composed of a double-porosity solid matrix infiltrated by liquid and gas. In this context, the term ‘double porosity’ pertains to the microstructural characteristic that allows the pore spaces in a continuum to be classified into two pore subspaces. We use the first law of thermodynamics to identify energy-conjugate variables and derive an expression for the ‘effective’, or constitutive, stress that is energy-conjugate to the rate of deformation of the solid matrix. The effective stress has the form , where σ is the total Cauchy stress tensor, B is the Biot coefficient, and is the mean fluid pressure weighted according to the local degrees of saturation and pore fractions. We identify other emerging energy-conjugate pairs relevant for constitutive modeling of double-porosity unsaturated continua, including the local suction versus degree of saturation pair and the pore volume fraction versus weighted pore pressure difference pair. Finally, we use the second law of thermodynamics to determine conditions for maximum plastic dissipation in the regime of inelastic deformation for the unsaturated two-porosity mixture.     

Deformation of thermosetting resins at impact rates of strain. Part 2: constitutive model with rejuvenation

The constitutive responses of three glassy thermoset polymers at impact rates of strain and slower, together with measurements of adiabatic heating, were reported earlier by the authors. The results are interpreted here in the context of a constitutive model proposed previously for amorphous polymers, expanded to incorporate strain-softening and the adiabatic heating deficit. In terms of the model, both features are a natural consequence of strain-induced evolution of the glass structure, as represented by Tool's “fictive temperature”—the phenomenon of structural rejuvenation. A representation is proposed for the evolution of fictive temperature with plastic strain, motivated by an approximate treatment of the kinetics of physical ageing/rejuvenation. Formulated in this manner, the model agrees reasonably well with experimental results across the wide range of strain rates of the previous experiments, 10⁻³ to , and across most of the range of strain to failure in compression. At the highest strains, however, an additional adiabatic heating deficit appears that is not predicted by the model, either suggesting the onset of structural breakdown possibly associated with the appearance of cracks or reflecting a need for better physical understanding of large deformations in glassy polymers.     

Cyclic behaviour of short glass fibre reinforced polyamide: Experimental study and constitutive equations

Polymer matrix composites are widely used in the automotive industry and undergo fatigue loadings. The investigation of the nonlinear cyclic behaviour of such materials is a required preliminary work for a confident fatigue design, but has not involved many publications in the literature. This paper presents an extensive experimental study conducted on a polyamide 66 reinforced with 35 wt% of short glass fibres (PA66 GF35), at room temperature. The material was tested in two conditions: dry-as-moulded (DAM) and at the equilibrium with air containing 50% of relative humidity (RH50).An exhaustive experimental campaign in tensile mode has been carried out, including various strain or stress rates, complex mechanical histories and local thermo-mechanical recordings. Such an extended database allowed us to highlight several complex physical phenomena: viscoelastic effects at different time scales, irrecoverable mechanisms, non-linear kinematic hardening, non-linear viscous flow rule, cyclic softening.Taking into account this advanced analysis, a constitutive model describing the cyclic behaviour is proposed. As the experimental database only includes uniaxial tensile tests, the general 3D anisotropic frame is reduced to an uniaxial model valid for a specific orientation distribution. The robust identification process is based on tests which enable the uncoupling between the underlined mechanical features. This strategy leads to a model which accurately predicts the cyclic behaviour of conditioned as well as dry materials under complex tensile loadings.     

Some asymptotic results concerning the buckling of a spherical shell of arbitrary thickness

For a spherical shell of arbitrary thickness which is subjected to an external hydrostatic pressure, symmetrical buckling takes place at a value of μ₁ which depends on and the mode number, where A₁ and A₂ are the undeformed inner and outer radii, and μ₁ is the ratio of the deformed inner radius to the undeformed inner radius. In the large mode number limit, we find that the dependence of μ₁ on has a boundary layer structure: it is a constant over almost the entire region of and decreases sharply from this constant value to unity as tends to unity (the thin-shell limit). Simple asymptotic expressions for the bifurcation condition are obtained. The classical result for thin shells is recovered directly from the equations of finite elasticity, and an asymptotic critical neutral curve (which envelops the neutral curves corresponding to different mode numbers) is obtained.     

Traction discontinuity problem in 3D BEM of the transversely isotropic material

The so-called edge problem, i.e., the traction discontinuity at edges or corners of the domain, is solved for the case of the transversely isotropic material being introduced into 3D BEM. The auxiliary relationship equations for the additional traction components are derived and the arbitrary inclination of the isotropic plane of the material is taken into account in the mean time. An example is given, and the comparison of the numerical result of the example with the theoritical solution is also illustrated.     

Nonlinear elasto-plastic model for dense granular flow

This work proposes a model for granular deformation that predicts the stress and velocity profiles in well-developed dense granular flows. Recent models for granular elasticity [Jiang, Y., Liu, M., 2003. Granular elasticity without the Coulomb condition. Phys. Rev. Lett. 91, 144301] and rate-sensitive fluid-like flow [Jop, P., Forterre, Y., Pouliquen, O., 2006. A constitutive law for dense granular flows. Nature 441, 727] are reformulated and combined into one universal elasto-plastic law, capable of predicting flowing regions and stagnant zones simultaneously in any arbitrary 3D flow geometry. The unification is performed by justifying and implementing a Kröner–Lee decomposition, with care taken to ensure certain continuum physical principles are necessarily upheld. The model is then numerically implemented in multiple geometries and results are compared to experiments and discrete simulations.     

Preparation of porous silica spheres with self-dispersing properties

A sol-gel procedure in a water/oil emulsion was introduced for the synthesis of porous silica spheres. Tetraethoxysilane was used as the silica source. The specific surface area and total pore volume of the product reached 772.3 m2/g and 0.663 cm3/g, respectively. The electrolyte washing process conferred a surface charge to the product, which displayed self-dispersal properties in water. The porous spheres have potential applications in the fields of drug delivery, controlled release capsules, indoor air pollutant scavengers, and hydrogen storage agents. The oil phase, which accounts for over 8O% of the chemical cost of the procedure, could largely be recycled by filtering, standing, and layering. The whole procedure is suitable for application as an industrial process.     

Presentation of the creep curves of solid biological materials by a simplified mathematical version of the generalized Kelvin-Voigt model

Published creep equations, derived from the Kelvin-Voigt model with a discrete retardation time spectrum were used to generate data for linear presentation of the creep function. It is shown that 2–3 term creep functions, containing 4–6 constants, could be reduced to an algebraic expression having two constants only. Consequently, entire creep curves of a variety of biological and food materials could be described by a single type of equation in the form of:P (t) = k ₀ + k₁ t + t/(k₂ + k₃t) whereP (t) is the creep parameter (compliance, strain or deformation),k ₀ a constant representing the instantaneous compliance,k ₁ a constant representing the steady-state flow, andk ₂ andk ₃ the characteristic constants of the creep function (t) obtained by regression from the linear relationship oft/ (t) vst.