Télécharger endoaxi1.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : endoaxi1.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Endoaxi1.dgibi: Jeux de données *
  5. * ------------------------------------ *
  6. * *
  7. *
  8. complet = faux;
  9. *
  10. *************************************************
  11. * pour calcul complet mettre complet à : vrai;
  12. *************************************************
  13. *******************************************************
  14. * *
  15. * TEST; MATERIAU ELASTO-PLASTIQUE ENDOMMAGEABLE *
  16. * *
  17. * EPROUVETTE EN TRACTION AVEC DEPLACEMENTS IMPOSES. *
  18. * MAILLAGE AXISYMETRIQUE. *
  19. * *
  20. * LES RESULTATS OBTENUS SONT COMPARES A LA SOLUTION *
  21. * ANALYTIQUE OBTENUE PAR ALGORITHME. *
  22. * ON TESTE - LA DEFORMATION PLASTIQUE CUMULEE *
  23. * (composante EPSE dans VARIABLES_INTERNES)*
  24. * - D, LA VARIABLE D'ENDOMMAGEMENT *
  25. * (DANS VARIABLES_INTERNES). *
  26. * - SIGMAzz ( dans CONTRAINTES ). *
  27. * *
  28. * Test effectue par Mr STRUB (cf rapport DMT/91 .558).*
  29. * Test modifie pour evolution 2/06/92. *
  30. *******************************************************
  31.  
  32. GRAPH = 'N';
  33.  
  34. OPTI ECHO 1 DIME 2 ELEM QUA8 MODE AXIS;
  35. opti trac bens;
  36. *------------------------------------------------------
  37. P1 = 0 0; P2 = 3E-3 0;
  38. P3 = 3E-3 30E-3; P4 = 0 30E-3;
  39. *
  40. L1 = P1 P2 DROIT 1 ;
  41. L2 = P2 P3 DROIT 3 ;
  42. L3 = P3 P4 DROIT 1 ;
  43. L4 = P4 P1 DROIT 3 ;
  44. *
  45. EPROU = L1 L2 L3 L4 DALLER PLAN ;
  46.  
  47. SI ( NEG GRAPH 'N');
  48. trac eprou;
  49. FINSI;
  50.  
  51. CL1 = BLOQ L1 UZ ;
  52. CL2 = BLOQ L3 UZ ;
  53. CL3 = BLOQ L4 UR ;
  54. CL = CL1 ET CL2 ET CL3 ;
  55. DEP = DEPI CL2 1E-3 ;
  56. *
  57. MO = MODE EPROU 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE'
  58. 'ENDOMMAGEABLE';
  59. *
  60. PROGEPS = PROG 0. .47285E-2 .76066E-2 .89226E-2
  61. 1.0533E-2 1.4011E-2 1.9096E-2 2.6081E-2 3.4654E-2
  62. 4.5828E-2 6.0813E-2 7.2835E-2 8.3328E-2 9.5526E-2
  63. 10.691E-2 11.808E-2 12.784E-2 13.891E-2 14.936E-2
  64. 15.700E-2 16.746E-2 17.720E-2 18.946E-2 19.860E-2
  65. 20.703E-2 21.607E-2 22.802E-2 23.685E-2 24.919E-2
  66. 26.083E-2 ;
  67. PROGSIG = PROG 0. 350.E06 370.59E06 379.18E06
  68. 387.37E06 400.52E06 414.22E06 427.94E06 440.69E06
  69. 453.63E06 467.26E06 476.30E06 483.24E06 490.46E06
  70. 496.56E06 502.06E06 506.53E06 511.28E06 515.51E06
  71. 518.45E06 522.29E06 525.71E06 529.81E06 532.73E06
  72. 535.33E06 538.03E06 541.47E06 543.92E06 547.23E06
  73. 550.24E06 ;
  74. CTRAC = EVOL MANU EPS PROGEPS TRAC PROGSIG ;
  75.  
  76. SI ( NEG GRAPH 'N');
  77. DESS CTRAC;
  78. FINSI;
  79.  
  80. MAT = MATE MO YOUN 74020.E06 NU 0.3 RHO 7800. ALPH 0. TALP 0. TREF 0.
  81. TRAC CTRAC EPSD 0.03 DC 0.23 EPSR 0.25 ;
  82. *
  83. * LISTE DES PAS DE CALCUL :
  84. * 0.13 MM INCREMENT ELASTIQUE
  85. * 0.16 MM INCREMENT PLASTIQUE
  86. si complet;
  87. LIS = PROG 0.13 0.16
  88. 0.5 1. 1.25 1.5 2. 2.5 3. 3.5 4. 4.5
  89. 5. 5.5 6. 6.5 7. 7.5 8.;
  90. sinon;
  91. LIS = prog 0.13 0.16 0.5 ;
  92. finsi;
  93. *
  94. LI1 = PROG 0. 30. ; LI2 = PROG 0. 30. ;
  95. EV = EVOL MANU T LI1 DEP(T) LI2 ;
  96. CHA = CHAR 'DIMP' DEP EV ;
  97. TAB = TABLE ;
  98.  
  99. TAB.'BLOCAGES_MECANIQUES' = CL;
  100. TAB.'CARACTERISTIQUES' = MAT;
  101. TAB.'CHARGEMENT' = CHA;
  102. TAB.'MODELE' = MO;
  103. TAB.'TEMPS_CALCULES' = LIS;
  104. TAB.'PRECISION' = 1E-08;
  105. *
  106. PASAPAS TAB ;
  107.  
  108. *******************************************************
  109. *************** TEST DE FONCTIONNEMENT ***********
  110. *******************************************************
  111.  
  112. dpc0 = extr TAB.VARIABLES_INTERNES.2 EPSE 1 1 1;
  113. dpc1 = extr TAB.VARIABLES_INTERNES.3 EPSE 1 1 1;
  114. dpc0t = 5.46e-4;dpc1t =1.12694E-02;
  115. si complet;
  116. dpc2 = extr TAB.VARIABLES_INTERNES.7 EPSE 1 1 1;
  117. dpc4 = extr TAB.VARIABLES_INTERNES.11 EPSE 1 1 1;
  118. dpc6 = extr TAB.VARIABLES_INTERNES.15 EPSE 1 1 1;
  119. dpc75 =extr TAB.VARIABLES_INTERNES.18 EPSE 1 1 1;
  120.  
  121. dpc2t = 6.03e-2; dpc4t = 12.65e-2; dpc6t = 19.29e-2;
  122. dpc75t = 24.27e-2;
  123. finsi;
  124. D0 = extr TAB.VARIABLES_INTERNES.2 'VHWD' 1 1 1;
  125. D1 = extr TAB.VARIABLES_INTERNES.3 'VHWD' 1 1 1;
  126. D0T = 0.; D1T= 0.0;
  127. si complet;
  128. D2 = extr TAB.VARIABLES_INTERNES.7 'VHWD' 1 1 1;
  129. D4 = extr TAB.VARIABLES_INTERNES.11 'VHWD' 1 1 1;
  130. D6 = extr TAB.VARIABLES_INTERNES.15 'VHWD' 1 1 1;
  131. D75 = extr TAB.VARIABLES_INTERNES.18 'VHWD' 1 1 1;
  132.  
  133. D2T= 0.0257; D4T= 0.0901; D6T= 0.1622; D75T= 0.2199;
  134. finsi;
  135.  
  136. SZZ0 = EXTR TAB.CONTRAINTES.2 SMZZ 1 1 1;
  137. SZZ1 = EXTR TAB.CONTRAINTES.3 SMZZ 1 1 1;
  138. SZZ0T = 354.3E6; SZZ1T = 4.07984E+08;
  139. si complet;
  140. SZZ2 = EXTR TAB.CONTRAINTES.7 SMZZ 1 1 1;
  141. SZZ4 = EXTR TAB.CONTRAINTES.11 SMZZ 1 1 1;
  142. SZZ6 = EXTR TAB.CONTRAINTES.15 SMZZ 1 1 1;
  143. SZZ75 = EXTR TAB.CONTRAINTES.18 SMZZ 1 1 1;
  144.  
  145. SZZ2T = 459.3E6; SZZ4T = 461.4E6; SZZ6T = 443.3E6;
  146. SZZ75T = 422.4E6;
  147. finsi;
  148.  
  149. EDPC0 = ABS((DPC0-DPC0T)/DPC0T);
  150. EDPC1 = ABS((DPC1-DPC1T)/DPC1T);
  151. si complet;
  152. EDPC2 = ABS((DPC2-DPC2T)/DPC2T);
  153. EDPC4 = ABS((DPC4-DPC4T)/DPC4T);
  154. EDPC6 = ABS((DPC6-DPC6T)/DPC6T);
  155. EDPC75 = ABS((DPC75-DPC75T)/DPC75T);
  156. finsi;
  157. ED0 = ABS((D0-D0T)/(D0T+0.001));
  158. ED1 = ABS((D1-D1T)/(D1T+0.001));
  159. si complet;
  160. ED2 = ABS((D2-D2T)/D2T);
  161. ED4 = ABS((D4-D4T)/D4T);
  162. ED6 = ABS((D6-D6T)/D6T);
  163. ED75 = ABS((D75-D75T)/D75T);
  164. finsi;
  165.  
  166. ESZZ0 = ABS((SZZ0-SZZ0T)/SZZ0T);
  167. ESZZ1 = ABS((SZZ1-SZZ1T)/SZZ1T);
  168. si complet;
  169. ESZZ2 = ABS((SZZ2-SZZ2T)/SZZ2T);
  170. ESZZ4 = ABS((SZZ4-SZZ4T)/SZZ4T);
  171. ESZZ6 = ABS((SZZ6-SZZ6T)/SZZ6T);
  172. ESZZ75 = ABS((SZZ75-SZZ75T)/SZZ75T);
  173.  
  174. finsi;
  175. si complet;
  176. LERR = PROG EDPC0 EDPC1 EDPC2 EDPC4 EDPC6 EDPC75 ED0
  177. ED1 ED2 ED4 ED6 ED75
  178. ESZZ0 ESZZ1 ESZZ2 ESZZ4 ESZZ6 ESZZ75;
  179. sinon;
  180. lerr = PROG EDPC0 EDPC1 ED0 ED1 ESZZ0 ESZZ1;
  181. finsi;
  182. list LERR;
  183.  
  184. ermax = maxi LERR;
  185. LIST ermax;
  186. temps;
  187. SI ( ERMAX <EG 0.05 );
  188. ERRE 0;
  189. SINON;
  190. ERRE 5;
  191. FINSI;
  192.  
  193. FIN;
  194. FIN;
  195.  
  196.  
  197.  
  198.  
  199.  
  200.  
  201.  
  202.  
  203.  
  204.  
  205.  
  206.  
  207.  
  208.  
  209.  
  210.  
  211.  
  212.  
  213.  

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