Numérotation des lignes :

C CALC3     SOURCE    CHAT      05/01/12    21:45:39     5004C CALC3.eso     SOURCE     bald     MAR 28/03/95     00:00:00        SUBROUTINE CALC3(KNN,KNT,KTN,KTT,SIG0,DU1,DU2,DV,NSTRS,     \$                   DELTA,RI0,DCON)CC-------------------------------------------------------------CC   ENTREESCC     KNN,KNT,KTN,KTT   = RAIDEURS DU JOINT 2D EQUIVALENTC     SIG0(NSTRS)       = CONTRAINTES INITIALES (AU PAS PRECEDENT)C     DU1               = INCREMENT DE LA U1C     DU2               = INCREMENT DE LA U2C     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 (3D)c       IF((ABS(SIG0(1)).LE.TOL1).AND.(ABS(SIG0(2)).LE.TOL1))THENcc    les deux tau sont nullesc          IF((ABS(DU1).LE.TOL2).AND.(ABS(DU2).LE.TOL2))THENcc       les deux incréments de déformation sont nulsc             DELTA=0.D0c             RI0(1,1)=0.D0             RI0(1,2)=0.D0             RI0(1,3)=0.D0c             RI0(2,1)=0.D0             RI0(2,2)=0.D0             RI0(2,3)=0.D0c             RI0(3,1)=0.D0             RI0(3,2)=0.D0             RI0(3,3)=KNNc          ELSEc             DELTA=SQRT(DU1**2+DU2**2)c             RI0(1,1)=KTT             RI0(1,2)=0.D0             RI0(1,3)=KTN*(DU1/SQRT(DU1**2+DU2**2))c             RI0(2,1)=0.D0             RI0(2,2)=KTT             RI0(2,3)=KTN*(DU2/SQRT(DU1**2+DU2**2))c             RI0(3,1)=KNT*(DU1/SQRT(DU1**2+DU2**2))             RI0(3,2)=KNT*(DU2/SQRT(DU1**2+DU2**2))             RI0(3,3)=KNNc          END IFc       ELSEc       DELTA=(SIG0(1)*DU1+SIG0(2)*DU2)/SQRT(SIG0(1)**2+SIG0(2)**2)c       RI0(1,1)=KTT*(1.D0+(SIG0(2)**2)/(SIG0(1)**2+SIG0(2)**2))       RI0(1,2)=-KTT*(SIG0(1)*SIG0(2))/(SIG0(1)**2+SIG0(2)**2)       RI0(1,3)=KTN*SIG0(1)/SQRT(SIG0(1)**2+SIG0(2)**2)c       RI0(2,1)=-KTT*(SIG0(2)*SIG0(1))/(SIG0(1)**2+SIG0(2)**2)       RI0(2,2)=KTT*(1.D0+(SIG0(1)**2)/(SIG0(1)**2+SIG0(2)**2))       RI0(2,3)=KTN*SIG0(2)/SQRT(SIG0(1)**2+SIG0(2)**2)c       RI0(3,1)=KNT*SIG0(1)/SQRT(SIG0(1)**2+SIG0(2)**2)       RI0(3,2)=KNT*SIG0(2)/SQRT(SIG0(1)**2+SIG0(2)**2)       RI0(3,3)=KNNc       END IFc       DCON(1)=RI0(1,1)*DU1+RI0(1,2)*DU2+RI0(1,3)*DV       DCON(2)=RI0(2,1)*DU1+RI0(2,2)*DU2+RI0(2,3)*DV       DCON(3)=RI0(3,1)*DU1+RI0(3,2)*DU2+RI0(3,3)*DVc       RETURNc       END

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