1 : C DYCHA5 SOURCE CHAT 05/01/12 23:02:11 5004
2 : C DYCHA4 SOURCE LAVARENN 96/08/05 21:17:45 2237
3 : SUBROUTINE DYCHA5(XNET,XTE,XXXN,XVITN,IDIM,IGP
4 : &,XPALB,NLIAB,I,IPALB,XFN,XFT,XPUS,iannul)
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 : * Calcul de la force de choc pour un choc elementaire *
13 : * avec amortissement pour la liaison ligne_cercle *
14 : * Param}tres: *
15 : * *
16 : * e XDEP Valeur du d{placement sur la normale *
17 : * e XVITN Valeur de la vitesse normale au plan tangent du contact*
18 : * e IGP Indice de glissement au pas pr{c{dent *
19 : * es XPALB Tableau contenant les param}tres de la liaison *
20 : * e NLIAB Nombre de liaisons sur la base B *
21 : * e I Num{ro de la liaison *
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 : * IBRAHIM PINTO, 05/97 *
30 : * *
31 : *--------------------------------------------------------------------*
32 :
33 :
34 :
35 : REAL*8 XNET,XTE(3),XXXN(3)
36 :
37 :
38 :
39 : DIMENSION XPALB(NLIAB,*)
40 : INTEGER IPALB(NLIAB,*)
41 : XFN = 0.D0
42 : XFT = 0.D0
43 : XPUS = 0.D0
44 : XRAY = XPALB(I,2)
45 : XRAIT = XPALB(I,5)
46 :
47 :
48 : ITYP =IPALB(I,1)
49 :
50 : *
51 : IF ((XNET.GE.XRAY).AND.(iannul.EQ.0)) THEN
52 : ID2 = 7 + IDIM
53 : ID3 = 7 + 2*IDIM
54 : ID4 = 7 + 3*IDIM
55 : ID6 = 7 + 5*IDIM
56 : ID7 = 7 + 6*IDIM
57 : XRAIN = XPALB(I,1)
58 : XAMON = XPALB(I,7)
59 : XGLIS = XPALB(I,3)
60 : XADHE = XPALB(I,4)
61 : XRAIT = XPALB(I,5)
62 : XAMOT = XPALB(I,6)
63 :
64 :
65 : *----dans le cas ou on suppose la normale de contact
66 : *----dans le plan du cercle on ne change pas xval
67 :
68 : XVAL = XNET - XRAY
69 :
70 :
71 : *----sinon il faut projeter sur la normale reelle de contact
72 : IF (ITYP.EQ.40) THEN
73 :
74 :
75 : XVAL=0.D0
76 : DO 4 ID=1,IDIM
77 : XVAL=XVAL+XTE(ID)*XXXN(ID)
78 : 4 CONTINUE
79 : XVAL=XVAL*(XNET-XRAY)/XNET
80 : ENDIF
81 :
82 :
83 : *----CALCUL DE LA FORCE DE CHOC NORMALE
84 :
85 : XFN = - XRAIN * XVAL - XAMON * XVITN
86 :
87 :
87 : CALL DYCHFR(IDIM,IGP,XPALB,NLIAB,I,ID2,ID3,ID4,ID6
88 : &,ID7,XGLIS,XADHE,XRAIT,XAMOT,XFN,XFT,XPUS)
89 :
90 : IF ( XFN.GT.0.D0) XFN = 0.D0
91 : ELSE
92 : IGP = -1
93 : ID4 = 7 + 3*IDIM
94 : DO 10 ID=1,IDIM
95 : XPALB(I,ID4+ID) = 0.D0
96 : 10 CONTINUE
97 : * end do
98 : ENDIF
99 : *
100 :
101 : RETURN
102 :
103 :
104 : END
105 :
106 :
107 :
108 :
© Cast3M 2003 - All rights reserved.
Disclaimer