Télécharger la plaquette Cast3M


fr en

Télécharger tpr15.eso

Retour à la liste

Numérotation des lignes : 

tpr15

  1. C TPR15     SOURCE    OF166741  23/12/04    21:15:14     11800          
  2.  
  3. C=======================================================================
  4. C=                            T P R 1 5                                =
  5. C=                            ---------                                =
  6. C=                                                                     =
  7. C=  Fonction :                                                         =
  8. C=  ----------                                                         =
  9. C=   Acquisition des caracteristiques d'integration aux noeuds, aux    =
  10. C=   points d'integration (Gauss) ou au centre de gravite pour         =
  11. C=   l'element thermique support PR15.                                 =
  12. C=                                                                     =
  13. C=  Parametres :  (E)=Entree  (S)=Sortie                               =
  14. C=  ------------                                                       =
  15. C=   ISUPTH  (E)   Indique le support d'integration demande            =
  16. C=   IPINTE  (S)   Pointeur sur le segment MINTE associe (ACTIF en S)  =
  17. C=======================================================================
  18.  
  19.       SUBROUTINE TPR15 (ISUPTH,IPINTE)
  20.  
  21.       IMPLICIT INTEGER(I-N)
  22.       IMPLICIT REAL*8 (A-H,O-Z)
  23.  
  24. -INC SMINTE
  25.  
  26.       PARAMETER ( O0=0.D0, O1=1.D0, O5=.5D0, O2=2.D0, O3=3.D0, O4=4.D0 )
  27.       PARAMETER ( S=0.774596669241483377035853079956479D0)
  28.       PARAMETER ( A=0.470142064105115089770441209513447D0)
  29.       PARAMETER ( B=0.101286507323456338800987361915124D0)
  30.       PARAMETER ( C=0.333333333333333333333333333333333D0)
  31.       PARAMETER ( W234=0.0661970763942530D0 )
  32.       PARAMETER ( W567=0.0629695902724135D0 )
  33.       PARAMETER ( W234HB=W234*5.D0/9.D0 )
  34.       PARAMETER ( W234C=W234*8.D0/9.D0 )
  35.       PARAMETER ( W567HB=W567*5.D0/9.D0 )
  36.       PARAMETER ( W567C=W567*8.D0/9.D0 )
  37.       PARAMETER ( W116=1.D0/16.D0 )
  38.       PARAMETER ( W01=0.1D0 )
  39.  
  40.       DIMENSION VT(6), DVXT(6),DVYT(6)
  41.  
  42.       NBNO = 15
  43.       IF      (ISUPTH.EQ.1) THEN
  44.         NBPGAU = 15
  45.       ELSE IF (ISUPTH.EQ.2) THEN
  46.         NBPGAU = 21
  47.       ELSE IF (ISUPTH.EQ.3) THEN
  48.         NBPGAU = 1
  49.       ELSE
  50.         RETURN
  51.       ENDIF
  52.       SEGINI,MINTE
  53.  
  54. C COORDONNEES INTRINSEQUES DES 15 NOEUDS
  55.       IF (ISUPTH.EQ.1) THEN
  56.         QSIGAU( 1) = O0
  57.         QSIGAU( 2) = O5
  58.         QSIGAU( 3) = O1
  59.         QSIGAU( 4) = O5
  60.         QSIGAU( 5) = O0
  61.         QSIGAU( 6) = O0
  62.         QSIGAU( 7) = O0
  63.         QSIGAU( 8) = O1
  64.         QSIGAU( 9) = O0
  65.         QSIGAU(10) = O0
  66.         QSIGAU(11) = O5
  67.         QSIGAU(12) = O1
  68.         QSIGAU(13) = O5
  69.         QSIGAU(14) = O0
  70.         QSIGAU(15) = O0
  71.  
  72.         ETAGAU( 1) = O0
  73.         ETAGAU( 2) = O0
  74.         ETAGAU( 3) = O0
  75.         ETAGAU( 4) = O5
  76.         ETAGAU( 5) = O1
  77.         ETAGAU( 6) = O5
  78.         ETAGAU( 7) = O0
  79.         ETAGAU( 8) = O0
  80.         ETAGAU( 9) = O1
  81.         ETAGAU(10) = O0
  82.         ETAGAU(11) = O0
  83.         ETAGAU(12) = O0
  84.         ETAGAU(13) = O5
  85.         ETAGAU(14) = O1
  86.         ETAGAU(15) = O5
  87.  
  88.         DZEGAU( 1) = -O1
  89.         DZEGAU( 2) = -O1
  90.         DZEGAU( 3) = -O1
  91.         DZEGAU( 4) = -O1
  92.         DZEGAU( 5) = -O1
  93.         DZEGAU( 6) = -O1
  94.         DZEGAU( 7) =  O0
  95.         DZEGAU( 8) =  O0
  96.         DZEGAU( 9) =  O0
  97.         DZEGAU(10) =  O1
  98.         DZEGAU(11) =  O1
  99.         DZEGAU(12) =  O1
  100.         DZEGAU(13) =  O1
  101.         DZEGAU(14) =  O1
  102.         DZEGAU(15) =  O1
  103.  
  104.         r_z = O1 / 15.D0
  105.         DO I = 1, NBPGAU
  106.           POIGAU(1) = r_z
  107.         ENDDO
  108.  
  109. C COORDONNEES ET POIDS DES 21 POINTS DE GAUSS
  110.       ELSE IF (ISUPTH.EQ.2) THEN
  111.         QSIGAU( 1) = C
  112.         QSIGAU( 2) = A
  113.         QSIGAU( 3) = O1-O2*A
  114.         QSIGAU( 4) = A
  115.         QSIGAU( 5) = B
  116.         QSIGAU( 6) = O1-O2*B
  117.         QSIGAU( 7) = B
  118.         QSIGAU( 8) = C
  119.         QSIGAU( 9) = A
  120.         QSIGAU(10) = O1-O2*A
  121.         QSIGAU(11) = A
  122.         QSIGAU(12) = B
  123.         QSIGAU(13) = O1-O2*B
  124.         QSIGAU(14) = B
  125.         QSIGAU(15) = C
  126.         QSIGAU(16) = A
  127.         QSIGAU(17) = O1-O2*A
  128.         QSIGAU(18) = A
  129.         QSIGAU(19) = B
  130.         QSIGAU(20) = O1-O2*B
  131.         QSIGAU(21) = B
  132.  
  133.         ETAGAU( 1) = C
  134.         ETAGAU( 2) = A
  135.         ETAGAU( 3) = A
  136.         ETAGAU( 4) = O1-O2*A
  137.         ETAGAU( 5) = B
  138.         ETAGAU( 6) = B
  139.         ETAGAU( 7) = O1-O2*B
  140.         ETAGAU( 8) = C
  141.         ETAGAU( 9) = A
  142.         ETAGAU(10) = A
  143.         ETAGAU(11) = O1-O2*A
  144.         ETAGAU(12) = B
  145.         ETAGAU(13) = B
  146.         ETAGAU(14) = O1-O2*B
  147.         ETAGAU(15) = C
  148.         ETAGAU(16) = A
  149.         ETAGAU(17) = A
  150.         ETAGAU(18) = O1-O2*A
  151.         ETAGAU(19) = B
  152.         ETAGAU(20) = B
  153.         ETAGAU(21) = O1-O2*B
  154.  
  155.         DZEGAU( 1) = -S
  156.         DZEGAU( 2) = -S
  157.         DZEGAU( 3) = -S
  158.         DZEGAU( 4) = -S
  159.         DZEGAU( 5) = -S
  160.         DZEGAU( 6) = -S
  161.         DZEGAU( 7) = -S
  162.         DZEGAU( 8) = O0
  163.         DZEGAU( 9) = O0
  164.         DZEGAU(10) = O0
  165.         DZEGAU(11) = O0
  166.         DZEGAU(12) = O0
  167.         DZEGAU(13) = O0
  168.         DZEGAU(14) = O0
  169.         DZEGAU(15) = +S
  170.         DZEGAU(16) = +S
  171.         DZEGAU(17) = +S
  172.         DZEGAU(18) = +S
  173.         DZEGAU(19) = +S
  174.         DZEGAU(20) = +S
  175.         DZEGAU(21) = +S
  176.  
  177.         POIGAU( 1) = W116
  178.         POIGAU( 2) = W234HB
  179.         POIGAU( 3) = W234HB
  180.         POIGAU( 4) = W234HB
  181.         POIGAU( 5) = W567HB
  182.         POIGAU( 6) = W567HB
  183.         POIGAU( 7) = W567HB
  184.         POIGAU( 8) = W01
  185.         POIGAU( 9) = W234C
  186.         POIGAU(10) = W234C
  187.         POIGAU(11) = W234C
  188.         POIGAU(12) = W567C
  189.         POIGAU(13) = W567C
  190.         POIGAU(14) = W567C
  191.         POIGAU(15) = W116
  192.         POIGAU(16) = W234HB
  193.         POIGAU(17) = W234HB
  194.         POIGAU(18) = W234HB
  195.         POIGAU(19) = W567HB
  196.         POIGAU(20) = W567HB
  197.         POIGAU(21) = W567HB
  198.  
  199. C COORDONNEES DU CENTRE DE GRAVITE
  200.       ELSE IF (ISUPTH.EQ.3) THEN
  201.         QSIGAU(1) = O0
  202.         ETAGAU(1) = O0
  203.         DZEGAU(1) = O0
  204.         POIGAU(1) = O1
  205.       ENDIF
  206.  
  207. C FONCTIONS DE FORME ET LEURS DERIVEES AU(X) NBPGAU POINT(S)
  208.       DO I = 1, NBPGAU
  209.         X = QSIGAU(I)
  210.         Y = ETAGAU(I)
  211.         Z = DZEGAU(I)
  212. C- FONCTIONS DE FORME DU T6 DE BASE
  213.         XY1   = O1-X-Y
  214.         VT(1) = O2*XY1*(O5-X-Y)
  215.         VT(2) = O4*X*XY1
  216.         VT(3) = O2*X*(X-O5)
  217.         VT(4) = O4*X*Y
  218.         VT(5) = O2*Y*(Y-O5)
  219.         VT(6) = O4*XY1*Y
  220. C- FONCTIONS DE FORME DU PRISME
  221.         ZM1 = (O1-Z)*O5
  222.         ZP1 = (O1+Z)*O5
  223.         ZM2 = O1-Z*Z
  224.  
  225.         SHPTOT(1, 7,I) = ZM2*XY1
  226.         SHPTOT(1, 8,I) = ZM2*X
  227.         SHPTOT(1, 9,I) = ZM2*Y
  228.  
  229.         SHPTOT(1, 1,I) = ZM1*VT(1)-SHPTOT(1,7,I)*O5
  230.         SHPTOT(1, 2,I) = ZM1*VT(2)
  231.         SHPTOT(1, 3,I) = ZM1*VT(3)-SHPTOT(1,8,I)*O5
  232.         SHPTOT(1, 4,I) = ZM1*VT(4)
  233.         SHPTOT(1, 5,I) = ZM1*VT(5)-SHPTOT(1,9,I)*O5
  234.         SHPTOT(1, 6,I) = ZM1*VT(6)
  235.  
  236.         SHPTOT(1,10,I) = ZP1*VT(1)-SHPTOT(1,7,I)*O5
  237.         SHPTOT(1,11,I) = ZP1*VT(2)
  238.         SHPTOT(1,12,I) = ZP1*VT(3)-SHPTOT(1,8,I)*O5
  239.         SHPTOT(1,13,I) = ZP1*VT(4)
  240.         SHPTOT(1,14,I) = ZP1*VT(5)-SHPTOT(1,9,I)*O5
  241.         SHPTOT(1,15,I) = ZP1*VT(6)
  242.  
  243. C- DERIVEES DES FONCTIONS DE FORME DU T6 DE BASE
  244.         DVXT(1) = O4*(X+Y)-O3
  245.         DVXT(2) = O4*(O1-X-X-Y)
  246.         DVXT(3) = O4*X-O1
  247.         DVXT(4) = O4*Y
  248.         DVXT(5) = O0
  249.         DVXT(6) = -O4*Y
  250.  
  251.         DVYT(1) = O4*(X+Y)-O3
  252.         DVYT(2) = -O4*X
  253.         DVYT(3) = O0
  254.         DVYT(4) = O4*X
  255.         DVYT(5) = O4*Y-O1
  256.         DVYT(6) = O4*(O1-X-Y-Y)
  257.  
  258. C- DERIVEES EN X DES FONCTIONS DE FORME
  259.         SHPTOT(2, 7,I) = -ZM2
  260.         SHPTOT(2, 8,I) =  ZM2
  261.         SHPTOT(2, 9,I) =  O0
  262.  
  263.         SHPTOT(2, 1,I) = ZM1*DVXT(1)-SHPTOT(2,7,I)*O5
  264.         SHPTOT(2, 2,I) = ZM1*DVXT(2)
  265.         SHPTOT(2, 3,I) = ZM1*DVXT(3)-SHPTOT(2,8,I)*O5
  266.         SHPTOT(2, 4,I) = ZM1*DVXT(4)
  267.         SHPTOT(2, 5,I) = ZM1*DVXT(5)-SHPTOT(2,9,I)*O5
  268.         SHPTOT(2, 6,I) = ZM1*DVXT(6)
  269.  
  270.         SHPTOT(2,10,I) = ZP1*DVXT(1)-SHPTOT(2,7,I)*O5
  271.         SHPTOT(2,11,I) = ZP1*DVXT(2)
  272.         SHPTOT(2,12,I) = ZP1*DVXT(3)-SHPTOT(2,8,I)*O5
  273.         SHPTOT(2,13,I) = ZP1*DVXT(4)
  274.         SHPTOT(2,14,I) = ZP1*DVXT(5)-SHPTOT(2,9,I)*O5
  275.         SHPTOT(2,15,I) = ZP1*DVXT(6)
  276.  
  277. C- DERIVEES EN Y DES FONCTIONS DE FORME
  278.         SHPTOT(3, 7,I) = -ZM2
  279.         SHPTOT(3, 8,I) = O0
  280.         SHPTOT(3, 9,I) = ZM2
  281.  
  282.         SHPTOT(3, 1,I) = ZM1*DVYT(1)-SHPTOT(3,7,I)*O5
  283.         SHPTOT(3, 2,I) = ZM1*DVYT(2)
  284.         SHPTOT(3, 3,I) = ZM1*DVYT(3)-SHPTOT(3,8,I)*O5
  285.         SHPTOT(3, 4,I) = ZM1*DVYT(4)
  286.         SHPTOT(3, 5,I) = ZM1*DVYT(5)-SHPTOT(3,9,I)*O5
  287.         SHPTOT(3, 6,I) = ZM1*DVYT(6)
  288.  
  289.         SHPTOT(3,10,I) = ZP1*DVYT(1)-SHPTOT(3,7,I)*O5
  290.         SHPTOT(3,11,I) = ZP1*DVYT(2)
  291.         SHPTOT(3,12,I) = ZP1*DVYT(3)-SHPTOT(3,8,I)*O5
  292.         SHPTOT(3,13,I) = ZP1*DVYT(4)
  293.         SHPTOT(3,14,I) = ZP1*DVYT(5)-SHPTOT(3,9,I)*O5
  294.         SHPTOT(3,15,I) = ZP1*DVYT(6)
  295.  
  296. C- DERIVEES EN Z DES FONCTIONS DE FORME
  297.         SHPTOT(4, 7,I) = -O2*Z*XY1
  298.         SHPTOT(4, 8,I) = -O2*Z*X
  299.         SHPTOT(4, 9,I) = -O2*Z*Y
  300.  
  301.         SHPTOT(4, 1,I) = -(VT(1)+SHPTOT(4,7,I))*O5
  302.         SHPTOT(4, 2,I) = -VT(2)*O5
  303.         SHPTOT(4, 3,I) = -(VT(3)+SHPTOT(4,8,I))*O5
  304.         SHPTOT(4, 4,I) = -VT(4)*O5
  305.         SHPTOT(4, 5,I) = -(VT(5)+SHPTOT(4,9,I))*O5
  306.         SHPTOT(4, 6,I) = -VT(6)*O5
  307.  
  308.         SHPTOT(4,10,I) = (VT(1)-SHPTOT(4,7,I))*O5
  309.         SHPTOT(4,11,I) = VT(2)*O5
  310.         SHPTOT(4,12,I) = (VT(3)-SHPTOT(4,8,I))*O5
  311.         SHPTOT(4,13,I) = VT(4)*O5
  312.         SHPTOT(4,14,I) = (VT(5)-SHPTOT(4,9,I))*O5
  313.         SHPTOT(4,15,I) = VT(6)*O5
  314.  
  315.       ENDDO
  316.  
  317.       IPINTE = MINTE
  318.  
  319. C      RETURN
  320.       END
  321.  
  322.  
  323.  

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