C DGCHAD SOURCE CHAT 05/01/12 22:48:29 5004 SUBROUTINE DGCHAD(IDIM,IGP,XPALB,NLIAB,I,I0,I1,I2,I3, + XADHE,XRAIT,XAMOT,XFN,XFT,XPUS,XVITG,XFORT,XPAS) IMPLICIT INTEGER(I-N) IMPLICIT REAL*8(A-H,O-Z) *---------------------------------------------------------------- * * * Op{rateur DYNE : algorithme de Fu - de Vogelaere * * ________________________________________________ * * * * LE PAS DE TEMPS DOIT ETRE CHOISI DE TELLE FACON QUE LE * * RESSORT TANGENTIEL NE PUISSE PAS PASSER DE LA COMPRESSION * * ET A L EXTENSION DURANT UN PAS DE TEMPS CE QUI EXCLUT LE * * PASSAGE BRUTAL D UN GLISSEMENT DANS UN SENS A UN GLISSEMENT * * DANS L AUTRE SENS * * * * * * * * * Param}tres: * * * * e IGP Indice de glissement au pas pr{c{dent * * es XPALB Tableau contenant les param}tres de la liaison * * e NLIAB Nombre de liaisons sur la base B * e I Num{ro de la liaison * * e I0 Pointeur de la vitesse tangentielle dans XPALB * * e I1 Pointeur de la vitesse tangentielle au pas pr{c{dent * * dans XPALB * * e I2 Pointeur de la force tangentielle au pas * * pr{c{dent dans XPALB * * e I3 Pointeur de la vitesse du glisseur au pas * pr{c{dent dans XPALB * * e XFN Valeur de la force de choc suivant la normale. * * e XGLIS Valeur du coefficient de glissement * * e XADHE Valeur du coefficient d'adh{rence * * e XRAIT Valeur de la raideur tangentielle * * e XAMOT Valeur de l'amortissement tangentiel * * s XFT Module de la force de choc projet{e sur le plan tangent* * s XPUS Valeur de la puissance d'usure instantan{e * * e XPAS Pas de temps * * * * * Auteurs Heidi NEDJAI Rene GIBERT * *--------------------------------------------------------------------* *CALCUL DES FORCES FT(n)=XFORT ET VG(n)=XVITG = 0 * *ON AVAIT ADHERENCE AU PAS (n-1) AVEC IGP=0 * *HYPOTHESE: ON CONTINUE A ADHERER AU PAS (n) * * * DIMENSION XPALB(NLIAB,*) DIMENSION XFORT (3) DIMENSION XVITG (3) DO 10 ID = 1,IDIM XVITG(ID) = 0.D0 10 CONTINUE *END DO * RECHERCHE DE XRAIT XPAS XAMOT PSXFT = 0.D0 DO 20 ID = 1,IDIM XFORT(ID) = XPALB(I,I2+ID) + - (XRAIT * XPAS + XAMOT) * XPALB(I,I0+ID) + - (XAMOT * (XPALB(I,I1+ID) - XPALB(I,I3+ID))) * PSXFT = PSXFT + XFORT(ID) * XFORT(ID) 20 CONTINUE XFT = SQRT(PSXFT) RETURN END