Télécharger la plaquette Cast3M


fr en

Télécharger rispl2.eso

Retour à la liste

Numérotation des lignes : 

rispl2

  1. C RISPL2    SOURCE    CHAT      05/01/13    03:00:54     5004
  2. C RISPL2.eso     SOURCE     bald     MAR 28/03/95     00:00:00
  3.         SUBROUTINE RISPL2(RI0,SIG0,DCON,DU,NSTRS,SGMT,DUp)
  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      DU                 = INCREMENT DE LA U
  13. C      NSRTS              = NOMBRE DE COMPOSANTES DE CONTRAINTES
  14. C      SGMT               = CONTRAINTE LIMITE NORMALE
  15. C
  16. C   SORTIES
  17. C
  18. C      DUp                = INCREMENT PLASTIQUE DE LA U
  19. C
  20. C-------------------------------------------------------------
  21. C
  22. c
  23.       IMPLICIT INTEGER(I-N)
  24.        IMPLICIT REAL*8(A-H,O-Z)
  25.        DIMENSION RI0(NSTRS,NSTRS), DCON(NSTRS), SIG0(NSTRS)
  26. c
  27.        TOL1=ABS(1.D-3*SGMT)
  28. c
  29.        DETERM=0.D0
  30.        DETERM=DETERM+RI0(1,1)*RI0(2,2)
  31.        DETERM=DETERM-RI0(1,2)*RI0(2,1)
  32. c
  33.        IF(ABS(SIG0(2)).LT.TOL1)THEN
  34. c
  35. c    le detérminant est nul
  36. c
  37.        DUe=0.D0
  38. c
  39.        ELSE
  40. c
  41.        DETER=0.D0
  42.        DETER=DETER+DCON(1)*RI0(2,2)
  43.        DETER=DETER-DCON(2)*RI0(1,2)
  44. c
  45.        DUe=DETER/DETERM
  46. c
  47.        END IF
  48. c
  49.        DUp=DU-DUe
  50. c
  51.        RETURN
  52. c
  53.        END
  54.  
  55.  
  56.  
  57.  

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