1. Mode d’emploi#

Il y a quatre étapes à suivre pour produire un maillage avec le superviseur GMSH:

  1. Définition de la géométrie;

  2. Définition des discrétisations;

  3. Création du maillage GMSH et des GROUP_MA et objets «Physical» associés;

  4. Importation du maillage GMSH dans Aster.

Exemple simple d’utilisation:

Dans l’exemple suivant, on utilise les fonctionnalités du superviseur pour générer le maillage d’une plaque rectangulaire:

Géométrie

from Utilitai.sup_gmsh import *

larg = 5.

H_beton = 3.

H_S1 = 4.

t_beton = 25.

prog_S1 = 1.1

On importe le module et on définit quelques paramètres.

# Geometrie

O = Point(0 , 0 )

A = Point(larg, 0 )

B = Point(larg, H_beton)

C = Point(0 , H_beton)

D = Point(0 , -H_S1 )

E = Point(larg, -H_S1 )

OA = Line(O,A)

AB = Line(A,B)

BC = Line(B,C)

OC = Line(O,C)

OD = Line(O,D)

DE = Line(D,E)

AE = Line(A,E)

S2 = Surface(OA,AB,BC,OC)

S1 = Surface(OD,DE,AE,OA)

On crée des points, des lignes entre les points et des surfaces à partir des lignes.

# Discretisation

OA.Transfinite(1)

BC.Transfinite(1)

DE.Transfinite(1)

N_beton = int(H_beton/t_beton + 0.5)

AB.Transfinite(N_beton)

OC.Transfinite(N_beton)

N_S1 = Progress(H_S1, r=prog_S1, h=t_beton)

OD.Transfinite(N_S1,prog_S1)

AE.Transfinite(N_S1,prog_S1)

S2.Transfinite()

S1.Transfinite()

On définit la discrétisation des lignes et des surfaces.

# Maillage

mesh = Mesh()

mesh.Physical(“FOND”,DE)

mesh.Physical(“LAT_G”,OC,OD)

mesh.Physical(“LAT_D”,AB,AE)

mesh.Physical(“INTERFAC”,OA)

mesh.Physical(“HAUT”,BC)

mesh.Physical(“S2”,S2)

mesh.Physical(“S1”,S1)

On crée l’objet maillage et on définit les groupes de mailles qui seront des GROUP_MA dans la SD maillage Aster et des « Physical» dans GMSH (ces derniers seront nommés \(\mathit{GM1}\) , \(\mathit{GM2}\) , etc…).

MA = mesh.LIRE_GMSH(

MODI_QUAD = “OUI”

)

Importation du maillage dans Aster: \(\mathit{MA}\) est un maillage Aster.