Télécharger la plaquette Cast3M


fr en

Télécharger donred.eso

Retour à la liste

Numérotation des lignes : 

donred

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

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