Télécharger la plaquette Cast3M


fr en

Télécharger dychpl.eso

Retour à la liste

Numérotation des lignes : 

dychpl

  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