Télécharger plas15.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : plas15.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * *
  5. GRAPH='N';
  6. SAUT PAGE ;
  7. *
  8. *******************************************************
  9. * TEST PLAS15
  10. *
  11. * Essai de compression simple d'un cube en beton
  12. * comportement élasto-plastique modèle OTTOSEN
  13. *
  14. * Un parallelogramme est soumis à un déplacement imposé
  15. * sur une de ses faces et fixé de manière isostatique.
  16. * On compare la courbe de réponse uniaxiale à la
  17. * solution exacte.
  18. *
  19. *******************************************************
  20.  
  21. TITRE 'PARALLELOGRAMME EN COMPRESSION UNIAXIALE ' ;
  22. OPTION ECHO 0 ;
  23. OPTION DIME 2 ELEM QUA4 MODE PLAN CONT ;
  24. *------------------------------------------------------
  25. * DEFINITION DE LA GEOMETRIE
  26. *------------------------------------------------------
  27. A1= 0. 0. ;
  28. B1= 1. 0. ;
  29. C1= 1. 1. ;
  30. D1= 0. 1. ;
  31. *
  32. L1 = A1 DROIT 1 B1 ;
  33. L2 = B1 DROIT 1 C1 ;
  34. L3 = C1 DROIT 1 D1 ;
  35. L4 = D1 DROIT 1 A1 ;
  36. LL = L1 ET L2 ET L3 ET L4 ;
  37. SU= SURF LL PLAN ;
  38. *------------------------------------------------------
  39. * DEFINITION DES CONDITIONS AUX LIMITES
  40. * ET DES DEPLACEMENTS IMPOSES
  41. *------------------------------------------------------
  42. CLX1 = BLOQ L4 UX;
  43. CLX2 = BLOQ L2 UX;
  44. CLY1 = BLOQ L1 UY;
  45. CLY2 = BLOQ L3 UY;
  46. CL1= CLX1 ET CLY1 ;
  47. CL2 = CLX2 ;
  48. DEP = DEPI CL2 -3.5E-4 ;
  49. LI2 = PROG 0. PAS 0.1 10. ;
  50. LI1 = LI2;
  51. EV = EVOL MANU T LI1 F(T) LI2 ;
  52.  
  53. CL = CL1 ET CL2;
  54. MO= MODE SU MECANIQUE ELASTIQUE PLASTIQUE 'OTTOSEN';
  55. TBET = TAILLE MO;
  56.  
  57. *
  58. *------------------------------------------------------
  59. * DEFINITION DES PROPRIETES MATERIELLES
  60. *------------------------------------------------------
  61. YG = 40000. ;
  62. RC = 10. ;
  63. MAT = ( MATE MO YOUN YG NU 0.21
  64. 'GFTR' 1.5D-4 'BTR' 0.25
  65. 'GS' ( 1.8/30.) 'LTR' 1. 'LCS' RC ) ET TBET ;
  66. *------------------------------------------
  67. * APPEL A PASAPAS
  68. * AU TEMPS T , LE DEPLACEMENT VAUDRA F(T)*D
  69. *------------------------------------------
  70. CHA1 = CHAR DIMP DEP EV ;
  71. LIS1 = PROG 0. PAS 0.05 4. ;
  72. TAB1 = TABLE ;
  73. TAB1.'CARACTERISTIQUES' = MAT;
  74. TAB1.'CHARGEMENT' = CHA1 ;
  75. TAB1.'MODELE' = MO;
  76. TAB1.'TEMPS_CALCULES' = LIS1;
  77. TAB1.'BLOCAGES_MECANIQUES' = CL;
  78. TAB1.FTOL=0.000000001;
  79. PASAPAS TAB1;
  80. *------------------------------------------
  81. * DEPOUILLEMENT ET
  82. * CALCUL DE LA SOLUTION ANALYTIQUE
  83. *
  84. *------------------------------------------
  85. *
  86. PGX = PROG 0. ;
  87. PGY = PROG 0. ;
  88. PGT = PROG 0. ;
  89. KE=( 4.*RC)/(3.*YG) ;
  90. KU=5.*KE ;
  91. DK = KU - KE ;
  92. EMAX=(RC/YG) + KE ;
  93. EULT= KU;
  94.  
  95. TT = TAB1 . CONTRAINTES ;
  96. NN = (DIME TT) - 1 ;
  97. *
  98. REPE BOU NN;
  99. SX = EXTR ( TAB1 . CONTRAINTES . (&BOU) )
  100. SMXX 1 1 1 ;
  101. SX = ABS ( SX );
  102. EPS = EPSI MO ( TAB1 . DEPLACEMENTS . (&BOU) ) ;
  103. E = ABS (EXTR EPS EPXX 1 1 1) ; D
  104. E = ABS E ;
  105. PGX = PGX ET ( PROG E ) ;
  106. PGY = PGY ET ( PROG SX ) ;
  107. SI ( E <EG (RC /(3. * YG)));
  108. S = YG * E ;
  109. SINON;
  110. SI ( E <EG EMAX) ;
  111. AA = 2. / (( KE * YG) ** 2 ) ;
  112. BB = (3./RC)+ ( 4. / ( KE * YG ) ) ;
  113. BB = BB - (( 4. * E ) / ((KE * YG) * KE ));
  114. CC =(( 2. * (E* E))/ (KE*KE)) - 1. ;
  115. CC = CC - (( 4. * E) /KE ) ;
  116. DD = (BB * BB ) - (( 4. * AA) * CC) ;
  117. S = ((-1.* BB) + (DD ** 0.5)) / ( 2. * AA ) ;
  118. SINON;
  119. SI ( E <EG EULT) ;
  120. DDE = E - KE ;
  121. AA = 1. / (( DK * YG) ** 2 ) ;
  122. BB = (1./RC) - (( 2. * DDE) / (( DK*DK) * YG ) ) ;
  123. CC =((DDE*DDE)/ (DK*DK)) - 1. ;
  124. DD = (BB * BB ) - (( 4. * AA) * CC) ;
  125. S = ((-1.* BB) + (DD ** 0.5)) / ( 2. * AA ) ;
  126. SINON;
  127. S=0.;
  128. FINSI;
  129. FINSI;
  130. FINSI;
  131. FINSI;
  132.  
  133. PGT = PGT ET ( PROG S ) ;
  134. FIN BOU ;
  135. EV1= EVOL ROUG MANU 'EPX' PGX 'SIGX' PGY ;
  136. EV2= EVOL TURQ MANU 'EPX' PGX 'STHE' PGT ;
  137. SI(NEG GRAPH 'N');
  138. DESS ( EV1 ET EV2 ) ;
  139. FINSI ;
  140. *
  141. TEMPS ;
  142. * CODE FONCTIONNEMENT
  143. ZZ =( MAXI ABS ( PGY - PGT ) ) / RC ;
  144. SI (ZZ <EG 2.E-2) ;
  145. MESS ' ';
  146. MESS ' ';
  147. MESS ' ';
  148. MESS ' COMPARAISON RESULTATS THEORIQUES ET CALCULES';
  149. MESS ' ';
  150. MESS ' JUSQU AU TEMPS 4. : PRECISION < A 2E-2 PRES';
  151. SINO;
  152. ERRE 5;
  153. FINSI;
  154. FIN;
  155.  
  156.  
  157.  
  158.  
  159.  
  160.  

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