Télécharger tpr15.eso

Retour à la liste

Numérotation des lignes :

tpr15
  1. C TPR15 SOURCE OF166741 23/12/04 21:15:14 11800
  2.  
  3. C=======================================================================
  4. C= T P R 1 5 =
  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 PR15. =
  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 TPR15 (ISUPTH,IPINTE)
  20.  
  21. IMPLICIT INTEGER(I-N)
  22. IMPLICIT REAL*8 (A-H,O-Z)
  23.  
  24. -INC SMINTE
  25.  
  26. PARAMETER ( O0=0.D0, O1=1.D0, O5=.5D0, O2=2.D0, O3=3.D0, O4=4.D0 )
  27. PARAMETER ( S=0.774596669241483377035853079956479D0)
  28. PARAMETER ( A=0.470142064105115089770441209513447D0)
  29. PARAMETER ( B=0.101286507323456338800987361915124D0)
  30. PARAMETER ( C=0.333333333333333333333333333333333D0)
  31. PARAMETER ( W234=0.0661970763942530D0 )
  32. PARAMETER ( W567=0.0629695902724135D0 )
  33. PARAMETER ( W234HB=W234*5.D0/9.D0 )
  34. PARAMETER ( W234C=W234*8.D0/9.D0 )
  35. PARAMETER ( W567HB=W567*5.D0/9.D0 )
  36. PARAMETER ( W567C=W567*8.D0/9.D0 )
  37. PARAMETER ( W116=1.D0/16.D0 )
  38. PARAMETER ( W01=0.1D0 )
  39.  
  40. DIMENSION VT(6), DVXT(6),DVYT(6)
  41.  
  42. NBNO = 15
  43. IF (ISUPTH.EQ.1) THEN
  44. NBPGAU = 15
  45. ELSE IF (ISUPTH.EQ.2) THEN
  46. NBPGAU = 21
  47. ELSE IF (ISUPTH.EQ.3) THEN
  48. NBPGAU = 1
  49. ELSE
  50. RETURN
  51. ENDIF
  52. SEGINI,MINTE
  53.  
  54. C COORDONNEES INTRINSEQUES DES 15 NOEUDS
  55. IF (ISUPTH.EQ.1) THEN
  56. QSIGAU( 1) = O0
  57. QSIGAU( 2) = O5
  58. QSIGAU( 3) = O1
  59. QSIGAU( 4) = O5
  60. QSIGAU( 5) = O0
  61. QSIGAU( 6) = O0
  62. QSIGAU( 7) = O0
  63. QSIGAU( 8) = O1
  64. QSIGAU( 9) = O0
  65. QSIGAU(10) = O0
  66. QSIGAU(11) = O5
  67. QSIGAU(12) = O1
  68. QSIGAU(13) = O5
  69. QSIGAU(14) = O0
  70. QSIGAU(15) = O0
  71.  
  72. ETAGAU( 1) = O0
  73. ETAGAU( 2) = O0
  74. ETAGAU( 3) = O0
  75. ETAGAU( 4) = O5
  76. ETAGAU( 5) = O1
  77. ETAGAU( 6) = O5
  78. ETAGAU( 7) = O0
  79. ETAGAU( 8) = O0
  80. ETAGAU( 9) = O1
  81. ETAGAU(10) = O0
  82. ETAGAU(11) = O0
  83. ETAGAU(12) = O0
  84. ETAGAU(13) = O5
  85. ETAGAU(14) = O1
  86. ETAGAU(15) = O5
  87.  
  88. DZEGAU( 1) = -O1
  89. DZEGAU( 2) = -O1
  90. DZEGAU( 3) = -O1
  91. DZEGAU( 4) = -O1
  92. DZEGAU( 5) = -O1
  93. DZEGAU( 6) = -O1
  94. DZEGAU( 7) = O0
  95. DZEGAU( 8) = O0
  96. DZEGAU( 9) = O0
  97. DZEGAU(10) = O1
  98. DZEGAU(11) = O1
  99. DZEGAU(12) = O1
  100. DZEGAU(13) = O1
  101. DZEGAU(14) = O1
  102. DZEGAU(15) = O1
  103.  
  104. r_z = O1 / 15.D0
  105. DO I = 1, NBPGAU
  106. POIGAU(1) = r_z
  107. ENDDO
  108.  
  109. C COORDONNEES ET POIDS DES 21 POINTS DE GAUSS
  110. ELSE IF (ISUPTH.EQ.2) THEN
  111. QSIGAU( 1) = C
  112. QSIGAU( 2) = A
  113. QSIGAU( 3) = O1-O2*A
  114. QSIGAU( 4) = A
  115. QSIGAU( 5) = B
  116. QSIGAU( 6) = O1-O2*B
  117. QSIGAU( 7) = B
  118. QSIGAU( 8) = C
  119. QSIGAU( 9) = A
  120. QSIGAU(10) = O1-O2*A
  121. QSIGAU(11) = A
  122. QSIGAU(12) = B
  123. QSIGAU(13) = O1-O2*B
  124. QSIGAU(14) = B
  125. QSIGAU(15) = C
  126. QSIGAU(16) = A
  127. QSIGAU(17) = O1-O2*A
  128. QSIGAU(18) = A
  129. QSIGAU(19) = B
  130. QSIGAU(20) = O1-O2*B
  131. QSIGAU(21) = B
  132.  
  133. ETAGAU( 1) = C
  134. ETAGAU( 2) = A
  135. ETAGAU( 3) = A
  136. ETAGAU( 4) = O1-O2*A
  137. ETAGAU( 5) = B
  138. ETAGAU( 6) = B
  139. ETAGAU( 7) = O1-O2*B
  140. ETAGAU( 8) = C
  141. ETAGAU( 9) = A
  142. ETAGAU(10) = A
  143. ETAGAU(11) = O1-O2*A
  144. ETAGAU(12) = B
  145. ETAGAU(13) = B
  146. ETAGAU(14) = O1-O2*B
  147. ETAGAU(15) = C
  148. ETAGAU(16) = A
  149. ETAGAU(17) = A
  150. ETAGAU(18) = O1-O2*A
  151. ETAGAU(19) = B
  152. ETAGAU(20) = B
  153. ETAGAU(21) = O1-O2*B
  154.  
  155. DZEGAU( 1) = -S
  156. DZEGAU( 2) = -S
  157. DZEGAU( 3) = -S
  158. DZEGAU( 4) = -S
  159. DZEGAU( 5) = -S
  160. DZEGAU( 6) = -S
  161. DZEGAU( 7) = -S
  162. DZEGAU( 8) = O0
  163. DZEGAU( 9) = O0
  164. DZEGAU(10) = O0
  165. DZEGAU(11) = O0
  166. DZEGAU(12) = O0
  167. DZEGAU(13) = O0
  168. DZEGAU(14) = O0
  169. DZEGAU(15) = +S
  170. DZEGAU(16) = +S
  171. DZEGAU(17) = +S
  172. DZEGAU(18) = +S
  173. DZEGAU(19) = +S
  174. DZEGAU(20) = +S
  175. DZEGAU(21) = +S
  176.  
  177. POIGAU( 1) = W116
  178. POIGAU( 2) = W234HB
  179. POIGAU( 3) = W234HB
  180. POIGAU( 4) = W234HB
  181. POIGAU( 5) = W567HB
  182. POIGAU( 6) = W567HB
  183. POIGAU( 7) = W567HB
  184. POIGAU( 8) = W01
  185. POIGAU( 9) = W234C
  186. POIGAU(10) = W234C
  187. POIGAU(11) = W234C
  188. POIGAU(12) = W567C
  189. POIGAU(13) = W567C
  190. POIGAU(14) = W567C
  191. POIGAU(15) = W116
  192. POIGAU(16) = W234HB
  193. POIGAU(17) = W234HB
  194. POIGAU(18) = W234HB
  195. POIGAU(19) = W567HB
  196. POIGAU(20) = W567HB
  197. POIGAU(21) = W567HB
  198.  
  199. C COORDONNEES DU CENTRE DE GRAVITE
  200. ELSE IF (ISUPTH.EQ.3) THEN
  201. QSIGAU(1) = O0
  202. ETAGAU(1) = O0
  203. DZEGAU(1) = O0
  204. POIGAU(1) = O1
  205. ENDIF
  206.  
  207. C FONCTIONS DE FORME ET LEURS DERIVEES AU(X) NBPGAU POINT(S)
  208. DO I = 1, NBPGAU
  209. X = QSIGAU(I)
  210. Y = ETAGAU(I)
  211. Z = DZEGAU(I)
  212. C- FONCTIONS DE FORME DU T6 DE BASE
  213. XY1 = O1-X-Y
  214. VT(1) = O2*XY1*(O5-X-Y)
  215. VT(2) = O4*X*XY1
  216. VT(3) = O2*X*(X-O5)
  217. VT(4) = O4*X*Y
  218. VT(5) = O2*Y*(Y-O5)
  219. VT(6) = O4*XY1*Y
  220. C- FONCTIONS DE FORME DU PRISME
  221. ZM1 = (O1-Z)*O5
  222. ZP1 = (O1+Z)*O5
  223. ZM2 = O1-Z*Z
  224.  
  225. SHPTOT(1, 7,I) = ZM2*XY1
  226. SHPTOT(1, 8,I) = ZM2*X
  227. SHPTOT(1, 9,I) = ZM2*Y
  228.  
  229. SHPTOT(1, 1,I) = ZM1*VT(1)-SHPTOT(1,7,I)*O5
  230. SHPTOT(1, 2,I) = ZM1*VT(2)
  231. SHPTOT(1, 3,I) = ZM1*VT(3)-SHPTOT(1,8,I)*O5
  232. SHPTOT(1, 4,I) = ZM1*VT(4)
  233. SHPTOT(1, 5,I) = ZM1*VT(5)-SHPTOT(1,9,I)*O5
  234. SHPTOT(1, 6,I) = ZM1*VT(6)
  235.  
  236. SHPTOT(1,10,I) = ZP1*VT(1)-SHPTOT(1,7,I)*O5
  237. SHPTOT(1,11,I) = ZP1*VT(2)
  238. SHPTOT(1,12,I) = ZP1*VT(3)-SHPTOT(1,8,I)*O5
  239. SHPTOT(1,13,I) = ZP1*VT(4)
  240. SHPTOT(1,14,I) = ZP1*VT(5)-SHPTOT(1,9,I)*O5
  241. SHPTOT(1,15,I) = ZP1*VT(6)
  242.  
  243. C- DERIVEES DES FONCTIONS DE FORME DU T6 DE BASE
  244. DVXT(1) = O4*(X+Y)-O3
  245. DVXT(2) = O4*(O1-X-X-Y)
  246. DVXT(3) = O4*X-O1
  247. DVXT(4) = O4*Y
  248. DVXT(5) = O0
  249. DVXT(6) = -O4*Y
  250.  
  251. DVYT(1) = O4*(X+Y)-O3
  252. DVYT(2) = -O4*X
  253. DVYT(3) = O0
  254. DVYT(4) = O4*X
  255. DVYT(5) = O4*Y-O1
  256. DVYT(6) = O4*(O1-X-Y-Y)
  257.  
  258. C- DERIVEES EN X DES FONCTIONS DE FORME
  259. SHPTOT(2, 7,I) = -ZM2
  260. SHPTOT(2, 8,I) = ZM2
  261. SHPTOT(2, 9,I) = O0
  262.  
  263. SHPTOT(2, 1,I) = ZM1*DVXT(1)-SHPTOT(2,7,I)*O5
  264. SHPTOT(2, 2,I) = ZM1*DVXT(2)
  265. SHPTOT(2, 3,I) = ZM1*DVXT(3)-SHPTOT(2,8,I)*O5
  266. SHPTOT(2, 4,I) = ZM1*DVXT(4)
  267. SHPTOT(2, 5,I) = ZM1*DVXT(5)-SHPTOT(2,9,I)*O5
  268. SHPTOT(2, 6,I) = ZM1*DVXT(6)
  269.  
  270. SHPTOT(2,10,I) = ZP1*DVXT(1)-SHPTOT(2,7,I)*O5
  271. SHPTOT(2,11,I) = ZP1*DVXT(2)
  272. SHPTOT(2,12,I) = ZP1*DVXT(3)-SHPTOT(2,8,I)*O5
  273. SHPTOT(2,13,I) = ZP1*DVXT(4)
  274. SHPTOT(2,14,I) = ZP1*DVXT(5)-SHPTOT(2,9,I)*O5
  275. SHPTOT(2,15,I) = ZP1*DVXT(6)
  276.  
  277. C- DERIVEES EN Y DES FONCTIONS DE FORME
  278. SHPTOT(3, 7,I) = -ZM2
  279. SHPTOT(3, 8,I) = O0
  280. SHPTOT(3, 9,I) = ZM2
  281.  
  282. SHPTOT(3, 1,I) = ZM1*DVYT(1)-SHPTOT(3,7,I)*O5
  283. SHPTOT(3, 2,I) = ZM1*DVYT(2)
  284. SHPTOT(3, 3,I) = ZM1*DVYT(3)-SHPTOT(3,8,I)*O5
  285. SHPTOT(3, 4,I) = ZM1*DVYT(4)
  286. SHPTOT(3, 5,I) = ZM1*DVYT(5)-SHPTOT(3,9,I)*O5
  287. SHPTOT(3, 6,I) = ZM1*DVYT(6)
  288.  
  289. SHPTOT(3,10,I) = ZP1*DVYT(1)-SHPTOT(3,7,I)*O5
  290. SHPTOT(3,11,I) = ZP1*DVYT(2)
  291. SHPTOT(3,12,I) = ZP1*DVYT(3)-SHPTOT(3,8,I)*O5
  292. SHPTOT(3,13,I) = ZP1*DVYT(4)
  293. SHPTOT(3,14,I) = ZP1*DVYT(5)-SHPTOT(3,9,I)*O5
  294. SHPTOT(3,15,I) = ZP1*DVYT(6)
  295.  
  296. C- DERIVEES EN Z DES FONCTIONS DE FORME
  297. SHPTOT(4, 7,I) = -O2*Z*XY1
  298. SHPTOT(4, 8,I) = -O2*Z*X
  299. SHPTOT(4, 9,I) = -O2*Z*Y
  300.  
  301. SHPTOT(4, 1,I) = -(VT(1)+SHPTOT(4,7,I))*O5
  302. SHPTOT(4, 2,I) = -VT(2)*O5
  303. SHPTOT(4, 3,I) = -(VT(3)+SHPTOT(4,8,I))*O5
  304. SHPTOT(4, 4,I) = -VT(4)*O5
  305. SHPTOT(4, 5,I) = -(VT(5)+SHPTOT(4,9,I))*O5
  306. SHPTOT(4, 6,I) = -VT(6)*O5
  307.  
  308. SHPTOT(4,10,I) = (VT(1)-SHPTOT(4,7,I))*O5
  309. SHPTOT(4,11,I) = VT(2)*O5
  310. SHPTOT(4,12,I) = (VT(3)-SHPTOT(4,8,I))*O5
  311. SHPTOT(4,13,I) = VT(4)*O5
  312. SHPTOT(4,14,I) = (VT(5)-SHPTOT(4,9,I))*O5
  313. SHPTOT(4,15,I) = VT(6)*O5
  314.  
  315. ENDDO
  316.  
  317. IPINTE = MINTE
  318.  
  319. C RETURN
  320. END
  321.  
  322.  
  323.  

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