devfxa
C DEVFXA SOURCE BP208322 18/12/20 21:15:31 10048 * SUBROUTINE DEVFXA(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 * * (supprime le 2018-12-14 par bp) * 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,*) c INTEGER IFEXA(*) * LOGICAL RIGIDE * * Boucle sur les points de chargement: * NPA = NPAS + 1 DO 100 I=1,NPFEXA FTOTA(I,1) = FEXA(I,NPA,1) FTOTA(I,2) = FEXA(I,NPA,2) 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,NPA,1,ID) FEXB(IP,2,ID) = FTEXB(IP,NPA,2,ID) 210 CONTINUE 200 CONTINUE ENDIF RETURN END
© Cast3M 2003 - Tous droits réservés.
Mentions légales