Télécharger la plaquette Cast3M


fr en

Télécharger shape.eso

Retour à la liste

Numérotation des lignes : 

shape

  1. C SHAPE     SOURCE    GOUNAND   21/06/02    21:17:38     11022          
  2.       SUBROUTINE SHAPE(QSI,ETA,DZE,IELE,SHP,IRET)
  3. C=======================================================================
  4. C
  5. C     FONCTIONS DE FORME ET LEUR DERIVEES
  6. C
  7. C=======================================================================
  8.       IMPLICIT INTEGER(I-N)
  9.       IMPLICIT REAL*8(A-H,O-Z)
  10.       PARAMETER(XZER=0.D0,UNDEMI=.5D0,UN=1.D0,DEUX=2.D0)
  11.       PARAMETER(TROIS=3.D0,QUATRE=4.D0)
  12.       DIMENSION SHP(6,*)
  13. C=======================================================================
  14. C     QSI ETA DZE COORDONNEES REDUITES DU POINT DE GAUSS
  15. C
  16. C     SHP FONCTIONS DE FORME SHP(1,INO) FONCTIONS DE FORME DU ND INO
  17. C                            SHP(2,INO) DERIVEES PR RPPRT A QSI
  18. C                            SHP(3,INO) DERIVEES PR RPPRT A ETA
  19. C                            SHP(4,INO) DERIVEES PR RPPRT A DZE
  20. C                            SHP(5,INO) DERIVEES PR RPPRT A 4
  21. C                               COORDONNEES BARYCENTRIQUES
  22. C     CAS DE DIMENSION 1
  23. C               2 OU 3 NOEUDS LIGNE
  24. C               2 NOEUDS POLYNOME DE HERMITE ET DERIVEES
  25. C     CAS DE DIMENSION 2
  26. C               3 OU 6 NOEUDS TRIANGLE
  27. C               7 NOEUDS TRIANGLE
  28. C               4 OU 8 NOEUDS QUADRILATERE
  29. C               9 NOEUDS QUADRILATERE
  30. C     CAS DE DIMENSION 3 FAIT DANS SHAPE3    ETIQUETTE 50
  31. C               4 OU 10 NOEUDS TETRAEDRE
  32. C               5 OU 13 NOEUDS PYRAMIDE
  33. C               6 OU 15 NOEUDS PRISME
  34. C               8 OU 20 NOEUDS CUBE
  35. C     CAS DE DIMENSION 3 QUAF FAIT DANS SHAPEQ    ETIQUETTE 51
  36. C               15 NOEUDS TETRAEDRE
  37. C               19 NOEUDS PYRAMIDE
  38. C               21 NOEUDS PRISME
  39. C               27 NOEUDS CUBE
  40. C=======================================================================
  41.       IRET=1
  42. C
  43. C    66 INDIQUE QUE L ELEMENT PAS ENCORE IMPLEMENTE
  44. C
  45. C           P1 S2 S3 T3 T4 T6 T7 Q4 Q5 Q8 Q9 R2 R3 C8 C20 P6 P15
  46.       GOTO ( 1, 2, 3, 4,66, 6,7, 8,66,10,10,12,13,50, 50,50, 50,
  47. C        L3 L4 L6 L8 MU T4 T10 P5 P13 ATT SUP RP3 LP6 LP8 POLY
  48.      1   18,19,66,66,66,50, 50,50, 50, 66, 66, 29, 30, 31, 66,
  49. C        CU27 PR21 TE15 PY19
  50.      1     51,51,51,51 ) IELE
  51.       if(iele.eq.260) then
  52.            write(6,*) ' on passe dans shape bizarre'
  53.       endif
  54. C
  55.   66  IRET=0
  56.       GOTO 666
  57. C
  58. C     FONCTIONS DE FORME MONODIMENSIONELLES
  59. C
  60. C
  61. C     FONCTION FORME 1 NOEUD
  62. C
  63.    1  CONTINUE
  64.       SHP(1,1)=UN
  65.       GOTO 666
  66. C     FONCTIONS DE FORME 2 NOEUDS
  67. C
  68.    2  CONTINUE
  69.       SHP(1,1)=(UN-QSI)/DEUX
  70.       SHP(1,2)=(UN+QSI)/DEUX
  71.       SHP(2,1)=-UNDEMI
  72.       SHP(2,2)= UNDEMI
  73.       GOTO 666
  74. C
  75. C     FONCTIONS DE FORME 3 NOEUDS
  76. C
  77.    3  CONTINUE
  78.       SHP(1,1)=(UN-QSI)/DEUX
  79.       SHP(1,3)=(UN+QSI)/DEUX
  80.       SHP(2,1)=-UNDEMI
  81.       SHP(2,3)= UNDEMI
  82.       SHP(1,2)=UN-QSI*QSI
  83.       SHP(2,2)=-DEUX*QSI
  84.       SHP(1,1)=SHP(1,1)-SHP(1,2)/DEUX
  85.       SHP(1,3)=SHP(1,3)-SHP(1,2)/DEUX
  86.       SHP(2,1)=SHP(2,1)-SHP(2,2)/DEUX
  87.       SHP(2,3)=SHP(2,3)-SHP(2,2)/DEUX
  88.       GOTO 666
  89. C
  90. C     POLYNOME DE HERMITE 2 NOEUDS 1 FCTN DE FORME ET SA DERIVEE
  91. C
  92. C104  SHP(1,1)=(UN-QSI)*(UN-QSI)*(DEUX+QSI)/QUATRE
  93. C     SHP(2,1)= TROIS*(QSI*QSI-UN)/QUATRE
  94. C     SHP(1,2)=(UN-QSI*QSI)*(UN-QSI)/QUATRE
  95. C     SHP(2,2)=(QSI-UN)*(UN+TROIS*QSI)/QUATRE
  96. C     SHP(1,3)=(UN+QSI)*(UN+QSI)*(DEUX-QSI)/QUATRE
  97. C     SHP(2,3)=TROIS*(UN-QSI*QSI)/QUATRE
  98. C     SHP(1,4)=(QSI*QSI-UN)*(UN+QSI)/QUATRE
  99. C     SHP(2,4)=(UN+QSI)*(TROIS*QSI-UN)/QUATRE
  100. C     GOTO 666
  101. C
  102. C     FONCTIONS DE FORME BIDIMENSIONNELLES
  103. C
  104.    4  CONTINUE
  105. C
  106. C     FONCTIONS DE FORME TRIANGLE A 3 NOEUDS
  107. C
  108.       AUX=UN-QSI-ETA
  109.       SHP(1,1)=AUX
  110.       SHP(1,2)=QSI
  111.       SHP(1,3)=ETA
  112.       SHP(2,1)=-UN
  113.       SHP(2,2)=UN
  114.       SHP(2,3)=XZER
  115.       SHP(3,1)=-UN
  116.       SHP(3,2)=XZER
  117.       SHP(3,3)=UN
  118.       GOTO 666
  119.    6  CONTINUE
  120. C
  121. C     FONCTIONS DE FORME TRIANGLE A 6 NOEUDS
  122. C
  123.       AUX=UN-QSI-ETA
  124.       SHP(1,1)=AUX
  125.       SHP(1,3)=QSI
  126.       SHP(1,5)=ETA
  127.       SHP(2,1)=-UN
  128.       SHP(2,3)=UN
  129.       SHP(2,5)=XZER
  130.       SHP(3,1)=-UN
  131.       SHP(3,3)=XZER
  132.       SHP(3,5)=UN
  133. C   AJOUTER NOEUD 2
  134.       SHP(1,2)=QUATRE*QSI*AUX
  135.       SHP(2,2)=QUATRE*(AUX-QSI)
  136.       SHP(3,2)=-QUATRE*QSI
  137.       DO 1010 I=1,3
  138.          SHP(I,1)=SHP(I,1)-SHP(I,2)/DEUX
  139.          SHP(I,3)=SHP(I,3)-SHP(I,2)/DEUX
  140.  1010 CONTINUE
  141. C   AJOUTER NOEUD 4
  142.       SHP(1,4)=QUATRE*QSI*ETA
  143.       SHP(2,4)=QUATRE*ETA
  144.       SHP(3,4)=QUATRE*QSI
  145.       DO 1020 I=1,3
  146.          SHP(I,3)=SHP(I,3)-SHP(I,4)/DEUX
  147.          SHP(I,5)=SHP(I,5)-SHP(I,4)/DEUX
  148.  1020 CONTINUE
  149. C   AJOUTER NOEUD 6
  150.       SHP(1,6)=QUATRE*ETA*AUX
  151.       SHP(2,6)=-QUATRE*ETA
  152.       SHP(3,6)=QUATRE*(AUX-ETA)
  153.       DO 1030 I=1,3
  154.          SHP(I,5)=SHP(I,5)-SHP(I,6)/DEUX
  155.          SHP(I,1)=SHP(I,1)-SHP(I,6)/DEUX
  156.  1030 CONTINUE
  157.       GOTO 666
  158.    7  CONTINUE
  159. C
  160. C     FONCTIONS DE FORME TRIANGLE A 7 NOEUDS
  161. C
  162.       AUX=UN-QSI-ETA
  163.       SHP(1,1)=AUX*(DEUX*AUX-UN)+TROIS*AUX*QSI*ETA
  164.       SHP(1,3)=QSI*(DEUX*QSI-UN)+TROIS*AUX*QSI*ETA
  165.       SHP(1,5)=ETA*(DEUX*ETA-UN)+TROIS*AUX*QSI*ETA
  166.  
  167.       SHP(1,2)=QUATRE*(AUX*QSI-TROIS*AUX*QSI*ETA)
  168.       SHP(1,4)=QUATRE*(QSI*ETA-TROIS*AUX*QSI*ETA)
  169.       SHP(1,6)=QUATRE*(AUX*ETA-TROIS*AUX*QSI*ETA)
  170.  
  171.       SHP(1,7)=TROIS*TROIS*TROIS*AUX*QSI*ETA
  172.  
  173. C     GRADIENTS DES FONCTIONS DE FORME A 7 NOEUDS
  174. C     NOEUD 1
  175.       SHP(2,1)=-3.D0*ETA*ETA+4.D0*QSI+7.D0*ETA-
  176.      * 6.D0*QSI*ETA-3.D0
  177.       SHP(3,1)=-3.D0*QSI*QSI+4.D0*ETA+7.D0*QSI-
  178.      * 6.D0*QSI*ETA-3.D0
  179. C     NOEUD 3
  180.       SHP(2,3)=-3.D0*ETA*ETA+4.D0*QSI+3.D0*ETA-
  181.      * 6.D0*QSI*ETA-1.D0
  182.       SHP(3,3)=-3.D0*QSI*QSI+3.D0*QSI-6.D0*QSI*ETA
  183. C     NOEUD 5
  184.       SHP(2,5)=-3.D0*ETA*ETA+3.D0*ETA-6.D0*QSI*ETA
  185.       SHP(3,5)=-3.D0*QSI*QSI+4.D0*ETA+3.D0*QSI-
  186.      * 6.D0*QSI*ETA-1.D0
  187. C     NOEUD 2
  188.       SHP(2,2)=QUATRE*TROIS*ETA*ETA-QUATRE*QUATRE*ETA-DEUX*QUATRE*QSI+
  189.      * DEUX*TROIS*QUATRE*QSI*ETA+QUATRE
  190.       SHP(3,2)=QUATRE*TROIS*QSI*QSI-QUATRE*QUATRE*QSI+
  191.      * DEUX*TROIS*QUATRE*QSI*ETA
  192. C     NOEUD 4
  193.       SHP(2,4)=QUATRE*TROIS*ETA*ETA-DEUX*QUATRE*ETA+
  194.      * DEUX*TROIS*QUATRE*QSI*ETA
  195.       SHP(3,4)=QUATRE*TROIS*QSI*QSI-DEUX*QUATRE*QSI+
  196.      * DEUX*TROIS*QUATRE*QSI*ETA
  197. C     NOEUD 6
  198.       SHP(2,6)=QUATRE*TROIS*ETA*ETA-QUATRE*QUATRE*ETA+
  199.      * DEUX*TROIS*QUATRE*QSI*ETA
  200.       SHP(3,6)=QUATRE*TROIS*QSI*QSI-QUATRE*QUATRE*QSI-DEUX*QUATRE*ETA+
  201.      * DEUX*TROIS*QUATRE*QSI*ETA+QUATRE
  202. C     NOEUD 7
  203.       SHP(2,7)=TROIS*TROIS*TROIS*ETA*(UN-DEUX*QSI-ETA)
  204.       SHP(3,7)=TROIS*TROIS*TROIS*QSI*(UN-QSI-DEUX*ETA)
  205.       GOTO 666
  206.    8  CONTINUE
  207. C
  208. C  FONCTIONS DE FORME QUADRANGLE  A 4 NOEUDS
  209. C
  210.       SHP(1,1)=(UN-QSI)*(UN-ETA)/QUATRE
  211.       SHP(1,2)=(UN+QSI)*(UN-ETA)/QUATRE
  212.       SHP(1,3)=(UN+QSI)*(UN+ETA)/QUATRE
  213.       SHP(1,4)=(UN-QSI)*(UN+ETA)/QUATRE
  214.       SHP(2,1)=(ETA-UN)/QUATRE
  215.       SHP(2,2)=-SHP(2,1)
  216.       SHP(2,3)=(ETA+UN)/QUATRE
  217.       SHP(2,4)=-SHP(2,3)
  218.       SHP(3,1)=(QSI-UN)/QUATRE
  219.       SHP(3,2)=-(UN+QSI)/QUATRE
  220.       SHP(3,3)=-SHP(3,2)
  221.       SHP(3,4)=-SHP(3,1)
  222.       GOTO 666
  223. C
  224. C     FONCTIONS DE FORME QUADRANGLE A 8 NOEUDS
  225. C
  226.   10  CONTINUE
  227.       SHP(1,1)=(UN-QSI)*(UN-ETA)/QUATRE
  228.       SHP(1,3)=(UN+QSI)*(UN-ETA)/QUATRE
  229.       SHP(1,5)=(UN+QSI)*(UN+ETA)/QUATRE
  230.       SHP(1,7)=(UN-QSI)*(UN+ETA)/QUATRE
  231.       SHP(2,1)=(ETA-UN)/QUATRE
  232.       SHP(2,3)=-SHP(2,1)
  233.       SHP(2,5)=(ETA+UN)/QUATRE
  234.       SHP(2,7)=-SHP(2,5)
  235.       SHP(3,1)=(QSI-UN)/QUATRE
  236.       SHP(3,3)=-(UN+QSI)/QUATRE
  237.       SHP(3,5)=-SHP(3,3)
  238.       SHP(3,7)=-SHP(3,1)
  239. C   AJOUTER NOEUD 2
  240.       SHP(1,2)=(UN-QSI*QSI)*(UN-ETA)/DEUX
  241.       SHP(2,2)=(ETA-UN)*QSI
  242.       SHP(3,2)=(QSI*QSI-1)/DEUX
  243.       DO 2010 I=1,3
  244.          SHP(I,1)=SHP(I,1)-SHP(I,2)/DEUX
  245.          SHP(I,3)=SHP(I,3)-SHP(I,2)/DEUX
  246.  2010 CONTINUE
  247. C   AJOUTER NOEUD 4
  248.       SHP(1,4)=(UN-ETA*ETA)*(UN+QSI)/DEUX
  249.       SHP(2,4)=(UN-ETA*ETA)/DEUX
  250.       SHP(3,4)=-ETA*(UN+QSI)
  251.       DO 2020 I=1,3
  252.          SHP(I,3)=SHP(I,3)-SHP(I,4)/DEUX
  253.          SHP(I,5)=SHP(I,5)-SHP(I,4)/DEUX
  254.  2020 CONTINUE
  255. C   AJOUTER NOEUD 6
  256.       SHP(1,6)=(UN-QSI*QSI)*(UN+ETA)/DEUX
  257.       SHP(2,6)=-QSI*(UN+ETA)
  258.       SHP(3,6)=(UN-QSI*QSI)/DEUX
  259.       DO 2030 I=1,3
  260.          SHP(I,5)=SHP(I,5)-SHP(I,6)/DEUX
  261.          SHP(I,7)=SHP(I,7)-SHP(I,6)/DEUX
  262.  2030 CONTINUE
  263. C   AJOUTER NOEUD 8
  264.       SHP(1,8)=(UN-ETA*ETA)*(UN-QSI)/DEUX
  265.       SHP(2,8)=(ETA*ETA-UN)/DEUX
  266.       SHP(3,8)=ETA*(QSI-UN)
  267.       DO 2040 I=1,3
  268.          SHP(I,7)=SHP(I,7)-SHP(I,8)/DEUX
  269.          SHP(I,1)=SHP(I,1)-SHP(I,8)/DEUX
  270.  2040 CONTINUE
  271. C
  272. C     FONCTIONS DE FORME QUADRANGLE A 9 NOEUDS
  273. C
  274.       IF (IELE.NE.11)  GOTO 666
  275.       SHP(1,9)=(UN-QSI*QSI)*(UN-ETA*ETA)
  276.       SHP(2,9)=DEUX*ETA*(QSI*QSI-UN)
  277.       SHP(3,9)=DEUX*QSI*(ETA*ETA-UN)
  278.       DO 2051 I=1,3
  279.          DO 20511 J=1,7,2
  280.             SHP(I,J)=SHP(I,J)+SHP(I,9)/QUATRE
  281.             SHP(I,J+1)=SHP(I,J+1)-SHP(I,9)/DEUX
  282. 20511    CONTINUE
  283.  2051 CONTINUE
  284.       GOTO 666
  285. C
  286. C     FONCTIONS DE FORME ELEMENTS DE RACCORD
  287. C            LIQUIDE - MASSIF LINEAIRE
  288. C
  289. C     FONCTIONS DE FORME RAC2
  290. C
  291.   12  CONTINUE
  292.       SHP(1,1)=(1.0-QSI)/DEUX
  293.       SHP(1,2)=(1.0+QSI)/DEUX
  294.       SHP(2,1)=-1.0/DEUX
  295.       SHP(2,2)= 1.0/DEUX
  296.       SHP(1,3)=SHP(1,2)
  297.       SHP(1,4)=SHP(1,1)
  298.       SHP(2,3)=SHP(2,2)
  299.       SHP(2,4)=SHP(2,1)
  300.       GOTO 666
  301. C
  302. C     FONCTIONS DE FORME RAC3
  303. C
  304.   13  CONTINUE
  305.       SHP(1,1)=(-QSI+QSI*QSI)*(UN-ETA)/QUATRE
  306.       SHP(1,2)=( UN -QSI*QSI)*(UN-ETA)/DEUX
  307.       SHP(1,3)=( QSI+QSI*QSI)*(UN-ETA)/QUATRE
  308.       SHP(1,4)=( QSI+QSI*QSI)*(UN+ETA)/QUATRE
  309.       SHP(1,5)=( UN -QSI*QSI)*(UN+ETA)/DEUX
  310.       SHP(1,6)=(-QSI+QSI*QSI)*(UN+ETA)/QUATRE
  311. C
  312.       SHP(2,1)=(-UN+DEUX*QSI)*(UN-ETA)/QUATRE
  313.       SHP(2,2)=(   -DEUX*QSI)*(UN-ETA)/DEUX
  314.       SHP(2,3)=( UN+DEUX*QSI)*(UN-ETA)/QUATRE
  315.       SHP(2,4)=( UN+DEUX*QSI)*(UN+ETA)/QUATRE
  316.       SHP(2,5)=(   -DEUX*QSI)*(UN+ETA)/DEUX
  317.       SHP(2,6)=(-UN+DEUX*QSI)*(UN+ETA)/QUATRE
  318. C
  319.       SHP(3,1)=-(-QSI+QSI*QSI)/QUATRE
  320.       SHP(3,2)=-( UN -QSI*QSI)/DEUX
  321.       SHP(3,3)=-( QSI+QSI*QSI)/QUATRE
  322.       SHP(3,4)= ( QSI+QSI*QSI)/QUATRE
  323.       SHP(3,5)= ( UN -QSI*QSI)/DEUX
  324.       SHP(3,6)= (-QSI+QSI*QSI)/QUATRE
  325.       GOTO 666
  326.   18  CONTINUE
  327. C
  328. C     FONCTIONS DE FORME LIA3
  329. C
  330.       AUX=UN-QSI-ETA
  331.       SHP(1,1)=AUX
  332.       SHP(1,2)=QSI
  333.       SHP(1,3)=ETA
  334.       SHP(2,1)=-UN
  335.       SHP(2,2)=UN
  336.       SHP(2,3)=XZER
  337.       SHP(3,1)=-UN
  338.       SHP(3,2)=XZER
  339.       SHP(3,3)=UN
  340.       SHP(1,4)=SHP(1,1)
  341.       SHP(1,5)=SHP(1,2)
  342.       SHP(1,6)=SHP(1,3)
  343.       SHP(2,4)=SHP(2,1)
  344.       SHP(2,5)=SHP(2,2)
  345.       SHP(2,6)=SHP(2,3)
  346.       SHP(3,4)=SHP(3,1)
  347.       SHP(3,5)=SHP(3,2)
  348.       SHP(3,6)=SHP(3,3)
  349.       GOTO 666
  350.   19  CONTINUE
  351. C
  352. C  FONCTIONS DE FORME LIA4
  353. C
  354.       SHP(1,1)=(UN-QSI)*(UN-ETA)/QUATRE
  355.       SHP(1,2)=(UN+QSI)*(UN-ETA)/QUATRE
  356.       SHP(1,3)=(UN+QSI)*(UN+ETA)/QUATRE
  357.       SHP(1,4)=(UN-QSI)*(UN+ETA)/QUATRE
  358.       SHP(2,1)=(ETA-UN)/QUATRE
  359.       SHP(2,2)=-SHP(2,1)
  360.       SHP(2,3)=(ETA+UN)/QUATRE
  361.       SHP(2,4)=-SHP(2,3)
  362.       SHP(3,1)=(QSI-UN)/QUATRE
  363.       SHP(3,2)=-(UN+QSI)/QUATRE
  364.       SHP(3,3)=-SHP(3,2)
  365.       SHP(3,4)=-SHP(3,1)
  366.       SHP(1,5)=SHP(1,1)
  367.       SHP(1,6)=SHP(1,2)
  368.       SHP(1,7)=SHP(1,3)
  369.       SHP(1,8)=SHP(1,4)
  370.       SHP(2,5)=SHP(2,1)
  371.       SHP(2,6)=SHP(2,2)
  372.       SHP(2,7)=SHP(2,3)
  373.       SHP(2,8)=SHP(2,4)
  374.       SHP(3,5)=SHP(3,1)
  375.       SHP(3,6)=SHP(3,2)
  376.       SHP(3,7)=SHP(3,3)
  377.       SHP(3,8)=SHP(3,4)
  378.       GOTO 666
  379.   29  CONTINUE
  380. C
  381. C     FONCTIONS DE FORME ELEMENTS DE RACCORD POREUX RAP3
  382. C
  383.       SHP(1,1)=(-QSI+QSI*QSI)/DEUX
  384.       SHP(1,2)=( UN -QSI*QSI)
  385.       SHP(1,3)=( QSI+QSI*QSI)/DEUX
  386.       SHP(1,4)=( QSI+QSI*QSI)/DEUX
  387.       SHP(1,5)=( UN -QSI*QSI)
  388.       SHP(1,6)=(-QSI+QSI*QSI)/DEUX
  389.       SHP(2,1)=(-UN+DEUX*QSI)/DEUX
  390.       SHP(2,2)=(   -DEUX*QSI)
  391.       SHP(2,3)=( UN+DEUX*QSI)/DEUX
  392.       SHP(2,4)=( UN+DEUX*QSI)/DEUX
  393.       SHP(2,5)=(   -DEUX*QSI)
  394.       SHP(2,6)=(-UN+DEUX*QSI)/DEUX
  395. c
  396.       SHP(1,7)= ( UN-QSI)/DEUX
  397.       SHP(1,8)= ( UN+QSI)/DEUX
  398.       SHP(1,9)= ( UN+QSI)/DEUX
  399.       SHP(1,10)=( UN-QSI)/DEUX
  400.       SHP(1,11)=( UN-QSI)/DEUX
  401.       SHP(1,12)=( UN+QSI)/DEUX
  402.       SHP(2,7)= -UN/DEUX
  403.       SHP(2,8)=  UN/DEUX
  404.       SHP(2,9)=  UN/DEUX
  405.       SHP(2,10)=-UN/DEUX
  406.       SHP(2,11)=-UN/DEUX
  407.       SHP(2,12)= UN/DEUX
  408.       GOTO 666
  409.   30  CONTINUE
  410. C
  411. C     FONCTIONS DE FORME ELEMENTS DE LIAISON POREUX RAP6
  412. C
  413.       SHP(1,1)= (UN-QSI-ETA)*(UN-DEUX*QSI-DEUX*ETA)
  414.       SHP(1,2)= QUATRE*QSI*(UN-QSI-ETA)
  415.       SHP(1,3)= QSI*(DEUX*QSI-UN)
  416.       SHP(1,4)= QUATRE*QSI*ETA
  417.       SHP(1,5)= ETA*(DEUX*ETA-UN)
  418.       SHP(1,6)= QUATRE*ETA*(UN-QSI-ETA)
  419.       SHP(2,1)=-(UN-DEUX*QSI-DEUX*ETA)-DEUX*(UN-QSI-ETA)
  420.       SHP(2,2)= QUATRE*(UN-DEUX*QSI-ETA)
  421.       SHP(2,3)= QUATRE*QSI-UN
  422.       SHP(2,4)= QUATRE*ETA
  423.       SHP(2,5)= XZER
  424.       SHP(2,6)=-QUATRE*ETA
  425.       SHP(3,1)=-(UN-DEUX*QSI-DEUX*ETA)-DEUX*(UN-QSI-ETA)
  426.       SHP(3,2)=-QUATRE*QSI
  427.       SHP(3,3)= XZER
  428.       SHP(3,4)= QUATRE*QSI
  429.       SHP(3,5)= QUATRE*ETA-UN
  430.       SHP(3,6)= QUATRE*(UN-QSI-DEUX*ETA)
  431.       DO 111 I=7,12
  432.       SHP(1,I)=SHP(1,I-6)
  433.       SHP(2,I)=SHP(2,I-6)
  434.       SHP(3,I)=SHP(3,I-6)
  435.   111 CONTINUE
  436. c
  437.       SHP(1,13)= UN-QSI-ETA
  438.       SHP(1,14)= QSI
  439.       SHP(1,15)= ETA
  440.       SHP(2,13)=-UN
  441.       SHP(2,14)= UN
  442.       SHP(2,15)= XZER
  443.       SHP(3,13)=-UN
  444.       SHP(3,14)= XZER
  445.       SHP(3,15)= UN
  446.       DO 222 I=16,18
  447.       SHP(1,I)=SHP(1,I-3)
  448.       SHP(2,I)=SHP(2,I-3)
  449.       SHP(3,I)=SHP(3,I-3)
  450.       SHP(1,I+3)=SHP(1,I-3)
  451.       SHP(2,I+3)=SHP(2,I-3)
  452.       SHP(3,I+3)=SHP(3,I-3)
  453.   222 CONTINUE
  454.       GOTO 666
  455.   31  CONTINUE
  456. C
  457. C     FONCTIONS DE FORME ELEMENTS DE LIAISON POREUX RAP8
  458. C
  459.       SHP(1,2)= (UN-QSI*QSI)*(UN-ETA)/DEUX
  460.       SHP(1,4)= (UN+QSI)*(UN-ETA*ETA)/DEUX
  461.       SHP(1,6)= (UN-QSI*QSI)*(UN+ETA)/DEUX
  462.       SHP(1,8)= (UN-QSI)*(UN-ETA*ETA)/DEUX
  463.       SHP(1,1)= (UN-QSI)*(UN-ETA)/QUATRE-UN/DEUX*(SHP(1,8)+SHP(1,2))
  464.       SHP(1,3)= (UN+QSI)*(UN-ETA)/QUATRE-UN/DEUX*(SHP(1,2)+SHP(1,4))
  465.       SHP(1,5)= (UN+QSI)*(UN+ETA)/QUATRE-UN/DEUX*(SHP(1,4)+SHP(1,6))
  466.       SHP(1,7)= (UN-QSI)*(UN+ETA)/QUATRE-UN/DEUX*(SHP(1,6)+SHP(1,8))
  467.       SHP(2,2)=-QSI*(UN-ETA)
  468.       SHP(2,4)= (UN-ETA*ETA)/DEUX
  469.       SHP(2,6)=-QSI*(UN+ETA)
  470.       SHP(2,8)=-(UN-ETA*ETA)/DEUX
  471.       SHP(2,1)=-(UN-ETA)/QUATRE-UN/DEUX*(SHP(2,8)+SHP(2,2))
  472.       SHP(2,3)= (UN-ETA)/QUATRE-UN/DEUX*(SHP(2,2)+SHP(2,4))
  473.       SHP(2,5)= (UN+ETA)/QUATRE-UN/DEUX*(SHP(2,4)+SHP(2,6))
  474.       SHP(2,7)=-(UN+ETA)/QUATRE-UN/DEUX*(SHP(2,6)+SHP(2,8))
  475.       SHP(3,2)=-(UN-QSI*QSI)/DEUX
  476.       SHP(3,4)=-ETA*(UN+QSI)
  477.       SHP(3,6)= (UN-QSI*QSI)/DEUX
  478.       SHP(3,8)=-ETA*(UN-QSI)
  479.       SHP(3,1)=-(UN-QSI)/QUATRE-UN/DEUX*(SHP(3,8)+SHP(3,2))
  480.       SHP(3,3)=-(UN+QSI)/QUATRE-UN/DEUX*(SHP(3,2)+SHP(3,4))
  481.       SHP(3,5)= (UN+QSI)/QUATRE-UN/DEUX*(SHP(3,4)+SHP(3,6))
  482.       SHP(3,7)= (UN-QSI)/QUATRE-UN/DEUX*(SHP(3,6)+SHP(3,8))
  483.       DO 333 I=9,16
  484.       SHP(1,I)=SHP(1,I-8)
  485.       SHP(2,I)=SHP(2,I-8)
  486.       SHP(3,I)=SHP(3,I-8)
  487.   333 CONTINUE
  488. c
  489.       SHP(1,17)= (UN-QSI)*(UN-ETA)/QUATRE
  490.       SHP(1,18)= (UN+QSI)*(UN-ETA)/QUATRE
  491.       SHP(1,19)= (UN+QSI)*(UN+ETA)/QUATRE
  492.       SHP(1,20)= (UN-QSI)*(UN+ETA)/QUATRE
  493.       SHP(2,17)=-(UN-ETA)/QUATRE
  494.       SHP(2,18)= (UN-ETA)/QUATRE
  495.       SHP(2,19)= (UN+ETA)/QUATRE
  496.       SHP(2,20)=-(UN+ETA)/QUATRE
  497.       SHP(3,17)=-(UN-QSI)/QUATRE
  498.       SHP(3,18)=-(UN+QSI)/QUATRE
  499.       SHP(3,19)= (UN+QSI)/QUATRE
  500.       SHP(3,20)= (UN-QSI)/QUATRE
  501.       DO 444 I=21,24
  502.       SHP(1,I)=SHP(1,I-4)
  503.       SHP(2,I)=SHP(2,I-4)
  504.       SHP(3,I)=SHP(3,I-4)
  505.       SHP(1,I+4)=SHP(1,I-4)
  506.       SHP(2,I+4)=SHP(2,I-4)
  507.       SHP(3,I+4)=SHP(3,I-4)
  508.   444 CONTINUE
  509.       GOTO 666
  510. C
  511. C     FONCTIONS DE FORME TRIDIMENSIONNELLES
  512. C
  513.   50  CALL SHAPE3(QSI,ETA,DZE,IELE,SHP,IRT1)
  514.       IRET=IRT1
  515.       GOTO 666
  516. C
  517. C     FONCTIONS DE FORME TRIDIMENSIONNELLES QUAF
  518. C
  519.  51   CALL SHAPEQ(QSI,ETA,DZE,IELE,SHP,IRT1)
  520.       IRET=IRT1
  521. C
  522.  666  CONTINUE
  523.       RETURN
  524.       END
  525.  
  526.  
  527.  
  528.  
  529.  
  530.  
  531.  
  532.  
  533.  
  534.  
  535.  

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