C DYCHE5    SOURCE    CHAT      05/01/12    23:04:59     5004
C DYCHE4    SOURCE    LAVARENN  96/08/05    21:18:27     2237
      SUBROUTINE DYCHE5(XNET,XTE,XXXN,IDIM,IGP,XPALB,NLIAB,I
     &,IPALB,XFN,XFT,XPUS,iannul)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
*--------------------------------------------------------------------*
*                                                                    *
*     Op{rateur DYNE : algorithme de Fu - de Vogelaere               *
*     ________________________________________________               *
*                                                                    *
*     Appel du module de frottement pour la liaison                  *
*     ligne_cercle          sans amortissement                       *
*                                                                    *
*     Param}tres:                                                    *
*                                                                    *
* e   XDEP    Valeur du d{placement sur la normale                   *
* e   IGP     Indice de glissement au pas pr{c{dent                  *
* es  XPALB   Tableau contenant les param}tres de la liaison         *
* e   NLIAB   Nombre de liaison sur la base B                        *
* e   I       Num{ro de la liaison                                   *
*  s  XFN     Valeur de la force de choc suivant la normale.         *
*  s  XFT     Module de la force de choc projetee sur le plan tangent*
*  s  XPUS    Valeur de la puissance d'usure instantanee             *
*                                                                    *
*                                                                    *
*     Auteur, date de cr{ation:                                      *
*                                                                    *
*     IBRAHIM PINTO, 05/97                                           *
*                                                                    *
*--------------------------------------------------------------------*


      REAL*8 XNET,XTE(3),XXXN(3)

      DIMENSION XPALB(NLIAB,*)
      INTEGER IPALB(NLIAB,*)
      XFN = 0.D0
      XFT = 0.D0
      XPUS = 0.D0
      XVAL=0.D0
      XRAY  = XPALB(I,2)

      ITYP =IPALB(I,1)


*
      IF ((XNET.GE.XRAY).AND.(iannul.EQ.0)) THEN
         ID2 = 6 + IDIM
         ID3 = 6 + 2*IDIM
         ID4 = 6 + 3*IDIM
         ID6 = 6 + 5*IDIM
         ID7 = 6 + 6*IDIM
         XRAIN = XPALB(I,1)
         XGLIS = XPALB(I,3)
         XADHE = XPALB(I,4)
         XRAIT = XPALB(I,5)
         XAMOT = XPALB(I,6)

*----dans le cas ou on suppose la normale de contact
*----dans le plan du cercle on ne change pas xval

         XVAL = XNET - XRAY

*----sinon il faut projeter sur la normale reelle de contact
         IF ( ITYP.EQ.39 ) THEN
          xval = 0.d0
            DO 4 ID=1,IDIM
               XVAL=XVAL+XTE(ID)*XXXN(ID)
  4         CONTINUE
            XVAL=XVAL*(XNET-XRAY)/XNET
         ENDIF
*----CALCUL DE LA FORCE DE CHOC NORMALE


            XFN = - XRAIN * XVAL

      CALL DYCHFR(IDIM,IGP,XPALB,NLIAB,I,ID2,ID3,ID4,ID6,ID7,XGLIS
     &,XADHE,XRAIT,XAMOT,XFN,XFT,XPUS)
      ELSE
         ID4 = 6 + 3*IDIM
         IGP = -1
         DO 10 ID=1,IDIM
            XPALB(I,ID4+ID) = 0.D0
 10         CONTINUE
*        end do
      ENDIF
*
      RETURN


      END