ConcreteCM
Developed and Implemented by:
Kristijan Kolozvari , California State University, Fullerton
Kutay Orakcal<span style=“color:black”>, Bogazici University, Istanbul, Turkey
John Wallace<span style=“color:black”>, Univeristy of California, Los Angeles
This command is used to construct a uniaxialMaterial ConcreteCM (Kolozvari et al., 2015), which is a uniaxial hysteretic constitutive model for concrete developed by Chang and Mander (1994). This model is a refined, rule-based, generalized, and non-dimensional constitutive model that allows calibration of the monotonic and hysteretic material modeling parameters, and can simulate the hysteretic behavior of confined and unconfined, ordinary and high-strength concrete, in both cyclic compression and tension (Figure 1). The model addresses important behavioral features, such as continuous hysteretic behavior under cyclic compression and tension, progressive stiffness degradation associated with smooth unloading and reloading curves at increasing strain values, and gradual crack closure effects. Details of the model are available in the report by Chang and Mander (1994).

The Chang and Mander (1994) model successfully generates continuous hysteretic stress-strain relationships with slope continuity for confined and unconfined concrete in both compression and tension. The compression envelope curve of the model is defined by the initial tangent slope, (E<sub class=“subscript”>c</sub>), the peak coordinate (<math></math>‘<sub class=“subscript”>c</sub>, f’<sub class=“subscript”>c</sub>), a parameter (r<sub class=“subscript”>c</sub>) from Tsai’s (1988) equation defining the shape of the envelope curve, and a parameter (<math></math><sup class=“superscript”>-</sup><sub class=“subscript”>cr</sub>) to define normalized (with respect to \(\epsilon\)‘<sub class=“subscript”>c</sub>) strain where the envelope curve starts following a straight line, until zero compressive stress is reached at the spalling strain, \(\epsilon\)<sub class=“subscript”>sp</sub>. These parameters can be controlled based on specific experimental results for a refined calibration of the compression envelope (Figure 2). Chang and Mander (1994) proposed empirical relationships for parameters E<sub class=“subscript”>c</sub>, <math></math>’<sub class=“subscript”>c</sub>, and r<sub class=“subscript”>c</sub> for unconfined concrete with compressive strength f’<sub class=“subscript”>c</sub>, based on review of previous research. Parameters f’<sub class=“subscript”>c</sub>, <math></math>’<sub class=“subscript”>c</sub>, E<sub class=“subscript”>c</sub>, r<sub class=“subscript”>c</sub>, and <math></math><sup class=“superscript”>-</sup><sub class=“subscript”>cr</sub> can also be calibrated to represent the stress-strain behavior of confined concrete in compression, to follow the constitutive relationships for confined concrete proposed by Mander et al (1988) or similar.

The shape of the tension envelope curve in the model is the same as that of the compression envelope; however, the tension envelope curve is shifted to a new origin that is based on the unloading strain from the compression envelope (Figure 2). As well, the strain ductility experienced previously on the compression envelope is also reflected on the tension envelope. The parameters associated with the tension envelope curve include the tensile strength of concrete (f<sub class=“subscript”>t</sub>), the monotonic strain at tensile strength ( \(\epsilon_\textrm{t}\) ), a parameter (r<sub class=“subscript”>t</sub>) from Tsai’s (1988) equation defining the shape of the tension envelope curve, and a parameter (<math></math><sup class=“superscript”>+</sup><sub class=“subscript”>cr</sub>) to define normalized (with respect to \(\epsilon\)<sub class=“subscript”>t</sub>) strain where the tension envelope curve starts following a straight line, until zero tensile stress is reached at a strain of \(\epsilon\)<sub class=“subscript”>crk</sub>. These parameters can also be controlled and calibrated based on specific experimental results or empirical relations proposed by other researchers (e.g., Belarbi and Hsu, 1994) to model the behavior of concrete in tension and the tension stiffening phenomenon. Concrete experiencing tension stiffening can be considered not to crack completely; that is, a large value for parameter <math></math><sup class=“superscript”>+</sup><sub class=“subscript”>cr</sub> (e.g., 10000) can be defined.
Source: /usr/local/cvs/OpenSees/SRC/material/uniaxial/
Input Format:
uniaxialMaterial ConcreteCM $mattag $fpcc $epcc $Ec $rc
$xcrn $ft $et $rt $xcrp <-GapClose $gap>
mattag
|
Unique uniaxialMaterial tag |
fpcc
|
Compressive strength (f’<sub class=“subscript”>c</sub>) |
epcc
|
Strain at compressive strength (<math></math>’<sub class=“subscript”>c</sub>) |
Ec
|
Initial tangent modulus (E<sub class=“subscript”>c</sub>) |
rc
|
Shape parameter in Tsai’s equation defined for compression (r<sub class=“subscript”>c</sub>) |
xcrn
|
Non-dimensional critical strain on compression envelope (<math></math><sup class=“superscript”>-</sup><sub class=“subscript”>cr</sub>, where the envelope curve starts following a straight line) |
ft
|
Tensile strength (f<sub class=“subscript”>t</sub>) |
et
|
Strain at tensile strength ( \(\epsilon_\textrm{t}\) ) |
rt
|
Shape parameter in Tsai’s equation defined for tension (r<sub class=“subscript”>t</sub>) |
xcrp
|
Non-dimensional critical strain on tension envelope (<math></math><sup class=“superscript”>+</sup><sub class=“subscript”>cr</sub>, where the envelope curve starts following a straight line - large value [e.g., 10000] recommended when tension stiffening is considered) |
<-GapClose $gap> |
gap = 0, less gradual gap closure (default); gap = 1, more gradual gap closure |
Example:
uniaxialMaterial ConcreteCM 1 -6.2 -0.0021 4500 7 1.035 0.30 0.00008 1.2 10000
Example of hysteretic stress-strain history generated by the model code is illustrated in Figure 3.

Discussion:
An optional input parameter gap is introduced in the ConcreteCM model implemented in OpenSees for providing the users with the opportunity to control the intensity of gap closure in the stress-strain behavior of concrete, which in-turn influences the level of pinching in the lateral load-displacement behavior of a RC wall. The original Chang and Mander (1994) model adopts a non-zero tangent stiffness at zero stress level upon unloading from the tension envelope, which is represented by gap = 1 in ConcreteCM. Using gap = 0 (default) produces less gradual gap closure, since it assumes zero tangent stiffness at zero stress level upon unloading from the tension envelope, and is suitable for most analyses. Figure 4 illustrates the effect of plastic stiffness upon unloading from tension envelope (E<sup class=“superscript”>+</sup><sub class=“subscript”>pl</sub>) on crack closure, i.e. use of more gradual (gap = 1) or less gradual (gap = 0) gap closure. The effect of parameter gap on predictions of flexural behavior of a RC wall is illustrated in Example 1 of MVLEM element.

Constitutive stress-strain concrete behavior is also implemented in OpenSees in uniaxialMaterial Cocnrete07. However, ConcreteCM incorporates sophisticated unloading/reloading rules defined originally by Chang and Mander (1994), as opposed to Concrete07 that adopts simplified hysteretic rules. Comparison between stress-strain response predicted using ConcreteCM and Concrete07 is shown in Figure 5.

References:
- Belarbi H. and Hsu T.C.C. (1994), “Constitutive Laws of Concrete in Tension and Reinforcing Bars Stiffened by Concrete”, ACI Structural Journal, V. 91, No. 4, pp. 465-474.
- Chang, G.A. and Mander, J.B. (1994), “Seismic Energy Based Fatigue Damage Analysis of Bridge Columns: Part I - Evaluation of Seismic Capacity”, NCEER Technical Report No. NCEER-94-0006, State University of New York, Buffalo.
- Kolozvari K., Orakcal K., and Wallace J. W. (2015). “Shear-Flexure Interaction Modeling of reinforced Concrete Structural Walls and Columns under Reversed Cyclic Loading”, Pacific Earthquake Engineering Research Center, University of California, Berkeley, PEER Report No. 2015/12
- Mander J.B., Priestley M.J.N., and Park R. (1988). “Theoretical Stress-Strain Model for Confined Concrete”, ASCE Journal of Structural Engineering, V. 114, No. 8, pp. 1804-1826.
- Orakcal K.(2004), “Nonlinear Modeling and Analysis of Slender Reinforced Concrete Walls”, PhD Dissertation, Department of Civil and Environmental Engineering, University of California, Los Angeles.
- Tsai W.T. (1988), “Uniaxial Compressional Stress-Strain Relation of Concrete”, ASCE Journal of Structural Engineering, V. 114, No. 9, pp. 2133-2136.
- Orakcal K.(2004), “Nonlinear Modeling and Analysis of Slender Reinforced Concrete Walls”, PhD Dissertation, Department of Civil and Environmental Engineering, University of California, Los Angeles.
- Mander J.B., Priestley M.J.N., and Park R. (1988). “Theoretical Stress-Strain Model for Confined Concrete”, ASCE Journal of Structural Engineering, V. 114, No. 8, pp. 1804-1826.
- Kolozvari K., Orakcal K., and Wallace J. W. (2015). “Shear-Flexure Interaction Modeling of reinforced Concrete Structural Walls and Columns under Reversed Cyclic Loading”, Pacific Earthquake Engineering Research Center, University of California, Berkeley, PEER Report No. 2015/12
- Chang, G.A. and Mander, J.B. (1994), “Seismic Energy Based Fatigue Damage Analysis of Bridge Columns: Part I - Evaluation of Seismic Capacity”, NCEER Technical Report No. NCEER-94-0006, State University of New York, Buffalo.