Télécharger dychpl2.eso

Retour à la liste

Numérotation des lignes :

  1. C DYCHPL2 SOURCE BP208322 19/02/25 21:15:49 10120
  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. * Operateur DYNE : algorithme de Fu - de Vogelaere *
  10. * ________________________________________________ *
  11. * *
  12. * Calcul de la force de choc pour un choc elementaire point-plan
  13. * elastoplastique, la force restant toujours positive
  14. * *
  15. * Parametres: *
  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 =2 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.2) THEN
  40. cbp 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. cbp ENDIF
  57. ELSE
  58. IF ((XVAL .GE. 0.D0).AND.(iannul.EQ.0)) THEN
  59. CALL LIRANG(XSEUIL,XORDON,XABSCI,I,NLIAB,NRG,XPENTE,NIP)
  60. XE = XABSCI(I,NRG) + XPENTE*(XSEUIL-XORDON(I,NRG))
  61. CALL LIRANG(XVAL,XABSCI,XORDON,I,NLIAB,NRG,XPENTE,NIP)
  62. XFLR = - (XORDON(I,NRG) + XPENTE*(XVAL-XABSCI(I,NRG)))
  63. DXFLP = XSEUIL - (ABS (XFLR))
  64. IF (DXFLP.LE.0.D0) THEN
  65. XDPLAS = XDEP - XJEU - XE
  66. XFLR = - XSEUIL
  67. ENDIF
  68. XFL = XFLR - XAMO*XVIT
  69. IF (XFL.GT.0.D0) XFL = 0.D0
  70. ELSE
  71. XFL = 0.D0
  72. ENDIF
  73. ENDIF
  74.  
  75. END
  76.  
  77.  
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  
  84.  

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