Télécharger plas9.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : plas9.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Plas9.dgibi: Jeux de données *
  5. * --------------------------------- *
  6. * *
  7. OPTION ECHO 0;
  8. * ajout de option epsilon lineaire pour la precision des test!
  9. OPTION epsilon lineaire;
  10.  
  11. GRAPH='N';
  12. SAUT PAGE ;
  13. *
  14. *******************************************************
  15. * TEST PLAS9
  16. *
  17. * Sortie du domaine élastique et phase plastique
  18. * (comportement élasto-plastique modèle DRUCKER-PRAGER
  19. * à écrouissage négatif).
  20. *
  21. * Un parallelépipède est soumis à un déplacement imposé
  22. * sur une de ses faces et fixé de manière isostatique.
  23. * On compare la courbe de réponse uniaxiale à la
  24. * solution exacte.
  25. *
  26. *******************************************************
  27.  
  28. TITRE 'PARALLELEPIPEDE EN TRACTION UNIAXIALE ' ;
  29. OPTION ECHO 0 ;
  30. OPTION DIME 3 ELEM CUB8 MODE TRID ;
  31. *------------------------------------------------------
  32. * DEFINITION DE LA GEOMETRIE
  33. *------------------------------------------------------
  34. A1= 0. 0. 0.;
  35. B1= 1.5 0. 0.;
  36. C1= 1.5 1.5 0.;
  37. D1= 0. 1.5 0.;
  38. A2= 0. 0. 1.;
  39. B2= 1.5 0. 1.;
  40. C2= 1.5 1.5 1.;
  41. D2= 0. 1.5 1.;
  42. *
  43. L1 = A1 DROIT 1 B1 ;
  44. L2 = B1 DROIT 1 C1 ;
  45. L3 = C1 DROIT 1 D1 ;
  46. L4 = D1 DROIT 1 A1 ;
  47. LL = L1 ET L2 ET L3 ET L4 ;
  48. SU1= SURF LL PLAN ;
  49. M1 = A2 DROIT 1 B2 ;
  50. M2 = B2 DROIT 1 C2 ;
  51. M3 = C2 DROIT 1 D2 ;
  52. M4 = D2 DROIT 1 A2 ;
  53. MM = M1 ET M2 ET M3 ET M4 ;
  54. SU2= SURF MM PLAN ;
  55. SU = SU1 VOLU 1 SU2 ;
  56. SI(NEG GRAPH 'N');
  57. TRAC (888 1000 2000 ) SU ;
  58. FINSI;
  59. *------------------------------------------------------
  60. * DEFINITION DES CONDITIONS AUX LIMITES
  61. * ET DES DEPLACEMENTS IMPOSES
  62. *------------------------------------------------------
  63. CL4 = (BLOQ L4 UX) ET ( BLOQ M4 UX) ;
  64. CL2 = (BLOQ L2 UX) ET ( BLOQ M2 UX) ;
  65. CL1 = (BLOQ L1 UY) ET (BLOQ M1 UY);
  66. CL5 = BLOQ SU1 UZ ;
  67. CL = CL1 ET CL2 ET CL4 ET CL5 ;
  68. *
  69. DEP= DEPI CL2 1.5E-4 ;
  70. *
  71. *------------------------------------------------------
  72. * DEFINITION DES PROPRIETES MATERIELLES
  73. *------------------------------------------------------
  74. MO1= MODE SU MECANIQUE ELASTIQUE PLASTIQUE
  75. DRUCKER_PRAGER CONS 'CONS1' ;
  76. MO2= MODE SU MECANIQUE ELASTIQUE CONS 'CONS2';
  77. MO = MO1 ET MO2 ;
  78. VALK = 2. ;
  79. VALALF =0.441 ;
  80. VALBE = 1.02 ;
  81. VALH = -10000.;
  82. YOYO = 40000. ;
  83. MAT1 = MATE MO1 YOUN YOYO NU 0.21
  84. ALFA VALALF BETA VALBE K VALK
  85. GAMM VALALF DELT VALBE ETA VALALF MU VALBE
  86. KL VALK H VALH ;
  87. MAT2 = MATE MO2 YOUN 2.D-3 NU 0. ;
  88. MAT = MAT1 ET MAT2 ;
  89. LI1 = PROG 0. PAS 1. 1000. ;
  90. LI2 = LI1 ;
  91. EV = EVOL MANU T LI1 F(T) LI2 ;
  92. *------------------------------------------
  93. * APPEL A PASAPAS
  94. * AU TEMPS T , LE DEPLACEMENT VAUDRA F(T)*D
  95. *------------------------------------------
  96. CHA1 = CHAR DIMP DEP EV ;
  97. LIS1 = PROG 0. PAS 0.1 4. ;
  98. TAB1 = TABLE ;
  99. TAB1.'CARACTERISTIQUES' = MAT;
  100. TAB1.'CHARGEMENT' = CHA1;
  101. TAB1.'MODELE' = MO;
  102. TAB1.'TEMPS_CALCULES' = LIS1;
  103. TAB1.'BLOCAGES_MECANIQUES' = CL;
  104. TAB1.FTOL=0.000001;
  105. PASAPAS TAB1;
  106. *------------------------------------------
  107. * DEPOUILLEMENT ET
  108. * CALCUL DE LA SOLUTION ANALYTIQUE
  109. *
  110. * ON VERIFIE EN GROS QUE LE CALCUL EST
  111. * ALLE JUSQU'AU BOUT
  112. *------------------------------------------
  113. *
  114. MU = (((2. * ( VALALF ** 2)) + ( VALBE ** 2) ) ** 0.5)
  115. / ( VALALF + VALBE) ;
  116. VALHMU = VALH * MU ;
  117. PGX = PROG 0. ;
  118. PGY = PROG 0. ;
  119. PGT = PROG 0. ;
  120. TT = TAB1 . CONTRAINTES ;
  121. NN = (DIME TT) - 1 ;
  122. SI (( TAB1 . TEMPS . NN) < 3.99) ;
  123. mess ' le temps est trop petit!!!!';
  124. ERRE 5 ;
  125. FINSI;
  126. *
  127. REPE BOU NN;
  128. *si ( &bou < 5); list ( tab1. contraintes . &bou);finsi;
  129. FF = EXTR ( redu TAB1 . CONTRAINTES . (&BOU) MO1)
  130. SMXX 1 1 1 ;
  131.  
  132. EPS = EPSI MO1 ( TAB1 . DEPLACEMENTS . (&BOU) ) ;
  133. GG = EXTR EPS EPXX 1 1 1 ;
  134. ETH = VALK / ( ( VALALF + VALBE ) * YOYO ) ;
  135. SI ( GG < ETH ) ;
  136. THE = YOYO * GG ;
  137. SINON;
  138. THE= ( VALK + ( VALHMU * GG )) / ( (VALALF + VALBE) +
  139. ( VALHMU / YOYO ) ) ;
  140. SI ( THE < 0. ) ; THE = 0. ; FINSI ;
  141. FINSI ;
  142. PGX = PGX ET ( PROG GG ) ;
  143. PGY = PGY ET ( PROG FF ) ;
  144. PGT = PGT ET ( PROG THE ) ;
  145. FIN BOU ;
  146. EV1= EVOL ROUG MANU ABS PGX ORD PGY ;
  147. EV2= EVOL TURQ MANU ABS PGX ORD PGT ;
  148. EV = EV1 ET EV2;
  149. SI(NEG GRAPH 'N');
  150. DESS EV;
  151. FINSI ;
  152. *
  153. TEMPS ;
  154. * CODE FONCTIONNEMENT
  155. RESI1=(MAXI ABS (PGY - PGT ))/ VALK ;
  156. SI(RESI1 &lt;EG 5E-5);
  157. MESS ' ';
  158. MESS ' ';
  159. MESS ' ';
  160. MESS ' ';
  161.  
  162. MESS ' COMPARAISON RESULTATS THEORIQUES ET CALCULES';
  163. MESS ' ';
  164. MESS ' JUSQU AU TEMPS 4. : PRECISION < A 5E-5 PRES';
  165. SINO;
  166. mess ' re_cal val_the precision ' (MAXI ABS(pgy - pgt)) valk resi1;
  167. ERRE 5;
  168. FINSI;
  169. FIN;
  170.  
  171.  
  172.  
  173.  
  174.  
  175.  
  176.  
  177.  
  178.  
  179.  
  180.  
  181.  
  182.  

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