1 : C DEVFXA SOURCE CHAT 05/01/12 22:46:00 5004
2 : SUBROUTINE DEVFXA(FEXA,IFEXA,FTOTA,NPFEXA,NA1,NPC1,NPAS,
3 : & FTEXB,FEXB,NPLB,IDIMB,RIGIDE)
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 : * Ajout des termes de forces ext{rieures libres. *
12 : * *
13 : * Param}tres: *
14 : * *
15 : * e FEXA Evolution des forces extérieures, en base A *
16 : * e IFEXA Numéro du mode correspondant au point de chargement *
17 : * es FTOTA Forces extérieures libres totalisées, en base A *
18 : * e NPFEXA Nombre de points de chargement *
19 : * e NA1 Nombre total d'inconnues dans la base A *
20 : * e NPC1 Nombre de pas de calcul *
21 : * e NPAS Num{ro du pas courant *
22 : * e FTEXB Evolution des forces extérieures sur base B, pour le *
23 : * calcul des moments *
24 : * es FEXB Forces extérieures sur base B, pour le calcul des *
25 : * moments * *
26 : * *
27 : * Auteur, date de création: *
28 : * *
29 : * Denis ROBERT-MOUGIN, le 22 mai 1989. *
30 : * *
31 : *--------------------------------------------------------------------*
32 : REAL*8 FEXA(NPFEXA,NPC1,*),FTOTA(NA1,*)
33 : REAL*8 FEXB(NPLB,2,*),FTEXB(NPLB,NPC1,2,*)
34 : INTEGER IFEXA(*)
35 : *
36 : LOGICAL RIGIDE
37 : *
38 : * Boucle sur les points de chargement:
39 : *
40 : NPA = NPAS + 1
41 : DO 100 I=1,NPFEXA
42 : FTOTA(I,1) = FEXA(I,NPA,1)
43 : FTOTA(I,2) = FEXA(I,NPA,2)
44 : 100 CONTINUE
45 : * END DO
46 : *
47 : * Dans le cas des corps rigides, boucle sur les chargements Base B
48 : IF (RIGIDE) THEN
49 : IF (IDIMB.EQ.6) THEN
50 : IDIM = 3
51 : ELSE
52 : IDIM = 2
53 : ENDIF
54 : DO 200 IP=1,NPLB
55 : DO 210 ID=1,IDIM
56 : FEXB(IP,1,ID) = FTEXB(IP,NPA,1,ID)
57 : FEXB(IP,2,ID) = FTEXB(IP,NPA,2,ID)
58 : 210 CONTINUE
59 : 200 CONTINUE
60 : ENDIF
61 : RETURN
62 : END
63 :
64 :
© Cast3M 2003 - All rights reserved.
Disclaimer