1 : C DYCHA3 SOURCE CHAT 05/01/12 23:00:43 5004
2 : SUBROUTINE DYCHA3(XDEP,XVITN,IDIM,IGP,XPALB,NLIAB,I,INTER
3 : & ,XFN,XFT,XPUS,iannul)
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 : * Calcul de la force de choc pour un choc {l{mentaire. *
12 : * *
13 : * Param}tres: *
14 : * *
15 : * e XDEP Valeur du d{placement sur la normale *
16 : * e XVITN Valeur de la vitesse normale au plan tangent du contact*
17 : * e IGP Indice de glissement au pas pr{c{dent *
18 : * es XPALB Tableau contenant les param}tres de la liaison *
19 : * e NLIAB Nombre de liaisons sur la base B *
20 : * e I Num{ro de la liaison *
21 : * e INTER Indice de liaison interieure / exterieure *
22 : * s XFN Valeur de la force de choc suivant la normale. *
23 : * s XFT Module de la force de choc projetee sur le plan tangent*
24 : * s XPUS Valeur de la puissance d'usure instantanee *
25 : * *
26 : * *
27 : * Auteur, date de cr{ation: *
28 : * *
29 : * Bertrand BEAUFILS, le 31 mai 1990. *
30 : * *
31 : *--------------------------------------------------------------------*
32 : DIMENSION XPALB(NLIAB,*)
33 : XFN = 0.D0
34 : XFT = 0.D0
35 : XPUS = 0.D0
36 : XJEU = XPALB(I,2)
37 : *
38 : XAD = ABS(XDEP)
39 : XAJ = ABS(XJEU)
40 : *
41 : IF ( (((XAD.LE.XAJ).AND.(INTER.EQ.0)).OR.
42 : & ((XAD.GE.XAJ).AND.(INTER.EQ.1))).AND.
43 : & (iannul.EQ.0)) THEN
44 : ID4 = 7 + 3*IDIM
45 : ID5 = 7 + 4*IDIM
46 : ID6 = 7 + 5*IDIM
47 : ID8 = 7 + 7*IDIM
48 : ID9 = 7 + 8*IDIM
49 : XRAIN = XPALB(I,1)
50 : XAMON = XPALB(I,7)
51 : XGLIS = XPALB(I,3)
52 : XADHE = XPALB(I,4)
53 : XRAIT = XPALB(I,5)
54 : XAMOT = XPALB(I,6)
55 : XVAL = XDEP - XJEU
56 : XFN = - XRAIN * XVAL - XAMON * XVITN
56 : CALL DYCHFR(IDIM,IGP,XPALB,NLIAB,I,ID4,ID5,ID6,ID8,ID9,XGLIS,
57 : & XADHE,XRAIT,XAMOT,XFN,XFT,XPUS)
58 : ELSE
59 : IGP = -1
60 : ID6 = 7 + 5*IDIM
61 : DO 10 ID=1,IDIM
62 : XPALB(I,ID6+ID) = 0.D0
63 : 10 CONTINUE
64 : * end do
65 : ENDIF
66 : *
67 : RETURN
68 : END
69 :
70 :
71 :
© Cast3M 2003 - All rights reserved.
Disclaimer