1 : C DYCHFR SOURCE CHAT 05/01/12 23:05:52 5004
2 : SUBROUTINE DYCHFR(IDIM,IGP,XPALB,NLIAB,I,I0,I1,I2,I3,I4,
3 : & XGLIS,XADHE,XRAIT,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 : * Gestion des phases glissement et adh{rence du frottement de *
12 : * COULOMB *
13 : * *
14 : * Param}tres: *
15 : * *
16 : * e IGP Indice de glissement au pas pr{c{dent *
17 : * es XPALB Tableau contenant les param}tres de la liaison *
18 : * e NLIAB Nombre de liaisons sur la base B *
19 : * e I Num{ro de la liaison *
20 : * e I0 Pointeur de la vitesse tangentielle dans XPALB *
21 : * e I1 Pointeur de la variation de d{placement par rapport a *
22 : * la position origine d'adh{rence, projet{e sur le plan *
23 : * tangent, dans XPALB *
24 : * e I2 Pointeur de la force tangentielle dans XPALB *
25 : * e I3 Pointeur de la vitesse tangentielle au pas pr{c{dent *
26 : * dans XPALB *
27 : * e I4 Pointeur de la force tangentielle @ la fin du glis- *
28 : * -sement dans XPALB *
29 : * e XFN Valeur de la force de choc suivant la normale. *
30 : * e XGLIS Valeur du coefficient de glissement *
31 : * e XADHE Valeur du coefficient d'adh{rence *
32 : * e XRAIT Valeur de la raideur tangentielle *
33 : * e XAMOT Valeur de l'amortissement tangentiel *
34 : * s XFT Module de la force de choc projet{e sur le plan tangent*
35 : * s XPUS Valeur de la puissance d'usure instantan{e *
36 : * *
37 : * *
38 : * Auteur, date de cr{ation: *
39 : * *
40 : * Bertrand BEAUFILS, le 31 mai 1990. *
41 : * I. POLITOPOULOS (bricolage pour le model d' Oden) ,juin 1993 *
42 : *--------------------------------------------------------------------*
43 : DIMENSION XPALB(NLIAB,*)
44 : XFT = 0.D0
45 : XPUS = 0.D0
46 : * modele d' Antunez
47 : IF (XRAIT.GT.0.D0) THEN
48 : IF (IGP.EQ.1) THEN
49 : *
50 : * Il y avait glissement au pas pr{cedent ; calcul du produit scalaire
51 : * de la vitesse et de la vitesse au pas pr{c{dent ; calcul du module
52 : * de la vitesse
53 : *
54 : PS = 0.D0
55 : DO 10 ID = 1,IDIM
56 : PS = PS + XPALB(I,I0+ID) * XPALB(I,I3+ID)
57 : 10 CONTINUE
58 : * end do
59 : IF (PS.GT.0.D0) THEN
60 : * Il y a glissement
60 : CALL DYCHGL(IDIM,IGP,XPALB,NLIAB,I,I0,I2,XGLIS,XFN,XFT,XPUS)
61 : ELSE
62 : * Il y a possibilit{ d'adh{rence
62 : CALL DYCHAD(IDIM,IGP,XPALB,NLIAB,I,I0,I1,I2,I4,XADHE,XRAIT,
63 : & XAMOT,XFN,XFT,XPUS)
64 : IF (IGP.EQ.1) THEN
65 : * En fait, il y a glissement.
65 : CALL DYCHGL(IDIM,IGP,XPALB,NLIAB,I,I0,I2,XGLIS,XFN,XFT,
66 : & XPUS)
67 : ELSE
68 : * Ca adh}re
69 : CONTINUE
70 : ENDIF
71 : ENDIF
72 : ELSE
73 : *
74 : * Il y avait adh{rence au pas pr{c{dent ; \a peut continuer
75 : *
75 : CALL DYCHAD(IDIM,IGP,XPALB,NLIAB,I,I0,I1,I2,I4,XADHE,XRAIT,
76 : & XAMOT,XFN,XFT,XPUS)
77 : IF (IGP.EQ.1) THEN
78 : * En fait, il y a glissement.
78 : CALL DYCHGL(IDIM,IGP,XPALB,NLIAB,I,I0,I2,XGLIS,XFN,XFT,XPUS)
79 : ELSE
80 : * Ca adh}re
81 : CONTINUE
82 : ENDIF
83 : ENDIF
84 : ELSE
85 : * modele d'Oden
86 : *
87 : IGP = 1
88 : EPSI = XAMOT
89 : * calcul du module de la vitesse
90 : *
91 : PSV = 0.D0
92 : DO 20 ID = 1,IDIM
93 : PSV = PSV + XPALB(I,I0+ID) * XPALB(I,I0+ID)
94 : 20 CONTINUE
95 :
96 : * end do
97 : VITT = SQRT(PSV)
98 : IF (VITT.GT.EPSI) THEN
99 : * Il y a glissement classique
99 : CALL DYCHGL(IDIM,IGP,XPALB,NLIAB,I,I0,I2,XGLIS,XFN,XFT,XPUS)
100 : ELSE
101 : * on est dans le domaine regularise
101 : CALL DYCGLR(IDIM,EPSI,XPALB,NLIAB,I,I0,I2,XGLIS,XFN,XFT,XPUS)
102 : ENDIF
103 : ENDIF
104 :
105 : RETURN
106 : END
107 :
© Cast3M 2003 - All rights reserved.
Disclaimer