Télécharger dycham2.eso

Retour à la liste

Numérotation des lignes :

dycham2
  1. C DYCHAM2 SOURCE BP208322 19/02/25 21:15:47 10120
  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. * Operateur DYNE *
  10. * ________________________________________________ *
  11. * *
  12. * Calcul de la force de choc pour un choc elementaire. *
  13. * Cas d'une LOI_DE_COMPORTEMENT elastique *
  14. * Parametres: *
  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 creation: *
  28. * *
  29. * Bruno FONTAINE, le 26/02/98 *
  30. * *
  31. *--------------------------------------------------------------------*
  32. REAL*8 XABSCI(NLIAB,*),XORDON(NLIAB,*)
  33. XFL = 0.D0
  34. *
  35. cbp if (iannul.eq.0) then
  36.  
  37. IF (IPERM.EQ.2) THEN
  38. cbp XVAL = abs(XDEP)
  39. XVAL = abs(XDEP-XJEU)
  40. CALL LIRANG(XVAL,XABSCI,XORDON,I,NLIAB,NRG,XPENTE,NIP)
  41. XFLA = (XORDON(I,NRG) + XPENTE*(XVAL-XABSCI(I,NRG)))
  42. cbp IF (XDEP.LT.0D0) XFLA = -XFLA
  43. IF (XDEP.LT.XJEU) XFLA = -XFLA
  44. XFL = -XFLA - XAMO * XVIT
  45. ELSE
  46. cbp IF (XDEP.GE.XJEU) THEN
  47. IF ((XDEP.GE.XJEU).AND.(iannul.EQ.0)) THEN
  48. XVAL = XDEP - XJEU
  49. CALL LIRANG(XVAL,XABSCI,XORDON,I,NLIAB,NRG,XPENTE,NIP)
  50. XFL = - (XORDON(I,NRG) + XPENTE*(XVAL-XABSCI(I,NRG)))
  51. & - XAMO * XVIT
  52. IF (XFL.GT.0.D0) XFL = 0.D0
  53. ENDIF
  54. ENDIF
  55.  
  56. cbp endif
  57. *
  58. *
  59. END
  60.  
  61.  
  62.  
  63.  
  64.  
  65.  
  66.  

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