Télécharger la plaquette Cast3M


fr en

Télécharger channeldie2.dgibi

Retour à la liste

Numérotation des lignes : 

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

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