Télécharger calc2.eso

Retour à la liste

Numérotation des lignes :

  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