Numérotation des lignes :

calc2
C CALC2     SOURCE    CHAT      05/01/12    21:45:37     5004C CALC2.eso     SOURCE     bald     MAR 28/03/95     00:00:00        SUBROUTINE CALC2(KNN,KNT,KTN,KTT,SIG0,DU,DV,NSTRS,     \$                   DELTA,RI0,DCON)CC-------------------------------------------------------------CC   ENTREESCC     KNN,KNT,KTN,KTT    = RAIDEURS DU JOINT 2DC     SIG0(NSTRS)        = CONTRAINTES INITIALES (AU PAS PRECEDENT)C     DU                 = INCREMENT DE LA UC     DV                 = INCREMENT DE LA VC     NSTRS              = NOMBRE DE COMPOSANTES DE CONTRAINTESCC   SORTIESCC     DELTA              = INCREMENT DE LA U EQUIVALENTEC     RI0(NSTRS,NSTRS)   = MATRICE DE RIGIDITE DU JOINTC     DCON(NSTRS)        = VECTEUR DES INCREMENTS DE CONTRAINTECC-------------------------------------------------------------Cc      IMPLICIT INTEGER(I-N)       IMPLICIT REAL*8(A-H,O-Z)       DIMENSION SIG0(NSTRS), RI0(NSTRS,NSTRS), DCON(NSTRS)       REAL*8 KNN, KNT, KTN, KTTc       TOL1=1.D-20       TOL2=1.D-8cc    calcul de la matrice Dt reale (2D)c       IF(ABS(SIG0(1)).LE.TOL1)THENcc    la tau est nullec          IF(ABS(DU).LE.TOL2)THENcc       l' incrément de déformation est nulc             DELTA=0.D0c             RI0(1,1)=0.D0             RI0(1,2)=0.D0c             RI0(2,1)=0.D0             RI0(2,2)=KNNc          ELSEc             DELTA=ABS(DU)c             RI0(1,1)=KTT             RI0(1,2)=KTN*(DU/ABS(DU))c             RI0(2,1)=KNT*(DU/ABS(DU))             RI0(2,2)=KNNc          END IFc       ELSEc       DELTA=DU*(SIG0(1)/ABS(SIG0(1)))c       RI0(1,1)=KTT       RI0(1,2)=KTN*(SIG0(1)/ABS(SIG0(1)))c       RI0(2,1)=KNT*(SIG0(1)/ABS(SIG0(1)))       RI0(2,2)=KNNc       END IFc       DCON(1)=RI0(1,1)*DU+RI0(1,2)*DV       DCON(2)=RI0(2,1)*DU+RI0(2,2)*DVc       RETURNc       END   

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