1 : C DEVEQ2 SOURCE CHAT 05/01/12 22:44:58 5004
2 : SUBROUTINE DEVEQ2(Q2,NA1,XK,XASM,XM,XDT,NPAS,FTOTA,FINERT)
3 : IMPLICIT INTEGER(I-N)
4 : IMPLICIT REAL*8(A-H,O-Z)
5 : *--------------------------------------------------------------------*
6 : * *
7 : * Op{rateur DYNE : algorithme de Fu - de Vogelaere *
8 : * ________________________________________________ *
9 : * *
10 : * Calcul des vitesses g{n{ralis{es pour le premier demi-pas *
11 : * de temps. *
12 : * *
13 : * Param}tres: *
14 : * *
15 : * es Q2(.,.) Vecteur des vitesses g{n{ralis{es *
16 : * e NA1 Nombre total d'inconnues dans la base A *
17 : * e XK Vecteur des raideurs g{n{ralis{es *
18 : * e XASM Vecteur des amortissements g{n{ralis{s *
19 : * e XM Vecteur des masses g{n{ralis{es *
20 : * e XDT Valeurs des pas de temps *
21 : * e NPAS Num{ro du pas de temps *
22 : * e FTOTA Forces ext{rieures totalis{es, sur la base A *
23 : * *
24 : * Auteur, date de cr{ation: *
25 : * *
26 : * Denis ROBERT-MOUGIN, le 19 mai 1989. *
27 : * *
28 : *--------------------------------------------------------------------*
29 : REAL*8 XK(*),XASM(*),XM(*),XDT(*),FTOTA(NA1,*),Q2(NA1,*)
30 : REAL*8 FINERT(NA1,*)
31 : *
32 : PDT = XDT(NPAS)
33 : PDTS4 = PDT / 4.D0
34 : *
35 : * . 4 . .
36 : * Q = ________ ( Q + h/4 ( G + G - A Q ) )
37 : * i,m+1/2 4 + A h i,m i,m i,m+1/2 i i,m
38 : * I
39 : *
40 : * Boucle sur les inconnues de la base A
41 : *
42 : DO 10 IA = 1,NA1
43 : UNSM2 = 1.D0 / ( XM(IA) - FINERT(IA,2) )
44 : UNSM3 = 1.D0 / ( XM(IA) - FINERT(IA,3) )
45 : AUX1 = 4.D0 / ( 4.D0 + ( XASM(IA) * UNSM2 * PDT ) )
46 : AUX2 = ( FTOTA(IA,3) * UNSM3 ) + ( FTOTA(IA,2) * UNSM2 )
47 : & - ( XASM(IA) * UNSM3 * Q2(IA,3) )
48 : AUX3 = Q2(IA,3) + ( PDTS4 * AUX2 )
49 : Q2(IA,2) = AUX1 * AUX3
50 : 10 CONTINUE
51 : * end do
52 : *
53 : END
54 :
© Cast3M 2003 - All rights reserved.
Disclaimer