4.10. Concrete

The concrete material model predicts the failure of brittle materials. Both cracking and crushing failure modes are accounted for. TB,CONCR accesses this material model, which is available with the reinforced concrete element SOLID65.

The criterion for failure of concrete due to a multiaxial stress state can be expressed in the form (Willam and Warnke([37])):

(4–304)

where:

F = a function (to be discussed) of the principal stress state (σxp, σyp, σzp)
S = failure surface (to be discussed) expressed in terms of principal stresses and five input parameters ft, fc, fcb, f1 and f2 defined in Table 4.4: Concrete Material Table
fc = uniaxial crushing strength
σxp, σyp, σzp = principal stresses in principal directions

If Equation 4–304 is satisfied, the material will crack or crush.

A total of five input strength parameters (each of which can be temperature dependent) are needed to define the failure surface as well as an ambient hydrostatic stress state. These are presented in Table 4.4: Concrete Material Table.

Table 4.4:  Concrete Material Table

(Input on TBDATA Commands with TB,CONCR)
LabelDescriptionConstant
ft Ultimate uniaxial tensile strength3
fc Ultimate uniaxial compressive strength4
fcb Ultimate biaxial compressive strength5
Ambient hydrostatic stress state6
f1 Ultimate compressive strength for a state of biaxial compression superimposed on hydrostatic stress state 7
f2 Ultimate compressive strength for a state of uniaxial compression superimposed on hydrostatic stress state 8

However, the failure surface can be specified with a minimum of two constants, ft and fc. The other three constants default to Willam and Warnke([37]):

(4–305)

(4–306)

(4–307)

However, these default values are valid only for stress states where the condition

(4–308)

(4–309)

is satisfied. Thus condition Equation 4–308 applies to stress situations with a low hydrostatic stress component. All five failure parameters should be specified when a large hydrostatic stress component is expected. If condition Equation 4–308 is not satisfied and the default values shown in Equation 4–305 thru Equation 4–307 are assumed, the strength of the concrete material may be incorrectly evaluated.

When the crushing capability is suppressed with fc = -1.0, the material cracks whenever a principal stress component exceeds ft.

Both the function F and the failure surface S are expressed in terms of principal stresses denoted as σ1, σ2, and σ3 where:

(4–310)

(4–311)

and σ1 σ2 σ3. The failure of concrete is categorized into four domains:

  1. 0 σ1 σ2 σ3 (compression - compression - compression)

  2. σ1 0 σ2 σ3 (tensile - compression - compression)

  3. σ1 σ2 0 σ3 (tensile - tensile - compression)

  4. σ1 σ2 σ3 0 (tensile - tensile - tensile)

In each domain, independent functions describe F and the failure surface S. The four functions describing the general function F are denoted as F1, F2, F3, and F4 while the functions describing S are denoted as S1, S2, S3, and S4. The functions Si (i = 1,4) have the properties that the surface they describe is continuous while the surface gradients are not continuous when any one of the principal stresses changes sign. The surface will be shown in Figure 4.29: 3-D Failure Surface in Principal Stress Space and Figure 4.31: Failure Surface in Principal Stress Space with Nearly Biaxial Stress. These functions are discussed in detail below for each domain.

4.10.1. The Domain (Compression - Compression - Compression)

0 σ1 σ2 σ3

In the compression - compression - compression regime, the failure criterion of Willam and Warnke([37]) is implemented. In this case, F takes the form

(4–312)

and S is defined as

(4–313)

Terms used to define S are:

(4–314)

(4–315)

(4–316)

(4–317)

σh is defined by Equation 4–309 and the undetermined coefficients a0, a1, a2, b0, b1, and b2 are discussed below.

This failure surface is shown as Figure 4.29: 3-D Failure Surface in Principal Stress Space. The angle of similarity η describes the relative magnitudes of the principal stresses. From Equation 4–314, η = 0° refers to any stress state such that σ3 = σ2 > σ1 (e.g. uniaxial compression, biaxial tension) while ξ = 60° for any stress state where σ32 = σ1 (e.g. uniaxial tension, biaxial compression). All other multiaxial stress states have angles of similarity such that 0° η 60°. When η = 0°, S1 Equation 4–313 equals r1 while if η = 60°, S1 equals r2. Therefore, the function r1 represents the failure surface of all stress states with η = 0°. The functions r1, r2 and the angle η are depicted on Figure 4.29: 3-D Failure Surface in Principal Stress Space.

Figure 4.29:  3-D Failure Surface in Principal Stress Space

3-D Failure Surface in Principal Stress Space

It may be seen that the cross-section of the failure plane has cyclic symmetry about each 120° sector of the octahedral plane due to the range 0° < η < 60° of the angle of similitude. The function r1 is determined by adjusting a0, a1, and a2 such that ft, fcb, and f1 all lie on the failure surface. The proper values for these coefficients are determined through solution of the simultaneous equations:

(4–318)

with

(4–319)

The function r2 is calculated by adjusting b0, b1, and b2 to satisfy the conditions:

(4–320)

ξ2 is defined by:

(4–321)

and ξ0 is the positive root of the equation

(4–322)

where a0, a1, and a2 are evaluated by Equation 4–318.

Since the failure surface must remain convex, the ratio r1 / r2 is restricted to the range

(4–323)

although the upper bound is not considered to be restrictive since r1 / r2 < 1 for most materials (Willam([36])). Also, the coefficients a0, a1, a2, b0, b1, and b2 must satisfy the conditions (Willam and Warnke([37])):

(4–324)

(4–325)

Therefore, the failure surface is closed and predicts failure under high hydrostatic pressure (ξ > ξ2). This closure of the failure surface has not been verified experimentally and it has been suggested that a von Mises type cylinder is a more valid failure surface for large compressive σh values (Willam([36])). Consequently, it is recommended that values of f1 and f2 are selected at a hydrostatic stress level in the vicinity of or above the expected maximum hydrostatic stress encountered in the structure.

Equation 4–322 expresses the condition that the failure surface has an apex at ξ = ξ0. A profile of r1 and r2 as a function of ξ is shown in Figure 4.30: A Profile of the Failure Surface.

Figure 4.30:  A Profile of the Failure Surface

A Profile of the Failure Surface

As a Function of ξα


The lower curve represents all stress states such that η = 0° while the upper curve represents stress states such that η = 60°. If the failure criterion is satisfied, the material is assumed to crush.

4.10.2. The Domain (Tension - Compression - Compression)

σ1 0 σ2 σ3

In the regime, F takes the form

(4–326)

and S is defined as

(4–327)

where cos η is defined by Equation 4–314 and

(4–328)

(4–329)

The coefficients a0, a1, a2, b0, b1, b2 are defined by Equation 4–318 and Equation 4–320 while

(4–330)

If the failure criterion is satisfied, cracking occurs in the plane perpendicular to principal stress σ1.

This domain can also crush. See (Willam and Warnke([37])) for details.

4.10.3. The Domain (Tension - Tension - Compression)

σ1 σ2 0 σ3

In the tension - tension - compression regime, F takes the form

(4–331)

and S is defined as

(4–332)

If the failure criterion for both i = 1, 2 is satisfied, cracking occurs in the planes perpendicular to principal stresses σ1 and σ2. If the failure criterion is satisfied only for i = 1, cracking occurs only in the plane perpendicular to principal stress σ1.

This domain can also crush. See (Willam and Warnke([37])) for details.

4.10.4. The Domain (Tension - Tension - Tension)

σ1 σ2 σ3 0

In the tension - tension - tension regimes, F takes the form

(4–333)

and S is defined as

(4–334)

If the failure criterion is satisfied in directions 1, 2, and 3, cracking occurs in the planes perpendicular to principal stresses σ1, σ2, and σ3.

If the failure criterion is satisfied in directions 1 and 2, cracking occurs in the plane perpendicular to principal stresses σ1 and σ2.

If the failure criterion is satisfied only in direction 1, cracking occurs in the plane perpendicular to principal stress σ1.

Figure 4.31:  Failure Surface in Principal Stress Space with Nearly Biaxial Stress

Failure Surface in Principal Stress Space with Nearly Biaxial Stress

Figure 4.31: Failure Surface in Principal Stress Space with Nearly Biaxial Stress represents the 3-D failure surface for states of stress that are biaxial or nearly biaxial. If the most significant nonzero principal stresses are in the σxp and σyp directions, the three surfaces presented are for σzp slightly greater than zero, σzp equal to zero, and σzp slightly less than zero. Although the three surfaces, shown as projections on the σxp - σyp plane, are nearly equivalent and the 3-D failure surface is continuous, the mode of material failure is a function of the sign of σzp. For example, if σxp and σyp are both negative and σzp is slightly positive, cracking would be predicted in a direction perpendicular to the σzp direction. However, if σzp is zero or slightly negative, the material is assumed to crush.


Release 18.2 - © ANSYS, Inc. All rights reserved.