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.  
  59.  

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