Télécharger la plaquette Cast3M


fr en

Télécharger ottovari_traction.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : ottovari_traction.dgibi
  2. ************************************************************************
  3. *                 Test du modele OTTOVARI en traction                  *
  4. ************************************************************************
  5.  
  6. ITRAC = FAUX ;
  7.  
  8. TITR 'Test en traction simple' ;
  9. OPTI 'DIME' 3 'ELEM' 'CUB8' 'MODE' 'TRID' ;
  10.  
  11. *------------------------------------------------------
  12. *           DEFINITION DE LA GEOMETRIE
  13. *------------------------------------------------------
  14. A1= 0.   0.   0.;
  15. B1= 1.5  0.   0.;
  16. C1= 1.5  1.5  0.;
  17. D1= 0.   1.5  0.;
  18. A2= 0.   0.   1.;
  19. B2= 1.5  0.   1.;
  20. C2= 1.5  1.5  1.;
  21. D2= 0.   1.5  1.;
  22. *
  23. L1 = A1 DROIT 1 B1  ;
  24. L2 = B1 DROIT 1 C1  ;
  25. L3 = C1 DROIT 1 D1  ;
  26. L4 = D1 DROIT 1 A1  ;
  27. LL = L1 ET L2 ET L3 ET L4 ;
  28. SU1= SURF LL 'PLAN' ;
  29. M1 = A2 DROIT 1 B2  ;
  30. M2 = B2 DROIT 1 C2  ;
  31. M3 = C2 DROIT 1 D2  ;
  32. M4 = D2 DROIT 1 A2  ;
  33. MM = M1 ET M2 ET M3 ET M4 ;
  34. SU2= SURF MM 'PLAN' ;
  35. SU = SU1 VOLU 1 SU2 ;
  36.  
  37. *------------------------------------------------------
  38. *           DEFINITION DES CONDITIONS AUX LIMITES
  39. *                ET DES DEPLACEMENTS IMPOSES
  40. *------------------------------------------------------
  41. CL4 = (BLOQ L4 'UX') ET (BLOQ M4 'UX') ;
  42. CL2 = (BLOQ L2 'UX') ET (BLOQ M2 'UX') ;
  43. CL1 = (BLOQ L1 'UY') ET (BLOQ M1 'UY') ;
  44. CL5 = BLOQ SU1 'UZ' ;
  45. CL  = CL1 ET CL2 ET CL4 ET CL5 ;
  46.  
  47. DEP = DEPI CL2 1.5E-4 ;
  48.  
  49. *------------------------------------------------------
  50. *           DEFINITION DES PROPRIETES MATERIELLES
  51. *------------------------------------------------------
  52. MO   = MODE SU 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' 'PLASTIQUE' 'OTTOVARI' 'CONS' 'BETON' ;
  53. TBET = TAILLE MO ;
  54.  
  55. YG  = 40000. ;
  56. RTR = 1.1 ;
  57. GF  = (RTR * 3.E-4) * 0.5 ;
  58. MA = MATE MO 'YOUN' YG 'NU' 0.21 'LTR' RTR 'GFTR' GF 'BTR' 0.2 'LCS' 20. ;
  59. MA = MA ET TBET ;
  60.  
  61. LI1 = PROG 0. 'PAS'  0.5 2.  'PAS'  0.5 18.;
  62. LI2 = PROG 0. 'PAS'  0.5 2.  'PAS' -0.5  1.
  63.               'PAS'  0.5 1.5 'PAS' -0.5  0.5 'PAS'  0.5  2.5
  64.               'PAS' -0.5 1.  'PAS'  0.5  5.  'PAS' -0.5 -1.;
  65. EV  = EVOL 'MANU' 'T' LI1 'F(T)' LI2 ;
  66.  
  67. *------------------------------------------
  68. *            APPEL A PASAPAS
  69. *------------------------------------------
  70. CHA1 = CHAR 'DIMP' DEP EV ;
  71. LIS1 = PROG 0. 'PAS' 0.1 18. ;
  72. TAB1 = TABLE ;
  73. TAB1.'MODELE' = MO ;
  74. TAB1.'CARACTERISTIQUES' = MA ;
  75. TAB1.'BLOCAGES_MECANIQUES' = CL ;
  76. TAB1.'CHARGEMENT' = CHA1 ;
  77. TAB1.'TEMPS_CALCULES' = LIS1 ;
  78. TAB1.FTOL = 0.000001 ;
  79. PASAPAS TAB1 ;
  80.  
  81. *------------------------------------------
  82. *          POST TRAITEMENT
  83. *------------------------------------------
  84. PGX = PROG 0. ;
  85. PGY = PROG 0. ;
  86. TT = TAB1 . CONTRAINTES ;
  87. NN = (DIME TT) - 1  ;
  88. REPE BOU NN ;
  89.   FF = EXTR (REDU TAB1 . CONTRAINTES . &BOU MO) 'SMXX' 1 1 1 ;
  90.   EPS = EPSI MO (TAB1 . DEPLACEMENTS . &BOU) ;
  91.   GG = EXTR EPS 'EPXX' 1 1 1 ;
  92.   PGX = PGX ET (PROG GG) ;
  93.   PGY = PGY ET (PROG FF) ;
  94. FIN BOU ;
  95. EV1 = EVOL 'ROUG' 'MANU' 'ABS' PGX 'ORD' PGY ;
  96. SI ITRAC ;
  97.   DESS EV1 'AXES' ;
  98. FINSI ;
  99.  
  100. FIN ;
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  

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