1 : C DEVPRO SOURCE CHAT 05/01/12 22:47:04 5004
2 : SUBROUTINE DEVPRO(XPHILB,FTOTB,FTOTBA,IBASB,INMSB,IPLSB,IORSB,
3 : & NSB,NPLSB,NA2,IDIMB,NPLB,NA1)
4 : IMPLICIT INTEGER(I-N)
5 : IMPLICIT REAL*8(A-H,O-Z)
6 : *--------------------------------------------------------------------*
7 : * *
8 : * Op{rateur DYNE : algorithme de Fu - de Vogelaere *
9 : * ________________________________________________ *
10 : * *
11 : * Projection des forces base r{elle sur base A. *
12 : * *
13 : * Param}tres: *
14 : * *
15 : * e XPHILB Tableau des vecteurs propres aux points de liaisons. *
16 : * e FTOTB Tableau des forces sur base B. *
17 : * s FTOTBA Tableau des forces base B projet{es sur base A. *
18 : * e IBASB Indique dans quelle sous base appartient le point de *
19 : * liaison. *
20 : * e INMSB Nombre de modes par sous base. *
21 : * e IORSB Donne l'indice du premier mode de la sous base dans *
22 : * l'ensemble des modes. *
23 : * e NSB Nombre total de sous base. *
24 : * e NPLSB Nombre total de points intervenant dans les liaisons *
25 : * d'une sous base. *
26 : * e NPLB Nombre total de points intervenant dans les liaisons. *
27 : * e IDIMB Nombre de ddl retenus. *
28 : * e NA1 Nombre total d'inconnues en base A. *
29 : * *
30 : * *
31 : * Auteur, date de cr{ation: *
32 : * *
33 : * Lionel VIVAN, le 22 Septembre 1989. *
34 : * *
35 : *--------------------------------------------------------------------*
36 : *
37 : INTEGER IBASB(*),INMSB(*),IPLSB(*),IORSB(*)
38 : REAL*8 XPHILB(NSB,NPLSB,NA2,*),FTOTB(NPLB,*),FTOTBA(*)
39 : *
40 : DO 10 IN = 1,NA1
41 : FTOTBA(IN) = 0.D0
42 : 10 CONTINUE
43 : * end do
44 : DO 20 IP = 1,NPLB
45 : ISB = IBASB(IP)
46 : IPLB = IPLSB(IP)
47 : NA3 = INMSB(ISB)
48 : INA2 = IORSB(ISB) - 1
49 : DO 30 IN = 1,NA3
50 : XRET = 0.D0
51 : DO 40 ID = 1,IDIMB
52 : XRET = XRET + XPHILB(ISB,IPLB,IN,ID) * FTOTB(IP,ID)
53 : 40 CONTINUE
54 : * end do
55 : IN2 = INA2 + IN
56 : FTOTBA(IN2) = FTOTBA(IN2) + XRET
57 : 30 CONTINUE
58 : * end do
59 : 20 CONTINUE
60 : * end do
61 : *
62 : END
63 :
© Cast3M 2003 - All rights reserved.
Disclaimer