Télécharger dou266.eso

Retour à la liste

Numérotation des lignes :

  1. C DOU266 SOURCE AM 12/02/27 23:39:37 7289
  2. SUBROUTINE DOU266(VELA,MATE,IFOU,LHOOK,DDHOOK,IRET)
  3. C
  4. C=======================================================================
  5. C
  6. C MATRICE DE HOOKE DES ELEMENTS DE ZONES COHESIVES ZCO2, ZCO3, ZCO4
  7. C
  8. C ENTREES
  9. C VELA() = materiau dans un tableau de travail
  10. C MATE = Nom du materiau
  11. C IFOU = num{ro d'harmonique de fourier: IFOUR de CCOPTIO
  12. C LHOOK = taille de la matrice de hooke
  13. C
  14. C SORTIES
  15. C DDHOOK(LHOOK,LHOOK) = matrice de hooke
  16. C IRET = 1 si option existante 0 SINON
  17. C
  18. C=======================================================================
  19. C
  20. IMPLICIT INTEGER(I-N)
  21. IMPLICIT REAL*8(A-H,O-Z)
  22. CHARACTER*8 MATE
  23. C
  24. DIMENSION VELA(*),DDHOOK(LHOOK,*)
  25. DATA DEUX/2.0D0/
  26. C
  27. IRET=1
  28. CALL ZERO(DDHOOK,LHOOK,LHOOK)
  29. C
  30. C MATERIAU ISOTROPE
  31. C
  32. IF(MATE.EQ.'ZONE_COH') THEN
  33. C
  34. C CAS TRIDIMENSIONNEL
  35. C
  36. IF(IFOU.EQ.2) THEN
  37. ZK1 =VELA(1)
  38. ZK2 =VELA(1)
  39. ZK3 =VELA(2)
  40. C
  41. DDHOOK(1,1)=ZK1
  42. DDHOOK(2,2)=ZK2
  43. DDHOOK(3,3)=ZK3
  44. C
  45. C CAS CONTRAINTES PLANES DANS LE PLAN (S1,N)
  46. C DEFORMATIONS PLANES DANS LE PLAN (S1,N)
  47. C AXISYMETRIE DANS LE PLAN (S1,N)
  48. C
  49. ELSE IF(IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.
  50. & IFOU.EQ.-1.OR.IFOU.EQ.0) THEN
  51. ZK1 =VELA(1)
  52. ZK2 =VELA(2)
  53. C
  54. DDHOOK(1,1)=ZK1
  55. DDHOOK(2,2)=ZK2
  56. C
  57. C OPTION NON DEFINIE
  58. C
  59. ELSE
  60. IRET=0
  61. ENDIF
  62. C
  63. C AUTRES CAS : NON ENCORE PREVUS ...!
  64. C
  65. ELSE
  66. IRET=0
  67. ENDIF
  68. C
  69. RETURN
  70. END
  71.  
  72.  
  73.  
  74.  
  75.  

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