1 : C DGCHFR SOURCE CHAT 05/01/12 22:48:33 5004
2 : SUBROUTINE DGCHFR(IDIM,IGP,XPALB,NLIAB,I,I0,I1,I2,I3,I4,
3 : + XGLIS,XADHE,XRAIT,XAMOT,XFN,XFT,XPUS,XPAS)
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 : * LE PAS DE TEMPS DOIT ETRE CHOISI DE TELLE FACON QUE LE *
12 : * RESSORT TANGENTIEL NE PUISSE PAS PASSER DE LA COMPRESSION *
13 : * ET A L EXTENSION DURANT UN PAS DE TEMPS CE QUI EXCLUT LE *
14 : * PASSAGE BRUTAL D UN GLISSEMENT DANS UN SENS A UN GLISSEMENT *
15 : * DANS L AUTRE SENS *
16 : * *
17 : * *
18 : * Gestion des phases glissement et adh{rence du frottement de *
19 : * COULOMB *
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 : DIMENSION XPALB(NLIAB,*)
46 : DIMENSION XFORT (3)
47 : DIMENSION XVITG (3)
48 : XFT = 0.D0
49 : XPUS = 0.D0
50 :
51 : * modele d' Antunez
52 : IF (IGP.EQ.1) THEN
53 : *
54 : *CAS1 IL Y AVAIT GLISSEMENT AU PAS PRECEDENT ON SUPPOSE QUE LE *
55 : *GLISSEMENT CONTINUE DANS LE MEME SENS *
55 : CALL DGCHGL(IDIM,IGP,XPALB,NLIAB,I,I0,I1,I2,I3,XGLIS,
56 : + XRAIT,XAMOT,XFN,XFT,XPUS,XVITG,XFORT,XPAS)
57 : *ON FAIT LE TEST SUR LES VITESSES DU GLISSEUR *
58 : *PRODUIT SCALAIRE DE LA VITESSE DU GLISSEUR AU PAS n ET AU *
59 : *PAS n-1 *
60 : *
61 : PSVITG = 0.D0
62 : DO 10 ID = 1,IDIM
63 : PSVITG = PSVITG + XVITG(ID) * XPALB(I,I3+ID)
64 : *
65 : 10 CONTINUE
66 : * end do
67 : IF (PSVITG.GT.0.D0) THEN
68 : *IL Y A GLISSEMENT
69 : IGP = 1
70 : *
71 : DO 11 ID = 1,IDIM
72 : XPALB(I,I2+ID) = XFORT(ID)
73 : XPALB(I,I3+ID) = XVITG(ID)
74 :
75 : 11 CONTINUE
76 : * end do
77 : ELSE
77 : CALL DGCHAD(IDIM,IGP,XPALB,NLIAB,I,I0,I1,I2,I3,XADHE,
78 : + XRAIT,XAMOT,XFN,XFT,XPUS,XVITG,XFORT,XPAS)
79 : IGP = 0
80 : DO 30 ID = 1,IDIM
81 : XPALB(I,I2+ID) = XFORT(ID)
82 : XPALB(I,I3+ID) = XVITG(ID)
83 : 30 CONTINUE
84 :
85 : ENDIF
86 : ELSE
87 : *(ELSE DE IGP=1 DU DEBUT)
88 : *CAS2 IL Y AVAIT ADHERENCE AU PAS PRECEDENT ,ON SUPPOSE QUE *
89 : *L ADHERENCE CONTINUE
90 :
90 : CALL DGCHAD(IDIM,IGP,XPALB,NLIAB,I,I0,I1,I2,I3,XADHE,
91 : + XRAIT,XAMOT,XFN,XFT,XPUS,XVITG,XFORT,XPAS)
92 : *
93 :
94 : *ON FAIT LE TEST
95 : IF (XFT.GT.(XADHE * ABS(XFN))) THEN
96 : *L ADHERENCE NE TIENT PAS
97 : IGP = 1
97 : CALL DGCHGL(IDIM,IGP,XPALB,NLIAB,I,I0,I1,I2,I3,XGLIS,
98 : + XRAIT,XAMOT,XFN,XFT,XPUS,XVITG,XFORT,XPAS)
99 :
100 : DO 40 ID = 1,IDIM
101 : XPALB(I,I2+ID) = XFORT(ID)
102 : XPALB(I,I3+ID) = XVITG(ID)
103 : ***ICI ON RAJOUTE BOUCLE 40
104 :
105 : 40 CONTINUE
106 : *End do
107 : ELSE
108 : IGP = 0
109 :
110 : DO 50 ID = 1,IDIM
111 : XPALB(I,I2+ID) = XFORT(ID)
112 : XPALB(I,I3+ID) = XVITG(ID)
113 : ***ICI ON RAJOUTE BOUCLE 50
114 :
115 : 50 CONTINUE
116 :
117 : endif
118 : endif
119 : end
120 :
121 :
122 :
123 :
124 :
125 :
126 :
127 :
128 :
© Cast3M 2003 - All rights reserved.
Disclaimer