Suiv.: 5.3.2 Construction du maillage 3D Sup.: 5.3 Un exemple en dimension 3 Préc.: 5.3 Un exemple en dimension 3 Index Table des matières

Analyse descendante et maillage de la section de départ

 
Figure: Coupe d'explication de la géométrie 

Le domaine est, topologiquement parlant, un cylindre. Par suite, la donnée du maillage de l'une de ses sections suffit à la construction du maillage 3D cherché.

La figure 5.4 montre une demie coupe dans le plan de symétrie du domaine. Elle permet de voir comment définir la montée du maillage 2D et, de ce fait, indique comment définir les différentes sections. On a donc successivement:

• la base : maillage 2D avec côte z mise à 0;
• sections 0 à 2 : translation de la base de vecteur (0.,0.,2.5);
• sections 2 à 4 : interpolation équirépartie entre la section 2 déjà créée et la section 4 définie maintenant comme image du maillage 2D via la fonction XYZ23;
• sections 4 à 5 : la section 5 est définie comme la dilatée de la section 4;
• sections 5 à 9 : translation de la section 5 de vecteur (0.,0.,4.) couplée avec une rotation de 5 degrés autour de l'axe des z;
• création du tableau ZINT de longueur 23;
• remplissage de ce tableau pour les indices 10 à 12;
• sections 9 à 12 : interpolation dirigée par la côte (tableau ci-dessus) entre la section 9 déjà créée et la section 12 définie maintenant comme image du maillage 2D via la fonction XYZ23;
• sections 12 à 16 : ces sections sont définies comme image du maillage 2D via la fonction XYZ23;
• sections 16 à 17 : la section 17 est définie comme la dilatée de la section 16;
• sections 17 à 20 : ces sections sont définies, de proche en proche, via la fonction XYZ33;
• sections 20 à 22 : la section 21 est définie comme la dilatée de la section 20, la section 22 est la dilatée de la section 21;

Les données permettant de créer le maillage 2D d'une section représentative (via le module APNOPO et son préprocesseur APNOXX). Le maillage obtenu est montré sur la figure 5.5.

'TEST_1_MA2D3D                                                       '
 COURBES
     1                                 $ IMPRE
 COURBE01(X,Y)=
       X*X+Y*Y-4.;                                                     
 COURBE02(X,Y)=
       X*X+Y*Y-25.;                                                    
FIN
'POINTS                                                              '
     1     5                            $ IMPRE NPOINT  $
$   NOP   NOREF(NOP)     X(NOP).            Y(NOP).  $
     1       1         -.200000E+01     0.000000E+00
     2       1         0.200000E+01     0.000000E+00
     3       2         0.500000E+01     0.000000E+00
     4       2         -.500000E+01     0.000000E+00
     5       0         0.000000E+00     0.000000E+00
'LIGNES                                                              '
     1     4                            $ IMPRE  NDLM   $
$   NOLIG NOELIG NEXTR1 NEXTR2 NOREFL NFFRON       RAISON $
       1      7      2      1      1     -3     0.100000E+01 
       5                                $ NOCE
       2      4      2      3      0      0     0.100000E+01
       3      7      3      4      2     -3     0.100000E+01
       5                                $ NOCE
       4      4      4      1      0      0     0.100000E+01
'QUAC                                                                '
     1     0     1     4    -1          $ IMPRE NIVEAU NUDSD NBRELI NS1L
 $ LISTE DES LIGNES DU CONTOUR :
     1     2     3     4
     7     1                            $ IMAX NQUAD
'SYMD                                                                '
     1     0     2                      $ IMPRE NIVEA1 NIVEA2 
     0     0                            $ NBNNF NBNNSD
  0.00000E+00  0.10000E+01  0.00000E+00 $ A. B. C.
'RECO                                                                '
    1    0    2    1  0.10000E-02    1  $ IMP NIV1 NIV2 NIV3 EPS IOPT
     0     0                            $ NBNNF NBNNSD
'SAUV                                                                '
     1     1     0                      $ IMPRE NINOPO NTNOPO
 fusee2d.nopo                          
 $ NOM FICHIER                          
'FIN                                                                 '

 
Figure: Le maillage 2D de départ