Suiv.: Test 5: Problème d'élasticité
Sup.: 5 Tests batch
Préc.: 5 Tests batch
Table des matières
Calculer la température pour un domaine bidimensionnel.
Figure: Le problème thermique bidimensionnel
Trouver u(x) tel que:
avec
Ce problème admet une solution exacte:
On utilise l'élément fini TRIA 2P3D (triangle P3 Hermite)
Etape 1: Maillage
Exécuter le programme: | D4NOXX | |
Modules utilisés: | APNOPO | ... mailleur 2D [Guide Modulef - 3] |
TRIHER | ... méthode de Voronoï | |
REGMA2 | ... régularisation | |
GIBBS | ... renumérotation | |
Bibliothèques: | NOP2, NOPO, UTSD, UTII | |
Fichier de données: | D4NOPO.D | |
Structure de données de sortie: | T4NOPO |
Etape 2: Interpolation [Guide Modulef - 4]
Exécuter le programme: | D4COXX |
Module utilisé: | COMACO |
Bibliothèques: | THER, COSD, UTSD, UTII |
Elément fini: | TRIA 2P3D |
Structures de données de sortie: | T4MAIL |
T4COOR |
Etape 3: Description des données [Guide Modulef - 4]
Appeler le préprocesseur: | FOMIXX | |
Fichier de données: | D4FOMI.D | ... de la bibliothèque TESD |
Structures de données de sortie: | T4FORC | |
T4MILI |
Etape 4: Calcul et assemblage des tableaux élémentaires, description des conditions aux limites et résolution
Exécuter le programme: | D4ARXX | |
Modules utilisés | CONDL1 | ... construction du pointeur sur les degrés de liberté associés à chaque noeud [Guide Modulef - 4] |
COBDC1 | ... description des conditions aux limites [Guide Modulef - 4] | |
PREPAC | ... allocation des pointeurs | |
THEASS | ... calcul et assemblage des tableaux élémentaires. | |
CLIMPC | ... introduction des conditions aux limites | |
CHOLPC | ... factorisation de Cholesky | |
DRCHPC | ... substitutions | |
Bibliothèques: | COSD, THER, RESB, RESO, RESR, UTSD, UTII | |
Structure de données de sortie: | T4B |
Remarques: L'élément fini TRIA 2P3D utilise le fichier en accès direct POBA. Le nom de ce fichier est POBA dans le test. Pour changer ce nom, il est nécessaire de modifier la variable NOMFIC du programme D4ARXX. Le fichier POBA est construit, une fois, quand la bande est chargée 9.6, par le programme POBAXX de PPAL.
La fonction utilisateur VALCLR permet de spécifier des
conditions aux limites imposées non constantes. Cette fonction se
trouve à la fin du programme D4ARXX dans le fichier
D4ARXX.F de la bibliothèque TEST.
Le nombre de mots nécessaires pour stocker la matrice est 5903.
Le nombre de mots nécessaires pour stocker le second membre est 365.
Le super-tableau M est dimensionné à 45 000 mots dans tous les programmes.
Dimension du | Temps | |
Etape | super-tableau LM | (DPS8 Multics) |
Maillage | 6 000 | |
Interpolation | 4 600 | 3.38 |
Préparation | 5 000 | 52.18 |
Solution | 10 000 | 8.54 |
La S.D. B, contenant la température et ses dérivées par rapport à x et y en tout noeud du maillage, est dans le fichier S4B.D de la bibliothèque TESD.
La suite de cette section donne l'impression totale ou partielle des diverses S.D. intervenant dans les différents programmes.
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& IMPRESSION DE LA S.D. NOPO DE NIVEAU 1 &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& TITRE : BENCHMARK MODULEF NUMERO 4 :TEST DE THERMIQUE DATE ET NOM UTILISATEUR : 05/04/91 george TYPE DE LA STRUCTURE DE DONNEES : NOPO NIVEAU ET NUMERO D'ETAT : 1 0 NOMBRE DE TABLEAUX ASSOCIES : 1 LE TABLEAU 1 : CF01 DE TYPE CARACTER A 10 MOTS CONTENU DE CE TABLEAU : COURBE ASSOCIEE : 1 NOM : CF01 TYPE : 4 VALEURS : 1 8COURBE01(X,Y)=X**2+(Y-1)**2-.16; TABLEAU N O P 2 ---------------- CARACTERISTIQUES DU MAILLAGE : DIMENSION DE L'ESPACE (NDIM ) : 2 MAXIMUM DES NUMEROS DE REFERENCE (NDSR ) : 2 MAXIMUM DES NUMEROS DE SOUS-DOMAINE (NDSD ) : 1 NOEUDS ET POINTS COINCIDENT PARTOUT (NCOPNP) : 1 NOMBRE D'ELEMENTS DU MAILLAGE (NE ) : 129 NOMBRE DE TRIANGLES (NTRI ) : 129 NOMBRE D'ELEMENTS FRONTALIERS (NEF ) : 57 NOMBRE DE NOEUDS (NOE ) : 81 NOMBRE DE NOEUDS PAR SEGMENT (HORS EXTREMITES) : 0 TYPE DES VALEURS DES COORDONNEES (NTYCOO) : REEL1MOT DIFFERENCE MAX + 1 ENTRE 2 NOEUDS D'UN ELEMENT : 13 NOMBRE D'ELEMENTS GROSSIERS (NBEGM ) : 0 NOMBRE DE MOTS DU TABLEAU NOP5 (LNOP5 ) : 1224 AXES DE REFERENCE X,Y,Z (NTACOO) : 1 TABLEAU N O P 4 ---------------- COORDONNEES DES POINTS ------------------------------------------------------------------------------- | POINT | X | Y | | POINT | X | Y | ------------------------------------------------------------------------------- | 1 | -1.00000 | 1.00000 | | 2 | -.880000 | 1.00000 | | 3 | -.760000 | 1.00000 | | 4 | -.927111 | 0.920056 | | 5 | -1.00000 | 0.875000 | | 6 | -.701447 | 0.913023 | | 7 | -.640000 | 1.00000 | | 8 | -1.00000 | 0.750000 | | 9 | -.808490 | 0.872889 | | 10 | -.906893 | 0.796814 | --------------------------------------- --------------------------------------- TABLEAU N O P 5 ---------------- DESCRIPTIF DES ELEMENTS | ELEMENT : 1 | TRIANGLE | S-D : 1 | NOMBRE DE MOTS ASSOCIES : 0 NOMBRE DE NOEUDS : 3 : 79 76 75 | ELEMENT : 2 | TRIANGLE | S-D : 1 | NOMBRE DE MOTS ASSOCIES : 4 NOMBRE DE NOEUDS : 3 : 79 75 78 REFERENCE DES SOMMETS : 0 0 2 | ELEMENT : 3 | TRIANGLE | S-D : 1 | NOMBRE DE MOTS ASSOCIES : 0 NOMBRE DE NOEUDS : 3 : 11 12 6 . . . &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& IMPRESSION DE LA S.D. MAIL DE NIVEAU 1 &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& TITRE : DATE ET NOM UTILISATEUR : 05/04/91 george TYPE DE LA STRUCTURE DE DONNEES : MAIL NIVEAU ET NUMERO D'ETAT : 1 0 NOMBRE DE TABLEAUX ASSOCIES : 0 TABLEAU M A I 2 ---------------- DIMENSION DE L'ESPACE (NDIM) : 2 NOMBRE D'INCONNUES DE LA FORMULATION VARIATIONNELLE (NINCFV) : 1 NOMBRE DE TYPES D'ELEMENTS (NTYELM) : 1 NOMBRE DE NUMEROS DE REFERENCES (NNR) : 2 NOMBRE DE TYPES DE NOEUDS (NTYNOE) : 1 NOMBRE DE TYPES DE POINTS (NTYPOI) : 1 NOMBRE DE MOTS DU TABLEAU MAI3 (LMAI3) : 36 LES POINTS NE SONT DEFINIS QUE PAR LEURS COORDONNEES NOEUDS ET POINTS COINCIDENT PARTOUT (NCOPNP) : 1 NOMBRE MAXIMUM DE MOTS POUR UN ELEMENT (NMMAEL) : 18 NOMBRE DE MOTS DU TABLEAU MAIL (LMAIL) : 1482 OPTION DE STOCKAGE (NOPFI) : 0 TABLEAU M A I 3 ---------------- NOMBRE D'ELEMENTS (NE) : 129 NOMBRE DE NOEUDS (NOE) : 81 NOMBRE DE POINTS (NP) : 81 NOMBRE DE SOUS-DOMAINES (NDSD) : 1 *** TABLEAU DES TYPES D'ELEMENTS *** -------------------------------------------------------------------------------- | S-D | NOEUD | SEGMENT | TRIANGLE | QUADRANGLE | -------------------------------------------------------------------------------- | | DROIT COURBE | DROIT COURBE | DROIT COURBE | DROIT COURBE | -------------------------------------------------------------------------------- | 1 | | | | |TRIA2P3D| | | | -------------------------------------------------------------------------------- -------------------------------------------------------------------------------- | S-D | TETRAEDRE | PENTAEDRE | HEXAEDRE | SUPER-ELEMENT | -------------------------------------------------------------------------------- | | DROIT COURBE | DROIT COURBE | DROIT COURBE | DROIT COURBE | -------------------------------------------------------------------------------- | 1 | | | | | | | | | -------------------------------------------------------------------------------- TABLEAUX M A I 5 A M A I 8 ------------------------------- DESCRIPTION DES TYPES DE NOEUDS NOMBRE D'INCONNUES VARIATIONNELLES DE CHAQUE TYPE DE NOEUD ---------------------------------------------------------- TYPE 1 : 1 TYPE DE NOEUD 1 : -------------------- TEMPERATURE (INCONNUE : 1) : VN TEMPERATURE (INCONNUE : 1) : DX TEMPERATURE (INCONNUE : 1) : DY TABLEAUX M A I B A M A I K ------------------------------- DESCRIPTION DES TYPES D'ELEMENTS TYPE D'ELEMENT 1 : ------------------- NOM DE L'ELEMENT : THERTRIA2P3D CODE DE L'ELEMENT : 100004 L'ELEMENT EST UN : TRIANGLE NOMBRE DE NOEUDS : 3 NOMBRE DE POINTS : 3 NOMBRE DE TABLEAUX ASSOCIES : 1 ------------------------------------------------------------------------------- | NOEUD | TYPE | COORDONNEES SUR L'ELEMENT UNITE | ------------------------------------------------------------------------------- | 1 | 1 | 0.0000000 0.0000000 | 2 | 1 | 1.000000 0.0000000 | 3 | 1 | 0.0000000 1.000000 ------------------------------------------------------------------------------- ------------------------------------------------------------------------------- | POINT | TYPE | COORDONNEES SUR L'ELEMENT UNITE | ------------------------------------------------------------------------------- | 1 | 1 | 0.0000000 0.0000000 | 2 | 1 | 1.000000 0.0000000 | 3 | 1 | 0.0000000 1.000000 ------------------------------------------------------------------------------- TABLEAU 1 ( TYPE ENTIER ) : 1 INDICE(S) DE VALEUR MAXI : 7 TABLEAU M A I L ---------------- DESCRIPTION DES ELEMENTS SI NOPTFI = 0 ------------------------------------------------------------------- | ELEMENT | TYPE | NBRE DE MOTS | S-D | TABLEAUX | NBRE DE NOEUDS | ------------------------------------------------------------------- | 1 | 1 | 0 | 1 | 0 | 3 | NUMERO DES NOEUDS : 79 76 75 ------------------------------------------------------------------- | 2 | 1 | 4 | 1 | 1 | 3 | NUMERO DES NOEUDS : 79 75 78 REFERENCE (SOMMETS) : 0 0 2 ------------------------------------------------------------------- | 3 | 1 | 0 | 1 | 0 | 3 | NUMERO DES NOEUDS : 11 12 6 ------------------------------------------------------------------- . . . &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& IMPRESSION DE LA S.D. COOR DE NIVEAU 1 &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& TITRE : DATE ET NOM UTILISATEUR : 05/04/91 george TYPE DE LA STRUCTURE DE DONNEES : COOR NIVEAU ET NUMERO D'ETAT : 1 0 NOMBRE DE TABLEAUX ASSOCIES : 0 TABLEAU C O O 2 --------------- TYPE DU TABLEAU COO4 (NTYT) : 2 NOMBRE DE SES INDICES (NINDI) : 2 DIMENSION DU DOMAINE (NDIM) : 2 VALEUR MAXIMALE DU DEUXIEME INDICE (M2) : 81 CODE DE LA SEGMENTATION (NCODS) : 1 NOMBRE DE BLOCS (NBLOC) : 1 TYPE DES AXES DES COORDONNEES (NTACOO) : 1 TABLEAU C O O 3 --------------- LISTE DES NUMERO DES COLONNES DE FIN DES PAGES ---------------------------------------------- PAGE 1 : 81 PAGE ---------------------------------------------- TABLEAU C O O 4 --------------- ------------------------------------------------------------------------------- | POINT | X | Y | | POINT | X | Y | ------------------------------------------------------------------------------- | 1 | -1.00000 | 1.00000 | | 2 | -.880000 | 1.00000 | | 3 | -.760000 | 1.00000 | | 4 | -.927111 | 0.920056 | | 5 | -1.00000 | 0.875000 | | 6 | -.701447 | 0.913023 | | 7 | -.640000 | 1.00000 | | 8 | -1.00000 | 0.750000 | | 9 | -.808490 | 0.872889 | | 10 | -.906893 | 0.796814 | --------------------------------------- --------------------------------------- . . &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& IMPRESSION DE LA S.D. FORC DE NIVEAU 1 &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& TITRE : DATE ET NOM UTILISATEUR : 31/01/94 pate TYPE DE LA STRUCTURE DE DONNEES : FORC NIVEAU ET NUMERO D'ETAT : 1 0 NOMBRE DE TABLEAUX ASSOCIES : 1 LE TABLEAU 1 : FOO1 DE TYPE REEL1MOT A 10 MOTS CONTENU DE CE TABLEAU : F OMEGA ET F GAMMA NOM : FOO1 TYPE : 2 VALEURS : -0.4000000E+01 -0.4000000E+01 -0.4000000E+01 -0.4000000E+01 -0.4000000E+01 -0.4000000E+01 -0.4000000E+01 -0.3200000E+00 -0.3200000E+00 -0.3200000E+00 TABLEAU F O R 2 ---------------- PRESENCE ( NOMBRE DE GROUPES ) OU ABSENCE ( 0 ) DE RENSE IGNEMENT -------------------------------------------------- | | | | | | SOUS-DOMAINES | ELEMENTS | | | | | -------------------------------------------------- | VOLUMES | 0 | 0 | | SURFACES | 1 | 0 | | LIGNES | 1 | 0 | | POINTS | 0 | 0 | -------------------------------------------------- TABLEAU F O R 4 ---------------- DESCRIPTION DES NUMEROS DE REFERENCE DES SURFACES NUMERO DE REFERENCE ( OU DE SOUS-DOMAINE ) : 1 TYPE DE L'ELEMENT ( > 0 : ISOTROPE, < 0 : ANISOTROPE ) : 1 TRAITEMENT PAR TABLEAU NOM DU TABLEAU : FOO1 ADRESSE DANS CELUI-CI : 1 TABLEAU F O R 5 ---------------- DESCRIPTION DES NUMEROS DE REFERENCE DES LIGNES NUMERO DE REFERENCE ( OU DE SOUS-DOMAINE ) : 1 TYPE DE L'ELEMENT ( > 0 : ISOTROPE, < 0 : ANISOTROPE ) : 1 TRAITEMENT PAR TABLEAU NOM DU TABLEAU : FOO1 ADRESSE DANS CELUI-CI : 8 &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& IMPRESSION DE LA S.D. MILI DE NIVEAU 1 &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& TITRE : DATE ET NOM UTILISATEUR : 31/01/94 pate TYPE DE LA STRUCTURE DE DONNEES : MILI NIVEAU ET NUMERO D'ETAT : 1 0 NOMBRE DE TABLEAUX ASSOCIES : 1 LE TABLEAU 1 : MII1 DE TYPE REEL1MOT A 2 MOTS CONTENU DE CE TABLEAU : K ET G NOM : MII1 TYPE : 2 VALEURS : 0.1000000E+01 0.3000000E+01 TABLEAU M I L 2 ---------------- PRESENCE ( NOMBRE DE GROUPES ) OU ABSENCE ( 0 ) DE RENSEIGNEMENT -------------------------------------------------- | | | | | | SOUS-DOMAINES | ELEMENTS | | | | | -------------------------------------------------- | VOLUMES | 0 | 0 | | SURFACES | 1 | 0 | | LIGNES | 1 | 0 | | POINTS | 0 | 0 | -------------------------------------------------- TABLEAU M I L 4 ---------------- DESCRIPTION DES NUMEROS DE REFERENCE DES SURFACES NUMERO DE REFERENCE ( OU DE SOUS-DOMAINE ) : 1 TYPE DE L'ELEMENT ( > 0 : ISOTROPE, < 0 : ANISOTROPE ) : 1 TRAITEMENT PAR TABLEAU NOM DU TABLEAU : MII1 ADRESSE DANS CELUI-CI : 1 TABLEAU M I L 5 ---------------- DESCRIPTION DES NUMEROS DE REFERENCE DES LIGNES NUMERO DE REFERENCE ( OU DE SOUS-DOMAINE ) : 1 TYPE DE L'ELEMENT ( > 0 : ISOTROPE, < 0 : ANISOTROPE ) : 1 TRAITEMENT PAR TABLEAU NOM DU TABLEAU : MII1 ADRESSE DANS CELUI-CI : 2 &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& IMPRESSION DE LA S.D. B DE NIVEAU 1 &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& TITRE : DATE ET NOM UTILISATEUR : 05/04/91 george TYPE DE LA STRUCTURE DE DONNEES : B NIVEAU ET NUMERO D'ETAT : 1 2 NOMBRE DE TABLEAUX ASSOCIES : 0 TABLEAU B 2 ------------ TYPE DU TABLEAU (NTYT) : 2 NOMBRE DE SES INDICES ET LEURS MAXI (NIND..) : 2 1 243 TRAITEMENT (1:PAGES DE MEME TAILLE,0:SINON) (NCOD) : 1 NOMBRE DE PAGES DU TABLEAU B4 (NBLOC) : 1 NOMBRE DE NOEUDS (NOE) : 81 NOMBRE DE TABLEAUX B4 DANS CETTE S.D. (NBBLOC) : 1 NOMBRE CONSTANT DE D.L. PAR NOEUD OU 0 (ND) : 3 NOMBRE DE D.L. OU LONGUEUR D'UNE PAGE DE B4 (NTDL) : 243 CODE DE STOCKAGE DE B (NCODSB) : -1 TABLEAU B 3 ------------ POINTEUR SUR LA DERNIERE LIGNE DE CHAQUE PAGE DE B 1 0 2 243 TABLEAU B4(NDSM,NTDL) NUMERO : 1 ------------------------------------- TABLEAU B4(NDSM,NTDL) PAGINE EN 1 OU PLUSIEURS PAGES ---------------------------------------------------- PAGE 1 ---------- ------------------------------------------------------ | NOEUD | DL | DL GLOBAL | LE (LES) CAS DE CHARGE | ------------------------------------------------------ | 1 | 1 | 1 | 1.00000 | 2 | 2 | -2.00000 | 3 | 3 | 0.000000 | 2 | 1 | 4 | 0.774400 | 2 | 5 | -1.76000 | 3 | 6 | 0.000000 | 3 | 1 | 7 | 0.577600 | 2 | 8 | -1.52000 | 3 | 9 | 0.000000 | 4 | 1 | 10 | 0.865922
Le dessin des isothermes est produit par le préprocesseur TRMCXX:
Appeler le préprocesseur: | TRMCXX |
Structures de données: | T4MAIL |
T4COOR | |
T4B |
Figure: Isothermes (dessin dû à TRMCXX)