Télécharger rispl3.eso

Retour à la liste

Numérotation des lignes :

  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