Télécharger la plaquette Cast3M


fr en

Télécharger ortho-coq4.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  ortho-coq4.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *------------------------------------------------------------------
  5. *
  6. *              TEST ORTHOTROPIE : PLAQUE EN FLEXION
  7. *   Michel BULIK (inspiré d'un test d'Alain MOAL (octobre 1996))
  8. *
  9. *   Il s'agit d'une plaque carrée, encastrée sur les bords et
  10. *   soumise à une pression uniforme. On modélise succésivement
  11. *   les cas suivants :
  12. *
  13. *   1. 1er axe d'orthotropie = Ox
  14. *      E_1 = 2.e+5
  15. *      E_2 = 2.e+4
  16. *
  17. *   2. 1er axe d'orthotropie = Oy
  18. *      E_1 = 2.e+4
  19. *      E_2 = 2.e+5
  20. *
  21. *   A cause de la symétrie du problème ces 2 cas sont équivalents
  22. *   aux 2 cas suivants :
  23. *
  24. *   3. 1er axe d'orthotropie = Oy
  25. *      E_1 = 2.e+5
  26. *      E_2 = 2.e+4
  27. *
  28. *   4. 1er axe d'orthotropie = Ox
  29. *      E_1 = 2.e+4
  30. *      E_2 = 2.e+5
  31. *
  32. *   Le test porte sur l'égalité des flèches au centre de la plaque
  33. *   dans les quatre cas.
  34. *
  35. *------------------------------------------------------------------
  36.  
  37. opti dime 3 mode trid elem qua4 ;
  39. graph = faux ;
  40.  
  41. *---- maillage
  42. p1 = -1.e3 -1.e3 0. ;
  43. p2 =  1.e3 -1.e3 0. ;
  44. p3 =  1.e3  1.e3 0. ;
  45. p4 = -1.e3  1.e3 0. ;
  46.  
  47. n1 = 4 ;
  48. n2 = 4 ;
  49.  
  50. l1 = p1 d n1 p2 ;
  51. l2 = p2 d n2 p3 ;
  52. l3 = p3 d n1 p4 ;
  53. l4 = p4 d n2 p1 ;
  54.  
  55. s0 = daller l1 l2 l3 l4 plan ;
  56.  
  57. *---- calcul
  58.  
  59. *---- cas 1
  60. youn1 =  200000. ;
  61. youn2 =   20000. ;
  62. youn0 = (youn1 + youn2)/2.;
  63. gcis0 = youn0 /(2.*(1.+0.3)) ;
  64. valnu12 = 0.3 ; comm '=> valnu21 = 0.03' ;
  65.  
  66. mod0 = MODE s0 mecanique elastique orthotrope coq4 ;
  67. mat0 = MATE mod0 direction (1. 0. 0.) parallele
  68.             YG1 youn1 YG2 youn2
  69.             NU12 valnu12
  70.             G12 gcis0
  71.             G23 gcis0
  72.             G13 gcis0
  73.             RHO 8.e-6 EPAI 50. ;
  74.  
  75. MOP = 'MODE' S0 'CHARGEMENT' 'PRESSION' 'COQ4' ;
  76. MAP = 'PRES' MOP 'PRES' -0.5 ;
  77. MODT = MOD0 'ET' MOP ;
  78. MATT = MAT0 ;
  79. pres0 = 'BSIG' MOP MAP ;
  80.  
  81.  
  82. rig0 = (rigi modT matT) et
  83.        (bloq depl rota (l1 et l2 et l3 et l4)) ;
  84.  
  85. dep0 = reso rig0 pres0 ;
  86.  
  87. *---- depouillement
  88. defo0 = defo s0 dep0 1. roug ;
  89. depz0 = exco uz dep0 ;
  90. fle1 = extr depz0 scal (s0 poin proc (0. 0. 0.)) ;
  91. mess 'fleche = ' fle1 ;
  92. titre 'Deplacement en z et deformee (fleche 1 = 'fle1' mm)';
  93. si graph ;
  94.    trace cach (1.e8 1.e8 1.e8) depz0 defo0 ;
  95. finsi ;
  96.  
  97. *---- cas 2
  98. youn1 =  20000. ;
  99. youn2 = 200000. ;
  100. youn0 = (youn1 + youn2)/2.;
  101. gcis0 = youn0 /(2.*(1.+0.3)) ;
  102. valnu12 = 0.03 ; comm '=> valnu21 = 0.3' ;
  103.  
  104. mat0 = MATE mod0 direction (1. 0. 0.) parallele
  105.             YG1 youn1 YG2 youn2
  106.             NU12 valnu12
  107.             G12 gcis0
  108.             G23 gcis0
  109.             G13 gcis0
  110.             RHO 8.e-6 EPAI 50. ;
  111. MODT = MOD0 'ET' MOP ;
  112. MATT = MAT0 ;
  113.  
  114. rig0 = (rigi modT matT) et
  115.        (bloq depl rota (l1 et l2 et l3 et l4)) ;
  116.  
  117. dep0 = reso rig0 pres0 ;
  118.  
  119. *---- depouillement
  120. defo0 = defo s0 dep0 1. roug ;
  121. depz0 = exco uz dep0 ;
  122. fle2 = extr depz0 scal (s0 poin proc (0. 0. 0.)) ;
  123. mess 'fleche = ' fle2 ;
  124. titre 'Deplacement en z et deformee (fleche 2 = 'fle2' mm)';
  125. si graph ;
  126.    trace cach (1.e8 1.e8 1.e8) depz0 defo0 ;
  127. finsi ;
  128.  
  129. *---- cas 3
  130. youn1 = 200000. ;
  131. youn2 =  20000. ;
  132. youn0 = (youn1 + youn2)/2.;
  133. gcis0 = youn0 /(2.*(1.+0.3)) ;
  134. valnu12 = 0.3 ; comm '=> valnu21 = 0.03' ;
  135.  
  136. mat0 = MATE mod0 direction (0. 1. 0.) parallele
  137.             YG1 youn1 YG2 youn2
  138.             NU12 valnu12
  139.             G12 gcis0
  140.             G23 gcis0
  141.             G13 gcis0
  142.             RHO 8.e-6 EPAI 50. ;
  143. MODT = MOD0 'ET' MOP ;
  144. MATT = MAT0 ;
  145. rig0 = (rigi modT matT) et
  146.        (bloq depl rota (l1 et l2 et l3 et l4)) ;
  147.  
  148. dep0 = reso rig0 pres0 ;
  149.  
  150. *---- depouillement
  151. defo0 = defo s0 dep0 1. roug ;
  152. depz0 = exco uz dep0 ;
  153. fle3 = extr depz0 scal (s0 poin proc (0. 0. 0.)) ;
  154. mess 'fleche = ' fle3 ;
  155. titre 'Deplacement en z et deformee (fleche 3 = 'fle3' mm)';
  156. si graph ;
  157.    trace cach (1.e8 1.e8 1.e8) depz0 defo0 ;
  158. finsi ;
  159.  
  160. *---- cas 4
  161. youn1 =  20000. ;
  162. youn2 = 200000. ;
  163. youn0 = (youn1 + youn2)/2.;
  164. gcis0 = youn0 /(2.*(1.+0.3)) ;
  165. valnu12 = 0.03 ; comm '=> valnu21 = 0.3' ;
  166.  
  167. mat0 = MATE mod0 direction (0. 1. 0.) parallele
  168.             YG1 youn1 YG2 youn2
  169.             NU12 valnu12
  170.             G12 gcis0
  171.             G23 gcis0
  172.             G13 gcis0
  173.             RHO 8.e-6 EPAI 50. ;
  174. MODT = MOD0 'ET' MOP ;
  175. MATT = MAT0 ;
  176. rig0 = (rigi modT matT) et
  177.        (bloq depl rota (l1 et l2 et l3 et l4)) ;
  178.  
  179. dep0 = reso rig0 pres0 ;
  180.  
  181. *---- depouillement
  182. defo0 = defo s0 dep0 1. roug ;
  183. depz0 = exco uz dep0 ;
  184. fle4 = extr depz0 scal (s0 poin proc (0. 0. 0.)) ;
  185. mess 'fleche = ' fle4 ;
  186. titre 'Deplacement en z et deformee (fleche 4 = 'fle4' mm)';
  187. si graph ;
  188.    trace cach (1.e8 1.e8 1.e8) depz0 defo0 ;
  189. finsi ;
  190.  
  191. *--------- test si c'est ok
  192. opti echo 0 ;
  194. diff12 = abs(fle2 - fle1) ;
  195. mess 'Différence entre les résultats 1 et 2 =' diff12 ;
  196. diff23 = abs(fle3 - fle2) ;
  197. mess 'Différence entre les résultats 2 et 3 =' diff23 ;
  198. diff34 = abs(fle4 - fle3) ;
  199. mess 'Différence entre les résultats 3 et 4 =' diff34 ;
  200. epsilo = 1.e-10 ;
  201. testok = (diff12 < epsilo) ET (diff23 < epsilo) ET (diff34 < epsilo) ;
  202. si (non testok) ;
  203.    erre 5 ;
  204. finsi ;
  205.  
  206. fin ;
  207.  
  208.  
  209.  
  210.  
  211.  
  212.  
  213.  
  214.  
  215.  
  216.  
  217.  

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