Numérotation des lignes :

flua
C FLUA      SOURCE    AM        09/12/08    21:15:14     6582      SUBROUTINE FLUA(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,     \$                       AN,BN,CN,NDI)CC_______________________________________________________________________CC SOUS PROGRAMME QUI CALCUL LES MATRICES DE DEFORMATIONC  DE FLUAGE PROPREC   D APRES LE MODELE AMORTISSEURC_______________________________________________________________________C      IMPLICIT INTEGER(I-N)      IMPLICIT REAL*8(A-H,O-Z)      DIMENSION VARF(NVARI),CMAT(NMAT),EPSFLU(6)      DIMENSION AN(*),BN(NDI,*),CN(NDI,*)C      CALL ZERO(AN,NSTRS,1)      CALL ZERO(BN,NDI,NSTRS)      CALL ZERO(CN,NDI,NSTRS)CC RECUPERATION DES VALEURS DES PARAMETRES MATERIAUC      VIDV = CMAT(23)      TDT  = CMAT(NMAT)      XNU  = CMAT(24)CC RECUPERATION DES VARIABLES INTERNES INITIALESCC    REMPLISSAGE DES VECTEURS DE DEFORMATION DE FLUAGEC      EPSFLU(1) = VARF(45)      EPSFLU(2) = VARF(46)      EPSFLU(3) = VARF(47)      EPSFLU(4) = VARF(48)      EPSFLU(5) = VARF(49)      EPSFLU(6) = VARF(50)CCC CONSTRUCTION DE LA MATRICEC      ADR = 0.D0      BDR = TDT/(2.D0*VIDV)      CDR = TDT/(2.D0*VIDV)C      DO 10 I=1,NSTRS          AN(I) = ADR * EPSFLU(I)  10  CONTINUEC      IF(IFOUR.EQ.-2) THENCC  - CONTRAINTES PLANESC         BN(1,1)=1.D0         BN(1,2)=-XNU         BN(2,1)=BN(1,2)         BN(2,2)=BN(1,1)         BN(3,3)=1.D0+XNU      ELSE IF((IFOUR.EQ.-1).OR.(IFOUR.EQ.0)) THENCC  - DEFORMATION PLANE OU AXISYMETRIQUEC         BN(1,1)=1.D0         BN(1,2)=-XNU         BN(1,3)=BN(1,2)         BN(2,1)=BN(1,2)         BN(2,2)=BN(1,1)         BN(2,3)=BN(1,2)         BN(3,1)=BN(1,2)         BN(3,2)=BN(1,2)         BN(3,3)=BN(1,1)         BN(4,4)=1.D0+XNU      ELSECC  - CALCUL TRIDIMENSIONELC         BN(1,1)=1.D0         BN(1,2)=-XNU         BN(1,3)=BN(1,2)         BN(2,1)=BN(1,2)         BN(2,2)=BN(1,1)         BN(2,3)=BN(1,2)         BN(3,1)=BN(1,2)         BN(3,2)=BN(1,2)         BN(3,3)=BN(1,1)         BN(4,4)=1.D0+XNU         BN(5,5)=BN(4,4)         BN(6,6)=BN(4,4)      ENDIFC      DO 20 I=1,NSTRS        DO 20 J=1,NSTRS          CN(I,J) = CDR * BN(I,J)          BN(I,J) = BDR * BN(I,J)  20  CONTINUEC      END   

© Cast3M 2003 - Tous droits réservés.
Mentions légales