Télécharger normbl.eso

Retour à la liste

Numérotation des lignes :

  1. C NORMBL SOURCE CHAT 05/01/13 02:00:18 5004
  2. SUBROUTINE NORMBL ( AAA , FT, ISGN1 , ISGN2 , NSURF, FFF)
  3. C Calcul du gradient pour les differents mecanismes
  4. C 13/06/95
  5. C
  6. IMPLICIT INTEGER(I-N)
  7. IMPLICIT REAL*8(A-H,O-Z)
  8. C
  9. C========================================================================
  10. C Cette routine calcule le gradient (ou la normale)
  11. C
  12. C ENTREES:
  13. C -------
  14. C AAA = Vecteurs contraintes egal a SIGT pour les deux surfaces
  15. C ISGN1 = 1 EN TRACTION -1 EN COMPRESSION
  16. C ISGN2 = -1 SI ON ECOULE SUR LA PLUS PETITE CONT. PRIN.
  17. C +1 SI ON ECOULE SUR LA PLUS GRANDE CONT. PRINC.
  18. C FT = Contrainte de reference pour le test de l'apex
  19. C SORTIES:
  20. C -------
  21. C FFF = Vecteur normal
  22. C
  23. C
  24. C========================================================================
  25. C
  26. PARAMETER(UNDEMI=.5D0, DEUX =2.D0 , PRECIS=1.D-4)
  27. C
  28. DIMENSION AAA(4) , FFF(4)
  29. C
  30. C========================================================================
  31. C Calcul de phi
  32. C====================================================================
  33. PHI = SQRT ((UNDEMI*(AAA(1)-AAA(2)))**2+AAA(4)**2)
  34. IF (PHI.LT.ABS(PRECIS*FT)) THEN
  35. C Cas de l'apex
  36. FFF(4) = 0.D0
  37. IF (NSURF.EQ.2) THEN
  38. FFF(1) = ISGN1*UNDEMI*SQRT(2.D0)
  39. FFF(2) = FFF(1)
  40. ELSE
  41. IF (AAA(1).GE.AAA(2)) THEN
  42. FFF(1) = ISGN1*(ISGN2 + 1.)*UNDEMI
  43. FFF(2) = ISGN1*(-ISGN2 + 1.)*UNDEMI
  44. ELSE
  45. FFF(1) = ISGN1*(-ISGN2 + 1.)*UNDEMI
  46. FFF(2) = ISGN1*(ISGN2 + 1.)*UNDEMI
  47. ENDIF
  48. ENDIF
  49. ELSE
  50. C Cas regulier
  51. FFF(1) = ISGN1*(ISGN2*UNDEMI/PHI*UNDEMI*(AAA(1)-AAA(2))
  52. & + UNDEMI)
  53. FFF(2) = ISGN1*(ISGN2*UNDEMI/PHI*UNDEMI*(AAA(2)-AAA(1))
  54. & + UNDEMI)
  55. FFF(4) = ISGN1*ISGN2*AAA(4)/PHI
  56. ENDIF
  57. C========================================================================
  58. C Fin de la routine
  59. C========================================================================
  60. RETURN
  61. END
  62.  
  63.  
  64.  

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