Modélisation A ============== Caractéristiques de la modélisation ----------------------------------- **AXIS-FOURIER (TRIA6)** .. image:: images/10000E0200001EE700000C67A8E67BB3FD0E328A.svg :width: 398 :height: 159 .. _RefImage_10000E0200001EE700000C67A8E67BB3FD0E328A.svg: Les axes de description du maillage sont :math:`x(r)` et :math:`y(z)` . Mode - Fourier : 1 :math:`T(A)=0.` .. csv-table:: ":math:`S=-3.`", "sur tout le domaine" ":math:`[\mathrm{BC}]` :", ":math:`\phi =2.`" ":math:`[\mathrm{CD}]` :", ":math:`h=2.` :math:`{T}_{\mathrm{ext}}=2.`" Caractéristiques du maillage ---------------------------- Nombre de nœuds : 25. Nombre de mailles et types : 8 TRIA6 Remarques --------- Le numéro du mode de Fourier n'affectant pas le chargement, le mot clé MODE_FOURIER n'est pas nécessaire dans la commande CALC_VECT_ELEM. L'utilisation de la commande CREA_CHAMP / ASSE n'est pas une recombinaison de Fourier mais une simple validation de ce mot clé. Valeurs testées --------------- .. csv-table:: "", "**Identification**", "**Référence**" ":math:`\theta =0`", ":math:`T(B)`", "1." "", ":math:`T(F)`", "0.25" "", ":math:`{\phi }_{r}(B)`", "–2" "", ":math:`{\phi }_{r}(F)`", "–1." "", ":math:`{\phi }_{\theta}(B)`", "1." "", ":math:`{\phi }_{\theta}(F)`", "0.5" "", ":math:`{\phi }_{z}(B)`", "0." "", ":math:`{\phi }_{z}(F)`", "0." "", "", "" ":math:`\theta =45`", ":math:`T(B)`", "0.7071" "", ":math:`T(F)`", "0.177" "", ":math:`{\phi }_{r}(B)`", "–1.414" "", ":math:`{\phi }_{r}(F)`", "–0.7071" "", ":math:`{\phi }_{\theta}(B)`", "–0.707" "", ":math:`{\phi }_{\theta}(F)`", "–0.3535" "", ":math:`{\phi }_{z}(B)`", "0." "", ":math:`{\phi }_{z}(F)`", "0." "", "", "" ":math:`\theta =135`", ":math:`T(B)`", "–0.707" "", ":math:`T(F)`", "–0.177" "", ":math:`{\phi }_{r}(B)`", "1.414" "", ":math:`{\phi }_{r}(F)`", "0.707" "", ":math:`{\phi }_{\theta}(B)`", "–0.707" "", ":math:`{\phi }_{\theta}(F)`", "–0.3535" "", ":math:`{\phi }_{z}(B)`", "0." "", ":math:`{\phi }_{z}(F)`", "0." Remarques --------- Les valeurs des flux aux nœuds sont moyennées sur les éléments contenant ce nœud. On remarque que la solution exacte n'est pas trouvée. Ceci est du au fait que l'intégration numérique de la matrice de rigidité thermique est approchée (formule à 3 points de GAUSS). Si on utilisait une formule à 6 points, on trouverait la solution exactement.