Télécharger normal.eso

Retour à la liste

Numérotation des lignes :

  1. C NORMAL SOURCE CHAT 05/01/13 02:00:11 5004
  2. SUBROUTINE NORMAL(A,B,C,SIGMA,ANOR,ALPHA,SELAS,IBOU,ITYP)
  3. C
  4. C CALCUL DE LA NORMNALE A LA SURFACE DE PLASTICITE
  5. C
  6. IMPLICIT INTEGER(I-N)
  7. IMPLICIT REAL*8(A-H,O-Z)
  8. DIMENSION SIGMA(*),ANOR(*)
  9. C CALCUL DE M*SIGMA
  10. DO 1001 IA=1,IBOU
  11. 1001 ANOR(IA)=0.
  12. GO TO (1,2,2,4,5,6,7,8,9,10,99,99,99,14),ITYP
  13. 1 CONTINUE
  14. ANOR(1)=A*SIGMA(1)-B*(SIGMA(2)+SIGMA(3))
  15. ANOR(2)=A*SIGMA(2)-B*(SIGMA(3)+SIGMA(1))
  16. ANOR(3)=A*SIGMA(3)-B*(SIGMA(1)+SIGMA(2))
  17. DO 1002 IA=1,3
  18. 1002 ANOR(IA+3)=C*SIGMA(3+IA)
  19. GO TO 1000
  20. 2 CONTINUE
  21. ANOR(1)=ALPHA*(SIGMA(1)*A-B*SIGMA(2))
  22. ANOR(2)=ALPHA*(SIGMA(2)*A-B*SIGMA(1))
  23. ANOR(4)= SIGMA(4)*A-B*SIGMA(5)
  24. ANOR(5)= SIGMA(5)*A-B*SIGMA(4)
  25. ANOR(3)=C*SIGMA(3)*ALPHA
  26. ANOR(6)=C*SIGMA(6)
  27. GO TO 1000
  28. 4 ANOR(3)=SIGMA(3)
  29. GO TO 1000
  30. 5 CONTINUE
  31. GO TO 1000
  32. 6 ANOR(1)=A*SIGMA(1)-B*SIGMA(2)
  33. ANOR(2)=A*SIGMA(2)-B*SIGMA(1)
  34. ANOR(4)=C*SIGMA(4)
  35. GO TO 1000
  36. 7 CONTINUE
  37. ANOR(1)=ALPHA*A*SIGMA(1)
  38. ANOR(4)= A*SIGMA(4)
  39. GO TO 1000
  40. 8 ANOR(1)=SIGMA(1)
  41. GO TO 1000
  42. C= Modes de calcul UNIDIMENSIONNELs (1D)
  43. 14 ANOR(1)=A*SIGMA(1)-B*(SIGMA(2)+SIGMA(3))
  44. ANOR(2)=A*SIGMA(2)-B*(SIGMA(3)+SIGMA(1))
  45. ANOR(3)=A*SIGMA(3)-B*(SIGMA(1)+SIGMA(2))
  46. C* GOTO 1000
  47. 9 CONTINUE
  48. 10 CONTINUE
  49. C NORMALISATION
  50. 1000 CONTINUE
  51. X=1./SELAS
  52. Y=0.
  53. DO 1003 IA=1,IBOU
  54. ANOR(IA)=ANOR(IA)*X
  55. 1003 Y=Y+ANOR(IA)*ANOR(IA)
  56. Y=1./SQRT(Y)
  57. DO 1004 IA=1,IBOU
  58. 1004 ANOR(IA)=ANOR(IA)*Y
  59. RETURN
  60. 99 WRITE(6,7999)
  61. 7999 FORMAT('0 NORMAL - CAS NON IMPLEMENTE ' /)
  62. RETURN
  63. END
  64.  
  65.  
  66.  
  67.  
  68.  

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