1 : C DYCHA41 SOURCE CHAT 05/01/12 23:01:19 5004
2 : C DYCHA41 SOURCE LAVARENN 96/08/05 21:17:45 2237
3 : SUBROUTINE DYCHA41(XDEP,XVITN,IDIM,IGP,XPALB,NLIAB,I,XFN,XFT
4 : & ,XPUS,XABSCI,XORDON,NIP,iannul)
5 : IMPLICIT INTEGER(I-N)
6 : IMPLICIT REAL*8(A-H,O-Z)
7 : *--------------------------------------------------------------------*
8 : * *
9 : * Op{rateur DYNE : *
10 : * ________________________________________________ *
11 : * *
12 : * Calcul de la force de choc pour un choc {l{mentaire. *
13 : * Cas d'une loi elastique *
14 : * *
15 : * Param}tres: *
16 : * *
17 : * e XDEP Valeur du d{placement sur la normale *
18 : * e XVITN Valeur de la vitesse normale au plan tangent du contact*
19 : * e IGP Indice de glissement au pas pr{c{dent *
20 : * es XPALB Tableau contenant les param}tres de la liaison *
21 : * e NLIAB Nombre de liaisons sur la base B *
22 : * e I Num{ro de la liaison *
23 : * e XABSCI Abscisses loi elastique *
24 : * e XORDON Ordonnees loi elastique *
25 : * s XFN Valeur de la force de choc suivant la normale. *
26 : * s XFT Module de la force de choc projetee sur le plan tangent*
27 : * s XPUS Valeur de la puissance d'usure instantanee *
28 : * *
29 : * *
30 : * Auteur, date de cr{ation: *
31 : * *
32 : * Bruno FONTAINE, le 26/02/98 *
33 : * *
34 : *--------------------------------------------------------------------*
35 : DIMENSION XPALB(NLIAB,*)
36 : REAL*8 XABSCI(NLIAB,*),XORDON(NLIAB,*)
37 : XFN = 0.D0
38 : XFT = 0.D0
39 : XPUS = 0.D0
40 : XJEU = XPALB(I,2)
41 : *
42 : IF ((XDEP.GE.XJEU).AND.(iannul.EQ.0)) THEN
43 : ID2 = 7 + IDIM
44 : ID3 = 7 + 2*IDIM
45 : ID4 = 7 + 3*IDIM
46 : ID6 = 7 + 5*IDIM
47 : ID7 = 7 + 6*IDIM
48 : XAMON = XPALB(I,7)
49 : XGLIS = XPALB(I,3)
50 : XADHE = XPALB(I,4)
51 : XRAIT = XPALB(I,5)
52 : XAMOT = XPALB(I,6)
53 : XVAL = XDEP - XJEU
53 : CALL LIRANG(XDEP,XABSCI,XORDON,I,NLIAB,NRG,XPENTE,NIP)
54 : XFN = - (XORDON(I,NRG) + XPENTE*(XDEP-XABSCI(I,NRG)))
55 : & - XAMON * XVITN
55 : CALL DYCHFR(IDIM,IGP,XPALB,NLIAB,I,ID2,ID3,ID4,ID6,ID7,XGLIS,
56 : & XADHE,XRAIT,XAMOT,XFN,XFT,XPUS)
57 : IF ( XFN.GT.0.D0) XFN = 0.D0
58 : ELSE
59 : IGP = -1
60 : ID4 = 7 + 3*IDIM
61 : DO 10 ID=1,IDIM
62 : XPALB(I,ID4+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