dyche4
C DYCHE4 SOURCE CHAT 05/01/12 23:04:34 5004 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 * * POINT_PLAN_FROTTEMENT 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: * * * * Bertrand BEAUFILS, le 5 juin 1990. * * * *--------------------------------------------------------------------* DIMENSION XPALB(NLIAB,*) XFN = 0.D0 XFT = 0.D0 XPUS = 0.D0 XJEU = XPALB(I,2) * IF ((XDEP.GE.XJEU).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) XVAL = XDEP - XJEU XFN = - XRAIN * XVAL & 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
© Cast3M 2003 - Tous droits réservés.
Mentions légales