Télécharger tte10.eso

Retour à la liste

Numérotation des lignes :

tte10
  1. C TTE10 SOURCE OF166741 23/12/04 21:15:22 11800
  2.  
  3. C=======================================================================
  4. C= T T E 1 0 =
  5. C= --------- =
  6. C= =
  7. C= Fonction : =
  8. C= ---------- =
  9. C= Acquisition des caracteristiques d'integration aux noeuds, aux =
  10. C= points d'integration (Gauss) ou au centre de gravite pour =
  11. C= l'element thermique support TE10. =
  12. C= =
  13. C= Parametres : (E)=Entree (S)=Sortie =
  14. C= ------------ =
  15. C= ISUPTH (E) Indique le support d'integration demande =
  16. C= IPINTE (S) Pointeur sur le segment MINTE associe (ACTIF en S) =
  17. C=======================================================================
  18.  
  19. SUBROUTINE TTE10 (ISUPTH,IPINTE)
  20.  
  21. IMPLICIT INTEGER(I-N)
  22. IMPLICIT REAL*8 (A-H,O-Z)
  23.  
  24. -INC SMINTE
  25.  
  26. PARAMETER ( XZER=0.D0 , DEUX=2.D0 , UN=1.D0 , QUATRE=4.D0 )
  27. PARAMETER ( UNDEMI=.5D0 )
  28. PARAMETER ( UNQUA=.25D0 )
  29. PARAMETER ( UNSIX=1.D0/6.D0 )
  30. PARAMETER ( OHUIT=-.80D0 )
  31. PARAMETER ( OQUCI= .45D0 )
  32. PARAMETER ( UN30=1.D0/30.D0 )
  33.  
  34. NBNO = 10
  35. IF (ISUPTH.EQ.1) THEN
  36. NBPGAU = 10
  37. ELSE IF (ISUPTH.EQ.2) THEN
  38. NBPGAU = 5
  39. ELSE IF (ISUPTH.EQ.3) THEN
  40. NBPGAU = 1
  41. ELSE
  42. RETURN
  43. ENDIF
  44. SEGINI,MINTE
  45.  
  46. C COORDONNEES ET POIDS DES 10 NOEUDS
  47. IF (ISUPTH.EQ.1) THEN
  48. QSIGAU( 2) = UNDEMI
  49. QSIGAU( 3) = UN
  50. QSIGAU( 4) = UNDEMI
  51. QSIGAU( 8) = UNDEMI
  52. ETAGAU( 4) = UNDEMI
  53. ETAGAU( 5) = UN
  54. ETAGAU( 6) = UNDEMI
  55. ETAGAU( 9) = UNDEMI
  56. DZEGAU( 7) = UNDEMI
  57. DZEGAU( 8) = UNDEMI
  58. DZEGAU( 9) = UNDEMI
  59. DZEGAU(10) = UN
  60. DO I = 1, NBPGAU
  61. POIGAU(I) = UN30
  62. ENDDO
  63.  
  64. C COORDONNEES ET POIDS DES 5 POINTS DE GAUSS
  65. ELSE IF (ISUPTH.EQ.2) THEN
  66.  
  67. QSIGAU(1) = UNQUA
  68. QSIGAU(2) = UNSIX
  69. QSIGAU(3) = UNSIX
  70. QSIGAU(4) = UNSIX
  71. QSIGAU(5) = UNDEMI
  72.  
  73. ETAGAU(1) = UNQUA
  74. ETAGAU(2) = UNSIX
  75. ETAGAU(3) = UNSIX
  76. ETAGAU(4) = UNDEMI
  77. ETAGAU(5) = UNSIX
  78.  
  79. DZEGAU(1) = UNQUA
  80. DZEGAU(2) = UNSIX
  81. DZEGAU(3) = UNDEMI
  82. DZEGAU(4) = UNSIX
  83. DZEGAU(5) = UNSIX
  84.  
  85. POIGAU(1) = OHUIT*UNSIX
  86. POIGAU(2) = OQUCI*UNSIX
  87. POIGAU(3) = OQUCI*UNSIX
  88. POIGAU(4) = OQUCI*UNSIX
  89. POIGAU(5) = OQUCI*UNSIX
  90.  
  91. C COORDONNEES ET POIDS DU CENTRE DE GRAVITE
  92. ELSE IF (ISUPTH.EQ.3) THEN
  93. QSIGAU(1) = XZer
  94. ETAGAU(1) = XZer
  95. DZEGAU(1) = XZer
  96. POIGAU(1) = UN/(DEUX+UN)
  97. ENDIF
  98.  
  99. C FONCTIONS DE FORME ET LEURS DERIVEES AU(X) NBPGAU POINT(S)
  100. DO I = 1, NBPGAU
  101. XI = QSIGAU(I)
  102. YI = ETAGAU(I)
  103. ZI = DZEGAU(I)
  104.  
  105. AUX = UN-XI-YI-ZI
  106. SHPTOT(1, 1,I) = AUX*(DEUX*AUX-UN)
  107. SHPTOT(1, 2,I) = QUATRE*XI*AUX
  108. SHPTOT(1, 3,I) = XI*(DEUX*XI-UN)
  109. SHPTOT(1, 4,I) = QUATRE*XI*YI
  110. SHPTOT(1, 5,I) = YI*(DEUX*YI-UN)
  111. SHPTOT(1, 6,I) = QUATRE*YI*AUX
  112. SHPTOT(1, 7,I) = QUATRE*ZI*AUX
  113. SHPTOT(1, 8,I) = QUATRE*XI*ZI
  114. SHPTOT(1, 9,I) = QUATRE*YI*ZI
  115. SHPTOT(1,10,I) = ZI*(DEUX*ZI-UN)
  116.  
  117. SHPTOT(2, 1,I) = UN-QUATRE*AUX
  118. SHPTOT(2, 2,I) = QUATRE*(AUX-XI)
  119. SHPTOT(2, 3,I) = QUATRE*XI-UN
  120. SHPTOT(2, 4,I) = QUATRE*YI
  121. SHPTOT(2, 5,I) = XZER
  122. SHPTOT(2, 6,I) = -QUATRE*YI
  123. SHPTOT(2, 7,I) = -QUATRE*ZI
  124. SHPTOT(2, 8,I) = QUATRE*ZI
  125. SHPTOT(2, 9,I) = XZER
  126. SHPTOT(2,10,I) = XZER
  127.  
  128. SHPTOT(3, 1,I) = UN-QUATRE*AUX
  129. SHPTOT(3, 2,I) = -QUATRE*XI
  130. SHPTOT(3, 3,I) = XZER
  131. SHPTOT(3, 4,I) = QUATRE*XI
  132. SHPTOT(3, 5,I) = QUATRE*YI-UN
  133. SHPTOT(3, 6,I) = QUATRE*(AUX-YI)
  134. SHPTOT(3, 7,I) = -QUATRE*ZI
  135. SHPTOT(3, 8,I) = XZER
  136. SHPTOT(3, 9,I) = QUATRE*ZI
  137. SHPTOT(3,10,I) = XZER
  138.  
  139. SHPTOT(4, 1,I) = UN-QUATRE*AUX
  140. SHPTOT(4, 2,I) = -QUATRE*XI
  141. SHPTOT(4, 3,I) = XZER
  142. SHPTOT(4, 4,I) = XZER
  143. SHPTOT(4, 5,I) = XZER
  144. SHPTOT(4, 6,I) = -QUATRE*YI
  145. SHPTOT(4, 7,I) = QUATRE*(AUX-ZI)
  146. SHPTOT(4, 8,I) = QUATRE*XI
  147. SHPTOT(4, 9,I) = QUATRE*YI
  148. SHPTOT(4,10,I) = QUATRE*ZI-UN
  149.  
  150. ENDDO
  151.  
  152. IPINTE = MINTE
  153.  
  154. c RETURN
  155. END
  156.  
  157.  
  158.  

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