Télécharger la plaquette Cast3M


fr en

Télécharger test_jointsoft.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  test_jointsoft.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ********************************
  5. *
  6. * Cas test joint_soft
  7. *
  8. *  Ancrage d'une barre d'acier 
  9. *   Essai push-pull
  10. *
  11. *  D. COMBESCURE - EMSI 1999
  12. *
  13. ********************************
  14. opti dime 2 elem seg2;
  15. *
  16. ncot = 1;
  17. *
  18. p0 = 0.  0.;
  19. p1 = 0.  1.;
  20. bar1 = p0 d ncot p1;
  21. *
  22. opti elem qua4;
  23. p01 = 0. 0.;
  24. p11 = 0. 1.;
  25. p02 = 1. 0.;
  26. p12 = 1. 1.;
  27. *
  28. beto = (tran (d ncot p01 p11) ncot p02);
  29. ligsup = d ncot p11 p12;
  30. ligext = d ncot p01 p11;
  31. elim 0.001 (beto et p01 et p11 et p02 et p12 et ligsup et ligext);
  32. *
  33. jointot = racc 1.d-4 ligext bar1 ;
  34. *
  35. * Acier
  36. *
  37. younbar = 200000.;
  38. *
  39. nubar = 0.;
  40. *
  41. Sbar = 3.14159*((10./2000)**2);
  42. *
  43. modbar1 = mode bar1 mecanique elastique barr;
  44. matbar1 = mate modbar1 YOUN younbar NU nubar SECT Sbar
  45.           RHO 0.;
  46. *
  47. *  Joint
  48. *
  49. * Courbe cisaillement
  50. *
  51. KN = (30000./0.05);
  52. FT = 2.1;
  53. *
  54. UXX1 = prog 0.  (6./KN) 1.d-3 3.d-3 10.5d-3 50.d-3;
  55. FXX1 = PROG 0.  6.0    13.0  13.0  6.0     6.0;
  56. * 
  57. tras = evol manu UXX1 FXX1;
  58. *
  59. UXX2 = prog 0. (FT/KN) (5.*(FT/KN)) (10.*(FT/KN)) (20.*(FT/KN));
  60. FXX2 = PROG 0.   FT      FT          (0.5*FT)        (0.5*FT) ;
  61. *
  62. trat = evol manu UXX2 FXX2;
  63. *
  64. PNOR = 0.;
  65. *
  66. *
  67. modj=mode jointot 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  68.                                    JOI2;
  69. matj='MATE' modj 'KS  ' KN 'KN  ' KN ;
  70. *modj=mode jointot 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
  71. *               'PLASTIQUE' 'JOINT_SOFT' JOI2;
  72. *matj='MATE' modj 'KS  ' KN 'KN  ' KN
  73. *                  'PNOR' 0.  'TRAS' TRAS 'TRAT' TRAT;
  74. *
  75. * BETON
  76. *
  77. YOUNBET = 30000.;
  78. nubet = 0.25;
  79. *
  80. modbet = mode beto 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' QUA4;
  81. matbet = mate modbet 'YOUN' younbet 'NU  ' nubet 'RHO ' 0.;
  82. *
  83. modtot = modbar1 et modj et modbet;
  84. mattot = matbar1 et matj et matbet;
  85. *
  86. BL0 = BLOQ DEPL ligsup;
  87. BL1 = BLOQ UY p0;
  88. BL2 = BLOQ UY p1;
  89. DEPP1 = DEPI (BL1 et BL2) (1.);
  90. Protemp = prog 0. 10. 20. 101.;
  91. ProHori = PROG 0. 0.01 -0.01 0.;
  92. Procalc = PROG 0. 0.001 pas 0.1 1. pas 0.25  10.;
  93. CharHori = CHAR DIMP DEPP1 (evol manu Protemp ProHori);
  94. *
  95. TAB1 = TABLE;
  96. ***TAB1 . 'MOVA'          = 'RIEN';
  97. TAB1 . 'TEMPS_CALCULES' = ProCalc ;
  98. *TAB1 . 'CONVERGENCE_FORCEE ' = 'FAUX';
  99. TAB1 .'MAXITERATION'  = 50;
  100. *TAB1.'PRECISION' = 1.D-5;
  101. TAB1 .'CHARGEMENT' = Charhori;
  102. TAB1 .'DYNAMIQUE' = faux;
  103. TAB1 .'MODELE' = modtot;
  104. TAB1 .'CARACTERISTIQUES'= mattot;
  105. TAB1 .'BLOCAGES_MECANIQUES' = BL0 ET BL1 et BL2;
  106. *
  107. PASAPAS TAB1;
  108. *
  109. *
  110. DEPT = TAB1.DEPLACEMENTS;
  111. CONT = TAB1.CONTRAINTES;
  112. TIME = INDEX DEPT;
  113. DIMT = DIME TIME;
  114. *
  115. TIMP = PROG 0.;
  116. DEPX = PROG 0.;
  117. GXYT = prog 0.;
  118. EFFX = prog 0.;
  119. *
  120. EFFXXX = PROG;
  121. REPETER BOU1 (DIMT - 1);
  122.   DEP1 = DEPT.(TIME.(&BOU1 + 1));
  123.   CON1 = CONT.(TIME.(&BOU1 + 1));
  124.   FOR1 = RESU (REAC (BL1 ET BL2) DEP1);
  125.   DEPX = DEPX ET (PROG  (extr DEP1 UY p1));
  126.   EFFX = EFFX ET (PROG (extr FOR1 FY p1));
  127.   TIMP = TIMP ET (PROG (TIME.(&BOU1 + 1)));
  128. FIN BOU1;
  129. *
  130. DXFX = EVOL MANU (DEPX) (EFFX);
  131. TIFX =  EVOL MANU  TIMP (EFFX);
  132. *
  133. *dreal = extr (ipol (prog 0.001) DEPX EFFX) 1;
  134. *dtheo =  extr (ipol (prog 0.001) UXX1 (Sbar*FXX1)) 1;
  135. *prec = abs (dreal - dtheo);
  136. *
  137. *xlim0 = 1.D-7;
  138. *
  139. def1 = defo (extr modtot mail) DEP1;
  140. *
  141. FIN;
  142. *
  143.  
  144.  
  145.  
  146.  
  147.  
  148.  
  149.  
  150.  
  151.  
  152.  
  153.  
  154.  
  155.  

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