Télécharger @FIS_3DS.procedur

Retour à la liste

Numérotation des lignes :

  1. * @FIS_3DS PROCEDUR CHAT 97/09/25 21:15:11 2851
  2. *
  3. *************************************************
  4. * VERSION 3.0 fournie par EDF/MMM *
  5. * Version du 15/06/96 *
  6. *************************************************
  7. *
  8. * parametres de la procedure principale :
  9. * ---------------------------------------
  10. *
  11. * c demi grand axe de l'ellipse
  12. * a demi petit axe de l'ellipse
  13. * rc0 rayon du tore
  14. * rc1 coefficient multiplicateur du parametre rc0
  15. * definissant l'epaisseur de la premiere couronne
  16. * de deraffinement
  17. * rc2 coefficient multiplicateur du parametre rc0
  18. * definissant l'epaisseur de la deuxieme couronne
  19. * de deraffinement
  20. * rc3 coefficient multiplicateur du parametre rc0
  21. * definissant l'epaisseur de la troisieme couronne
  22. * de deraffinement (si ndt=2)
  23. * nc nombre de couronnes
  24. * ns nombre de secteurs sur 90 degres
  25. * nt nombre de divisions sur un quart d'ellipse
  26. * eps demi-angle d'ouverture de la fissure (degres)
  27. * lo longueur du bloc
  28. * to largeur du bloc
  29. * ho hauteur du bloc
  30. * ndt nombre de couronnes de deraffinement (1 ou 2)
  31. * nsdt Nombre de secteurs sur 90 degres au niveau des
  32. * couronnes de deraffinement des tranches (2 ou 4)
  33. * beta impose le decoupage le long de la generatrice
  34. * alpha impose l'angle des differentes tranches
  35. * xl impose le nombre d'elements pour la prolongation
  36. * du bloc initial suivant l'axe x (longueur)
  37. * xt impose le nombre d'elements pour la prolongation
  38. * du bloc initial suivant l'axe y (largeur)
  39. * xh impose le nombre d'elements pour la prolongation
  40. * du bloc initial suivant l'axe -z (hauteur)
  41. * d1 'oui' si deplacement des noeuds intermediares du
  42. * boudin
  43. * d2 'oui' si maillage 'rayonnant' sur la face
  44. * exterieure au boudin torique
  45. * On doit avoir d2 = 'non' dans le cas de l'utilisation
  46. * de la procedure plaque epaisse fissuree
  47. debproc @FIS_3DS c*flottant a*flottant lo*flottant
  48. to*flottant ho*flottant
  49. nt*entier nc*entier ns*entier
  50. rc0*flottant rc1*flottant rc2*flottant rc3*flottant
  51. alpha*flottant
  52. ndt*entier nsdt*entier
  53. xl*entier xt*entier xh*entier ;
  54. *
  55. * DEBUT PARAMETRES UTILISATEUR FIXES
  56. *
  57. eps = 0.1 ;
  58. beta = 1. ;
  59. d1 = 'oui' ;
  60. d2 = 'non' ;
  61. *
  62. * FIN PARAMETRES UTILISATEUR FIXES
  63. *
  64. si (c < a) ;
  65. saut ligne ;
  66. mess '================================================== ' ;
  67. mess ' ATTENTION !!! ' ;
  68. saut ligne ;
  69. mess ' Le premier parametre (demi grand axe) doit etre ' ;
  70. mess ' superieur aux deuxieme parametre (demi petit axe)' ;
  71. mess ' Demi grand axe : c =' c ;
  72. mess ' Demi petit axe : a =' a ;
  73. saut ligne ;
  74. mess ' ------ ON SORT DE LA PROCEDURE ----- ' ;
  75. mess '================================================== ' ;
  76. quit @FIS_3DS ;
  77. finsi ;
  78. *
  79. asurc = a / c ;
  80. si ((asurc < .09) ou (asurc > 1.01)) ;
  81. saut ligne ;
  82. mess '================================================== ' ;
  83. mess ' ATTENTION !!! ' ;
  84. saut ligne ;
  85. mess ' L"allongement de l"ellipse : ' ;
  86. mess ' (demi petit axe)/(demi grand axe) doit etre ' ;
  87. mess ' compris entre 0,1 et 1 ' ;
  88. mess ' (a / c) =' asurc ;
  89. saut ligne ;
  90. mess ' ------ ON SORT DE LA PROCEDURE ----- ' ;
  91. mess '================================================== ' ;
  92. quit @FIS_3DS ;
  93. finsi ;
  94. *
  95. asurto = a / to ;
  96. si ((asurto < .0999) ou (asurto > .8001)) ;
  97. saut ligne ;
  98. mess '================================================== ' ;
  99. mess ' ATTENTION !!! ' ;
  100. saut ligne ;
  101. mess ' Le rapport a/to (demi-petit axe / largeur du bloc)' ;
  102. mess ' doit etre compris entre 0,1 et 0,8 ' ;
  103. mess ' (a/to) =' a/to ;
  104. saut ligne ;
  105. mess ' ------ ON SORT DE LA PROCEDURE ----- ' ;
  106. mess '================================================== ' ;
  107. quit @FIS_3DS ;
  108. finsi ;
  109. *
  110. si ((nc < 1) ou (nc > 6)) ;
  111. saut ligne ;
  112. mess '================================================== ' ;
  113. mess ' ATTENTION !!! ' ;
  114. saut ligne ;
  115. mess ' Le nombre de couronnes du "boudin torique" ' ;
  116. mess ' doit etre compris entre 1 et 6 ' ;
  117. mess ' nc =' nc ;
  118. saut ligne ;
  119. mess ' ------ ON SORT DE LA PROCEDURE ----- ' ;
  120. mess '================================================== ' ;
  121. quit @FIS_3DS ;
  122. finsi ;
  123. *
  124. si ((ns < 2) ou (ns > 9)) ;
  125. saut ligne ;
  126. mess '================================================== ' ;
  127. mess ' ATTENTION !!! ' ;
  128. saut ligne ;
  129. mess ' Le nombre de secteurs du "boudin torique" ' ;
  130. mess ' doit etre compris entre 2 et 9 (sur 90 degres) ' ;
  131. mess ' ns =' ns ;
  132. saut ligne ;
  133. mess ' ------ ON SORT DE LA PROCEDURE ----- ' ;
  134. mess '================================================== ' ;
  135. quit @FIS_3DS ;
  136. finsi ;
  137. *
  138. si ((rc1 < .5) ou (rc1 > 3.)) ;
  139. saut ligne ;
  140. mess '================================================== ' ;
  141. mess ' ATTENTION !!! ' ;
  142. saut ligne ;
  143. mess ' Le parametre rc1 doit etre compris entre 0,5 et 3 ' ;
  144. mess ' rc1 =' rc1 ;
  145. saut ligne ;
  146. mess ' ------ ON SORT DE LA PROCEDURE ----- ' ;
  147. mess '================================================== ' ;
  148. quit @FIS_3DS ;
  149. finsi ;
  150. *
  151. si ((rc2 < .5) ou (rc2 > 5.)) ;
  152. saut ligne ;
  153. mess '================================================== ' ;
  154. mess ' ATTENTION !!! ' ;
  155. saut ligne ;
  156. mess ' Le parametre rc2 doit etre compris entre 0,5 et 5 ' ;
  157. mess ' rc2 =' rc2 ;
  158. saut ligne ;
  159. mess ' ------ ON SORT DE LA PROCEDURE ----- ' ;
  160. mess '================================================== ' ;
  161. quit @FIS_3DS ;
  162. finsi ;
  163. *
  164. si ((rc3 < .5) ou (rc3 > 5.)) ;
  165. saut ligne ;
  166. mess '================================================== ' ;
  167. mess ' ATTENTION !!! ' ;
  168. saut ligne ;
  169. mess ' Le parametre rc3 doit etre compris entre 0,5 et 5 ' ;
  170. mess ' rc3 =' rc3 ;
  171. saut ligne ;
  172. mess ' ------ ON SORT DE LA PROCEDURE ----- ' ;
  173. mess '================================================== ' ;
  174. quit @FIS_3DS ;
  175. finsi ;
  176. *
  177. vtot lff levreinf sar slat sinf sav_s ssup_s slat_s boudin ep3
  178. = @FIS_3 c a lo to ho eps nt nc ns rc0 rc1 rc2 rc3
  179. beta alpha ndt nsdt xl xt xh d1 d2 ;
  180. *
  181. finproc vtot lff levreinf sar slat sinf sav_s ssup_s
  182. slat_s boudin ep3 ;
  183.  
  184.  

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