1 : C DYCGLR SOURCE CHAT 05/01/12 23:00:14 5004
2 : C DYCGLR SOURCE BEAU 90/07/26 21:16:02
3 : SUBROUTINE DYCGLR(IDIM,EPSI,XPALB,NLIAB,I,I0,I1,XGLIS,XFN,XFT,
4 : & XPUS)
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 glissement *
13 : * *
14 : * Param}tres: *
15 : * e EPSI Vitesse tangentielle de regularisation d' Oden * *
16 : * es XPALB Tableau contenant les param}tres de la liaison *
17 : * e NLIAB Nombre de liaisons sur la base B *
18 : * e I Num{ro de la liaison *
19 : * e I0 Pointeur de la vitesse tangentielle dans XPALB *
20 : * e I1 Pointeur de la force tangentielle dans XPALB *
21 : * e XFN Valeur de la force de choc suivant la normale. *
22 : * e XGLIS Valeur du coefficient de glissement *
23 : * s XFT Module de la force de choc projet{e sur le plan tangent*
24 : * s XPUS Valeur de la puissance d'usure instantan{e *
25 : * *
26 : * *
27 : * Auteur, date de cr{ation: *
28 : * *
29 : * I. Politopoulos ,Juin 1993 *
30 : * *
31 : *--------------------------------------------------------------------*
32 : DIMENSION XPALB(NLIAB,*)
33 : XFTT = - XGLIS * ABS(XFN)
34 : *
35 : PSVITT = 0.D0
36 : DO 10 ID = 1,IDIM
37 : PSVITT = PSVITT + XPALB(I,I0+ID) * XPALB(I,I0+ID)
38 : 10 CONTINUE
39 : * end do
40 : VITT = SQRT(PSVITT)
41 : PSXFT = 0.D0
42 : DO 11 ID = 1,IDIM
43 : XPALB(I,I1+ID) = XFTT *
44 : & (2- (VITT/EPSI)) * XPALB(I,I0+ID)/EPSI
45 : PSXFT = PSXFT + XPALB(I,I1+ID) * XPALB(I,I1+ID)
46 : 11 CONTINUE
47 : * end do
48 : XFT = SQRT(PSXFT)
49 : XPUS = ABS(XFN * VITT) * (VITT/EPSI)
50 :
51 :
52 : RETURN
53 : END
54 :
© Cast3M 2003 - All rights reserved.
Disclaimer