Procedures | Cast3M

   1 : * PAS_EPTH  PROCEDUR  CHAT      13/03/21    21:15:07     7735           
   2 : 'DEBP' PAS_EPTH   PRECED*TABLE ;
   3 : *
   4 : * cette procedure calcule de la deformation thermique de reference 
   5 : * et  calcule de la deformation thermique initiale
   6 : * doit se faire une fois en début de calcul (initiale ou reprise )
   7 : * car maillage a pu changer. dans ce cas on croise les doigts en esperant 
   8 : * que les redu fonctionnent!
   9 : * résutalts dans  WTAB.'ETREF'  WTAB.'ETHER2'  WTAB.'MSRTHER2'
  10 :  WTAB=PRECED.'WTABLE'; 
  11 : * mess ' pas_epth mat0';list wtab.'MAT0';
  12 :  MATMECA= 'REDU'   WTAB.'MAT0'  WTAB.'MO_TOT'  ;            
  13 :  WTAB.'ETREF' = 'EPTH' WTAB.'MO_TOT' MATMECA
  14 :            (WTAB.'TEMPERATURE_REFERENCE' - WTAB.'TALPHA_REFERENCE') ;           
  15 :  TST = WTAB.'TET1' - WTAB.'TALPHA_REFERENCE'; 
  16 :  'OUBL' MATMECA;
  17 :  MATMECA= 'REDU'   WTAB.'MAT1'  WTAB.'MO_TOT'  ; 
  18 :  ETT = 'EPTH' WTAB.'MO_TOT' MATMECA TST ;  
  19 :  'OUBL' MATMECA;                  
  20 : * on met le resultat dans ether2 car on fera dans la boucle ether1=ethetr2  
  21 :  WTAB.'ETHER2' = ETT - WTAB.'ETREF';  
  22 : *           
  23 : *                                                                       
  24 :       'SI'  WTAB.'POR1' ;                                                      
  25 :           MAT_POR1 = REDU WTAB.'MAT1' WTAB.'MOD_POR' ;
  26 :           ZALM   = ('EXCO' 'ALPM' MAT_POR1 'MSR0')
  27 :                  'CHAN' 'TYPE' 'SCALAIRE' ;                             
  28 :           ZBI    = ('EXCO' 'COB' MAT_POR1 'MSR0' )
  29 :                  'CHAN' 'TYPE' 'SCALAIRE';                              
  30 :           ZALPH  =  ( 'EXCO' 'ALPH' MAT_POR1  'MSR0')
  31 :                          'CHAN' 'TYPE' 'SCALAIRE'    ;         
  32 : 
  33 :           ZFAC = 3. * ( ZALM - ( ZBI * ZALPH ) );                       
  34 :           ZFAC = ZFAC  'CHAN' 'TYPE' 'SCALAIRE';                        
  35 : *                                                                       
  36 :           ZDT = 'REDU'  ( WTAB.'TET1' -
  37 :                    WTAB.'TEMPERATURE_REFERENCE') WTAB.'MAI_POR';             
  38 :           ZDDT = 'CHAN' 'CHAM' ZDT WTAB.'MOD_POR' 'STRESSES' 'SCALAIRE';        
  39 :           ZFAC = 'CHAN' 'STRESSES' WTAB.'MOD_POR' ZFAC ;                        
  40 :                                                                         
  41 :           MSRTT  = ZFAC  * ZDDT ;                                       
  42 :           WTAB.'MSRTHER2' = MSRTT ;                                   
  43 :           'DETR' ZALM ; 'DETR' ZBI ; 'DETR' ZALPH ;                     
  44 :           'DETR' ZFAC ; 'DETR' ZDT ; 'DETR' ZDDT ;                      
  45 :       'FINSI';                                                          
  46 : *                                                                               
  47 : 'FINPROC';
  48 :  
  49 :  
  50 :  
  51 :  
  52 :  
  53 :