Télécharger fluaendo.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : fluaendo.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. complet = faux ;
  5. *
  6. *************************************************
  7. * pour calcul complet mettre complet à : vrai;
  8. *************************************************
  9. OPTI ECHO 1 DIME 2 ELEM QUA8 MODE AXIS;
  10. *
  11. *-----------------------------------------------------------------
  12. * MAILLAGE AXISYMETRIQUE EPROUVETTE CYLINDRIQUE
  13. * MATERIAU VISCO-PLASTIQUE ENDOMMAGEABLE DEPENDANT DE LA
  14. * =============
  15. * TEMPERATURE POUR N M KK A R EVOL
  16. *-----------------------------------------------------------------
  17. P1 = 0 0; P2 = 3E-3 0; P3 = 3E-3 30E-3; P4 = 0 30E-3;
  18. *
  19. L1 = P1 P2 DROIT 1 ;
  20. L2 = P2 P3 DROIT 1 ;
  21. L3 = P3 P4 DROIT 1 ;
  22. L4 = P4 P1 DROIT 1 ;
  23. *
  24. * mesh
  25. *
  26. EPROU = L1 L2 L3 L4 DALLER PLAN ;
  27. *
  28. * boundary conditions
  29. *
  30. CL1 = BLOQ L1 UZ ;
  31. CL3 = BLOQ L4 UR ;
  32. CL = CL1 ET CL3 ;
  33. *
  34. * MODE defines the behavior of the material and the finite element
  35. * formulation
  36. *
  37. MO = MODE EPROU MECANIQUE ELASTIQUE VISCOPLASTIQUE VISCODOMMAGE;
  38. *
  39. * where is the load applied ?
  40. *
  41. pres=pression mass mo -137. l3 ;
  42. *
  43. * Definition des coefficients variables avec la temperature
  44. *
  45. PROGKSI1 = PROG 0. 2.E5 ;
  46. PROGK1 = PROG 15. 15. ;
  47. CTRAC1 = EVOL MANU KSI PROGKSI1 K PROGK1 ;
  48. *
  49. PROGKSI2 = PROG 0. 2.E5 ;
  50. PROGK2 = PROG 75. 75. ;
  51. CTRAC2 = EVOL MANU KSI PROGKSI2 K PROGK2 ;
  52. *
  53. CTRAC = NUAGE 'T'*'FLOTTANT' 'EVOL'*'EVOLUTION'
  54. 1000. CTRAC1 1050. CTRAC2 ;
  55. *
  56. * lors de l'ecoulement :
  57. * interpolation des coefficients avec la temperature
  58. * car les listes de temperatures ne sont pas identiques pour N
  59. * et les autre coef.
  60. EVN = EVOL MANU 'T' (PROG 1000 1025 1050 ) 'N '
  61. (PROG 9 7 5);
  62. EVM = EVOL MANU 'T' (PROG 1000 1050) 'M ' (PROG 8 4 );
  63. EVKK = EVOL MANU 'T' (PROG 1000 1050) 'KK ' (PROG 2200 800 );
  64. EVR = EVOL MANU 'T' (PROG 1000 1050) 'R ' (PROG 7 3 );
  65. EVA = EVOL MANU 'T' (PROG 1000 1050) 'A ' (PROG 639 1261 );
  66. *
  67. * Definition du materiau variable
  68. *
  69. MATVAR = MATE MO YOUN 150000. NU 0.3 rho 7800. alph 0.
  70. N EVN M EVM KK EVKK ALP1 0. BLP1 0.
  71. R EVR A EVA EVOL CTRAC SMAX 0. ;
  72. *
  73. * Definition des cartes de temperature
  74. *
  75. TEMP0 = MANU CHPO EPROU 1 'T' 1000. ;
  76. TEMP1 = MANU CHPO EPROU 1 'T' 1050. ;
  77. *
  78. * Definition du chargement
  79. *
  80. LI1 = PROG 0. 3.E7 ; LI2 = PROG 1. 1. ;
  81. EV = EVOL MANU T LI1 LOAD LI2 ;
  82. CHA1 = CHAR 'MECA' PRES EV ;
  83. TEMPSs = TABLE; TEMPE = TABLE;
  84. TEMPSs.0 = 0. ;TEMPSs.1 = 6000.;
  85. TEMPE.0 = TEMP0; TEMPE.1 = TEMP1;
  86. CHA2 = CHAR 'T' TEMPSs TEMPE;
  87. CHA = CHA1 ET CHA2;
  88.  
  89. TAB = TABLE ;
  90. TAB.'BLOCAGES_MECANIQUES' = CL;
  91. TAB.'CARACTERISTIQUES' = MATVAR;
  92. TAB.'MODELE' = MO;
  93. TAB.'CHARGEMENT' = CHA;
  94. *
  95. si complet;
  96. LIS = PROG 0. pas 1.E-1 1. 1.5 2 pas 1 10 15 20 pas 10 100
  97. 150 200 pas 100 800. 850. pas 50 1100 pas 20 1400 ;
  98. sinon;
  99. LIS = PROG 0. pas 1.E-1 1. 1.5 2 ;
  100. finsi;
  101. *
  102. TAB.'TEMPS_CALCULES' = LIS;
  103. PASAPAS TAB ;
  104. *
  105. *----------------------------------------------------------------
  106. * CONTROLE DES RESULTATS AVEC DE LA SOLUTION DE REFERENCE
  107. * OBTENUE PAR ALGORITHME ou par code
  108. *----------------------------------------------------------------
  109. *
  110. si complet;
  111. REF_D = 5.87E-2 ;
  112. REF_P = 1.36E-3 ;
  113. sinon;
  114. REF_D = 3.98156E-05;
  115. REF_P = 1.51511E-05;
  116. finsi;
  117. *
  118. *tind=index (tab.resucont); ntind=dime tind ;
  119. *T=tind.ntind;
  120. *
  121. D = EXTR ( PECHE TAB VARIABLES_INTERNES) 'VHWD' 1 1 1;
  122. P = EXTR ( PECHE TAB VARIABLES_INTERNES) 'EPSE' 1 1 1;
  123. *
  124. list D; list P;
  125. errd = ABS ((REF_D - D) / REF_D) ;
  126. errp = ABS ((REF_P - P) / REF_P) ;
  127. *
  128. err = MAXI (prog errd errp) ;
  129. *
  130. temps;
  131. SI ( ERR <EG 0.05 );
  132. ERRE 0;
  133. SINON;
  134. ERRE 5;
  135. FINSI;
  136. *
  137. FIN;
  138.  
  139.  
  140.  
  141.  
  142.  
  143.  
  144.  
  145.  
  146.  
  147.  
  148.  
  149.  
  150.  

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