Comportements métallo-mécaniques ================================ Pour le comportement métallurgique (cf. [:ref:`R4.04.01 `]), trois lois de comportement sont disponibles : une loi caractéristique des transformations métallurgiques de l’acier, de l'acier avec phases de revenus et une loi caractéristique des alliages de zirconium. L’acier peut comporter (au plus) sept phases métallurgiques différentes: * quatre phases froides (notées :math:`\alpha_\alpha`) que sont la ferrite, la perlite, la bainite et la martensite * une phase chaude, l'austénite (noté :math:`\gamma`) * deux phases froides de revenu pour la bainite et la martensite Le zircaloy peut comporter (au plus) trois phases métallurgiques différentes : * Deux phases froides: une phase pure :math:`\alpha` et une phase mélangée :math:`\alpha` * Une phase chaude :math:`\beta` Pour le comportement mécanique avec la prise en compte des transformations métallurgiques, il existe deux modèles. Le premier modèle (cf. [:ref:`R4.04.02 `]) est utilisable pour l’acier et pour le Zircaloy. On choisit le matériau souhaité en activant, dans le mot-clef facteur :ref:`COMPORTEMENT ` des opérateurs non-linéaires, le mot-clé RELATION_KIT qui vaut ACIER ou ZIRC. Les différentes relations relatives à ce modèle sont identiques pour ces deux matériaux (on traite les mêmes phénomènes) mais le nombre de phases en présence est différent. Le second modèle (cf. [:ref:`R4.04.05 `]) est uniquement disponible pour le Zircaloy (``RELATION_KIT=’ZIRC’``)et correspond au mot clé META_LEMA_ANI sous :ref:`COMPORTEMENT `. Mot clé facteur META_ACIER ------------------------------ Paramètres à renseigner pour la métallurgie de l’acier. Paramètres pour les changements de phases ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ TRC est un concept de type table produit par l’opérateur [:ref:`DEFI_TRC `] et contenant l’ensemble des informations fournies par les diagrammes TRC (Transformation en Refroidissement Continu) de l’acier considéré. AR3 est la température quasi-statique de début de décomposition de l’austénite au refroidissement. Remarque Il faut être cohérent entre cette valeur et celle correspondante à la courbe de refroidissement donnée par la commande [:ref:`DEFI_TRC `] . Au-delà d'une différence de 10°C, une alarme sera émise et on considère que les résultats pourraient être faux. ALPHA est le coefficient :math:`\alpha` de la loi de Koïstinen-Marbürger exprimant la quantité de martensite formée en fonction de la température .. math:: {Z}_{\alpha_{m}}=1-\exp(\alpha ({M}_{s}-T)) MSO est la température de début de transformation :math:`{M}_{s,0}` martensitique lorsque celle-ci est totale. Dans ce cas :math:`{M}_{s}={M}_{s0}`. AC1 est la température quasi-statique de début de transformation en austénite au chauffage. AC3 est la température quasi-statique de fin de transformation en austénite. TAUX_1 est la valeur de la fonction "retard" (cf. [:ref:`R4.04.01 `]) :math:`\tau_1(T)` intervenant dans le modèle de transformation austénitique à la température AC1. TAUX_3 est la valeur de la fonction "retard" (cf. [:ref:`R4.04.01 `]) :math:`\tau_3(T)` intervenant dans le modèle de transformation austénitique à la température AC3. L’évolution de la proportion d’austénite est alors définie par : .. math:: \dot{Z}=\frac{Z-{Z}_{\text{eq}}\left(T\right)}{\tau \left(T\right)} avec :math:`{Z}_{\text{eq}}\left(T\right)` qui a donc la forme de [:ref:`Courbe `]. .. figure:: ../../../man_u/u4/u4.43.01/images/Forme1.gif :name: u4.43.01-ZEQForme Courbe :math:`{Z}_{\text{eq}}\left(T\right)` et :math:`\tau \left(T\right)` qui a la forme de [:ref:`Courbe `]. .. figure:: ../../../man_u/u4/u4.43.01/images/Forme2.gif :name: u4.43.01-TAUForme Courbe :math:`\tau \left(T\right)` Paramètres pour la taille de grain ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Les quatre paramètres suivants entraînent le calcul de taille de grains s’ils sont renseignés: .. math:: \frac{d}{{dt}}(d)=\frac{1}{\lambda}\left (\frac{1}{d}-\frac{1}{{d}_{\lim}} \right) avec :math:`\lambda ={\lambda}_{0}\cdot \exp \left(\frac{{Q}_{\text{app}}}{R T} \right )` et :math:`{d}_{\lim}={d}_{10}\exp \left (-\frac{{W}_{\text{app}}}{R T} \right )` . LAMBDA0 est le paramètre matériau :math:`{\lambda}_{0}` QSR_K est le paramètre d'énergie d’activation :math:`\frac{{Q}_{\text{app}}}{R T}` D10 est le paramètre matériau :math:`{d}_{10}` WSR_K est le paramètre d'énergie d’activation :math:`\frac{{W}_{\text{app}}}{R T}` Mot clé facteur META_ACIER_REVENU ------------------------------------- Paramètres à renseigner pour la métallurgie de l’acier avec ses phases de revenu. Paramètres pour les changements de phases ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ On considère une loi de type Johnson-Mehl-Avrami pour décrire l’évolution du pourcentage de phase revenue :math:`{\tilde{Z}}^{r}_{\alpha_{\alpha}}` en fonction du temps de maintien :math:`t^{*}` à une température :math:`T^{m}` pour les deux phases de revenu (bainite et martensite). .. math:: {\tilde{Z}}^{r}_{\alpha_{\alpha}}(t^{*})=\left ( 1 - \exp\left\lbrace -b ({t^{*}})^n \right\rbrace \right ) TEMP est la température de revenu :math:`T^{r}` et TEMP_MAINTIEN est la température de maintien :math:`T^{m}` intervenant dans l'identification des paramètres. BAINITE_B et BAINITE_N sont les coefficients :math:`b` et :math:`n` pour le cas de la bainite :math:`{\alpha_{\alpha}}=b`. MARTENSITE_B et MARTENSITE_N sont les coefficients :math:`b` et :math:`n` pour le cas de la martensite :math:`{\alpha_{\alpha}}=m`. Mot clé facteur META_ZIRC ----------------------------- Paramètres à renseigner pour la métallurgie du Zircaloy (cf. [:ref:`R4.04.04 `]). .. code-block:: text META_ZIRC = _F( ◆ TDEQ = float, ◆ N = float, ◆ K = float, ◆ T1C = float, ◆ T2C = float, ◆ AC = float, ◆ M = float, ◆ QSR_K = float, ◆ T1R = float, ◆ T2R = float, ◆ AR = float, ◆ BR = float, ) On rapppelle que le matériau contient deux phases: * :math:`\alpha` : phase à froid hexagonale compacte * :math:`\beta` : phase à chaud cubique centrée TDEQ est la température de début de transformation :math:`\alpha \iff \beta` à l’équilibre. N et K sont des paramètres matériau relatifs au modèle donnant la proportion de :math:`\beta` en fonction de la température, à l’équilibre. T1C est la température de début de transformation :math:`\alpha` en :math:`\beta` au chauffage. T2C est un paramètre matériau intervenant dans le calcul de la température de début de transformation :math:`\alpha` en :math:`\beta` au chauffage. AC et M sont des paramètres matériau intervenant dans le modèle d’évolution de :math:`\beta` au chauffage. T2R est un paramètre matériau intervenant dans le calcul de la température de début de transformation :math:`\beta` en :math:`\alpha` au refroidissement. AR et BR sont des paramètres matériau intervenant dans le modèle d’évolution de :math:`\beta` au refroidissement. QSR_K est la constante d’Arrhénius exprimé en degré Kelvin. Mot-clé facteur DURT_META ----------------------------- Définition des caractéristiques relatives au calcul de dureté associée à la métallurgie des aciers. La dureté est calculée en utilisant une loi de mélange linéaire sur la micro-dureté des constituants : Le modèle s’écrit :math:`H_{V}=\sum_{k}{Z}_{k}{H_{V}}_{k}` avec :math:`H_{V}` la dureté (ici Vickers par exemple) du point polyphasé, :math:`{Z}_{k}` la proportion de la phase :math:`k` et :math:`{H_{V}}_{k}` la dureté de la phase :math:`k` . F1_DURT est la micro-dureté de la phase à froid :math:`\mathrm{F1}` (ferrite pour l’acier). F2_DURT est la micro-dureté de la phase à froid :math:`\mathrm{F2}` (perlite pour l’acier). F3_DURT est la micro-dureté de la phase à froid :math:`\mathrm{F3}` (bainite pour l’acier). F4_DURT est la micro-dureté de la phase à froid :math:`\mathrm{F4}` (martensite pour l’acier). C_DURT est la micro-dureté pour la phase à chaud (austénite pour l’acier). F3_REVENU_DURT est la micro-dureté de la phase *revenue* à froid :math:`\mathrm{F3}` (bainite pour l’acier). F4_REVENU_DURT est la micro-dureté de la phase *revenue* à froid :math:`\mathrm{F4}` (martensite pour l’acier). Mots clés facteur ELAS_META, ELAS_META_FO ----------------------------------------- Définition des caractéristiques élastiques, de dilatation et de limites d’élasticité pour la modélisation d’un matériau subissant des transformations métallurgiques (voir [:ref:`R4.04.02 `] ou [:ref:`R4.04.05 `]). Ces coefficients peuvent être soient constants par rapport à la température ELAS_META, soient dépendre de la température ELAS_META_FO (paramètre ’TEMP’). Certains coefficients dépendent de la structure métallurgique (paramètre ’META’). **Remarque:** Concernant le modèle META_LEMA_ANI, la dilatation thermique s’écrit classiquement sans distinction des phases. Par conséquent, les mots clés ‘C_ALPHA’, ‘PHASE_REFE’ et ‘EPSF_EPSC_TREF’ sont obligatoires mais non pris en compte dans les équations. Seul le coefficient de dilatation ‘F_ALPHA’ est considéré. Ce modèle est une loi sans seuil donc les limites d’élasticité et la loi des mélanges ne sont pas utiles. **Remarque:** Concernant les autres modèles, pour un acier on renseigne au maximum 5 limites d’élasticité, pour le Zircaloy on en renseigne au maximum trois. Opérandes ~~~~~~~~~ .. code-block:: text E = young Module de Young, identique pour toutes les phases métallurgiques. .. code-block:: text NU = nu Coefficient de Poisson, identique pour toutes les phases métallurgiques. .. code-block:: text F_ALPHA = fal Coefficient de dilatation thermique moyen des phases froides. .. code-block:: text C_ALPHA = cal Coefficient de dilatation thermique moyen de la phase chaude. .. code-block:: text PHASE_REFE = / ’CHAUD’ / ’FROID’ Choix de la phase métallurgique de référence (phase chaude ou phase froide). En effet, pour définir la déformation thermique nulle, il faut définir la température de référence :math:`{T}_{\mathrm{ref}}` (définie dans AFFE_MATERIAU) et la phase métallurgique de référence, de sorte que la déformation thermique soit considérée nulle à :math:`{T}_{\mathrm{ref}}` et dans l’état métallurgique de référence. .. code-block:: text EPSF_EPSC_TREF = deltae Déformation de la phase non de référence par rapport à la phase de référence à la température :math:`{T}_{\mathrm{ref}}` : traduit la différence de compacité entre les structures cristallographiques cubiques à faces centrées (type austénitique) et cubiques centrées (type ferritique). .. code-block:: text TEMP_DEF_ALPHA = Tda Température par rapport à laquelle on définit le coefficient de dilatation. Dans le cas où C_ALPHA est une fonction, cet opérande est obligatoire. .. code-block:: text PRECISION = eps Ce réel indique avec quelle précision une température :math:`T` est proche de la température de référence (cf. [:ref:`§3.1.4 <§3.1.4>`]). .. code-block:: text F1_SY = F1sy Limite d’élasticité de la phase froide 1 pour un comportement plastique. .. code-block:: text F2_SY = F2sy Limite d’élasticité de la phase froide 2 pour un comportement plastique. .. code-block:: text F3_SY = F3sy Limite d’élasticité de la phase froide 3 pour un comportement plastique. .. code-block:: text F4_SY = F4sy Limite d’élasticité de la phase froide 4 pour un comportement plastique. .. code-block:: text C_SY = Fsy Limite d’élasticité de la phase chaude pour un comportement plastique. .. code-block:: text SY_MELANGE = f Fonction utilisée pour la loi de mélange sur la limite d’élasticité du matériau multiphasé pour un comportement plastique. :math:`{\sigma}_{y}=\left(1-f\left(z\right)\right){\sigma}_{y}^{\gamma}+f\left(z\right){\sigma}_{y}^{\alpha}` .. code-block:: text F1_S_VP = F1svp Limite d’élasticité de la phase froide 1 pour un comportement visqueux. .. code-block:: text F2_S_VP = F2svp Limite d’élasticité de la phase froide 2 pour un comportement visqueux. .. code-block:: text F3_S_VP = F3svp Limite d’élasticité de la phase froide 3 pour un comportement visqueux. .. code-block:: text F4_S_VP = F4svp Limite d’élasticité de la phase froide 4 pour un comportement visqueux. .. code-block:: text C_S_VP = Csvp Limite d’élasticité de la phase chaude pour un comportement visqueux. .. code-block:: text S_VP_MELANGE = Svp Fonction utilisée pour la loi de mélange sur la limite d’élasticité du matériau multiphasé pour un comportement visqueux. :math:`{\sigma}_{c}=\left(1-f\left(z\right)\right){\sigma}_{c}^{\gamma}+f\left(z\right){\sigma}_{c}^{\alpha}` .. _RefNumPara__120386_1070254933: Mot clé facteur META_ECRO_LINE ------------------------------ .. code-block:: text Définition de cinq modules d’écrouissage utilisés dans la modélisation du phénomène d’écrouissage isotrope linéaire d’un matériau subissant des changements de phases métallurgiques (voir [:ref:`R4.04.02 `]). Ces modules dépendent de la température. Opérandes ~~~~~~~~~ .. code-block:: text F1_D_SIGM_EPSI = dsde1 Pente de la courbe de traction pour la phase froide 1. .. code-block:: text F2_D_SIGM_EPSI = dsde2 Pente de la courbe de traction pour la phase froide 2. .. code-block:: text F3_D_SIGM_EPSI = dsde3 Pente de la courbe de traction pour la phase froide 3. .. code-block:: text F4_D_SIGM_EPSI = dsde4 Pente de la courbe de traction pour la phase froide 4. .. code-block:: text C_D_SIGM_EPSI = dsdec Pente de la courbe de traction pour la phase chaude. .. image:: images/10000000000002D6000001E90DCCCD948845BAFA.png :width: 3.3508in :height: 2.3102in .. _RefImage_10000000000002D6000001E90DCCCD948845BAFA.png: Le module de Young :math:`E` est à préciser par les mots-clés META_ELAS ou META_ELAS_FO. .. _RefNumPara__120388_1070254933: Mot clé facteur META_TRACTION ----------------------------- Définition de cinq courbes de traction utilisées dans la modélisation du phénomène d’écrouissage isotrope non linéaire d’un matériau subissant des changements de phases métallurgiques (voir [:ref:`R4.04.02 `]). Les courbes de traction peuvent éventuellement dépendre de la température. Opérandes ~~~~~~~~~ .. code-block:: text SIGM_F1 = r_p1 .. code-block:: text Courbe écrouissage isotrope :math:`R`en fonction de la déformation plastique cumulée :math:`p`pour la phase froide 1 . .. code-block:: text SIGM_F2 = r_p2 .. code-block:: text Courbe écrouissage isotrope :math:`R`en fonction de la déformation plastique cumulée :math:`p`pour la phase froide 2 . .. code-block:: text SIGM_F3 = r_p3 .. code-block:: text Courbe écrouissage isotrope :math:`R`en fonction de la déformation plastique cumulée :math:`p`pour la phase froide 3 . .. code-block:: text SIGM_F4 = r_p4 .. code-block:: text Courbe écrouissage isotrope :math:`R`en fonction de la déformation plastique cumulée :math:`p`pour la phase froide 4 . .. code-block:: text SIGM_C = r_p c .. code-block:: text Courbe écrouissage isotrope :math:`R`en fonction de la déformation plastique cumulée :math:`p`pour la phase chaude. **Remarque:** *Attention il ne s’agit pas de la courbe* :math:`\sigma` *fonction de* :math:`\varepsilon` *mais de la courbe* :math:`R` *fonction de* :math:`p` \*. On passe de l’une à l’autre en effectuant les calculs suivants:* :math:`R=\sigma -\mathrm{limite}d’\mathrm{élasticité},p=\varepsilon –(\sigma /E).` .. _RefNumPara__120390_1070254933: Mot clé facteur META_VISC_FO ---------------------------- Définition des paramètres visqueux de la loi de comportement viscoplastique avec prise en compte de la métallurgie (voir [:ref:`R4.04.02 `]). Le modèle viscoplastique de type Norton-Hoff comporte 5 paramètres; les paramètres classique :math:`\eta` , :math:`n` de la loi d’écoulement en puissance, la limite élastique d’écoulement visqueuse, les paramètres :math:`C` et :math:`m` relatifs à la restauration d’écrouissage d’origine visqueuse. Ces paramètres dépendent de la température et de la structure métallurgique. Les paramètres limites d’élasticité sont définis dans ELAS_META. Opérandes F1_ETA/F2_ETA/F3_ETA/F4_ETA/C_ETA ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ .. code-block:: text F1_ETA = eta1 Paramètre :math:`\eta` de la loi d’écoulement visco-plastique, pour la phase froide1. .. code-block:: text F2_ETA = eta2 Paramètre :math:`\eta` de la loi d’écoulement visco-plastique, pour la phase froide2. .. code-block:: text F3_ETA = eta3 Paramètre :math:`\eta` de la loi d’écoulement visco-plastique, pour la phase froide3. .. code-block:: text F4_ETA = eta4 Paramètre :math:`\eta` de la loi d’écoulement visco-plastique, pour la phase froide4. .. code-block:: text C_ETA = etac Paramètre :math:`\eta` de la loi d’écoulement visco-plastique, pour la phase à chaud. Opérandes F1_N/F2_N/F3_N/F4_N/C_N ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ .. code-block:: text F1_N = n1 Paramètre :math:`n` de la loi d’écoulement visco-plastique, pour la phase froide 1. .. code-block:: text F2_N = n2 Paramètre :math:`n` de la loi d’écoulement visco-plastique, pour la phase froide 2. .. code-block:: text F3_N = n3 Paramètre :math:`n` de la loi d’écoulement visco-plastique, pour la phase froide 3. .. code-block:: text F4_N = n4 Paramètre :math:`n` de la loi d’écoulement visco-plastique, pour la phase froide 4. .. code-block:: text C_N = n5 Paramètre n de la loi d’écoulement visco-plastique, pour la phase à chaud. Opérandes F1_C/F2_C/F3_C/F4_C/C_C ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ .. code-block:: text F1_C = C1 Paramètre :math:`C` relatif à la restauration d’écrouissage d’origine visqueuse, pour la phase froide 1. .. code-block:: text F2_C = C2 Paramètre :math:`C` relatif à la restauration d’écrouissage d’origine visqueuse, pour la phase froide 2. .. code-block:: text F3_C = C3 Paramètre :math:`C` relatif à la restauration d’écrouissage d’origine visqueuse, pour la phase froide 3. .. code-block:: text F4_C = C4 Paramètre :math:`C` relatif à la restauration d’écrouissage d’origine visqueuse, pour la phase froide 4. .. code-block:: text C_C = C5 Paramètre :math:`C` relatif à la restauration d’écrouissage d’origine visqueuse, pour la phase à chaud. Opérandes F1_M/F2_M/F3_M/F4_M/C_M ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ .. code-block:: text F1_M = m1 Paramètre :math:`m` relatif à la restauration d’écrouissage d’origine visqueuse, pour la phase froide 1. .. code-block:: text F2_M = m2 Paramètre :math:`m` relatif à la restauration d’écrouissage d’origine visqueuse, pour la phase froide 2. .. code-block:: text F3_M = m3 Paramètre :math:`m` relatif à la restauration d’écrouissage d’origine visqueuse, pour la phase froide 3. .. code-block:: text F4_M = m4 Paramètre :math:`m` relatif à la restauration d’écrouissage d’origine visqueuse, pour la phase froide 4. .. code-block:: text C_M = m5 Paramètre :math:`m` relatif à la restauration d’écrouissage d’origine visqueuse, pour la phase à chaud. .. _RefNumPara__120392_1070254933: Mot clé facteur META_PT ----------------------- Définition des caractéristiques utilisées dans la modélisation de la plasticité de transformation d’un matériau qui subit des changements de phases métallurgiques (voir [:ref:`R4.04.02 `]). Le modèle est le suivant : :math:`\Delta {\varepsilon}^{\mathrm{pt}}=\frac{3}{2}\sigma \sum_{i=1}^{i=4}{K}_{i}{F}_{i}^{'}({Z}_{i})\langle \Delta {Z}_{i}\rangle` Opérandes ~~~~~~~~~ .. code-block:: text F1_K = Kf, F2_K = Kp, F3_K = Kb, F4_K = Km Constantes :math:`{K}_{i}` utilisées dans le modèle de plasticité de transformation, pour les différentes phases à froid. Pour l’acier = phase ferritique, perlitique, bainitique et martensitique. .. code-block:: text F1_D_F_META=F’f, F2_D_F_META=F’p, F3_D_F_META=F’b, F4_D_F_META=F’m, Fonctions :math:`{F}_{i}^{'}` utilisées dans le modèle de plasticité de transformation, pour les différentes phases à froid. Pour l’acier : phase ferritique, perlitique, bainitique et martensitique. .. _RefNumPara__120394_1070254933: Mot clé facteur META_RE ----------------------- Définition des caractéristiques utilisées dans la modélisation du phénomène de restauration d’écrouissage d’un matériau qui subit des changements de phases métallurgiques (voir [:ref:`R4.04.02 `]). Opérandes ~~~~~~~~~ .. code-block:: text C_F1_THETA=Tgf, C_F2_THETA=Tgp, C_F3_THETA=Tgb, C_F4_THETA=Tgm Constantes caractérisant le taux d’écrouissage transmis lors de la transformation de la phase à chaud C en phase à froid. Pour l’acier; transformation de l’austénite en ferrite, perlite, bainite et martensite. Ainsi, :math:`\theta =0` correspond à une restauration totale et :math:`\theta =1` à une transmission totale de l’écrouissage. .. code-block:: text F1_C_THETA=Tfg, F2_C_THETA=Tpg, F3_C_THETA=Tbg, F4_C_THETA=Tmg Constantes caractérisant le taux d’écrouissage transmis lors de la transformation des phases à froid en phase à chaud. Pour l’acier; transformation de la ferrite, de la perlite, de la bainite et de la martensite en austénite. Ainsi, :math:`\theta =0` correspond à une restauration totale et :math:`\theta =1` à une transmission totale de l’écrouissage. .. _RefNumPara__120396_1070254933: Mot clé META_LEMA_ANI --------------------- Définition des paramètres de la loi META_LEMA_ANI (cf. [:ref:`R4.04.05 `]), élasto-visqueuse sans seuil avec un comportement anisotrope. Brièvement, le modèle s’écrit soit dans le repère cylindrique :math:`(r,\theta ,z)` , soit dans le repère cartésien :math:`(\mathit{Ox},\mathit{Oy},\mathit{Oz})` : Partition des déformations : :math:`\epsilon ={\epsilon}^{e}+\alpha \Delta T\mathrm{Id}+{\epsilon}^{v}` Loi d’écoulement de la déformation visqueuse : :math:`\dot{{\epsilon}^{v}}=\dot{p}\frac{M:\sigma }{{\sigma}_{\mathit{eq}}}` Critère de Hill: :math:`{\sigma}_{\mathit{eq}}=\sqrt{\sigma :M:\sigma }` Matrice de Hill :math:`\mathrm{M}` *en coordonnées cylindriques* : :math:`{\underline{M}}_{(r,\theta ,z)}=\left[\begin{array}{cccccc}{M}_{\mathit{rrrr}}& {M}_{\mathit{rr}\theta \theta }& {M}_{\mathit{rrzz}}& 0& 0& 0\\ {M}_{\mathit{rr}\theta \theta }& {M}_{\theta \theta \theta \theta }& {M}_{\theta \theta zz}& 0& 0& 0\\ {M}_{\mathit{rrzz}}& {M}_{\theta \theta zz}& {M}_{\mathit{zzzz}}& 0& 0& 0\\ 0& 0& 0& {M}_{r\theta r\theta }& 0& 0\\ 0& 0& 0& 0& {M}_{\mathit{rzrz}}& 0\\ 0& 0& 0& 0& 0& {M}_{\theta z\theta z}\end{array}\right]` :math:`\text{avec}\lbrace \begin{array}{ccc}{M}_{\mathit{rrrr}}+{M}_{\mathit{rr}\theta \theta }+{M}_{\mathit{rrzz}}& =& 0\\ {M}_{\mathit{rr}\theta \theta }+{M}_{\theta \theta \theta \theta }+{M}_{\theta \theta zz}& =& 0\\ {M}_{\mathit{rrzz}}+{M}_{\theta \theta zz}+{M}_{\mathit{zzzz}}& =& 0\end{array}` ou bien, en coordonnées cartésiennes: :math:`{\underline{M}}_{(x,y,z)}=\left[\begin{array}{cccccc}{M}_{\mathit{xxxx}}& {M}_{\mathit{xxyy}}& {M}_{\mathit{xxzz}}& 0& 0& 0\\ {M}_{\mathit{xxyy}}& {M}_{\mathit{yyyy}}& {M}_{\mathit{yyzz}}& 0& 0& 0\\ {M}_{\mathit{xxzz}}& {M}_{\mathit{yyzz}}& {M}_{\mathit{zzzz}}& 0& 0& 0\\ 0& 0& 0& {M}_{\mathit{xyxy}}& 0& 0\\ 0& 0& 0& 0& {M}_{\mathit{xzxz}}& 0\\ 0& 0& 0& 0& 0& {M}_{\mathit{yzyz}}\end{array}\right]` :math:`\text{avec}\lbrace \begin{array}{ccc}{M}_{\mathit{xxxx}}+{M}_{\mathit{xxyy}}+{M}_{\mathit{xxzz}}& =& 0\\ {M}_{\mathit{xxyy}}+{M}_{\mathit{yyyy}}+{M}_{\mathit{yyzz}}& =& 0\\ {M}_{\mathit{xxzz}}+{M}_{\mathit{yyzz}}+{M}_{\mathit{zzzz}}& =& 0\end{array}` Loi des mélanges sur la matrice :math:`\mathrm{M}` : :math:`\begin{array}{c}M=\lbrace \begin{array}{ccc}{M}^{c}& \text{si}& 0.00⩽{Z}_{f}⩽0.01\\ {M}^{2}={Z}_{f}{M}^{1}+(1-{Z}_{f}){M}^{c}& \text{si}& 0.01⩽{Z}_{f}⩽0.99\\ {M}^{1}& \text{si}& 0.99⩽{Z}_{f}⩽1.00\end{array}\\ {Z}_{f}={Z}_{1}+{Z}_{2};\text{}{Z}_{c}={Z}_{3}=1-{Z}_{f}\end{array}` Vitesse de déformation équivalente: :math:`\dot{p}={\left(\frac{{\sigma}_{\mathit{eq}}}{{\mathit{ap}}^{m}}\right)}^{n}{e}^{-Q/T}` ou de manière équivalente: :math:`{\sigma}_{\mathit{eq}}=\underset{\text{contrainte visqueuse}{\sigma}_{v}}{\underset{⏟}{a{\left({e}^{Q/T}\right)}^{1/n}{p}^{m}{\dot{p}}^{1/n}}}={\sigma}_{v}` Loi des mélanges sur la contrainte visqueuse :math:`{\sigma}_{v}` : :math:`{\sigma}_{\mathit{eq}}={\sigma}_{v}=\sum_{i=1}^{3}{f}_{i}({Z}_{\alpha}){\sigma}_{\mathit{vi}}` avec :math:`{\sigma}_{\mathit{vi}}={a}_{i}{\left({e}^{{Q}_{i}/T}\right)}^{1/{n}_{i}}{p}^{{m}_{i}}{\dot{p}}^{1/{n}_{i}}` **Remarque** : dans le cas isotrope, on a: :math:`\begin{array}{}{M}_{\mathrm{rrrr}}={M}_{\theta \theta \theta \theta }={M}_{\mathrm{zzzz}}=1\\ {M}_{\mathrm{r\theta r\theta }}={M}_{\mathrm{rzrz}}={M}_{\theta z\theta z}=0.75\end{array}` :math:`\begin{array}{c}{M}_{\mathit{xxxx}}={M}_{\mathit{yyyy}}={M}_{\mathit{zzzz}}=1\\ {M}_{\mathit{xyxy}}={M}_{\mathit{xzxz}}={M}_{\mathit{yzyz}}=0.75\end{array}` Le choix du type de coordonnées (cylindriques ou cartésiennes) se fait respectivement par la détection du mot-clé F_MRR_RR ou du mot-clé F_MXX_XX. Opérandes ~~~~~~~~~ Le tableau ci-dessous résume les correspondances entre les symboles des équations et les mots clés. .. csv-table:: "Symbole dans les équations", "Mot clé" ":math:`\mathrm{a1}` , :math:`\mathrm{a2}` , :math:`\mathrm{a3}`", "’F1_A’, ’F2_A’, ’C_A’" ":math:`\mathrm{m1}` , :math:`\mathrm{m2}` , :math:`\mathrm{m3}`", "’F1_M’, ’F2_M’, ’C_M’" ":math:`\mathrm{n1}` , :math:`\mathrm{n2}` , :math:`\mathrm{n3}` :math:``", "’F1_N’, ’F2_N’, ’C_N’" ":math:`\mathrm{Q1}` , :math:`\mathrm{Q2}` , :math:`\mathrm{Q3}`", "’F1_Q’, ’F2_Q’, ’C_Q’" La matrice de Hill est connue soit pour la phase froide (1) ’F_M**_**’, soit pour la phase chaude (3) ’C_M**_**’. **Remarque :** *Les coefficients* ’F1_Q’, ’F2_Q’ *et* ’C_Q’ *sont en degré Kelvin.* .. _RefNumPara__121343_1070254933: