Télécharger dychpl2.eso

Retour à la liste

Numérotation des lignes :

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

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