Télécharger courbu.eso

Retour à la liste

Numérotation des lignes :

  1. C COURBU SOURCE GOUNAND 14/05/28 21:15:04 8056
  2. SUBROUTINE COURBU(SIG,FT,ISGN1,ISGN2,NSURF,FFF)
  3. C Calcul de la derivee seconde pour le calcul de la matrice tangente
  4. C 21/08/95
  5. C
  6. IMPLICIT INTEGER(I-N)
  7. IMPLICIT REAL*8(A-H,O-Z)
  8. C
  9. C========================================================================
  10. C
  11. C ENTREES:
  12. C -------
  13. C SIG = Vecteurs contraintes egal a SIGT pour les deux surfaces
  14. C ISGN1 = 1 EN TRACTION -1 EN COMPRESSION
  15. C ISGN2 = -1 SI ON ECOULE SUR LA PLUS PETITE CONT. PRIN.
  16. C +1 SI ON ECOULE SUR LA PLUS GRANDE CONT. PRINC.
  17. C FT = Contrainte de reference pour le test de l'apex
  18. C SORTIES:
  19. C -------
  20. C FFF = Derivee seconde (matrice 3*3)
  21. C
  22. C
  23. C========================================================================
  24. C
  25. PARAMETER(XZER=0.D0,UNDEMI=.5D0,UNQUAR=.25D0,UN=1.D0,
  26. & DEUX =2.D0,PRECIS=1.D-4)
  27. C
  28. DIMENSION SIG(4),XPP(3,3),FFF(3,3),XPSIG(3,1),XPSIGT(1,3)
  29. C
  30. CALL ZERO(XPP,3,3)
  31. CALL ZERO(FFF,3,3)
  32. XPP(1,1) = UNDEMI
  33. XPP(1,2) = (-1.D0)*UNDEMI
  34. XPP(2,1) = (-1.D0)*UNDEMI
  35. XPP(2,2) = UNDEMI
  36. XPP(3,3) = DEUX
  37. SIG(3) = SIG(4)
  38. C========================================================================
  39. C Calcul de phi
  40. C====================================================================
  41. PHI = SQRT ((UNDEMI*(SIG(1)-SIG(2)))**2+SIG(3)**2)
  42. IF (PHI.LT.ABS(PRECIS*FT)) THEN
  43. C Cas de l'apex
  44. IF (NSURF.EQ.2) THEN
  45. PHI = (ISGN2*ISGN1)*(UN/ABS(PRECIS*FT))
  46. CALL PLMATR(XPP,UNDEMI*PHI,FFF,XZER,3,3,FFF)
  47. ELSE
  48. CALL ZERO(FFF,3,3)
  49. ENDIF
  50. ELSE
  51. C Cas regulier
  52. CALL MATVE1(XPP,SIG,3,3,XPSIG,2)
  53. CALL MATVE1(XPP,SIG,3,3,XPSIGT,1)
  54. CALL MATMAT(XPSIG,XPSIGT,3,1,3,FFF)
  55. PHI = (ISGN1*ISGN2)*(UN/PHI)
  56. CALL PLMATR(XPP,UNDEMI*PHI,FFF,-UNQUAR*(PHI**3),3,3,FFF)
  57. ENDIF
  58. C========================================================================
  59. C Fin de la routine
  60. C========================================================================
  61. RETURN
  62. END
  63.  
  64.  
  65.  
  66.  
  67.  
  68.  

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