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

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