1 : C DYNE16 SOURCE CHAT 05/01/12 23:14:35 5004
2 : SUBROUTINE DYNE16(Q1,Q2,Q3,NA1,FTOTA,XPTB,NPLB,IDIMB,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 : * D{calage ou mise @ z{ro dans les tableaux pour le pas suivant. *
11 : * *
12 : * Param}tres: *
13 : * *
14 : * es Q1( , ) Vecteur des d{placements g{n{ralis{s *
15 : * es Q2( , ) Vecteur des vitesses g{n{ralis{es *
16 : * es Q3( , ) Vecteur des acc{l{rations g{n{ralis{es *
17 : * e NA1 Nombre total d'inconnues sur la base A *
18 : * es FTOTA Forces ext{rieures totalis{es, sur la base A *
19 : * es XPTB Tableau des d{placements des points de liaisons *
20 : * e NPLB Nombre maxi de points de liaisons *
21 : * e IDIMB Nombre maxi de ddl par liaisons *
22 : * *
23 : * Auteur, date de cr{ation: *
24 : * *
25 : * Denis ROBERT-MOUGIN, le 16 juin 1989. *
26 : * *
27 : *--------------------------------------------------------------------*
28 : REAL*8 Q1(NA1,*),Q2(NA1,*),Q3(NA1,*),FTOTA(NA1,*)
29 : REAL*8 XPTB(NPLB,4,*),FINERT(NA1,*)
30 : *
31 : PARAMETER ( ZERO=0.D0 )
32 : *
33 : * Boucle sur les inconnues de la base A
34 : *
35 : DO 10 I=1,NA1
36 : Q1(I,3) = Q1(I,1)
37 : Q1(I,4) = Q1(I,2)
38 : Q1(I,1) = ZERO
39 : Q1(I,2) = ZERO
40 : *
41 : Q2(I,3) = Q2(I,1)
42 : Q2(I,4) = Q2(I,2)
43 : Q2(I,1) = ZERO
44 : Q2(I,2) = ZERO
45 : *
46 : Q3(I,3) = Q3(I,1)
47 : Q3(I,4) = Q3(I,2)
48 : Q3(I,1) = ZERO
49 : Q3(I,2) = ZERO
50 : *
51 : FTOTA(I,3) = FTOTA(I,1)
52 : FTOTA(I,4) = FTOTA(I,2)
53 : FTOTA(I,1) = ZERO
54 : FTOTA(I,2) = ZERO
55 : *
56 : FINERT(I,3) = FINERT(I,1)
57 : FINERT(I,4) = FINERT(I,2)
58 : FINERT(I,1) = ZERO
59 : FINERT(I,2) = ZERO
60 : 10 CONTINUE
61 : * end do
62 : DO 20 IP = 1,NPLB
63 : DO 22 ID = 1,IDIMB
64 : XPTB(IP,3,ID) = XPTB(IP,1,ID)
65 : XPTB(IP,4,ID) = XPTB(IP,2,ID)
66 : XPTB(IP,1,ID) = ZERO
67 : XPTB(IP,2,ID) = ZERO
68 : 22 CONTINUE
69 : * end do
70 : 20 CONTINUE
71 : * end do
72 : *
73 : END
74 :
© Cast3M 2003 - All rights reserved.
Disclaimer