Specific characteristics of concrete ==== .. _RefNumPara__120069_1070254933: Keyword factor THER_HYDR ---- Allows you to define the behavior associated with the hydration of concrete. The hydration of concrete is a phenomenon that is accompanied by a release of heat depending on the temperature [:ref:`R7.01.12 `]. :math:`\begin{array}{cc}\begin{array}{}\frac{d\beta }{\mathrm{dt}}+\text{div}q=Q\frac{d\xi (T)}{\mathrm{dt}}+s\\ q=-\lambda \mathrm{grad}T\end{array}\}& \text{éq 7.1-1}\\ \frac{d\xi }{\mathrm{dt}}=\mathrm{AFF}(\xi ,T)& \text{éq 7.1-2}\end{array}` Operands LAMBDA/BETA ~~~~ .. code-block:: text LAMBDA = normal Isotropic thermal conductivity as a function of temperature. .. code-block:: text BETA = beta Volume enthalpy as a function of temperature. The extensions are at least linear, the volume enthalpy can be defined as the integral of the volume heat. Operand AFFINITE ~~~~ .. code-block:: text AFFINITE = AFF Depending on the degree of hydration and the temperature. In general, we use: :math:`\text{AFF}(\xi ,T)=A(\xi )\mathrm{exp}(-\frac{{E}_{a}}{\mathrm{RT}})` with :math:`\text{QSR\_K}=\frac{{E}_{a}}{R}` the Arrhenius constant expressed in degrees Kelvin, and A determined by a calorimetric test of concrete (a function of the quantity HYDR). Operand CHAL_HYDR ~~~~ .. code-block:: text CHAL_HYDR = Q Heat released per hydration unit (assumed to be constant), this function depends on the type of concrete. .. _RefNumPara__120071_1070254933: Keyword factor SECH_GRANGER ---- Definition of the parameters characterizing the diffusion coefficient :math:`D\left(C,T\right)` used in the nonlinear drying equation proposed by Granger (cf. [:ref:`R7.01.12 `]). These characteristics are constants, while the diffusion coefficient depends on the calculation variable, i.e. the current :math:`C` water concentration, (as thermal conductivity depended on temperature). A/B operands/ QSR_K/TEMP_0_C ~~~~ These coefficients make it possible to express the diffusion coefficient in its form most commonly used in the literature and proposed by Granger: :math:`D\left(C,T\right)=a\text{.}{e}^{\left(b\text{.}C\right)}\frac{T}{{T}_{0}}{e}^{\left[-\frac{Q}{R}\left(\frac{1}{T}-\frac{1}{{T}_{0}}\right)\right]}` .. code-block:: text A= a Diffusion coefficient varying from :math:`0.5{10}^{-13}` and :math:`{2.10}^{-13}{m}^{2}/s` for concrete. .. code-block:: text B= b Coefficient of the order of :math:`0.05` for concrete. .. code-block:: text QSR_K = QsR QsR is generally :math:`4700.K`. (:math:`R` is the ideal gas constant). .. code-block:: text TEMP_0_C = T0 Reference temperature in the Arrhenius law. Reference temperature :math:`\mathit{T0}` is **in degrees Celsius**, and converted to Kelvin during resolution. .. _RefNumPara__120073_1070254933: Keyword factor SECH_MENSI ---- Definition of the parameters characterizing the diffusion coefficient used in the nonlinear drying equation proposed by Mensi (cf. [:ref:`R7.01.12 `]). These characteristics are constants, while the diffusion coefficient depends on the calculation variable, i.e. the current :math:`C` water concentration, (as thermal conductivity depended on temperature). It is a simplified formulation of the general case, constituting Mensi's law. A/B operands These coefficients make it possible to express the diffusion coefficient according to Mensi's law: :math:`\xi =\frac{{\epsilon }_{r}^{n-1}-{\epsilon }_{y}^{n-1}}{{\epsilon }_{y}^{n}-{\epsilon }_{r}^{n-1}}` :math:`D\left(C\right)=a\text{.}{e}^{\left(b\text{.}C\right)}` A= a Diffusion coefficient varying from :math:`0.5{.10}^{-13}` and :math:`{2.10}^{-13}{m}^{2}/s` for concrete. .. code-block:: text B= b Coefficient of the order of :math:`0.05` for concrete. .. _RefNumPara__120075_1070254933: Keyword factor SECH_BAZANT ---- Definition of the parameters characterizing the diffusion coefficient used in the nonlinear drying equation proposed by Bazant (confer [:ref:`R7.01.12 `]). These characteristics are constants, while the diffusion coefficient depends on the calculation variable, i.e. the current :math:`C` water concentration, (as thermal conductivity depended on temperature). This formulation constitutes Bazant's law. Operands D1/ ALPHA_BAZANT /N/ FONC_DESORP ~~~~ These coefficients make it possible to express the diffusion coefficient according to Bazant's law: :math:`D(h)={d}_{1}(\alpha +\frac{1-\alpha }{1+{(\frac{1-h}{1-0\text{.}\text{75}})}^{n}})` where :math:`h` is the degree of hydration, linked to the water concentration by the desorption curve. .. code-block:: text D1 = d1 Diffusion coefficient which is of the order of :math:`{3.10}^{-13}{m}^{2}/s` for concrete. .. code-block:: text ALPHA_BAZANT = alpha Coefficient varying from :math:`0.025` to :math:`0.1` for concrete. .. code-block:: text N = n Exponent of the order of 6 for concrete. .. code-block:: text FONC_DESORP = desorp Desorption curve, allowing to go from water concentration to hydration degree :math:`h`. **Important note:** desorp *is a function of the calculation variable,* :math:`C` * *, the water concentration, which is assimilated for the resolution at a temperature, of type* *' TEMP ' .* .. _RefNumPara__120077_1070254933: Keyword factor SECH_NAPPE ---- The diffusion coefficient, characterizing the non-linear drying equation, is expressed using a sheet, a tabulated function of the water concentration, a calculation variable, and of the temperature, an auxiliary calculation variable, given in the form of a data structure of the evol_ther type. For the resolution of drying by operator THER_NON_LINE, the water concentration is assimilated to a temperature, of the 'TEMP' type. For the consistency of the data, the parameters of the table, i.e. the calculation variable and the auxiliary variable, cannot be of the same type. A new type of variable has been added in DEFI_NAPPE, the "type of temperature calculated prior to drying", 'TSEC', which actually corresponds to a temperature. Operand FONCTION ~~~~ The diffusion coefficient is expressed using a tabulated function of the parameters :math:`C` and :math:`T\text{.}` .. code-block:: text FONCTION = fund_name Tablecloth name. .. _RefNumPara__120079_1070254933: Keyword factor PINTO_MENEGOTTO ---- Definitions of the coefficients of the cyclic elastoplasticity behavior relationship of steel reinforcements in reinforced concrete according to the Pinto-Menegotto model (cf. [:ref:`R5.03.09 `]). The initial traction curve (start of loading) is defined by: * :math:`\sigma =E\mathrm{.}\epsilon` as :math:`\sigma \le {\sigma }_{y}`; :math:`E` defined as ELAS * :math:`\sigma ={\sigma }_{y}` for :math:`\frac{{\sigma }_{y}}{E}\le \epsilon \le {\epsilon }_{h}` * :math:`\sigma \mathrm{=}{\sigma }_{u}\mathrm{-}({\sigma }_{u}\mathrm{-}{\sigma }_{y}){(\frac{{\varepsilon }_{u}\mathrm{-}\varepsilon }{{\varepsilon }_{u}\mathrm{-}{\varepsilon }_{h}})}^{4}` for :math:`{\epsilon }_{h}\le \epsilon <{\epsilon }_{u}` (:math:`\varepsilon` can't go past :math:`{\varepsilon }_{u}`) Curve :math:`s=f(e)` in the :math:`{n}^{\mathit{ième}}` cycle is defined by: :math:`{\sigma }_{L}^{\text{*}}=b{\varepsilon }_{L}^{\text{*}}+(\frac{1-b}{{(1+{({\varepsilon }_{L}^{\text{*}})}^{R})}^{1/R}}){\varepsilon }_{L}^{\text{*}}` with :math:`R={R}_{0}-\frac{{a}_{1}\xi }{{a}_{2}+\xi }` and :math:`b=\frac{{E}_{h}}{E}`, :math:`{E}_{h}`: asymptotic work hardening slope where :math:`{\epsilon }^{\text{*}}` is defined by: :math:`{\varepsilon }^{\text{*}}=\frac{\varepsilon -{\varepsilon }_{r}^{n-1}}{{\varepsilon }_{y}^{n}-{\varepsilon }_{r}^{n-1}}`. where :math:`{\sigma }^{\text{*}}` is defined by: :math:`{\sigma }^{\text{*}}=\frac{\sigma -{\sigma }_{r}^{n-1}}{{\sigma }_{y}^{n}-{\sigma }_{r}^{n-1}}`. Quantity :math:`{\epsilon }_{y}^{n}` is deduced from cycle :math:`n-1` by: :math:`\begin{array}{}{\varepsilon }_{y}^{n}={\varepsilon }_{r}^{n-1}+\frac{{\sigma }_{y}^{n}-{\sigma }_{r}^{n-1}}{E}\\ {\sigma }_{y}^{n}={\sigma }_{y}^{n-1}\text{.}\mathrm{sign}({\varepsilon }_{y}^{n-1}-{\varepsilon }_{r}^{n-1})+{\varepsilon }_{H}({\varepsilon }_{r}^{n-1}-{\varepsilon }_{y}^{n-1})\end{array}` Variable :math:`\xi` is defined by: where :math:`{\epsilon }_{r}^{n-1}` represents the deformation reached at the end of the :math:`n-1` th half-cycle and :math:`{\epsilon }_{y}^{n-1},{\epsilon }_{y}^{n}` represent the deformations at the end of linearity of the half-cycles and :math:`n`. :math:`b` represents either the value provided by the user (keyword EP_SUR_E) or, failing that: :math:`b=\frac{{E}_{H}}{E}\text{avec}{E}_{H}=\frac{{\sigma }_{u}-{\sigma }_{y}}{{\varepsilon }_{u}-\frac{{\sigma }_{y}}{E}}` In case of buckling, (if :math:`L/D>5`): * in compression we replace :math:`b` by :math:`{b}_{c}=a\left(5.0-L/D\right){e}^{\left(b{\xi }^{\text{'}}\frac{E}{{\sigma }_{y}-{\sigma }_{\mathrm{\infty }}}\right)}` * In traction, we calculate a new slope :math:`{E}_{r}=E\left({a}_{5}+\left(1.0-{a}_{5}\right){e}^{\left(-{a}_{6}\left({\epsilon }_{r}^{n-1}-{\epsilon }_{y}^{n-1}\right)\right)}\right)` with :math:`{a}_{5}=1+\frac{5-L/D}{7.5}`. :math:`{\xi }^{\text{'}}` represents the biggest "plastic excursion" during loading: :math:`{\xi }^{\text{'}}=\underset{n}{\mathit{max}}\left({\epsilon }_{r}^{n}-{\epsilon }_{y}^{n}\right)` and :math:`{\sigma }_{\infty }=4\frac{{\sigma }_{y}}{L/D}` In the case of buckling, we add to :math:`{\sigma }_{y}^{n}` the value :math:`{\sigma }_{s}^{\text{*}}={\gamma }_{s}bE\frac{b-{b}_{c}}{1-{b}_{c}}` with :math:`{\gamma }_{s}=\frac{11-L/D}{10({e}^{\frac{\mathrm{cL}}{D}}-1)}`. Operands ~~~~ .. code-block:: text SY = sign Initial elasticity limit, noted :math:`{\sigma }_{y}` in the equations. .. code-block:: text EPSI_ULTM = epsu, noted :math:`{\varepsilon }_{u}` in the equations. Ultimate deformation. SIGM_ULTM = sigmu, noted :math:`{\sigma }_{u}` in the equations. Ultimate constraint. .. code-block:: text ◊ ELAN = L/D Slender bar (>5: buckling). EPSP_HARD = epsh, noted :math:`{\varepsilon }_{h}` in the equations. Deformation corresponding to the end of the plastic bearing. .. code-block:: text ◊ EP_SUR_E = b Hardening slope to Young's modulus ratio (if no value is given, we take :math:`b=\frac{{E}_{H}}{E}`). .. code-block:: text A 1_PM = a1 Coefficient defining the traction curve of the model. .. code-block:: text A 2_PM = a2 Coefficient defining the traction curve of the model. .. code-block:: text A 6_PM = a6 Coefficient defining the tensile curve of the model in case of buckling. .. code-block:: text C_PM = c used in :math:`{\gamma }_{s}` Coefficient defining the tensile curve of the model in case of buckling. .. code-block:: text A_PM = a Coefficient defining the tensile curve of the model in case of buckling. .. code-block:: text R_PM = Coefficient :math:`{R}_{O}` (20. by default). The Young's modulus E and the thermal expansion coefficient ALPHA should be specified by the keywords ELAS or ELAS_FO. .. _RefNumPara__120081_1070254933: Tags factor BPEL_BETON, BPEL_ACIER ---- Definition of the characteristics involved in the behavior model of prestressed cables in the regulatory framework of BPEL [:ref:`R7.01.02 `]. The linear elastic characteristics of the concrete material and the steel material must be simultaneously defined under the keyword ELAS. Operands ~~~~ Behavior: BPEL_BETON Keyword factor for the definition of the characteristic parameters of the concrete material that are involved in the estimation of tension losses along prestress cables. This factor keyword can only be used in conjunction with the factor keyword ELAS. .. code-block:: text PERT_FLUA = xflu Flat rate of tension loss due to concrete creep, compared to the initial tension. :math:`\mathrm{\Delta }{F}_{\mathit{flu}}={x}_{\mathit{flu}}\mathrm{.}{F}_{0}` where :math:`{F}_{0}` refers to the initial voltage defined by DEFI_CABLE_BP. [U4.42.04] The value by default is 0: in this case, tension losses due to concrete creep are not taken into account. Attention, this value will not be affected by entering the relaxation coefficient R_J in DEFI_CABLE_BP. The value xflu must therefore take this effect into account (multiplication by :math:`r(t)=\frac{t}{t+9{r}_{m}}`, :math:`t` corresponding to the date on which we want to estimate the state of the structure and :math:`{r}_{m}` the average radius). .. code-block:: text PERT_RETR = xret Flat rate of tension loss due to concrete shrinkage, compared to the initial tension. :math:`\mathrm{\Delta }{F}_{\mathit{ret}}={x}_{\mathit{ret}}\mathrm{.}{F}_{0}` where :math:`{F}_{0}` refers to the initial voltage. The value by default is 0: in this case, tension losses due to concrete shrinkage are not taken into account. Attention, this value will not be affected by entering the relaxation coefficient :math:`{R}_{J}` in DEFI_CABLE_BP. The value xret must therefore take this effect into account (multiplication by :math:`r(t)=\frac{t}{t+9{r}_{m}}`, :math:`t` corresponding to the date on which we want to estimate the state of the structure and :math:`{r}_{m}` the average radius). Behavior: BPEL_ACIER Keyword factor for the definition of the characteristic parameters of the steel material that are involved in the estimation of tension losses along prestress cables. This factor keyword can only be used in conjunction with the factor keyword ELAS. .. code-block:: text RELAX_1000 = rh1000 Relaxation of steel at 1000 hours, expressed in%. The value by default is 0: in this case, the tension losses due to relaxation of the steel are not taken into account. .. code-block:: text MU0_RELAX = mu0 Dimensional relaxation coefficient of prestressed steel. The value by default is 0. .. code-block:: text F_PRG = fprg Guaranteed stress of the maximum load at break (following the BPEL) If tension losses due to relaxation of steel are taken into account (RELAX_1000 given by a non-zero value), it is mandatory to fill in the F_PRG operand with a non-zero value. .. code-block:: text FROT_COURB = f Coefficient of friction between the cable and the sheath, partly curved, in :math:`{\mathit{rad}}^{–1}`. The value by default is 0. .. code-block:: text FROT_LINE = phi Coefficient of friction per unit length :math:`{m}^{-1}`, on the right side. The value by default is 0. **Note:** *FROT_LINE = FROT_COURB* PERT_LIGNE. * .. _RefNumPara__120083_1070254933: Tags factor ETCC_BETON, ETCC_ACIER ---- Definition of the characteristics involved in the behavior model of prestressed cables, in the regulatory framework of ETCC [:ref:`R7.01.02 `]. The linear elastic characteristics of the concrete material and the steel material must be simultaneously defined under the keyword ELAS. Operands ~~~~ Behavior: ETCC_BETON Keyword factor to indicate in order to be able to calculate the tension in the cables according to the formulas in ETCC. No information is required. This factor keyword can only be used in conjunction with the factor keyword ELAS. Behavior: ETCC_ACIER Keyword factor for the definition of the characteristic parameters of the steel material that are involved in the estimation of tension losses along prestress cables. This factor keyword can only be used in conjunction with the factor keyword ELAS. .. code-block:: text RELAX_1000 = rh1000 Relaxation of steel at 1000 hours, expressed in%. The value by default is 0: in this case, the tension losses due to relaxation of the steel are not taken into account. .. code-block:: text F_PRG = fprg Guaranteed stress of the maximum load at break (following ETCC). If tension losses due to relaxation of steel are taken into account (RELAX_1000 given by a non-zero value), it is mandatory to fill in the F_PRG operand with a non-zero value. .. code-block:: text COEF_FROT = f Coefficient of friction between the cable and its partly curved sheath, in :math:`{\mathit{rad}}^{–1}`. The value by default is 0. .. code-block:: text PERT_LIGNE = phi Online loss coefficient in :math:`{\mathit{rad.m}}^{-1}`. The value by default is 0. **Note:** *PERT_LIGNE = FROT_LINE/FROT_COURB.* .. _RefNumPara__120085_1070254933: Keyword factor BETON_DOUBLE_DP ---- The 3D behavior model developed in *Code_Aster* is formulated in the framework of thermo-plasticity, for the description of the non-linear behavior of concrete, in tension, and in compression, taking into account the irreversible variations in the thermal and mechanical characteristics of concrete, which are particularly sensitive to high temperatures [:ref:`R7.01.03 `]. Functions may depend on the following command variables: 'TEMP', 'INST', 'HYDR', 'SECH'. BETON_DOUBLE_DP allows you to define all the characteristics associated with the law of behavior with the double criteria of Drücker Prager. In addition to these characteristics, the elasticity module, the Poisson's ratio, and the thermal expansion coefficient :math:`\alpha`, as well as the endogenous shrinkage and desiccation shrinkage coefficients, must be defined under the keyword ELAS for the real coefficients, or ELAS_FO, for the coefficients defined by functions, or sheets. All the characteristics of the model, :math:`\text{(E, nu,}\alpha ,f\text{'}c,f\text{'}t,\beta ,\text{Gc,Gt)}` of the [function] type, may depend on one or two variables including temperature, hydration and drying. When they depend on temperature, they are a function of the maximum temperature reached during the loading history :math:`\theta`, which is stored in memory for each Gauss point, as an internal variable. This makes it possible to take into account the irreversible variations in these characteristics at high temperature. F_C /F_T/COEF_BIAX operands ~~~~ .. code-block:: text F_C= f'c Uniaxial compression strength :math:`f{\text{'}}_{c}`. .. code-block:: text F_T= f't Uniaxial tensile strength :math:`f{\text{'}}_{t}`. .. code-block:: text COEF_BIAX = beta The ratio of biaxial compressive strength to uniaxial compressive strength :math:`\beta`. Operands ENER_COMP_RUPT/ENER_TRAC_RUPT/COEF_ELAS_COMP ~~~~ .. code-block:: text ENER_COMP_RUPT = GC The breaking energy in compression :math:`{G}_{c}` *,* .. code-block:: text ENER_TRAC_RUPT = Gt Breakthrough energy in traction :math:`{G}_{t}`. .. code-block:: text COEF_ELAS_COMP = phi The compressive elastic limit, given by a proportionality coefficient as a percentage of the strength at peak :math:`{f}_{c}^{\text{'}}(\theta )`, is generally of the order of 30% for standard concrete. It is important to note that this parameter is a real and not a function. Operands LONG_CARA ~~~~ This operand makes it possible to overload the characteristic length calculated automatically, for each mesh, according to its dimensions (from its surface in 2D, from its volume in 3D). The characteristic length calculated automatically makes it possible, when the fineness of the mesh varies from one calculation to another, to maintain stable results while avoiding localization phenomena. This length, calculated automatically or given by the user, leads to the value of the ultimate tensile work hardening according to the formula (for linear post-peak work hardening): :math:`{\kappa }_{u}\left(\theta \right)=\frac{2\text{.}{G}_{t}\left(\theta \right)}{{l}_{c}\text{.}{f}_{t}^{\prime }\left(\theta \right)}` In the particular case of a mesh containing adjacent cells whose dimensions are very different, the ultimate work-outs of the model BETON_DOUBLE_DP calculated on the basis of the characteristic length of the cells are therefore very different, which can cause convergence problems or lead to a state of constraints that are not very physical. (This characteristic length is calculated from the volume of the current mesh). For this reason, it is proposed to give the user the possibility to define an average length that overrides the characteristic length calculated for each mesh. The default value of*Code_Aster* is the characteristic length calculated for each mesh. Choosing an arbitrary and identical length for all cells can also cause convergence difficulties. The best solution consists in creating a mesh whose variations in the dimensions of the cells respect the direction of variation in the stress field, and to use the characteristic length calculated automatically according to the size of the cells. Overloading by LONG_CARA should be reserved for specific cases, when the user cannot freely intervene on the mesh. In the case where the user defines the characteristic length in the material, he will choose a :math:`({G}_{t},\text{LONG\_CARA})` torque such that :math:`\frac{2\text{.}{G}_{t}\left(\theta \right)}{{l}_{c}\text{.}{f}_{t}^{\prime }\left(\theta \right)}` is worth the value he wants for the ultimate work hardening in traction :math:`{\kappa }_{u}`. (The usual value of the deformation associated with the ultimate tensile work hardening of an average concrete is :math:`5.E\mathrm{-}4`). Operands ECRO_COMP_P_PIC/ECRO_TRAC_P_PIC ~~~~ The parameters used to define the compression and tension softening curve are optional, and have values by default. .. code-block:: text ECRO_COMP_P_PIC =/'LINEAIRE' /'PARABOLE' Shape of the post-peak curve in text-type compression, which can take the values' LINEAIRE 'and' PARABOLE '. The nonlinear curve is then of the parabolic type. .. code-block:: text ECRO_TRAC_P_PIC =/'LINEAIRE' /'EXPONENT' Shape of the post-peak curve in text-type tension, which can take the values' LINEAIRE 'and' EXPONENT '. The nonlinear curve is then of the exponential type. .. _RefNumPara__120087_1070254933: Tag factor BETON_GRANGER, V_BETON_GRANGER ---- Definition of material parameters for the Granger viscoelastic model, modeling the natural flow of concrete. There are 2 behavioral relationships: the first (BETON_GRANGER) does not take into account the phenomenon of aging but models the humidity effect. The second (BETON_GRANGER_V) takes into account the effects of aging and humidity (Cf. [:ref:`R7.01.01 `]). In 1D and under the action of a constant stress :math:`{\sigma }_{0}`, the creep deformation as a function of time and time of loading :math:`{t}_{c}` is written as: :math:`{\epsilon }_{\mathit{fl}}\left(t\right)=J\left(t,{t}_{c},\right)\cdot {\sigma }_{0}` :math:`(0.5{10}^{-4}<{\epsilon }_{d0}<1.5{10}^{-4})` The creep function :math:`J\left(t,{t}_{c}\right)` is equal to: :math:`J(t,{t}_{c})=k\left(a\left({t}_{c}\right)\right)\cdot \sum _{s=0}^{n}{J}_{s}\left(1-\mathit{exp}\left(\frac{t-{t}_{c}}{{\mathrm{\tau }}_{s}}\right)\right)` where :math:`a` is the age of the material, i.e. how long it took for the concrete to be put in place. It is an internal variable in the model. The :math:`k\left(a\left({t}_{c}\right)\right)` function that appears in :math:`J\left(t,{t}_{c}\right)` is used to model aging because it introduces the direct dependence on the loading time. For example, we can use curve CEB which models aging due to hydration: :math:`k\left(a\right)=\frac{{28}^{0.2}+0.1}{{a}^{0.2}+0.1}` with :math:`a` in days Without aging this function is constant and is equal to :math:`1`, the creep function depends in this case only on the time elapsed as soon as charging :math:`t-{t}_{c}` was carried out. Hygrometry is taken into account through an equivalent constraint :math:`S=h\cdot \sigma`, :math:`h` being the relative humidity of the material. It is therefore necessary to fill in the isothermal desorption curve :math:`c` which makes it possible to go from the water content :math:`C` to :math:`h` *:* :math:`h={c}^{-1}\left(C\right)`. This curve is provided to Code_Aster under the keyword ELAS_FO (FONC_DESORP, see § :ref:`4.1.8 `). **Note:** *This behavior may be associated with drying and thermohydration withdrawals defined by the operands* *K_DESSIC,* *B_ENDOGE* *and* ALPHAsous the key word* *ELAS_FO.* For law BETON_GRANGER, the parameters of the law should be entered under the keyword: BETON_GRANGER. For law BETON_GRANGER_V, you must also enter the keyword BETON_GRANGER but you must add the keyword V_BETON_GRANGER for parameters specific to the aging law. The internal variables of the law of behavior are described in [:ref:`R7.01.01 `]. Operands for clean creep ~~~~ .. code-block:: text 8 material coefficients of the creep function, homogeneous at one time. J1 = J1 ... ... J8 = J8 .. code-block:: text 8 delay times of the creep function. TAUX_1 =tau1 ... ... TAUX_8 =tau8 Operands for aging ~~~~ .. code-block:: text FONC_V = k (a) Aging function. .. _RefNumPara__120089_1070254933: Tag factor MAZARS, MAZARS_FO ---- The Mazars behavior model is a model of damaging elastic behavior used to describe the softening behavior of concrete. It distinguishes between traction and compression behavior, but uses only one scalar damage variable (confer [:ref:`R7.01.08 `]). The Mazars model implanted corresponds to the 2012 version, i.e. to the reformulation improving bi-compression and pure shear behavior. The parameters can be a function of temperature, hydration and drying. In this case you should use MAZARS_FO. In the case of this law of behavior, the parameters depend on the maximum temperature reached during the entire loading history. Functions may depend on the following command variables: 'TEMP', 'HYDR', 'SECH'. MAZARS (or MAZARS_FO) allows you to define all the characteristics associated with the Mazars behavior model. In addition to these characteristics, elastic constants should be defined under the keyword ELAS for real coefficients or ELAS_FO for temperature-dependent coefficients. Operands: EPSD0/EPSC0/EPST0 /AC/ AT/ AT/ BC/ BT/ k ~~~~ ♦/epsd0 = epsd0 Deformation damage threshold. ♦/ePSc0 = epsc0 epst0 = epst0 Threshold of damage in deformation, compression and traction. If :math:`{\mathrm{\epsilon }}_{d0}` is given :math:`{\mathrm{\epsilon }}_{t0}` and :math:`{\mathrm{\epsilon }}_{c0}` are calculated by: :math:`\begin{array}{c}{\mathrm{\epsilon }}_{t0}={\mathrm{\epsilon }}_{d0}\\ {\mathrm{\epsilon }}_{c0}=\frac{{\mathrm{\epsilon }}_{d0}}{\mathrm{\nu }\ast \sqrt{2}}\end{array}` ♦ AC = ac Coefficient allowing to fix the shape of the post-peak curve under compression. Introduce a horizontal asymptote which is the :math:`\varepsilon` axis for :math:`\mathrm{Ac}=1` and the horizontal for going through the peak for :math:`\mathrm{Ac}=0` (generally). ♦ AT = at Coefficient for fixing the shape of the post-peak curve under traction. Introduce a horizontal asymptote which is the :math:`\varepsilon` axis for :math:`\mathrm{Ac}=1` and the horizontal through the peak for :math:`\mathrm{Ac}=0` (usually :math:`0.7<\mathrm{At}<1`). ♦ bc=BC Coefficient allowing to fix the shape of the post-peak curve under compression. Depending on its value, it may correspond to a sudden drop in stress (:math:`\mathrm{BC}<{10}^{4}`) or a preliminary phase of stress increase followed by a more or less rapid decrease (generally :math:`{10}^{3}<\mathrm{Bc}<2.{10}^{3}`). ♦ BT=bt Coefficient for fixing the shape of the post-peak curve under traction. Depending on its value, it may correspond to a sudden drop in stress (:math:`\mathrm{BC}<{10}^{4}`) or a preliminary phase of stress increase followed by a more or less rapid decrease (generally :math:`{10}^{\mathrm{4 }}<\mathrm{Bt}<{10}^{5}`). ♦ k=K Parameter introducing a horizontal asymptote in pure shear. It is between 0 and 1. Recommended value :math:`\mathrm{0,7}`. Operand CHI ~~~~ ♦ CHI =chi As part of the coupling BETON_UMLV with the law of MAZARS. The chi parameter makes it possible to define the importance of the coupling: :math:`\mathrm{CHI}=0`: no coupling, :math:`\mathrm{CHI}=1`: total coupling. Total coupling leads to premature onset of concrete damage. The recommended value is in the range: :math:`[0.4;0.7]`. Operand SIGM_LIM, EPSI_LIM ~~~~ ◊ SIGM_LIM =sglim Definition of limit stress. ◊ EPSI_LIM =eplim Definition of limit deformation. The sigm_ LIM and espi_ LIM operands make it possible to define the stress and deformation limits that correspond to the service and ultimate limit states, classically used during studies in civil engineering. The sigm_ LIM and espi_ LIM operands make it possible to define the stress and deformation limits that correspond to the service and ultimate limit states, classically used during studies in civil engineering. **Note:** These terminals are mandatory when using the mazars_ UNIL behavior (confer [:ref:`R7.01.08 `] Damage model from MAZARS, [:external:ref:`U4.42.07 `] DEFI_MATER_GC). *In other cases they are not taken into account*. .. _RefNumPara__120094_1070254933: Keyword BETON_UMLV ---- The creep law UMLV assumes a total decoupling between spherical and deviatory components: the deformations induced by spherical stresses are purely spherical and the deformations induced by deviatory stresses are purely deviatory [:ref:`R7.01.06 `]. Moreover, the natural creep deformation is assumed to be proportional to the internal relative humidity: Spherical part: :math:`{\epsilon }^{s}=h\cdot f\left({\sigma }^{s}\right)` and, deviatory part: :math:`\underline{\underline{{\epsilon }^{d}}}=h\cdot f\left(\stackrel{~}{\underline{\underline{\sigma }}}\right)` where :math:`h` refers to internal relative humidity. Behavioral model BETON_UMLV is a non-aging viscoelastic model developed in partnership with the University of Marne-la-Vallée to describe the natural creep of concrete. It is particularly suitable for multiaxial configurations by not assuming the value of the creep Poisson's ratio. Spherical stresses are at the origin of the migration of water absorbed at the interfaces between hydrates at the level of macro-porosity and absorbed within the micro-porosity in the capillary porosity. The diffusion of inter-lamellar water from the hydrate pores to the capillary porosity takes place irreversibly. The total spherical creep deformation is therefore written as the sum of a reversible part and an irreversible part: :math:`{\varepsilon }^{\mathrm{fs}}=\underset{\begin{array}{}\mathrm{partie}\\ \mathrm{réversible}\end{array}}{\underset{\underbrace{}}{{\varepsilon }_{r}^{\mathrm{fs}}}}+\underset{\begin{array}{}\mathrm{partie}\\ \mathrm{irréversible}\end{array}}{\underset{\underbrace{}}{{\varepsilon }_{i}^{\mathrm{fs}}}}` The process of spherical creep deformation is governed by the following system of coupled equations: :math:`\{\begin{array}{c}{\dot{\varepsilon }}^{\mathrm{fs}}=\frac{1}{{\eta }_{r}^{s}}\cdot \left[h\cdot {\sigma }^{s}-{k}_{r}^{s}\cdot {\varepsilon }_{r}^{\mathrm{fs}}\right]-{\dot{\varepsilon }}_{i}^{\mathrm{fs}}\\ {\dot{\varepsilon }}_{i}^{\mathrm{fs}}=\frac{1}{{\eta }_{i}^{s}}{\langle \left[{k}_{r}^{s}\cdot {\varepsilon }^{\mathrm{fs}}-({k}_{r}^{s}+{k}_{i}^{s})\cdot {\varepsilon }_{i}^{\mathrm{fs}}\right]-\left[h{\sigma }^{s}-{k}_{r}^{s}\cdot {\varepsilon }_{r}^{\mathrm{fs}}\right]\rangle }^{+}\end{array}` where :math:`{k}_{r}^{s}` refers to the apparent stiffness associated with the skeleton formed by hydrate blocks at the mesoscopic scale; :math:`{\eta }_{r}^{s}` the apparent viscosity associated with the diffusion mechanism within the capillary porosity; :math:`{k}_{i}^{s}` refers to the apparent stiffness intrinsically associated with hydrates at the microscopic scale and the :math:`{\eta }_{i}^{s}` apparent viscosity associated with the mechanism of interfoliar diffusion. (The square brackets:math: `{⟨⟩}} ^ {\ text {+}}}` refer to the Mac Cauley operator: :math: `{⟨x⟩}} ^ {\ text {+}}} =\ frac {1} {1} {2}\ left (x+|x|\ right)} =\ frac {1} {1} {2}\ left (x+|x|\ right)`) Deviatory stresses are at the origin of a sliding mechanism (or mechanism of virtual dislocation) of sheets of CSH in nano-porosity. Under deviatoric stress, creep takes place at constant volume. Moreover, creep law UMLV assumes the isotropy of deviatoric creep. Phenomenologically, the sliding mechanism includes a reversible viscoelastic contribution of water strongly adsorbed to the sheets of CSH and an irreversible viscous contribution of free water: :math:`\underset{\begin{array}{c}\mathit{déformation}\\ \mathit{déviatorique}\\ \mathit{totale}\end{array}}{\underset{\underbrace{}}{{\underline{\underline{\varepsilon }}}^{\mathit{fd}}}}\mathrm{=}\underset{\begin{array}{c}\mathit{contribution}\\ \mathit{eau}\\ \mathit{absorbée}\end{array}}{\underset{\underbrace{}}{{\underline{\underline{\varepsilon }}}_{\text{}r}^{\mathit{fd}}}}+\underset{\begin{array}{c}\mathit{contribution}\\ \mathit{eau}\\ \mathit{libre}\end{array}}{\underset{\underbrace{}}{{\underline{\underline{\varepsilon }}}_{\text{}i}^{\mathit{fd}}}}` The*th* main component of total deviatoric deformation is governed by the following system of equations: :math:`{\tilde{\dot{\sigma }}}^{j}(1+\frac{{\eta }_{r}^{d}}{{\eta }_{i}^{d}})+\frac{{k}_{r}^{d}}{{\eta }_{i}^{d}}{\tilde{\sigma }}^{j}\mathrm{=}{\eta }_{r}^{d}{\ddot{\varepsilon }}^{d,j}+{k}_{r}^{d}{\dot{\varepsilon }}^{d,j}` where :math:`{k}_{r}^{d}` refers to the stiffness associated with the capacity of the absorbed water to transmit charges (*load bearing water*); :math:`{\eta }_{r}^{d}` the viscosity associated with the water adsorbed by the hydrate sheets and :math:`{\eta }_{i}^{d}` refers to the viscosity associated with free water. Operand ~~~~ .. code-block:: text K_RS = K_RS :math:`{k}_{r}^{s}` apparent stiffness associated with the skeleton formed by hydrate blocks at the mesoscopic scale. .. code-block:: text K_IS = K_IS .. code-block:: text :math:`{k}_{i}^{s}` apparent stiffness intrinsically associated with hydrates at the microscopic scale. .. code-block:: text K_RD = K_RD .. code-block:: text :math:`{k}_{r}^{d}` stiffness associated with the capacity of adsorbed water to transmit loads (load bearing water). .. code-block:: text ETA_RS = ETA_RS .. code-block:: text :math:`{\eta }_{r}^{s}` apparent viscosity associated with the diffusion mechanism within the capillary porosity. .. code-block:: text ETA_IS = ETA_IS .. code-block:: text :math:`{\eta }_{i}^{s}` apparent viscosity associated with the interlamellar diffusion mechanism. .. code-block:: text ETA_RD = ETA_RD .. code-block:: text :math:`{\eta }_{r}^{d}` viscosity associated with water absorbed by hydrate sheets. ETA_FD = ETA_FD makes it possible to take into account desiccation creep according to Bazant's law. .. code-block:: text Note: *The desorption curve giving the hygrometry* :math:`h` *as a function of the water concentration* :math:`C` *must be entered under the keyword* *ELAS_FO.* .. _RefNumPara__120096_1070254933: Keyword factor BETON_ECRO_LINE ---- Definition of a linear work-hardening curve taking into account confinement in the specific case of concrete. In order to improve compression behavior, a reversibility threshold is defined ([:ref:`R7.01.04 `] model ENDO_ISOT_BETON). Operands ~~~~ .. code-block:: text D_SIGM_EPSI = dsde (AND) Slope of the traction curve. .. code-block:: text SYT = sigt Maximum stress in simple traction. .. code-block:: text SYC = SIGC Maximum stress in simple compression (it does not exist for a Poisson's ratio :math:`\nu =0`, in this case :math:`\mathrm{SYC}` is not specified) Young's module :math:`E` should be specified by the keywords ELAS or ELAS_FO. .. _RefNumPara__120334_1070254933: Keyword factor ENDO_ORTH_BETON ---- Definition of the parameters of the law of behavior ENDO_ORTH_BETON, allowing to describe the anisotropy induced by concrete damage, as well as the unilateral effects [:ref:`R7.01.09 `]. Refer to documents [:ref:`R7.01.09 `] and [:external:ref:`V6.04.176 `] for the precise meaning of the parameters and the identification procedure. Operand ALPHA ~~~~ Coupling constant between the evolution of tensile damage and that of compression damage. It should be taken between :math:`0` and :math:`1`, rather close to :math:`1`. The value by default is :math:`0.9`. Operands K0/K1/K2 ~~~~ K0 = k0 A constant part of the threshold function. Allows you to calibrate the height of the peak under traction. K1 = k1 Parameter of the threshold function used to increase the compression threshold. K2 = k2 Parameter for controlling the shape of the rupture envelope for biaxial tests. The value by default is :math:`{7.10}^{-4}`. Operands ECROB/ECROD ~~~~ .. code-block:: text ECROB = ecrob Term for blocked energy (equivalent to work hardening energy) relating to the evolution of tensile damage. It allows you to control the shape of the peak under traction. .. code-block:: text ECROD = record Term for blocked energy (equivalent to work hardening energy) relating to the evolution of compression damage. It makes it possible to control the shape of the peak under compression. The Young's modulus :math:`E` and the Poisson's ratio :math:`\nu` are to be specified by the keywords ELAS or ELAS_FO. .. _RefNumPara__120336_1070254933: Keywords factor ENDO_SCALAIRE/ENDO_SCALAIRE_FO ---- Definition of the parameters of the law of behavior ENDO_SCALAIRE [:ref:`R5.03.25 `], which describes the brittle elastic rupture of a homogeneous isotropic material. This law is only available for damage gradient modeling GRAD_VARI. K, P, M operand ~~~~ These are the internal parameters of the model that define work hardening, see [:ref:`R5.03.25 `]: :math:`k` refers to an energy density :math:`\mathrm{Pa}`, :math:`k` and :math:`m` are dimensionless parameters. :math:`k` and :math:`m` can be adjusted using the non-local scale :math:`D` (approximately the half-width of the location bandwidth) and the following macroscopic parameters: :math:`E` the Young's modulus, :math:`{G}_{f}` the cracking energy and :math:`{f}_{t}` the value of the peak stress in simple tension. The readjustment relationships are then written as: :math:`k=\frac{3{G}_{f}}{4D};\phantom{\rule{4em}{0ex}}m=\frac{3E{G}_{f}}{2{f}_{t}^{2}D};\phantom{\rule{4em}{0ex}}c=\frac{3}{8}D{G}_{f}` where :math:`c` is the parameter specified by NON_LOCAL = _F (C_GRAD_VARI = c), which also depends on the macroscopic response. As for the parameter :math:`p`, greater than 1, it controls the curvature of the post-peak response. Operands C_COMP, C_VOLU ~~~~ These are the internal, dimensionless parameters of the model that define the shape of the load surface (with one exception), see [:ref:`R5.03.25 `]. The default values make it possible to find the energy model (symmetric) for which the load surface corresponds to a level line of the elastic energy density (ellipsoid of rotation around the :math:`(\mathrm{1,1}\mathrm{,1})` axis which is centered at the beginning of coordinates). In the more general case, the non-centered ellipsoidal load surface (always the :math:`\mathrm{1,1}\mathrm{,1}` axis) can be defined by three parameters that are more accessible to measurement: :math:`{f}_{t}` the value of the stress at the peak in simple tension, :math:`{f}_{c}` the value of the stress at the peak in simple compression, and :math:`\tau` the value of the stress at the peak in pure shear. The adjustment relationships are as follows: :math:`{c}_{\mathit{comp}}=\frac{1+\nu }{1-2\nu }\frac{({f}_{c}-{f}_{t})\tau \sqrt{3}}{2{f}_{t}{f}_{c}};\phantom{\rule{6em}{0ex}}{c}_{\mathit{volu}}=\frac{2(1+\nu )}{1-2\nu }\left[{\left(\frac{({f}_{c}+{f}_{t})\tau \sqrt{3}}{2{f}_{t}{f}_{c}}\right)}^{2}-1\right]` Operands COEF_RIGI_MINI ~~~~ .. code-block:: text COEF_RIGI_MINI This is the parameter for regularizing the matrix tangent to the break, to avoid zero pivots if the cracking were to cut the part into several pieces not maintained by the boundary conditions. It does not depend on control variables. The Young's modulus :math:`E` and the Poisson's ratio :math:`\nu` are to be specified by the keywords ELAS or ELAS_FO. The non-location parameter is entered under the keyword C_GRAD_VARI behind the keyword factor NON_LOCAL. It is linked to the macroscopic parameters by: .. _RefNumPara__120346_1070254933: Keyword factor ENDO_FISS_EXP/ENDO_FISS_EXP_FO ---- Definition of the parameters of the law of behavior ENDO_FISS_EXP [:ref:`R5.03.27 `], which describes the brittle elastic rupture of a homogeneous isotropic material. This law is only available for damage gradient modeling GRAD_VARI. Operand K, M, P, Q ~~~~ These are the internal parameters of the model that define work hardening, see [:ref:`R5.03.27 `]. Their identification is supported by the DEFI_MATER_GC [:external:ref:`U4.42.07 `] command, using quantities that are experimentally accessible. Operands TAU, SIG0, BETA ~~~~ These are the internal parameters of the model that define the shape of the load surface (with one exception), see [:ref:`R5.03.27 `]. It is based on the Von Mises stress and the exponential stress tensor and compares well to experimental results on concrete under biaxial loading. The parameter BETA is more numerical in nature and only has the advantage of making the elasticity domain limited, including for hydrostatic compressions; the default value fulfills this function well, without affecting the shape of the domain in the areas of interest. Again, the DEFI_MATER_GC [:external:ref:`U4.42.07 `] command makes it possible to identify these parameters from quantities that are experimentally accessible (limits in tension and compression). Operand REST_RIGIDITE ~~~~ Stiffness restoration is active for compression deformation directions. To avoid an abrupt change in regime during the transition from traction to compression, a function S regulates the jump in stiffness, see [:ref:`R5.03.27 `]. The parameter REST_RIGIDITE, which is positive, controls this regularization; it corresponds to the coefficient gamma of the function S'. A value of 0 leads to not restoring stiffness (i.e. the model is without restoring stiffness) while a very large value amounts to virtually eliminating regularization. The DEFI_MATER_GC [:external:ref:`U4.42.07 `] command makes it easy to quantify this parameter by stipulating what proportion of the stiffness is restored for a deformation level corresponding to the initial compression threshold. Operands COEF_RIGI_MINI ~~~~ This is the parameter for regularizing the matrix tangent to the break, to avoid zero pivots if the cracking were to cut the part into several pieces not maintained by the boundary conditions. It does not depend on control variables. The Young's modulus :math:`E` and the Poisson's ratio :math:`\nu` are to be specified by the keywords ELAS or ELAS_FO. The non-locality parameter is entered under the keyword C_GRAD_VARI behind the keyword factor NON_LOCAL; the command DEFI_MATER_GC [:external:ref:`U4.42.07 `] makes it possible to identify it from quantities that are experimentally accessible. .. _RefNumPara__120348_1070254933: Keyword factor GLRC_DM ---- This keyword factor allows you to define the parameters of the GLRC_DM law of behavior. It is a model of global damage to a reinforced concrete slab formulated in terms of generalized deformation/stress relationships (membrane extension, membrane flexure and force, bending moment), cf. [:ref:`R7.01,32 `]. Operands ~~~~ .. code-block:: text NYT = Not Membrane force of the damage threshold under simple traction of a reinforced concrete slab (unit of force per length). .. code-block:: text NYC = Nc Membrane force of the "damage" threshold (end of linearity of the compression curve) in simple compression of a reinforced concrete slab (unit of force per length). .. code-block:: text MYF = Mf Flexing moment of the damage threshold in simple bending of a reinforced concrete slab (unit of force). .. code-block:: text GAMMA_T = GMT Relative damaging slope compared to the elastic slope under simple traction (:math:`0<{\gamma }_{\mathrm{MT}}<1`). .. code-block:: text GAMMA_C = GMC Relative damaging slope compared to the elastic slope in simple compression (:math:`0<{\gamma }_{\mathit{MC}}<1`). GAMMA_F = GMF Relative damaging slope compared to the elastic slope in simple bending (:math:`0<{\gamma }_{F}<1`). .. code-block:: text ALPHA_C = Alfc Parameter for modulating the compression damage function to introduce a decoupling of the traction and compression thresholds and inducing a curvature of the compression curve. The membrane damage function is written as: :math:`{\xi }_{m}(x,{d}_{\mathrm{1,}}{d}_{2})=\frac{1}{2}((\frac{1+{\gamma }_{\mathrm{mt}}{d}_{1}}{1+{d}_{1}}+\frac{1+{\gamma }_{\mathrm{mt}}{d}_{2}}{1+{d}_{2}})H(x)+(\frac{{\alpha }_{c}+{\gamma }_{\mathrm{mc}}{d}_{1}}{{\alpha }_{c}+{d}_{1}}+\frac{{\alpha }_{c}+{\gamma }_{\mathrm{mc}}{d}_{2}}{{\alpha }_{c}+{d}_{2}})H(-x))` We can refer to the reference documentation [:ref:`R7.01.32 `] section §3.2.4 where a summary of the identification of model parameters is presented. .. _RefNumPara__120350_1070254933: Keyword factor DHRC ---- This factor keyword makes it possible to define the parameters of the DHRC law of behavior. It is a global damage model of a reinforced concrete slab formulated using a homogenization method, in terms of generalized deformation/stress relationships (membrane extension, membrane flexure and force, flexing moment) and including internal state variables of damage and sliding at the steel-concrete interface, see [:external:ref:`R7.01.36 `]. The :math:`258` parameters of the law to be identified, by homogenization and by the method of least squares on different damage values, correspond to: * to the parameters controlling the components of the tensors of damageable elastic stiffness :math:`A`, of coupling of generalized deformations and slips :math:`B` and of stored energy in sliding :math:`C`, for which no analytical expression is available; * to macroscopic threshold parameters which are linked to microscopic threshold parameters. Operands ~~~~ .. code-block:: text NYD = nyd List of the two damage thresholds :math:`{G}^{\zeta ,\mathit{crit}}` in simple traction of reinforced concrete .. code-block:: text SCRIT = script List of the four sliding thresholds :math:`{\Sigma }_{\alpha }^{\zeta ,\mathit{crit}}` equivalent steel-concrete .. code-block:: text AA_C = Alpha_AC Parameters (:math:`42`) :math:`{\alpha }^{\mathrm{Ac}}` of the dependencies in component damage variables (:math:`21` supra-diagonal terms) of the tensor of order :math:`4` symmetric :math:`A` membrane-flexure of the plate, in the compression domain, in the reinforcement coordinate system :math:`(x,y)`, in Voigt notations, in Voigt notations, identified by homogenization and by the method of least squares on various values of :math:`{D}_{\rho }`, in the upper zone (:math:`1`) then in the lower zone (:math:`2`): :math:`{A}_{\beta \delta \tau \upsilon }^{\rho }({D}_{\rho })={A}_{\beta \delta \tau \upsilon }^{0}\frac{{\alpha }_{\beta \delta \tau \upsilon }^{\mathrm{Ac}\rho }+{\gamma }_{\beta \delta \tau \upsilon }^{\mathrm{Ac}\rho }{D}_{\rho }}{{\alpha }_{\beta \delta \tau \upsilon }^{\mathrm{Ac}\rho }+{D}_{\rho }}` **Note:** *In current practice, due to the isotropy of concrete and the orientation of the steels along the axis* :math:`(x,y)` *, we will have:* AAC131 = 1., AAC161 = 1., AAC231 = 1., = 1., = 1., = 1., = 1. AAC261 AAC341 AAC351 , AAC461 = 1., AAC561 = 1.; AAC132 = 1.; = 1.; = 1.; = 1.; = 1. AAC162 AAC232 AAC262 , AAC342 = 1., AAC352 = 1., AAC462 = 1., AAC562 = 1.. .. code-block:: text AA_T = Alpha_AT Parameters (:math:`42`) :math:`{\alpha }^{\mathrm{At}}` of the dependencies in component damage variables (:math:`21` supra-diagonal terms) of the tensor of order :math:`4` symmetric :math:`A` membrane-flexure of the plate, in the tensile domain, in the frame coordinate system :math:`(x,y)`, in Voigt notations, in Voigt notations, identified by homogenization and by the method of least squares on various values of :math:`{D}_{\rho }`, in the upper zone (:math:`1`) then in the lower zone (:math:`2`): :math:`{A}_{\beta \delta \tau \upsilon }^{\rho }({D}_{\rho })={A}_{\beta \delta \tau \upsilon }^{0}\frac{{\alpha }_{\beta \delta \tau \upsilon }^{\mathrm{At}\rho }+{\gamma }_{\beta \delta \tau \upsilon }^{\mathrm{At}\rho }{D}_{\rho }}{{\alpha }_{\beta \delta \tau \upsilon }^{\mathrm{At}\rho }+{D}_{\rho }}` **Note:** *In current practice, because of the isotropy of concrete and the orientation of the steels along the axis* :math:`(x,y)` *, we will take:* AAT131 = 1., AAT161 = 1., AAT231 = 1., = 1., = 1., = 1., = 1. AAT261 AAT341 AAT351 , AAT461 = 1., AAT561 = 1.; AAT132 = 1.; = 1.; = 1.; = 1.; = 1. AAT162 AAT232 AAT262 , AAT342 = 1., AAT352 = 1., AAT462 = 1., AAT562 = 1.. .. code-block:: text GA_C = gamma_AC Parameters (:math:`42`) :math:`{\gamma }^{\mathrm{Ac}}` of the dependencies in component damage variables (:math:`21` supra-diagonal terms) of the tensor of order :math:`4` symmetric :math:`A` membrane-flexure of the plate, in the compression domain, in the reinforcement coordinate system :math:`(x,y)`, in Voigt notations, in Voigt notations, identified by homogenization and by the method of least squares on various values of :math:`{D}_{\rho }`, in the upper zone (:math:`1`) then in the lower zone (:math:`2`). **Note:** *In current practice, because of the isotropy of concrete and the orientation of the steels along the axis* :math:`(x,y)` *, we will take:* GAC131 = 1., GAC161 = 1., GAC231 = 1., = 1., = 1., = 1., = 1. GAC261 GAC341 GAC351 , GAC461 = 1., GAC561 = 1.; GAC132 = 1.; = 1.; = 1.; = 1.; = 1. GAC162 GAC232 GAC262 , GAC342 = 1., GAC352 = 1., GAC462 = 1., GAC562 = 1.. .. code-block:: text GA_T = gamma_AT Parameters (:math:`42`) :math:`{\gamma }^{\mathrm{At}}` of the dependencies in damage variables of the components (:math:`21` supra-diagonal terms) of the tensor of order :math:`4` symmetric :math:`A` membrane-flexure of the plate, in the tensile domain, in the frame coordinate system :math:`(x,y)`, for sliding in the upper grid (:math:`1`) or lower ( :math:`2`), in Voigt notations, identified by homogenization and by the least squares method on various values of :math:`{D}_{\rho }`, in the upper zone (:math:`1`) then the lower zone (:math:`2`). **Note:** *In current practice, because of the isotropy of concrete and the orientation of the steels along the axis* :math:`(x,y)` *, we will take:* GAT131 = 1., GAT161 = 1., GAT231 = 1., = 1., = 1., = 1., = 1. GAT261 GAT341 GAT351 , GAT461 = 1., GAT561 = 1.; GAT132 = 1.; = 1.; = 1.; = 1.; = 1. GAT162 GAT232 GAT262 , GAT342 = 1., GAT352 = 1., GAT462 = 1., GAT562 = 1.. .. code-block:: text AB = alpha_B Parameters (:math:`24`) :math:`{\alpha }^{B}` of the dependencies in component damage variables (:math:`24` supra-diagonal terms) of the symmetric 3rd-order tensor :math:`B` of membrane-flexion-sliding coupling of the plate, in the frame coordinate system :math:`(x,y)`, for sliding in the upper grid (:math:`1`) then lower grid () then lower (:math:`2`) grid, in Voigt notations , identified by homogenization and by the least squares method on various values of :math:`{D}_{\rho }`: with: :math:`{B}_{\beta \delta \zeta }^{m\rho \pi }({D}_{\rho })=\frac{{\gamma }_{\beta \delta \zeta }^{Bm\rho \pi }{D}_{\rho }}{{\alpha }_{\beta \delta \zeta }^{Bm\rho \pi }+{D}_{\rho }}`; :math:`{B}_{\beta \delta \zeta }^{f\rho \pi }({D}_{\rho })=\frac{{\gamma }_{\beta \delta \zeta }^{Bf\rho \pi }{D}_{\rho }}{{\alpha }_{\beta \delta \zeta }^{Bf\rho \pi }+{D}_{\rho }}` **Note:** *In current practice, because of the isotropy of concrete and the orientation of the steels along the axis* :math:`(x,y)` *, we will take:* AB311 = 1., AB321 = 1., AB611 = 1., AB621 = 1. , AB312 = 1., AB322 = 1., AB612 = 1., AB622 = 1. .. code-block:: text GB = gamma_B Parameters (24) :math:`{\gamma }^{B}` of the dependencies in component damage variables (:math:`24` supra-diagonal terms) of the tensor of order :math:`3` symmetric :math:`B` of membrane-flexion-sliding coupling of the plate, in the frame coordinate system :math:`(x,y)`, for sliding in the upper grid (:math:`1`) then lower () then lower (:math:`2`) grid sliding, in Voigt notations , identified by homogenization and by the least squares method on various values of :math:`{D}_{\rho }`. **Note:** *In current practice, because of the isotropy of concrete and the orientation of the steels along the axis* :math:`(x,y)` *, we will take:* GB311 = 0, GB321 = 0, GB611 = 0, = 0, GB621 = 0, GB312 = 0, GB322 = 0, GB612 = 0, GB622 = 0. .. code-block:: text C0 = C0 Components (:math:`6` non-zero supra-diagonal terms) of the :math:`2` order symmetric tensor of the free energy of free sliding steel-concrete :math:`{C}^{0}` of the plate before damage, according to the directions of the slides considered, in the frame coordinate system :math:`(x,y)`, in the frame coordinate system, in an upper grid (:math:`1`) or lower grid () or lower (:math:`2`) grid, in Voigt notations, identified by homogenization: :math:`(\begin{array}{cccc}{C}_{\mathrm{xx}}^{01}& {C}_{\mathrm{yx}}^{01}& 0& 0\\ \text{}& {C}_{\mathrm{yy}}^{01}& 0& 0\\ \text{}& \text{}& {C}_{\mathrm{xx}}^{02}& {C}_{\mathrm{yx}}^{02}\\ \text{}& \text{}& \text{}& {C}_{\mathrm{yy}}^{02}\end{array})` **Note:** *In current practice, because of the isotropy of concrete and the orientation of the steels along the axis* :math:`(x,y)` *, we will have:* C0211= 0, C0212= 0. .. code-block:: text AC = alpha_C Parameters (:math:`6`) :math:`{\alpha }^{C}` of the dependencies in damage variables of the components (:math:`6` supra-diagonal terms) of the symmetric tensor :math:`C` identified by homogenization and by the least squares method on different values of :math:`{D}_{\rho }`: :math:`{C}_{\beta \delta }^{\rho }({D}_{\rho })={C}_{\beta \delta }^{0\rho }\frac{{\alpha }_{\beta \delta }^{C\rho }+{\gamma }_{\beta \delta }^{C\rho }{D}_{\rho }}{{\alpha }_{\beta \delta }^{C\rho }+{D}_{\rho }}` **Note:** In current practice, because of the isotropy of concrete and the orientation of the steels along the :math:`(x,y)` axes, such as: C0211= 0, C0212= 0, we will take: AC211 = 1., AC212 = 1., = 1.. .. code-block:: text GC = gamma_C Parameters (:math:`6`) :math:`{\gamma }^{C}` of the dependencies in damage variables of the components (:math:`6` supra-diagonal terms) of the symmetric tensor :math:`C` identified by homogenization and by the least squares method on different values of :math:`{D}_{\rho }`. **Note:** In current practice, because of the isotropy of concrete and the orientation of the steels along the :math:`(x,y)` axes, such as: C0211= 0, C0212= 0, we will take: GC211 = 1., GC212 = 1., = 1.. Refer to the reference documentation [:external:ref:`R7.01.36 `] where a summary of the identification of model parameters is presented. .. _RefNumPara__120352_1070254933: Keyword factor BETON_REGLE_PR ---- This keyword is used to define the material parameters used by behavior BETON_REGLE_PR ("Parabole-Rectangle" rule). This behavior can only be used in 2D (plane stresses or plane deformations) or in shells (models DKT, COQUE_3D) (see for example the ssnop129a test). It is reduced to a one-dimensional behavior, which is written, in each of the main directions of the 2D deformation tensor: * In traction: :math:`\{\begin{array}{ccc}\sigma \text{=}E\varepsilon & \text{si}& 0<\epsilon <\frac{{\sigma }_{y}^{t}}{E}\\ \sigma \text{=}{\sigma }_{y}^{t}+{E}_{T}(\varepsilon -\frac{{\sigma }_{y}^{t}}{E})& \text{si}& \frac{{\sigma }_{y}^{t}}{E}<\varepsilon <\frac{{\sigma }_{y}^{t}}{E}(1-\frac{E}{{E}_{T}})\\ \sigma \text{=}0& \text{sinon}& \end{array}` * In compression: :math:`\{\begin{array}{ccc}\sigma ={\sigma }_{y}^{c}\left[1-{(1-\frac{\varepsilon }{{\varepsilon }_{c}})}^{n}\right]& \text{si}& \varepsilon >{\varepsilon }_{c}\\ \sigma ={\sigma }_{y}^{c}& \text{sinon}& \end{array}` Operands ~~~~ .. code-block:: text DSIGM_EPSI = And Post-peak tangent module in traction :math:`{E}_{t}` (negative). .. code-block:: text SYT = Syt Ultimate tensile stress :math:`{\sigma }_{y}^{t}`. .. code-block:: text SYC = Syc Ultimate compression stress :math:`{\sigma }_{y}^{c}`. It must be given positive. .. code-block:: text EPSC = Epsc Ultimate compression deformation :math:`{\varepsilon }_{c}`. It must be given positive. .. code-block:: text N = n Exponent of the law of work hardening in compression. .. _RefNumPara__120354_1070254933: Keyword JOINT_BA ---- These parameters define the nonlinear behavior model of the steel-concrete bond used for the fine calculation of reinforced concrete structures where the prediction of cracks and the redistribution of stresses in concrete are very important. Available for analyses under the effect of monotonic and cyclic loads, the model is written in the framework of the thermodynamic formulation of irreversible processes. It makes it possible to take into account the damage to the shear interface, in combination with the effects of crack friction, as well as irreversible deformations. The [:ref:`R7.01.21 `] document describes the corresponding details. This model should be used with 2D "joint" elements [:ref:`R3.06.09 `]. Steel reinforcements can be modelled with plane (QUAD4) or one-dimensional (BARRE) elements. **Note:** *Taking into account the effect of thermal loading is not possible for the moment.* Operands ~~~~ .. code-block:: text HPEN = HPEN Parameter of penetration between surfaces by crushing concrete. We check that :math:`\mathrm{HPEN}>0.`. GTT = GTT Link stiffness module. We check that :math:`{\text{G}}_{\mathit{beton}}\mathrm{\le }\text{GTT}\mathrm{\le }{\text{G}}_{\mathit{acier}}\mathrm{.}` .. code-block:: text GAMD0 = Gam0 Perfect adherence threshold or elastic deformation limit. We check that :math:`1.E-4<\mathrm{Gam0}<1.E-2`. .. code-block:: text AD1 = ad1 Parameter for the evolution of damage in region 1 (transition from small deformations to large landslides). .. code-block:: text We check that :math:`1.E-1<\mathrm{AD1}<1.E+1`. BD1 = bd1 Power parameter describing the evolution of the damage variable in region 1 (transition from small deformations to large slides). We check that :math:`\mathrm{BD1}<1.E-1`. GAMD2 = Gam2 Threshold of major landslides. We check that :math:`1.E-4<\mathrm{Gam2}<1.E+0`. .. code-block:: text AD2 = ad2 Parameter for the evolution of damage in region 2 (maximum bond strength and friction degradation). .. code-block:: text We check that :math:`\mathrm{AD2}<1.E-6`. .. code-block:: text BD2 = bd2 Power parameter describing the evolution of the damage variable in region 2 (maximum bond strength and friction degradation). .. code-block:: text We check that :math:`\mathrm{BD2}<1.E-1`. .. code-block:: text VIFROT = vifrot Material parameter describing the influence of crack friction. .. code-block:: text We check that :math:`\mathrm{VIFROT}<0.0E+0`. .. code-block:: text FA = alpha Material parameter related to kinematic work hardening by friction of cracks. .. code-block:: text We check that :math:`\mathrm{FA}<0.0E+0`. .. code-block:: text FC = c Parameter describing the influence of confinement on bond strength. .. code-block:: text We check that :math:`\mathrm{FC}<0.0E+0`. .. code-block:: text EPSTR0 = EPSN Threshold of elastic deformation in the normal direction before rupture. We check that :math:`1.E-4<\mathrm{EPSN}<1.E+0`. .. code-block:: text ADN = DNA Parameter of damage in the normal direction by opening the crack. .. code-block:: text We check that :math:`\mathrm{ADN}<1.E-10`. .. code-block:: text BDN = good Power parameter describing the evolution of the damage variable in the normal direction. .. code-block:: text We check that :math:`\mathrm{BDN}<1.E-1`. .. _RefNumPara__120356_1070254933: Keyword BETON_RAG ---- This model is used to estimate the long-term behavior of structures affected by the alkali-granulate reaction [R7.01.26]. It makes it possible to assess deformations and anisotropic damage (cracking) of affected structures. It includes a Rankine criterion in tension and a Drücker-Prager criterion in compression. Both criteria are associated with a law of evolution leading to softening behavior. Operands ~~~~ Operands linked to the damage model ^^^^ ♦ ENDO_MC Parameter of fragility of concrete under compression. ♦ ENDO_MT Parameter of fragility of concrete under tension. ♦ ENDO_SIGUC Equivalent stress of concrete in compression. The unit of this parameter is consistent with a constraint. ♦ ENDO_SIGUT Equivalent tensile stress of concrete. The unit of this parameter is consistent with a constraint. ♦ ENDO_DRUPRA This term is a characteristic of the compression criterion. It corresponds to the angle in **radians** of the Drucker Prager criterion. Operands related to the creep model ^^^^ ♦ FLUA_SPH_KR Deferred compressibility module for reversible spherical creep. The unit of this parameter is consistent with a constraint. ♦ FLUA_SPH_KI Deferred compressibility module for irreversible spherical creep. The unit of this parameter is consistent with a constraint. ♦ FLUA_SPH_NR Viscosity coefficient for reversible spherical creep. The unit of this parameter is homogeneous to a*time constraint. ♦ FLUA_SPH_NI Viscosity coefficient for irreversible spherical creep. The unit of this parameter is homogeneous to a*time constraint. ♦ FLUA_DEV_KR Deferred compressibility module for reversible deviatoric creep. The unit of this parameter is consistent with a constraint. ♦ FLUA_DEV_KI Deferred compressibility module for irreversible deviatoric creep. The unit of this parameter is consistent with a constraint. ♦ FLUA_DEV_NR Viscosity coefficient for reversible deviatoric creep. The unit of this parameter is homogeneous to a*time constraint. ♦ FLUA_DEV_NI Viscosity coefficient for irreversible deviatoric creep. The unit of this parameter is homogeneous to a*time constraint. Operands related to endogenous withdrawal and gel formation ^^^^ ♦ PW_A, PW_B For unsaturated media, these coefficients are those of Van Genuchten. ♦ GEL_MG Gel elasticity module. The unit of this parameter is consistent with a constraint. ♦ GEL_BG 'Biot' coefficient of the gel. ♦ GEL_VG Maximum volume of gel that can be created by the chemical reaction; it corresponds to the theoretical volume of gel created per unit volume of concrete maintained under saturated conditions for an infinite amount of time. ♦ GEL_A0 Advancement at which the initial connected porosity is filled. ♦ GEL_ALPHA0 Advancement kinetics identified in Tref. ♦ GEL_TREF Reference temperature for GEL_ALPHA0, in the chemical progress equation. The unit of this parameter is homogeneous at degrees Celsius :math:`[°C]` and is converted to degrees Kelvin :math:`[K]` during resolution (addition of :math:`273.15[°C]`). ♦ GEL_EAR It is the ratio between the activation energy of the reaction and the ideal gas constant. The unit of this parameter is homogeneous at degrees Kelvin :math:`[K]`. The universal ideal gas constant is the product of the Avogadro number and the Boltzmann constant. This product is worth :math:`R\phantom{\rule{2em}{0ex}}=\phantom{\rule{2em}{0ex}}\mathrm{8,3144621}\phantom{\rule{2em}{0ex}}[J].{[\mathit{mol}]}^{-1}.{[K]}^{-1}`. ♦ GEL_SR0 Saturation threshold at which the evolution of the chemical reaction becomes possible. ♦ RAG_EPSI0 Viscous deformation for RAG. .. _RefNumPara__120358_1070254933: Keyword BETON_BURGER ---- The creep model BETON_BURGER assumes a decomposition between spherical and deviatory components: the deformations induced by spherical stresses are purely spherical and the deformations induced by deviatory stresses are purely deviatory [:ref:`R7.01.35 `]. Moreover, the natural creep deformation is assumed to be proportional to the internal relative humidity: Spherical part: :math:`{\epsilon }^{s}=h\cdot f\left({\sigma }^{s}\right)` and, deviatory part: :math:`\underline{\underline{{\epsilon }^{d}}}=h\cdot f\left(\stackrel{~}{\underline{\underline{\sigma }}}\right)` where :math:`h` refers to internal relative humidity. Behaviour model BETON_BURGER is a model based on the BETON_UMLV [:ref:`R7.01.06 `] model for describing the natural creep of concrete. It is particularly suitable for multiaxial configurations by not assuming the value of the creep Poisson's ratio. The changes made focus on taking into account a consolidation of creep translated by a non-linear term on the long-term behavior of the model. In addition, the parts spherical and deviatory are now constructed identically, leaving the possibility of controlling the apparent Poisson's ratio of creep. The spherical and deviatory parts are described by equivalent rheological chains, the so-called Bürger chain. This model was initially built using a Kelvin Voigt stage (reversible part) coupled in series to a Maxwell body (irreversible part). The model also makes it possible to take into account the effect of temperature on creep deformations via an Arrhenius law. The elastic mechanical characteristics E and NU must be defined in parallel under the keyword ELAS. Values under the ELAS keywords will be compared to the values entered under the BETON_BURGER keyword. If they are different, a fatal error will be issued. If the ELASn keyword is not entered, code_aster will do it automatically by taking the elastic characteristics of the BETON_BURGER keyword. Operands ~~~~ .. code-block:: text YoungModulus = E :math:`E` Young's module. This operand is mandatory due to the use of Mfront .. code-block:: text FishRatio = NU :math:`\mathrm{\nu }` Poisson's ratio. This operand is mandatory due to the use of Mfront .. code-block:: text K_RS = K_RS :math:`{k}_{r}^{s}` apparent stiffness associated with the reversible spherical part of the creep deformations .. code-block:: text K_RD = K_RD .. code-block:: text :math:`{k}_{r}^{d}` apparent stiffness associated with the reversible deviatoric portion of creep deformations ETA_RS = ETA_RS .. code-block:: text :math:`{\eta }_{r}^{s}` apparent viscosity associated with reversible spherical deformations .. code-block:: text ETA_IS = ETA_IS .. code-block:: text :math:`{\mathrm{\eta }}_{i}^{s}` apparent viscosity associated with irreversible spherical deformations .. code-block:: text ETA_RD = ETA_RD .. code-block:: text :math:`{\eta }_{r}^{d}` viscosity at reversible deviatory deformations .. code-block:: text ETA_ID = ETA_ID .. code-block:: text :math:`{\eta }_{i}^{d}` viscosity with irreversible deviatory deformations .. code-block:: text KAPPA = KAPPA :math:`\kappa` term affecting the long-term viscosity (:math:`{\eta }_{i}^{s}` and :math:`{\eta }_{i}^{d}`) of the material QSR_K = EAC/R :math:`{E}_{\mathit{ac}}/R` is generally equivalent to :math:`4700.K`. (:math:`R` is the ideal gas constant). .. code-block:: text TEMP_0_C = T0 Reference temperature in the Arrhenius law. Reference temperature :math:`{T}_{0}` is **in degrees Celsius**, and converted to Kelvin during resolution. .. code-block:: text ETA_FD = ETA_FD makes it possible to take into account desiccation creep according to Bazant's law. .. code-block:: text Note: *The desorption curve giving the hygrometry* :math:`h` *as a function of the water concentration* :math:`C` *must be entered under the keyword* *ELAS_FO.* .. _RefNumPara__120366_1070254933: Keyword FLUA_PORO_BETON ---- This model is used to estimate the behavior of structures subjected to delayed creep deformations [R7.01.30]. Deferred deformations of the solid skeleton are called natural creep deformations. They can be permanent or reversible, depending on the underlying phenomena. In the model, the constitutive equations of natural creep are always used coupled with a plastic deformation model that ensures the compatibility of the stress field with the resistance criteria. In practice, the proper creep model involves two usual rheological modules: a Kelvin-Voigt stage for reversible visco-elastic creep and a plastic Maxwell stage for permanent creep. The effects of intraporeous pressure are considered through the poromechanical framework. The units of the material parameters are given in R7.01.30. Operands ~~~~ Operands related to mechanical behavior ^^^^ ♦ AND Tensile strength. ♦ RC Compressive strength. ♦ DELT Coefficient for taking into account confinement. Operands related to the creep model ^^^^ ♦ TREF Reference temperature for creep. ♦ TAUK Characteristic time of the Kelvin module. ♦ YKSY Ratio: Kelvin Stiffness/Young's Modulus ♦ TAUM Characteristic time of the Maxwell module. ♦ EKFL Characteristic deformation of the creep potential. ♦ XFLU Non-linearity coefficient. ♦ NRJM Creep activation energy. ♦ DFMX Creep damage. Operands related to shrinkage/creep in desiccation ^^^^ ♦ PORO Porosity of the material. ◊ MSHR Van Genuchten module. ◊ MVGN Van Genuchten exhibitor. ◊ BSHR Biot coefficient of capillary pressure. ◊ SFLD Characteristic desiccation stress. Operands related to thermal damage ^^^^ ♦ DT80 Damage at 80°C. ♦ TSTH Threshold temperature for thermal damage. Hydration-related operands ^^^^ ♦ HYDR Advancement of hydration. ♦ HYDS Threshold for solidification of progress. .. _RefNumPara__120371_1070254933: Keyword ENDO_PORO_BETON ---- This model is used to estimate the behavior of structures subjected to mechanical (in tension, compression or shear) or thermal damage [R7.01.30]. Concrete cracking is described using a non-linear model combining plastic criteria and damages, all anisotropic. These damages form the link between total stress and effective stress in the sense of damage theory. In traction, the Rankine criteria make it possible to reproduce localized structural macro-cracks. In compression, the Drucker-Prager criterion makes it possible to establish shear damage, in connection with the expansion resulting from the non-associated flow based on the Drucker-Prager criterion. The units of the material parameters are given in R7.01.30. Operands ~~~~ Operands related to mechanical behavior ^^^^ ♦ AND Tensile strength. ♦ EPT Deformation at peak traction. ♦ RC Compressive strength. ♦ EPC Deformation at the peak of compression. ♦ DELT Coefficient for taking into account confinement. ♦ BETA Dilatance for the non-associated plastic Drucker-Prager flow. ♦ REF Crack closure stress. ♦ EKDC Characteristic plastic deformation of Drucker-Prager damage. ♦ GFT Tensile cracking energy. ♦ GFR Energy to close tensile cracks. ♦ DIM3 Size of the element in the 3rd direction (for cases AXIS or 2D). Operands related to thermal damage ^^^^ ♦ DT80 Damage at 80°C. ♦ TSTH Threshold temperature for thermal damage. Hydration-related operands ^^^^ ♦ HYDR Advancement of hydration. ♦ HYDS Threshold for solidification of progress. .. _RefNumPara__120374_1070254933: Keyword FLUA_ENDO_PORO ---- This model is used to estimate the behavior of structures subjected to delayed creep deformations and mechanical (in tension, compression or shear) or thermal damage in a coupled manner [:ref:`R7.01.30 `]. This law takes up the laws of behavior FLUA_PORO_BETON and ENDO_PORO_BETON by adding couplings between physical phenomena. The units of the material parameters are given in R7.01.30. Operands ~~~~ Operands related to mechanical behavior ^^^^ ♦ AND Tensile strength. ♦ EPT Deformation at peak traction. ♦ RC Compressive strength. ♦ EPC Deformation at the peak of compression. ♦ DELT Coefficient for taking into account confinement. ♦ BETA Dilatance for the non-associated plastic Drucker-Prager flow. ♦ REF Crack closure stress. ♦ EKDC Characteristic plastic deformation of Drucker-Prager damage. ♦ GFT Tensile cracking energy. ♦ GFR Energy to close tensile cracks. ♦ DIM3 Size of the element in the 3rd direction (for cases AXIS or 2D). Operands related to the creep model ^^^^ ♦ TREF Reference temperature for creep. ♦ TAUK Characteristic time of the Kelvin module. ♦ YKSY Ratio: Kelvin Stiffness/Young's Modulus ♦ TAUM Characteristic time of the Maxwell module. ♦ EKFL Characteristic deformation of the creep potential. ♦ XFLU Non-linearity coefficient. ♦ NRJM Creep activation energy. ♦ DFMX Creep damage. Operands related to shrinkage/creep in desiccation ^^^^ ♦ PORO Porosity of the material. ◊ MSHR Van Genuchten module. ◊ MVGN Van Genuchten exhibitor. ◊ BSHR Biot coefficient of capillary pressure. ◊ SFLD Characteristic desiccation stress. Operands related to thermal damage ^^^^ ♦ DT80 Damage at 80°C. ♦ TSTH Threshold temperature for thermal damage. Hydration-related operands ^^^^ ♦ HYDR Advancement of hydration. ♦ HYDS Threshold for solidification of progress. .. _RefNumPara__120376_1070254933: Keyword RGI_BETON ---- This model is used to estimate the behavior of structures subjected to internal swelling reactions (RGI) whether they come from the alkali-aggregate reaction (RAG) or from the DEF (Deferred Ettringite Formation) resulting from the RSI (Internal Sulfatic Reaction) [:ref:`R7.01.30 `]. The volume of expansive product created and calculated by chemical models, allows the calculation of intraperous pressure which, combined with external loading, makes it possible to evaluate anisotropic diffuse cracking using anisotropic cracking criteria and work hardening laws. It should be noted that all the phenomena and equations of law FLUA_ENDO_BETON and therefore a fortiori of FLUA_PORO_BETON and ENDO_PORO_BETON are reused in RGI_BETON. In fact, the correct consideration of swelling phenomena requires that creep and damage be taken into account in a coupled manner. The units of the material parameters are given in R7.01.30. Operands ~~~~ Operands related to the alkali-aggregate reaction (RAG) ^^^^ ♦ VRAG Maximum volume ratio of RAS. ♦ TAUG Characteristic time of RAS. ♦ SRSG Water saturation threshold to activate RAS. ♦ TRAG Reference temperature at which the test is carried out. ♦ NRJG Activation energy for the kinetics of RAS. Operands related to the internal sulphatic reaction (RSI) ^^^^ ♦ VDEF Maximum volume ratio of DEF that initializes the quantities of aluminum and sulfate. ♦ TDEF Characteristic precipitation time of DEF. ♦ SSAD Sulfates/Aluminum Ratio. ♦ SRSD Characteristic degree of saturation for precipitation. ♦ SNAK Reference alkali concentration. ♦ NRJP Activation energy for fixation and precipitation. ♦ NRJD Activation energy for dissolution. ♦ TTDD Reference threshold temperature. ♦ TTRD Reference temperature for characteristic time TDEF. ♦ TDID Characteristic time for the dissolution of sulfoaluminates. ♦ TFID Characteristic time for the fixation of aluminates in hydro-garnets. ♦ CNAD Alkaline concentration ♦ CNAK Concentration of characteristic alkalis controlling the dissolution of primary sulfoaluminates. ♦ CNAB Alkaline concentration blocking the precipitation of DEF. ♦ EXND Exponent of the dissolution-temperature coupling. ♦ EXMD Exponent of the precipitation-fixation coupling of aluminums. ♦ TTKF Threshold temperature for fixing aluminums in hydro-garnets. ♦ NRJF Fixing energy of aluminiums in hydro-garnets. Operands related to internal swelling reactions (RGI) ^^^^ ♦ VVRG Volume of voids connected to RGI creation sites. ♦ KRGI Stress concentration factor. ♦ HRGI Ratio of the plastic work hardening modulus and the Young's modulus for the plastic criterion of RGI. ♦ EKRG Characteristic deformation for the damages of RGI. ◊ ALAT Minimal chemical progress beyond which the volume of gel is effective. ◊ HPEV Multiplying coefficient of the hardening of the tensile strength for high swelling of RGI (> 0.4%). ♦ KGEL Stiffness of the product of RGI considered. Operands related to mechanical behavior ^^^^ ♦ AND Tensile strength. ♦ EPT Deformation at peak traction. ♦ RC Compressive strength. ♦ EPC Deformation at the peak of compression. ♦ DELT Coefficient for taking into account confinement. ♦ BETA Dilatance for the non-associated plastic Drucker-Prager flow. ♦ REF Crack closure stress. ♦ EKDC Characteristic plastic deformation of Drucker-Prager damage. ♦ GFT Tensile cracking energy. ♦ GFR Energy to close tensile cracks. ♦ DIM3 Size of the element in the 3rd direction (for cases AXIS or 2D). Operands related to the creep model ^^^^ ♦ TREF Reference temperature for creep. ♦ TAUK Characteristic time of the Kelvin module. ♦ YKSY Ratio: Kelvin Stiffness/Young's Modulus ♦ TAUM Characteristic time of the Maxwell module. ♦ EKFL Characteristic deformation of the creep potential. ♦ XFLU Non-linearity coefficient. ♦ NRJM Creep activation energy. ♦ DFMX Creep damage. Operands related to shrinkage/creep in desiccation ^^^^ ♦ PORO Porosity of the material. ◊ MSHR Van Genuchten module. ◊ MVGN Van Genuchten exhibitor. ◊ BSHR Biot coefficient of capillary pressure. ◊ SFLD Characteristic desiccation stress. Operands related to thermal damage ^^^^ ♦ DT80 Damage at 80°C. ♦ TSTH Threshold temperature for thermal damage. Hydration-related operands ^^^^ ♦ HYDR Advancement of hydration. ♦ HYDS Threshold for solidification of progress. Keyword RGI_BETON_BA ---- This model is an extension of law RGI_BETON integrating reinforcements to model the behavior of reinforced concrete. It uses the same parameters as well as additional parameters for defining the reinforcements. The definitions of the parameters common to RGI_BETON are not given again. Operands ~~~~ Global operands ^^^^ ♦ NREN Number of frames (between 1 and 5) ♦ YOUM Young's modulus of concrete alone. ♦ NUM Poisson's ratio of concrete alone. Operands linked to a given frame ^^^^ The following operands exist for i from 1 to 5. Parameters with no by default values are required if NREN is greater than or equal to i. ◊ ROAi Density of the frame i. ◊ Ei Young's modulus of the frame i. ◊ SYi Limit of elasticity of the frame i. ◊ TYRi Stress at the steel/concrete interface of the reinforcement i. ◊ VRi1 Coordinate x of the direction of the frame i. ◊ VRi2 Coordinate y of the reinforcement direction i. ◊ VRi3 The z coordinate of the reinforcement direction i. ◊ D_SIGM_EPSI i Armature work hardening module i. ◊ TAUKi Reversible characteristic time for the Kelvin creep of the frame i. ◊ TAUMi Irreversible characteristic time for Maxwell creep of the frame i. ◊ EKRi Characteristic deformation for the relaxation of the frame i. ◊ SKRi Characteristic constraint for setting EKRi. (The value by default is in MPa) ◊ ATRi Reference value for the thermal activation coefficient ◊ CTMi Coupling coefficient for the impact of stresses on activation energy ◊ XFLi Exhibitor for thermal activation of relaxation ◊ PREi Initial prestress (imposed as initial stress at the first step) ◊ TTRi Reference temperature for setting TAUMi (°C) ◊ XNRi Nonlinear relaxation coefficient ◊ MUSi Threshold beyond which thermal activation becomes a function of the charge rate ◊ YKYi Reversible relaxation reduction coefficient rated at TTRi .. _RefNumPara__120360_1070254933: Keyword factor ENDO_LOCA_EXP/ENDO_LOCA_EXP_FO ---- Definition of the parameters of the law of behavior ENDO_LOCA_EXP [:ref:`R7.01.42 `], which describes the homogeneous cracking of a concrete structure via an isotropic local quasibrittle damage model. Operands KAPPA and P ~~~~ These are the internal parameters of the model that define work hardening, see [:ref:`R7.01.42 `]. Their identification is supported by the DEFI_MATER_GC [:external:ref:`U4.42.07 `] command, using quantities that are more experimentally accessible. Operands SIGC, SIG0, BETA0 These are the internal parameters of the model that define the shape of the load surface (with one exception), see [:ref:`R7.01.42 `]. It is based on the Von Mises stress and the exponential stress tensor and compares well to experimental results on concrete under biaxial loading. The parameter BETA0 is more numerical in nature and only has the advantage of making the elasticity domain limited, including for hydrostatic compressions; the default value fulfills this function well, without affecting the shape of the domain in the areas of interest. Again, the DEFI_MATER_GC [:external:ref:`U4.42.07 `] command makes it possible to identify these parameters from quantities that are experimentally accessible (limits in tension and compression). Operand REST_RIGIDITE ~~~~ Stiffness restoration is active for compression deformation directions. To avoid an abrupt change in regime during the transition from traction to compression, a function S regulates the jump in stiffness, see [:ref:`R7.01.42 `]. The parameter REST_RIGIDITE, which is positive, controls this regularization; it corresponds to the coefficient gamma of the function S'. A value of 0 leads to not restoring stiffness (i.e. the model is without restoring stiffness) while a very large value amounts to virtually eliminating regularization. The DEFI_MATER_GC [:external:ref:`U4.42.07 `] command makes it easy to quantify this parameter by stipulating what proportion of the stiffness is restored for a deformation level corresponding to the initial compression threshold. .. _RefNumPara__121345_1070254933: