Procedures | Cast3M

   1 : * CHAMINT   PROCEDUR  GBM       04/12/22    21:15:02     4973           
   2 : *//////////////////////////////////////////////////////////////////
   3 : * integre le chargement sur le temps, sort un chargement
   4 : 'DEBPROC'  CHAMINT CH1*'CHARGEMENT' ;
   5 :   ncha = 'DIME' CH1 ;
   6 :   xstart = -1.D15;
   7 : * boucle sur les sous chargements
   8 :   'REPETER' bouc ncha ;
   9 : *    on extrait la nieme evolution
  10 :      lstp = 'EXTRAIRE' CH1 'EVOL' &bouc ;
  11 :      px = 'EXTRAIRE'  lstp ABSC ;
  12 :      py = 'EXTRAIRE'  lstp ORDO ;
  13 :      xdeb = 'EXTRAIRE'  px 1 ;
  14 :      xfin = 'EXTRAIRE'  px ('DIME'  px) ; 
  15 :      ydeb = 'EXTRAIRE'  py 1 ;
  16 :      yfin = 'EXTRAIRE'  py ('DIME'  py) ;
  17 :      n = 'DIME'  px ;
  18 :      xx0 =  xdeb ;
  19 :      yy0 = ydeb ;
  20 :      pry = 0.D0 ;
  21 :      p_prim = 'PROG'  pry ;
  22 :      chsp = 'EXTRAIRE' CH1 'CHAM' &bouc ;
  23 :      SI (maxi (abs chsp) > 1.D-30);
  24 :        i = 2 ;
  25 :        'REPETER' bb (n - 1) ;
  26 :         xx1 = 'EXTRAIRE' px i ;
  27 :         dx = xx1 - xx0 ;
  28 :         yy1 = 'EXTRAIRE' py i ;
  29 : *       integration trapeze
  30 :         pry = (yy1 + yy0) * dx * 0.5D0 + pry ;
  31 :         p_prim = p_prim 'ET' ('PROG' pry);
  32 :         'SI' ((('ABS' (yy1)) > 1.D-30) 'ET' (xx1 > xstart));
  33 : *           'MESSAGE' 'CEST INFERIEUR' yy1;
  34 :            xstart = xx1;
  35 :         'FINSI';
  36 :         xx0 = xx1 ;
  37 :         yy0 = yy1 ;
  38 :         i = i + 1 ;
  39 :         'FIN'  bb ;
  40 : *       On insere au nouveau chargement
  41 :         ev2 = 'EVOL'  manu px p_prim ;
  42 :         'SI' (&bouc 'EGA' 1) ;
  43 :            CH2 = 'CHARGEMENT' chsp ev2 ;
  44 :         'SINON' ;
  45 :            CH2 = CH2 'ET' ('CHARGEMENT' chsp ev2) ;
  46 :         'FINSI' ;
  47 :      SINON ;
  48 :         'SI' (&bouc 'EGA' 1) ;
  49 :            CH2 = 'CHARGEMENT' chsp lstp ;
  50 :         'SINON' ;
  51 :            CH2 = CH2 'ET' ('CHARGEMENT' chsp lstp) ;
  52 :         'FINSI' ;
  53 :      FINSI;
  54 :   'FIN' bouc;
  55 : *  mess 'SALUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUT' xstart;
  56 : 'FINPROC' CH2 xstart;
  57 : *\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\
  58 :