Procedures | Cast3M

   1 : * VITETFOR  PROCEDUR  AM        96/09/09    21:15:09     2269           
   2 : 'DEBPROC' VITETFOR   MTA*'TABLE' ETAB12*'TABLE'
   3 :                      FCORU*'CHPOINT' ;
   4 : *
   5 : * MTA   table de la rigidite associee  au contact
   6 : *
   7 :   FLAG = VRAI ;
   8 :   MAN = ETAB12.'FMAN' ;
   9 :   MAN = 'CHAN' 'ATTRIBUT' MAN 'NATURE' 'DISCRET' ;
  10 : *
  11 :   MDI = 'DIME' MTA;
  12 :  'SI' ( 'NEG' MDI 0) ;
  13 :     A = MTA.7 ;
  14 :     LISEA = MTA.6;
  15 : 
  16 : 
  17 :     MALI = 'MAXI' LISEA;
  18 :    'SI' ( 'EGA' MALI 1 ) ;
  19 :       GEO1 = 'EXTR' A 'MAIL' 'MULT' ;
  20 : *
  21 : *  on recupere quelques quantites necessaires
  22 : *
  23 :       H = ETAB12.'DT' ;
  24 :       H2 = H * H ;
  25 : 
  26 : *------
  27 : *
  28 : *  calcul de l'operateur
  29 : *
  30 :       ATOT = ETAB12.'MASSE'  'ET' A ;
  31 : *
  32 : *  calcul du 1-er second membre et de COE1
  33 : *
  34 :       FU =  ( -0.25 * H2 ) * MAN;
  35 : 
  36 :       VDT =  ( H * ETAB12.'VITI')  +  ETAB12.'DEPI' ;
  37 :       VDT = 'CHAN' 'ATTRIBUT' VDT 'NATURE' 'DIFFUS' ;
  38 :       FL = 'EXCO' ( A * VDT) 'FLX' 'FLX'  ;
  39 :       FL = FL - ( 'REDU' FCORU GEO1 ) ;
  40 :       FT = FL 'ET' FU ;
  41 :       COE1 = 'RESO'  ATOT FT 'NOUNIL' ;
  42 : *
  43 : *  calcul du 2-eme second membre et de COE2
  44 : *
  45 :       FL =  'EXCO' ( A * ETAB12.'VITI' ) 'FLX' 'FLX' ;
  46 :       FU =  (  -0.5 * H ) * MAN;
  47 :       FT = FL 'ET' FU ;
  48 :       COE2 = 'RESO'  ATOT FT 'NOUNIL' ;
  49 : *
  50 : *   calcul des vecteurs VADD et FADD
  51 : *
  52 :       DTL =  (( 2. / H ) * COE1 ) - COE2 ;
  53 :       DTLL = 'REAC' A  DTL ;
  54 :       VADD =  'RESO' ETAB12.'MASSE'  DTLL  ;
  55 :       FADD = (-4./ H2 ) *  ( 'REAC' A COE1) ;
  56 :       'SI' ( 'EXIS' ETAB12 'AMORTISSEMENT' ) ;
  57 :            FADD = FADD - (ETAB12.'AMORTISSEMENT' * VADD );
  58 :       'FINSI' ;
  59 :        FLAG = FAUX ;
  60 :    'FINSI';
  61 :  'FINSI';
  62 : *
  63 : 'SI' FLAG ;
  64 :   VADD = 0. * ETAB12.'VITI' ;
  65 :   FADD = 0. * MAN ;
  66 : 'FINSI' ;
  67 : VADD =  'ENLE' VADD 'LX' ;
  68 : *
  69 : 
  70 : 'FINPROC' VADD FADD ;
  71 : 
  72 : 
  73 :  
  74 :