1 : C DGCHAD SOURCE CHAT 05/01/12 22:48:29 5004
2 :
3 : SUBROUTINE DGCHAD(IDIM,IGP,XPALB,NLIAB,I,I0,I1,I2,I3,
4 : + XADHE,XRAIT,XAMOT,XFN,XFT,XPUS,XVITG,XFORT,XPAS)
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 : * LE PAS DE TEMPS DOIT ETRE CHOISI DE TELLE FACON QUE LE *
13 : * RESSORT TANGENTIEL NE PUISSE PAS PASSER DE LA COMPRESSION *
14 : * ET A L EXTENSION DURANT UN PAS DE TEMPS CE QUI EXCLUT LE *
15 : * PASSAGE BRUTAL D UN GLISSEMENT DANS UN SENS A UN GLISSEMENT *
16 : * DANS L AUTRE SENS
17 : * *
18 : * *
19 : * *
20 : * *
21 : * Param}tres: *
22 : * *
23 : * e IGP Indice de glissement au pas pr{c{dent *
24 : * es XPALB Tableau contenant les param}tres de la liaison *
25 : * e NLIAB Nombre de liaisons sur la base B
26 : * e I Num{ro de la liaison *
27 : * e I0 Pointeur de la vitesse tangentielle dans XPALB *
28 : * e I1 Pointeur de la vitesse tangentielle au pas pr{c{dent *
29 : * dans XPALB *
30 : * e I2 Pointeur de la force tangentielle au pas *
31 : * pr{c{dent dans XPALB *
32 : * e I3 Pointeur de la vitesse du glisseur au pas
33 : * pr{c{dent dans XPALB *
34 : * e XFN Valeur de la force de choc suivant la normale. *
35 : * e XGLIS Valeur du coefficient de glissement *
36 : * e XADHE Valeur du coefficient d'adh{rence *
37 : * e XRAIT Valeur de la raideur tangentielle *
38 : * e XAMOT Valeur de l'amortissement tangentiel *
39 : * s XFT Module de la force de choc projet{e sur le plan tangent*
40 : * s XPUS Valeur de la puissance d'usure instantan{e *
41 : * e XPAS Pas de temps * *
42 : * *
43 : * Auteurs Heidi NEDJAI Rene GIBERT *
44 : *--------------------------------------------------------------------*
45 : *CALCUL DES FORCES FT(n)=XFORT ET VG(n)=XVITG = 0 *
46 : *ON AVAIT ADHERENCE AU PAS (n-1) AVEC IGP=0 *
47 : *HYPOTHESE: ON CONTINUE A ADHERER AU PAS (n) *
48 : *
49 : *
50 : DIMENSION XPALB(NLIAB,*)
51 : DIMENSION XFORT (3)
52 : DIMENSION XVITG (3)
53 : DO 10 ID = 1,IDIM
54 : XVITG(ID) = 0.D0
55 : 10 CONTINUE
56 : *END DO
57 : * RECHERCHE DE XRAIT XPAS XAMOT
58 :
59 : PSXFT = 0.D0
60 : DO 20 ID = 1,IDIM
61 : XFORT(ID) = XPALB(I,I2+ID)
62 : + - (XRAIT * XPAS + XAMOT) * XPALB(I,I0+ID)
63 : + - (XAMOT * (XPALB(I,I1+ID) - XPALB(I,I3+ID)))
64 : *
65 : PSXFT = PSXFT + XFORT(ID) * XFORT(ID)
66 :
67 : 20 CONTINUE
68 : XFT = SQRT(PSXFT)
69 :
70 : RETURN
71 : END
72 :
© Cast3M 2003 - All rights reserved.
Disclaimer