Télécharger dychpl.eso

Retour à la liste

Numérotation des lignes :

  1. C DYCHPL SOURCE CHAT 05/01/12 23:06:57 5004
  2. SUBROUTINE DYCHPL(XDEP,XVIT,XDPLAS,XRAID,XJEU,XAMO,XSEUIL,XFL,
  3. & IPERM,iannul)
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8(A-H,O-Z)
  6. *--------------------------------------------------------------------*
  7. * *
  8. * Op{rateur DYNE : algorithme de Fu - de Vogelaere *
  9. * ________________________________________________ *
  10. * *
  11. * Calcul de la force de choc pour un choc {l{mentaire point-plan
  12. * elastoplastique, la force restant toujours positive
  13. * *
  14. * Param}tres: *
  15. * *
  16. * e XDEP Valeur du d{placement. *
  17. * e XVIT Valeur de la vitesse normale
  18. * e XRAID Valeur de la raideur. *
  19. * e XJEU Valeur de la limite de force elastique (positive)
  20. * e XAMO Valeur de l'amortissement
  21. * e XSEUIL Valeur de la force plastique
  22. * e IPERM =1 si liaison permanente *
  23. * es XDPLAS Valeur du deplacement plastique (du pas precedent
  24. * en entree, du pas actuel en sortie)
  25. * s XFL Valeur de la force de choc. *
  26. * *
  27. * *
  28. * Auteur, date de cr{ation: *
  29. * *
  30. * I. POLITOPOULOS 19/09/1994
  31. * *
  32. *--------------------------------------------------------------------*
  33.  
  34. XFL = 0.D0
  35. XVAL = XDEP - XJEU - XDPLAS
  36.  
  37. IF (IPERM.EQ.1) THEN
  38. IF (iannul.EQ.0) THEN
  39. XE = XSEUIL/XRAID
  40. XFLR = XRAID * XVAL
  41. IF (XFLR.LT.0) XCRIT=-XSEUIL
  42. IF (XFLR.GE.0) XCRIT= XSEUIL
  43.  
  44. DXFLP = XFLR - XCRIT
  45. IF (ABS(XFLR).LE.ABS(XCRIT)) THEN
  46. XFL = XFLR
  47. ELSE
  48. XFL = XCRIT
  49. XDPLAS = XDEP - XJEU - XCRIT/XRAID
  50. ENDIF
  51. XFL = - XFL - XAMO*XVIT
  52.  
  53. ENDIF
  54. ELSE
  55. IF ((XVAL .GE. 0.D0).AND.(iannul.EQ.0)) THEN
  56. XE = XSEUIL/XRAID
  57. XFLR = - XRAID * XVAL
  58. DXFLP = XSEUIL - (ABS (XFLR))
  59. IF (DXFLP.LE.0.D0) THEN
  60. XDPLAS = XDEP - XJEU - XE
  61. XFLR = - XSEUIL
  62. ENDIF
  63. XFL = XFLR - XAMO*XVIT
  64. IF (XFL.GT.0.D0) XFL = 0.D0
  65. ELSE
  66. XFL = 0.D0
  67. ENDIF
  68. ENDIF
  69.  
  70. END
  71.  
  72.  
  73.  
  74.  
  75.  
  76.  

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