1 : C DYCHAM2 SOURCE CHAT 05/01/12 23:03:08 5004
2 : C DYCHAM2 SOURCE LAVARENN 96/08/05 21:17:56 2237
3 : SUBROUTINE DYCHAM2(XDEP,XVIT,XJEU,XAMO,XABSCI,XORDON,NIP,
4 : & NLIAB,I,XFL,IPERM,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 elementaire. *
13 : * Cas d'une loi elastique *
14 : * Param}tres: *
15 : * *
16 : * e XDEP Valeur du d{placement. *
17 : * e XVIT Valeur de la vitesse. *
18 : * e XJEU Valeur du jeu. *
19 : * e XAMO Valeur de l'amortissement. *
20 : * e IPERM =1 si liaison permanente *
21 : * e XABSCI Abscisses loi elastique *
22 : * e XORDON Ordonnees loi elastique *
23 : * *
24 : * s XFL Valeur de la force de choc. *
25 : * *
26 : * *
27 : * Auteur, date de cr{ation: *
28 : * *
29 : * Bruno FONTAINE, le 26/02/98 *
30 : * *
31 : *--------------------------------------------------------------------*
32 : REAL*8 XABSCI(NLIAB,*),XORDON(NLIAB,*)
33 : XFL = 0.D0
34 : *
35 : if (iannul.eq.0) then
36 :
37 : IF (IPERM.EQ.1) THEN
38 : XVAL = abs(XDEP)
38 : CALL LIRANG(XVAL,XABSCI,XORDON,I,NLIAB,NRG,XPENTE,NIP)
39 : XFLA = (XORDON(I,NRG) + XPENTE*(XVAL-XABSCI(I,NRG)))
40 : if (xdep.lt.0d0) xfla = -xfla
41 : xfl = -xfla - XAMO * XVIT
42 : ELSE
43 : IF (XDEP.GE.XJEU) THEN
44 : XVAL = XDEP - XJEU
44 : CALL LIRANG(XVAL,XABSCI,XORDON,I,NLIAB,NRG,XPENTE,NIP)
45 : XFL = - (XORDON(I,NRG) + XPENTE*(XVAL-XABSCI(I,NRG)))
46 : & - XAMO * XVIT
47 : IF ( XFL.GT.0.D0) XFL = 0.D0
48 : ENDIF
49 : ENDIF
50 :
51 : endif
52 : *
53 : *
54 : END
55 :
56 :
57 :
58 :
59 :
© Cast3M 2003 - All rights reserved.
Disclaimer