Télécharger la plaquette Cast3M


fr en

Télécharger plexus1.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  plexus1.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * CAS TEST DU 92/01/16    PROVENANCE : PLA2
  5.  
  6. *******************************************
  7. *  Test plexus1.dgibi: jeux de données    *
  8. *    ---------------------------------    *
  9. *******************************************
  10.  
  11. * SI GRAPH = N PAS DE GRAPHIQUE AFFICHE
  12. * SINON SI GRAPH DIFFERENT DE N TOUS
  13. * LES GRAPHIQUES SONT AFFICHES
  14.  
  15. GRAPH = 'N' ;
  16.  
  17. SAUT PAGE;
  18. SI (NEG GRAPH 'N') ;
  19.   OPTI ECHO 1   ;
  20.   OPTI TRAC X ;
  21. SINO ;
  22.   OPTI ECHO 0   ;
  23. FINSI ;
  24.  
  25. SAUT PAGE;
  26.  
  27. *******************************************
  28. *                PLEXUS1                  *
  29. *                                         *
  30. *       CAS TEST DE L INTERFACE           *
  31. *        CASTEM 2000 - PLEXUS             *
  32. *                                         *
  33. * Cet exemple permet de tester la         *
  34. *  relecture de la bande CAST3M cree par  *
  35. *  PLEXUS                                 *
  36. *******************************************
  37.  
  38. * repertoire des fichiers "divers"
  39. DIVERS = VENV 'CASTEM_DIVERS';
  40. *
  41. OPTI REST FORMAT
  42.       ('CHAINE' DIVERS '/plexus1.couplage') ;
  43.  
  44. TITRE 'TEST NTERFACE CAST3M - PLEXUS';
  45.  
  46. *-Restitution de la bande cree par PLEXUS
  47. TEMPS ;
  48. REST FORMAT ;
  49.  
  50.  
  51.  
  52. SI (NEG GRAPH 'N');
  53.    TITR 'MAILLAGE TOUT DANS PLEXUS1.DGIBI';
  54.    TRAC 'QUAL' TOUT  ;
  55. FINSI;
  56.  
  57. *    TESTS DES CHPO CREES PAR PLEXUS
  58. *    ===============================
  59. ND5  = NOEU 5  ;
  60. ND11 = NOEU 11 ;
  61. *
  62. * creation d un chpo de deplacements et
  63. * de vitesses ( 2 et 3eme sauvegarde)
  64.  
  65. DEPl2 = TABPLEX . 2 . DEPL ;
  66. VITE3 = TABPLEX . 3 . VITE ;
  67.  
  68. SI (NEG GRAPH 'N');
  69.    TITR 'CHAMP DE DEPLACEMENT au pas 2' ;
  70.    TRAC DEPl2 TOUT  ;
  71.    TITR 'CHAMP DE VITESSE au pas 3' ;
  72.    TRAC VITE3 TOUT  ;
  73. FINSI;
  74.  
  75. * extraction du deplacement suivant r pour 
  76. * le noeud no 5 et de la vitesse suivant z 
  77. * du noeud 11
  78.  
  79. FD5R  = EXTR DEPL2 ND5  UR ;
  80. FV11Z = EXTR VITE3 ND11 UZ ;
  81.  
  82. *    TESTS DES MCHAML CREES PAR PLEXUS
  83. *    =================================
  84.  
  85. CONT3 = TABPLEX . 3 . CONT ;
  86. EPST2 = TABPLEX . 2 . EPTS ;
  87.  
  88. * on extrait du champ de contrainte et du 
  89. * champ des defo. totales la composante S_2
  90. *(resp. D_1) pour laquelle on change le nom
  91. * de la composante.
  92.  
  93. CHN22 = EXCO S_2 CONT3 N22 ;
  94. CHD11 = EXCO D_1 EPST2 D11 ;
  95.  
  96. * on cree un objet de type MODE , un objet 
  97. * de type MATE associes a la 1 ere zone du 
  98. * maillage ZONE_1
  99.  
  100. MODL1=MODE ZONE_1 MECANIQUE ELASTIQUE COQ2;
  101. CAR1 =MATE MODL1 'YOUN' 200.E3 'NU' 0. 
  102.                  'EPAI' 2.    'RHO' 7.8E-9;
  103.  
  104. * on reduit le nouveau chamelem sur la 
  105. * premiere zone du maillage ZONE_1
  106.  
  107. CHN22R=REDU CHN22 ZONE_1;
  108.  
  109. * on extrait la composante N22 du 1er 
  110. * sous-paquet du 1er element pour le 1er 
  111. * point de gauss (dans notre cas le 1er 
  112. * noeud : mchaml au noeuds en sortie de 
  113. * PLEXUS)
  114.  
  115. N22_1 = EXTR CHN22R N22 1 1 1 ;
  116. TEMPS ;
  117. VAL1 = 1.5492E-5 ;
  118. VAL2 = 335.39 ;
  119. VAL3 = 2.2031 ;
  120. RESI1 = 100*(ABS((FD5R  - VAL1 )/ VAL1 ));
  121. RESI2 = 100*(ABS((FV11Z - VAL2 )/ VAL2 ));
  122. RESI3 = 100*(ABS((N22_1 - VAL3 )/ VAL3 ));
  123.  
  124. MESS 'deplacement theorique :' VAL1 'mm';
  125. MESS 'deplacement calculee  :' FD5R 'mm';
  126. MESS '    Soit un ecart de : ' RESI1 '%';
  127. SAUTER 1 LIGNES ;
  128.  
  129. MESS 'vitesse theorique:' VAL2    'mm/s';
  130. MESS 'vitesse calculee :' FV11Z    'm/s';
  131. MESS '    Soit un ecart de : ' RESI2 '%';
  132. SAUTER 1 LIGNES ;
  133.  
  134. MESS 'contrainte theorique:' VAL3  'MPa';
  135. MESS 'contrainte calculee :' N22_1 'MPa';
  136. MESS '    Soit un ecart de : ' RESI3 '%';
  137. SAUTER 1 LIGNES ;
  138.  
  139.  
  140. *           CODE FONCTIONNEMENT
  141. *           -------------------
  142.  
  143. SI((RESI1 <EG 1E-2) ET (RESI2 <EG 1E-2) ET
  144.    (RESI3 <EG 1E-2));
  145.   ERRE  0;
  146. SINO;
  147.   ERRE  5;
  148. FINSI;
  149. FIN;
  150.  
  151.  
  152.  
  153.  
  154.  
  155.  
  156.  

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