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

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