Télécharger rigi_ic_2d.dgibi

Retour à la liste

Numérotation des lignes :

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

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