Procedures | Cast3M

   1 : * M_DAMPIN  PROCEDUR  AM        96/12/23    21:47:16     2448           
   2 : *********************************************************************
   3 : 'DEBPROC' m_dampin tmod*'TABLE' stiff*'RIGIDITE' damp*'LISTREEL';
   4 : *********************************************************************
   5 : * matrice de damping modal dans l'espace physique
   6 : * formulation en rigidite
   7 : *********************************************************************
   8 : *
   9 : nmod=('DIME' tmod) - 2;
  10 : mmod= 'DIME' damp;
  11 : 'SI' ('OU' ('NEG' mmod nmod) ('NEG' tmod.'SOUSTYPE' 'BASE_DE_MODES'));
  12 :   'MESS' 'tache impossible';
  13 : 'SINON';
  14 :   j=0; 'REPE' loopm mmod; j=j+1;
  15 :     frj=tmod.j .'FREQUENCE';
  16 :     mdj=tmod.j .'DEFORMEE_MODALE';
  17 :     maj=tmod.j .'MASSE_GENERALISEE';
  18 :     omj=2*Pi*frj;
  19 :     kdj=2 * (extr damp j) / 100 / (omj**3) / maj;
  20 :     daj='XXT' (stiff*mdj) kdj;
  21 :     si (ega j 1); dat=daj; sinon; dat=dat et daj; finsi;
  22 :   'FIN' loopm;
  23 : 'FINSI';
  24 : *
  25 : 'FINPROC' dat;
  26 :