1.5.6 Exemple

Modulef
Suiv.: 1.5.7 Limitations et Bogues Sup.: Le préprocesseur EMC² Préc.: L'application EDIT_MESH Index Table des matières

Dans cette section nous montrons comment mailler un carré unité troué. Il faut commencer par exécuter le programme emc2 et entrer le numéro du "device" graphique de F3D et de plus répondre à quelques questions supplémentaires (dépend de l'implémentation de F3D).

Construction du carré unité troué

Cliquer dans l'item POINT du menu construction
Cliquer dans l'item PT_XY du menu désignation
Taper au clavier 0=0=1=0=0=1=1=1=.5=.5=, 5 points s'affichent au milieu de l'écran car l'échelle est égale à 1
Cliquer dans l'item VOIR_TOUT pour bien voir les 5 points. Le 5ème point sera le centre du trou.
Cliquer dans l'item SEGMENT du menu construction (en haut, celui du bas est l'item SEGMENT du menu de désignation)
Cliquer dans l'item POINT du menu désignation (en bas)
Cliquer dans la fenêtre graphique sur les points extremités des 4 segments à créer après chaque couple de cliques, un segment est construit; donc il faut cliquer 8 fois par exemple près des 8 points suivant (0,0), (1,0), (1,0), (1,1), (1,1), (0,1), (0,1), (0,0) dans cet ordre.
Cliquer dans l'item CERCLE du menu construction (en haut) puis cliquer dans l'item CENTRE
Cliquer dans la fenêtre graphique près un point (.5,.5) (on est toujours en désignation de point)
Cliquer dans l'item RAYON
Taper au clavier 1/4= (un cercle de centre (.5,.5) et de rayon 1/4 s'affiche), attention le cercle comme la droite n'est qu'un support de construction donc il faut le transformer en arc
Cliquer dans l'item ARC du menu pour construire des arcs
Cliquer dans l'item IDEM
Cliquer dans l'item CERCLE du menu désignation (pour désigner le cercle)
Cliquer dans la fenêtre graphique (comme il n'y a qu'un cercle on ne peut que désigner celui-ci). Le cercle ne servant plus, on le détruit.
Cliquer dans l'item DETRUIRE du menu général (à gauche)
Cliquer dans la fenêtre graphique (le cercle et l'arc disparaissent car il sont superposés). Attention on est toujours en mode de destruction.
Cliquer dans l'item RAFRAICHIR pour réafficher (l'arc est toujours la)

On a fini la construction de la géométrie (figure 1.5.6).

On peut passer à l'application PREP_MESH.

Figure 1.29: 8.125

La géométrie du domaine carré troué

Définition de la discrétisation sur la frontière

Cliquer dans l'item PREP_MESH du menu général
Cliquer dans l'item NB_INTERVAL du menu de l'application (en haut)
Taper au clavier 4=
TOUT du menu de l'application (en haut) (tous les éléments sont découpés en 4 intervalles)
Taper au clavier 12=
Puis cliquer dans la fenêtre graphique près de l'arc qui se découpe en 12 intervalles). On n'a pas besoin de changer de mode de désignation.
Cliquer dans l'item SAUVER du menu général
Taper le nom de fichier carre_troue(CR), on a créé le fichier carre_troue.emc2_bd

Figure 1.30: 8.125

La discrétisation des lignes du domaine carré troué

le fichier carre_troue.emc2_bd :

'-- TYPE N BD(1) BD(2) BD(3) BD(4) BD(5) NBNODE RAISON NUREFG NUREFD NUREF1G NUREF1D NUREF2G NUREF2D FIS'
'DROITE'      1 -1.000000 0.0000000 1.000000 0.0000000 0.0000000 0.0000000 0 0.0000000 0 0 0 0 0 0 F
'DROITE'      2 -1.000000 1.000000 0.0000000 0.0000000 0.0000000 0.0000000 0 0.0000000 0 0 0 0 0 0 F
'POINT'       3 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0 1.000000 0 0 0 0 0 0 F
'POINT'       4 0.0000000 1.000000 0.0000000 0.0000000 0.0000000 0.0000000 0 1.000000 0 0 0 0 0 0 F
'POINT'       5 0.0000000 0.0000000 1.000000 0.0000000 0.0000000 0.0000000 0 1.000000 0 0 0 0 0 0 F
'POINT'       6 0.0000000 1.000000 1.000000 0.0000000 0.0000000 0.0000000 0 1.000000 0 0 0 0 0 0 F
'POINT'       7 0.0000000 0.5000000 0.5000000 0.0000000 0.0000000 0.0000000 0 1.000000 0 0 0 0 0 0 F
'SEGMENT'     8 -3.000000 0.0000000 0.0000000 1.000000 0.0000000 0.0000000 5 1.000000 0 0 0 0 0 0 F
'SEGMENT'     9 -3.000000 1.000000 0.0000000 1.000000 1.000000 0.0000000 5 1.000000 0 0 0 0 0 0 F
'SEGMENT'     10 -3.000000 1.000000 1.000000 0.0000000 1.000000 0.0000000 5 1.000000 0 0 0 0 0 0 F
'SEGMENT'     11 -3.000000 0.0000000 1.000000 0.0000000 0.0000000 0.0000000 5 1.000000 0 0 0 0 0 0 F
'CERCLE'      12 0.2500000 0.5000000 0.5000000 0.0000000 0.0000000 0.0000000 2 1.000000 0 0 0 0 0 0 F
'ARC'         13 -2.000000 0.5000000 0.5000000 0.7500000 0.5000000 6.283185 13 1.000000 0 0 0 0 0 0 F
'MASQUE' 0 -0.4632312 1.465731 -5.0000010E-02 1.052500 0.0000000 0.0000000 0 0. 0 0 0  0 0 0  F
'RAYON' 0 0.2500000 0. 0. 0. 0. 0. 0 0. 0 0 0  0 0 0  F
'ANGLE' 0 0.0000000 0. 0. 0. 0. 0. 0 0. 0 0 0  0 0 0  F
'DISTANCE' 0 1.000000 0. 0. 0. 0. 0. 0 0. 0 0 0  0 0 0  F
'RAPPORT' 0 1.000000 0. 0. 0. 0. 0. 0 0. 0 0 0  0 0 0  F
'NOMBRE' 0 0.0000000 0. 0. 0. 0. 0. 0 0. 0 0 0  0 0 0  F
'ECHELLE' 0 14.32432 0. 0. 0. 0. 0. 0 0. 0 0 0  0 0 0  F
'NUREF' 0 0.0000000 0. 0. 0. 0. 0. 0 0. 0 0 0  0 0 0  F
'NB_INTERVALS' 0 12.00000 0. 0. 0. 0. 0. 0 0. 0 0 0  0 0 0  F
'RAISON' 0 1.000000 0. 0. 0. 0. 0. 0 0. 0 0 0  0 0 0  F
'NUDSD' 0 1.000000 0. 0. 0. 0. 0. 0 0. 0 0 0  0 0 0  F

On a fini la préparation un maillage (figure 1.5.6)

Réalisation et écriture du maillage

Cliquer dans l'item EDIT_MESH du menu général pour activer l'application
[4] EDIT_MESH puis taper 4 (CR) au clavier pour prendre les options par défaut du mailleur. Le maillage s'affiche et le trou est maillé. On va supprimer le maillage du trou:
Cliquer dans l'item S_DOM du menu désignation.
Cliquer dans l'item SUPPRIMER du menu edit_mesh (en haut).
Cliquer dans le cercle maillé, le maillage du trou disparait.
Cliquer dans l'item SAUVER puis taper le type de sauvegarde par exemple
[4] am_fmt(CR) puis le nom du fichier carre_troue(CR), un fichier carre_troue.am_fmt est créé.
Cliquer dans l'item FINIR pour sortir d'EMC²

Figure 1.31: 8.125

Le maillage final du domaine carré troué

le fichier carre_troue.am_fmt :

 44 60  -- nbs,nbt
      1      2     29        25     24     29
     24     31     29         3     43     30
     15     16     31        22     34     23
     26     25     30        23     34     31
     32     27     43        28     27     32
     22     33     34        21     40     38
      4     32     43        34     15     31
     33     14     34        21     33     22
      5      6     39        20     42     40
     28     32     39        41     42     19
      7      8     36        18     36     37
     28     39     35        34     14     15
     12     13     38         6     35     39
     28     35     17        43     27     26
     38     13     14        23     31     24
      7     36     17        31     16     29
     11     12     40         3     30      2
     37      8     44        30     25     29
     41     10     42        30     29      2
     44     10     41        29     16      1
     35      6      7        17     35      7
     36      8     37        17     36     18
     41     37     44        18     37     19
     38     14     33        21     38     33
     39     32      4         5     39      4
     40     12     38        20     40     21
     42     10     11        37     41     19
     42     11     40        19     42     20
     43     26     30         4     43      3
     44      8      9        10     44      9
   0.000000E+00   0.000000E+00      2.500000E-01   0.000000E+00
   5.000000E-01   0.000000E+00      7.500000E-01   0.000000E+00
   1.000000E+00   0.000000E+00      1.000000E+00   2.500000E-01
   1.000000E+00   5.000000E-01      1.000000E+00   7.500000E-01
   1.000000E+00   1.000000E+00      7.500000E-01   1.000000E+00
   5.000000E-01   1.000000E+00      2.500000E-01   1.000000E+00
   0.000000E+00   1.000000E+00      0.000000E+00   7.500000E-01
   0.000000E+00   5.000000E-01      0.000000E+00   2.500000E-01
   7.500000E-01   5.000000E-01      7.165064E-01   6.250000E-01
   6.250000E-01   7.165064E-01      5.000000E-01   7.500000E-01
   3.750000E-01   7.165064E-01      2.834937E-01   6.250000E-01
   2.500000E-01   5.000001E-01      2.834936E-01   3.750001E-01
   3.749999E-01   2.834937E-01      4.999999E-01   2.500000E-01
   6.249999E-01   2.834936E-01      7.165062E-01   3.749999E-01
   2.077567E-01   2.075703E-01      4.017089E-01   1.482612E-01
   1.483122E-01   4.021604E-01      7.024530E-01   2.009149E-01
   2.021900E-01   7.009841E-01      1.471843E-01   5.790281E-01
   8.609530E-01   3.646002E-01      8.519265E-01   6.222693E-01
   7.943417E-01   7.355210E-01      2.050139E-01   8.383228E-01
   8.379703E-01   1.980158E-01      4.015485E-01   8.609668E-01
   7.210978E-01   8.353896E-01      5.831956E-01   8.603261E-01
   5.796704E-01   1.470356E-01      8.526899E-01   8.640727E-01
       0       0       0       0       0       0       0       0       0       0
       0       0       0       0       0       0       0       0       0       0
       0       0       0       0       0       0       0       0       0       0
       0       0       0       0       0       0       0       0       0       0
       0       0       0       0       0       0       0       0       0       0
       0       0       0       0       0       0       0       0       0       0
       0       0       0       0       0       0       0       0       0       0
       0       0       0       0       0       0       0       0       0       0
       0       0       0       0       0       0       0       0       0       0
       0       0       0       0       0       0       0       0       0       0
       0       0       0       0

Une aile NACA0012

Introduction

La paramétrisation du naca0012 est pour

Pour normaliser le naca0012, on utilise le changement de variables:

Avec le petit programme fortran suivant:

      PROGRAM NACA12
      DOUBLE PRECISION X,Y,C,XX,YY
      PARAMETER (C=1.008930411365D0)
      INTEGER NB,I
C --------  ENTRER LE NOMBRE D'INTERVALLES -----------
      READ *,NB
      DO I=0,NB
C ------------------------------------------------
C      POUR QUE LA PROGRESSION DES POINTS GENERE  |
C        SOIT PLUS DENSE AU VOISINAGE DE 0        |
C                4                                |
C      ON PREND X  PLUTOT QUE X                   |
C -------------------------------------------------
       XX= (DBLE(I)/DBLE(NB))**4
       X = C*XX
       Y = 5*.12*( 0.2969*SQRT(X) -0.126*X
     +            -0.3516*X**2    +0.2843*X**3 -0.1015*X**4)
       YY = Y / C
       PRINT*, XX,YY
      ENDDO
      END

on crée le fichier "naca12.21pts" en l'utilisant avec 20 intervalles.

le fichier "naca12.21pts" :

 0.0000000000000000 0.0000000000000000
 6.2500000000000000E-06 4.4290213310617220E-04
 1.0000000000000000E-04 1.7659364374154260E-03
 5.0625000000000000E-04 3.9520447457293480E-03
 1.6000000000000000E-03 6.9724900930098020E-03
 3.9062500000000000E-03 1.0785816129632120E-02
 8.1000000000000000E-03 1.5335254409279880E-02
 1.5006250000000000E-02 2.0543538719326360E-02
 2.5600000000000000E-02 2.6304013971036800E-02
 4.1006250000000000E-02 3.2467169756843240E-02
 6.2500000000000000E-02 3.8822480697810590E-02
 9.1506250000000000E-02 4.5076893778196040E-02
 0.1296000000000000 5.0833559338836810E-02
 0.1785062500000000 5.5577241169031020E-02
 0.2401000000000000 5.8675384962268570E-02
 0.3164062500000000 5.9403968668221290E-02
 0.4096000000000000 5.7000809116986590E-02
 0.5220062500000000 5.0728415111155510E-02
 0.6561000000000000 3.9881091560527040E-02
 0.8145062500000000 2.3576580331835940E-02
 1.000000000000000 2.0060198232592660E-14

La construction de la géométrie

Cliquer dans l'item SPLINE du menu construction
Taper < au clavier ou cliquer dans l'item Fi_XY du menu désignation puis entrer le nom du fichier contenant les points du NACA0012 en tapant naca12.21pts(CR), la spline du naca s'affiche à l'écran
Cliquer dans l'item VOIR_TOUT pour mieux voir la spline
On veut mailler l'extérieur du naca0012 pour cela il faut définir la frontière correspondant à l'infini que l'on placera à -5 en amont, +11 en aval et 8 orthogonalement. On construit les 3 points (-5,0), (0,8), (11,0)
- Cliquer dans l'item POINT du menu construction (en haut)
- Entrer au clavier x-5=0=0=8=11=0=
- Cliquer dans l'item VOIR_TOUT (à droite)
On voit les 3 points. On va construire l'arc passant par ces 3 points:
- Cliquer dans l'item ARC du menu construction (en haut)
- Cliquer dans l'item POINT du menu désignation (en bas)
- Cliquer successivement près des points (11,0), (0,8), (-5,0)
Pour finir le domaine il suffit de construire les 2 segments allant de "l'infini" au naca.
- Cliquer dans l'item SEGMENT du menu construction (en haut)
- Cliquer près du point (-5,0)
- Cliquer dans l'item PRECEDENTE du menu gestion de l'écran (à droite). On utilise le contexte d'affichage précédant qui met le naca plein écran.
- Cliquer près du point (0,0) (en dessous à gauche pour être sûr de désigner le premier point).
- Cliquer dans l'item SUIVANTE du menu gestion de l'écran (à droite), on voit tous les objets.
- Cliquer près du point (11,0)
Figure 1.32: 8.125

Toute la géométrie du naca0012 et de l'infini

Figure 1.33: 8.125

Zoom autour du naca0012

On a fini théoriquement la construction du domaine géométrique (figures 5 et 5).

Discrétisation des contours

On passe dans l'application PREP_MESH pour définir la discrétisation des contours.

Cliquer dans l'item PREP_MESH
La discrétisation des lignes se fait comme suit :
- Cliquer dans l'item NB_INTERVAL En vue de définir le nombre d'intervalles
- Taper 26=
- Cliquer près de l'arc (attention, on est en mode de désignation QLCONQ donc il ne faut pas cliquer près du point (8,0)). L'arc est discrétisé en 26 intervalles.
- Taper 39=
- Cliquer près du segment de gauche.
- Taper 30=.
- Cliquer près du segment de droite.
- Taper 50=.
- Cliquer dans l'item SPLINE du menu désignation .
- Cliquer dans la fenêtre graphique (il n'y a qu'une spline donc on la désignera).
- Cliquer dans l'item RAISON pour changer la répartion de point sur les lignes
- Taper 1.1=
- Cliquer près de la spline.
- Cliquer dans l'item SEGMENT du menu désignation .
- Cliquer près du segment de droite
- Taper 1/1.2=
- Cliquer près du segment de gauche
  
  Figure 1.34: 8.125
  
  Vu de la discrétisation de tous les contours
  
  Figure 1.35: 8.125
  
  Zoom de la discrétisation des contours autour du naca

On a fini la définition de la dicrétisation des lignes (figures 2 et 2).

Réalisation et édition du maillage

On crée le maillage en passant dans l'application EDIT_MESH.

Cliquer dans l'item EDIT_MESH puis taper quatre fois (CR) pour répondre par défaut aux 4 questions. Un maillage s'affiche (figure 1).

Figure 1.36: 8.125

Maillage du demi naca0012
On définit des références en vue de la prise en compte des conditions aux limites : 1 sur l'infini , 2 sur le NACA, 3 au bord de fuite.
Cliquer dans l'item MODIF_REF du menu l'application
- Taper 1= la référence courante est 1
- Cliquer dans l'item SOMMET du menu désignation pour passer en mode de désignation de sommet
- Cliquer dans la maillage près d'un sommet du maillage qui est sur l'arc et non une extrémité puis cliquer près des 2 extremités toujours dans le maillage. Tous les sommets de l'arc ont pour référence 1.
- Taper 2= la référence courante est 2
- On fait un zoom autour du naca
  - Cliquer dans l'item ZOOM + du menu gestion de l'écran
  - Cliquer dans l'item P_TABL du menu désignation pour passer en désignation de point curseur
  - Cliquer 2 fois en désignant les 2 coins opposés dans un rectangle fictif contenant le naca. Si le zoom n'est pas assez fort recliquer dans l'item ZOOM + puis entre 2 autres coins opposés Pour avoir un zoom autour de début du naca (coordonnées (0,0)):
    - Cliquer dans l'item C_MASQUE puis taper x0=0= pour centrer la fenêtre graphique autour du point (0,0) (Le x est le raccourci de POINT_XY)
    - Cliquer dans l'item ECHELLE puis taper 1000= pour avoir une échelle de 1000
    (figure 2)
    
    Figure 1.37: 8.125
    
    Zoom de 1000 du Maillage du demi naca0012
- Cliquer dans l'item SOMMET du menu désignation pour repasser mode de désignation de sommet
- En cas d'erreur (l'item MODIF_REF n'est plus marqué) pour continuer, il suffit de recliquer dans celui-ci et de taper 2= pour redéfinir la référence courante.
- Cliquer dans le maillage près d'un sommet du maillage qui est sur le naca et non une extremité, puis cliquer dans le maillage près de l'extrémité aval (de coordonnées (0,0)).
- Cliquer dans l'item C_MASQUE puis taper x0=1= pour centrer la fenêtre graphique autour du point (0,1) qui est le bord de fuite.
- Taper 3= la référence courante est 3
- Cliquer dans le maillage près du bord de fuite.
On a fini l'entrée des références. On sauve le maillage du demi naca.
Cliquer dans l'item SAUVER , taper mesh(CR) puis naca12_5(CR). Le fichier
[4] naca12_5.mesh est créé.
on symétrise le maillage pour obtenir un maillage complet, on commence par changer le point de vue pour voir tout le naca.
Cliquer dans l'item C_MASQUE puis taper x1/2=0= pour centrer la fenêtre graphique autour du point (,0)
Cliquer dans l'item ECHELLE puis taper 10= pour avoir une échelle de 10
Cliquer dans l'item SYMETRIE
- Cliquer dans l'item SEGMENT du menu désignation
- Cliquer près d'un des 2 segments
- Cliquer dans l'item TOUT , tout le maillage est symétrisé, la partie symétrisée est affichée en traits pointillés (figure 6).
  
  Figure 1.38: 8.125
  
  Zoom de 10 du maillage du naca0012 symétrisé
- Cliquer dans l'item EVALUER pour réaliser effectivement le maillage total (figure 6).
  
  Figure 1.39: 8.125
  
  Zoom de 10 du maillage du naca0012 complet
On sauve le maillage du naca.
Cliquer dans l'item SAUVER , taper mesh(CR), taper naca12(CR) le fichier
[4] naca12.mesh est créé

Suiv.: 1.5.7 Limitations et Bogues Sup.: Le préprocesseur EMC² Préc.: L'application EDIT_MESH Index Table des matières