Télécharger elas20.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : elas20.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Elas20.dgibi: Jeux de données *
  5. * ---------------------------------- *
  6. * *
  7. SAUT PAGE ;
  8. *****************************************************
  9. * *
  10. * TEST TUYAU DROIT ET COUDE SOUS PRESSION *
  11. * *
  12. * description : *
  13. * ----------- *
  14. * *
  15. * Test en statique lineaire, chargement pression. *
  16. * *
  17. * La structure est modelisee en elements TUYAU : *
  18. * - 2 parties droites de longueur 3.048m maillées *
  19. * par 5 éléments chacune. *
  20. * - 1 coude entre ces deux parties droites de *
  21. * rayon de courbure égal à 0.61m, maillé par *
  22. * 10 éléments. *
  23. * *
  24. * Propriétés des matériaux : *
  25. * - Module d'Young : YOUN1 = 2.e11 *
  26. * - Coefficient de Poisson : NU1 = 0.3 *
  27. * *
  28. * Chargement : *
  29. * - Pression : PRES1 = 15 Mpa *
  30. * *
  31. * Conditions limites : *
  32. * - Encastrement a l'extrémité PA *
  33. * *
  34. *****************************************************
  35. OPTI ECHO 0 ;
  36. *
  37. OPTI DIME 3 ELEM SEG2 ;
  38. *
  39. GRAPH = 'N' ;
  40. *
  41. *****************************************************
  42. * GEOMETRIE : MAILLAGE *
  43. *****************************************************
  44. *
  45. OEIL = 0. 0. 10. ;
  46. *
  47. LDRO1 = 3.048 ;
  48. RCO1 = 0.61 ;
  49. LTOT1 = LDRO1 + RCO1 ;
  50. *
  51. PA = 0. 0. 0. ;
  52. PAP = 0. LDRO1 0. ;
  53. PBP = RCO1 LTOT1 0. ;
  54. PB = LTOT1 LTOT1 0. ;
  55. *
  56. LIG1 = PA D 5 PAP ;
  57. LIG2 = PBP D 5 PB ;
  58. LCOU1 = CERC 10 PAP ( RCO1 LDRO1 0. ) PBP ;
  59. *
  60. LIGNE = LIG1 ET LCOU1 ET LIG2 ;
  61. *
  62. TASS LIGNE ;
  63. ELIM LIGNE 1.E-5 ;
  64. *
  65. 'SI' ( 'NEG' GRAPH 'N' ) ;
  66. 'TITRE' 'TUYPRES : MAILLAGE' ;
  67. TRAC LIGNE OEIL QUAL ;
  68. 'FINSI' ;
  69. *
  70. *****************************************************
  71. * MODELES *
  72. *****************************************************
  73. *
  74. OBLI1 = MODE LIG1 MECANIQUE ELASTIQUE TUYAU ;
  75. OBLI2 = MODE LIG2 MECANIQUE ELASTIQUE TUYAU ;
  76. OBCO1 = MODE LCOU1 MECANIQUE ELASTIQUE TUYAU ;
  77. OBLIGN = OBLI1 ET OBLI2 ET OBCO1 ;
  78. *
  79. MOP1 = 'MODE' LIG1 'CHARGEMENT' 'PRESSION' TUYAU ;
  80. MOP2 = 'MODE' LIG2 'CHARGEMENT' 'PRESSION' TUYAU ;
  81. MOP3 = 'MODE' LCOU1 'CHARGEMENT' 'PRESSION' TUYAU ;
  82. MOPT = MOP1 'ET' MOP2 'ET' MOP3 ;
  83. *
  84. *****************************************************
  85. * VALEUR DE LA PRESSION *
  86. *****************************************************
  87. *
  88. PRES1 = 15.E6 ;
  89. *
  90. *****************************************************
  91. * MATERIAUX *
  92. *****************************************************
  93. *
  94. MATLI1 = 'MATE' ( OBLI1 ET OBLI2 ) 'YOUN' 2.1E11 'NU' 0.3 ;
  95. MATCO1 = 'MATE' OBCO1 'YOUN' 2.1E11 'NU' 0.3 ;
  96. MATLIGN = MATLI1 ET MATCO1 ;
  97. *
  98. MAP1 = 'MATE' MOP1 'PRES' PRES1 'EPAI' 26.187E-3 'RAYO' 203.2E-3;
  99. MAP2 = 'MATE' MOP2 'PRES' PRES1 'EPAI' 26.187E-3 'RAYO' 203.2E-3;
  100. MAP3 = 'MATE' MOP3 'PRES' PRES1 'EPAI' 26.187E-3 'RAYO' 203.2E-3
  101. 'RACO' 0.61;
  102. MAPT = MAP1 'ET' MAP2 'ET' MAP3 ;
  103. *
  104. *****************************************************
  105. * CARACTERISTIQUES *
  106. *****************************************************
  107. *
  108. CARLI1 = 'CARA' OBLI1 'EPAI' 26.187E-3 'RAYO' 203.2E-3
  109. 'VECT' ( 0. 0. -1. ) 'PRES' PRES1 ;
  110. CARLI2 = 'CARA' OBLI2 'EPAI' 26.187E-3 'RAYO' 203.2E-3
  111. 'VECT' ( 0. 0. -1. ) 'PRES' PRES1 ;
  112. CARCO1 = 'CARA' OBCO1 'EPAI' 26.187E-3 'RAYO' 203.2E-3
  113. 'RACO' 0.61 'VECT' ( 0. 0. -1. ) 'PRES' PRES1;
  114. CARLIGN = CARLI1 ET CARLI2 ET CARCO1 ;
  115. *
  116. CARTOT = MATLIGN ET CARLIGN ;
  117. *
  118. *****************************************************
  119. * CONDITIONS LIMITES *
  120. *****************************************************
  121. *
  122. RIGBLOQ = 'BLOQ' 'DEPL' 'ROTA' PA ;
  123. *
  124. *****************************************************
  125. * CALCUL DE LA RIGIDITE *
  126. *****************************************************
  127. *
  128. RIGLIGN = 'RIGI' OBLIGN CARTOT ;
  129. *
  130. RIGTOT = RIGLIGN ET RIGBLOQ ;
  131. *
  132. *****************************************************
  133. * CHARGEMENT *
  134. *****************************************************
  135. *
  136. PRZ = 'CHAN' 'TYPE' ('EXCO' MAPT 'PRES') 'CONTRAINTES';
  137. FPRES1 = 'PRES' MOPT PRZ MAPT ;
  138. *
  139. *****************************************************
  140. * RESOLUTION *
  141. *****************************************************
  142. *
  143. DEPL1 = 'RESO' RIGTOT FPRES1 ;
  144. *
  145. *****************************************************
  146. * EXTRACTION DES RESULTATS *
  147. *****************************************************
  148. *
  149. VXB1 = 0.822E-3 ;
  150. UXB1 = 'EXTR' DEPL1 'UX' PB ;
  151. ERG1 = 100. * ( 'ABS' (( VXB1 - UXB1 ) / VXB1 )) ;
  152. *
  153. VYB1 = VXB1 ;
  154. UYB1 = 'EXTR' DEPL1 'UY' PB ;
  155. ERG2 = 100. * ( 'ABS' (( VYB1 - UYB1 ) / VYB1 )) ;
  156. *
  157. VXBP1 = 0.137E-3 ;
  158. UXBP1 = 'EXTR' DEPL1 'UX' PBP ;
  159. ERG3 = 100. * ( 'ABS' (( VXBP1 - UXBP1 ) / VXBP1 )) ;
  160. *
  161. VYBP1 = 0.822E-3 ;
  162. UYBP1 = 'EXTR' DEPL1 'UY' PBP ;
  163. ERG4 = 100. * ( 'ABS' (( VYBP1 - UYBP1 ) / VYBP1 )) ;
  164. *
  165. VYAP1 = 0.685E-3 ;
  166. UYAP1 = 'EXTR' DEPL1 'UY' PAP ;
  167. ERG5 = 100. * ( 'ABS' (( VYAP1 - UYAP1 ) / VYAP1 )) ;
  168. *
  169. *****************************************************
  170. * AFFICHAGE DES RESULTATS *
  171. *****************************************************
  172. *
  173. 'MESS' 'RESULTATS : TEST TUYPRES' ;
  174. 'MESS' '------------------------' ;
  175. *
  176. 'MESS' ' deplacement ux (m) du point PB' ;
  177. 'MESS' ' ------------------------------' ;
  178. 'MESS' ' Valeur theorique : ' VXB1 'm' ;
  179. 'MESS' ' Valeur calculee : ' UXB1 'm' ;
  180. 'MESS' ' Ecart : ' ERG1 '%' ;
  181. *
  182. 'MESS' ' deplacement uy (m) du point PB' ;
  183. 'MESS' ' ------------------------------' ;
  184. 'MESS' ' Valeur theorique : ' VYB1 'm' ;
  185. 'MESS' ' Valeur calculee : ' UYB1 'm' ;
  186. 'MESS' ' Ecart : ' ERG2 '%' ;
  187. *
  188. 'MESS' ' deplacement ux (m) du point PBP' ;
  189. 'MESS' ' -------------------------------' ;
  190. 'MESS' ' Valeur theorique : ' VXBP1 'm' ;
  191. 'MESS' ' Valeur calculee : ' UXBP1 'm' ;
  192. 'MESS' ' Ecart : ' ERG3 '%' ;
  193. *
  194. 'MESS' ' deplacement uy (m) du point PBP' ;
  195. 'MESS' ' -------------------------------' ;
  196. 'MESS' ' Valeur theorique : ' VYBP1 'm' ;
  197. 'MESS' ' Valeur calculee : ' UYBP1 'm' ;
  198. 'MESS' ' Ecart : ' ERG4 '%' ;
  199. *
  200. 'MESS' ' deplacement uy (m) du point PAP' ;
  201. 'MESS' ' -------------------------------' ;
  202. 'MESS' ' Valeur theorique : ' VYAP1 'm' ;
  203. 'MESS' ' Valeur calculee : ' UYAP1 'm' ;
  204. 'MESS' ' Ecart : ' ERG5 '%' ;
  205. *
  206. ERGMAX = 'MAXI' ( 'PROG' ERG1 ERG2 ERG3 ERG4 ERG5 ) ;
  207. *
  208. 'SI' ( ERGMAX <EG 5.) ;
  209. ERRE 0 ;
  210. 'SINON' ;
  211. ERRE 5 ;
  212. 'FINSI' ;
  213. *
  214. TEMPS ;
  215. *
  216. FIN ;
  217.  
  218.  
  219.  
  220.  
  221.  
  222.  
  223.  
  224.  
  225.  
  226.  

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