SOLID65 allows the presence of four different materials within each element; one matrix material (e.g. concrete) and a maximum of three independent reinforcing materials. The concrete material is capable of directional integration point cracking and crushing besides incorporating plastic and creep behavior. The reinforcement (which also incorporates creep and plasticity) has uniaxial stiffness only and is assumed to be smeared throughout the element. Directional orientation is accomplished through user specified angles.

The stress-strain matrix [D] used for this element is defined as:

(1–99)

where:

Nr = number of reinforcing materials (maximum of three, all reinforcement is ignored if M1 is zero. Also, if M1, M2, or M3 equals the concrete material number, the reinforcement with that material number is ignored)

[Dc] = stress-strain matrix for concrete, defined by Equation 1–100

[Dr]i = stress-strain matrix for reinforcement i, defined by Equation 1–101

M1, M2, M3 = material numbers associated of reinforcement (input as MAT1, MAT2, and MAT3 on R command)

The matrix [D^c] is derived by specializing and inverting the orthotropic stress-strain relations defined by Equation 2–4 to the case of an isotropic material or

(1–100)

where:

The orientation of the reinforcement i within an element is depicted in Figure 1.10: Reinforcement Orientation. The element coordinate system is denoted by (X, Y, Z) and describes the coordinate system for reinforcement type i. The stress-strain matrix with respect to each coordinate system has the form

(1–101)

where:

It may be seen that the only nonzero stress component is , the axial stress in the direction of reinforcement type i. The reinforcement direction is related to element coordinates X, Y, Z through

(1–102)

where:

θi = angle between the projection of the axis on XY plane and the X axis (input as THETA1, THETA2, and THETA3 on R command)

φi = angle between the axis and the XY plane (input as PHI1, PHI2, and PHI3 on R command)

= direction cosines between axis and element X, Y, Z axes

Figure 1.10: Reinforcement Orientation

Since the reinforcement material matrix is defined in coordinates aligned in the direction of reinforcement orientation, it is necessary to construct a transformation of the form

(1–103)

in order to express the material behavior of the reinforcement in global coordinates. The form of this transformation by Schnobrich is:

(1–104)

where the coefficients a_ij are defined as

(1–105)

The vector is defined by Equation 1–102 while and are unit vectors mutually orthogonal to thus defining a Cartesian coordinate referring to reinforcement directions. If the operations presented by Equation 1–103 are performed substituting Equation 1–101 and Equation 1–104, the resulting reinforcement material matrix in element coordinates takes the form

(1–106)

where:

Therefore, the only direction cosines used in [D^R]_i involve the uniquely defined unit vector .

The matrix material (concrete) is capable of plasticity, creep, cracking and crushing. The concrete material model with its cracking and crushing capabilities is discussed in Concrete (with SOLID65). This material model predicts either elastic behavior, cracking behavior or crushing behavior. If elastic behavior is predicted, the concrete is treated as a linear elastic material (discussed above). If cracking or crushing behavior is predicted, the elastic, stress-strain matrix is adjusted as discussed below for each failure mode.

The presence of a crack at an integration point is represented through modification of the stress-strain relations by introducing a plane of weakness in a direction normal to the crack face. Also, a shear transfer coefficient β_t (constant C₁ with TB,CONCR) is introduced which represents a shear strength reduction factor for those subsequent loads which induce sliding (shear) across the crack face. The stress-strain relations for a material that has cracked in one direction only become:

(1–107)

where the superscript ck signifies that the stress strain relations refer to a coordinate system parallel to principal stress directions with the x^ck axis perpendicular to the crack face. If KEYOPT(7) = 0, R^t = 0.0. If KEYOPT(7) = 1, R^t is the slope (secant modulus) as defined in the figure below. R^t works with adaptive descent and diminishes to 0.0 as the solution converges.

Figure 1.11: Strength of Cracked Condition

where:

If the crack closes, then all compressive stresses normal to the crack plane are transmitted across the crack and only a shear transfer coefficient β_c (constant C₂ with TB,CONCR) for a closed crack is introduced. Then can be expressed as

(1–108)

The stress-strain relations for concrete that has cracked in two directions are:

(1–109)

If both directions reclose,

(1–110)

The stress-strain relations for concrete that has cracked in all three directions are:

(1–111)

If all three cracks reclose, Equation 1–110 is followed. In total there are 16 possible combinations of crack arrangement and appropriate changes in stress-strain relationships incorporated in SOLID65. A note is output if 1 >β_c >β_t >0 are not true.

The transformation of to element coordinates has the form

(1–112)

where [T^ck] has a form identical to Equation 1–104 and the three columns of [A] in Equation 1–105 are now the principal direction vectors.

The open or closed status of integration point cracking is based on a strain value called the crack strain. For the case of a possible crack in the x direction, this strain is evaluated as

(1–113)

where:

The vector {ε^ck} is computed by:

(1–114)

where:

{ε'}, in turn, is defined as:

(1–115)

where:

n = substep number

{Δεn} = total strain increment (based on {Δun}, the displacement increment over the substep)

If is less than zero, the associated crack is assumed to be closed.

If is greater than or equal to zero, the associated crack is assumed to be open. When cracking first occurs at an integration point, the crack is assumed to be open for the next iteration.

If the material at an integration point fails in uniaxial, biaxial, or triaxial compression, the material is assumed to crush at that point. In SOLID65, crushing is defined as the complete deterioration of the structural integrity of the material (e.g. material spalling). Under conditions where crushing has occurred, material strength is assumed to have degraded to an extent such that the contribution to the stiffness of an element at the integration point in question can be ignored.