Télécharger la plaquette Cast3M


fr en

Télécharger vitetfor.procedur

Retour à la liste

Numérotation des lignes : 

  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.  

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