10x
005925
2025-01-21

TS according to Quast and modified Quast

Local Anisotropic Tension-Stiffening Augmentation

When using a concrete damage model, tension stiffening is represented by the local anisotropic tension-stiffening augmentation [4]. With an anisotropic damage model, this method replaces the previous, globally acting approach: instead of describing the tension zone of the cross-section once via a single stress-strain curve valid for the entire cross-section, it superimposes a reinforcement-induced tensile stress on the basic material response of the concrete at each concrete integration point. The correction is calculated locally (per integration point) and anisotropically (separately in both principal strain directions); the angle between the reinforcement direction and the respective principal tension direction is explicitly taken into account.

Info

Scope by element type: For surface elements (shells, slabs) with a layered cross-section, the evaluation is performed fully anisotropically as described, separately in both principal strain directions. The method is also available for member elements, but with limitations regarding the angle reference: since only one governing (axial) strain direction exists, the distinction between two principal directions—and thus the angle factor aj,k between reinforcement and principal tension direction in the true sense—does not apply; the reinforcement is treated as parallel to the member axis. For solid elements, the method is not available. Likewise, the explicit reinforcement simulation using the member type "reinforcement bar" is currently not included; only the longitudinal reinforcement in surface reinforcement and member reinforcement is considered.

Info

Scope by material model: When an anisotropic concrete damage model is used ("Anisotropic Damage"), tension stiffening is represented via the local anisotropic tension-stiffening augmentation as described above; for this combination, this replaces the previous, globally acting approach. For the isotropic damage model ("Isotropic Damage", Isotropic Damage), instead of the angle factor aj,k, a projection of the two current principal tensile strains onto the reinforcement direction is used, with the same geometric weights cos2(Δθ) or sin2(Δθ); the result is normalized to the larger of the two principal tensile strains and yields a single activation value that is used equally for both principal directions. The angle factor is thus not eliminated entirely but converted into a strain-based projection, since an isotropic damage model itself provides no information about which direction has actually cracked. The advantages of the new approach, including the improved convergence behavior, are retained.

Activation Factor

At each concrete integration point, the principal strains ε1, ε2 and the associated principal strain angle θp are determined; only tension directions (εk > 0) are relevant for tension stiffening. Each reinforcement layer j contributes to the activation factor αk of the principal direction k via three multiplicative partial factors:

  • gj (capacity factor): decreases with increasing reinforcement strain and becomes zero when the yield strain of the reinforcement is reached.
  • aj,k (angle factor): for the anisotropic damage model, 1 for parallel and 0 for perpendicular alignment of reinforcement and principal tension direction, and the complement for the respective other principal direction. For the isotropic damage model, a different calculation applies without complement formation, yielding the same value for both principal directions, which can still be non-zero at exactly 90° under biaxial tension; see the banner "Scope by material model" above.
  • pj (position factor): 1 within the concrete cover zone of the considered reinforcement layer, decreasing outside (strain compatibility between concrete and reinforcement).

The contributions of all reinforcement layers are combined via a smoothed p-norm aggregator and limited to [0, 1]:

Info

If there is no reinforcement, or no direction-compatible reinforcement, in the governing principal tension direction (αk = 0), the tension-stiffening correction is omitted entirely for this direction; the basic material then provides the unchanged tensile behavior.

Target Stress Curve: Two Variants

The target stress in principal direction k is obtained as:

Two variants are available for the underlying target stress curve σ̂(ε), based on the two following subsections:

Target stress curve: two variants
Variant Target stress curve based on Behavior Recommendation
according to Quast Section "Quast" enforcing: can both add and remove tensile stress in order to drive the mean response to the target value for calibration to a specified mean tensile stress-strain behavior
according to Quast, modified Section "Modified Approach According to Quast" purely additive (Δσk ≥ 0): adds tensile stress only if the target value exceeds the current concrete stress default setting (numerically more robust)

Notes for Existing Models

Important

Deviations from older program versions are intentional. The local anisotropic tension-stiffening augmentation completely replaces the previous, globally acting approach; this used the maximum reinforcement strain across all reinforcement layers as the sole activation driver for the entire cross-section, regardless of the direction and position of the reinforcement, which could lead to physically unjustified activation, for example on the compression side of a cross-section subjected to bending. The new augmentation specifically suppresses activation where it is physically unjustified. Deviating results compared to older calculations are therefore to be expected, especially for multi-layer or multi-directional reinforcement as well as for shrinkage-dominated states. For simple bending without shrinkage, both approaches provide qualitatively comparable results. In addition, the method does not model discrete cracks, crack spacings, or crack width; crack width determination continues to be performed via a separate, code-compliant verification step.

Quast

This model for capturing the contribution of the concrete in tension between the cracks is based on a defined stress-strain curve of the concrete in the tension zone (parabola-rectangle diagram).

The basic assumptions of the approach by Quast can be summarized as follows:

  • Full contribution of the concrete in tension up to reaching the crack strain εcr or the calculated concrete tensile strength fct,R.
  • Reduced stiffening contribution of the concrete in the tension zone according to the existing concrete strain.
  • No application of tension stiffening after the yielding of the governing reinforcement bar.

In summary, this means that the calculated tensile strength fct,R is not a fixed value but relates to the existing strain in the governing steel (tension) fiber. The maximum tensile strength fct,R decreases linearly to zero from the defined crack strain εcr until the yield strain of the reinforcement steel is reached in the governing steel fiber. This is achieved by the stress-strain curve shown in the image below in the tension zone of the concrete (parabola-rectangle diagram) and the determination of a reduction factor VMB (stiffening contribution of the concrete).

The next image schematically shows the stress states for increasing loading with tension stiffening.

The stress-strain curve in the tension zone can be described by the following equations:

for 0 < ε < εcr for ε > εcr

  • The fullness of the parabola in the first section can be controlled by the exponent nPR.
  • The exponent should be adjusted such that the transition from the compression to the tension zone occurs with the same modulus of elasticity as much as possible.

To determine the reduction factor VMB, the strain at the most heavily tensioned steel fiber is used. The position of the reference point is illustrated in the following graphic.

The reduction factor VMB decreases with increasing steel strain. In the diagram for the factor VMB (see image below), it can be seen that the factor VMB is reduced exactly to zero at the onset of yielding of the reinforcement.

The curve of the reduction factor VMB in State II (ε > εcr) can be controlled via the exponent nVMB.

  • According to Pfeiffer [2], values of nVMB = 1 (linear) to nVMB = 2 (parabola) are empirical values for components subjected to bending stress.
  • Quast [3] uses the exponent nVMB = 1 (linear) in his model and thereby achieves good agreement in the recalculation of column tests.
  • According to Pfeiffer [2], pure tension tests can be represented with nVMB = 2 with acceptable agreement.

The assumption of a parabola-rectangle diagram for the cracked concrete tension zone is to be regarded as a computational aid. At first glance, there are large differences compared to experimentally determined stress-strain curves on the tension side of plain concrete.

When considering the existing stresses in the reinforced concrete cross-section subjected to bending stress, it becomes apparent that the parabola-rectangle diagram is indeed better able to describe the strains and stresses on average. In a bending beam, a concrete body forms between two cracks. This acts as a kind of wall into which tensile forces are gradually re-introduced by the reinforcement. As a result, a very irregular distribution of stress and strain develops. On average, however, a strain plane with a parabola-rectangle course can be constructed, by means of which the mean curvature can be captured.

For the Quast model, the design values to be used were proposed as follows

  • for the tensile strength fct,R
  • for the crack strain εcr,R

The design value for the tensile strength fct,R is thus smaller than the specifications of the Eurocode. This is justified by the description of the stress-strain curve and the determination of the reduction factor VMB, in which the assumed tensile stress and the resulting tensile force are only gradually reduced after the tensile strain is exceeded. At a strain of 2 ⋅ εcr, an effective tensile stress of about 0.95 ⋅ fct,R still results. Thus, the reduction of the stiffness can be well predicted for bending stress. For pure tensile stress, the values given above for fct,R are too low. Here, according to Pfeiffer [2], the values from EC 2 should be used for the design value of the tensile strength. The values recommended by Quast [3] for fct,R = 1/20 ⋅ fcm can be achieved by assuming 60% of the tensile strengths specified in EC 2. When assuming fct,R = 0.6 ⋅ fctm, on the one hand, cracking of the cross-section is predicted too early. On the other hand, however, a reduction of the tensile strength under permanent load (approx. 70%) or a temporarily higher load (e.g., the short-term application of the rare action combination), which leads to a damaged tension zone, is already taken into account. The individual design values for the tension zone of the concrete can be described as follows:

Modified Approach According to Quast

The modified approach according to Quast uses the same definition of the calculated tensile strength fct,R and the crack strain εcr,R as the approach according to Quast (see previous section). The difference lies in the manner of application to the concrete stress, not in the target curve itself:

  • With the approach according to Quast, the concrete stress in the tension zone is directly set to the value of the target curvec = σ̂(ε), see target curve in section "Quast"); compared to the stress of the unmodified basic material model, this can mean both an increase and a reduction.
  • With the modified approach, the target curve is instead purely additively superimposed on the concrete stress σc,Basis(ε) provided by the basic material model (e.g., damage model).

The concrete stress can thus only be increased by the tension stiffening, never reduced. This additive procedure is numerically more robust than the enforcing variant according to Quast, since it does not lead to an abrupt stiffness jump when the basic material stress already exceeds the target curve.

Info

The modified approach according to Quast provides the additive target stress curve of the local anisotropic tension-stiffening augmentation (see chapter introduction) and is recommended there as the default setting.


References
Parent Chapter