1 : C DYCHAD SOURCE CHAT 05/01/12 23:02:39 5004
2 : SUBROUTINE DYCHAD(IDIM,IGP,XPALB,NLIAB,I,I0,I1,I2,I3,XADHE,XRAIT,
3 : & XAMOT,XFN,XFT,XPUS)
4 : IMPLICIT INTEGER(I-N)
5 : IMPLICIT REAL*8(A-H,O-Z)
6 : *--------------------------------------------------------------------*
7 : * *
8 : * Op{rateur DYNE : algorithme de Fu - de Vogelaere *
9 : * ________________________________________________ *
10 : * *
11 : * Calcul de la force d'adh{rence et test de sa validit{ *
12 : * *
13 : * Param}tres: *
14 : * *
15 : * es XPALB Tableau contenant les param}tres de la liaison *
16 : * e NLIAB Nombre de liaisons sur la base B *
17 : * e I Num{ro de la liaison *
18 : * e I0 Pointeur de la vitesse tangentielle dans XPALB *
19 : * e I1 Pointeur de la variation de d{placement par rapport a *
20 : * la position origine d'adh{rence, projet{e sur le plan *
21 : * tangent, dans XPALB *
22 : * e I2 Pointeur de la force tangentielle dans XPALB *
23 : * e I3 Pointeur de la force tangentielle @ la fin du glis- *
24 : * -sement dans XPALB *
25 : * e XFN Valeur de la force de choc suivant la normale. *
26 : * e XADHE Valeur du coefficient d'adh{rence *
27 : * e XRAIT Valeur de la raideur tangentielle *
28 : * e XAMOT Valeur de l'amortissement tangentiel *
29 : * es IGP Indice de glissement au pas pr{c{dent *
30 : * s XFT Module de la force de choc projet{e sur le plan tangent*
31 : * s XPUS Valeur de la puissance d'usure instantan{e *
32 : * *
33 : * *
34 : * Auteur, date de cr{ation: *
35 : * *
36 : * Bertrand BEAUFILS, le 31 mai 1990. *
37 : * *
38 : *--------------------------------------------------------------------*
39 : DIMENSION XPALB(NLIAB,*)
40 : * Calcul de la force tangentielle introduite par l'oscillateur
41 : * {quivalent
42 : PSXFT = 0.D0
43 : DO 10 ID = 1,IDIM
44 : XPALB(I,I2+ID) = XPALB(I,I3+ID) - XRAIT * XPALB(I,I1+ID)
45 : & - XAMOT * XPALB(I,I0+ID)
46 : PSXFT = PSXFT + XPALB(I,I2+ID) * XPALB(I,I2+ID)
47 : 10 CONTINUE
48 : * end do
49 : XFT = SQRT(PSXFT)
50 : *
51 : IF (XFT.GT.( XADHE * ABS(XFN) ) ) THEN
52 : * L'adh{rence ne tient pas
53 : IGP = 1
54 : ELSE
55 : * Ca colle
56 : XPUS = 0.D0
57 : IGP = 0
58 : ENDIF
59 : *
60 : RETURN
61 : END
62 :
© Cast3M 2003 - All rights reserved.
Disclaimer