Procedures | Cast3M

   1 : * M_DAMP_K  PROCEDUR  AM        96/12/23    21:47:23     2448           
   2 : *********************************************************************
   3 : 'DEBPROC' m_damp_K tmod*'TABLE' stiff*'RIGIDITE' damp*'LISTREEL';
   4 : *********************************************************************
   5 : * matrice de damping modal + proportionel en K 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 :   frm=tmod.mmod .'FREQUENCE';
  15 :   omm=2*Pi*frm;
  16 :   ksm=('EXTR' damp mmod) / 100;
  17 :   a1=2 * ksm / omm;
  18 :   dat=a1 * stiff;
  19 :   'SI' (mmod > 1);
  20 :     j=mmod; 'REPE' loopm (mmod-1); j=j-1;
  21 :       frj=tmod.j .'FREQUENCE';
  22 :       mdj=tmod.j .'DEFORMEE_MODALE';
  23 :       maj=tmod.j .'MASSE_GENERALISEE';
  24 :       omj=2*Pi*frj;
  25 :       ksj=('EXTR' damp j) / 100;
  26 :       kbj=ksj - (ksm * omj / omm);
  27 :       kdj=2 * kbj / (omj**3) / maj;
  28 :       daj='XXT' (stiff*mdj) kdj;
  29 :       dat=dat 'ET' daj;
  30 :     'FIN' loopm;
  31 :   'FINSI';
  32 : 'FINSI';
  33 : *
  34 : 'FINPROC' dat;
  35 :