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 ;
  38. graph = faux ;
  39.  
  40. *---- maillage
  41. p1 = -1.e3 -1.e3 0. ;
  42. p2 = 1.e3 -1.e3 0. ;
  43. p3 = 1.e3 1.e3 0. ;
  44. p4 = -1.e3 1.e3 0. ;
  45.  
  46. n1 = 4 ;
  47. n2 = 4 ;
  48.  
  49. l1 = p1 d n1 p2 ;
  50. l2 = p2 d n2 p3 ;
  51. l3 = p3 d n1 p4 ;
  52. l4 = p4 d n2 p1 ;
  53.  
  54. s0 = daller l1 l2 l3 l4 plan ;
  55.  
  56. *---- calcul
  57.  
  58. *---- cas 1
  59. youn1 = 200000. ;
  60. youn2 = 20000. ;
  61. youn0 = (youn1 + youn2)/2.;
  62. gcis0 = youn0 /(2.*(1.+0.3)) ;
  63. valnu12 = 0.3 ; comm '=> valnu21 = 0.03' ;
  64.  
  65. mod0 = MODE s0 mecanique elastique orthotrope coq4 ;
  66. mat0 = MATE mod0 direction (1. 0. 0.) parallele
  67. YG1 youn1 YG2 youn2
  68. NU12 valnu12
  69. G12 gcis0
  70. G23 gcis0
  71. G13 gcis0
  72. RHO 8.e-6 EPAI 50. ;
  73.  
  74. MOP = 'MODE' S0 'CHARGEMENT' 'PRESSION' 'COQ4' ;
  75. MAP = 'MATE' MOP 'PRES' -0.5 ;
  76. MODT = MOD0 'ET' MOP ;
  77. MATT = MAT0 'ET' MAP ;
  78. pres0 = 'PRES' MOP MAP ;
  79.  
  80.  
  81. rig0 = (rigi modT matT) et
  82. (bloq depl rota (l1 et l2 et l3 et l4)) ;
  83.  
  84. dep0 = reso rig0 pres0 ;
  85.  
  86. *---- depouillement
  87. defo0 = defo s0 dep0 1. roug ;
  88. depz0 = exco uz dep0 ;
  89. fle1 = extr depz0 scal (s0 poin proc (0. 0. 0.)) ;
  90. mess 'fleche = ' fle1 ;
  91. titre 'Deplacement en z et deformee (fleche 1 = 'fle1' mm)';
  92. si graph ;
  93. trace cach (1.e8 1.e8 1.e8) depz0 defo0 ;
  94. finsi ;
  95.  
  96. *---- cas 2
  97. youn1 = 20000. ;
  98. youn2 = 200000. ;
  99. youn0 = (youn1 + youn2)/2.;
  100. gcis0 = youn0 /(2.*(1.+0.3)) ;
  101. valnu12 = 0.03 ; comm '=> valnu21 = 0.3' ;
  102.  
  103. mat0 = MATE mod0 direction (1. 0. 0.) parallele
  104. YG1 youn1 YG2 youn2
  105. NU12 valnu12
  106. G12 gcis0
  107. G23 gcis0
  108. G13 gcis0
  109. RHO 8.e-6 EPAI 50. ;
  110. MODT = MOD0 'ET' MOP ;
  111. MATT = MAT0 'ET' MAP ;
  112.  
  113. rig0 = (rigi modT matT) et
  114. (bloq depl rota (l1 et l2 et l3 et l4)) ;
  115.  
  116. dep0 = reso rig0 pres0 ;
  117.  
  118. *---- depouillement
  119. defo0 = defo s0 dep0 1. roug ;
  120. depz0 = exco uz dep0 ;
  121. fle2 = extr depz0 scal (s0 poin proc (0. 0. 0.)) ;
  122. mess 'fleche = ' fle2 ;
  123. titre 'Deplacement en z et deformee (fleche 2 = 'fle2' mm)';
  124. si graph ;
  125. trace cach (1.e8 1.e8 1.e8) depz0 defo0 ;
  126. finsi ;
  127.  
  128. *---- cas 3
  129. youn1 = 200000. ;
  130. youn2 = 20000. ;
  131. youn0 = (youn1 + youn2)/2.;
  132. gcis0 = youn0 /(2.*(1.+0.3)) ;
  133. valnu12 = 0.3 ; comm '=> valnu21 = 0.03' ;
  134.  
  135. mat0 = MATE mod0 direction (0. 1. 0.) parallele
  136. YG1 youn1 YG2 youn2
  137. NU12 valnu12
  138. G12 gcis0
  139. G23 gcis0
  140. G13 gcis0
  141. RHO 8.e-6 EPAI 50. ;
  142. MODT = MOD0 'ET' MOP ;
  143. MATT = MAT0 'ET' MAP ;
  144. rig0 = (rigi modT matT) et
  145. (bloq depl rota (l1 et l2 et l3 et l4)) ;
  146.  
  147. dep0 = reso rig0 pres0 ;
  148.  
  149. *---- depouillement
  150. defo0 = defo s0 dep0 1. roug ;
  151. depz0 = exco uz dep0 ;
  152. fle3 = extr depz0 scal (s0 poin proc (0. 0. 0.)) ;
  153. mess 'fleche = ' fle3 ;
  154. titre 'Deplacement en z et deformee (fleche 3 = 'fle3' mm)';
  155. si graph ;
  156. trace cach (1.e8 1.e8 1.e8) depz0 defo0 ;
  157. finsi ;
  158.  
  159. *---- cas 4
  160. youn1 = 20000. ;
  161. youn2 = 200000. ;
  162. youn0 = (youn1 + youn2)/2.;
  163. gcis0 = youn0 /(2.*(1.+0.3)) ;
  164. valnu12 = 0.03 ; comm '=> valnu21 = 0.3' ;
  165.  
  166. mat0 = MATE mod0 direction (0. 1. 0.) parallele
  167. YG1 youn1 YG2 youn2
  168. NU12 valnu12
  169. G12 gcis0
  170. G23 gcis0
  171. G13 gcis0
  172. RHO 8.e-6 EPAI 50. ;
  173. MODT = MOD0 'ET' MOP ;
  174. MATT = MAT0 'ET' MAP ;
  175. rig0 = (rigi modT matT) et
  176. (bloq depl rota (l1 et l2 et l3 et l4)) ;
  177.  
  178. dep0 = reso rig0 pres0 ;
  179.  
  180. *---- depouillement
  181. defo0 = defo s0 dep0 1. roug ;
  182. depz0 = exco uz dep0 ;
  183. fle4 = extr depz0 scal (s0 poin proc (0. 0. 0.)) ;
  184. mess 'fleche = ' fle4 ;
  185. titre 'Deplacement en z et deformee (fleche 4 = 'fle4' mm)';
  186. si graph ;
  187. trace cach (1.e8 1.e8 1.e8) depz0 defo0 ;
  188. finsi ;
  189.  
  190. *--------- test si c'est ok
  191. opti echo 0 ;
  192. diff12 = abs(fle2 - fle1) ;
  193. mess 'Différence entre les résultats 1 et 2 =' diff12 ;
  194. diff23 = abs(fle3 - fle2) ;
  195. mess 'Différence entre les résultats 2 et 3 =' diff23 ;
  196. diff34 = abs(fle4 - fle3) ;
  197. mess 'Différence entre les résultats 3 et 4 =' diff34 ;
  198. epsilo = 1.e-10 ;
  199. testok = (diff12 < epsilo) ET (diff23 < epsilo) ET (diff34 < epsilo) ;
  200. si (non testok) ;
  201. erre 5 ;
  202. finsi ;
  203.  
  204. fin ;
  205.  
  206.  
  207.  
  208.  
  209.  
  210.  
  211.  
  212.  
  213.  

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