Télécharger endoaxi3.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : endoaxi3.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Endoaxi3.dgibi: Jeux de données *
  5. * ------------------------------------ *
  6. * *
  7. complet = faux;
  8. *
  9. *************************************************
  10. * pour calcul complet mettre complet à : vrai;
  11. ***************************************************
  12. *******************************************************
  13. * *
  14. * TEST; MATERIAU ELASTO-PLASTIQUE ENDOMMAGEABLE *
  15. * MATERIAU DEPENDANT DE LA TEMPERATURE : *
  16. * COURBE DE TRACTION *
  17. * COEFFICIENT DE DILATATION *
  18. * THERMIQUE *
  19. * TEMPERATURE DE REFERENCE = 30°C *
  20. * TEMPERATURE DE REFERENCE POUR COEF. DE DILATATION *
  21. * THERMIQUES = 10°C *
  22. * TEMPERATURE INITIALE = 50°C *
  23. * *
  24. * EPROUVETTE EN TRACTION AVEC DEPLACEMENTS IMPOSES *
  25. * + CHARGEMENT THERMIQUE *
  26. * MAILLAGE AXISYMETRIQUE. *
  27. * *
  28. * LES RESULTATS OBTENUS SONT COMPARES A LA SOLUTION *
  29. * ANALYTIQUE OBTENUE PAR ALGORITHME. *
  30. * ON TESTE - LA DEFORMATION PLASTIQUE CUMULEE *
  31. * (composante EPSE dans VARIABLES_INTERNES)*
  32. * - D, LA VARIABLE D'ENDOMMAGEMENT *
  33. * DANS VARIABLES_INTERNES . *
  34. * - SIGMAzz ( DANS CONTRAINTES ). *
  35. ******************************************************* *
  36.  
  37.  
  38. OPTI ECHO 0 DIME 2 ELEM QUA8 MODE AXIS;
  39. *------------------------------------------------------
  40. P1 = 0 0; P2 = 3E-3 0;
  41. P3 = 3E-3 30E-3; P4 = 0 30E-3;
  42. *
  43. L1 = P1 P2 DROIT 1 ;
  44. L2 = P2 P3 DROIT 1 ;
  45. L3 = P3 P4 DROIT 1 ;
  46. L4 = P4 P1 DROIT 1 ;
  47. *
  48. EPROU = L1 L2 L3 L4 DALLER PLAN ;
  49.  
  50. CL1 = BLOQ L1 UZ ;
  51. CL2 = BLOQ L3 UZ ;
  52. CL3 = BLOQ L4 UR ;
  53. CL = CL1 ET CL2 ET CL3 ;
  54. DEP = DEPI CL2 -1E-3 ;
  55. *
  56. MO = MODE EPROU 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE'
  57. 'ENDOMMAGEABLE';
  58. *
  59. * module tangent valant 783MPa
  60. *-----------------------------
  61. PROGEPS1 = PROG 0. (350./74020.) 0.3 0.6 ;
  62. PROGSIG1 = PROG 0. 350.E06 581.1976E06 816.0976E06 ;
  63. CTRAC1 = EVOL MANU EPS PROGEPS1 SIG PROGSIG1 ;
  64. *
  65. * module tangent valant 491MPa
  66. *-----------------------------
  67. PROGEPS2 = PROG 0. (350./74020.) 0.3 0.6 ;
  68. PROGSIG2 = PROG 0. 350.E06 494.9783E06 642.2783E06 ;
  69. CTRAC2 = EVOL MANU EPS PROGEPS2 SIG PROGSIG2 ;
  70. *
  71. * module tangent valant 199MPa
  72. *-----------------------------
  73. PROGEPS3 = PROG 0. (350./74020.) 0.3 0.6 ;
  74. PROGSIG3 = PROG 0. 350.E06 408.7590E06 468.4590E06 ;
  75. CTRAC3 = EVOL MANU EPS PROGEPS3 SIG PROGSIG3 ;
  76. *
  77. * module tangent variable
  78. *------------------------
  79. CTRAC = NUAGE 'T'*'FLOTTANT' 'TRAC'*'EVOLUTION'
  80. 0. CTRAC1 100. CTRAC2 200. CTRAC3 ;
  81. *
  82. * coefficient de dilatation thermique vaiable
  83. *--------------------------------------------
  84. *
  85. PROGTEMP = PROG 0. 100. 200. ;
  86. PROGALFA = PROG 0. .25E-3 .5E-3 ;
  87. EVALPH = EVOL MANU 'T' PROGTEMP 'ALPH' PROGALFA ;
  88. *
  89. MATvar = MATE MO YOUN 74020.E06 NU 0.3 RHO 7800.
  90. ALPH EVALPH TRAC CTRAC EPSD 0. DC 0.23 EPSR 0.25 ;
  91. *
  92. ************** DEFINITION DU CHARGEMENT ***************
  93. *
  94. LI1 = PROG 0. 1 31. ; LI2 = PROG 0. 0. 30. ;
  95. EV = EVOL MANU T LI1 DEP(T) LI2 ;
  96. CHA1 = CHAR 'DIMP' DEP EV ;
  97. TEMPREF = MANU CHPO EPROU 1 'T' 30. ;
  98. TEMP0 = MANU CHPO EPROU 1 'T' 50. ;
  99. TEMP200 = MANU CHPO EPROU 1 'T' 200. ;
  100. t0 = 0. ; T1 = 1. ; t20 = 20. ;
  101. TEMP = TABLE; TEMPERA = TABLE;
  102. TEMP.0 = t0 ; TEMP.1 = t1 ; TEMP.2 = t20 ;
  103. TEMPERA.0 = TEMPREF;TEMPERA.1 = TEMP0;
  104. TEMPERA.2 = TEMP200;
  105. CHA2 = CHAR 'T' TEMP TEMPERA;
  106. CHA = CHA1 ET CHA2;
  107. *
  108. * LISTE DES PAS DE CALCUL :
  109. * 1.035 MM PREMIER INCREMENT PLASTIQUE
  110. *
  111. si complet ;
  112. LIS = PROG 1.035 1.1 pas 0.1 7;
  113. sinon;
  114. LIS = PROG 1.035 1.1 pas 0.1 1.6;
  115. finsi;
  116.  
  117. *
  118. TAB = TABLE ;
  119. TAB.'PRECISION' = 1E-3;
  120. TAB.'BLOCAGES_MECANIQUES' = CL;
  121. TAB.'CARACTERISTIQUES' = MATvar;
  122. TAB.'CHARGEMENT' = CHA;
  123. TAB.'TEMPS_CALCULES' = LIS;
  124. TAB.'MODELE' = MO;
  125. TAB.'TALPHA_REFERENCE' = 10. ;
  126. TAB.'TEMPERATURE_REFERENCE' = TEMPREF ;
  127. PASAPAS TAB ;
  128. *
  129. *------------------------------------------------------
  130. * CONTROLE DES RESULTATS AVEC DE LA SOLUTION
  131. * DE REFERENCE OBTENUE PAR ALGORITHME
  132. *------------------------------------------------------
  133. *
  134. si complet;
  135. TA = 4.; TB = 7.;
  136. REF_D = PROG 7.5E-2 17.7E-2 ;
  137. REF_P = PROG 10.3E-2 20.9E-2 ;
  138. REF_S = (PROG -376.9 -369.1)*1E6 ;
  139. SINON;
  140. TA = 1.2; TB = 1.6;
  141. REF_D = PROG 3.52404E-03 1.24384E-02 ;
  142. REF_P = PROG 5.75205E-03 1.97372E-02 ;
  143. REF_S = PROG -3.52473E+08 -3.57838E+08;
  144. finsi;
  145. *
  146. D1=EXTR (PECHE TAB VARIABLES_INTERNES TA) 'VHWD' 1 1 1;
  147. D2=EXTR (PECHE TAB VARIABLES_INTERNES TB) 'VHWD' 1 1 1;
  148. *
  149. P1=EXTR (PECHE TAB VARIABLES_INTERNES TA) 'EPSE' 1 1 1;
  150. P2=EXTR (PECHE TAB VARIABLES_INTERNES TB) 'EPSE' 1 1 1;
  151. *
  152. S1 = EXTR (PECHE TAB CONTRAINTES TA) 'SMZZ' 1 1 1;
  153. S2 = EXTR (PECHE TAB CONTRAINTES TB) 'SMZZ' 1 1 1;
  154. *
  155. EF_D = PROG D1 D2 ;
  156. EF_P = PROG P1 P2 ;
  157. EF_S = PROG S1 S2 ;
  158. *
  159. PROG_ERD = MAXI (ABS ((EF_D - REF_D) / REF_D)) ;
  160. PROG_ERP = MAXI (ABS ((EF_P - REF_P) / REF_P)) ;
  161. PROG_ERS = MAXI (ABS ((EF_S - REF_S) / REF_S)) ;
  162. *
  163. ERR = MAXI (PROG PROG_ERD PROG_ERP PROG_ERS) ;
  164. *
  165. temps;
  166. SI ( ERR <EG 0.05 );
  167. ERRE 0;
  168. SINON;
  169. ERRE 5;
  170. FINSI;
  171. *
  172. FIN;
  173.  
  174.  
  175.  
  176.  
  177.  
  178.  
  179.  
  180.  
  181.  
  182.  
  183.  
  184.  
  185.  
  186.  
  187.  
  188.  

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