Télécharger testIC20.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier :
  2. ************************************************************************
  3. ************************************************************************
  4. * Test TestIC20.dgibi: Jeux de données *
  5. * ------------------------------------ *
  6. * *
  7. *
  8. * FICHIER GIBIANE POUR TESTER LES ELEMENTS
  9. * INCOMPRESSIBLES volumiques quadratiques
  10. *
  11. * Maillage d'un cube avec les
  12. * éléments IC
  13. *
  14. * On soumet le cube à un champ de contrainte
  15. * hydrostatique (contribution: alpha=1) et à un champ
  16. * de contrainte déviatorique (contribution: beta=0.001)
  17. *
  18. * Les résultats théoriques sont déterminés dans le cas
  19. * incompressible (coefficient de Poisson proche de 0.5)
  20. *
  21. * cf rapport DMT sur ces elements
  22. *
  23. *
  24. *-----MAILLAGE-----------
  25. *
  26. opti dime 2 mode plan defo elem qua8 ;
  27. opti dime 2 mode plan defo elem tri6 ;
  28. *
  29. p1 = 0. 0. ;
  30. p2 = 100. 0. ;
  31. p3 = 100. 100. ;
  32. p4 = 0. 100. ;
  33. *
  34. nm2 = 10;
  35. li1 = d nm2 p1 p2 ;
  36. li2 = d nm2 p2 p3 ;
  37. li3 = d nm2 p3 p4 ;
  38. li4 = d nm2 p4 p1 ;
  39. *
  40. sur1 = dall li1 li2 li3 li4 'PLAN' ;
  41. *
  42. opti dime 3 elem cu20 ;
  43. opti dime 3 elem pr15 ;
  44. e_z = 0. 0. 100. ;
  45. vol1 = volu sur1 tran nm2 e_z ;
  46. *sur2 = changer (face 2 vol1) tri6 ;
  47. sur2 = face 2 vol1;
  48. *sur3 = changer (face 3 vol1) tri6 ;
  49. sur3 = face 3 vol1 ;
  50. svol1 = enve vol1 ;
  51. *svol1 = changer svol1 tri6 ;
  52. *opti dime 3 elem te10 ;
  53. *vol1 = volu svol1 ;
  54. *elim 1.e-4 svol1 sur3;
  55. *elim 1.e-4 svol1 sur2;
  56. *elim 1.e-4 svol1 sur1;
  57.  
  58. ph3 = point vol1 proc (p3 plus e_z) ;
  59. *
  60. *-----MODELE MECANIQUE------------
  61. *
  62. * Valeur du coefficient de Poisson
  63. *
  64. nu1 = 0.499 ;
  65. *
  66. * Courbe de traction :
  67. lep1 = prog 0. 2.e-3 1. ;
  68. lsm1 = prog 0. 420.e6 ((1.-2.e-3)*1.e9 + 420.e6) ;
  69. trac1 = evol manu lep1 lsm1 ;
  70. *dess trac1 titr
  71. *'Courbe de traction (Young = 210 GPa, pente plastique = 1 GPa)' nclk;
  72. *
  73. *mod1 = modele vol1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' IC20 ;
  74. *mod2 = modele vol1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  75. * 'PLASTIQUE' 'ISOTROPE' ic20 ;
  76. mod1 = modele vol1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' bbar ;
  77. mod2 = modele vol1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  78. 'PLASTIQUE' 'ISOTROPE' bbar ;
  79. *mod1 = modele vol1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' IC10 ;
  80. *mod2 = modele vol1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  81. * 'PLASTIQUE' 'ISOTROPE' IC10 ;
  82. mat1 = mate mod1 'YOUNG' 2.e11 'NU' nu1 ;
  83. mat2 = mate mod2 'YOUNG' 210.e9 'NU' nu1 trac trac1 ;
  84. *
  85. *
  86. *-----CONDITIONS AUX LIMITES------------
  87. *
  88. bl1 = bloq sur1 'UZ' ;
  89. bl2 = bloq li1 'UY' ;
  90. bl3 = bloq p1 'UX' 'UY' ;
  91. bl4 = bloq li4 'UX' ;
  92. *
  93. *
  94. *----------------CHARGEMENT-----------------
  95. *
  96. * Coefficients de contribution de chargement
  97. * (alpha-->hydrostatique)
  98. * (beta-->déviatorique)
  99. *
  100. alpha = 1. ;
  101. beta = 1.e-3 ;
  102. *
  103. *-----Chargement déviatorique------------
  104. *
  105. *-----sur li1-----
  106. *
  107. chpx1 = coor 3 sur3 ;
  108. chpresy1 = pres 'MASS' mod1 (chpx1 * (-2.)*1.e5) ;
  109. chpres1 = beta * chpresy1 ;
  110. vecfor1 = vect chpres1 1.E-7 'FX' 'FY' 'FZ' rouge ;
  111. *
  112. *
  113. *-----Chargement hydrostatique------------
  114. *
  115. *-----svol1 ----
  116. *
  117. presy1 = pres 'MASS' mod1 (1.e5) svol1 ;
  118. pres1 = alpha * presy1 ;
  119. vfor1 = vect pres1 5.E-7 'FX' 'FY' 'FZ' vert ;
  120. *
  121. *trac vol1 cach ;
  122. *trac vfor1 svol1 cach ;
  123. *
  124. *---------CALCUL DES DEPLACEMENTS --------------------
  125. *
  126. *
  127. * Solution analytique au point Ph3: uzan.
  128. *
  129. x3 = coor 1 ph3 ;
  130. z3 = coor 3 ph3 ;
  131. mu = 1.e11 / ( 1 + nu1 ) ;
  132. uzan = ( 2 * x3 * z3 ) / ( 2 * mu ) * 1.e5 * beta ;
  133. *
  134. * Solution Castem au point P3: uycas.
  135. rigi1 = rigi mod1 mat1 ;
  136. *
  137. deptot = reso (rigi1 et bl1 et bl2 et bl3 et bl4) pres1 ;
  138. *deptot = reso (rigi1 et bl1 et bl3) pres1 ;
  139. *( pres2 et pres3 et pres4 et pres1 et (beta *chforx3)) ;
  140. *(chpres2 et chpres3 et chpres4 et chpres1 et
  141. uzcas = extr deptot 'UZ' ph3 ;
  142. *list (redu deptot sur2) ;
  143. *
  144. * Calcul de l erreur sur le déplacement.
  145. *
  146. erry = abs (( uzan - uzcas ) / ( uzan )) * 100 ;
  147. *
  148. * Test de précision sur le calcul du déplacement uy
  149. * (<1%).
  150. *
  151. si ( erry >eg 1. ) ;
  152. * erreur 5 ;
  153. mess 'erry' erry ;
  154. finsi ;
  155. *
  156. *
  157. *-----------CALCUL CONTRAINTES-------------------------
  158. *
  159. *
  160. * Solution analytique: pression = alpha*1.E5.
  161. *
  162. * Solution Castem : pression maximale = maxpres1.
  163. *
  164. sig1 = sigma mod1 deptot mat1 ;
  165. sxx = exco 'SMXX' sig1 'P' ;
  166. syy = exco 'SMYY' sig1 'P' ;
  167. szz = exco 'SMZZ' sig1 'P' ;
  168. sp = (sxx + syy + szz) / 3. ;
  169. maxpres1 = maxi sp ;
  170. *
  171. *trac sig1 mod1 ;
  172. * Calcul de l erreur sur la pression.
  173. *
  174. errpress = abs ((maxpres1 + (alpha * 1.E5)) /
  175. (alpha * 1.E5)) * 100. ;
  176. *
  177. * Test de précision sur le calcul de la pression (<1%).
  178. *
  179. mess 'errpress' errpress maxpres1;
  180. si ( errpress >eg 1. ) ;
  181. erreur 5 ;
  182. finsi ;
  183. *
  184. ev_ic = evol manu (prog 0. 1. ) (prog 0. 10.) ;
  185. ch_hy = char meca ev_ic pres1 ;
  186.  
  187. t_ic = table ;
  188. t_ic . modele = mod2 ;
  189. t_ic . caracteristiques = mat2 ;
  190. t_ic . blocages_mecaniques = bl1 et bl2 et bl3 ;
  191. t_ic . chargement = ch_hy ;
  192. t_ic . temps_sauves = prog 1. ;
  193. pasapas t_ic ;
  194. s_ic = t_ic . contraintes . 1 ;
  195. *trac s_ic mod2 ;
  196. *list (redu t_ic . deplacements . 1 sur2) ;
  197. *
  198.  
  199. fin;
  200.  
  201.  
  202.  
  203.  
  204.  
  205.  
  206.  
  207.  
  208.  

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