Télécharger la plaquette Cast3M


fr en

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.   'TALP' 0. 'TREF' 1000.
  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.  
  151.  
  152.  

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