1 : C DYCHE4 SOURCE CHAT 05/01/12 23:04:34 5004
2 : SUBROUTINE DYCHE4(XDEP,IDIM,IGP,XPALB,NLIAB,I,XFN,XFT,XPUS,iannul)
3 : IMPLICIT INTEGER(I-N)
4 : IMPLICIT REAL*8(A-H,O-Z)
5 : *--------------------------------------------------------------------*
6 : * *
7 : * Op{rateur DYNE : algorithme de Fu - de Vogelaere *
8 : * ________________________________________________ *
9 : * *
10 : * Appel du module de frottement pour la liaison *
11 : * POINT_PLAN_FROTTEMENT sans amortissement *
12 : * *
13 : * Param}tres: *
14 : * *
15 : * e XDEP Valeur du d{placement sur la normale *
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 liaison sur la base B *
19 : * e I Num{ro de la liaison *
20 : * s XFN Valeur de la force de choc suivant la normale. *
21 : * s XFT Module de la force de choc projetee sur le plan tangent*
22 : * s XPUS Valeur de la puissance d'usure instantanee *
23 : * *
24 : * *
25 : * Auteur, date de cr{ation: *
26 : * *
27 : * Bertrand BEAUFILS, le 5 juin 1990. *
28 : * *
29 : *--------------------------------------------------------------------*
30 : DIMENSION XPALB(NLIAB,*)
31 : XFN = 0.D0
32 : XFT = 0.D0
33 : XPUS = 0.D0
34 : XJEU = XPALB(I,2)
35 : *
36 : IF ((XDEP.GE.XJEU).AND.(iannul.EQ.0)) THEN
37 : ID2 = 6 + IDIM
38 : ID3 = 6 + 2*IDIM
39 : ID4 = 6 + 3*IDIM
40 : ID6 = 6 + 5*IDIM
41 : ID7 = 6 + 6*IDIM
42 : XRAIN = XPALB(I,1)
43 : XGLIS = XPALB(I,3)
44 : XADHE = XPALB(I,4)
45 : XRAIT = XPALB(I,5)
46 : XAMOT = XPALB(I,6)
47 : XVAL = XDEP - XJEU
48 : XFN = - XRAIN * XVAL
48 : CALL DYCHFR(IDIM,IGP,XPALB,NLIAB,I,ID2,ID3,ID4,ID6,ID7,XGLIS,
49 : & XADHE,XRAIT,XAMOT,XFN,XFT,XPUS)
50 : ELSE
51 : ID4 = 6 + 3*IDIM
52 : IGP = -1
53 : DO 10 ID=1,IDIM
54 : XPALB(I,ID4+ID) = 0.D0
55 : 10 CONTINUE
56 : * end do
57 : ENDIF
58 : *
59 : RETURN
60 : END
61 :
62 :
© Cast3M 2003 - All rights reserved.
Disclaimer