Télécharger la plaquette Cast3M


fr en

Télécharger chamint.procedur

Retour à la liste

Numérotation des lignes : 

  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.  

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