Télécharger dohmas.eso

Retour à la liste

Numérotation des lignes :

  1. C DOHMAS SOURCE CHAT 05/01/12 22:53:45 5004
  2.  
  3. C=======================================================================
  4. C= D O H M A S =
  5. C= ----------- =
  6. C= =
  7. C= Fonction : =
  8. C= ---------- =
  9. C= Calcul de la matrice de Hooke (elasticite lineaire) dans le cas =
  10. C= des materiaux isotropes =
  11. C= =
  12. C= Parametres : (E)=Entree (S)=Sortie =
  13. C= ------------ =
  14. C= VELA (E) Valeurs des parametres elastiques du materiau =
  15. C= MATE (E) Chaine contenant le type du materiau =
  16. C= IFOU (E) Mode de calcul (cf. IFOUR de CCOPTIO) =
  17. C= LHOOK (S) Dimension de la matrice de Hooke =
  18. C= KCAS (E) =1 si on veut la matrice de Hooke associee au mode =
  19. C= de calcul (non necessairement inversible) =
  20. C= =2 si on veut la matrice de Hooke complete pour =
  21. C= l'inverser par la suite (matrice de souplesse) =
  22. C= DDHOOK (S) Matrice de Hooke dans le repere GLOBAL =
  23. C= IRET (S) =1 si option existante , =0 sinon =
  24. C=======================================================================
  25.  
  26. SUBROUTINE DOHMAS (VELA,MATE,IFOU,LHOOK,KCAS,DDHOOK,IRET)
  27.  
  28. IMPLICIT INTEGER(I-N)
  29. IMPLICIT REAL*8 (A-H,O-Z)
  30.  
  31. PARAMETER (XUnDemi=0.5,XUn=1.)
  32.  
  33. DIMENSION VELA(*),DDHOOK(LHOOK,*)
  34. CHARACTER*(*) MATE
  35.  
  36. IRET=1
  37. C= ERREUR si le materiau n'est pas massif ISOTROPE ou POREUX (isotrope)
  38. IF (MATE.NE.'ISOTROPE'.AND.MATE.NE.'POREUX ') THEN
  39. IRET=0
  40. RETURN
  41. ENDIF
  42.  
  43. C= Initialisation de la matrice de Hooke
  44. CALL ZERO(DDHOOK,LHOOK,LHOOK)
  45.  
  46. C= Recuperation des module d'Young et coefficient de Poisson
  47. YOU=VELA(1)
  48. XNU=VELA(2)
  49.  
  50. C= -> Contraintes planes avec KCAS=1
  51. IF (IFOU.EQ.-2.AND.KCAS.EQ.1) THEN
  52. AUX=YOU/(XUn-XNU*XNU)
  53. AUX1=AUX*XNU
  54. DDHOOK(1,1)=AUX
  55. DDHOOK(1,2)=AUX1
  56. DDHOOK(2,1)=AUX1
  57. DDHOOK(2,2)=AUX
  58. DDHOOK(4,4)=XUnDemi*YOU/(XUn+XNU)
  59.  
  60. C= -> Deformations planes, Axisymetrie et Contraintes planes avec KCAS=2
  61. ELSE IF (IFOU.EQ.-1.OR.IFOU.EQ.0.OR.IFOU.EQ.-3.OR.
  62. . (IFOU.EQ.-2.AND.KCAS.EQ.2)) THEN
  63. AUX0=YOU/((XUn+XNU)*(XUn-XNU-XNU))
  64. AUX=AUX0*(XUn-XNU)
  65. AUX1=AUX0*XNU
  66. DDHOOK(1,1)=AUX
  67. DDHOOK(2,1)=AUX1
  68. DDHOOK(3,1)=AUX1
  69. DDHOOK(1,2)=AUX1
  70. DDHOOK(2,2)=AUX
  71. DDHOOK(3,2)=AUX1
  72. DDHOOK(1,3)=AUX1
  73. DDHOOK(2,3)=AUX1
  74. DDHOOK(3,3)=AUX
  75. DDHOOK(4,4)=XUnDemi*YOU/(XUn+XNU)
  76.  
  77. C= -> Serie de Fourier et 3D
  78. ELSE IF (IFOU.EQ.1.OR.IFOU.EQ.2) THEN
  79. AUX0=YOU/((XUn+XNU)*(XUn-XNU-XNU))
  80. AUX=AUX0*(XUn-XNU)
  81. AUX1=AUX0*XNU
  82. GEGE=XUnDemi*YOU/(XUn+XNU)
  83. DDHOOK(1,1)=AUX
  84. DDHOOK(2,1)=AUX1
  85. DDHOOK(3,1)=AUX1
  86. DDHOOK(1,2)=AUX1
  87. DDHOOK(2,2)=AUX
  88. DDHOOK(3,2)=AUX1
  89. DDHOOK(1,3)=AUX1
  90. DDHOOK(2,3)=AUX1
  91. DDHOOK(3,3)=AUX
  92. DDHOOK(4,4)=GEGE
  93. DDHOOK(5,5)=GEGE
  94. DDHOOK(6,6)=GEGE
  95.  
  96. C= -> Modes de calcul 1D
  97. ELSE IF (IFOU.GE.3.AND.IFOU.LE.15) THEN
  98. IF (IFOU.EQ.6.AND.KCAS.EQ.1) THEN
  99. DDHOOK(1,1)=YOU
  100. ELSE IF ((IFOU.EQ.4.OR.IFOU.EQ.8).AND.KCAS.EQ.1) THEN
  101. AUX=YOU/(XUn-XNU*XNU)
  102. AUX1=AUX*XNU
  103. DDHOOK(1,1)=AUX
  104. DDHOOK(2,1)=AUX1
  105. DDHOOK(1,2)=AUX1
  106. DDHOOK(2,2)=AUX
  107. ELSE IF ((IFOU.EQ.5.OR.IFOU.EQ.10.OR.IFOU.EQ.13).AND.
  108. . KCAS.EQ.1) THEN
  109. AUX=YOU/(XUn-XNU*XNU)
  110. AUX1=AUX*XNU
  111. DDHOOK(1,1)=AUX
  112. DDHOOK(3,1)=AUX1
  113. DDHOOK(1,3)=AUX1
  114. DDHOOK(3,3)=AUX
  115. ELSE
  116. AUX0=YOU/((XUn+XNU)*(XUn-XNU-XNU))
  117. AUX=AUX0*(XUn-XNU)
  118. AUX1=AUX0*XNU
  119. DDHOOK(1,1)=AUX
  120. DDHOOK(2,1)=AUX1
  121. DDHOOK(3,1)=AUX1
  122. DDHOOK(1,2)=AUX1
  123. DDHOOK(2,2)=AUX
  124. DDHOOK(3,2)=AUX1
  125. DDHOOK(1,3)=AUX1
  126. DDHOOK(2,3)=AUX1
  127. DDHOOK(3,3)=AUX
  128. ENDIF
  129.  
  130. C= -> Cas non prevus actuellement
  131. ELSE
  132. IRET=0
  133. ENDIF
  134.  
  135. RETURN
  136. END
  137.  
  138.  
  139.  
  140.  

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