Télécharger la plaquette Cast3M


fr en

Télécharger rigi_ic_2d.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : rigi_ic_2d.dgibi
  2. ************************************************************************
  3. *            Petit test sur les elements 2D incompressibles            *
  4. *                      (ICT3, ICQ4, ICT6, ICQ8)                        *
  5. * Verification du calcul de la rigidite                                *
  6. * Calculs en 2D : deformations et contraintes planes, axisymmetrie     *
  7. * Sur un essai de traction simple, en effort impose, on verifie qu'il  *
  8. * est bien retrouve a l'encastrement                                   *
  9. ************************************************************************
  10.  
  11.  
  12. ** Options generales
  13. OPTI 'DIME' 2 'ECHO' -1 ;
  15. leg = MOTS 'TRI3' 'QUA4' 'TRI6' 'QUA8' ;
  16. lef = MOTS 'ICT3' 'ICQ4' 'ICT6' 'ICQ8' ;
  17. yo = 200.E9 ;
  18. nu = 0.3 ;
  19. pri = 1.E3 ;
  20. errmax = 0. ;
  21. tol1 = 1.E-10 ;
  22.  
  23.  
  24. ** Points du maillage
  25. p1 = 0. 0. ;
  26. p2 = 1. 0. ;
  27. p3 = 1. 2. ;
  28. p4 = 0. 2. ;
  29.  
  30.  
  31. ** Valeurs theoriques de l'effort de reaction (F = pression * section)
  32. fdp = pri * 1. ;
  33. fcp = pri * 1. ;
  34. fax = pri * pi * 1. ;
  35.  
  36.  
  37. ** Boucle sur les types d'elements finis
  38. n1 = DIME lef ;
  39. REPE b1 n1 ;
  40.   typeg = EXTR leg &b1 ;
  41.   typef = EXTR lef &b1 ;
  42.   MESS ;
  43.   MESS 'Element fini ' typef ;
  44.   MESS '-----------------' ;
  45.   MESS 'Hypothese    | F calculee     | F theorique    | Ecart relatif';
  46. * maillage
  47.   OPTI 'ELEM' typeg ;
  48.   lb = DROI 1 p1 p2 ;
  49.   ld = DROI 1 p2 p3 ;
  50.   lh = DROI 1 p3 p4 ;
  51.   lg = DROI 1 p4 p1 ;
  52.   s1 = DALL lb ld lh lg ;
  53. * boucle sur les hypotheses 2D
  54.   lf = PROG ;
  55.   REPE b2 3 ;
  56. *   matrices de blocage
  57.     SI (EGA &b2 1) ;
  58.       moh = 'def.  planes ' ;
  59.       fth = fdp ;
  60.       OPTI 'MODE' 'PLAN' 'DEFO' ;
  61.       bl1 = BLOQ 'UY' lb ;
  62.       bl2 = BLOQ 'UX' lg ;
  63.     FINSI ;
  64.     SI (EGA &b2 2) ;
  65.       moh = 'cont. planes ' ;
  66.       fth = fcp ;
  67.       OPTI 'MODE' 'PLAN' 'CONT' ;
  68.       bl1 = BLOQ 'UY' lb ;
  69.       bl2 = BLOQ 'UX' lg ;
  70.     FINSI ;
  71.     SI (EGA &b2 3) ;
  72.       moh = 'axisymetrie  ' ;
  73.       fth = fax ;
  74.       OPTI 'MODE' 'AXIS' ;
  75.       bl1 = BLOQ 'UZ' lb ;
  76.       bl2 = BLOQ 'UR' lg ;
  77.     FINSI ;
  78.     bl = bl1 ET bl2 ;
  79. *   matrice de rigidite
  80.     mo = MODE s1 'MECANIQUE' typef ;
  81.     ma = MATE mo 'YOUN' yo 'NU' nu ;
  82.     r = RIGI mo ma ;
  83. *   second membre
  84.     f = PRES 'MASS' mo lh (-1. * pri) ;
  85. *   resolution et calcul de l'effort de reaction a la base
  86.     u = RESO (r ET bl) f ;
  87.     f1 = MAXI 'ABS' (RESU (REAC u bl1)) ;
  88. *   calcul de l'ecart relatif
  89.     SI (EGA &b2 1) ;
  90.       err1 = ABS ((f1 - fdp) / fdp) ;
  91.     FINSI ;
  92.     SI (EGA &b2 2) ;
  93.       err1 = ABS ((f1 - fcp) / fcp) ;
  94.     FINSI ;
  95.     SI (EGA &b2 3) ;
  96.       err1 = ABS ((f1 - fax) / fax) ;
  97.     FINSI ;
  98.     MESS moh '|' f1 '|' fth '|' err1 ;
  99.     errmax = MAXI (PROG errmax err1) ;
  100.   FIN b2 ;
  101. FIN b1 ;
  102.  
  103.  
  104. ** Fin du cas test
  105. SAUT 1 'LIGN' ;
  106. SI (errmax > tol1) ;
  107.   MESS 'ECHEC  DU CAS TEST' ;
  108.   MESS 'Ecart relatif max. = ' errmax ;
  109.   ERRE 5 ;
  110. SINON ;
  111.   MESS 'SUCCES DU CAS TEST' ;
  112. FINSI ;
  113.  
  114.  
  115. FIN ;
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  

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