Télécharger la plaquette Cast3M


fr en

Télécharger dycham.eso

Retour à la liste

Numérotation des lignes : 

dycham

  1. C DYCHAM    SOURCE    BP208322  20/09/18    21:16:18     10718          
  2.       SUBROUTINE DYCHAM(XDEP,XVIT,XRAID,XJEU,XAMO,ETA,
  3.      &                  XFL,DFDX,DFDV,IPERM,iannul)
  4.       IMPLICIT INTEGER(I-N)
  5.       IMPLICIT REAL*8(A-H,O-Z)
  6. *--------------------------------------------------------------------*
  7. *     operateurs DYn* : calcul de la force d'un CHoc AMorti          *
  8. *--------------------------------------------------------------------*
  9. *                                                                    *
  10. *     Parametres:                                                    *
  11. *                                                                    *
  12. * e   XDEP    Valeur du deplacement.                                 *
  13. * e   XVIT    Valeur de la vitesse.                                  *
  14. * e   XRAID   Valeur de la raideur.                                  *
  15. * e   XJEU    Valeur du jeu.                                         *
  16. * e   XAMO    Valeur de l'amortissement.                             *
  17. * e   ETA     Parametre de regularisation.                           *
  18. * e   IPERM   =2 si liaison permanente                               *
  19. *  s  XFL     Valeur de la force de choc.                            *
  20. *  s  DFDX    Derivee de la force de choc.                           *
  21. *                                                                    *
  22. *--------------------------------------------------------------------*
  23.  
  24. *     Contact permanent (<=> ressort+amortisseur localise)
  25.       IF (IPERM.EQ.2) THEN
  26. cbp         IF (iannul.EQ.0) THEN
  27.             XVAL = XDEP - XJEU
  28.             XFL  = -XRAID*XVAL - XAMO*XVIT
  29.             DFDX = -XRAID
  30.             DFDV = -XAMO
  31. cbp         ENDIF
  32.  
  33. *     Contact strict      
  34.       ELSEIF (ETA.EQ.0.D0) THEN
  35.          IF ((XDEP.GE.XJEU).AND.(iannul.EQ.0)) THEN
  36.             XVAL = XDEP - XJEU
  37.             XFL  = -XRAID*XVAL - XAMO*XVIT
  38.             DFDX = -XRAID
  39.             DFDV = -XAMO
  40.             IF ( XFL.GT.0.D0) THEN 
  41.                XFL = 0.D0
  42.                DFDX = 0.D0 
  43.                DFDV = 0.D0
  44.             ENDIF
  45.          ELSE
  46.             XFL = 0.D0
  47.             DFDX = 0.D0
  48.             DFDV = 0.D0
  49.          ENDIF
  50.  
  51. *     Contact regularise
  52.       ELSE
  53.          XVAL = XDEP - XJEU
  54.          XFLR = XRAID*XVAL+XAMO*XVIT
  55.          XAUX = SQRT(XFLR**2 + 4.*XRAID*ETA)
  56.          XFL  = -0.5*(XFLR  + XAUX)      
  57.          DFDX = -0.5*XRAID*(1.+XFLR/XAUX)
  58.          DFDV = -0.5*XAMO *(1.+XFLR/XAUX)
  59.       ENDIF
  60. *
  61. *
  62.       END
  63.  
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  

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