Télécharger la plaquette Cast3M


fr en

Télécharger calc2.eso

Retour à la liste

Numérotation des lignes : 

calc2

  1. C CALC2     SOURCE    CHAT      05/01/12    21:45:37     5004
  2. C CALC2.eso     SOURCE     bald     MAR 28/03/95     00:00:00
  3.         SUBROUTINE CALC2(KNN,KNT,KTN,KTT,SIG0,DU,DV,NSTRS,
  4.      $                   DELTA,RI0,DCON)
  5. C
  6. C-------------------------------------------------------------
  7. C
  8. C   ENTREES
  9. C
  10. C     KNN,KNT,KTN,KTT    = RAIDEURS DU JOINT 2D
  11. C     SIG0(NSTRS)        = CONTRAINTES INITIALES (AU PAS PRECEDENT)
  12. C     DU                 = INCREMENT DE LA U
  13. C     DV                 = INCREMENT DE LA V
  14. C     NSTRS              = NOMBRE DE COMPOSANTES DE CONTRAINTES
  15. C
  16. C   SORTIES
  17. C
  18. C     DELTA              = INCREMENT DE LA U EQUIVALENTE
  19. C     RI0(NSTRS,NSTRS)   = MATRICE DE RIGIDITE DU JOINT
  20. C     DCON(NSTRS)        = VECTEUR DES INCREMENTS DE CONTRAINTE
  21. C
  22. C-------------------------------------------------------------
  23. C
  24. c
  25.       IMPLICIT INTEGER(I-N)
  26.        IMPLICIT REAL*8(A-H,O-Z)
  27.        DIMENSION SIG0(NSTRS), RI0(NSTRS,NSTRS), DCON(NSTRS)
  28.        REAL*8 KNN, KNT, KTN, KTT
  29. c
  30.        TOL1=1.D-20
  31.        TOL2=1.D-8
  32. c
  33. c    calcul de la matrice Dt reale (2D)
  34. c
  35.        IF(ABS(SIG0(1)).LE.TOL1)THEN
  36. c
  37. c    la tau est nulle
  38. c
  39.           IF(ABS(DU).LE.TOL2)THEN
  40. c
  41. c       l' incrément de déformation est nul
  42. c
  43.              DELTA=0.D0
  44. c
  45.              RI0(1,1)=0.D0
  46.              RI0(1,2)=0.D0
  47. c
  48.              RI0(2,1)=0.D0
  49.              RI0(2,2)=KNN
  50. c
  51.           ELSE
  52. c
  53.              DELTA=ABS(DU)
  54. c
  55.              RI0(1,1)=KTT
  56.              RI0(1,2)=KTN*(DU/ABS(DU))
  57. c
  58.              RI0(2,1)=KNT*(DU/ABS(DU))
  59.              RI0(2,2)=KNN
  60. c
  61.           END IF
  62. c
  63.        ELSE
  64. c
  65.        DELTA=DU*(SIG0(1)/ABS(SIG0(1)))
  66. c
  67.        RI0(1,1)=KTT
  68.        RI0(1,2)=KTN*(SIG0(1)/ABS(SIG0(1)))
  69. c
  70.        RI0(2,1)=KNT*(SIG0(1)/ABS(SIG0(1)))
  71.        RI0(2,2)=KNN
  72. c
  73.        END IF
  74. c
  75.        DCON(1)=RI0(1,1)*DU+RI0(1,2)*DV
  76.        DCON(2)=RI0(2,1)*DU+RI0(2,2)*DV
  77. c
  78.        RETURN
  79. c
  80.        END
  81.  
  82.  
  83.  

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