Télécharger la plaquette Cast3M


fr en

Télécharger channeldie1.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : channeldie1.dgibi
  2. *----------------------------------------------------------------------*
  3. *                  C H A N N E L D I E 1 . D G I B I                   *
  4. *----------------------------------------------------------------------*
  5. * Objet :                                                              
  6. * -------
  7. *
  8. *    Cas-test de validation des elements BBAR pour les CUB8 et PRI6.
  9. *
  10. *    Calcul elastique de la compression d'un lopin de metal place dans
  11. * une matrice l'empechant de se dilater lateralement (lopin coince dans
  12. * un "canal", d'ou "channel die").
  13. *
  14. *    Validation du calcul par comparaison a la solution analytique.
  15. *
  16. *    Un 2nd chargement en deplacement impose sur la face superieure,
  17. * non uniforme en espace (varaition quadratique), permet de valider
  18. * l'integration dans les elements en verifiant que la pression est 
  19. * est bien sous-integree (champ constant dans l'element), ainsi que
  20. * l'equilibre (F-Bsig).
  21. *
  22. *    Verification du fonctionnement de l'operateur MASSE.
  23. *
  24. *    Verification et validation de l'operateur KSIG en deplacements.
  25. *
  26. * Description :
  27. * -------------
  28. *
  29. * Type de calcul : Mecanique Elastique
  30. * Mode de calcul : 3D
  31. * Type d'element : CUB8, PRI6
  32. * Chargement     : Deplacement impose
  33. *
  34. *----------------------------------------------------------------------*
  35. opti dime 3 mode TRID elem CUB8 ;
  36.  
  37. * Pour affichages, mettre ig1 a VRAI :
  38. ig1      = faux ;
  39.  
  40. *------------------------------ MAILLAGE ------------------------------*
  41.  
  42. * Points :
  43. p1       = 0. 0. 0. ;
  44. p2       = 100. 0. 0. ;
  45. p3       = 100. 100. 0. ;
  46. p4       = 0. 100. 0. ;
  47. p5       = 0. 0. 100. ;
  48.  
  49. * Contour base :
  50. ne1      = 5 ;
  51. li1      = d ne1 p1 p2 ;
  52. li2      = d ne1 p2 p3 ;
  53. li3      = d ne1 p3 p4 ;
  54. li4      = d ne1 p4 p1 ;
  55.  
  56. * Surface et Volume :
  57. sur1     = dall li1 li2 li3 li4 'PLAN' ;
  58. sur2     = (sur1 plus p2) chan tri3 ;
  59. sur1     = sur1 et sur2 ;
  60. elim sur1 1.e-3 ;
  61. vol1     = volu tran sur1 p5 ne1 ;
  62.  
  63. si ig1 ;
  64.   trac cach vol1 titr 'Maillage lopin' ;
  65. fins ;
  66.  
  67. *--------------------- MODELE / CARACTERISTIQUES ----------------------*
  68.  
  69. * Valeurs modules d'elasticite :
  70. ym1      = 1.5e11 ;
  71. nu1      = 0.499 ;
  72.  
  73. * Modele et caracteristiques mecaniques :
  74. mod1     = modele vol1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' 'BBAR' ;
  75. mat1     = mate mod1 'YOUNG' 1.5e11 'NU' nu1 'RHO' 7.6e3 ;
  76.  
  77. * Matrice de raideur :
  78. rig1     = rigi mod1 mat1 ;
  79. mas1     = mass mod1 mat1 ;
  80.  
  81. *----------------------- CONDITIONS AUX LIMITES -----------------------*
  82.  
  83. * Definition des points d'interet :
  84. ptx0     = (vol1 coor 1) poin mini ;
  85. pty0     = (vol1 coor 2) poin mini ;
  86. ptz0     = sur1 ;
  87. ptz1     = vol1 face 2 ;
  88. clx0     = bloq ux ptx0 ;
  89. cly0     = bloq uy pty0 ;
  90. clz0     = bloq uz ptz0 ;
  91. clz1     = bloq uz ptz1 ;
  92. cl0      = clx0 et cly0 et clz0 et clz1 ;
  93.  
  94. * Affichages points CL :
  95. si ig1 ;
  96.   trac ((ptx0 coul roug) et (aret vol1)) titr 'Points Ux = 0' ;
  97.   trac ((pty0 coul roug) et (aret vol1)) titr 'Points Uy = 0' ;
  98.   trac ((ptz0 coul roug) et (ptz1 coul rouge) et (aret vol1)) titr 'Points Uz = 0 (surfaces bloquees par matrice)' ;
  99. fins ;
  100.  
  101. *----------------------------- CHARGEMENT -----------------------------*
  102.  
  103. * Deplacement impose essai channel die :
  104. Uy1      = -0.1 ;
  105. pty1     = (vol1 coor 2) poin maxi ;
  106. cly1     = bloq uy pty1 ;
  107. dcl1     = depi cly1 Uy1 ;
  108.  
  109. * Deplacement impose non homogene :
  110. ptx1     = (vol1 coor 1) poin maxi ;
  111. clx1     = bloq ux ptx1 ;
  112. chy1     = ptz1 coor 2 ;
  113. Uz1      = (((-0.01 * (chy1 - 100.)) ** 2) / -100.) nomc uz ;
  114. dcl2     = depi clz1 Uz1 ;
  115.  
  116. si ig1 ;
  117.   trac ((pty1 coul vert) et (aret vol1)) titr 'Points deplacement Uy impose essai channel die' ;
  118.   trac Uz1 ptz1 titr '2nd chargement : deplacement non homogene' ;
  119. fins ;
  120.  
  121. *---------------------------- DEPLACEMENTS ----------------------------*
  122.  
  123. * Solution analytique au point P3:
  124. UxAna1   = -1. * nu1 / (1. - nu1) * Uy1 ;
  125.  
  126. * Solution Castem au point P3: uycas.
  127. deptot   = reso (rig1 et cl0 et cly1) dcl1 ;
  128. dep2     = reso (rig1 et cl0 et cly1 et clx1) dcl2 ;
  129.  
  130. * deplacement Ux au point P3 :
  131. UxSim1   = extr deptot 'UX' p3 ;
  132.  
  133. * Deformee :
  134. def0     = defo deptot (aret vol1) 0. blan ; 
  135. def1     = defo deptot (aret vol1) 100 roug ; 
  136. mot1     = chai format '(F7.4)' 'deformee (vue de dessus) : Ux =' UxSim1 ;
  137.  
  138. si ig1 ;
  139.   trac (0 0 1.e6) (def0 et def1) titr mot1 ;
  140. fins ;
  141.  
  142. opti oeil (1.e6 -1.e6 0.8e6) ;
  143.  
  144. * Calcul de l'erreur sur le déplacement.
  145. err1     = abs (( UxAna1 - UxSim1 ) / ( UxAna1 )) ;
  146.  
  147. *---------------------------- CONTRAINTES -----------------------------*
  148.  
  149. * Solution analytique: 
  150. syyana1  = ym1 / (1. - (nu1 * nu1)) * uy1 / 100. ;
  151. szzana1  = nu1 * syyana1 ;
  152. prana1   = syyana1 + szzana1 / 3. ;
  153.  
  154. * Solution Castem : pression maximale = maxpres1.
  155. sig1     = sigma mod1 deptot mat1 ;
  156.  
  157. * Validation KSIG avec BBAR :
  158. ksg1     = ksigm mod1 sig1 ;
  159. depto2   = reso (rig1 et ksg1 et cl0 et cly1) dcl1 ;
  160. err6     = maxi abs (deptot - depto2) avec (mots ux uy uz) ;
  161.  
  162. * Affichage :
  163. si ig1 ;
  164.   trac sig1 mod1 titr 'Champ de contrainte' ;
  165. fins ;
  166.  
  167. * Calcul de la pression :
  168. syySim1  = exco 'SMYY' sig1 'P' ;
  169. szzSim1  = exco 'SMZZ' sig1 'P' ;
  170. prSim1   = (syySim1 + szzSim1) / 3. ;
  171. syySim1  = mini syySim1 ;
  172. szzSim1  = mini szzSim1 ;
  173. prSim1   = mini prSim1 ;
  174.  
  175. * Calcul de l'erreur sur les contraintes :
  176. err2     = ((syyana1 - syySim1) / syyana1) + ((szzana1 - szzSim1) / szzana1) + ((prana1 - prSim1) / prana1) ;
  177. err2     = abs err2 ;
  178.  
  179. * Verification integration pression chargement non homogene :
  180. sig2     = sigm mod1 mat1 dep2 ;
  181. p2 I2 I3 = inva mod1 sig2 ;
  182. * Dans prismes :
  183. el1      = (vol1 elem pri6) elem 1 ;
  184. model1   = redu mod1 el1 ;
  185. pmax1 pmin1 = (maxi (redu p2 model1)) (mini (redu p2 model1)) ;
  186. err3     = (maxi abs (pmax1 - pmin1)) / (maxi abs pmax1) ;
  187. * Dans hexaedres :
  188. el2      = (vol1 elem cub8) elem 1 ;
  189. model2   = redu mod1 el2 ;
  190. pmax2 pmin2 = (maxi (redu p2 model2)) (mini (redu p2 model2)) ;
  191. err4     = (maxi abs (pmax2 - pmin2)) / (maxi abs pmax2) ;
  192.  
  193. * Verification equilibre chargement non homogene :
  194. bsg2     = bsig mod1 sig2 ;
  195. rea2     = reac (cl0 et cly1 et clx1) dep2 ;
  196. res2     = rea2 - bsg2 ;
  197. err5     = (maxi abs res2) / (maxi abs rea2) ;
  198.  
  199. * Affichages solution chargement non homogene :
  200. si ig1 ;
  201.   trac (defo vol1 dep2 2.e3) cach titr ' Deformee chargement non homogene' ;
  202.   trac sig2 mod1 titr ' Contraintes chargement non homogene' ;
  203.  
  204. * Affichage valeur champs aux points d'integration :
  205.   list (redu sig2 model1) ;
  206.   list (redu p2   model1) ;
  207.   list (redu sig2 model2) ;
  208.   list (redu p2   model2) ;
  209. fins ;
  210.  
  211. *----------------------------- VALIDATION -----------------------------*
  212.  
  213. * Affichages :
  214. opti echo 0 ; saut 1 lign ;
  215. mess (chai format '(F6.2)' ' > Erreur relative deplacement au point P3 =' err1) ;
  216. mess (chai format '(F6.2)' ' > Erreur relative contraintes et pression =' err2) ;
  217. opti echo 1 ; saut 1 lign ;
  218.  
  219. * Test de validation:
  220. err0     = maxi abs (prog err1 err2 err3 err4 err5 err6) ;
  221. prec0    = (vale prec) ** 0.5 * 10. ;
  222. list err0 ;
  223. list err1 ;
  224. list err2 ;
  225. list err3 ;
  226. list err4 ;
  227. list err5 ;
  228. list err6 ;
  229. list prec0 ;
  230. si ( err0 >eg prec0 ) ;
  231.   erreur 5 ;
  232. finsi ;
  233.  
  234. fin;
  235.  
  236.  
  237.  
  238.  
  239.  
  240.  
  241.  

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