Télécharger lirang.eso

Retour à la liste

Numérotation des lignes :

  1. C LIRANG SOURCE CHAT 05/01/13 01:20:21 5004
  2. SUBROUTINE LIRANG(XTOT,XABSCI,XORDON,I,NLIAB,NRG,XPENTE,NIP)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-H,O-Z)
  5. *--------------------------------------------------------------------*
  6. * *
  7. * Operateur DYNE : algorithme de Fu - de Vogelaere *
  8. * ________________________________________________ *
  9. * *
  10. * Repere l'indice NRG du tableau XABSCI representant la liste de *
  11. * reels des abscisses de la loi de comportement elasto-plastique.*
  12. * Pour une valeur XTOT de l'abscisse, cet indice sera tel que : *
  13. * XABSCI(I,NRG)<= XTOT <= XABSCI(I,NRG+1) *
  14. * *
  15. * Parametres: *
  16. * *
  17. * e XTOT Valeur donnee de l'abscisse *
  18. * e XABSCI Tableau contenant les abscisses de la loi plastique *
  19. * pour la liaison point-point-plastique *
  20. * e I numero de la liaison. *
  21. * pour la liaison point-point-plastique *
  22. * s NRG Rang du point repere *
  23. * * *
  24. * Auteur, date de creation: *
  25. * *
  26. * Nicolas WECXSTEEN modifications 04/96 point-point-..-plastique *
  27. * *
  28. *--------------------------------------------------------------------*
  29. *
  30. REAL*8 XABSCI(NLIAB,*),XORDON(NLIAB,*)
  31. *
  32.  
  33.  
  34. NRG = 1
  35. XVAL = XABSCI(I,2)
  36. 10 IF (XTOT.GE.XVAL) THEN
  37. NRG = NRG + 1
  38. if (nrg.eq.NIP) then
  39. XPENTE = (XORDON(I,NRG)-XORDON(I,NRG-1))/
  40. & (XABSCI(I,NRG)-XABSCI(I,NRG-1))
  41. go to 20
  42. endif
  43.  
  44. XVAL = XABSCI(I,NRG+1)
  45. GOTO 10
  46. endif
  47.  
  48. XPENTE = (XORDON(I,NRG+1)-XORDON(I,NRG))/
  49. & (XABSCI(I,NRG+1)-XABSCI(I,NRG))
  50.  
  51. 20 continue
  52.  
  53. END
  54.  
  55.  
  56.  
  57.  

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