Télécharger dycham2.eso

Retour à la liste

Numérotation des lignes :

  1. C DYCHAM2 SOURCE CHAT 05/01/12 23:03:08 5004
  2. C DYCHAM2 SOURCE LAVARENN 96/08/05 21:17:56 2237
  3. SUBROUTINE DYCHAM2(XDEP,XVIT,XJEU,XAMO,XABSCI,XORDON,NIP,
  4. & NLIAB,I,XFL,IPERM,iannul)
  5. IMPLICIT INTEGER(I-N)
  6. IMPLICIT REAL*8(A-H,O-Z)
  7. *--------------------------------------------------------------------*
  8. * *
  9. * Op{rateur DYNE *
  10. * ________________________________________________ *
  11. * *
  12. * Calcul de la force de choc pour un choc elementaire. *
  13. * Cas d'une loi elastique *
  14. * Param}tres: *
  15. * *
  16. * e XDEP Valeur du d{placement. *
  17. * e XVIT Valeur de la vitesse. *
  18. * e XJEU Valeur du jeu. *
  19. * e XAMO Valeur de l'amortissement. *
  20. * e IPERM =1 si liaison permanente *
  21. * e XABSCI Abscisses loi elastique *
  22. * e XORDON Ordonnees loi elastique *
  23. * *
  24. * s XFL Valeur de la force de choc. *
  25. * *
  26. * *
  27. * Auteur, date de cr{ation: *
  28. * *
  29. * Bruno FONTAINE, le 26/02/98 *
  30. * *
  31. *--------------------------------------------------------------------*
  32. REAL*8 XABSCI(NLIAB,*),XORDON(NLIAB,*)
  33. XFL = 0.D0
  34. *
  35. if (iannul.eq.0) then
  36.  
  37. IF (IPERM.EQ.1) THEN
  38. XVAL = abs(XDEP)
  39. CALL LIRANG(XVAL,XABSCI,XORDON,I,NLIAB,NRG,XPENTE,NIP)
  40. XFLA = (XORDON(I,NRG) + XPENTE*(XVAL-XABSCI(I,NRG)))
  41. if (xdep.lt.0d0) xfla = -xfla
  42. xfl = -xfla - XAMO * XVIT
  43. ELSE
  44. IF (XDEP.GE.XJEU) THEN
  45. XVAL = XDEP - XJEU
  46. CALL LIRANG(XVAL,XABSCI,XORDON,I,NLIAB,NRG,XPENTE,NIP)
  47. XFL = - (XORDON(I,NRG) + XPENTE*(XVAL-XABSCI(I,NRG)))
  48. & - XAMO * XVIT
  49. IF ( XFL.GT.0.D0) XFL = 0.D0
  50. ENDIF
  51. ENDIF
  52.  
  53. endif
  54. *
  55. *
  56. END
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  

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