Télécharger la plaquette Cast3M


fr en

Télécharger dycha5.eso

Retour à la liste

Numérotation des lignes : 

dycha5

  1. C DYCHA5    SOURCE    CHAT      05/01/12    23:02:11     5004
  2. C DYCHA4    SOURCE    LAVARENN  96/08/05    21:17:45     2237
  3.       SUBROUTINE DYCHA5(XNET,XTE,XXXN,XVITN,IDIM,IGP
  4.      &,XPALB,NLIAB,I,IPALB,XFN,XFT,XPUS,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 elementaire            *
  13. *     avec amortissement pour la liaison ligne_cercle                *
  14. *     Param}tres:                                                    *
  15. *                                                                    *
  16. * e   XDEP    Valeur du d{placement sur la normale                   *
  17. * e   XVITN   Valeur de la vitesse normale au plan tangent du contact*
  18. * e   IGP     Indice de glissement au pas pr{c{dent                  *
  19. * es  XPALB   Tableau contenant les param}tres de la liaison         *
  20. * e   NLIAB   Nombre de liaisons sur la base B                       *
  21. * e   I       Num{ro de la liaison                                   *
  22. *  s  XFN     Valeur de la force de choc suivant la normale.         *
  23. *  s  XFT     Module de la force de choc projetee sur le plan tangent*
  24. *  s  XPUS    Valeur de la puissance d'usure instantanee             *
  25. *                                                                    *
  26. *                                                                    *
  27. *     Auteur, date de cr{ation:                                      *
  28. *                                                                    *
  29. *    IBRAHIM PINTO, 05/97                                            *
  30. *                                                                    *
  31. *--------------------------------------------------------------------*
  32.  
  33.  
  34.  
  35.       REAL*8 XNET,XTE(3),XXXN(3)
  36.  
  37.  
  38.  
  39.       DIMENSION XPALB(NLIAB,*)
  40.       INTEGER IPALB(NLIAB,*)
  41.       XFN = 0.D0
  42.       XFT = 0.D0
  43.       XPUS = 0.D0
  44.       XRAY  = XPALB(I,2)
  45.       XRAIT = XPALB(I,5)
  46.  
  47.  
  48.       ITYP =IPALB(I,1)
  49.  
  50. *
  51.       IF ((XNET.GE.XRAY).AND.(iannul.EQ.0)) THEN
  52.          ID2 = 7 + IDIM
  53.          ID3 = 7 + 2*IDIM
  54.          ID4 = 7 + 3*IDIM
  55.          ID6 = 7 + 5*IDIM
  56.          ID7 = 7 + 6*IDIM
  57.          XRAIN = XPALB(I,1)
  58.          XAMON = XPALB(I,7)
  59.          XGLIS = XPALB(I,3)
  60.          XADHE = XPALB(I,4)
  61.          XRAIT = XPALB(I,5)
  62.          XAMOT = XPALB(I,6)
  63.  
  64.  
  65. *----dans le cas ou on suppose la normale de contact
  66. *----dans le plan du cercle on ne change pas xval
  67.  
  68.          XVAL = XNET - XRAY
  69.  
  70.  
  71. *----sinon il faut projeter sur la normale reelle de contact
  72.          IF (ITYP.EQ.40) THEN
  73.  
  74.  
  75.             XVAL=0.D0
  76.             DO 4 ID=1,IDIM
  77.               XVAL=XVAL+XTE(ID)*XXXN(ID)
  78.  4         CONTINUE
  79.             XVAL=XVAL*(XNET-XRAY)/XNET
  80.          ENDIF
  81.  
  82.  
  83. *----CALCUL DE LA FORCE DE CHOC NORMALE
  84.  
  85.          XFN = - XRAIN * XVAL - XAMON * XVITN
  86.  
  87.  
  88.        CALL DYCHFR(IDIM,IGP,XPALB,NLIAB,I,ID2,ID3,ID4,ID6
  89.      &,ID7,XGLIS,XADHE,XRAIT,XAMOT,XFN,XFT,XPUS)
  90.  
  91.          IF ( XFN.GT.0.D0) XFN = 0.D0
  92.        ELSE
  93.          IGP = -1
  94.          ID4 = 7 + 3*IDIM
  95.          DO 10 ID=1,IDIM
  96.             XPALB(I,ID4+ID) = 0.D0
  97.  10         CONTINUE
  98. *        end do
  99.          ENDIF
  100. *
  101.  
  102.       RETURN
  103.  
  104.  
  105.       END
  106.  
  107.  
  108.  
  109.  
  110.  

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