Bigoni–Piccolroaz yield criterion

From HandWiki
Short description: Yield criterion

The Bigoni–Piccolroaz yield criterion is a yielding model, based on a phenomenological approach, capable of describing the mechanical behavior of a broad class of pressure-sensitive granular materials such as soil, concrete, porous metals and ceramics.

General concepts

The hardening process during hydrostatic powder compaction described with the Bigoni & Piccolroaz yield surface.

The idea behind the Bigoni-Piccolroaz criterion is that of deriving a function capable of transitioning between the yield surfaces typical of different classes of materials only by changing the function parameters.[1] The reason for this kind of implementation lies in the fact that the materials towards which the model is targeted undergo consistent changes during manufacturing and working conditions. The typical example is that of the hardening of a power specimen by compaction and sintering during which the material changes from granular to dense.

The Bigoni-Piccolroaz yielding criterion can be represented in the Haigh–Westergaard stress space as a convex smooth surface and in fact the criterion itself is based on the mathematical definition of the surface in the above-mentioned space as a proper interpolation of experimental points.

Mathematical formulation

Bigoni-Piccolroaz yield surface in the 3D space of principal stresses

The Bigoni-Piccolroaz yield surface is thought as a direct interpolation of experimental data. This criterion represents a smooth and convex surface, which is closed both in hydrostatic tension and compression and has a drop-like shape, particularly suited to describe frictional and granular materials. This criterion has also been generalized to the case of surfaces with corners.[2]

Design principles

Since the whole idea of the model is to tailor a function to experimental data, the authors have defined a certain group of features as desirable, even if not essential, among those:

  • smoothness of the surface;
  • possibility of changing the shape and thus the interpolation on a broad class of experimental data for different materials;
  • possibility to represent known criteria with limit set of parameters;
  • convexity of the surface.

Parametric function

Representation of coordinates [math]\displaystyle{ p, q }[/math] and [math]\displaystyle{ \theta }[/math] of the BP yield surface in Haigh Westergaard stress space.
Deviatoric sections of a BP yield surface.
Meridian sections of a BP yield surface.

The Bigoni–Piccolroaz yield criterion is a seven-parameter surface defined as:

[math]\displaystyle{ f(p,q,\theta) = F(p) + \frac{q}{g(\theta)} = 0, }[/math]

where p, q and [math]\displaystyle{ \theta }[/math] are invariants dependent on the stress tensor, while [math]\displaystyle{ F(p) }[/math] is the "meridian" function:

[math]\displaystyle{ F(p) = \left\{ \begin{array}{ll} -M p_c \sqrt{(\phi - \phi^m)[2(1 - \alpha)\phi + \alpha]}, & \phi \in [0,1], \\ +\infty, & \phi \notin [0,1], \end{array} \right. }[/math]
[math]\displaystyle{ \phi = \frac{p + c}{p_c + c}, }[/math]

describing the pressure-sensitivity and [math]\displaystyle{ g(\theta) }[/math] is the "deviatoric" function:[3]

[math]\displaystyle{ g(\theta) = \frac{1}{\cos[\beta \frac{\pi}{6} - \frac{1}{3} \cos^{-1}(\gamma \cos 3\theta)]}, }[/math]

describing the Lode-dependence of yielding.

The mathematical definitions of the parameters [math]\displaystyle{ p, q }[/math] and [math]\displaystyle{ \theta }[/math] are:

[math]\displaystyle{ p = -\frac{tr ~ \sigma }{3}, \quad q = \sqrt{3\textit{J}_2}, \quad \theta = \frac{1}{3} cos^{-1} \left( \frac{3\sqrt{3}}{2} \frac{J_3}{J_2^{3/2}} \right) }[/math]

With:

[math]\displaystyle{ J_2 = \frac{1}{2} \textbf{S} \cdot \textbf{S}, \quad \textbf{S} = \sigma - \frac{tr \ \sigma}{3} \textbf{I}, \quad J_3 = \frac{1}{3} tr \ \textbf{S}^3 }[/math]

Where [math]\displaystyle{ \textbf{S} }[/math] is the deviatoric stress, [math]\displaystyle{ \textbf{I} }[/math] is the identity tensor, [math]\displaystyle{ \sigma }[/math] is the stress tensor and the dot indicates the scalar product.

A better understanding of those important parameters can be grasped by using their geometrical representation in the Haigh–Westergaard stress space.[4] Considering the tern of principal stresses and the deviatoric plane [math]\displaystyle{ \pi }[/math], orthogonal to the trisector of the first quadrant and passing through the origin of the coordinate system, the tern [math]\displaystyle{ p, q }[/math] and [math]\displaystyle{ \theta }[/math] unequivocally represents a point in the space acting as a cylindrical coordinate system with the trisector as an axis:

  • [math]\displaystyle{ - \sqrt{3} ~p }[/math] is the distance of the point from the deviatoric plane [math]\displaystyle{ \pi }[/math];
  • [math]\displaystyle{ \sqrt{2/3}~ q }[/math] is the distance from the trisector;
  • [math]\displaystyle{ \theta }[/math] represents the angle between the projections of [math]\displaystyle{ q }[/math] and the axis [math]\displaystyle{ \sigma_1 }[/math] on the deviatoric plane [math]\displaystyle{ \pi }[/math].

The usage of p and q instead of the correct cylindrical coordinates [math]\displaystyle{ \xi }[/math] and [math]\displaystyle{ \rho }[/math]:

[math]\displaystyle{ \xi = \frac{1}{\sqrt{3}} tr(\sigma)= - \sqrt{3} ~p, \quad \rho = \sqrt{2 ~ J_2} =\sqrt{\frac{2}{3}} q }[/math]

is justified by the easier physical interpretation:[5] p is the hydro-static pressure on the material point, q is the Von Mises equivalent stress.

The described yield function corresponds to the yield surface:

[math]\displaystyle{ q = -f(p) g(\theta), \quad p \in [-c,p_c], \quad \theta \in [0,\pi/3] }[/math]

which makes explicit the relation between the two functions [math]\displaystyle{ f(p) }[/math] and [math]\displaystyle{ g(\theta) }[/math] and the shape of the meridian and deviatoric sections, respectively.

The seven, non-negative material parameters:

[math]\displaystyle{ \underbrace{M \gt 0,~ p_c \gt 0,~ c \geq 0,~ 0 \lt \alpha \lt 2,~ m \gt 1}_{\mbox{defining}~\displaystyle{F(p)}},~~~ \underbrace{0\leq \beta \leq 2,~ 0 \leq \gamma \lt 1}_{\mbox{defining}~\displaystyle{g(\theta)}}, }[/math]

define the shape of the meridian and deviatoric sections. In particular, some of the parameters are easily relatable to mechanical properties: [math]\displaystyle{ M }[/math] controls the pressure sensitivity, [math]\displaystyle{ c }[/math] and [math]\displaystyle{ p_c }[/math] are the yield strength under isotropic conditions of tension and compression. The other parameters define the shape of the surface when intersected by the meridian and deviatoric planes: [math]\displaystyle{ \alpha }[/math] and [math]\displaystyle{ m }[/math] define the meridian section, [math]\displaystyle{ \beta }[/math] and [math]\displaystyle{ \gamma }[/math] define the deviatoric section.

Related yielding criteria

Having been designed to allow consistent changes in the surface shape in the Haigh–Westergaard stress space, the Bigoni-Piccolroaz yield surface can be used as a generalized formulation for several criteria,[6] such as the well known von Mises, Tresca, Mohr–Coulomb.

See also

External links

The Bigoni-Piccolroaz yield surface is a powerful instrument for the characterization of granular materials and it arises great interest in the field of the definition of constitutive models for ceramics, rock and soil which is a task of fundamental importance for better design of products using these materials.

References

  1. Bigoni, D. and Piccolroaz, A., (2004), Yield criteria for quasibrittle and frictional materials, International Journal of Solids and Structures 41, 2855–2878.
  2. Piccolroaz, A. and Bigoni, D. (2009), Yield criteria for quasibrittle and frictional materials: a generalization to surfaces with corners, International Journal of Solids and Structures 46, 3587–3596.
  3. Podgórski, J. (1984). Limit state condition and the dissipation function for isotropic materials. Archives of Mechanics, 36 (3), pp. 323–342.
  4. Niels Saabye Ottosen, Matti Ristinmaa 2005 The Mechanics of constitutive modeling, Elsevier pp. 150-153
  5. Neto Peric Owen 2008 Computational methods for plasticity, Wiley pp. 157-165
  6. Bigoni, D. Nonlinear Solid Mechanics: Bifurcation Theory and Material Instability. Cambridge University Press, 2012 . ISBN:9781107025417.