Télécharger @mesu.procedur

Retour à la liste

Numérotation des lignes :

  1. * @MESU PROCEDUR FANDEUR 14/10/10 21:15:01 8178
  2. *-------------------------------------------------
  3. ****** PROCEDURE @MESU ******
  4. *-------------------------------------------------
  5. *
  6. *
  7. * CETTE PROCEDURE A ETE MISE GRACIEUSEMENT
  8. * A DISPOSOTION DE LA COMMUNAUTE CASTEM2000
  9. * PAR M. LIBEYRE ( CEA/DSM/DRFC )
  10. *
  11. * TEL : ( 33 1 ) 42 25 46 03
  12. *
  13. *-------------------------------------------------
  14. *
  15. * MESURE D'UN MAILLAGE (LONGUEUR, AIRE OU VOLUME)
  16. *
  17. *--------------------------------------------------------
  18. DEBPROC @MESU GEO*MAILLAGE ;
  19. LTY = GEO ELEM TYPE ;
  20. ITCAR=FAUX;
  21. * mess 'liste des types d elements :' ;
  22. * list LTY ;
  23. NTY = DIME LTY ;
  24. I = 0 ;
  25. E1 = FAUX ;
  26. E2 = FAUX ;
  27. E3 = FAUX ;
  28. REPETER BOUC NTY ;
  29. I = I + 1 ;
  30. ITY = EXTR LTY I ;
  31. EE1 = (EGA ITY SEG2) OU (EGA ITY SEG3) ;
  32. EE2 = (EGA ITY TRI3) OU (EGA ITY TRI6) OU (EGA ITY QUA4)
  33. OU (EGA ITY QUA8) ;
  34. EE3 = (EGA ITY CU20) OU (EGA ITY PRI6) OU (EGA ITY PR15)
  35. OU (EGA ITY PYR5) OU (EGA ITY PY13) OU (EGA ITY TET4)
  36. OU (EGA ITY TE10) ;
  37. E1 = E1 OU EE1 ;
  38. E2 = E2 OU EE2 ;
  39. E3 = E3 OU EE3 ;
  40. FIN BOUC ;
  41. SI E1;
  42. * MESS 'MAILLAGE DE TYPE LIGNE' ;
  43. FINSI ;
  44. SI E2;
  45. * MESS 'MAILLAGE DE TYPE SURFACE' ;
  46. FINSI ;
  47. SI E3;
  48. * MESS 'MAILLAGE DE TYPE MASSIF' ;
  49. FINSI ;
  50. SI ((E1 ET E2) OU (E1 ET E3) OU (E2 ET E3)) ;
  51. MESS 'LE MAILLAGE DOIT ETRE D UN TYPE UNIQUE : LIGNE, '
  52. 'SURFACE OU MASSIF' ;
  53. P1 =0. ;
  54. SINON ;
  55. SI E1 ;
  56. IMA = 1 ;
  57. FINSI ;
  58. SI E2 ;
  59. IMA = 2 ;
  60. FINSI ;
  61. SI E3 ;
  62. IMA = 3 ;
  63. FINSI ;
  64. * MESS 'DIMENSION DE GEO :';
  65. * LIST IMA ;
  66. SI (EGA IMA 1) ;
  67. * MESS 'CALCUL DE LA LONGUEUR' ;
  68. GEO1 = CHANGER GEO SEG2 ;
  69. MODL1 = MODE GEO1 MECANIQUE ELASTIQUE BARR ;
  70. MAT1 = MATE MODL1 'RHO' 1. ;
  71. CAR1 = CARA MODL1 'SECT' 1. ;
  72. ITCAR=VRAI;
  73. FINSI ;
  74. SI (EGA IMA 2) ;
  75. * MESS 'CALCUL DE L AIRE' ;
  76. I = 0 ;
  77. REPETER BOUC2 NTY ;
  78. I = I + 1 ;
  79. GEOT = GEO ELEM (EXTR LTY I) ;
  80. SI (EGA I 1) ;
  81. GEO1 = CHANGER GEOT TRI3 ;
  82. SINON ;
  83. GEO1 = GEO1 ET (CHANGER GEOT TRI3) ;
  84. FINSI ;
  85. FIN BOUC2 ;
  86. IDI = VALEUR DIME ;
  87. SI (EGA IDI 2) ;
  88. MODL1 = 'MODE' GEO1 MECANIQUE ELASTIQUE TRI3 ;
  89. SINON ;
  90. MODL1 = 'MODE' GEO1 MECANIQUE ELASTIQUE COQ3 ;
  91. CAR1 = 'CARA' MODL1 'EPAI' 1. ;
  92. ITCAR=VRAI;
  93. FINSI ;
  94. MAT1 = 'MATE' MODL1 'RHO' 1. ;
  95. FINSI ;
  96. SI (EGA IMA 3) ;
  97. * MESS 'CALCUL DU VOLUME' ;
  98. I = 0 ;
  99. REPETER BOUC3 NTY ;
  100. I = I + 1 ;
  101. GEOT = GEO ELEM (EXTR LTY I) ;
  102. SI (EGA I 1) ;
  103. GEO1 = CHANGER GEOT TET4 ;
  104. SINON ;
  105. GEO1 = GEO1 ET (CHANGER GEOT TET4) ;
  106. FINSI ;
  107. FIN BOUC3 ;
  108. MODL1 = 'MODE' GEO1 MECANIQUE ELASTIQUE 'TET4' ;
  109. MAT1 = 'MATE' MODL1 'RHO' 1. ;
  110. FINSI ;
  111. SI (IMA < 3) ;
  112. SI ITCAR ;
  113. MASS1 = MASSE MODL1 (MAT1 ET CAR1) ;
  114. SINON ;
  115. MASS1 = MASSE MODL1 MAT1 ;
  116. FINSI ;
  117. SINON ;
  118. MASS1 = MASSE MODL1 MAT1 ;
  119. FINSI ;
  120. GEO2 = CHANGER GEO POI1 ;
  121. GAM1 = MANU CHPO GEO2 1 UX 1. ;
  122. POIP = MASS1 * GAM1 ;
  123. P1 = @TOTAL POIP GEO2 FX ;
  124. FINSI;
  125. FINPROC P1 ;
  126.  
  127.  

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