Télécharger la plaquette Cast3M


fr en

Télécharger rispl3.eso

Retour à la liste

Numérotation des lignes : 

rispl3

  1. C RISPL3    SOURCE    CHAT      05/01/13    03:00:58     5004
  2. C RISPL3.eso     SOURCE     bald     MAR 28/03/95     00:00:00
  3.         SUBROUTINE RISPL3(RI0,SIG0,DCON,DU1,DU2,NSTRS,SGMT,DU1p,DU2p)
  4. C
  5. C-------------------------------------------------------------
  6. C
  7. C   ENTREES
  8. C
  9. C      RI0(NSTRS,NSTRS)   = MATRICE DE RIGIDITE ELASTIQUE DU JOINT
  10. C      SIG0(NSTRS)        = CONTRAINTES INITIALES (AU PAS PRECEDENT)
  11. C      DCON(NSTRS)        = VECTEUR DES INCREMENTS DE CONTRAINTE
  12. C      DU1                = INCREMENT DE LA U1
  13. C      DU2                = INCREMENT DE LA U2
  14. C      NSRTS              = NOMBRE DE COMPOSANTES DE CONTRAINTES
  15. C      SGMT               = CONTRAINTE LIMITE NORMALE
  16. C
  17. C   SORTIES
  18. C
  19. C      DU1p               = INCREMENT PLASTIQUE DE LA DU1
  20. C      DU2p               = INCREMENT PLASTIQUE DE LA DU2
  21. C
  22. C-------------------------------------------------------------
  23. C
  24. c
  25.       IMPLICIT INTEGER(I-N)
  26.        IMPLICIT REAL*8(A-H,O-Z)
  27.        DIMENSION RI0(NSTRS,NSTRS), DCON(NSTRS), SIG0(NSTRS)
  28. c
  29.        TOL1=ABS(1.D-3*SGMT)
  30. c
  31.        DETERM=0.D0
  32.        DETERM=DETERM+RI0(1,1)*RI0(2,2)*RI0(3,3)
  33.        DETERM=DETERM+RI0(1,2)*RI0(2,3)*RI0(3,1)
  34.        DETERM=DETERM+RI0(1,3)*RI0(2,1)*RI0(3,2)
  35.        DETERM=DETERM-RI0(1,3)*RI0(2,2)*RI0(3,1)
  36.        DETERM=DETERM-RI0(1,2)*RI0(2,1)*RI0(3,3)
  37.        DETERM=DETERM-RI0(1,1)*RI0(3,2)*RI0(2,3)
  38. c
  39.        IF(ABS(SIG0(3)).LT.TOL1)THEN
  40. c
  41. c    le detérminant est nul
  42. c
  43.        DU1e=0.D0
  44.        DU2e=0.D0
  45. c
  46.        ELSE
  47. c
  48.        DETER1=0.D0
  49.        DETER1=DETER1+DCON(1)*RI0(2,2)*RI0(3,3)
  50.        DETER1=DETER1+RI0(1,2)*RI0(2,3)*DCON(3)
  51.        DETER1=DETER1+RI0(1,3)*DCON(2)*RI0(3,2)
  52.        DETER1=DETER1-RI0(1,3)*RI0(2,2)*DCON(3)
  53.        DETER1=DETER1-RI0(1,2)*DCON(2)*RI0(3,3)
  54.        DETER1=DETER1-DCON(1)*RI0(3,2)*RI0(2,3)
  55. c
  56.        DETER2=0.D0
  57.        DETER2=DETER2+RI0(1,1)*DCON(2)*RI0(3,3)
  58.        DETER2=DETER2+DCON(1)*RI0(2,3)*RI0(3,1)
  59.        DETER2=DETER2+RI0(1,3)*RI0(2,1)*DCON(3)
  60.        DETER2=DETER2-RI0(1,3)*DCON(2)*RI0(3,1)
  61.        DETER2=DETER2-DCON(1)*RI0(2,1)*RI0(3,3)
  62.        DETER2=DETER2-RI0(1,1)*DCON(3)*RI0(2,3)
  63. c
  64.        DU1e=DETER1/DETERM
  65.        DU2e=DETER2/DETERM
  66. c
  67.        END IF
  68. c
  69.        DU1p=DU1-DU1e
  70.        DU2p=DU2-DU2e
  71. c
  72.        RETURN
  73. c
  74.        END
  75.  
  76.  
  77.  
  78.  

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