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

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