Télécharger diaxyz.eso

Retour à la liste

Numérotation des lignes :

  1. C DIAXYZ SOURCE CHAT 05/01/12 22:50:02 5004
  2. SUBROUTINE DIAXYZ(ITYP,SIGT,SIGEL,VECPRO,IBO,A,B,C,B1,B2)
  3. C
  4. C EFFECTUE LA TRANSFORMATION DE SIGAMA DIAGONALISE EN
  5. C SIGAM GLOBAL
  6. C
  7. C ENTREE SIGT
  8. C SORTIE SIGEL
  9. C ITYP TYPE DE L'ELEMENT
  10. C
  11. IMPLICIT INTEGER(I-N)
  12. IMPLICIT REAL*8(A-H,O-Z)
  13. DIMENSION SIGT(*),SIGEL(*),VECPRO(*)
  14. C
  15. GO TO (70,71,71,73,74,75,76,77,78,79,80,81,71),ITYP
  16. C
  17. C MODIF DU 05-01-96
  18. C POUR ITYP=13 (COQUES EPAISSES) ON PASSAIT ICI (AU 83)
  19. C DESORMAIS ON VA AU 71 COMME LES CALCULS EN CONTRAINTE PLANE
  20. C
  21. C 83 CONTINUE
  22. C COQUES AVEC CISAILLEMENTS TRANSVERSES
  23. C SIGEL(1)=B*(SIGT(1)-SIGT(2))
  24. C SIGEL(2)=B*(SIGT(1)+SIGT(2))
  25. C SIGEL(3)=SIGT(3)
  26. C SIGEL(4)=SIGT(4)
  27. C SIGEL(5)=SIGT(5)
  28. C SIGEL(6)=SIGT(6)
  29. C GO TO 740
  30. 81 CONTINUE
  31. C TUYAUX
  32. 80 CONTINUE
  33. C POUTRES
  34. DO 810 I=1,6
  35. 810 SIGEL(I)=SIGT(I)
  36. GO TO 740
  37. 79 CONTINUE
  38. GO TO 740
  39. 78 CONTINUE
  40. GO TO 740
  41. 77 SIGEL(1)=SIGT(1)
  42. GO TO 740
  43. 76 DO 760 IA=1,4
  44. 760 SIGEL(IA)=SIGT(IA)
  45. GO TO 740
  46. 75 CONTINUE
  47. SIGEL(1)=B*(SIGT(1)-SIGT(2))
  48. SIGEL(2)=B*(SIGT(1)+SIGT(2))
  49. SIGEL(4)=SIGT(4)
  50. SIGEL(3)=SIGT(3)
  51. GO TO 740
  52. 74 CONTINUE
  53. DO 7401 IB=1,3
  54. SIGEL(IB)=0.
  55. DO 7402 IC=1,3
  56. IP=3*(IC-1)+IB
  57. 7402 SIGEL(IB)=SIGEL(IB)+VECPRO(IP)*SIGT(IC)
  58. 7401 SIGEL(3+IB)=SIGT(3+IB)
  59. GO TO 740
  60. 71 CONTINUE
  61. B1=B
  62. B2=B
  63. DO 710 I=1,IBO
  64. IP=(I-1)*3
  65. SIGEL(IP+1)=B1*SIGT(IP+1)-B2*SIGT(IP+2)
  66. SIGEL(IP+2)=B2*SIGT(IP+1)+B1*SIGT(IP+2)
  67. 710 SIGEL(3+IP)=SIGT(3+IP)
  68. GO TO 740
  69. 73 CONTINUE
  70. SIGEL(3)=SIGT(3)
  71. GO TO 740
  72. 70 CONTINUE
  73. SIGT(1)=SIGT(1)*A
  74. SIGT(2)=SIGT(2)*B
  75. SIGT(3)=SIGT(3)*C
  76. SIGEL(1)=SIGT(1)+SIGT(2)+SIGT(3)
  77. SIGEL(2)=SIGT(1)-SIGT(2)+SIGT(3)
  78. SIGEL(3)=SIGT(1)-2.*SIGT(3)
  79. SIGEL(4)=SIGT(4)
  80. SIGEL(5)=SIGT(5)
  81. SIGEL(6)=SIGT(6)
  82. 740 CONTINUE
  83. RETURN
  84. END
  85.  
  86.  
  87.  

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