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

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