Télécharger rispl2.eso

Retour à la liste

Numérotation des lignes :

  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