Télécharger chrp3d.eso

Retour à la liste

Numérotation des lignes :

  1. C CHRP3D SOURCE CHAT 07/10/22 21:15:32 5921
  2. SUBROUTINE CHRP3D(PPP,SIEPIN,SIEPOO,IOP)
  3. IMPLICIT REAL*8(A-H,o-Z)
  4. IMPLICIT INTEGER (I-N)
  5. *
  6. * ------------------------------------------------------------------
  7. *
  8. * CHGT DE REPERE POUR LES CONTRAINTES ET LES DEFORMATIONS
  9. * ECRITES SOUS FORME VECTORIELLE EN 3D
  10. * H. BUNG 02-93
  11. * ------------------------------------------------------------------
  12. *
  13. * PPP : MATRICE DE PASSAGE
  14. * SIEPIN : VECTEUR DE CONTRAINTES OU DE DEFORMATIONS EN ENTREE
  15. * SIEPOO : VECTEUR DE CONTRAINTES OU DE DEFORMATIONS EN SORTIE
  16. * IOP ENTREE SORTIE
  17. * = 0 SIG_G SIG_L
  18. * = 1 SIG_L SIG_G
  19. * = 2 EPS_G EPS_L
  20. * = 3 EPS_L EPS_G
  21. *
  22. *-- RAPPEL
  23. * [SIG] = [ S_11, S_22, S_33, S_12, S_23, S_13 ]
  24. * [EPS] = [ E_11, E_22, E_33, 2*E_12, 2*E_23, 2*E_13 ]
  25. *
  26. *
  27. * IMPLICIT NONE
  28. *
  29. *
  30. *--- VARIABLES GLOBALES
  31. * REAL *8 PPP(3,3),SIEPIN(*),SIEPOO(*)
  32. * INTEGER IOP
  33. DIMENSION PPP(3,*),SIEPIN(*),SIEPOO(*)
  34. *
  35. *--- VARIABLES LOCALES
  36.  
  37. DIMENSION AIN(9),AOO(9)
  38. DIMENSION LLM(9),LLV(6),LDEF(6)
  39. *
  40. DATA LLM/1,4,6, 4,2,5, 6,5,3/
  41. DATA LLV/1,5,9, 2,6,3/
  42. DATA LDEF/2,3, 4,6, 7,8/
  43. *
  44. *
  45. *---- PASSAGE VERS LA FORME MATRICIELLE
  46. DO IK=1,9
  47. AIN(IK) = SIEPIN(LLM(IK))
  48. END DO
  49. *
  50. IF(IOP.GE.2) THEN
  51. *--- TENSEUR DE DEFORMATIONS
  52. DO IK=1,6
  53. IKK = LDEF(IK)
  54. AIN(IKK) = 0.5D0*AIN(IKK)
  55. ENDDO
  56. ENDIF
  57. *
  58. KOP=MOD(IOP,2)
  59. IF(KOP.EQ.0)THEN
  60. * ENTREE : GLOBAL --> SORTIE : LOCAL
  61. CALL DR3GLM1(PPP,AIN,AOO)
  62. ENDIF
  63. IF(KOP.EQ.1) THEN
  64. * ENTREE : LOCAL --> SORTIE : GLOBAL
  65. CALL DR3GLM2(PPP,AOO,AIN)
  66. ENDIF
  67. IF(KOP.NE.0.AND.KOP.NE.1)THEN
  68. * STOP
  69. STOP 'CHR3D_1'
  70. ENDIF
  71. C SELECT CASE(KOP)
  72. C CASE (0)
  73. C* ENTREE : GLOBAL --> SORTIE : LOCAL
  74. C CALL DR3GLM1(PPP,AIN,AOO)
  75. C CASE (1)
  76. C* ENTREE : LOCAL --> SORTIE : GLOBAL
  77. C CALL DR3GLM2(PPP,AOO,AIN)
  78. C CASE DEFAULT
  79. C* STOP
  80. C STOP 'CHR3D_1'
  81. C END SELECT
  82. *
  83. *---- PASSAGE VERS LA FORME VECTORIELLE
  84. DO IK=1,6
  85. SIEPOO(IK) = AOO(LLV(IK))
  86. END DO
  87. *
  88. *---- POUR LES DEFORMATIONS
  89. IF(IOP.GE.2) THEN
  90. SIEPOO(4) = 2.D0*SIEPOO(4)
  91. SIEPOO(5) = 2.D0*SIEPOO(5)
  92. SIEPOO(6) = 2.D0*SIEPOO(6)
  93. ENDIF
  94. *
  95. END
  96.  
  97.  
  98.  

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