1 : C DGCHGL SOURCE CHAT 05/01/12 22:48:36 5004
2 :
3 : SUBROUTINE DGCHGL(IDIM,IGP,XPALB,NLIAB,I,I0,I1,I2,I3,
4 : + XGLIS,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 *
46 : *ON AVAIT GLISSEMENT AU PAS (n-1) AVEC IGP=1 *
47 : *HYPOTHESE: ON CONTINUE A GLISSER DANS LE MEME SENS *
48 : *
49 : *
50 : DIMENSION XPALB(NLIAB,*)
51 : DIMENSION XFORT (3)
52 : DIMENSION XVITG (3)
53 : XFT = - XGLIS * ABS(XFN)
54 : IGP = 1
55 : *
56 :
57 : *
58 : *PSFT = PRODUIT SCALAIRE DE LA FORCE PAR ELLE MEME *
59 : PSFT = 0.D0
60 : DO 30 ID = 1,IDIM
61 : PSFT = PSFT + XPALB(I,I2+ID) * XPALB(I,I2+ID)
62 : 30 CONTINUE
63 : *END DO
64 : *
65 : PSFT = SQRT(PSFT)
66 : *
67 : *CALCUL DE LA FORCE DE FROTTEMENT XFORT(n)
68 : DO 40 ID = 1,IDIM
69 : XFORT(ID) = XGLIS * ABS(XFN) * XPALB(I,I2+ID) / PSFT
70 :
71 : 40 CONTINUE
72 : *END DO
73 : *
74 : *CALCUL DE LA VITESSE DU GLISSEUR VG(n)
75 : DO 50 ID = 1,IDIM
76 : XVITG(ID) = XPALB(I,I0+ID) + ((XFORT(ID)
77 : + - XPALB(I,I2+ID)) / (XRAIT * XPAS + XAMOT))
78 : + - (XAMOT * (XPALB(I,I1+ID)
79 : + - XPALB(I,I3+ID)) / (XRAIT * XPAS + XAMOT))
80 : *
81 :
82 : 50 CONTINUE
83 : *END DO
84 :
85 :
86 : RETURN
87 : END
88 :
89 :
90 :
© Cast3M 2003 - All rights reserved.
Disclaimer