Télécharger donred.eso

Retour à la liste

Numérotation des lignes :

  1. C DONRED SOURCE BP208322 17/03/20 21:15:01 9359
  2. SUBROUTINE DONRED(NBG,IELE,MELE,NNO,NPINT,IPT,IRET)
  3. C=======================================================================
  4. C
  5. C CALCULE LES COORDONNEES REDUITES ET LES POIDS D'INTEGRATION
  6. C DES POINTS DE GAUSS D'UN ELEMENT
  7. C
  8. C ENTREES
  9. C NBG = NOMBRE DE POINT DE GAUSS
  10. C IELE = NUMEROTATION DES ELEMENTS DANS NOMS(32)
  11. C (VOIR CCGEOME)
  12. C MELE = NUMEROTATION DANS NOMTP
  13. C NNO = NOMBRE DE NOEUDS DE L'ELEMENT NOMS(IELE)
  14. C NPINT = NOMBRE DE POINTS D'INTEGRATION DONS LE CAS DES
  15. C ELEMENTS COQUES INTEGRES
  16. C SORTIES
  17. C IPT = POINTEUR SUR MINTE
  18. C IRET = 1 OU 0 SUIVANT QUE IELE ET NBG COMPATIBLES OU PAS
  19. C
  20. C PAS DE MESSAGES D'ERREUR ECRIT DANS CETTE ROUTINE
  21. C
  22. C QSI(NBG) ETA(NBG) DZE(NBG) = COORDONNEES REDUITES
  23. C
  24. C POI(NBG) = POIDS DE L'INTEGRATION
  25. C C A D VOLUME DE L ELEMENT DE BASE
  26. C FOIS CONTRIBUTION DU POINT DE GAUSS
  27. C
  28. C ON NE FAIT LA DISTINCTION QU'ENTRE SEGMENTS
  29. C TRIANGLES
  30. C CARRES
  31. C TETRAEDRES
  32. C PRISMES
  33. C PYRAMIDES
  34. C CUBES
  35. C POLYGONES
  36. C
  37. *C PPJ CHANGEMENT DES REGLES D'INTEGRATION POUR LES JOINTS
  38. C=======================================================================
  39. IMPLICIT INTEGER(I-N)
  40. IMPLICIT REAL*8(A-H,O-Z)
  41. -INC CCREEL
  42. -INC SMINTE
  43. C
  44. C QUELQUES CONSTANTES EN VRAC
  45. C
  46. DATA HUIT/8.D0/
  47. DATA DEUX,TROIS,QUATRE/2.D0,3.D0,4.D0/
  48. DATA UNQUA,UNDEMI,UN/.25D0,.5D0,1.D0/
  49. DATA UNTIER/.33333333333333333D0/
  50. DATA DEUTIE/.66666666666666666D0/
  51. DATA UNCINQ/.2D0/
  52. DATA UNSIX/.166666666666666666D0/
  53. DATA UN24/.4166666666666666666D-1/
  54. C
  55. C INTEGRATION 2 POINTS LIGNE (ordre 3)
  56. C
  57. DATA X577/.577350269189626D0/
  58. C
  59. C INTEGRATION 3 POINTS LIGNE (ordre 5)
  60. C
  61. DATA P555/.555555555555555555D0/
  62. DATA P888/.888888888888888888D0/
  63. DATA X774/.774596669241483D0/
  64. C
  65. C INTEGRATION 4 POINTS LIGNE (ordre 7)
  66. C
  67. DATA X861/.861136311594053D0/
  68. DATA X339/.339981043584856D0/
  69. DATA P347/.347854845137454D0/
  70. DATA P652/.652145154862546D0/
  71. C
  72. C INTEGRATION 4 POINTS TRIANGLE (ordre 3)
  73. C attention : poids negatif !
  74. C P260=25./96. P281=-27/96.
  75. DATA P260/.2604166666666666D0/
  76. DATA P281/-.28125D0/
  77. C
  78. C INTEGRATION 6 POINTS TRIANGLE (ordre 4?)
  79. C
  80. DATA X044/.44594849091596D0/
  81. DATA X009/.09157621350977D0/
  82. DATA X111/.11169079483900D0/
  83. DATA X005/.05497587182766D0/
  84. C
  85. C INTEGRATION 7 POINTS TRIANGLE (ordre 5)
  86. C
  87. DATA X101/.1012865073235D0/
  88. DATA X797/.7974269853531D0/
  89. DATA X470/.4701420641051D0/
  90. DATA X059/.0597158717898D0/
  91. DATA P125/.1259391805448D0/
  92. DATA P132/.1323941527885D0/
  93. C
  94. C INTEGRATION 4 POINTS TETRAEDRE (ordre 2)
  95. C
  96. DATA X138/.1381966011250105D0/
  97. DATA X585/.5854101966249684D0/
  98. C
  99. C INTEGRATION 11 POINTS TETRAEDRE (ordre 4)
  100. C tiré de [DHATT, FEM, p371] (attention : poids negatif !)
  101. C XONZ14=11/14 XUN14=1/14. P01315=-74/5625 P00762=343/45000
  102. DATA XONZ14/0.785714285714285698D0/
  103. DATA XUN14 /0.071428571428571426D0/
  104. DATA P01315/-.013155555555555555D0/
  105. DATA P00762/0.007622222222222222D0/
  106. DATA P00124/0.001244444444444444D0/
  107. C
  108. C INTEGRATION 15 POINTS TETRAEDRE (ordre 5)
  109. C tiré de [DHATT, FEM, p371]
  110. c X3197 = (7.+sqrt(5.))/34. X0919 = (7.-sqrt(5.))/34.
  111. c X7240 = (13.+3.*sqrt(5.)))/34. X0406 = (13.-3.*sqrt(5.)))/34.
  112. c X0563 = (5.-sqrt(15.))/20. X4436 = (5.-sqrt(15.))/20.
  113. c P0197 = 112./5670. P0088 = 5./567.
  114. c P0119=(2665.+14.sqrt(15.))/226800. P0115=(2665.-14.sqrt(15.))/226800.
  115. DATA X3197/0.319793627829629890D0/
  116. DATA X0919/0.091971078052723032D0/
  117. DATA X7240/0.724086765841830959D0/
  118. DATA X0406/0.040619116511110234D0/
  119. DATA X0563/0.056350832689629149D0/
  120. DATA X4436/0.443649167310370796D0/
  121. DATA P0197/0.019753086419753086D0/
  122. DATA P0119/0.011989513963169770D0/
  123. DATA P0115/0.011511367871045397D0/
  124. DATA P0088/0.008818342151675485D0/
  125. C
  126. NBPGAU=NBG
  127. NBNO=NNO
  128. SEGINI MINTE
  129. IPT=MINTE
  130. IRET=1
  131. DO 100 I=1,NBG
  132. QSIGAU(I)=XZERO
  133. ETAGAU(I)=XZERO
  134. DZEGAU(I)=XZERO
  135. POIGAU(I)=XZERO
  136. 100 CONTINUE
  137. *
  138. * PETIT DEBRANCHEMENT POUR ELEMENT TUYO
  139. IF(MELE.EQ.96) GO TO 96
  140. * ET COQUES INTEGREES (sauf pour la matrice masse)
  141. IF(NPINT.NE.0.AND.NBG.NE.3)GOTO 95
  142. * ET ELEMENT TR6H
  143. IF(MELE.EQ.157)GOTO 1570
  144. C
  145. C 66 INDIQUE QUE L ELEMENT N EST PAS ENCORE IMPLEMENTE
  146. C DANS CE CAS IRET EST MIS A XZERO
  147. C
  148. C P1 S2 S3 T3 T4 T6 T7 Q4 Q5 Q8 Q9 R2 R3 C8 C20 P6 P15
  149. GOTO ( 1, 2, 2, 4, 4, 4, 4, 8, 8, 8, 8, 2,12,14, 14,16, 16,
  150. C L3 L4 L6 L8 MU T4 T10 P5 P13 ATT SUP RP3 LP6 LP8
  151. 1 4, 8,66,66,66,23, 23,25, 25, 66, 66, 2, 4, 8,
  152. C POLY
  153. 2 30),IELE
  154. C
  155. 66 IRET=0
  156. IPT=0
  157. GOTO 666
  158.  
  159. C=======================================================================
  160. C ELEMENT POI1
  161. C=======================================================================
  162. C
  163. 1 CONTINUE
  164. POIGAU(1)=UN
  165. GOTO 666
  166. C
  167. C=======================================================================
  168. C ELEMENT LIGNE LONGUEUR =2
  169. C=======================================================================
  170. C
  171. 2 IF (NBG.NE.1) GOTO 220
  172. C
  173. C UN POINT D INTEGRATION (ordre 1)
  174. C
  175. POIGAU(1)=DEUX
  176. GOTO 666
  177. 220 IF (NBG.NE.2) GOTO 230
  178. C
  179. C 2 POINTS D INTEGRATION (ordre 3)
  180. C
  181. POIGAU(1)= UN
  182. POIGAU(2)= UN
  183. QSIGAU(1)=-X577
  184. QSIGAU(2)= X577
  185. * les contraintes et les deformations des él. poutre et tuyau sont
  186. * calculées analytiquement aux noeuds. On change les coordonées des
  187. * points de Gauss afin d'avoir les bonnes valeurs quand on fait "CHANGER"
  188. *
  189. * ON INTEGRE AU NOEUD (LOBATTO) POUR L'ELEMENT DE JOINT 85 ET 170 (JOI2 ET JGI2)
  190. * (--> ordre 2)
  191. * IF(MELE.EQ.29.OR.MELE.EQ.42.OR.MELE.EQ.85) THEN
  192. IF(MELE.EQ.29.OR.MELE.EQ.42) THEN
  193. QSIGAU(1)=-1.D0
  194. QSIGAU(2)= 1.D0
  195. ENDIF
  196. GOTO 666
  197. C
  198. C 3 POINTS D INTEGRATION (ordre 5)
  199. C
  200. 230 IF (NBG.NE.3) GOTO 240
  201. POIGAU(1)= P555
  202. POIGAU(2)= P888
  203. POIGAU(3)= P555
  204. QSIGAU(1)=-X774
  205. QSIGAU(3)= X774
  206. GOTO 666
  207. C
  208. C 4 POINTS D INTEGRATION (ordre 7)
  209. C
  210. 240 IF(NBG.NE.4) GOTO 66
  211. QSIGAU(1)=-X861
  212. QSIGAU(2)=-X339
  213. QSIGAU(3)= X339
  214. QSIGAU(4)= X861
  215. POIGAU(1)= P347
  216. POIGAU(2)= P652
  217. POIGAU(3)= P652
  218. POIGAU(4)= P347
  219. GOTO 666
  220. C
  221. C=======================================================================
  222. C ELEMENT TRIANGLE SURFACE=1/2
  223. C=======================================================================
  224. C
  225. 4 CONTINUE
  226. C
  227. C 1 POINTS D INTEGRATION (ordre 1)
  228. C
  229. IF (NBG.NE.1) GOTO 430
  230. QSIGAU(1)=UNTIER
  231. ETAGAU(1)=UNTIER
  232. POIGAU(1)=UNDEMI
  233. GOTO 666
  234. C
  235. C 3 POINTS D INTEGRATION (ordre 3)
  236. C [Zienkiewicz, p.165] ou [Bathe, p.467]
  237. C
  238. 430 IF(NBG.NE.3) GOTO 440
  239. IF(MELE.EQ.56.OR.MELE.EQ.128) THEN
  240. QSIGAU(1)=UNSIX
  241. ETAGAU(1)=UNSIX
  242. QSIGAU(2)=UNSIX
  243. ETAGAU(2)=DEUTIE
  244. QSIGAU(3)=DEUTIE
  245. ETAGAU(3)=UNSIX
  246. C
  247. C CAS DU JOT3, DU JCT3 ET DU JGT3 ON SE PLACE AUX SOMMETS
  248. C
  249. ELSEIF(MELE.EQ.87)THEN
  250. QSIGAU(2)=UN
  251. ETAGAU(3)=UN
  252. C
  253. ELSE
  254. QSIGAU(1)=UNDEMI
  255. ETAGAU(2)=UNDEMI
  256. QSIGAU(3)=UNDEMI
  257. ETAGAU(3)=UNDEMI
  258. ENDIF
  259. POIGAU(1)=UNSIX
  260. POIGAU(2)=UNSIX
  261. POIGAU(3)=UNSIX
  262. GOTO 666
  263. C
  264. C 4 POINTS D INTEGRATION (ordre 3)
  265. C [Zienkiewicz, p.165] (attention : poids negatif !)
  266. C
  267. 440 IF (NBG.NE.4) GOTO 460
  268. QSIGAU(1)=UNTIER
  269. QSIGAU(2)=UNCINQ
  270. QSIGAU(3)=UNCINQ*TROIS
  271. QSIGAU(4)=UNCINQ
  272. ETAGAU(1)=UNTIER
  273. ETAGAU(2)=UNCINQ
  274. ETAGAU(3)=UNCINQ
  275. ETAGAU(4)=UNCINQ*TROIS
  276. POIGAU(1)=P281
  277. POIGAU(2)=P260
  278. POIGAU(3)=P260
  279. POIGAU(4)=P260
  280. GOTO 666
  281. C
  282. C 6 POINTS D INTEGRATION COQUE EPAISSE
  283. C
  284. 460 IF(NBG.NE.6) GOTO 470
  285. QSIGAU(1)=UNSIX
  286. ETAGAU(1)=UNSIX
  287. QSIGAU(2)=UNSIX
  288. ETAGAU(2)=DEUTIE
  289. QSIGAU(3)=DEUTIE
  290. ETAGAU(3)=UNSIX
  291. DZEGAU(1)=-X577
  292. DZEGAU(2)=-X577
  293. DZEGAU(3)=-X577
  294. QSIGAU(4)=UNSIX
  295. ETAGAU(4)=UNSIX
  296. QSIGAU(5)=UNSIX
  297. ETAGAU(5)=DEUTIE
  298. QSIGAU(6)=DEUTIE
  299. ETAGAU(6)=UNSIX
  300. DZEGAU(4)= X577
  301. DZEGAU(5)= X577
  302. DZEGAU(6)= X577
  303. POIGAU(1)=UNSIX
  304. POIGAU(2)=UNSIX
  305. POIGAU(3)=UNSIX
  306. POIGAU(4)=UNSIX
  307. POIGAU(5)=UNSIX
  308. POIGAU(6)=UNSIX
  309. GOTO 666
  310. C
  311. C 7 POINTS D INTEGRATION (ordre 5)
  312. C [Zienkiewicz, p.165] ou BATHE P280
  313. 470 IF(NBG.NE.7) GOTO 777
  314. QSIGAU(1)=X101
  315. QSIGAU(2)=X797
  316. QSIGAU(3)=X101
  317. QSIGAU(4)=X470
  318. QSIGAU(5)=X470
  319. QSIGAU(6)=X059
  320. QSIGAU(7)=UNTIER
  321. ETAGAU(1)=QSIGAU(1)
  322. ETAGAU(2)=QSIGAU(1)
  323. ETAGAU(3)=QSIGAU(2)
  324. ETAGAU(4)=QSIGAU(6)
  325. ETAGAU(5)=QSIGAU(4)
  326. ETAGAU(6)=QSIGAU(4)
  327. ETAGAU(7)=QSIGAU(7)
  328. POIGAU(1)=P125*UNDEMI
  329. POIGAU(2)=P125*UNDEMI
  330. POIGAU(3)=P125*UNDEMI
  331. POIGAU(4)=P132*UNDEMI
  332. POIGAU(5)=P132*UNDEMI
  333. POIGAU(6)=P132*UNDEMI
  334. POIGAU(7)=.1125D0
  335. GOTO 666
  336. C
  337. C PLUSIEURS NAPPES DE 3 POINTS D INTEGRATION
  338. C
  339. 777 IF(NBG.NE.9.AND.NBG.NE.12.AND.NBG.NE.15.AND.NBG.NE.18.AND.
  340. 1 NBG.NE.21) GOTO 66
  341. NNAPPE=NBG/3
  342. POIDS=UN/DBLE(NBG)
  343. DO 7771 IA=1,NNAPPE
  344. NP=IA*3-2
  345. POIGAU(NP )=POIDS
  346. POIGAU(NP+1)=POIDS
  347. POIGAU(NP+2)=POIDS
  348. QSIGAU(NP )= UNSIX
  349. ETAGAU(NP )= UNSIX
  350. QSIGAU(NP+1)= UNSIX
  351. ETAGAU(NP+1)= DEUTIE
  352. QSIGAU(NP+2)= DEUTIE
  353. ETAGAU(NP+2)= UNSIX
  354. XXXX = ((DEUX*DBLE(IA)-UN)/DBLE(NNAPPE))-UN
  355. DZEGAU(NP )= XXXX
  356. DZEGAU(NP+1)= XXXX
  357. DZEGAU(NP+2)= XXXX
  358. 7771 CONTINUE
  359. GOTO 666
  360. C
  361. C=======================================================================
  362. C ELEMENT CARRE SURFACE=4
  363. C=======================================================================
  364. C
  365. 8 IF(NBG.NE.1) GOTO 84
  366. C
  367. C UN POINT D INTEGRATION (ordre 1)
  368. C
  369. POIGAU(1)=QUATRE
  370. GOTO 666
  371. 84 IF(NBG.NE.4) GOTO 85
  372. C
  373. C 4 POINTS D INTEGRATION (ordre 3)
  374. C
  375. QSIGAU(1)=-X577
  376. QSIGAU(2)= X577
  377. QSIGAU(3)= X577
  378. QSIGAU(4)=-X577
  379. ETAGAU(1)=-X577
  380. ETAGAU(2)=-X577
  381. ETAGAU(3)= X577
  382. ETAGAU(4)= X577
  383. C
  384. C CAS DU JOI4, JCI4 ET JGI4 ON SE PLACE AUX SOMMETS
  385. C
  386. IF(MELE.EQ.88)THEN
  387. QSIGAU(1)=-UN
  388. QSIGAU(2)= UN
  389. QSIGAU(3)= UN
  390. QSIGAU(4)=-UN
  391. ETAGAU(1)=-UN
  392. ETAGAU(2)=-UN
  393. ETAGAU(3)= UN
  394. ETAGAU(4)= UN
  395. ENDIF
  396. C
  397. POIGAU(1)=UN
  398. POIGAU(2)=UN
  399. POIGAU(3)=UN
  400. POIGAU(4)=UN
  401. GOTO 666
  402. 85 IF(NBG.NE.5) GOTO 89
  403. C
  404. C 5 POINTS D INTEGRATION
  405. C
  406. QSIGAU(1)=-X577
  407. QSIGAU(2)= X577
  408. QSIGAU(3)= X577
  409. QSIGAU(4)=-X577
  410. ETAGAU(1)=-X577
  411. ETAGAU(2)=-X577
  412. ETAGAU(3)= X577
  413. ETAGAU(4)= X577
  414. POIGAU(1)=UN
  415. POIGAU(2)=UN
  416. POIGAU(3)=UN
  417. POIGAU(4)=UN
  418. POIGAU(5)=QUATRE
  419. GOTO 666
  420. C
  421. C 9 POINTS D INTEGRATION (ordre 5)
  422. C
  423. 89 IF (NBG.NE.9) GOTO 816
  424. QSIGAU(1)=-X774
  425. QSIGAU(3)= X774
  426. QSIGAU(4)=-X774
  427. QSIGAU(6)= X774
  428. QSIGAU(7)=-X774
  429. QSIGAU(9)= X774
  430. ETAGAU(1)=-X774
  431. ETAGAU(2)=-X774
  432. ETAGAU(3)=-X774
  433. ETAGAU(7)= X774
  434. ETAGAU(8)= X774
  435. ETAGAU(9)= X774
  436. POIGAU(1)=P555*P555
  437. POIGAU(3)=POIGAU(1)
  438. POIGAU(7)=POIGAU(1)
  439. POIGAU(9)=POIGAU(1)
  440. POIGAU(2)=P888*P555
  441. POIGAU(4)=POIGAU(2)
  442. POIGAU(6)=POIGAU(2)
  443. POIGAU(8)=POIGAU(2)
  444. POIGAU(5)=P888*P888
  445. GOTO 666
  446. 816 IF(NBG.NE.16) GOTO 88
  447. C
  448. C CARRE 16 POINTS D INTEGRATION (ordre 7)
  449. C
  450. DO 8161 IA=1,13,4
  451. QSIGAU(IA )=-X861
  452. QSIGAU(IA+1 )=-X339
  453. QSIGAU(IA+2 )= X339
  454. QSIGAU(IA+3 )= X861
  455. POIGAU(IA )= P347
  456. POIGAU(IA+1 )= P652
  457. POIGAU(IA+2 )= P652
  458. POIGAU(IA+3 )= P347
  459. 8161 CONTINUE
  460. DO 8162 IA=1,4
  461. ETAGAU(IA )=-X861
  462. ETAGAU(IA+4 )=-X339
  463. ETAGAU(IA+8 )= X339
  464. ETAGAU(IA+12)= X861
  465. POIGAU(IA )= P347*POIGAU(IA )
  466. POIGAU(IA+4 )= P652*POIGAU(IA+4 )
  467. POIGAU(IA+8 )= P652*POIGAU(IA+8 )
  468. POIGAU(IA+12)= P347*POIGAU(IA+12)
  469. 8162 CONTINUE
  470. GOTO 666
  471. C
  472. C ELEMENT CARRE 8 POINTS D INTEGRATION POUR COQ8
  473. C
  474. 88 IF(NBG.NE.8) GOTO 888
  475. QSIGAU(1)=-X577
  476. QSIGAU(2)= X577
  477. QSIGAU(3)= X577
  478. QSIGAU(4)=-X577
  479. C
  480. ETAGAU(1)=-X577
  481. ETAGAU(2)=-X577
  482. ETAGAU(3)= X577
  483. ETAGAU(4)= X577
  484. C
  485. QSIGAU(5)=-X577
  486. QSIGAU(6)= X577
  487. QSIGAU(7)= X577
  488. QSIGAU(8)=-X577
  489. C
  490. ETAGAU(5)=-X577
  491. ETAGAU(6)=-X577
  492. ETAGAU(7)= X577
  493. ETAGAU(8)= X577
  494. C
  495. DO 881 IA=1,4
  496. POIGAU(IA )= UN
  497. POIGAU(IA+4)= UN
  498. DZEGAU(IA )=-X577
  499. DZEGAU(IA+4)= X577
  500. 881 CONTINUE
  501. GOTO 666
  502. C
  503. C PLUSIEURS NAPPES DE POINTS DE 4 POINTS D INTEGRATION
  504. C
  505. 888 IF(NBG.NE.12.AND.NBG.NE.16.AND.NBG.NE.20.AND.NBG.NE.24.AND.
  506. 1 NBG.NE.28) GOTO 66
  507. NNAPPE=NBG/4
  508. POIDS=HUIT/DBLE(NBG)
  509. DO 8881 IA=1,NNAPPE
  510. NP=IA*4-3
  511. POIGAU(NP )=POIDS
  512. POIGAU(NP+1)=POIDS
  513. POIGAU(NP+2)=POIDS
  514. POIGAU(NP+3)=POIDS
  515. QSIGAU(NP )=-X577
  516. QSIGAU(NP+1)= X577
  517. QSIGAU(NP+2)= X577
  518. QSIGAU(NP+3)=-X577
  519. ETAGAU(NP )=-X577
  520. ETAGAU(NP+1)=-X577
  521. ETAGAU(NP+2)= X577
  522. ETAGAU(NP+3)= X577
  523. XXXX = ((DEUX*DBLE(IA)-UN)/DBLE(NNAPPE))-UN
  524. DZEGAU(NP )= XXXX
  525. DZEGAU(NP+1)= XXXX
  526. DZEGAU(NP+2)= XXXX
  527. DZEGAU(NP+3)= XXXX
  528. 8881 CONTINUE
  529. GOTO 666
  530. C
  531. C=======================================================================
  532. C ELEMENT RACCORD 2
  533. C=======================================================================
  534. C
  535. 12 CONTINUE
  536. C
  537. C UN POINT D INTEGRATION
  538. C
  539. IF(NBG.NE.1) GOTO 122
  540. POIGAU(1)=DEUX
  541. GOTO 666
  542. 122 IF (NBG.NE.2) GOTO 123
  543. C
  544. C 2 POINTS D INTEGRATION
  545. C
  546. C QUE LES ELEMENTS RACCORDS
  547. C
  548. POIGAU(1)= UN
  549. POIGAU(2)= UN
  550. QSIGAU(1)=-X577
  551. QSIGAU(2)= X577
  552. ETAGAU(1)=-UN
  553. ETAGAU(2)=-UN
  554. GOTO 666
  555. C
  556. C 3 POINTS D INTEGRATION
  557. C
  558. C QUE LE LINESPRING
  559. C
  560. 123 IF(NBG.NE.3) GOTO 124
  561. QSIGAU(1)=-X774
  562. QSIGAU(3)= X774
  563. ETAGAU(1)=-UN
  564. ETAGAU(2)=-UN
  565. ETAGAU(3)=-UN
  566. POIGAU(1)= P555
  567. POIGAU(2)= P888
  568. POIGAU(3)= P555
  569. *
  570. * on integre au noeud (Lobatto) pour l'element de joint 86 (JOI3)
  571. *
  572. IF(MELE.EQ.86) THEN
  573. ETAGAU(1)=-UN
  574. ETAGAU(2)=-UN
  575. ETAGAU(3)=-UN
  576. QSIGAU(1)=-1.D0
  577. QSIGAU(3)= 1.D0
  578. POIGAU(1)= UNTIER
  579. POIGAU(2)= 2*DEUTIE
  580. POIGAU(3)= UNTIER
  581. ENDIF
  582.  
  583. GOTO 666
  584. C
  585. C 4 POINTS D INTEGRATION
  586. C
  587. 124 IF(NBG.NE.4) GOTO 66
  588. QSIGAU(1)=-X861
  589. QSIGAU(2)=-X339
  590. QSIGAU(3)= X339
  591. QSIGAU(4)= X861
  592. POIGAU(1)= P347
  593. POIGAU(2)= P652
  594. POIGAU(3)= P652
  595. POIGAU(4)= P347
  596. GOTO 666
  597. C
  598. C=======================================================================
  599. C ELEMENT CUBE VOLUME=8
  600. C=======================================================================
  601. C
  602. 14 IF(NBG.NE.1) GOTO 148
  603. C
  604. C UN POINT D INTEGRATION (ordre 1)
  605. C
  606. POIGAU(1)=HUIT
  607. GOTO 666
  608. 148 IF (NBG.NE.8) GOTO 1427
  609. C
  610. C 8 POINTS D INTEGRATION (ordre 3)
  611. C
  612. QSIGAU(1)=-X577
  613. QSIGAU(2)= X577
  614. QSIGAU(3)= X577
  615. QSIGAU(4)=-X577
  616. QSIGAU(5)=-X577
  617. QSIGAU(6)= X577
  618. QSIGAU(7)= X577
  619. QSIGAU(8)=-X577
  620. ETAGAU(1)=-X577
  621. ETAGAU(2)=-X577
  622. ETAGAU(3)= X577
  623. ETAGAU(4)= X577
  624. ETAGAU(5)=-X577
  625. ETAGAU(6)=-X577
  626. ETAGAU(7)= X577
  627. ETAGAU(8)= X577
  628. DO 141 IA=1,4
  629. POIGAU(IA )= UN
  630. POIGAU(IA+4)= UN
  631. DZEGAU(IA )=-X577
  632. DZEGAU(IA+4)= X577
  633. 141 CONTINUE
  634. GOTO 666
  635. C
  636. C 27 POINTS D INTEGRATION (ordre 5)
  637. C
  638. 1427 IF(NBG.NE.27) GOTO 1440
  639. DO 151 IA=1,25,3
  640. QSIGAU(IA )=-X774
  641. QSIGAU(IA+2)= X774
  642. 151 CONTINUE
  643. DO 152 IA=1,9
  644. DZEGAU(IA )=-X774
  645. DZEGAU(IA+18)= X774
  646. 152 CONTINUE
  647. DO 153 IA=1,3
  648. ETAGAU(IA )=-X774
  649. ETAGAU(IA+6 )= X774
  650. ETAGAU(IA+9 )=-X774
  651. ETAGAU(IA+15)= X774
  652. ETAGAU(IA+18)=-X774
  653. ETAGAU(IA+24)= X774
  654. 153 CONTINUE
  655. POIGAU(1)=P555*P555
  656. POIGAU(3)=POIGAU(1)
  657. POIGAU(7)=POIGAU(1)
  658. POIGAU(9)=POIGAU(1)
  659. POIGAU(2)=P888*P555
  660. POIGAU(4)=POIGAU(2)
  661. POIGAU(6)=POIGAU(2)
  662. POIGAU(8)=POIGAU(2)
  663. POIGAU(5)=P888*P888
  664. DO 154 IA=1,9
  665. XX=POIGAU(IA)
  666. POIGAU(IA )=XX*P555
  667. POIGAU(IA+9 )=XX*P888
  668. POIGAU(IA+18)=POIGAU(IA)
  669. 154 CONTINUE
  670. GO TO 666
  671. 1440 IF(NBG.NE.5) go to 66
  672. * special shb8
  673. QSIGAU(1) = -0.906179845938664D0
  674. QSIGAU(2) = -0.538469310105683D0
  675. QSIGAU(3) = 0.D0
  676. QSIGAU(4) = 0.538469310105683D0
  677. QSIGAU(5) = 0.906179845938664D0
  678. C
  679. POIGAU(1) = 0.236926885056189D0
  680. POIGAU(2) = 0.478628670499366D0
  681. POIGAU(3) = 0.568888888888889D0
  682. POIGAU(4) = 0.478628670499366D0
  683. POIGAU(5) = 0.236926885056189D0
  684. C
  685. GOTO 666
  686. C
  687. C=======================================================================
  688. C ELEMENT PRISME VOLUME=1
  689. C=======================================================================
  690. C
  691. C
  692. C UN POINT D INTEGRATION
  693. C
  694. 16 IF(NBG.NE.1) GOTO 162
  695. POIGAU(1)=UN
  696. QSIGAU(1)=UNTIER
  697. ETAGAU(1)=UNTIER
  698. DZEGAU(1)=0.
  699. GOTO 666
  700. C
  701. C 2 POINTS D INTEGRATION (ordre 1/3)
  702. C dans le plan (QSI,ETA), idem triangle a 1 pt de G --> ordre 1
  703. C selon DZE, idem segment a 2 ptG -> ordre 3
  704. C
  705. 162 IF(NBG.NE.2) GO TO 166
  706. QSIGAU(1)=UNTIER
  707. QSIGAU(2)=UNTIER
  708. ETAGAU(1)=UNTIER
  709. ETAGAU(2)=UNTIER
  710. DZEGAU(1)=-X577
  711. DZEGAU(2)= X577
  712. POIGAU(1)=UNDEMI
  713. POIGAU(2)=UNDEMI
  714. GOTO 666
  715. C
  716. C 6 POINTS D INTEGRATION (ordre 2/3)
  717. C dans le plan (QSI,ETA), idem triangle a 3 pt de G --> ordre 2
  718. C selon DZE, idem segment a 2 ptG -> ordre 3
  719. c
  720. 166 IF(NBG.NE.6) GO TO 168
  721. QSIGAU(1)=UNSIX
  722. ETAGAU(1)=UNSIX
  723. QSIGAU(2)=DEUTIE
  724. ETAGAU(2)=UNSIX
  725. QSIGAU(3)=UNSIX
  726. ETAGAU(3)=DEUTIE
  727. QSIGAU(4)=UNSIX
  728. ETAGAU(4)=UNSIX
  729. QSIGAU(5)=DEUTIE
  730. ETAGAU(5)=UNSIX
  731. QSIGAU(6)=UNSIX
  732. ETAGAU(6)=DEUTIE
  733. DO IA=1,3
  734. DZEGAU(IA )=-X577
  735. DZEGAU(IA+3)= X577
  736. ENDDO
  737. DO IA=1,6
  738. POIGAU(IA)=UNSIX
  739. ENDDO
  740. GOTO 666
  741. C
  742. C 8 POINTS D INTEGRATION (ordre 3/3)
  743. C (attention : poids negatif !)
  744. 168 IF(NBG.NE.8) GO TO 169
  745. QSIGAU(1)=UNTIER
  746. QSIGAU(2)=UNCINQ
  747. QSIGAU(3)=UNCINQ*TROIS
  748. QSIGAU(4)=UNCINQ
  749. ETAGAU(1)=UNTIER
  750. ETAGAU(2)=UNCINQ
  751. ETAGAU(3)=UNCINQ
  752. ETAGAU(4)=UNCINQ*TROIS
  753. POIGAU(1)=P281
  754. POIGAU(2)=P260
  755. POIGAU(3)=P260
  756. POIGAU(4)=P260
  757. DO IA=1,4
  758. QSIGAU(IA+4)=QSIGAU(IA)
  759. ETAGAU(IA+4)=ETAGAU(IA)
  760. POIGAU(IA+4)=POIGAU(IA)
  761. DZEGAU(IA) =-X577
  762. DZEGAU(IA+4)= X577
  763. ENDDO
  764. GOTO 666
  765. C
  766. C 9 POINTS D INTEGRATION
  767. C
  768. 169 IF (NBG.NE.9 ) GOTO 1612
  769. DO IA=1,7,3
  770. QSIGAU(IA )=UNDEMI
  771. ETAGAU(IA+1)=UNDEMI
  772. QSIGAU(IA+2)=UNDEMI
  773. ETAGAU(IA+2)=UNDEMI
  774. ENDDO
  775. DO IA=1,3
  776. DZEGAU(IA )=-X774
  777. POIGAU(IA )=P555*UNSIX
  778. POIGAU(IA+3)=P888*UNSIX
  779. POIGAU(IA+6)=P555*UNSIX
  780. DZEGAU(IA+6)= X774
  781. ENDDO
  782. GOTO 666
  783. C
  784. C 12 POINTS D INTEGRATION (ordre 3/5)
  785. C dans le plan (QSI,ETA), idem triangle a 4 pt de G que
  786. C [Zienkiewicz, p.165] -> ordre 3
  787. C selon DZE, idem segment a 3 ptG -> ordre 5
  788. 1612 IF (NBG.NE.12 ) GOTO 1621
  789. DO 171 IA=2,10,4
  790. QSIGAU(IA)=UNCINQ
  791. ETAGAU(IA)=UNCINQ
  792. 171 CONTINUE
  793. DO 172 IA=3,11,4
  794. QSIGAU(IA )=UNCINQ*TROIS
  795. ETAGAU(IA )=UNCINQ
  796. QSIGAU(IA+1)=UNCINQ
  797. ETAGAU(IA+1)=UNCINQ*TROIS
  798. 172 CONTINUE
  799. DO 173 IA=1,9,4
  800. QSIGAU(IA)=UNTIER
  801. ETAGAU(IA)=UNTIER
  802. 173 CONTINUE
  803. DO 174 IA=1,4
  804. DZEGAU(IA )=-X774
  805. DZEGAU(IA+8)= X774
  806. 174 CONTINUE
  807. DO 175 IA=2,4
  808. POIGAU(IA )=P260*P555
  809. POIGAU(IA+4)=P260*P888
  810. POIGAU(IA+8)=P260*P555
  811. 175 CONTINUE
  812. POIGAU(1)=P281*P555
  813. POIGAU(5)=P281*P888
  814. POIGAU(9)=P281*P555
  815. GO TO 666
  816. C
  817. C 21 POINTS D INTEGRATION (ordre 5)
  818. C dans le plan (QSI,ETA), idem triangle a 7 pt de G que
  819. C [Zienkiewicz, p.165] -> ordre 5
  820. C selon DZE, segment a 3 ptG -> ordre 5
  821. 1621 IF(NBG.NE.21) GO TO 66
  822. QSIGAU(1)=X101
  823. QSIGAU(2)=X797
  824. QSIGAU(3)=X101
  825. QSIGAU(4)=X470
  826. QSIGAU(5)=X470
  827. QSIGAU(6)=X059
  828. QSIGAU(7)=UNTIER
  829. ETAGAU(1)=QSIGAU(1)
  830. ETAGAU(2)=QSIGAU(1)
  831. ETAGAU(3)=QSIGAU(2)
  832. ETAGAU(4)=QSIGAU(6)
  833. ETAGAU(5)=QSIGAU(4)
  834. ETAGAU(6)=QSIGAU(4)
  835. ETAGAU(7)=QSIGAU(7)
  836. POIGAU(1)=P125*UNDEMI
  837. POIGAU(2)=P125*UNDEMI
  838. POIGAU(3)=P125*UNDEMI
  839. POIGAU(4)=P132*UNDEMI
  840. POIGAU(5)=P132*UNDEMI
  841. POIGAU(6)=P132*UNDEMI
  842. POIGAU(7)=.1125D0
  843. do 176 iu=1,7
  844. dzegau(iu+14)=X774
  845. dzegau(iu+7)=0.d0
  846. 176 dzegau(iu)=-X774
  847. do 177 iu=1,7
  848. poigau(iu+7)= poigau(iu)*P888
  849. poigau(iu)=poigau(iu)*p555
  850. poigau (iu+14)=poigau(iu)
  851. 177 continue
  852. do 178 iu=1,7
  853. qsigau(iu+7)=qsigau(iu)
  854. qsigau(iu+14)=qsigau(iu)
  855. etagau(iu+7)= etagau(iu)
  856. etagau(iu+14)=etagau(iu)
  857. 178 continue
  858. go to 666
  859. C
  860. C=======================================================================
  861. C ELEMENT LIAISON 3
  862. C=======================================================================
  863. C
  864. 18 CONTINUE
  865. C
  866. C 1 POINTS D INTEGRATION
  867. C
  868. IF (NBG.NE.1) GOTO 183
  869. QSIGAU(1)=UNTIER
  870. ETAGAU(1)=UNTIER
  871. POIGAU(1)=UNDEMI
  872. GOTO 666
  873. C
  874. C 3 POINTS D INTEGRATION ZIENK P201
  875. C
  876. 183 IF(NBG.NE.3) GOTO 184
  877. QSIGAU(1)=UNDEMI
  878. ETAGAU(2)=UNDEMI
  879. QSIGAU(3)=UNDEMI
  880. ETAGAU(3)=UNDEMI
  881. POIGAU(1)=UNSIX
  882. POIGAU(2)=UNSIX
  883. POIGAU(3)=UNSIX
  884. GOTO 666
  885. C
  886. C 4 POINTS D INTEGRATION ZIENK P201
  887. C
  888. 184 IF (NBG.NE.4) GOTO 187
  889. QSIGAU(1)=UNTIER
  890. QSIGAU(2)=UNCINQ
  891. QSIGAU(3)=UNCINQ*TROIS
  892. QSIGAU(4)=UNCINQ
  893. ETAGAU(1)=UNTIER
  894. ETAGAU(2)=UNCINQ
  895. ETAGAU(3)=UNCINQ
  896. ETAGAU(4)=UNCINQ*TROIS
  897. POIGAU(1)=-27.D0/96.D0
  898. POIGAU(2)= 25.D0/96.D0
  899. POIGAU(3)=POIGAU(2)
  900. POIGAU(4)=POIGAU(2)
  901. GOTO 666
  902. C
  903. C 7 POINTS D INTEGRATION BATHE P280
  904. C
  905. 187 IF(NBG.NE.7) GOTO 66
  906. QSIGAU(1)=X101
  907. QSIGAU(2)=X797
  908. QSIGAU(3)=X101
  909. QSIGAU(4)=X470
  910. QSIGAU(5)=X470
  911. QSIGAU(6)=X059
  912. QSIGAU(7)=UNTIER
  913. ETAGAU(1)=QSIGAU(1)
  914. ETAGAU(2)=QSIGAU(1)
  915. ETAGAU(3)=QSIGAU(2)
  916. ETAGAU(4)=QSIGAU(6)
  917. ETAGAU(5)=QSIGAU(4)
  918. ETAGAU(6)=QSIGAU(4)
  919. ETAGAU(7)=QSIGAU(7)
  920. POIGAU(1)=P125*UNDEMI
  921. POIGAU(2)=P125*UNDEMI
  922. POIGAU(3)=P125*UNDEMI
  923. POIGAU(4)=P132*UNDEMI
  924. POIGAU(5)=P132*UNDEMI
  925. POIGAU(6)=P132*UNDEMI
  926. POIGAU(7)=.1125D0
  927. GOTO 666
  928. C
  929. C=======================================================================
  930. C ELEMENT LIAISON 4
  931. C=======================================================================
  932. C
  933. 19 IF(NBG.NE.1) GOTO 194
  934. C
  935. C UN POINT D INTEGRATION
  936. C
  937. POIGAU(1)=QUATRE
  938. GOTO 666
  939. 194 IF(NBG.NE.4) GOTO 199
  940. C
  941. C 4 POINTS D INTEGRATION
  942. C
  943. QSIGAU(1)=-X577
  944. QSIGAU(2)= X577
  945. QSIGAU(3)= X577
  946. QSIGAU(4)=-X577
  947. ETAGAU(1)=-X577
  948. ETAGAU(2)=-X577
  949. ETAGAU(3)= X577
  950. ETAGAU(4)= X577
  951. POIGAU(1)=UN
  952. POIGAU(2)=UN
  953. POIGAU(3)=UN
  954. POIGAU(4)=UN
  955. GOTO 666
  956. C
  957. C 9 POINTS D INTEGRATION
  958. C
  959. 199 IF (NBG.NE.9 ) GOTO 1916
  960. QSIGAU(1)=-X774
  961. QSIGAU(3)= X774
  962. QSIGAU(4)=-X774
  963. QSIGAU(6)= X774
  964. QSIGAU(7)=-X774
  965. QSIGAU(9)= X774
  966. ETAGAU(1)=-X774
  967. ETAGAU(2)=-X774
  968. ETAGAU(3)=-X774
  969. ETAGAU(7)= X774
  970. ETAGAU(8)= X774
  971. ETAGAU(9)= X774
  972. POIGAU(1)=P555*P555
  973. POIGAU(3)=POIGAU(1)
  974. POIGAU(7)=POIGAU(1)
  975. POIGAU(9)=POIGAU(1)
  976. POIGAU(2)=P888*P555
  977. POIGAU(4)=POIGAU(2)
  978. POIGAU(6)=POIGAU(2)
  979. POIGAU(8)=POIGAU(2)
  980. POIGAU(5)=P888*P888
  981. GOTO 666
  982. 1916 IF(NBG.NE.16) GOTO 66
  983. C
  984. C CARRE 16 POINTS D INTEGRATION
  985. C
  986. DO 19161 IA=1,13,4
  987. QSIGAU(IA )=-X861
  988. QSIGAU(IA+1 )=-X339
  989. QSIGAU(IA+2 )= X339
  990. QSIGAU(IA+3 )= X861
  991. POIGAU(IA )= P347
  992. POIGAU(IA+1 )= P652
  993. POIGAU(IA+2 )= P652
  994. POIGAU(IA+3 )= P347
  995. 19161 CONTINUE
  996. DO 19162 IA=1,4
  997. ETAGAU(IA )=-X861
  998. ETAGAU(IA+4 )=-X339
  999. ETAGAU(IA+8 )= X339
  1000. ETAGAU(IA+12)= X861
  1001. POIGAU(IA )= P347*POIGAU(IA )
  1002. POIGAU(IA+4 )= P652*POIGAU(IA+4 )
  1003. POIGAU(IA+8 )= P652*POIGAU(IA+8 )
  1004. POIGAU(IA+12)= P347*POIGAU(IA+12)
  1005. 19162 CONTINUE
  1006. GOTO 666
  1007. C
  1008. C=======================================================================
  1009. C ELEMENT TETRAEDRE VOLUME=1/6
  1010. C=======================================================================
  1011. C
  1012. 23 IF (NBG.NE.1) GOTO 234
  1013. C
  1014. C 1 POINT D INTEGRATION (ordre 1)
  1015. C
  1016. QSIGAU(1)=UNQUA
  1017. ETAGAU(1)=UNQUA
  1018. DZEGAU(1)=UNQUA
  1019. POIGAU(1)=UNSIX
  1020. GOTO 666
  1021. C
  1022. C 4 POINTS D INTEGRATION (ordre 2)
  1023. c [Zienkiewicz, p.166]
  1024. C
  1025. 234 IF (NBG.NE.4) GOTO 235
  1026. QSIGAU(1)=X138
  1027. QSIGAU(2)=X585
  1028. QSIGAU(3)=X138
  1029. QSIGAU(4)=X138
  1030. ETAGAU(1)=X138
  1031. ETAGAU(2)=X138
  1032. ETAGAU(3)=X585
  1033. ETAGAU(4)=X138
  1034. DZEGAU(1)=X138
  1035. DZEGAU(2)=X138
  1036. DZEGAU(3)=X138
  1037. DZEGAU(4)=X585
  1038. POIGAU(1)=UN24
  1039. POIGAU(2)=UN24
  1040. POIGAU(3)=UN24
  1041. POIGAU(4)=UN24
  1042. GOTO 666
  1043. C
  1044. C 5 POINTS D INTEGRATION (ordre 3)
  1045. c [Zienkiewicz, p.166] (attention : poids negatif !)
  1046. C
  1047. 235 IF(NBG.NE.5) GOTO 236
  1048. QSIGAU(1)=UNQUA
  1049. QSIGAU(2)=UNSIX
  1050. QSIGAU(3)=UNSIX
  1051. QSIGAU(4)=UNSIX
  1052. QSIGAU(5)=UNDEMI
  1053. ETAGAU(1)=UNQUA
  1054. ETAGAU(2)=UNSIX
  1055. ETAGAU(3)=UNSIX
  1056. ETAGAU(4)=UNDEMI
  1057. ETAGAU(5)=UNSIX
  1058. DZEGAU(1)=UNQUA
  1059. DZEGAU(2)=UNSIX
  1060. DZEGAU(3)=UNDEMI
  1061. DZEGAU(4)=UNSIX
  1062. DZEGAU(5)=UNSIX
  1063. POIGAU(1)=-.80D0*UNSIX
  1064. POIGAU(2)= .45D0*UNSIX
  1065. POIGAU(3)= .45D0*UNSIX
  1066. POIGAU(4)= .45D0*UNSIX
  1067. POIGAU(5)= .45D0*UNSIX
  1068. GOTO 666
  1069. C
  1070. C 11 POINTS D INTEGRATION (ordre 4)
  1071. C [DHATT, FEM, p371] (attention : poids negatif !)
  1072. C
  1073. 236 IF(NBG.NE.11) GOTO 237
  1074. * 1er point = le centre
  1075. QSIGAU(1)=UNQUA
  1076. ETAGAU(1)=UNQUA
  1077. DZEGAU(1)=UNQUA
  1078. POIGAU(1)=P01315
  1079. * points 2 a 5
  1080. QSIGAU(2)=XONZ14
  1081. ETAGAU(3)=XONZ14
  1082. DZEGAU(4)=XONZ14
  1083. ETAGAU(2)=XUN14
  1084. DZEGAU(2)=XUN14
  1085. QSIGAU(3)=XUN14
  1086. DZEGAU(3)=XUN14
  1087. QSIGAU(4)=XUN14
  1088. ETAGAU(4)=XUN14
  1089. QSIGAU(5)=XUN14
  1090. ETAGAU(5)=XUN14
  1091. DZEGAU(5)=XUN14
  1092. POIGAU(2)=P00762
  1093. POIGAU(3)=P00762
  1094. POIGAU(4)=P00762
  1095. POIGAU(5)=P00762
  1096. * points 6 a 11
  1097. QSIGAU( 6)=UNQUA
  1098. ETAGAU( 6)=UNQUA
  1099. DZEGAU( 6)=UNQUA
  1100. QSIGAU( 7)=UNQUA
  1101. DZEGAU( 7)=UNQUA
  1102. QSIGAU( 8)=UNQUA
  1103. ETAGAU( 9)=UNQUA
  1104. DZEGAU( 9)=UNQUA
  1105. DZEGAU(10)=UNQUA
  1106. ETAGAU(11)=UNQUA
  1107. ETAGAU( 7)=UNSIX
  1108. ETAGAU( 8)=UNSIX
  1109. DZEGAU( 8)=UNSIX
  1110. QSIGAU( 9)=UNSIX
  1111. QSIGAU(10)=UNSIX
  1112. ETAGAU(10)=UNSIX
  1113. QSIGAU(11)=UNSIX
  1114. DZEGAU(11)=UNSIX
  1115. POIGAU(6)=P00124
  1116. POIGAU(7)=P00124
  1117. POIGAU(8)=P00124
  1118. POIGAU(9)=P00124
  1119. POIGAU(10)=P00124
  1120. POIGAU(11)=P00124
  1121. GOTO 666
  1122. C
  1123. C 15 POINTS D INTEGRATION (ordre 5)
  1124. C [DHATT, FEM, p371]
  1125. C
  1126. 237 IF(NBG.NE.15) GOTO 66
  1127. * 1er point = le centre
  1128. QSIGAU(1)=UNQUA
  1129. ETAGAU(1)=UNQUA
  1130. DZEGAU(1)=UNQUA
  1131. POIGAU(1)=P0197
  1132. * points 2 a 5
  1133. QSIGAU(2)=X3197
  1134. ETAGAU(2)=X3197
  1135. DZEGAU(2)=X3197
  1136. c QSIGAU(3)=X7240
  1137. QSIGAU(3)=X0406
  1138. ETAGAU(3)=X3197
  1139. DZEGAU(3)=X3197
  1140. QSIGAU(4)=X3197
  1141. c ETAGAU(4)=X7240
  1142. ETAGAU(4)=X0406
  1143. DZEGAU(4)=X3197
  1144. QSIGAU(5)=X3197
  1145. ETAGAU(5)=X3197
  1146. c DZEGAU(5)=X7240
  1147. DZEGAU(5)=X0406
  1148. c POIGAU(2)=P0119
  1149. c POIGAU(3)=P0119
  1150. c POIGAU(4)=P0119
  1151. c POIGAU(5)=P0119
  1152. POIGAU(2)=P0115
  1153. POIGAU(3)=P0115
  1154. POIGAU(4)=P0115
  1155. POIGAU(5)=P0115
  1156. * points 6 a 9
  1157. QSIGAU(6)=X0919
  1158. ETAGAU(6)=X0919
  1159. DZEGAU(6)=X0919
  1160. c QSIGAU(7)=X0406
  1161. QSIGAU(7)=X7240
  1162. ETAGAU(7)=X0919
  1163. DZEGAU(7)=X0919
  1164. QSIGAU(8)=X0919
  1165. c ETAGAU(8)=X0406
  1166. ETAGAU(8)=X7240
  1167. DZEGAU(8)=X0919
  1168. QSIGAU(9)=X0919
  1169. ETAGAU(9)=X0919
  1170. c DZEGAU(9)=X0406
  1171. DZEGAU(9)=X7240
  1172. c POIGAU(6)=P0115
  1173. c POIGAU(7)=P0115
  1174. c POIGAU(8)=P0115
  1175. c POIGAU(9)=P0115
  1176. POIGAU(6)=P0119
  1177. POIGAU(7)=P0119
  1178. POIGAU(8)=P0119
  1179. POIGAU(9)=P0119
  1180. * points 10 a 15
  1181. QSIGAU(10)=X4436
  1182. ETAGAU(10)=X0563
  1183. DZEGAU(10)=X0563
  1184. QSIGAU(11)=X0563
  1185. ETAGAU(11)=X4436
  1186. DZEGAU(11)=X0563
  1187. QSIGAU(12)=X0563
  1188. ETAGAU(12)=X0563
  1189. DZEGAU(12)=X4436
  1190. QSIGAU(13)=X0563
  1191. ETAGAU(13)=X4436
  1192. DZEGAU(13)=X4436
  1193. QSIGAU(14)=X4436
  1194. ETAGAU(14)=X0563
  1195. DZEGAU(14)=X4436
  1196. QSIGAU(15)=X4436
  1197. ETAGAU(15)=X4436
  1198. DZEGAU(15)=X0563
  1199. POIGAU(10)=P0088
  1200. POIGAU(11)=P0088
  1201. POIGAU(12)=P0088
  1202. POIGAU(13)=P0088
  1203. POIGAU(14)=P0088
  1204. POIGAU(15)=P0088
  1205. GOTO 666
  1206. C
  1207. C=======================================================================
  1208. C ELEMENT PYRAMIDE VOLUME=2/3
  1209. C=======================================================================
  1210. C
  1211. 25 IF(NBG.NE.1) GOTO 255
  1212. C
  1213. C UN POINT D INTEGRATION (ordre 1)
  1214. C
  1215. DZEGAU(1)=UNQUA
  1216. POIGAU(1)=DEUTIE
  1217. GOTO 666
  1218. C
  1219. C 5 POINTS D INTEGRATION (ordre 2?)
  1220. C regle d'integration retrouvee par [Kubatko et al.,FEAD, 2013]
  1221. 255 IF(NBG.NE.5) GOTO 256
  1222. AUX= 0.5D0
  1223. POI1=0.1333333333333333D0
  1224. POI2= 0.1333333333333333D0
  1225. H1 = 0.1531754163448146D0
  1226. H2 = 0.6372983346207416D0
  1227. QSIGAU(1)= AUX
  1228. QSIGAU(2)= XZERO
  1229. QSIGAU(3)=-AUX
  1230. QSIGAU(4)= XZERO
  1231. QSIGAU(5)= XZERO
  1232. ETAGAU(1)= XZERO
  1233. ETAGAU(2)= AUX
  1234. ETAGAU(3)= XZERO
  1235. ETAGAU(4)=-AUX
  1236. ETAGAU(5)= XZERO
  1237. DO 251 IA=1,4
  1238. DZEGAU(IA)= H1
  1239. POIGAU(IA)= POI1
  1240. 251 CONTINUE
  1241. DZEGAU(5)= H2
  1242. POIGAU(5)= POI2
  1243. GOTO 666
  1244. C
  1245. C 6 POINTS D INTEGRATION (???)
  1246. C
  1247. 256 IF (NBG.NE.6) GOTO 2510
  1248. AUX =.5702963741068025D0
  1249. POI1=.1024890634400000D0
  1250. POI2=.1100000000000000D0
  1251. POI3=.1467104129066667D0
  1252. H1 =.1666666666666667D0
  1253. H2 =.8063183038464675D-1
  1254. H3 =.6098484849057127D0
  1255. QSIGAU(1)= AUX
  1256. QSIGAU(2)= XZERO
  1257. QSIGAU(3)=-AUX
  1258. QSIGAU(4)= XZERO
  1259. QSIGAU(5)= XZERO
  1260. QSIGAU(6)= XZERO
  1261. ETAGAU(1)= XZERO
  1262. ETAGAU(2)= AUX
  1263. ETAGAU(3)= XZERO
  1264. ETAGAU(4)=-AUX
  1265. ETAGAU(5)= XZERO
  1266. ETAGAU(6)= XZERO
  1267. DO 2516 IA=1,4
  1268. DZEGAU(IA)= H1
  1269. POIGAU(IA)= POI1
  1270. 2516 CONTINUE
  1271. DZEGAU(5)= H2
  1272. DZEGAU(6)= H3
  1273. POIGAU(5)= POI2
  1274. POIGAU(6)= POI3
  1275. GOTO 666
  1276. c
  1277. C rem bp : pour 10 et 17 points, cf. [Kubatko et al.,FEAD, 2013]
  1278. C
  1279. C 10 POINTS D INTEGRATION (ordre 4?)
  1280. C
  1281. 2510 IF (NBG.NE.10) GOTO 2517
  1282. * points 1 a 4
  1283. QSIGAU(1)=-0.657966997126504D0
  1284. QSIGAU(2)= 0.657966997126504D0
  1285. ETAGAU(3)=-0.657966997126504D0
  1286. ETAGAU(4)= 0.657966997126504D0
  1287. DO IA=1,4
  1288. DZEGAU(IA)=0.0392482838982998D0
  1289. POIGAU(IA)=0.0423460604458543D0
  1290. ENDDO
  1291. * points 5
  1292. DZEGAU(5)=0.12513695309945289D0
  1293. POIGAU(5)=0.13792226838670674D0
  1294. * points 6 a 9
  1295. QSIGAU(6)=-0.325290778196272556D0
  1296. QSIGAU(7)= 0.325290778196272556D0
  1297. QSIGAU(8)= 0.325290778196272556D0
  1298. QSIGAU(9)=-0.325290778196272556D0
  1299. ETAGAU(6)=-0.325290778196272556D0
  1300. ETAGAU(7)=-0.325290778196272556D0
  1301. ETAGAU(8)= 0.325290778196272556D0
  1302. ETAGAU(9)= 0.325290778196272556D0
  1303. DO IA=6,9
  1304. DZEGAU(IA)=0.32238414957239719D0
  1305. POIGAU(IA)=0.07088305859532487D0
  1306. ENDDO
  1307. * points 10
  1308. DZEGAU(10)=0.677232788884158832D0
  1309. POIGAU(10)=0.07582792211524662D0
  1310. GOTO 666
  1311. C
  1312. C 17 POINTS D INTEGRATION (ordre 5?)
  1313. C
  1314. 2517 IF (NBG.NE.17) GOTO 258
  1315. * points 1 a 4
  1316. QSIGAU(1)=-0.3196022296916385D0
  1317. ETAGAU(1)=-0.3196022296916385D0
  1318. QSIGAU(2)=0.3196022296916385D0
  1319. ETAGAU(2)=-0.3196022296916385D0
  1320. QSIGAU(3)=0.3196022296916385D0
  1321. ETAGAU(3)=0.3196022296916385D0
  1322. QSIGAU(4)=-0.3196022296916385D0
  1323. ETAGAU(4)=0.3196022296916385D0
  1324. DO IA=1,4
  1325. DZEGAU(IA)=0.0634862039888108D0
  1326. POIGAU(IA)=0.0576400632904598D0
  1327. ENDDO
  1328. * points 5 a 8
  1329. QSIGAU(5)=0.7980415331554783D0
  1330. QSIGAU(6)=-0.7980415331554783D0
  1331. ETAGAU(7)=0.7980415331554783D0
  1332. ETAGAU(8)=-0.7980415331554783D0
  1333. DO IA=5,8
  1334. DZEGAU(IA)=0.0954161323042548D0
  1335. POIGAU(IA)=0.0181053421609235D0
  1336. ENDDO
  1337. * points 9 a 12
  1338. QSIGAU(9)=-0.2341903767177774D0
  1339. ETAGAU(9)=-0.2341903767177774D0
  1340. QSIGAU(10)=0.2341903767177774D0
  1341. ETAGAU(10)=-0.2341903767177774D0
  1342. QSIGAU(11)=0.2341903767177774D0
  1343. ETAGAU(11)=0.2341903767177774D0
  1344. QSIGAU(12)=-0.2341903767177774D0
  1345. ETAGAU(12)=0.2341903767177774D0
  1346. DO IA=9,12
  1347. DZEGAU(IA)=0.3174282137715947D0
  1348. POIGAU(IA)=0.0639133011727061D0
  1349. ENDDO
  1350. * points 13 a 16
  1351. QSIGAU(13)=0.4433429394393417D0
  1352. QSIGAU(14)=-0.4433429394393417D0
  1353. ETAGAU(15)=0.4433429394393417D0
  1354. ETAGAU(16)=-0.4433429394393417D0
  1355. DO IA=13,16
  1356. DZEGAU(IA)=0.4914145363689955D0
  1357. POIGAU(IA)=0.0153472930269660D0
  1358. ENDDO
  1359. * point 17
  1360. DZEGAU(17)=0.7246590783754809D0
  1361. POIGAU(17)=0.0466426680624461D0
  1362. GOTO 666
  1363. C
  1364. C rem bp : pour 8 et 27 points, methode decrite dans [Bedrosian, 1992]
  1365. c elle-meme basee sur [Hammer, Marlow et Stroud, 1957]
  1366. C valeurs generees automatiquement par un dgibi bien ecrit ;)
  1367. C
  1368. C 8 POINTS D INTEGRATION (ordre 3)
  1369. C
  1370. 258 IF (NBG.NE.8) GOTO 2527
  1371. QSIGAU(1)=-.2631840555697141D0
  1372. DZEGAU(1)=0.5441518440112245D0
  1373. POIGAU(1)=0.0503929410399128D0
  1374. ETAGAU(2)=0.2631840555697141D0
  1375. DZEGAU(2)=0.5441518440112245D0
  1376. POIGAU(2)=0.0503929410399128D0
  1377. ETAGAU(3)=-.2631840555697141D0
  1378. DZEGAU(3)=0.5441518440112245D0
  1379. POIGAU(3)=0.0503929410399128D0
  1380. QSIGAU(4)=0.2631840555697141D0
  1381. DZEGAU(4)=0.5441518440112245D0
  1382. POIGAU(4)=0.0503929410399128D0
  1383. QSIGAU(5)=-.5066163033497884D0
  1384. DZEGAU(5)=0.1225148226554399D0
  1385. POIGAU(5)=0.1162737256267542D0
  1386. ETAGAU(6)=0.5066163033497884D0
  1387. DZEGAU(6)=0.1225148226554399D0
  1388. POIGAU(6)=0.1162737256267542D0
  1389. ETAGAU(7)=-.5066163033497884D0
  1390. DZEGAU(7)=0.1225148226554399D0
  1391. POIGAU(7)=0.1162737256267542D0
  1392. QSIGAU(8)=0.5066163033497884D0
  1393. DZEGAU(8)=0.1225148226554399D0
  1394. POIGAU(8)=0.1162737256267542D0
  1395. GOTO 666
  1396. C
  1397. C 27 POINTS D INTEGRATION (ordre 5)
  1398. C
  1399. 2527 IF(NBG.NE.27) GO TO 66
  1400. c rem bp : on a remplace l'ancienne methode qui semblait erronee!
  1401. QSIGAU(1)=-.2285043056539678D0
  1402. DZEGAU(1)=0.7050022098884978D0
  1403. POIGAU(1)=0.0046220220692254D0
  1404. QSIGAU(2)=-.1142521528269839D0
  1405. ETAGAU(2)=0.1142521528269839D0
  1406. DZEGAU(2)=0.7050022098884978D0
  1407. POIGAU(2)=0.0073952353107607D0
  1408. ETAGAU(3)=0.2285043056539678D0
  1409. DZEGAU(3)=0.7050022098884978D0
  1410. POIGAU(3)=0.0046220220692254D0
  1411. QSIGAU(4)=-.1142521528269839D0
  1412. ETAGAU(4)=-.1142521528269839D0
  1413. DZEGAU(4)=0.7050022098884978D0
  1414. POIGAU(4)=0.0073952353107607D0
  1415. DZEGAU(5)=0.7050022098884978D0
  1416. POIGAU(5)=0.0118323764972171D0
  1417. QSIGAU(6)=0.1142521528269839D0
  1418. ETAGAU(6)=0.1142521528269839D0
  1419. DZEGAU(6)=0.7050022098884978D0
  1420. POIGAU(6)=0.0073952353107607D0
  1421. ETAGAU(7)=-.2285043056539678D0
  1422. DZEGAU(7)=0.7050022098884978D0
  1423. POIGAU(7)=0.0046220220692254D0
  1424. QSIGAU(8)=0.1142521528269839D0
  1425. ETAGAU(8)=-.1142521528269839D0
  1426. DZEGAU(8)=0.7050022098884978D0
  1427. POIGAU(8)=0.0073952353107607D0
  1428. QSIGAU(9)=0.2285043056539678D0
  1429. DZEGAU(9)=0.7050022098884978D0
  1430. POIGAU(9)=0.0046220220692254D0
  1431. QSIGAU(10)=-.5058087078539260D0
  1432. DZEGAU(10)=0.3470037660383507D0
  1433. POIGAU(10)=0.0225688687129423D0
  1434. QSIGAU(11)=-.2529043539269630D0
  1435. ETAGAU(11)=0.2529043539269630D0
  1436. DZEGAU(11)=0.3470037660383507D0
  1437. POIGAU(11)=0.0361101899407077D0
  1438. ETAGAU(12)=0.5058087078539260D0
  1439. DZEGAU(12)=0.3470037660383507D0
  1440. POIGAU(12)=0.0225688687129423D0
  1441. QSIGAU(13)=-.2529043539269630D0
  1442. ETAGAU(13)=-.2529043539269630D0
  1443. DZEGAU(13)=0.3470037660383507D0
  1444. POIGAU(13)=0.0361101899407077D0
  1445. DZEGAU(14)=0.3470037660383507D0
  1446. POIGAU(14)=0.0577763039051324D0
  1447. QSIGAU(15)=0.2529043539269630D0
  1448. ETAGAU(15)=0.2529043539269630D0
  1449. DZEGAU(15)=0.3470037660383507D0
  1450. POIGAU(15)=0.0361101899407077D0
  1451. ETAGAU(16)=-.5058087078539260D0
  1452. DZEGAU(16)=0.3470037660383507D0
  1453. POIGAU(16)=0.0225688687129423D0
  1454. QSIGAU(17)=0.2529043539269630D0
  1455. ETAGAU(17)=-.2529043539269630D0
  1456. DZEGAU(17)=0.3470037660383507D0
  1457. POIGAU(17)=0.0361101899407077D0
  1458. QSIGAU(18)=0.5058087078539260D0
  1459. DZEGAU(18)=0.3470037660383507D0
  1460. POIGAU(18)=0.0225688687129423D0
  1461. QSIGAU(19)=-.7180557413198901D0
  1462. DZEGAU(19)=0.0729940240731481D0
  1463. POIGAU(19)=0.0242494384359393D0
  1464. QSIGAU(20)=-.3590278706599451D0
  1465. ETAGAU(20)=0.3590278706599451D0
  1466. DZEGAU(20)=0.0729940240731481D0
  1467. POIGAU(20)=0.0387991014975029D0
  1468. ETAGAU(21)=0.7180557413198901D0
  1469. DZEGAU(21)=0.0729940240731481D0
  1470. POIGAU(21)=0.0242494384359393D0
  1471. QSIGAU(22)=-.3590278706599451D0
  1472. ETAGAU(22)=-.3590278706599451D0
  1473. DZEGAU(22)=0.0729940240731481D0
  1474. POIGAU(22)=0.0387991014975029D0
  1475. DZEGAU(23)=0.0729940240731481D0
  1476. POIGAU(23)=0.0620785623960047D0
  1477. QSIGAU(24)=0.3590278706599451D0
  1478. ETAGAU(24)=0.3590278706599451D0
  1479. DZEGAU(24)=0.0729940240731481D0
  1480. POIGAU(24)=0.0387991014975029D0
  1481. ETAGAU(25)=-.7180557413198901D0
  1482. DZEGAU(25)=0.0729940240731481D0
  1483. POIGAU(25)=0.0242494384359393D0
  1484. QSIGAU(26)=0.3590278706599451D0
  1485. ETAGAU(26)=-.3590278706599451D0
  1486. DZEGAU(26)=0.0729940240731481D0
  1487. POIGAU(26)=0.0387991014975029D0
  1488. QSIGAU(27)=0.7180557413198901D0
  1489. DZEGAU(27)=0.0729940240731481D0
  1490. POIGAU(27)=0.0242494384359393D0
  1491. GOTO 666
  1492. C
  1493. 30 CONTINUE
  1494. C
  1495. C=======================================================================
  1496. C ELEMENT POLYGONE A NNO COTES
  1497. C=======================================================================
  1498. C
  1499. POIDS = COS(XPI/NNO) * SIN(XPI/NNO)
  1500. C
  1501. IF ( NBG .EQ.1) THEN
  1502. C
  1503. QSIGAU(1) = XZERO
  1504. ETAGAU(1) = XZERO
  1505. POIGAU(1) = POIDS * NNO
  1506. C
  1507. ELSE
  1508. C
  1509. SCOS2 = XZERO
  1510. SCOS4 = XZERO
  1511. SQSI2 = XZERO
  1512. SQSI4 = XZERO
  1513. TGPSN = TAN(XPI/NNO)
  1514. TG2PSN = TGPSN * TGPSN
  1515. C
  1516. DO 300, IA=1, NNO
  1517. C
  1518. ALPHA = (2*IA - 1) * XPI / NNO
  1519. BETA = 2 * (IA - 1) * XPI / NNO
  1520. CIA = COS(ALPHA)
  1521. CIA2 = CIA*CIA
  1522. CIB = COS(BETA)
  1523. CIB2 = CIB*CIB
  1524. SIB = SIN(BETA)
  1525. SIB2 = SIB*SIB
  1526. C
  1527. C Calcul de la somme des abscisses carré et puissance 4
  1528. C des points de GAUSS
  1529. C
  1530. SCOS2 = SCOS2+CIA2
  1531. SCOS4 = SCOS4+CIA2*CIA2
  1532. C
  1533. C Calcul de la somme des integrales sur chaque triangle
  1534. C de QSI ** 2
  1535. C
  1536. SQSI2 = SQSI2 + 3 * CIB2 + TG2PSN * SIB2
  1537. C
  1538. C Calcul de la somme des integrales sur chaque triangle
  1539. C de QSI ** 4
  1540. C
  1541. SQSI4 = SQSI4 + CIB2*CIB2 + 2*TG2PSN*CIB2*SIB2 +
  1542. & TG2PSN*TG2PSN*SIB2*SIB2/5
  1543. C
  1544. 300 CONTINUE
  1545. C
  1546. SOM2 = SQSI2 * COS (XPI/NNO) ** 4 * TGPSN / 6
  1547. SOM4 = SQSI4 * COS (XPI/NNO) ** 6 * TGPSN / 3
  1548. C
  1549. IF (NBG .EQ. NNO) THEN
  1550. C
  1551. C Recherche du poids des points de GAUSS N points
  1552. C
  1553. PE = POIDS
  1554. C
  1555. C Distance relative des points de GAUSS au centre du polygome
  1556. C
  1557. DIST = SQRT ( SOM2/(POIDS*SCOS2))
  1558. C
  1559. IF (NNO.EQ.5) DIST=DIST * 1.0080508208396880D0
  1560. IF (NNO.EQ.6) DIST=DIST * 1.0352424358375543D0
  1561. IF (NNO.EQ.7) DIST=DIST * 1.04449199849566D0
  1562. IF (NNO.EQ.8) DIST= 0.75741767161339812D0
  1563. IF (NNO.EQ.9) DIST=DIST * 1.13045302559206D0
  1564. IF (NNO.EQ.10) DIST= DIST * 1.151017710991479D0
  1565. IF (NNO.EQ.12) DIST= DIST * 1.11993200981186D0
  1566.  
  1567. ELSEIF (NBG.EQ.NNO+1) THEN
  1568. C
  1569. C Recherche du poids des points de GAUSS N+1 points
  1570. C
  1571. PE = SOM2*SOM2*SCOS4/(SOM4*SCOS2*SCOS2)
  1572. PC = NNO*(POIDS-PE)
  1573. C
  1574. C Distance relative des points de GAUSS au centre du polynome
  1575. C
  1576. DIST = SQRT (SOM4*SCOS2/(SOM2*SCOS4))
  1577. C
  1578. C ABCSISSE ET POIDS DU POINT CENTRAL
  1579. C
  1580. QSIGAU(NBG) = XZERO
  1581. ETAGAU(NBG) = XZERO
  1582. POIGAU(NBG) = PC
  1583. C
  1584. ELSE
  1585. GOTO 66
  1586. ENDIF
  1587. C
  1588. DO 305, IA = 1, NNO
  1589. C
  1590. ALPHA = (2*IA - 1) * XPI / NNO
  1591. C
  1592. QSIGAU(IA) = DIST * COS(ALPHA)
  1593. ETAGAU(IA) = DIST * SIN(ALPHA)
  1594. POIGAU(IA) = PE
  1595. C
  1596. 305 CONTINUE
  1597. C
  1598. ENDIF
  1599. GOTO 666
  1600. C
  1601. 95 CONTINUE
  1602. C
  1603. C=======================================================================
  1604. C ELEMENTS DE COQUE INTEGRES
  1605. C=======================================================================
  1606. C
  1607. IF(MELE.EQ.28.OR.MELE.EQ.93) THEN
  1608. NNAPPE=NBG/3
  1609. POIDS=UNSIX
  1610. DO 7772 IA=1,NNAPPE
  1611. NP=IA*3-2
  1612. POIGAU(NP )=POIDS
  1613. POIGAU(NP+1)=POIDS
  1614. POIGAU(NP+2)=POIDS
  1615. QSIGAU(NP )= UNDEMI
  1616. ETAGAU(NP+1)= UNDEMI
  1617. QSIGAU(NP+2)= UNDEMI
  1618. ETAGAU(NP+2)= UNDEMI
  1619. XXXX = DEUX*((DBLE(IA)-UN)/(DBLE(NNAPPE)-UN))-UN
  1620. DZEGAU(NP )= XXXX
  1621. DZEGAU(NP+1)= XXXX
  1622. DZEGAU(NP+2)= XXXX
  1623. 7772 CONTINUE
  1624. CALL SIMPSN(NNAPPE,3,1,POIGAU)
  1625. ELSE
  1626. GOTO 66
  1627. ENDIF
  1628. GOTO 666
  1629. C
  1630. 96 CONTINUE
  1631. C
  1632. C ELEMENT TUYO LONGUEUR = 2
  1633. C
  1634. IF(NBG.EQ.1) THEN
  1635. QSIGAU(1)= 0.D0
  1636. POIGAU(1)= 2.D0
  1637. ELSE
  1638. NTETA= NBG/2
  1639. DO 961 IN=1,NTETA
  1640. QSIGAU(IN)= -X577
  1641. QSIGAU(IN+NTETA)= X577
  1642. 961 CONTINUE
  1643. *
  1644. DTETA=2.D0*XPI/DBLE(NTETA)
  1645. TETA=0.D0
  1646. DO 962 IN=1,NTETA
  1647. ETAGAU(IN)=TETA
  1648. ETAGAU(IN+NTETA)=TETA
  1649. POIGAU(IN)=DTETA
  1650. POIGAU(IN+NTETA)=DTETA
  1651. TETA=TETA+DTETA
  1652. 962 CONTINUE
  1653. ENDIF
  1654. GO TO 666
  1655. *
  1656. 1570 CONTINUE
  1657. C
  1658. C=======================================================================
  1659. C ELEMENT HOMOGENE TRH6
  1660. C=======================================================================
  1661. C
  1662. IF(MELE.NE.157) GOTO 66
  1663. QSIGAU(1)=X044
  1664. ETAGAU(1)=X044
  1665. QSIGAU(2)=UN-(DEUX*X044)
  1666. ETAGAU(2)=X044
  1667. QSIGAU(3)=X044
  1668. ETAGAU(3)=UN-(DEUX*X044)
  1669. QSIGAU(4)=X009
  1670. ETAGAU(4)=X009
  1671. QSIGAU(5)=UN-(DEUX*X009)
  1672. ETAGAU(5)=X009
  1673. QSIGAU(6)=X009
  1674. ETAGAU(6)=UN-(DEUX*X009)
  1675. POIGAU(1)=X111
  1676. POIGAU(2)=X111
  1677. POIGAU(3)=X111
  1678. POIGAU(4)=X005
  1679. POIGAU(5)=X005
  1680. POIGAU(6)=X005
  1681. GOTO 666
  1682.  
  1683. *
  1684. 666 CONTINUE
  1685. IF(IRET.EQ.0) SEGSUP MINTE
  1686. IF(IRET.EQ.1) SEGDES MINTE
  1687. RETURN
  1688. END
  1689.  
  1690.  
  1691.  
  1692.  
  1693.  

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