Suiv.: pent3r1d Sup.: 6 Fiches techniques Préc.: hexa3q2d Index Table des matières

hexa3q2c

DESCRIPTION DE L'ELEMENT

• Dimension de l'espace : 3, élément de volume
• Inconnues : , et , les déplacements suivant x, y et z
• Formulation variationnelle : en 'déplacements'
• Numéro de code : 200 018
• numéro de code géométrique : 7 (hexaèdre)
• Interpolation : Q2-Lagrange isoparamétrique
• Nombre de points : 20
• Nombre de types de points : 1
• Les points sont définis par leurs seules coordonnées : NCACAR=0
• Nombre de noeuds : 20
• Nombre de types de noeuds : 1
• Type des noeuds : les sommets et les milieux des arêtes, 3 d.l. par noeud, 'VN' pour et
• Points et noeuds sont confondus : NCOPNP=1
• Fichier POBA : nécessaire
• Blocage (S.D. BDCL) : (noeud) = valeur
• Sous-programmes utilisateurs : oui

 
Figure: Elément HEXA 3Q2C 

Coordonnées sur l'élément de référence:

INTEGRATION NUMERIQUE

-

Intégration de volume: Intégration numérique aux 27 points de Gauss.

• npiv=27
• voir HEXA 3Q2D
• ibidem

-

Intégration de surface: Intégration numérique aux 9 points de Gauss de la face.

• npis=9
• ibidem
• ibidem

LES DONNEES A FOURNIR (en double précision)

-

S.D. MILI:

• la matrice de masse , partie
  • la chaleur massique supposée constante par élément.
• la matrice de rigidité , partie

  le paramètre NOTEL permet de spécifier la nature du problème: isotrope (> 0), anisotrope (< 0).

  • matériau homogène isotrope : k
  • matériau anisotrope :
• la matrice de rigidité , partie
  • g(k) le coefficient de transfert aux npis points d'intégration
• le tableau des flux : comme pour pour si elle n'est pas calculée, rien sinon.

-

S.D. FORC: NDSM seconds membres.

• terme de volume , partie
  • pour i=1,NDSM; k=1,npiv; la composante j du cas de charge i au point d'intégration k.
• terme de surface , partie
  • pour i=1,NDSM; k=1,npis; la composante j du cas de charge i au point d'intégration k.

Remarques

• Pour utiliser les sous-programmes utilisateurs, on fixera la valeur de NOTEL à la valeur ci-dessus en ajoutant 20 en isotrope (en retranchant 20 en orthotrope et anisotrope). Cette valeur est demandée lors de la construction des S.D. MILI ou FORC (cf [Guide Modulef - 4]).
• Le fichier POBA ets organisé comme suit:
  • tableau 3Q25 (double précision, 2274 valeurs) : NDIM=3,NBPOLY=20, NPIV=27,1,1,1,
  • tableau 2Q25 (double précision, 249 valeurs) : idem pour les faces quadrangulaires avec donc NDIM=2,NBPOLY=8,NPIS=9
  • tableau 7F25 (double précision)

  • tableau 3T2I (simple précision, 38 valeurs): NPI=27, NPIS=9, (les NPI poids de la formule d intégration numérique), (les NPIS poids sur le carré ,en double précision)
  • tableau 2T2D (simple précision, 144 valeurs): valeurs des dérivées des 8 polynômes de base du quadrangle Q2 aux NPIS points d intégration numérique
  • tableau 2T2D (simple précision, 1080 valeurs): coefficients des polynômes de base de Q2(R3) aux points d'intégration sur les faces de l'hexaèdre