C D2VFXA SOURCE CHAT 05/01/12 22:35:22 5004 C DEVFXA SOURCE LAVARENN 96/10/30 21:16:44 2349 SUBROUTINE D2VFXA(FEXA,IFEXA,FTOTA,NPFEXA,NA1,NPC1,NPAS, & FTEXB,FEXB,NPLB,IDIMB,RIGIDE) IMPLICIT INTEGER(I-N) IMPLICIT REAL*8(A-H,O-Z) *--------------------------------------------------------------------* * * * Op{rateur DYNE : algorithme de Fu - de Vogelaere * * ________________________________________________ * * * * Ajout des termes de forces ext{rieures libres. * * * * Param}tres: * * * * e FEXA Evolution des forces extérieures, en base A * * e IFEXA Numéro du mode correspondant au point de chargement * * es FTOTA Forces extérieures libres totalisées, en base A * * e NPFEXA Nombre de points de chargement * * e NA1 Nombre total d'inconnues dans la base A * * e NPC1 Nombre de pas de calcul * * e NPAS Num{ro du pas courant * * e FTEXB Evolution des forces extérieures sur base B, pour le * * calcul des moments * * es FEXB Forces extérieures sur base B, pour le calcul des * * moments * * * * * Auteur, date de création: * * * * Denis ROBERT-MOUGIN, le 22 mai 1989. * * * *--------------------------------------------------------------------* REAL*8 FEXA(NPFEXA,NPC1,*),FTOTA(NA1,*) REAL*8 FEXB(NPLB,2,*),FTEXB(NPLB,NPC1,2,*) INTEGER IFEXA(*) * LOGICAL RIGIDE * * Boucle sur les points de chargement: * DO 100 I=1,NPFEXA FTOTA(I,1) = FEXA(I,NPAS,1) 100 CONTINUE * END DO * * Dans le cas des corps rigides, boucle sur les chargements Base B IF (RIGIDE) THEN IF (IDIMB.EQ.6) THEN IDIM = 3 ELSE IDIM = 2 ENDIF DO 200 IP=1,NPLB DO 210 ID=1,IDIM FEXB(IP,1,ID) = FTEXB(IP,NPAS,1,ID) 210 CONTINUE 200 CONTINUE ENDIF RETURN END