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. OPTI REST FORMAT
  38. '/u/castem/divers/plexus1.couplage' ;
  39.  
  40. TITRE 'TEST NTERFACE CAST3M - PLEXUS';
  41.  
  42. *-Restitution de la bande cree par PLEXUS
  43. TEMPS ;
  44. REST FORMAT ;
  45.  
  46.  
  47.  
  48. SI (NEG GRAPH 'N');
  49. TITR 'MAILLAGE TOUT DANS PLEXUS1.DGIBI';
  50. TRAC 'QUAL' TOUT ;
  51. FINSI;
  52.  
  53. * TESTS DES CHPO CREES PAR PLEXUS
  54. * ===============================
  55. ND5 = NOEU 5 ;
  56. ND11 = NOEU 11 ;
  57. *
  58. * creation d un chpo de deplacements et
  59. * de vitesses ( 2 et 3eme sauvegarde)
  60.  
  61. DEPl2 = TABPLEX . 2 . DEPL ;
  62. VITE3 = TABPLEX . 3 . VITE ;
  63.  
  64. SI (NEG GRAPH 'N');
  65. TITR 'CHAMP DE DEPLACEMENT au pas 2' ;
  66. TRAC DEPl2 TOUT ;
  67. TITR 'CHAMP DE VITESSE au pas 3' ;
  68. TRAC VITE3 TOUT ;
  69. FINSI;
  70.  
  71. * extraction du deplacement suivant r pour
  72. * le noeud no 5 et de la vitesse suivant z
  73. * du noeud 11
  74.  
  75. FD5R = EXTR DEPL2 ND5 UR ;
  76. FV11Z = EXTR VITE3 ND11 UZ ;
  77.  
  78. * TESTS DES MCHAML CREES PAR PLEXUS
  79. * =================================
  80.  
  81. CONT3 = TABPLEX . 3 . CONT ;
  82. EPST2 = TABPLEX . 2 . EPTS ;
  83.  
  84. * on extrait du champ de contrainte et du
  85. * champ des defo. totales la composante S_2
  86. *(resp. D_1) pour laquelle on change le nom
  87. * de la composante.
  88.  
  89. CHN22 = EXCO S_2 CONT3 N22 ;
  90. CHD11 = EXCO D_1 EPST2 D11 ;
  91.  
  92. * on cree un objet de type MODE , un objet
  93. * de type MATE associes a la 1 ere zone du
  94. * maillage ZONE_1
  95.  
  96. MODL1=MODE ZONE_1 MECANIQUE ELASTIQUE COQ2;
  97. CAR1 =MATE MODL1 'YOUN' 200.E3 'NU' 0.
  98. 'EPAI' 2. 'RHO' 7.8E-9;
  99.  
  100. * on reduit le nouveau chamelem sur la
  101. * premiere zone du maillage ZONE_1
  102.  
  103. CHN22R=REDU CHN22 ZONE_1;
  104.  
  105. * on extrait la composante N22 du 1er
  106. * sous-paquet du 1er element pour le 1er
  107. * point de gauss (dans notre cas le 1er
  108. * noeud : mchaml au noeuds en sortie de
  109. * PLEXUS)
  110.  
  111. N22_1 = EXTR CHN22R N22 1 1 1 ;
  112. TEMPS ;
  113. VAL1 = 1.5492E-5 ;
  114. VAL2 = 335.39 ;
  115. VAL3 = 2.2031 ;
  116. RESI1 = 100*(ABS((FD5R - VAL1 )/ VAL1 ));
  117. RESI2 = 100*(ABS((FV11Z - VAL2 )/ VAL2 ));
  118. RESI3 = 100*(ABS((N22_1 - VAL3 )/ VAL3 ));
  119.  
  120. MESS 'deplacement theorique :' VAL1 'mm';
  121. MESS 'deplacement calculee :' FD5R 'mm';
  122. MESS ' Soit un ecart de : ' RESI1 '%';
  123. SAUTER 1 LIGNES ;
  124.  
  125. MESS 'vitesse theorique:' VAL2 'mm/s';
  126. MESS 'vitesse calculee :' FV11Z 'm/s';
  127. MESS ' Soit un ecart de : ' RESI2 '%';
  128. SAUTER 1 LIGNES ;
  129.  
  130. MESS 'contrainte theorique:' VAL3 'MPa';
  131. MESS 'contrainte calculee :' N22_1 'MPa';
  132. MESS ' Soit un ecart de : ' RESI3 '%';
  133. SAUTER 1 LIGNES ;
  134.  
  135.  
  136. * CODE FONCTIONNEMENT
  137. * -------------------
  138.  
  139. SI((RESI1 <EG 1E-2) ET (RESI2 <EG 1E-2) ET
  140. (RESI3 <EG 1E-2));
  141. ERRE 0;
  142. SINO;
  143. ERRE 5;
  144. FINSI;
  145. FIN;
  146.  
  147.  
  148.  
  149.  
  150.  

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