Mazars damage mechanics constitutive law for concrete
This law is only used for mechanical analysis for concrete solids undergoing large strains.
Prepro: LMAZAR.F
| Plane stress state | YES |
| Plane strain state | YES |
| Axisymmetric state | YES |
| 3D state | YES |
| Generalized plane state | NO |
| Line 1 (2I5, 60A1) | |
|---|---|
| IL | Law number |
| ITYPE | 588 |
| COMMENT | Any comment (up to 60 characters) that will be reproduced on the output listing |
| Line 1 (5I5) | |
|---|---|
| NINTV | (=1 par défaut) |
| ISOL | = 0 : use of total stresses in the constitutive law |
| ≠ 0 : use of effective stresses in the constitutive law : see Appendix 8 | |
| INOLO | Non-local approach index (0 = local approach, 1 = non-local approach) |
| ITEMOIN | Constitutive matrix index (only for 2D analysis) |
| ITEMP | Temperature index (only for 2D analysis) |
| = 0 : No temperature influence | |
| = 1 : Dotreppe's formulation (only for 2D analysis) | |
| = 2 : Eurocode 2 formulation (only for 2D analysis) | |
| = 3 : Bakker-Stabler formulation (only for 2D analysis) | |
| Line 1 (3G10.0) | |
|---|---|
| E | YOUNG's elastic modulus |
| ANU | POISSON's ratio. |
| FT | Tensile strength |
| Line 2 (5G10.0) | |
| ACOM | 1st parameter for damage due to compresssion |
| BCOM | 2nd parameter for damage due to compresssion |
| ATRA | 1st parameter for damage due to tension |
| BTRA | 2nd parameter for damage due to tension |
| BETA | Parameter for damage due to shear |
| Line 3 (4G10.0) | |
| DIV | Time integration parameter |
| RHOS | Specific mass |
| RMAX | Maximum radius (Non-local method) |
| ALC | Internal length (Non-local method) |
= 6 : for the 3-D state
= 4 : for the other cases.
The stresses are the components of CAUCHY stress tensor in global (X,Y,Z) coordinates.
For the 3-D state :
| SIG(1) | $\sigma_{xx}$ |
| SIG(2) | $\sigma_{yy}$ |
| SIG(3) | $\sigma_{zz}$ |
| SIG(4) | $\sigma_{xy}$ |
| SIG(5) | $\sigma_{xz}$ |
| SIG(6) | $\sigma_{yz}$ |
For the other cases :
| SIG(1) | $\sigma_{xx}$ |
| SIG(2) | $\sigma_{yy}$ |
| SIG(3) | $\sigma_{xy}$ |
| SIG(4) | $\sigma_{zz}$ |
= 17 : for 2D plane strain or axisymmetric analysis
= 16 : for 3D analysis
2D plane strain or axisymmetric analysis :
| Q(1) | THICK : Thickness (axisymmetric or plane stress state) |
| Q(2) | D : Damage |
| Q(3) | DC : Damage due to compression |
| Q(4) | DT : Damage due to tension |
| Q(5) | EPSTIL : Equivalent deformation |
| Q(6) | RHOS : Actualised specific mass |
| Q(7) | DXXT : Strain computed at the previous step |
| Q(8) | DYYT : Strain computed at the previous step |
| Q(9) | DZZT : Strain computed at the previous step |
| Q(10) | DXYT : Strain computed at the previous step |
| Q(11) | FLAG : Flag for increasing damage (=1- or constant damage (=0) |
| Q(12) | EPSB1 |
| Q(13) | EQVIL1 |
| Q(14) | ESPNL : Non-local equivalent strain |
| Q(15) | OMEGA : Area associated to integration point |
| Q(16) | EPSD0 damage threshold |
| Q(17) | Thermal damage G |
3D analysis :
| Q(1) | THICK : Thickness (axisymmetric or plane stress state) |
| Q(2) | D : Damage |
| Q(3) | DC : Damage due to compression |
| Q(4) | DT : Damage due to tension |
| Q(5) | EPSTIL : Equivalent deformation |
| Q(6) | RHOS : Actualised specific mass |
| Q(7) | DXXT : Strain computed at the previous step |
| Q(8) | DYYT : Strain computed at the previous step |
| Q(9) | DZZT : Strain computed at the previous step |
| Q(10) | DXYT : Strain computed at the previous step |
| Q(11) | DXZT : Strain computed at the previous step |
| Q(12) | DYZT : Strain computed at the previous step |
| Q(13) | FLAG : Flag for increasing damage (=1- or constant damage (=0) |
| Q(14) | EPSB1 |
| Q(15) | EQVIL1 |
| Q(16) | EPSD0 damage threshold |