Télécharger endoaxi2.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : endoaxi2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Endoaxi2.dgibi: Jeux de données *
  5. * ------------------------------------ *
  6. * *
  7.  
  8. *******************************************************
  9. * *
  10. * TEST; MATERIAU ELASTO-PLASTIQUE ENDOMMAGEABLE *
  11. * MATERIAU DEPENDANT DE LA TEMPERATURE : *
  12. * COURBE DE TRACTION *
  13. * COEFFICIENT DE DILATATION *
  14. * THERMIQUE *
  15. * *
  16. * EPROUVETTE EN TRACTION AVEC DEPLACEMENTS IMPOSES *
  17. * + CHARGEMENT THERMIQUE *
  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.  
  29.  
  30. OPTI ECHO 0 DIME 2 ELEM QUA8 MODE AXIS;
  31. *------------------------------------------------------
  32. P1 = 0 0; P2 = 3E-3 0;
  33. P3 = 3E-3 30E-3; P4 = 0 30E-3;
  34. *
  35. L1 = P1 P2 DROIT 1 ;
  36. L2 = P2 P3 DROIT 1 ;
  37. L3 = P3 P4 DROIT 1 ;
  38. L4 = P4 P1 DROIT 1 ;
  39. *
  40. EPROU = L1 L2 L3 L4 DALLER PLAN ;
  41.  
  42. CL1 = BLOQ L1 UZ ;
  43. CL2 = BLOQ L3 UZ ;
  44. CL3 = BLOQ L4 UR ;
  45. CL = CL1 ET CL2 ET CL3 ;
  46. DEP = DEPI CL2 1E-3 ;
  47. *
  48. MO = MODE EPROU 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE'
  49. 'ENDOMMAGEABLE';
  50. *
  51. * module tangent valant 783MPa
  52. *-----------------------------
  53. PROGEPS1 = PROG 0. (350./74020.) 0.3 0.6 ;
  54. PROGSIG1 = PROG 0. 350.E06 581.1976E06 816.0976E06 ;
  55. CTRAC1 = EVOL MANU EPS PROGEPS1 SIG PROGSIG1 ;
  56. *
  57. * module tangent valant 491MPa
  58. *-----------------------------
  59. PROGEPS2 = PROG 0. (350./74020.) 0.3 0.6 ;
  60. PROGSIG2 = PROG 0. 350.E06 494.9783E06 642.2783E06 ;
  61. CTRAC2 = EVOL MANU EPS PROGEPS2 SIG PROGSIG2 ;
  62. *
  63. * module tangent valant 199MPa
  64. *-----------------------------
  65. PROGEPS3 = PROG 0. (350./74020.) 0.3 0.6 ;
  66. PROGSIG3 = PROG 0. 350.E06 408.7590E06 468.4590E06 ;
  67. CTRAC3 = EVOL MANU EPS PROGEPS3 SIG PROGSIG3 ;
  68. *
  69. * module tangent variable
  70. *------------------------
  71. CTRAC = NUAGE 'T'*'FLOTTANT' 'TRAC'*'EVOLUTION'
  72. 0. CTRAC1 100. CTRAC2 200. CTRAC3 ;
  73. *
  74. * coefficient de dilatation thermique vaiable
  75. *--------------------------------------------
  76. *
  77. PROGTEMP = PROG 0. 100. 200. ;
  78. PROGALFA = PROG 0. 1E-3 2E-3 ;
  79. EVALPH = EVOL MANU 'T' PROGTEMP 'ALPH' PROGALFA ;
  80. *
  81. * MAT = MATE MO YOUN 74020.E06 NU 0.3 RHO 7800.
  82. * ALPH 1.E-3 TRAC CTRAC1 EPSD 0. DC 0.23 EPSR 0.25 ;
  83. *
  84. MATvar = MATE MO YOUN 74020.E06 NU 0.3 RHO 7800.
  85. ALPH EVALPH TRAC CTRAC EPSD 0. DC 0.23 EPSR 0.25 ;
  86. *
  87. ************** DEFINITION DU CHARGEMENT ***************
  88. *
  89. LI1 = PROG 0. 30. ; LI2 = PROG 0. 30. ;
  90. EV = EVOL MANU T LI1 DEP(T) LI2 ;
  91. CHA1 = CHAR 'DIMP' DEP EV ;
  92. TEMP0 = MANU CHPO EPROU 1 'T' 0. ;
  93. TEMP200 = MANU CHPO EPROU 1 'T' 200. ;
  94. t0 = 0. ; t20 = 20. ;
  95. TEMP = TABLE; TEMPERA = TABLE;
  96. TEMP.0 = t0 ; TEMP.1 = t20 ;
  97. TEMPERA.0 = TEMP0; TEMPERA.1 = TEMP200;
  98. CHA2 = CHAR 'T' TEMP TEMPERA;
  99. CHA = CHA1 ET CHA2;
  100. *
  101. * LISTE DES PAS DE CALCUL :
  102. * 0.14 MM INCREMENT ELASTIQUE
  103. * 0.15 MM INCREMENT PLASTIQUE
  104. LIS = PROG 0.14 0.15 0.2 0.4 0.7 1 pas .3
  105. 3.1 3.3 pas .2 4.1 ;
  106. *
  107. TAB = TABLE ;
  108. TAB.'PRECISION' = 1E-3;
  109. TAB.'BLOCAGES_MECANIQUES' = CL;
  110. TAB.'CARACTERISTIQUES' = MATvar;
  111. TAB.'CHARGEMENT' = CHA;
  112. TAB.'TEMPS_CALCULES' = LIS;
  113. TAB.'MODELE' = MO;
  114. PASAPAS TAB ;
  115. *
  116. *
  117. *------------------------------------------------------
  118. * CONTROLE DES RESULTATS AVEC DE LA SOLUTION
  119. * DE REFERENCE OBTENUE PAR ALGORITHME
  120. *------------------------------------------------------
  121. *
  122. REF_D = PROG 2.197E-2 3.417E-2 6.337E-2 8.131E-2 ;
  123. REF_P = PROG 3.646E-2 5.443E-2 9.296E-2 1.143E-1 ;
  124. REF_S = (PROG 368.9 376.0 386.4 389.2)*1E6 ;
  125. *
  126. D1 = EXTR (TAB.VARIABLES_INTERNES.7) 'VHWD' 1 1 1;
  127. D2 = EXTR (TAB.VARIABLES_INTERNES.9) 'VHWD' 1 1 1;
  128. D3 = EXTR (TAB.VARIABLES_INTERNES.14) 'VHWD' 1 1 1;
  129. D4 = EXTR (TAB.VARIABLES_INTERNES.18) 'VHWD' 1 1 1;
  130. *
  131. P1 = EXTR (TAB.VARIABLES_INTERNES.7) 'EPSE' 1 1 1;
  132. P2 = EXTR (TAB.VARIABLES_INTERNES.9) 'EPSE' 1 1 1;
  133. P3 = EXTR (TAB.VARIABLES_INTERNES.14) 'EPSE' 1 1 1;
  134. P4 = EXTR (TAB.VARIABLES_INTERNES.18) 'EPSE' 1 1 1;
  135. *
  136. S1 = EXTR (TAB.CONTRAINTES.7) 'SMZZ' 1 1 1;
  137. S2 = EXTR (TAB.CONTRAINTES.9) 'SMZZ' 1 1 1;
  138. S3 = EXTR (TAB.CONTRAINTES.14) 'SMZZ' 1 1 1;
  139. S4 = EXTR (TAB.CONTRAINTES.18) 'SMZZ' 1 1 1;
  140. *
  141. EF_D = PROG D1 D2 D3 D4 ;
  142. EF_P = PROG P1 P2 P3 P4 ;
  143. EF_S = PROG S1 S2 S3 S4 ;
  144. *
  145. PROG_ERD = MAXI (ABS ((EF_D - REF_D) / REF_D)) ;
  146. PROG_ERP = MAXI (ABS ((EF_P - REF_P) / REF_P)) ;
  147. PROG_ERS = MAXI (ABS ((EF_S - REF_S) / REF_S)) ;
  148. *
  149. ERR = MAXI (PROG PROG_ERD PROG_ERP PROG_ERS) ;
  150. *
  151. temps;
  152. SI ( ERR <EG 0.05 );
  153. ERRE 0;
  154. SINON;
  155. ERRE 5;
  156. FINSI;
  157. *
  158. FIN;
  159.  
  160.  
  161.  
  162.  
  163.  
  164.  
  165.  
  166.  
  167.  
  168.  
  169.  
  170.  
  171.  

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