Télécharger la plaquette Cast3M


fr en

Télécharger explorer.procedur

Retour à la liste

Numérotation des lignes : 

  1. * EXPLORER  PROCEDUR  PASCAL    22/04/14    21:15:02     11341          
  2. DEBPROC  EXPLORER ;
  3. *=======================================================================
  4. *
  5. * PROCEDURE DONT LE BUT EST D'EXPLORER LES RESULTATS ISSUS DE :
  6. * - PASAPAS 
  7. * - VIBR (ajout bp, 2016-02)
  8. *
  9. * SYNTAXE :
  10. * EXPLORER MATAB | (LMOC1) | (TOPT1) ;
  11. *                |  'CHAR' | 
  12. * avec :
  13. * -MATAB : table de resultats issues de PASAPAS ou de VIBR
  14. * -CHAR  : mot-cle si l'on souhaite visualiser le chargement
  15. * -LMOC1 : mots-cles definissant les actions a produire
  16. * -TOPT1 : table d'options 
  17. *
  18. *=======================================================================
  19.  
  20. MESS '==============   DEBUT DE LA PROCEDURE EXPLORER   ==============';
  21.  
  22. *----------------------------------------------------------------------*
  23. *          RECUP et VERIF DES DONNEES D'ENTREE                         *
  24. *----------------------------------------------------------------------*
  25.  
  26. ARGU MATAB*'TABLE' ;
  27. ARGU MOT1/'MOT' ;
  28. ARGU LMOC1/'LISTMOTS' ;
  29. ARGU TOPT1/'TABLE' ;
  30.  
  31. * quelques logiques :
  32. FLMOT1 = EXIS MOT1;
  33. FLMOC1 = EXIS LMOC1;
  34. FLOPT1 = EXIS TOPT1;
  35. FLTRACX1 = EGA (VALE 'TRAC') 'X';
  36. FLTRACX2 = EGA (VALE 'TRAC') 'OPEN';
  37. FLTRACX = FLTRACX1 'OU' FLTRACX2;
  38. * Le choix est interactif si pas de mot-cles et trace X
  39. FLCHOI = FLTRACX ET (NON FLMOC1);
  40.  
  41.  
  42. * Quel est le type de resultats a traiter?
  43. SI (EXIS MATAB 'SOUSTYPE');
  44.   TYPTAB = MOT MATAB . 'SOUSTYPE';
  45. SINON;
  46.   MESS '*  SOUSTYPE DE LA TABLE INCONNU :';
  47.   MESS '*  => ON SUPPOSE QU IL S AGIT D UNE TABLE PASAPAS';
  48.   TYPTAB = MOT 'PASAPAS';
  49. FINSI;
  50.  
  51.  
  52. *----------------------------------------------------------------------*
  53. *                   CAS D'UNE TABLE ISSUE DE VIBR                      *
  54. *----------------------------------------------------------------------*
  55.  
  56. SI (    (EGA TYPTAB 'BASE_MODALE') 
  57.      OU (EGA TYPTAB 'LIAISONS_STATIQUES')
  58.      OU (EGA TYPTAB 'BASE_DE_MODES') );
  59. *  on teste l'existence des arguments     
  60.    SI (NON FLMOC1);
  61.      LMOC1 = MOTS 'TABL' 'DEFO' 'DEF0' 'MAIL';
  62.    FINSI;
  63.    SI FLOPT1; 
  64.      POSTVIBR MATAB LMOC1 TOPT1;
  65.    SINON;
  66.      POSTVIBR MATAB LMOC1;
  67.    FINSI;
  68.    QUITTER EXPLORER;
  69. FINSI;
  70.  
  71.  
  72. * ON VERIFIE LE TYPE AVANT DE POURSUIVRE
  73. SI (NEG TYPTAB 'PASAPAS');
  74.   MESS '* TABLE DE SOUSTYPE PASAPAS, BASE_MODALE OU LIAISONS_STATIQUES';
  75.   MESS '* ATTENDUE EN ENTREE D EXPLORER !';
  76.   ERRE 21;
  77. FINSI;
  78.  
  79.  
  80. *----------------------------------------------------------------------*
  81. *                   CAS D'UNE TABLE ISSUE DE PASAPAS                   *
  82. *----------------------------------------------------------------------*
  83.  
  84. * Si MATAB ne contient pas de resultat ==> proposer option CHAR :      *
  85. IDEPL1   = EXIS MATAB 'DEPLACEMENTS' ;
  86. ISIGM1   = EXIS MATAB 'CONTRAINTES' ;
  87. ITHER1   = EXIS MATAB 'TEMPERATURES' ;
  88. SI IDEPL1 ;
  89.   IDEPL1   = (DIME MATAB . 'DEPLACEMENTS') '>' 1 ;
  90. FINSI ;
  91. SI ISIGM1 ;
  92.   ISIGM1   = (DIME MATAB . 'CONTRAINTES' ) '>' 1 ;
  93. FINSI ;
  94. SI ITHER1 ;
  95.   ITHER1   = (DIME MATAB . 'TEMPERATURES') '>' 1 ;
  96. FINSI ;
  97. IRESU1   = (IDEPL1 ET ISIGM1) 'OU' ITHER1 ;
  98. *
  99. MESG00   = ' ' ;
  100. SI (NON IRESU1) ;
  101.   SI (NON (EXIS MOT1)) ;
  102.     MOT1     = MOT 'CHAR' ;
  103.     MESG00   = ' PAS DE RESULTAT => option CHAR : ' ;
  104.   FINSI ;
  105. FINSI ;
  106.  
  107. ICHGMOD1 = exis MATAB 'MODELES' ;
  108. ICHGMAT1 = FAUX;
  109. 'SI' ('EXIS' MATAB 'CHARGEMENT');
  110.    ICHGMAT1 = exis (MATAB . 'CHARGEMENT') 'MATE' ;
  111. 'FINSI';
  112.  
  113. *----------------------------------------------------------------------*
  114. *              SOUS-CAS DU CHARGEMENT D'UNE TABLE PASAPAS              *
  115. *----------------------------------------------------------------------*
  116.  
  117. SI (EGA MOT1 'CHAR') ;
  118.  
  119. * Le choix est interactif si pas de table d'options et trace X
  120.   FLCHOI = FLTRACX ET (NON FLOPT1);
  121.  
  122.   MOD1     = MATAB . 'MODELE' ;
  123.   CGT0     = MATAB . 'CHARGEMENT' ;
  124.   MAIL1    = EXTR MOD1 'MAIL' ;
  125. *
  126. *  Boucle principale :
  127.   LCHAR1   = EXTR CGT0 'COMP' ;
  128.   REPE BPRIN1 ;
  129. *
  130. *   Choix du chargement a vualiser :
  131.     SI FLCHOI;
  132. *     en interactif
  133.       TRAC 'NCLK' MAIL1 ;
  134.       MESG0    = CHAI MESG00 'Choisissez le chargement a visualiser ';
  135.       MRET1    = MENU MESG0 LCHAR1 ;
  136.       MESG00   = ' ' ;
  137.       SI (EGA MRET1 'Quitter') 
  138.         QUITTER BPRIN1 ;
  139.       FINSI ;
  140.     SINON;
  141. *     par defaut
  142.       MRET1 = EXTR LCHAR1 1;
  143. *     par une table d option
  144.       SI FLOPT1; 
  145.         SI (EXIS TOPT1 'TYPE');
  146.            MRET1 = TOPT1 . 'TYPE';
  147.         FINSI;
  148.       FINSI;
  149.     FINSI;
  150.     CGI1     = EXTR CGT0 MRET1 ;
  151.  
  152. *   Petit travail si chargement est defini par 2 tables
  153.     OPTI 'ERRE' 'CONT' ;
  154.     EVI1     = EXTR CGI1 'EVOL' ;
  155.     OPTI 'ERRE' 'NORM' ;
  156.     SI (EGA (TYPE EVI1) 'ANNULE') ;
  157.       TX1 TY1  = EXTR CGI1 MRET1 'TABL' ;
  158.       LI1      = PROG ;
  159.       REPE BX1 (DIME TX1) ;
  160.         LI1     = LI1 ET (PROG TX1 . (&BX1 - 1)) ;
  161.       FIN BX1 ;
  162.     SINON ;
  163.       LI1      = EXTR EVI1 'ABSC' ;
  164.     FINSI ;
  165. *
  166. *   Traitement CHPOINT / MCHAML
  167.     I2       = 1 ;
  168.     TI2      = EXTR LI1 I2 ;
  169.     CH2      = TIRE CGI1 MRET1 TI2 ;
  170.     SI (EGA (TYPE CH2 ) 'MCHAML') ;
  171.       MM1      = MOD1 ;
  172.     SINON ;
  173.       SI (EGA (TYPE CH2 ) 'CHPOINT') ;
  174.         MM1      = MAIL1 ;
  175.       SINON ;
  176.         ERRE ' On attend un MCHAML ou un CHPOINT ' ;
  177.       FINSI ;
  178.     FINSI ;
  179.  
  180. *   composante
  181.     LCOMP1   = EXTR CH2  'COMP' ;
  182.     LCOMP2   = LCOMP1 ;
  183.     MCOMP1   = EXTR LCOMP1 1 ;
  184.     SI FLOPT1;
  185.       SI (EXIS TOPT1 'COMP');
  186.          MCOMP1 = TOPT1 . 'COMP';
  187.       FINSI;
  188.     FINSI;
  189.  
  190. *   Boucle sur les instants du chargement 
  191.     NLI1     = DIME LI1 ;
  192.     REPE BINST1 ;
  193. *
  194. *     Traitement de l'echelle d'isovaleur 
  195.       XMIN1    = 'MINI' CH2  'AVEC' ('MOTS' MCOMP1) ;
  196.       XMAX1    = 'MAXI' CH2  'AVEC' ('MOTS' MCOMP1) ;
  197.       REPE BLEG1 NLI1 ;
  198.         TI3      = EXTR LI1 &BLEG1 ;
  199.         CHI3     = TIRE CGI1 MRET1 TI3 ;
  200.         XMINI1   = 'MINI' CHI3 'AVEC' ('MOTS' MCOMP1) ;
  201.         XMAXI1   = 'MAXI' CHI3 'AVEC' ('MOTS' MCOMP1) ;
  202.         SI (XMINI1 '<' XMIN1) ; XMIN1 = XMINI1 ; FINSI ;
  203.         SI (XMAXI1 '>' XMAX1) ; XMAX1 = XMAXI1 ; FINSI ;
  204.       FIN BLEG1 ;
  205.       XINT1    = XMAX1 - XMIN1 ;
  206.       XPAS1    = XINT1 / 18. ;
  207.       SI (XPAS1 NEG 0.) ;
  208.         EE1      = ENTI 'INFE' (('LOG' XPAS1) / ('LOG' 10.)) ;
  209.         YPAS1    = ENTI 'INFE' (XPAS1 / (10. ** EE1)) ;
  210.         YPAS1    = YPAS1 * (10. ** EE1) ;
  211.         NN1      = ENTI 'INFE' (XMIN1 / YPAS1) ;
  212.         YMIN1    = NN1 * YPAS1 ;
  213.         NN1      = ENTI 'INFE' (XINT1 / YPAS1) ;
  214.         YMAX1    = YMIN1 + ((NN1 + 1) * YPAS1) ;
  215.         LISO1    = PROG YMIN1 'PAS' YPAS1 YMAX1 ;
  216.       SINON ;
  217.         LISO1    = TEXT ' ' ;
  218.       FINSI ;
  219.  
  220. *     Cas non-interactif :
  221.       SI (NON FLCHOI);
  222. *       on itere betement
  223.         TI2 = EXTR LI1 I2 ;
  224.         CH2 = TIRE CGI1 MRET1 TI2 ;
  225.         CH2 = 'EXCO' CH2  MCOMP1 MCOMP1 ;
  226.         MESG1  = CHAI ' Chargement ' MRET1 ' a l instant '
  227.                 (CHAI 'FORMAT' '(E10.2)' TI2) ;
  228.         TRAC CH2  MM1 'TITR' MESG1 LISO1 ;
  229. *       on itere betement
  230.         I2  = I2 + 1 ;
  231.         SI (I2 > NLI1); QUITTER BINST1; FINSI;
  232.         ITERER BINST1;
  233.       FINSI;
  234. *
  235. *     Choix des instants a vualiser 
  236.       TRAC 'NCLK' CH2  MM1 LISO1 ;
  237.  
  238.       MESG1    = CHAI ' Chargement ' MRET1 ' a l instant '
  239.                 (CHAI 'FORMAT' '(E10.2)' TI2) ;
  240.       MRET2    = MENU MESG1 ' Visualiser ' ' Composante '
  241.                    ' <-- ' ' --> ' ' No instant' ;
  242. *
  243. *     On quitte pour boucle principale
  244.       SI (EGA MRET2 'Quitter') ;
  245.         QUITTER BINST1 ;
  246.       FINSI ;
  247. *
  248. *     On visualise le champ a l'instant demande
  249.       SI (EGA MRET2 ' Visualiser ') ;
  250.         TRAC CH2  MM1 'TITR' MESG1 LISO1 ;
  251.       FINSI ;
  252. *
  253. *     Traitement choix instant suivant / precedent
  254.       SI (EGA MRET2 ' <-- ') ;
  255.         I2        = I2 - 1 ;
  256.         I2        = I2 - ((I2 - NLI1) / NLI1 * NLI1) ;
  257.       FINSI ;
  258.       SI (EGA MRET2 ' --> ') ;
  259.         I2        = I2 + 1 - (I2 / NLI1 * NLI1) ;
  260.       FINSI ;
  261. *
  262. *     Traitement choix No instant 
  263.       SI (EGA MRET2 ' No instant') ;
  264.         MESG2     = CHAI
  265.           ' Entrez le No du pas de temps a visualiser (1 a '
  266.           'FORMAT' '(I)' NLI1 ') : ' ;
  267.         REPE BSAIS1 ;
  268.           I2        = SAIS MESG2 'ENTIER' ;
  269.           SI ((I2 '<' 0) 'OU' (I2 '>' NLI1)) ;
  270.             MESG2     = CHAI ' Nombre inacceptable !' MESG2 ;
  271.             'ITER' BSAIS1 ;
  272.           FINSI ;
  273.           QUITTER BSAIS1 ;
  274.         FIN BSAIS1 ;
  275.       FINSI ;
  276. *
  277. *     Champ pour nouvel instant choisi
  278.       TI2     = EXTR LI1 I2 ;
  279.       CH2     = TIRE CGI1 MRET1 TI2 ;
  280.       CH2     = 'REDU' CH2  MM1 ;
  281.  
  282. *     Traitement des options d'affichage
  283. *      SI (EGA MRET2 ' Options ') ;
  284. *        TRAC 'NCLK' CH2  MM1 'TITR' MESG1 ;
  285. *        MRET3   = MENU ' Choisissez une option d affichage '
  286. *                  ' Composante ' ;
  287. *      FINSI ;
  288. *
  289. *     Traitement choix composante
  290.       SI (EGA MRET2 ' Composante ') ;
  291.         MRET3   = MENU ' Choisissez la composante a afficher '
  292.           (LCOMP1 ET ('MOTS' 'Tout')) ;
  293.         SI (EXIS LCOMP1 MRET3) ;
  294.           LCOMP2  = 'MOTS' MRET3 ;
  295.         FINSI ;
  296.         SI (EGA MRET3 'Tout') ;
  297.           LCOMP2  = LCOMP1 ;
  298.         FINSI ;
  299.       FINSI ;
  300.       SI (NEG LCOMP1 LCOMP2) ;
  301.         CH2     = 'EXCO' CH2  LCOMP2 LCOMP2 ;
  302.       FINSI ;
  303.       MCOMP1   = EXTR LCOMP2 1 ;
  304. *List Resu CH2  ;
  305. *
  306. *     Traitement choix vecteur
  307. *      SI (EGA MRET3 ' Vecteur ') ;
  308. *        TYP3     = EXTR CH2  'TYPE' ;
  309. *        SI (EGA (TYPE TYP3) 'LISTMOTS') ;
  310. *          MESG1     = CHAI
  311. *            ' On attend un CHPOINT !' ;
  312. *        FINSI ;
  313. *        List TYP3 ;
  314. *        SI ((EGA TYP3 'FORCES') 'OU' (EGA TYP3 'DEPLACEM')) ;
  315. *        SINON ;
  316. *          MESG1     = CHAI
  317. *            ' On attend un CHPOINT de FORCES ou DEPLACEMENTS !' ;
  318. *          TRAC 'NCLK' CH2  MM1 'TITR' MESG1 ;
  319. *        FINSI ;
  320. *      FINSI ;
  321. *
  322.     FIN BINST1 ;
  323.  
  324.     SI (NON FLCHOI); QUITTER BPRIN1; FINSI;
  325.  
  326.   FIN BPRIN1 ;
  327.   QUITTER EXPLORER ;
  328.  
  329. FINSI ;
  330.  
  331.  
  332. *----------------------------------------------------------------------*
  333. *              SOUS-CAS DES RESULTATS D'UNE TABLE PASAPAS              *
  334. *----------------------------------------------------------------------*
  335.  
  336. * INITIALISATIONS DES NOMS DES COMPOSANTES DES EFFORTS
  337. leforces = TEXT 'FX FY FZ';
  338. SI (EGA (VALE DIME) 2); 
  339.    leforces = TEXT 'FX FY';
  340. FINSI;
  341. SI (EGA (VALE 'MODE') 'AXIS');
  342.    leforces = TEXT 'FR FZ';
  343. FINSI;
  344. SI (EGA (VALE 'MODE') 'FOUR');
  345.    leforces = TEXT 'FR UT FZ';
  346. FINSI;
  347.  
  348. * INITIALISATIONS ------------------------------------------------------
  349.  
  350. * ON COMMENCE PAR TRACER LE MAILLAGE lorsqu'on est en interactif
  351. si ICHGMOD1 ;
  352.   MOD0     = MATAB . 'MODELES' . 0 ;
  353. sino ;
  354.   MOD0     = MATAB . 'MODELE' ;
  355. fins ;
  356. lemail   = EXTR MOD0 'MAIL' ;
  357. lemail22 = EXTR MOD0 'MAIL' 'FROTTEMENT';
  358. * enlever les elements de frottements car les multiplicateurs de lagrange de frottement ne suivent pas la
  359. * deformee
  360. lemail = lemail diff lemail22;
  361. * suite a la demande de CB
  362. SI (EXIS MATAB 'RIGIDITE_CONSTANTE');
  363.   SI (EGA (TYPE MATAB . 'RIGIDITE_CONSTANTE') 'RIGIDITE');
  364.     lemail = lemail et (EXTR MATAB . 'RIGIDITE_CONSTANTE' 'MAIL');
  365.   FINSI;
  366. FINSI;
  367. SI FLCHOI; TRAC 'NCLK' 'CACH' lemail; FINSI;
  368.  
  369. * ON VERIFIE L'EXISTENCE DES CHAMPS -> lopt 
  370. lopt = MOTS;
  371. SI (EXIS matab 'DEPLACEMENTS');
  372.    tdpl = MATAB . 'DEPLACEMENTS';
  373.    lopt = MOTS 'Depl';
  374. FINSI;
  375. SI ( EXIS matab 'CONTRAINTES');
  376.    lopt = lopt ET (MOTS 'Cont');
  377.    tctr = MATAB . 'CONTRAINTES';
  378.    modmec = EXTR  MOD0 'FORMULATION' 'MECANIQUE';
  379.    matmec = REDU MATAB . WTABLE . 'CARACTERISTIQUES' modmec;
  380. FINSI;
  381. SI ( EXIS matab 'VARIABLES_INTERNES');
  382.    tvar = MATAB . 'VARIABLES_INTERNES';
  383.    lopt = lopt ET (MOTS 'Var_');
  384. FINSI;
  385. SI ( EXIS matab 'REACTIONS');
  386.    trea = MATAB . 'REACTIONS';
  387.    lopt = lopt ET (MOTS 'Reac');
  388. FINSI;
  389. SI ( EXIS matab 'TEMPERATURES');
  390.    lopt = lopt ET (MOTS 'Temp');
  391.    modthe = EXTR  MOD0 'FORMULATION' 'THERMIQUE';
  392.    maithe = EXTR modthe 'MAILLAGE';
  393.    tctem = MATAB . 'TEMPERATURES';
  394. FINSI;
  395. lopt = lopt ET (MOTS 'Quitter');
  396.  
  397. * ON INDEXE LA TABLE DES TEMPS (DES PAS DE CALCUL)
  398. ttem = MATAB . 'TEMPS';
  399. tindex = INDEX ttem;
  400. ltindex = DIME tindex - 1;
  401. ind=0;d1=d1;
  402.  
  403. * VALEUR PAR DEFAUT DES LOGIQUES
  404. evolu = FAUX;
  405. defor = FAUX;
  406. ampdef = TEXT ' ';
  407. contr = FAUX;
  408. react = FAUX;
  409. tempe = FAUX;
  410. evo_esp = FAUX;
  411. evo_tem = FAUX;
  412. depla  = FAUX;
  413. var_int= FAUX;
  414. quit   = FAUX;
  415.  
  416. * En cas de modele de contact, valeur amplitude deformee = 1 par defaut
  417. si (exis MOD0 form contact) ;
  418.   ampdef = 1. ;
  419. fins ;
  420.  
  421. * BOUCLE PRINCIPALE ----------------------------------------------------
  422.  
  423. REPE bouc2;
  424.  
  425. SI  (EXIS matab 'DEPLACEMENTS');
  426.   SI FLCHOI;
  427. *   Par defaut :
  428.     si (&bouc2 ega 1) ;
  429.       defor = VRAI ;
  430.       contr = VRAI ;
  431.     fins ;
  432. *   Choix interactif :
  433.     SI ( EXIS matab 'TEMPERATURES');
  434.       defoP contP reacP evolP tempP = defor contr react evolu tempe ;
  435.       defor contr react evolu tempe quit= CHOI 'cochez les options désirées'
  436.       defor contr react evolu tempe quit;
  437.       si (contP et tempe) ; contr = faux ; fins ;
  438.       si (contr et tempP) ; contr = vrai ; fins ;
  439.       si (react et tempe) ; react = faux ; tempe = faux ; iter bouc2 ; fins ;
  440.     SINON ;
  441.       defor contr react evolu quit= CHOI 'cochez les options désirées'
  442.       defor contr react evolu quit;
  443.     FINSI ;
  444.   SINON;
  445.     SI FLMOC1;
  446. *     Choix par un listmots en argument
  447.       defor = EXIS LMOC1 'DEPL';     
  448.       contr = EXIS LMOC1 'CONT';
  449.       react = EXIS LMOC1 'REAC';
  450.       evolu = EXIS LMOC1 'EVOL'; 
  451.       var_int = EXIS LMOC1 'VAR';
  452.     SINON;
  453. *     Choix par defaut
  454.       defor = VRAI;
  455.       contr = VRAI;
  456.     FINSI;
  457.     SI FLOPT1;
  458.       SI (EXIS TOPT1 'AMPL');
  459.         ampdef = TOPT1 . 'AMPL';
  460.       FINSI;
  461.     FINSI;
  462.   FINSI;
  463. SINON;
  464.   SI ( EXIS matab 'TEMPERATURES');
  465.     SI FLCHOI;
  466.       tempe evolu quit = CHOI 'cochez les options désirées'
  467.       tempe evolu quit;
  468.     SINON;
  469.       SI FLMOC1;
  470.         tempe = EXIS LMOC1 'TEMP'; 
  471.       SINON;
  472.         tempe = VRAI;
  473.       FINSI;
  474.     FINSI;
  475.   SINON ;
  476.     MESS '* LA TABLE NE CONTIENT NI DEPLACEMENTS NI TEMPERATURES !';
  477.     QUITTER EXPLORER;    
  478.   FINSI; 
  479. FINSI;
  480.  
  481. * CHOIX = QUIT  -->  ON SORT
  482. SI quit ; QUITTER EXPLORER; FINSI;
  483. * 1 seule iteration si pas de choix interactif
  484. SI ((NON FLCHOI) ET (EGA &bouc2 2)); QUITTER EXPLORER; FINSI;
  485.  
  486. * AUCUN CHOIX = QUIT  -->  ON NE FAIT QUE REPROPOSER UN NOUVEAUX CHOIX
  487. SI ((NON defor) ET (NON contr) ET (NON react) 
  488.       ET ( NON evolu) ET (NON tempe));       
  489.    ITERER bouc2;
  490. FINSI;
  491.  
  492. * CHOIX = EVOLU
  493. SI evolu;
  494.  
  495.   SI FLCHOI;
  496. *   Choix interactif
  497.     ret = MENU '  Faites votre choix' 'evo_espace' 'evo_temps' ;
  498.     SI ( EGA ret 'Quitter' ) ; QUITTER EXPLORER; FINSI;
  499.     SI ( EGA ret 'evo_espace') ; 
  500.       evo_esp=VRAI; evo_tem=FAUX;
  501.     SINON;
  502.       evo_esp=FAUX;
  503.       SI ( EGA ret 'evo_temps' ); evo_tem=VRAI; FINSI;
  504.     FINSI;
  505.   SINON;
  506. *   Choix par une table d'option en argument
  507.     SI (NON FLOPT1);
  508.       MESS '* Une table d options est necessaire en argument ';
  509.       MESS '* dans le cas du trace d evolutions';
  510.       ERRE 21;
  511.     FINSI;
  512.     SI (NON (EXIS TOPT1 'EVOL'));
  513.       MESS '* Une table d options avec l indice EVOL est necessaire ';
  514.       MESS '* en argument dans le cas du trace d evolutions';
  515.       ERRE 21;
  516.     FINSI;
  517.     SI (NON (EXIS TOPT1 . 'EVOL' 'TYPE'));
  518.       MESS '* Une table d options avec l indice EVOL . TYPE est ' ;
  519.       MESS '* necessaire en argument dans le cas du trace d evolutions';
  520.       ERRE 21;
  521.     FINSI;
  522.     evo_esp = EGA  TOPT1 . 'EVOL' . 'TYPE'  'ESPA';
  523.     evo_tem = EGA  TOPT1 . 'EVOL' . 'TYPE'  'TEMP';
  524. *   logiques deja affectes + haut
  525.     depla = defor;
  526.     lret = mots ;
  527.     si depla; ret = mot 'Depl'; lret = lret et (mots ret); finsi;
  528.     si contr; ret = mot 'Cont'; lret = lret et (mots ret); finsi;
  529.     si react; ret = mot 'Reac'; lret = lret et (mots ret); finsi;
  530.     si tempe; ret = mot 'Temp'; lret = lret et (mots ret); finsi;
  531.     si var_int; ret = mot 'Var_'; lret = lret et (mots ret); finsi;
  532.     SI ((dime lret) NEG 1);
  533.        MESS '* Champ concerne par l evolution ambigu parmi :';
  534.        LIST lret;
  535.        MESS '* Veuillez fournir un LISTMOT plus explicite';
  536.        ERRE 21;
  537.     FINSI;
  538.   FINSI;
  539.  
  540. * EVOLUTION SPATIALE  
  541.   SI evo_esp;
  542. *    evolution spatiale sur quelle ligne? et de quelle grandeur?
  543.      SI FLCHOI;
  544. *       Choix interactif
  545.         depla=FAUX; contr=FAUX; react=FAUX; tempe=FAUX; var_int=faux;
  546.         List *maillage ligne ;
  547.         mai = SAIS ' Nom de la ligne suivant laquelle tracer (voir liste): '
  548.              'MAILLAGE' ;
  549.         ret = MENU 'Choisissez le champ concerne' lopt;
  550. *       MESS ' retour de menu ' ret;      
  551.         SI ( EGA ret 'Quitter' ); QUITTER EXPLORER; FINSI;    
  552.         SI ( EGA ret 'Depl') ; depla=VRAI; FINSI;
  553.         SI ( EGA ret 'Cont') ; contr=VRAI; FINSI;     
  554.         SI ( EGA ret 'Reac') ; react=VRAI; FINSI;     
  555.         SI ( EGA ret 'Temp') ; tempe=VRAI; FINSI;
  556.         SI ( EGA ret 'Var_') ; var_int=VRAI;FINSI;
  557.      SINON;
  558. *       Choix par une table d'option en argument
  559.         SI (EXIS TOPT1 . 'EVOL' 'LIGN');
  560.           mai=TOPT1 . 'EVOL' . 'LIGN';
  561.         SINON;
  562.           MESS '* Une table d options avec l indice EVOL . LIGN est ';
  563.           MESS '* necessaire en argument dans le cas du trace ';
  564.           MESS '* d evolutions spatiales';
  565.           ERRE 21;
  566.         FINSI;
  567.      FINSI;
  568.      SI ( EGA ret 'Depl') ; 
  569.        ICHPO=VRAI; mot1 = MOT 'DEPLACEMENTS'; 
  570.      FINSI;
  571.      SI ( EGA ret 'Cont') ; 
  572.        ICHPO=FAUX; mot1 = MOT 'CONTRAINTES'; 
  573.      FINSI;        
  574.      SI ( EGA ret 'Reac') ; 
  575.        ICHPO=VRAI;  mot1 = MOT 'REACTIONS';
  576.      FINSI;      
  577.      SI ( EGA ret 'Temp') ; 
  578.        ICHPO=VRAI;  mot1 = MOT 'TEMPERATURES';
  579.      FINSI;
  580.      SI ( EGA ret 'Var_') ; 
  581.        ICHPO=FAUX;  mot1 = MOT 'VARIABLES_INTERNES';
  582.      FINSI;                
  583.      lcomp= EXTR   (matab . mot1 . 1)  'COMP' ; 
  584.      SI FLCHOI;
  585. *      Choix interactif
  586.        compo = MENU  'Composantes?' lcomp;
  587.      SINON;
  588.        SI (EXIS TOPT1 . 'EVOL' 'COMP');
  589. *        Choix par une table d'option en argument
  590.          compo = TOPT1 . 'EVOL' . 'COMP';
  591.        SINON;
  592. *        Choix par defaut
  593.          compo = EXTR lcomp 1;
  594.        FINSI;
  595.      FINSI;                
  596.      cp= ttem . ind;   
  597.      SI ( EGA compo 'Quitter' ); QUITTER EXPLORER; FINSI;
  598. *    titre du trace
  599.      SI FLCHOI;
  600. *      Choix interactif
  601.        tit= SAIS 'Donner un titre pour le dessin (entrer quotes) :' 
  602.                  'MOT';
  603.      SINON;
  604.        SI (EXIS TOPT1 . 'EVOL' 'TITR');
  605. *        Choix par une table d'option en argument
  606.          tit = TOPT1 . 'EVOL' . 'TITR';
  607.        SINON;
  608. *        Choix par defaut
  609.          tit = CHAI 'evolution spatiale de' ' ' compo;
  610.        FINSI;
  611.      FINSI;
  612.      tcp = ttem . ind;
  613.      ch = CHAI tit ' pour le temps' ' ' tcp;
  614.      TITRE ch;
  615. *    trace initial
  616.      SI ( depla 'OU' tempe 'OU' react );
  617.        SI depla; mata=MATAB . 'DEPLACEMENTS';FINSI;
  618.        SI tempe; mata=MATAB . 'TEMPERATURES';FINSI;
  619.        SI react; mata=MATAB . 'REACTIONS'   ;FINSI;
  620.        tcp= ttem . ind;
  621.  
  622. *      Pour l'indice 0 les reactions sont 'VIDES' ==> Test sur l'existence de la composante
  623.        SI ('EXIS' mata . ind compo);
  624.          ev = EVOL 'CHPO'  mata . ind mai compo;
  625.        SINON;
  626.          ev ='VIDE' 'EVOLUTION';
  627.        FINS;
  628.        ichpo=vrai;     
  629. *        dess ev;
  630.      SINON ;
  631.        SI var_int; mata = MATAB . 'VARIABLES_INTERNES' ; FINSI;
  632.        SI contr  ; mata = MATAB . 'CONTRAINTES'        ; FINSI;
  633.        rdt = CHAN 'CHPO' mata . ind MOD0 ;
  634.        ichpo=FAUX;
  635.        ev =  EVOL 'CHPO'  rdt  mai compo;     
  636.      FINSI;   
  637. *    boucle second niveau sur les pas de temps ......     
  638.      REPE bouc77 ; 
  639.        SI FLCHOI;
  640. *        Choix interactif
  641.          dess ev;
  642.          cha= CHAI 'pas courant :' ind; 
  643.          ret= MENU cha 'options' 'n°' 'retracer' 'precedent' 'suivant'; 
  644.          SI (EGA ret 'options') ; QUITTER bouc77  ; FINSI;
  645.          SI (EGA ret 'precedent') ; ind = ind - 1; FINSI;
  646.          SI (EGA ret 'suivant')   ; ind = ind + 1; FINSI;
  647.          SI (EGA ret 'retracer')  ; ind = ind    ; FINSI; 
  648.          SI (EGA ret 'Quitter')   ; QUITTER bouc2; FINSI;
  649.          SI (EGA ret 'n°'); ind=SAIS 'Numéro du pas : ' 'ENTIER'; FINS;
  650.          SI (ind > ltindex); ind = 0      ; FINSI;
  651.          SI (ind < 0);       ind = ltindex; FINSI;
  652.        FINSI;        
  653.        tcp = ttem . ind;
  654.        ch = CHAI tit ' pour le temps' ' ' tcp;
  655.        TITRE ch; 
  656.        SI  ichpo;
  657. *         Pour l'indice 0 les reactions sont 'VIDES' ==> Test sur l'existence de la composante
  658.           SI ('EXIS' mata . ind compo);
  659.             ev = EVOL 'CHPO'  mata . ind mai compo;
  660.           SINON;
  661.             ev ='VIDE' 'EVOLUTION';
  662.           FINS;
  663.        SINON;
  664.           rdt= CHAN 'CHPO' mata . ind MOD0 ;
  665.           ev =  EVOL 'CHPO'  rdt  mai compo; 
  666.        FINSI;        
  667.        si (&bouc77 ega 1) ; dess ev ; fins ;
  668.        SI (NON FLCHOI);
  669. *        non interactif : on trace tous les pas de temps
  670.          ind = ind + 1;
  671.          SI (ind > ltindex); QUITTER bouc77; FINSI;
  672.        FINSI;               
  673.      FIN bouc77; 
  674. *    fin de la boucle second niveau sur les pas de temps ......     
  675.   FINSI;
  676.  
  677. * EVOLUTION TEMPORELLE  
  678.   SI evo_tem;
  679. *    evolution temporelle en quel point? et de quelle grandeur?
  680.      SI FLCHOI;
  681. *      Choix interactif
  682.        depla=FAUX;contr=FAUX; react=FAUX; tempe=FAUX;var_int=faux; 
  683. *        ret  = MENU 'Choisissez le champ concerne' 
  684. *               'depla' 'contr' 'react' 'tempe' 'var_int'; 
  685.        ret  = MENU 'Choisissez le champ concerne' lopt;
  686.        SI ( EGA ret 'Quitter' ); QUITTER explorer; FINSI;     
  687.      FINSI;
  688.      SI ( EGA ret 'Depl') ; ichpo=VRAI;
  689.            mot1 = MOT 'DEPLACEMENTS';FINSI;
  690.      SI ( EGA ret 'Cont') ; ichpo=FAUX;
  691.            mot1 = MOT 'CONTRAINTES'; FINSI;        
  692.      SI ( EGA ret 'Reac') ; ichpo=VRAI;
  693.            mot1 = MOT 'REACTIONS';FINSI;      
  694.      SI ( EGA ret 'Temp') ; ichpo=VRAI;
  695.            mot1 = MOT 'TEMPERATURES';FINSI;
  696.      SI ( EGA ret 'Var_') ; ichpo=FAUX;
  697.           mot1 = MOT 'VARIABLES_INTERNES';FINSI; 
  698.      lcomp= EXTR   matab . mot1 . 1 'COMP'; 
  699.      SI FLCHOI;
  700. *      Choix interactif
  701.        compo=  MENU  'Composantes?' lcomp;
  702.        SI ( EGA compo 'Quitter' ); QUITTER EXPLORER; FINSI;         
  703.        SI ichpo;
  704.           List *point ;
  705.           P1= SAIS 'Point concerne (voir liste)? : '  'POINT'; 
  706.        SINON;
  707.           N1=SAIS 'Numero de la zone concernee : '  'ENTIER'; 
  708.           N2=SAIS 'Numero de l element concerne : '  'ENTIER';
  709.           N3=SAIS 'Numero du point de gauss concerne : '  'ENTIER';
  710.        FINSI;
  711.      SINON;
  712.        SI (EXIS TOPT1 . 'EVOL' 'COMP');
  713. *        Choix par une table d'option en argument
  714.          compo = TOPT1 . 'EVOL' . 'COMP';
  715.        SINON;
  716. *        Choix par defaut
  717.          compo = EXTR lcomp 1;
  718.        FINSI;
  719.        SI ichpo;
  720.          SI (EXIS TOPT1 . 'EVOL' 'POIN'); p1 = TOPT1 . 'EVOL' . 'POIN'; 
  721.          SINON; MESS '* Le point concerne doit etre fourni a lindice' 
  722.                 'EVOL . POIN de la table d option';
  723.                 ERRE 21;
  724.          FINSI;
  725.        SINON;
  726.          SI (EXIS TOPT1 . 'EVOL' 'ZONE'); n1 = TOPT1 . 'EVOL' . 'ZONE'; 
  727.          SINON; n1 = 1;
  728.          FINSI;
  729.          SI (EXIS TOPT1 . 'EVOL' 'ELEM'); n2 = TOPT1 . 'EVOL' . 'ELEM'; 
  730.          SINON; n2 = 1;
  731.          FINSI;
  732.          SI (EXIS TOPT1 . 'EVOL' 'PTG' ); n3 = TOPT1 . 'EVOL' . 'PTG' ; 
  733.          SINON; n3 = 1;
  734.          FINSI;
  735.        FINSI;
  736.      FINSI;                
  737.      SI ichpo; ev1 = evol 'TEMP' matab mot1 compo p1;
  738.      SINON;    ev1 = evol 'TEMP' matab mot1 compo n1 n2 n3;
  739.      FINSI;
  740.      dess ev1;    
  741.   FINSI;   
  742.  
  743. * CHOIX = AUTRE QUE EVOLU
  744. SINON;
  745.  
  746. * boucle second niveau sur les pas de temps ......     
  747.   REPE bouc1;
  748.  
  749.     cha = CHAI 'pas courant :' ind;
  750.     SI FLCHOI;
  751. *     Choix interactif
  752.       ret = MENU cha 'options' 'n°' 'retracer' 'precedent' 'suivant';
  753.       SI (EGA ret 'options') ; QUITTER bouc1  ; FINSI;
  754.       SI (EGA ret 'precedent') ; ind = ind - 1; FINSI;
  755.       SI (EGA ret 'suivant')   ; ind = ind + 1; FINSI;
  756.       SI (EGA ret 'retracer')  ; ind = ind    ; FINSI;
  757.       SI (EGA ret 'Quitter')   ; QUITTER bouc2; FINSI;
  758.       SI (EGA ret 'n°') ; ind = SAIS 'Numéro du pas : ' 'ENTIER'; FINSI;
  759.     FINSI;
  760.  
  761. *   indice du pas    
  762.     SI (ind > ltindex); ind = 0      ; FINSI;
  763.     SI (ind < 0);       ind = ltindex; FINSI;
  764.  
  765. *   Si le modele evolue au cours du calcul :
  766.     SI ICHGMOD1 ;
  767.       SI ( EXIS matab 'CONTRAINTES');
  768.         modmec = EXTR  (MATAB . 'MODELES'  . ind) 'FORMULATION' 'MECANIQUE';
  769.             lemail = extr modmec 'MAIL' ;
  770. *  enlever les contacts
  771.         SI ICHGMAT1 ;
  772.           mati1  = tire (MATAB . 'CHARGEMENT') 'MATE' (MATAB . 'TEMPS' . ind) ;
  773.         SINON ;
  774.           mati1  = matmec ;
  775.         FINSI ;
  776.         matmec = REDU mati1 modmec;
  777.       FINSI ;
  778.       SI tempe ;
  779.         modthe = EXTR  (MATAB . 'MODELES'  . ind) 'FORMULATION' 'THERMIQUE';
  780.         modthe = EXTR  modthe 'MATERIAU' 'CONDUCTION' ;
  781.         maithe = EXTR  modthe 'MAIL' ;
  782.       FINSI ;
  783.     FINSI ;
  784.  
  785. *   objet a tracer
  786.     SI defor; dpl = DEFO lemail tdpl.ind ampdef;FINSI;
  787.     SI contr; 
  788.       ctr = tctr.ind; 
  789.       moctr = mot 'contraintes';
  790.     SINON; 
  791.       SI var_int;
  792.         ctr = tvar.ind;       
  793.         moctr = mot 'variables internes';
  794.       FINSI; 
  795.     FINSI;
  796.     SI react; vec = VECT trea.ind leforces;FINSI;
  797.     SI (defor ET react);dpv = DEFO lemail tdpl.ind ampdef vec;FINSI;
  798.  
  799. *   tracer
  800.     SI (defor ET (NON contr) ET (NON react) et (non tempe));
  801.       TITRE 'deforme au temps ' MATAB . 'TEMPS' . ind;
  802.       TRAC lemail dpl;
  803.     FINSI;
  804.     SI ((NON defor) ET contr ET (NON react));
  805.       TITRE moctr ' au temps ' MATAB . 'TEMPS' . ind;
  806.       TRAC ctr modmec matmec ;
  807.     FINSI;
  808.     SI (defor ET contr ET (NON react));
  809.       TITRE moctr ' au temps ' MATAB . 'TEMPS' . ind;
  810.       TRAC ctr modmec matmec dpl;
  811.     FINSI;
  812.     SI ((NON defor) ET (NON contr) ET react);
  813.       TITRE 'reactions au temps ' MATAB . 'TEMPS' . ind;
  814.       TRAC lemail vec;
  815.     FINSI;
  816.     SI (defor ET (NON contr) ET react);
  817.       TITRE 'deforme au temps ' MATAB . 'TEMPS' . ind;
  818.       TRAC dpv lemail;
  819.     FINSI;
  820.     SI ((NON defor) ET contr ET react);
  821.       TITRE moctr ' au temps ' MATAB . 'TEMPS' . ind;
  822.       TRAC ctr modmec matmec lemail vec;
  823.     FINSI;
  824.     SI (defor ET contr ET react);
  825.       TITRE 'contraintes au temps ' MATAB . 'TEMPS' . ind;
  826.       TRAC ctr modmec matmec  dpv ;
  827.     FINSI;
  828.     SI ((NON defor) et tempe) ;
  829.       tctemp=tctem.ind;  
  830.       TITRE ' Temperatures au temps ' MATAB . 'TEMPS' . ind;
  831.       TRAC tctemp  maithe ;
  832.     FINSI; 
  833.     SI (defor et tempe) ;
  834.       tctemp=tctem.ind;  
  835.       TITRE ' Temperatures au temps ' MATAB . 'TEMPS' . ind;
  836.       TRAC tctemp  dpl ;
  837.     FINSI; 
  838.  
  839.     SI (NON FLCHOI);
  840.       ind = ind + 1;
  841.       SI (ind > ltindex); QUITTER bouc1; FINSI;
  842.     FINSI;
  843.  
  844.   FIN bouc1;
  845. * boucle second niveau sur les pas de temps ......  
  846.  
  847. FINSI;
  848. * FIN DISTINCTION CHOIX = EVOLU / AUTRE QUE EVOLU
  849.  
  850. FIN bouc2;
  851. * FIN BOUCLE PRINCIPALE ------------------------------------------------
  852.  
  853.  
  854. FINPROC  ;
  855.  
  856.  
  857.  
  858.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales