Télécharger dychpl.eso

Retour à la liste

Numérotation des lignes :

  1. C DYCHPL SOURCE BP208322 19/02/25 21:15:48 10120
  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. * Operateur DYNE : algorithme de Fu - de Vogelaere *
  9. * ________________________________________________ *
  10. * *
  11. * Calcul de la force de choc pour un choc elementaire point-plan
  12. * elastoplastique, la force restant toujours positive
  13. * *
  14. * Parametres: *
  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 =2 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 creation: *
  29. * *
  30. * I. POLITOPOULOS 19/09/1994
  31. * *
  32. *--------------------------------------------------------------------*
  33.  
  34. XFL = 0.D0
  35. XVAL = XDEP - XJEU - XDPLAS
  36.  
  37. IF (IPERM.EQ.2) THEN
  38. cbp 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. cbp ENDIF
  53. ELSE
  54. IF ((XVAL .GE. 0.D0).AND.(iannul.EQ.0)) THEN
  55. XE = XSEUIL/XRAID
  56. XFLR = - XRAID * XVAL
  57. DXFLP = XSEUIL - (ABS (XFLR))
  58. IF (DXFLP.LE.0.D0) THEN
  59. XDPLAS = XDEP - XJEU - XE
  60. XFLR = - XSEUIL
  61. ENDIF
  62. XFL = XFLR - XAMO*XVIT
  63. IF (XFL.GT.0.D0) XFL = 0.D0
  64. ELSE
  65. XFL = 0.D0
  66. ENDIF
  67. ENDIF
  68.  
  69. END
  70.  
  71.  
  72.  
  73.  
  74.  
  75.  
  76.  

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