Télécharger la plaquette Cast3M


fr en

Télécharger ortho-cu20.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  ortho-cu20.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *------------------------------------------------------------------
  5. *
  6. *     TEST ORTHOTROPIE : PLAQUE EN FLEXION - MAILLAGE VOLUMIQUE
  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. *   Le troisième axe d'orthotropie est toujours parallèle à l'axe
  22. *   Oz.
  23. *
  24. *   Le test porte sur l'égalité des flèches au centre de la plaque
  25. *   dans les deux cas (physiquement équivalents).
  26. *
  27. *------------------------------------------------------------------
  28.  
  29. opti dime 3 mode trid elem cu20 ;
  31. graph = faux ;
  32.  
  33. *---- maillage
  34. p1 = -1.e3 -1.e3 0. ;
  35. p2 =  1.e3 -1.e3 0. ;
  36. p3 =  1.e3  1.e3 0. ;
  37. p4 = -1.e3  1.e3 0. ;
  38.  
  39. n1 = 4 ;
  40. n2 = 4 ;
  41. n3 = 1 ;
  42.  
  43. l1 = p1 d n1 p2 ;
  44. l2 = p2 d n2 p3 ;
  45. l3 = p3 d n1 p4 ;
  46. l4 = p4 d n2 p1 ;
  47.  
  48. s0 = daller l1 l2 l3 l4 plan ;
  49. s1 = s0 moin (0. 0. 50.) ;
  50.  
  51. scl1 = l1 tran n3 (0. 0. -50.) ;
  52. scl2 = l2 tran n3 (0. 0. -50.) ;
  53. scl3 = l3 tran n3 (0. 0. -50.) ;
  54. scl4 = l4 tran n3 (0. 0. -50.) ;
  55. scl0 = scl1 et scl2 et scl3 et scl4 ;
  56.  
  57. v0 = s1 volu n3 s0 ;
  58. elim 0.1 scl0 v0 ;
  59.  
  60. *---- calcul
  61.  
  62. *---- cas 1
  63. youn1 =  200000. ;
  64. youn2 =   20000. ;
  65. youn3 =    2000. ;
  66. youn0 = (youn1 + youn2 + youn3)/3.;
  67. gcis0 = youn0 /(2.*(1.+0.3)) ;
  68. valnu12 = 0.3 ; comm '=> valnu21 = 0.03 ' ;
  69. valnu23 = 0.3 ; comm '=> valnu32 = 0.03 ' ;
  70. valnu13 = 0.3 ; comm '=> valnu31 = 0.003' ;
  71.  
  72. mod0 = MODE v0 mecanique elastique orthotrope  ;
  73. mat0 = MATE mod0 direction (1. 0. 0.) (0. 1. 0.) parallele
  74.             YG1 youn1 YG2 youn2 YG3 youn3
  75.             NU12 valnu12 NU23 valnu23 NU13 valnu13
  76.             G12 gcis0
  77.             G23 gcis0
  78.             G13 gcis0
  79.             RHO 8.e-6 ;
  80. pres0 = pres mass mod0 0.5 s0 ;
  81.  
  82. rig0 = (rigi mod0 mat0) et (bloq depl scl0) ;
  83.  
  84. dep0 = reso rig0 pres0 ;
  85.  
  86. *---- depouillement
  87. defo0 = defo v0 dep0 1. roug ;
  88. depz0 = exco uz dep0 ;
  89. fle1 = extr depz0 scal (v0 poin proc (0. 0. 0.)) ;
  90. titre 'Deplacement en z et deformee (fleche 1 = 'fle1' mm)';
  91. si graph ;
  92.    trace cach (1.e8 1.e8 1.e8) depz0 defo0 ;
  93. finsi ;
  94.  
  95. *---- cas 2
  96. youn1 =  20000. ;
  97. youn2 = 200000. ;
  98. youn3 =   2000. ;
  99. youn0 = (youn1 + youn2 + youn3)/3.;
  100. gcis0 = youn0 /(2.*(1.+0.3)) ;
  101. valnu12 = 0.03 ; comm '=> valnu21 = 0.3  ' ;
  102. valnu23 = 0.3  ; comm '=> valnu32 = 0.003' ;
  103. valnu13 = 0.3  ; comm '=> valnu31 = 0.03 ' ;
  104.  
  105. mat0 = MATE mod0 direction (0. 1. 0.) (1. 0. 0.) parallele
  106.             YG1 youn1 YG2 youn2 YG3 youn3
  107.             NU12 valnu12 NU23 valnu23 NU13 valnu13
  108.             G12 gcis0
  109.             G23 gcis0
  110.             G13 gcis0
  111.             RHO 8.e-6 ;
  112.  
  113. rig0 = (rigi mod0 mat0) et (bloq depl scl0) ;
  114.  
  115. dep0 = reso rig0 pres0 ;
  116.  
  117. *---- depouillement
  118. defo0 = defo v0 dep0 1. roug ;
  119. depz0 = exco uz dep0 ;
  120. fle2 = extr depz0 scal (v0 poin proc (0. 0. 0.)) ;
  121. titre 'Deplacement en z et deformee (fleche 2 = 'fle2' mm)';
  122. si graph ;
  123.    trace cach (1.e8 1.e8 1.e8) depz0 defo0 ;
  124. finsi ;
  125.  
  126. *---------------- test si c'est ok
  127. ladiff = abs(fle2 - fle1) ;
  128. opti echo 0 ;
  130. mess 'Différence entre les 2 résultats =' ladiff ;
  131. si(ladiff > 1.e-10) ;
  132.    erre 5 ;
  133. finsi ;
  134. fin ;
  135.  
  136.  
  137.  
  138.  
  139.  
  140.  
  141.  

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