Télécharger la plaquette Cast3M


fr en

Télécharger shape3.eso

Retour à la liste

Numérotation des lignes : 

shape3

  1. C SHAPE3    SOURCE    GOUNAND   21/06/02    21:17:39     11022          
  2.       SUBROUTINE SHAPE3(QSI,ETA,DZE,IELE,SHP,IRET)
  3. C=======================================================================
  4. C
  5. C     FONCTIONS DE FORME TRIDIMENSIONNELLES
  6. C        EBERSOLT
  7. C=======================================================================
  8.       IMPLICIT INTEGER(I-N)
  9.       IMPLICIT REAL*8(A-H,O-Z)
  10.       DIMENSION SHP(6,*)
  11.       DATA XZER/0.D0/,UNDEMI/.5D0/,UN/1.D0/,DEUX/2.D0/,QUATRE/4.D0/
  12.       DATA CINQ/5.D0/,HUIT/8.D0/,TROIS/3.D0/
  13.       DATA UNQUA/.25D0/
  14. C=======================================================================
  15. C     QSI ETA DZE COORDONNEES REDUITES DU POINT DE GAUSS
  16. C
  17. C     SHP FONCTIONS DE FORME SHP(1,INO) FONCTIONS DE FORME DU ND INO
  18. C                            SHP(2,INO) DERIVEES PR RPPRT A QSI
  19. C                            SHP(3,INO) DERIVEES PR RPPRT A ETA
  20. C                            SHP(4,INO) DERIVEES PR RPPRT A DZE
  21. C                            SHP(5,INO) DERIVEES PR RPPRT A  4
  22. C                              DERIVEES BARYCENTRIQUE
  23. C     CAS DE DIMENSION 3
  24. C               4 OU 10 NOEUDS  TETRAEDRE
  25. C               5 OU 13 NOEUDS  PYRAMIDE
  26. C               6 OU 15 NOEUDS  PRISME
  27. C               8 OU 20 NOEUDS  CUBE
  28. C=======================================================================
  29. C
  30.       IRET=1
  31. C           P1 S2 S3 T3 T4 T6 T7 Q4 Q5 Q8 Q9 R2 R3 C8 C20 P6 P15
  32.       GOTO (66,66,66,66,66,66,66,66,66,66,66,66,66,14, 15,16, 17,
  33. C        L3 L4 L6 L8 MU T4 T10 P5 P13
  34.      1   66,66,66,66,66,23, 24,25, 30) IELE
  35. C
  36.   66  IRET=0
  37.       GOTO 666
  38. C
  39.   16  CONTINUE
  40. C
  41. C     FONCTIONS DE FORME PRISMATIQUES  A 6 NOEUDS
  42. C
  43.       AUX=UN-QSI-ETA
  44.       AUX1=(UN-DZE)/DEUX
  45.       AUX2=(UN+DZE)/DEUX
  46.       SHP(1,1)=AUX*AUX1
  47.       SHP(1,2)=QSI*AUX1
  48.       SHP(1,3)=ETA*AUX1
  49.       SHP(2,1)=-AUX1
  50.       SHP(2,2)=AUX1
  51.       SHP(2,3)=XZER
  52.       SHP(3,1)=-AUX1
  53.       SHP(3,2)=XZER
  54.       SHP(3,3)=AUX1
  55.       SHP(4,1)=-AUX/DEUX
  56.       SHP(4,2)=-QSI/DEUX
  57.       SHP(4,3)=-ETA/DEUX
  58.       SHP(1,4)=AUX*AUX2
  59.       SHP(1,5)=QSI*AUX2
  60.       SHP(1,6)=ETA*AUX2
  61.       SHP(2,4)=-AUX2
  62.       SHP(2,5)=AUX2
  63.       SHP(2,6)=XZER
  64.       SHP(3,4)=-AUX2
  65.       SHP(3,5)=XZER
  66.       SHP(3,6)=AUX2
  67.       SHP(4,4)=AUX/DEUX
  68.       SHP(4,5)=QSI/DEUX
  69.       SHP(4,6)=ETA/DEUX
  70.       GOTO 666
  71.   17  CONTINUE
  72. C
  73. C     FONCTION DE FORME PRISMATIQUES A 15 NOEUDS
  74. C
  75.       AUX=UN-QSI-ETA
  76.       DZEM=(UN-DZE)
  77.       DZEP=(UN+DZE)
  78.       PAUX =QSI-UNDEMI*DZE-UN
  79.       PAUX1=ETA-UNDEMI*DZE-UN
  80.       PAUX2=QSI+ETA+UNDEMI*DZE
  81.       PAUX3=QSI+UNDEMI*DZE-UN
  82.       PAUX4=ETA+UNDEMI*DZE-UN
  83.       PAUX5=QSI+ETA-UNDEMI*DZE
  84.       SHP(1,1)=-AUX*DZEM*PAUX2
  85.       SHP(1,2)=DEUX*QSI*AUX*DZEM
  86.       SHP(1,3)=QSI*DZEM*PAUX
  87.       SHP(2,1)=(PAUX2-AUX)*DZEM
  88.       SHP(2,2)=DEUX*(AUX-QSI)*DZEM
  89.       SHP(2,3)=(PAUX+QSI)*DZEM
  90.       SHP(3,1)=(PAUX2-AUX)*DZEM
  91.       SHP(3,2)=-DEUX*QSI*DZEM
  92.       SHP(3,3)=XZER
  93.       SHP(4,1)=AUX*(PAUX2-UNDEMI*DZEM)
  94.       SHP(4,2)=-DEUX*QSI*AUX
  95.       SHP(4,3)=-QSI*(UNDEMI*DZEM+PAUX)
  96.       SHP(1,4)=DEUX*QSI*ETA*DZEM
  97.       SHP(1,5)=ETA*DZEM*PAUX1
  98.       SHP(1,6)=DEUX*ETA*AUX*DZEM
  99.       SHP(2,4)=DEUX*ETA*DZEM
  100.       SHP(2,5)=XZER
  101.       SHP(2,6)=-DEUX*ETA*DZEM
  102.       SHP(3,4)=DEUX*QSI*DZEM
  103.       SHP(3,5)=(ETA+PAUX1)*DZEM
  104.       SHP(3,6)=DEUX*(AUX-ETA)*DZEM
  105.       SHP(4,4)=-DEUX*QSI*ETA
  106.       SHP(4,5)=-ETA*(PAUX1+UNDEMI*DZEM)
  107.       SHP(4,6)=-DEUX*ETA*AUX
  108.       SHP(1,7)=AUX*DZEM*DZEP
  109.       SHP(1,8)=QSI*DZEM*DZEP
  110.       SHP(1,9)=ETA*DZEM*DZEP
  111.       SHP(2,7)=-DZEM*DZEP
  112.       SHP(2,8)=DZEM*DZEP
  113.       SHP(2,9)=XZER
  114.       SHP(3,7)=-DZEM*DZEP
  115.       SHP(3,8)=XZER
  116.       SHP(3,9)=DZEM*DZEP
  117.       SHP(4,7)=-DEUX*DZE*AUX
  118.       SHP(4,8)=-DEUX*DZE*QSI
  119.       SHP(4,9)=-DEUX*DZE*ETA
  120.       SHP(1,10)=-AUX*DZEP*PAUX5
  121.       SHP(1,11)=DEUX*QSI*AUX*DZEP
  122.       SHP(1,12)=QSI*DZEP*PAUX3
  123.       SHP(2,10)=(PAUX5-AUX)*DZEP
  124.       SHP(2,11)=DEUX*(AUX-QSI)*DZEP
  125.       SHP(2,12)=(PAUX3+QSI)*DZEP
  126.       SHP(3,10)=(PAUX5-AUX)*DZEP
  127.       SHP(3,11)=-DEUX*QSI*DZEP
  128.       SHP(3,12)=XZER
  129.       SHP(4,10)=AUX*(UNDEMI*DZEP-PAUX5)
  130.       SHP(4,11)=DEUX*QSI*AUX
  131.       SHP(4,12)=QSI*(PAUX3+UNDEMI*DZEP)
  132.       SHP(1,13)=DEUX*QSI*ETA*DZEP
  133.       SHP(1,14)=ETA*DZEP*PAUX4
  134.       SHP(1,15)=DEUX*ETA*AUX*DZEP
  135.       SHP(2,13)=DEUX*ETA*DZEP
  136.       SHP(2,14)=XZER
  137.       SHP(2,15)=-DEUX*ETA*DZEP
  138.       SHP(3,13)=DEUX*QSI*DZEP
  139.       SHP(3,14)=(ETA+PAUX4)*DZEP
  140.       SHP(3,15)=DEUX*(AUX-ETA)*DZEP
  141.       SHP(4,13)=DEUX*QSI*ETA
  142.       SHP(4,14)=ETA*(PAUX4+UNDEMI*DZEP)
  143.       SHP(4,15)=DEUX*ETA*AUX
  144.       GOTO 666
  145.   23  CONTINUE
  146. C
  147. C     FONCTIONS DE FORME TETRAEDRE A 4 NOEUDS
  148. C
  149.       AUX=UN-QSI-ETA-DZE
  150.       SHP(1,1)=AUX
  151.       SHP(1,2)=QSI
  152.       SHP(1,3)=ETA
  153.       SHP(1,4)=DZE
  154.       SHP(2,1)=-UN
  155.       SHP(2,2)= UN
  156.       SHP(2,3)= XZER
  157.       SHP(2,4)= XZER
  158.       SHP(3,1)=-UN
  159.       SHP(3,2)= XZER
  160.       SHP(3,3)= UN
  161.       SHP(3,4)= XZER
  162.       SHP(4,1)=-UN
  163.       SHP(4,2)= XZER
  164.       SHP(4,3)= XZER
  165.       SHP(4,4)= UN
  166.       GOTO 666
  167.   24  CONTINUE
  168. C
  169. C     FONCTIONS DE FORME TETRAEDRE A 10 NOEUDS
  170. C
  171.       AUX=UN-QSI-ETA-DZE
  172.       SHP(1,1)=AUX*(DEUX*AUX-UN)
  173.       SHP(1,2)=QUATRE*QSI*AUX
  174.       SHP(1,3)=QSI*(DEUX*QSI-UN)
  175.       SHP(1,4)=QUATRE*QSI*ETA
  176.       SHP(1,5)=ETA*(DEUX*ETA-UN)
  177.       SHP(1,6)=QUATRE*ETA*AUX
  178.       SHP(1,7)=QUATRE*DZE*AUX
  179.       SHP(1,8)=QUATRE*QSI*DZE
  180.       SHP(1,9)=QUATRE*ETA*DZE
  181.       SHP(1,10)=DZE*(DEUX*DZE-UN)
  182.       SHP(2,1)=UN-QUATRE*AUX
  183.       SHP(2,2)=QUATRE*(AUX-QSI)
  184.       SHP(2,3)=QUATRE*QSI-UN
  185.       SHP(2,4)=QUATRE*ETA
  186.       SHP(2,5)=XZER
  187.       SHP(2,6)=-QUATRE*ETA
  188.       SHP(2,7)=-QUATRE*DZE
  189.       SHP(2,8)= QUATRE*DZE
  190.       SHP(2,9)= XZER
  191.       SHP(2,10)=XZER
  192.       SHP(3,1)= UN-QUATRE*AUX
  193.       SHP(3,2)= -QUATRE*QSI
  194.       SHP(3,3)= XZER
  195.       SHP(3,4)= QUATRE*QSI
  196.       SHP(3,5)= QUATRE*ETA-UN
  197.       SHP(3,6)= QUATRE*(AUX-ETA)
  198.       SHP(3,7)=-QUATRE*DZE
  199.       SHP(3,8)= XZER
  200.       SHP(3,9)= QUATRE*DZE
  201.       SHP(3,10)=XZER
  202.       SHP(4,1)= UN-QUATRE*AUX
  203.       SHP(4,2)= -QUATRE*QSI
  204.       SHP(4,3)= XZER
  205.       SHP(4,4)= XZER
  206.       SHP(4,5)= XZER
  207.       SHP(4,6)= -QUATRE*ETA
  208.       SHP(4,7)=QUATRE*(AUX-DZE)
  209.       SHP(4,8)=QUATRE*QSI
  210.       SHP(4,9)=QUATRE*ETA
  211.       SHP(4,10)=QUATRE*DZE-UN
  212.       GOTO 666
  213.   14  CONTINUE
  214. C
  215. C     FONCTIONS DE FORME CUBE A 8 NOEUDS
  216. C
  217.       QSIM= UN-QSI
  218.       QSIP= UN+QSI
  219.       ETAM= UN-ETA
  220.       ETAP= UN+ETA
  221.       DZEM= UN-DZE
  222.       DZEP= UN+DZE
  223.       SHP(1,1)=QSIM*ETAM*DZEM/HUIT
  224.       SHP(1,2)=QSIP*ETAM*DZEM/HUIT
  225.       SHP(1,3)=QSIP*ETAP*DZEM/HUIT
  226.       SHP(1,4)=QSIM*ETAP*DZEM/HUIT
  227.       SHP(1,5)=QSIM*ETAM*DZEP/HUIT
  228.       SHP(1,6)=QSIP*ETAM*DZEP/HUIT
  229.       SHP(1,7)=QSIP*ETAP*DZEP/HUIT
  230.       SHP(1,8)=QSIM*ETAP*DZEP/HUIT
  231.       SHP(2,1)=    -ETAM*DZEM/HUIT
  232.       SHP(2,2)=     ETAM*DZEM/HUIT
  233.       SHP(2,3)=     ETAP*DZEM/HUIT
  234.       SHP(2,4)=    -ETAP*DZEM/HUIT
  235.       SHP(2,5)=    -ETAM*DZEP/HUIT
  236.       SHP(2,6)=     ETAM*DZEP/HUIT
  237.       SHP(2,7)=     ETAP*DZEP/HUIT
  238.       SHP(2,8)=    -ETAP*DZEP/HUIT
  239.       SHP(3,1)=-QSIM*DZEM/HUIT
  240.       SHP(3,2)=-QSIP*DZEM/HUIT
  241.       SHP(3,3)= QSIP*DZEM/HUIT
  242.       SHP(3,4)= QSIM*DZEM/HUIT
  243.       SHP(3,5)=-QSIM*DZEP/HUIT
  244.       SHP(3,6)=-QSIP*DZEP/HUIT
  245.       SHP(3,7)= QSIP*DZEP/HUIT
  246.       SHP(3,8)= QSIM*DZEP/HUIT
  247.       SHP(4,1)=-QSIM*ETAM/HUIT
  248.       SHP(4,2)=-QSIP*ETAM/HUIT
  249.       SHP(4,3)=-QSIP*ETAP/HUIT
  250.       SHP(4,4)=-QSIM*ETAP/HUIT
  251.       SHP(4,5)= QSIM*ETAM/HUIT
  252.       SHP(4,6)= QSIP*ETAM/HUIT
  253.       SHP(4,7)= QSIP*ETAP/HUIT
  254.       SHP(4,8)= QSIM*ETAP/HUIT
  255.       GOTO 666
  256.   15  CONTINUE
  257. C
  258. C     FONCTIONS DE FORME CUBE A 20 NOEUDS
  259. C
  260.       QSIM=(UN-QSI)
  261.       QSIP=(UN+QSI)
  262.       ETAM=(UN-ETA)
  263.       ETAP=(UN+ETA)
  264.       DZEM=(UN-DZE)
  265.       DZEP=(UN+DZE)
  266. C     NOEUDS AUX SOMMETS 1 3 5 7 13 15 17 19
  267.       SHP(1,1 )=QSIM*ETAM*DZEM*(QSIM+ETAM+DZEM-CINQ)/HUIT
  268.       SHP(1,3 )=QSIP*ETAM*DZEM*(QSIP+ETAM+DZEM-CINQ)/HUIT
  269.       SHP(1,5 )=QSIP*ETAP*DZEM*(QSIP+ETAP+DZEM-CINQ)/HUIT
  270.       SHP(1,7 )=QSIM*ETAP*DZEM*(QSIM+ETAP+DZEM-CINQ)/HUIT
  271.       SHP(1,13)=QSIM*ETAM*DZEP*(QSIM+ETAM+DZEP-CINQ)/HUIT
  272.       SHP(1,15)=QSIP*ETAM*DZEP*(QSIP+ETAM+DZEP-CINQ)/HUIT
  273.       SHP(1,17)=QSIP*ETAP*DZEP*(QSIP+ETAP+DZEP-CINQ)/HUIT
  274.       SHP(1,19)=QSIM*ETAP*DZEP*(QSIM+ETAP+DZEP-CINQ)/HUIT
  275.       SHP(2,1 )=    -ETAM*DZEM*(DEUX*QSIM+ETAM+DZEM-CINQ)/HUIT
  276.       SHP(2,3 )=     ETAM*DZEM*(DEUX*QSIP+ETAM+DZEM-CINQ)/HUIT
  277.       SHP(2,5 )=     ETAP*DZEM*(DEUX*QSIP+ETAP+DZEM-CINQ)/HUIT
  278.       SHP(2,7 )=    -ETAP*DZEM*(DEUX*QSIM+ETAP+DZEM-CINQ)/HUIT
  279.       SHP(2,13)=    -ETAM*DZEP*(DEUX*QSIM+ETAM+DZEP-CINQ)/HUIT
  280.       SHP(2,15)=     ETAM*DZEP*(DEUX*QSIP+ETAM+DZEP-CINQ)/HUIT
  281.       SHP(2,17)=     ETAP*DZEP*(DEUX*QSIP+ETAP+DZEP-CINQ)/HUIT
  282.       SHP(2,19)=    -ETAP*DZEP*(DEUX*QSIM+ETAP+DZEP-CINQ)/HUIT
  283.       SHP(3,1 )=    -QSIM*DZEM*(QSIM+DEUX*ETAM+DZEM-CINQ)/HUIT
  284.       SHP(3,3 )=    -QSIP*DZEM*(QSIP+DEUX*ETAM+DZEM-CINQ)/HUIT
  285.       SHP(3,5 )=     QSIP*DZEM*(QSIP+DEUX*ETAP+DZEM-CINQ)/HUIT
  286.       SHP(3,7 )=     QSIM*DZEM*(QSIM+DEUX*ETAP+DZEM-CINQ)/HUIT
  287.       SHP(3,13)=    -QSIM*DZEP*(QSIM+DEUX*ETAM+DZEP-CINQ)/HUIT
  288.       SHP(3,15)=    -QSIP*DZEP*(QSIP+DEUX*ETAM+DZEP-CINQ)/HUIT
  289.       SHP(3,17)=     QSIP*DZEP*(QSIP+DEUX*ETAP+DZEP-CINQ)/HUIT
  290.       SHP(3,19)=     QSIM*DZEP*(QSIM+DEUX*ETAP+DZEP-CINQ)/HUIT
  291.       SHP(4,1 )=    -QSIM*ETAM*(QSIM+ETAM+DEUX*DZEM-CINQ)/HUIT
  292.       SHP(4,3 )=    -QSIP*ETAM*(QSIP+ETAM+DEUX*DZEM-CINQ)/HUIT
  293.       SHP(4,5 )=    -QSIP*ETAP*(QSIP+ETAP+DEUX*DZEM-CINQ)/HUIT
  294.       SHP(4,7 )=    -QSIM*ETAP*(QSIM+ETAP+DEUX*DZEM-CINQ)/HUIT
  295.       SHP(4,13)=     QSIM*ETAM*(QSIM+ETAM+DEUX*DZEP-CINQ)/HUIT
  296.       SHP(4,15)=     QSIP*ETAM*(QSIP+ETAM+DEUX*DZEP-CINQ)/HUIT
  297.       SHP(4,17)=     QSIP*ETAP*(QSIP+ETAP+DEUX*DZEP-CINQ)/HUIT
  298.       SHP(4,19)=     QSIM*ETAP*(QSIM+ETAP+DEUX*DZEP-CINQ)/HUIT
  299. C     NOEUDS SUR LES COTES PARALLELES A L AXE QSI 2 6 14 18
  300.       SHP(1,2 )=UNQUA*(UN-QSI*QSI)*ETAM*DZEM
  301.       SHP(1,6 )=UNQUA*(UN-QSI*QSI)*ETAP*DZEM
  302.       SHP(1,14)=UNQUA*(UN-QSI*QSI)*ETAM*DZEP
  303.       SHP(1,18)=UNQUA*(UN-QSI*QSI)*ETAP*DZEP
  304.       SHP(2,2 )=-UNDEMI*QSI*ETAM*DZEM
  305.       SHP(2,6 )=-UNDEMI*QSI*ETAP*DZEM
  306.       SHP(2,14)=-UNDEMI*QSI*ETAM*DZEP
  307.       SHP(2,18)=-UNDEMI*QSI*ETAP*DZEP
  308.       SHP(3,2 )=-UNQUA*(UN-QSI*QSI)*DZEM
  309.       SHP(3,6 )= UNQUA*(UN-QSI*QSI)*DZEM
  310.       SHP(3,14)=-UNQUA*(UN-QSI*QSI)*DZEP
  311.       SHP(3,18)= UNQUA*(UN-QSI*QSI)*DZEP
  312.       SHP(4,2 )=-UNQUA*(UN-QSI*QSI)*ETAM
  313.       SHP(4,6 )=-UNQUA*(UN-QSI*QSI)*ETAP
  314.       SHP(4,14)= UNQUA*(UN-QSI*QSI)*ETAM
  315.       SHP(4,18)= UNQUA*(UN-QSI*QSI)*ETAP
  316. C     NOEUDS SUR LES COTES PARALELLES A L AXE ETA 4 8 16 20
  317.       SHP(1,4 )=UNQUA*QSIP*(UN-ETA*ETA)*DZEM
  318.       SHP(1,8 )=UNQUA*QSIM*(UN-ETA*ETA)*DZEM
  319.       SHP(1,16)=UNQUA*QSIP*(UN-ETA*ETA)*DZEP
  320.       SHP(1,20)=UNQUA*QSIM*(UN-ETA*ETA)*DZEP
  321.       SHP(2,4 )= UNQUA*(UN-ETA*ETA)*DZEM
  322.       SHP(2,8 )=-UNQUA*(UN-ETA*ETA)*DZEM
  323.       SHP(2,16)= UNQUA*(UN-ETA*ETA)*DZEP
  324.       SHP(2,20)=-UNQUA*(UN-ETA*ETA)*DZEP
  325.       SHP(3,4 )=-UNDEMI*QSIP*ETA*DZEM
  326.       SHP(3,8 )=-UNDEMI*QSIM*ETA*DZEM
  327.       SHP(3,16)=-UNDEMI*QSIP*ETA*DZEP
  328.       SHP(3,20)=-UNDEMI*QSIM*ETA*DZEP
  329.       SHP(4,4 )=-UNQUA*QSIP*(UN-ETA*ETA)
  330.       SHP(4,8 )=-UNQUA*QSIM*(UN-ETA*ETA)
  331.       SHP(4,16)= UNQUA*QSIP*(UN-ETA*ETA)
  332.       SHP(4,20)= UNQUA*QSIM*(UN-ETA*ETA)
  333. C     NOEUDS SUR LES COTES PARALELLES A L AXE DZE 9 10 11 12
  334.       SHP(1,9 )= UNQUA*QSIM*ETAM*(UN-DZE*DZE)
  335.       SHP(1,10)= UNQUA*QSIP*ETAM*(UN-DZE*DZE)
  336.       SHP(1,11)= UNQUA*QSIP*ETAP*(UN-DZE*DZE)
  337.       SHP(1,12)= UNQUA*QSIM*ETAP*(UN-DZE*DZE)
  338.       SHP(2,9 )=-UNQUA*ETAM*(UN-DZE*DZE)
  339.       SHP(2,10)= UNQUA*ETAM*(UN-DZE*DZE)
  340.       SHP(2,11)= UNQUA*ETAP*(UN-DZE*DZE)
  341.       SHP(2,12)=-UNQUA*ETAP*(UN-DZE*DZE)
  342.       SHP(3,9 )=-UNQUA*QSIM*(UN-DZE*DZE)
  343.       SHP(3,10)=-UNQUA*QSIP*(UN-DZE*DZE)
  344.       SHP(3,11)= UNQUA*QSIP*(UN-DZE*DZE)
  345.       SHP(3,12)= UNQUA*QSIM*(UN-DZE*DZE)
  346.       SHP(4,9 )=-UNDEMI*QSIM*ETAM*DZE
  347.       SHP(4,10)=-UNDEMI*QSIP*ETAM*DZE
  348.       SHP(4,11)=-UNDEMI*QSIP*ETAP*DZE
  349.       SHP(4,12)=-UNDEMI*QSIM*ETAP*DZE
  350.       GOTO 666
  351.   25  CONTINUE
  352. C
  353. C     FONCTIONS DE FORME PYRAMIDE A 5 NOEUDS
  354. C
  355.       DZEM =UN-DZE
  356.       DZEM4=(UN-DZE)*QUATRE
  357.       AUX = QSI+ETA+DZE-UN
  358.       AUX1=-QSI+ETA+DZE-UN
  359.       AUX2=-QSI-ETA+DZE-UN
  360.       AUX3= QSI-ETA+DZE-UN
  361.       IF(ABS(DZE-UN).LT.1.D-6) THEN
  362.       DO 251 IB=1,5
  363.          DO 2511 IA=1,3
  364.             SHP(IA,IB)= XZER
  365.  2511    CONTINUE
  366.  251  CONTINUE
  367.       SHP(1,5)= UN
  368.       SHP(2,1)= UNDEMI
  369.       SHP(2,3)=-UNDEMI
  370.       SHP(3,2)= UNDEMI
  371.       SHP(3,4)=-UNDEMI
  372.       DO 252 IA=1,4
  373.       SHP(4,IA)=-UNDEMI
  374.   252 CONTINUE
  375.       SHP(4,5)= UN
  376.       GO TO 666
  377.       ENDIF
  378.       SHP(1,1)=AUX1*AUX2/DZEM4
  379.       SHP(1,2)=AUX2*AUX3/DZEM4
  380.       SHP(1,3)=AUX *AUX3/DZEM4
  381.       SHP(1,4)=AUX *AUX1/DZEM4
  382.       SHP(1,5)=DZE
  383.       SHP(2,1)=(-AUX1-AUX2)/DZEM4
  384.       SHP(2,2)=( AUX2-AUX3)/DZEM4
  385.       SHP(2,3)=( AUX +AUX3)/DZEM4
  386.       SHP(2,4)=( AUX1-AUX )/DZEM4
  387.       SHP(2,5)=XZER
  388.       SHP(3,1)=( AUX2-AUX1)/DZEM4
  389.       SHP(3,2)=(-AUX2-AUX3)/DZEM4
  390.       SHP(3,3)=( AUX3-AUX )/DZEM4
  391.       SHP(3,4)=( AUX +AUX1)/DZEM4
  392.       SHP(3,5)=XZER
  393.       SHP(4,1)=(AUX1+AUX2+AUX1*AUX2/DZEM)/DZEM4
  394.       SHP(4,2)=(AUX2+AUX3+AUX2*AUX3/DZEM)/DZEM4
  395.       SHP(4,3)=(AUX3+AUX +AUX3*AUX /DZEM)/DZEM4
  396.       SHP(4,4)=( AUX+AUX1+ AUX*AUX1/DZEM)/DZEM4
  397.       SHP(4,5)=UN
  398.       GOTO 666
  399.   30  CONTINUE
  400. C
  401. C     FONCTIONS DE FORME PYRAMIDE A 13 NOEUDS
  402. C
  403.       DZEM=UN-DZE
  404.       DZEM2=(UN-DZE)*DEUX
  405.       PAUX = QSI+ETA+DZE-UN
  406.       PAUX1=-QSI+ETA+DZE-UN
  407.       PAUX2=-QSI-ETA+DZE-UN
  408.       PAUX3= QSI-ETA+DZE-UN
  409.       QAUX = QSI-UNDEMI
  410.       QAUX1= ETA-UNDEMI
  411.       QAUX2=-QSI-UNDEMI
  412.       QAUX3=-ETA-UNDEMI
  413.       IF(ABS(DZE-UN).LT.1.D-6) THEN
  414.       DO 301 IB=1,13
  415.          DO 3011 IA=1,3
  416.             SHP(IA,IB)= XZER
  417.  3011    CONTINUE
  418.   301 CONTINUE
  419.       SHP(1,13)= UN
  420.       SHP(2,1 )=-UNDEMI
  421.       SHP(2,5 )= UNDEMI
  422.       SHP(2,9 )= DEUX
  423.       SHP(2,11)=-DEUX
  424.       SHP(3,3 )=-UNDEMI
  425.       SHP(3,7 )= UNDEMI
  426.       SHP(3,10)= DEUX
  427.       SHP(3,12)=-DEUX
  428.       DO 302 IA=1,7,2
  429.       SHP(4,IA)=   UN/QUATRE
  430.       SHP(4,IA+1)= XZER
  431.   302 CONTINUE
  432.       DO 303 IA=9,12
  433.       SHP(4,IA)=-UN
  434.   303 CONTINUE
  435.       SHP(4,13)= TROIS
  436.       GO TO 666
  437.       ENDIF
  438.       SHP(1,1 )= PAUX1*PAUX2*QAUX /DZEM2
  439.       SHP(1,2 )=-PAUX1*PAUX2*PAUX3/DZEM2
  440.       SHP(1,3 )= PAUX2*PAUX3*QAUX1/DZEM2
  441.       SHP(1,4 )=-PAUX2*PAUX3*PAUX /DZEM2
  442.       SHP(1,5 )= PAUX3*PAUX *QAUX2/DZEM2
  443.       SHP(1,6 )=-PAUX3*PAUX *PAUX1/DZEM2
  444.       SHP(1,7 )= PAUX *PAUX1*QAUX3/DZEM2
  445.       SHP(1,8 )=-PAUX *PAUX1*PAUX2/DZEM2
  446.       SHP(1,9 )= DZE*PAUX1*PAUX2/DZEM
  447.       SHP(1,10)= DZE*PAUX2*PAUX3/DZEM
  448.       SHP(1,11)= DZE*PAUX3*PAUX /DZEM
  449.       SHP(1,12)= DZE*PAUX *PAUX1/DZEM
  450.       SHP(1,13)= DEUX*DZE*(DZE-UNDEMI)
  451. C
  452.       SHP(2,1 )=(PAUX1*PAUX2-(PAUX1+PAUX2)*QAUX)/DZEM2
  453.       SHP(2,2 )=(PAUX2*PAUX3+PAUX1*PAUX3-PAUX1*PAUX2)/DZEM2
  454.       SHP(2,3 )=(PAUX2-PAUX3)*QAUX1/DZEM2
  455.       SHP(2,4 )=(PAUX3*PAUX-PAUX2*PAUX-PAUX2*PAUX3)/DZEM2
  456.       SHP(2,5 )=((PAUX+PAUX3)*QAUX2-PAUX3*PAUX)/DZEM2
  457.       SHP(2,6 )=(PAUX3*PAUX-PAUX*PAUX1-PAUX3*PAUX1)/DZEM2
  458.       SHP(2,7 )=(PAUX1-PAUX)*QAUX3/DZEM2
  459.       SHP(2,8 )=(PAUX*PAUX2+PAUX*PAUX1-PAUX1*PAUX2)/DZEM2
  460.       SHP(2,9 )=(-PAUX2-PAUX1)*DZE/DZEM
  461.       SHP(2,10)=(PAUX2-PAUX3)*DZE/DZEM
  462.       SHP(2,11)=(PAUX+PAUX3)*DZE/DZEM
  463.       SHP(2,12)=(PAUX1-PAUX)*DZE/DZEM
  464.       SHP(2,13)=XZER
  465. C
  466.       SHP(3,1 )=(PAUX2-PAUX1)*QAUX/DZEM2
  467.       SHP(3,2 )=(PAUX1*PAUX3+PAUX1*PAUX2-PAUX2*PAUX3)/DZEM2
  468.       SHP(3,3 )=(PAUX2*PAUX3-(PAUX2+PAUX3)*QAUX1)/DZEM2
  469.       SHP(3,4 )=(PAUX3*PAUX+PAUX2*PAUX-PAUX2*PAUX3)/DZEM2
  470.       SHP(3,5 )=(PAUX3-PAUX)*QAUX2/DZEM2
  471.       SHP(3,6 )=(PAUX*PAUX1-PAUX3*PAUX1-PAUX3*PAUX)/DZEM2
  472.       SHP(3,7 )=((PAUX1+PAUX)*QAUX3-PAUX*PAUX1)/DZEM2
  473.       SHP(3,8 )=(PAUX*PAUX1-PAUX1*PAUX2-PAUX*PAUX2)/DZEM2
  474.       SHP(3,9 )=(PAUX2-PAUX1)*DZE/DZEM
  475.       SHP(3,10)=(-PAUX3-PAUX2)*DZE/DZEM
  476.       SHP(3,11)=(PAUX3-PAUX)*DZE/DZEM
  477.       SHP(3,12)=(PAUX+PAUX1)*DZE/DZEM
  478.       SHP(3,13)=XZER
  479. C
  480.       SHP(4,1 )=(PAUX1+PAUX2+PAUX1*PAUX2/DZEM)*QAUX/DZEM2
  481.       SHP(4,2 )=-UN*(PAUX2*PAUX3+PAUX1*PAUX3+PAUX1*PAUX2+
  482.      +           PAUX1*PAUX2*PAUX3/DZEM)/DZEM2
  483.       SHP(4,3 )=(PAUX2+PAUX3+PAUX2*PAUX3/DZEM)*QAUX1/DZEM2
  484.       SHP(4,4 )=-UN*(PAUX3*PAUX +PAUX2*PAUX+PAUX2*PAUX3+
  485.      +           PAUX2*PAUX3*PAUX/DZEM)/DZEM2
  486.       SHP(4,5 )=(PAUX+PAUX3+PAUX*PAUX3/DZEM)*QAUX2/DZEM2
  487.       SHP(4,6 )=-UN*(PAUX*PAUX1+PAUX3*PAUX1+PAUX3*PAUX+
  488.      +           PAUX3*PAUX*PAUX1/DZEM)/DZEM2
  489.       SHP(4,7 )=(PAUX1+PAUX+PAUX*PAUX1/DZEM)*QAUX3/DZEM2
  490.       SHP(4,8 )=-UN*(PAUX1*PAUX2+PAUX*PAUX2+PAUX*PAUX1+
  491.      +           PAUX*PAUX1*PAUX2/DZEM)/DZEM2
  492.       SHP(4,9 )=PAUX1*PAUX2/DZEM/DZEM+(PAUX2+PAUX1)*DZE/DZEM
  493.       SHP(4,10)=PAUX2*PAUX3/DZEM/DZEM+(PAUX3+PAUX2)*DZE/DZEM
  494.       SHP(4,11)=PAUX3*PAUX /DZEM/DZEM+(PAUX +PAUX3)*DZE/DZEM
  495.       SHP(4,12)=PAUX *PAUX1/DZEM/DZEM+(PAUX +PAUX1)*DZE/DZEM
  496.       SHP(4,13)=QUATRE*DZE-UN
  497.       GO TO 666
  498. C
  499.  666  CONTINUE
  500.       RETURN
  501.       END
  502.  
  503.  
  504.  

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