Télécharger la plaquette Cast3M


fr en

Télécharger tcu20.eso

Retour à la liste

Numérotation des lignes : 

tcu20

  1. C TCU20     SOURCE    OF166741  23/12/04    21:15:11     11800          
  2.  
  3. C=======================================================================
  4. C=                            T C U 2 0                                =
  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 CU20.                                 =
  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 TCU20 (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=0.5D0, O25=0.25D0, O125=0.125D0,
  27.      &            O2=2.D0, O8=8.D0 )
  28.       PARAMETER ( X=0.774596669241483377035853079956479D0 )
  29.  
  30.       NBNO = 20
  31.       IF      (ISUPTH.EQ.1) THEN
  32.         NBPGAU = 20
  33.       ELSE IF (ISUPTH.EQ.2) THEN
  34.         NBPGAU = 27
  35.       ELSE IF (ISUPTH.EQ.3) THEN
  36.         NBPGAU = 1
  37.       ELSE
  38.         CALL ERREUR(5)
  39.         RETURN
  40.       ENDIF
  41.       SEGINI,MINTE
  42.  
  43. C COORDONNEES DES 20 NOEUDS
  44.       IF (ISUPTH.EQ.1) THEN
  45.         QSIGAU( 1) = -O1
  46.         QSIGAU( 2) =  O0
  47.         QSIGAU( 3) =  O1
  48.         QSIGAU( 4) =  O1
  49.         QSIGAU( 5) =  O1
  50.         QSIGAU( 6) =  O0
  51.         QSIGAU( 7) = -O1
  52.         QSIGAU( 8) = -O1
  53.         QSIGAU( 9) = -O1
  54.         QSIGAU(10) =  O1
  55.         QSIGAU(11) =  O1
  56.         QSIGAU(12) = -O1
  57.         QSIGAU(13) = -O1
  58.         QSIGAU(14) =  O0
  59.         QSIGAU(15) =  O1
  60.         QSIGAU(16) =  O1
  61.         QSIGAU(17) =  O1
  62.         QSIGAU(18) =  O0
  63.         QSIGAU(19) = -O1
  64.         QSIGAU(20) = -O1
  65.  
  66.         ETAGAU( 1) = -O1
  67.         ETAGAU( 2) = -O1
  68.         ETAGAU( 3) = -O1
  69.         ETAGAU( 4) =  O0
  70.         ETAGAU( 5) =  O1
  71.         ETAGAU( 6) =  O1
  72.         ETAGAU( 7) =  O1
  73.         ETAGAU( 8) =  O0
  74.         ETAGAU( 9) = -O1
  75.         ETAGAU(10) = -O1
  76.         ETAGAU(11) =  O1
  77.         ETAGAU(12) =  O1
  78.         ETAGAU(13) = -O1
  79.         ETAGAU(14) = -O1
  80.         ETAGAU(15) = -O1
  81.         ETAGAU(16) =  O0
  82.         ETAGAU(17) =  O1
  83.         ETAGAU(18) =  O1
  84.         ETAGAU(19) =  O1
  85.         ETAGAU(20) =  O0
  86.  
  87.         DZEGAU( 1) = -O1
  88.         DZEGAU( 2) = -O1
  89.         DZEGAU( 3) = -O1
  90.         DZEGAU( 4) = -O1
  91.         DZEGAU( 5) = -O1
  92.         DZEGAU( 6) = -O1
  93.         DZEGAU( 7) = -O1
  94.         DZEGAU( 8) = -O1
  95.         DZEGAU( 9) =  O0
  96.         DZEGAU(10) =  O0
  97.         DZEGAU(11) =  O0
  98.         DZEGAU(12) =  O0
  99.         DZEGAU(13) =  O1
  100.         DZEGAU(14) =  O1
  101.         DZEGAU(15) =  O1
  102.         DZEGAU(16) =  O1
  103.         DZEGAU(17) =  O1
  104.         DZEGAU(18) =  O1
  105.         DZEGAU(19) =  O1
  106.         DZEGAU(20) =  O1
  107. C Poids :
  108.         r_z = O8/20.D0
  109.         DO I = 1, NBPGAU
  110.           POIGAU(I) = r_z
  111.         ENDDO
  112.  
  113. C COORDONNEES DES 27 POINTS DE GAUSS
  114.       ELSE IF (ISUPTH.EQ.2) THEN
  115.         QSIGAU( 1) = -X
  116.         QSIGAU( 2) =  O0
  117.         QSIGAU( 3) =  X
  118.         QSIGAU( 4) =  X
  119.         QSIGAU( 5) =  X
  120.         QSIGAU( 6) =  O0
  121.         QSIGAU( 7) = -X
  122.         QSIGAU( 8) = -X
  123.         QSIGAU( 9) =  O0
  124.         QSIGAU(10) = -X
  125.         QSIGAU(11) =  O0
  126.         QSIGAU(12) =  X
  127.         QSIGAU(13) =  X
  128.         QSIGAU(14) =  X
  129.         QSIGAU(15) =  O0
  130.         QSIGAU(16) = -X
  131.         QSIGAU(17) = -X
  132.         QSIGAU(18) =  O0
  133.         QSIGAU(19) = -X
  134.         QSIGAU(20) =  O0
  135.         QSIGAU(21) =  X
  136.         QSIGAU(22) =  X
  137.         QSIGAU(23) =  X
  138.         QSIGAU(24) =  O0
  139.         QSIGAU(25) = -X
  140.         QSIGAU(26) = -X
  141.         QSIGAU(27) =  O0
  142.  
  143.         ETAGAU( 1) = -X
  144.         ETAGAU( 2) = -X
  145.         ETAGAU( 3) = -X
  146.         ETAGAU( 4) =  O0
  147.         ETAGAU( 5) =  X
  148.         ETAGAU( 6) =  X
  149.         ETAGAU( 7) =  X
  150.         ETAGAU( 8) =  O0
  151.         ETAGAU( 9) =  O0
  152.         ETAGAU(10) = -X
  153.         ETAGAU(11) = -X
  154.         ETAGAU(12) = -X
  155.         ETAGAU(13) =  O0
  156.         ETAGAU(14) =  X
  157.         ETAGAU(15) =  X
  158.         ETAGAU(16) =  X
  159.         ETAGAU(17) =  O0
  160.         ETAGAU(18) =  O0
  161.         ETAGAU(19) = -X
  162.         ETAGAU(20) = -X
  163.         ETAGAU(21) = -X
  164.         ETAGAU(22) =  O0
  165.         ETAGAU(23) =  X
  166.         ETAGAU(24) =  X
  167.         ETAGAU(25) =  X
  168.         ETAGAU(26) =  O0
  169.         ETAGAU(27) =  O0
  170.  
  171.         DZEGAU( 1) = -X
  172.         DZEGAU( 2) = -X
  173.         DZEGAU( 3) = -X
  174.         DZEGAU( 4) = -X
  175.         DZEGAU( 5) = -X
  176.         DZEGAU( 6) = -X
  177.         DZEGAU( 7) = -X
  178.         DZEGAU( 8) = -X
  179.         DZEGAU( 9) = -X
  180.         DZEGAU(10) =  O0
  181.         DZEGAU(11) =  O0
  182.         DZEGAU(12) =  O0
  183.         DZEGAU(13) =  O0
  184.         DZEGAU(14) =  O0
  185.         DZEGAU(15) =  O0
  186.         DZEGAU(16) =  O0
  187.         DZEGAU(17) =  O0
  188.         DZEGAU(18) =  O0
  189.         DZEGAU(19) =  X
  190.         DZEGAU(20) =  X
  191.         DZEGAU(21) =  X
  192.         DZEGAU(22) =  X
  193.         DZEGAU(23) =  X
  194.         DZEGAU(24) =  X
  195.         DZEGAU(25) =  X
  196.         DZEGAU(26) =  X
  197.         DZEGAU(27) =  X
  198. C POIDS
  199. C 5/9 POUR UNE COORDONNEE EGALE A + OU - X
  200. C 8/9 POUR UNE COORDONNEE NULLE
  201.         W5   = 125.D0/729.D0
  202.         W8   = 200.D0/729.D0
  203.         W88  = 320.D0/729.D0
  204.         W888 = 512.D0/729.D0
  205.         POIGAU( 1) = W5
  206.         POIGAU( 2) = W8
  207.         POIGAU( 3) = W5
  208.         POIGAU( 4) = W8
  209.         POIGAU( 5) = W5
  210.         POIGAU( 6) = W8
  211.         POIGAU( 7) = W5
  212.         POIGAU( 8) = W8
  213.         POIGAU( 9) = W88
  214.  
  215.         POIGAU(10) = W8
  216.         POIGAU(11) = W88
  217.         POIGAU(12) = W8
  218.         POIGAU(13) = W88
  219.         POIGAU(14) = W8
  220.         POIGAU(15) = W88
  221.         POIGAU(16) = W8
  222.         POIGAU(17) = W88
  223.         POIGAU(18) = W888
  224.  
  225.         POIGAU(19) = W5
  226.         POIGAU(20) = W8
  227.         POIGAU(21) = W5
  228.         POIGAU(22) = W8
  229.         POIGAU(23) = W5
  230.         POIGAU(24) = W8
  231.         POIGAU(25) = W5
  232.         POIGAU(26) = W8
  233.         POIGAU(27) = W88
  234.  
  235. C COORDONNEES DU CENTRE DE GRAVITE
  236.       ELSE IF (ISUPTH.EQ.3) THEN
  237.         QSIGAU(1) = O0
  238.         ETAGAU(1) = O0
  239.         DZEGAU(1) = O0
  240.         POIGAU(1) = O8
  241.       ENDIF
  242.  
  243. C FONCTIONS DE FORME ET LEURS DERIVEES AU(X) NBPGAU POINT(S)
  244.       DO I = 1, NBPGAU
  245.         XI = QSIGAU(I)
  246.         YI = ETAGAU(I)
  247.         ZI = DZEGAU(I)
  248. C--- FONCTIONS DE FORME
  249.         SHPTOT(1, 1,I) = -O125*(O1-XI)*(O1-YI)*(O1-ZI)*( XI+YI+ZI+O2)
  250.         SHPTOT(1, 2,I) =  O25 *(O1-XI*XI)*(O1-YI  )*(O1-ZI)
  251.         SHPTOT(1, 3,I) = -O125*(O1+XI)*(O1-YI)*(O1-ZI)*(-XI+YI+ZI+O2)
  252.         SHPTOT(1, 4,I) =  O25 *(O1+XI  )*(O1-YI*YI)*(O1-ZI)
  253.         SHPTOT(1, 5,I) = -O125*(O1+XI)*(O1+YI)*(O1-ZI)*(-XI-YI+ZI+O2)
  254.         SHPTOT(1, 6,I) =  O25 *(O1-XI*XI)*(O1+YI  )*(O1-ZI)
  255.         SHPTOT(1, 7,I) = -O125*(O1-XI)*(O1+YI)*(O1-ZI)*( XI-YI+ZI+O2)
  256.         SHPTOT(1, 8,I) =  O25 *(O1-XI  )*(O1-YI*YI)*(O1-ZI)
  257.         SHPTOT(1, 9,I) =  O25 *(O1-XI)*(O1-YI)*(O1-ZI*ZI)
  258.         SHPTOT(1,10,I) =  O25 *(O1+XI)*(O1-YI)*(O1-ZI*ZI)
  259.         SHPTOT(1,11,I) =  O25 *(O1+XI)*(O1+YI)*(O1-ZI*ZI)
  260.         SHPTOT(1,12,I) =  O25 *(O1-XI)*(O1+YI)*(O1-ZI*ZI)
  261.         SHPTOT(1,13,I) = -O125*(O1-XI)*(O1-YI)*(O1+ZI)*( XI+YI-ZI+O2)
  262.         SHPTOT(1,14,I) =  O25 *(O1-XI*XI)*(O1-YI  )*(O1+ZI)
  263.         SHPTOT(1,15,I) = -O125*(O1+XI)*(O1-YI)*(O1+ZI)*(-XI+YI-ZI+O2)
  264.         SHPTOT(1,16,I) =  O25 *(O1+XI  )*(O1-YI*YI)*(O1+ZI)
  265.         SHPTOT(1,17,I) = -O125*(O1+XI)*(O1+YI)*(O1+ZI)*(-XI-YI-ZI+O2)
  266.         SHPTOT(1,18,I) =  O25 *(O1-XI*XI)*(O1+YI  )*(O1+ZI)
  267.         SHPTOT(1,19,I) = -O125*(O1-XI)*(O1+YI)*(O1+ZI)*( XI-YI-ZI+O2)
  268.         SHPTOT(1,20,I) =  O25 *(O1-XI  )*(O1-YI*YI)*(O1+ZI)
  269. C---  DERIVEES EN X DES FONCTIONS DE FORME
  270.         SHPTOT(2, 1,I) =  O125*(O1-YI)*(O1-ZI)*( O2*XI+YI+ZI+O1)
  271.         SHPTOT(2, 2,I) = -O5  *XI*(O1-YI  )*(O1-ZI)
  272.         SHPTOT(2, 3,I) = -O125*(O1-YI)*(O1-ZI)*(-O2*XI+YI+ZI+O1)
  273.         SHPTOT(2, 4,I) =  O25 *  (O1-YI*YI)*(O1-ZI)
  274.         SHPTOT(2, 5,I) = -O125*(O1+YI)*(O1-ZI)*(-O2*XI-YI+ZI+O1)
  275.         SHPTOT(2, 6,I) = -O5  *XI*(O1+YI  )*(O1-ZI)
  276.         SHPTOT(2, 7,I) =  O125*(O1+YI)*(O1-ZI)*( O2*XI-YI+ZI+O1)
  277.         SHPTOT(2, 8,I) = -O25 *  (O1-YI*YI)*(O1-ZI)
  278.         SHPTOT(2, 9,I) = -O25 *(O1-YI)*(O1-ZI*ZI)
  279.         SHPTOT(2,10,I) =  O25 *(O1-YI)*(O1-ZI*ZI)
  280.         SHPTOT(2,11,I) =  O25 *(O1+YI)*(O1-ZI*ZI)
  281.         SHPTOT(2,12,I) = -O25 *(O1+YI)*(O1-ZI*ZI)
  282.         SHPTOT(2,13,I) =  O125*(O1-YI)*(O1+ZI)*( O2*XI+YI-ZI+O1)
  283.         SHPTOT(2,14,I) = -O5  *XI*(O1-YI  )*(O1+ZI)
  284.         SHPTOT(2,15,I) = -O125*(O1-YI)*(O1+ZI)*(-O2*XI+YI-ZI+O1)
  285.         SHPTOT(2,16,I) =  O25 *  (O1-YI*YI)*(O1+ZI)
  286.         SHPTOT(2,17,I) = -O125*(O1+YI)*(O1+ZI)*(-O2*XI-YI-ZI+O1)
  287.         SHPTOT(2,18,I) = -O5  *XI*(O1+YI  )*(O1+ZI)
  288.         SHPTOT(2,19,I) =  O125*(O1+YI)*(O1+ZI)*( O2*XI-YI-ZI+O1)
  289.         SHPTOT(2,20,I) = -O25 *  (O1-YI*YI)*(O1+ZI)
  290. C---  DERIVEES EN Y DES FONCTIONS DE FORME
  291.         SHPTOT(3, 1,I) =  O125*(O1-XI)*(O1-ZI)*( XI+O2*YI+ZI+O1)
  292.         SHPTOT(3, 2,I) = -O25 *(O1-XI*XI)  *(O1-ZI)
  293.         SHPTOT(3, 3,I) =  O125*(O1+XI)*(O1-ZI)*(-XI+O2*YI+ZI+O1)
  294.         SHPTOT(3, 4,I) = -O5  *(O1+  XI)*YI*(O1-ZI)
  295.         SHPTOT(3, 5,I) = -O125*(O1+XI)*(O1-ZI)*(-XI-O2*YI+ZI+O1)
  296.         SHPTOT(3, 6,I) =  O25 *(O1-XI*XI)  *(O1-ZI)
  297.         SHPTOT(3, 7,I) = -O125*(O1-XI)*(O1-ZI)*( XI-O2*YI+ZI+O1)
  298.         SHPTOT(3, 8,I) = -O5  *(O1-  XI)*YI*(O1-ZI)
  299.         SHPTOT(3, 9,I) = -O25 *(O1-XI)*(O1-ZI*ZI)
  300.         SHPTOT(3,10,I) = -O25 *(O1+XI)*(O1-ZI*ZI)
  301.         SHPTOT(3,11,I) =  O25 *(O1+XI)*(O1-ZI*ZI)
  302.         SHPTOT(3,12,I) =  O25 *(O1-XI)*(O1-ZI*ZI)
  303.         SHPTOT(3,13,I) =  O125*(O1-XI)*(O1+ZI)*( XI+O2*YI-ZI+O1)
  304.         SHPTOT(3,14,I) = -O25 *(O1-XI*XI)  *(O1+ZI)
  305.         SHPTOT(3,15,I) =  O125*(O1+XI)*(O1+ZI)*(-XI+O2*YI-ZI+O1)
  306.         SHPTOT(3,16,I) = -O5  *(O1+  XI)*YI*(O1+ZI)
  307.         SHPTOT(3,17,I) = -O125*(O1+XI)*(O1+ZI)*(-XI-O2*YI-ZI+O1)
  308.         SHPTOT(3,18,I) =  O25 *(O1-XI*XI)  *(O1+ZI)
  309.         SHPTOT(3,19,I) = -O125*(O1-XI)*(O1+ZI)*( XI-O2*YI-ZI+O1)
  310.         SHPTOT(3,20,I) = -O5  *(O1-  XI)*YI*(O1+ZI)
  311. C---  DERIVEES EN Z DES FONCTIONS DE FORME
  312.         SHPTOT(4, 1,I) =  O125*(O1-XI)*(O1-YI)*( XI+YI+O2*ZI+O1)
  313.         SHPTOT(4, 2,I) = -O25 *(O1-XI*XI)*(O1-  YI)
  314.         SHPTOT(4, 3,I) =  O125*(O1+XI)*(O1-YI)*(-XI+YI+O2*ZI+O1)
  315.         SHPTOT(4, 4,I) = -O25 *(O1+  XI)*(O1-YI*YI)
  316.         SHPTOT(4, 5,I) =  O125*(O1+XI)*(O1+YI)*(-XI-YI+O2*ZI+O1)
  317.         SHPTOT(4, 6,I) = -O25 *(O1-XI*XI)*(O1+  YI)
  318.         SHPTOT(4, 7,I) =  O125*(O1-XI)*(O1+YI)*( XI-YI+O2*ZI+O1)
  319.         SHPTOT(4, 8,I) = -O25 *(O1-  XI)*(O1-YI*YI)
  320.         SHPTOT(4, 9,I) = -O5 *(O1-XI)*(O1-YI)*ZI
  321.         SHPTOT(4,10,I) = -O5 *(O1+XI)*(O1-YI)*ZI
  322.         SHPTOT(4,11,I) = -O5 *(O1+XI)*(O1+YI)*ZI
  323.         SHPTOT(4,12,I) = -O5 *(O1-XI)*(O1+YI)*ZI
  324.         SHPTOT(4,13,I) = -O125*(O1-XI)*(O1-YI)*( XI+YI-O2*ZI+O1)
  325.         SHPTOT(4,14,I) =  O25 *(O1-XI*XI)*(O1-  YI)
  326.         SHPTOT(4,15,I) = -O125*(O1+XI)*(O1-YI)*(-XI+YI-O2*ZI+O1)
  327.         SHPTOT(4,16,I) =  O25 *(O1+  XI)*(O1-YI*YI)
  328.         SHPTOT(4,17,I) = -O125*(O1+XI)*(O1+YI)*(-XI-YI-O2*ZI+O1)
  329.         SHPTOT(4,18,I) =  O25 *(O1-XI*XI)*(O1+  YI)
  330.         SHPTOT(4,19,I) = -O125*(O1-XI)*(O1+YI)*( XI-YI-O2*ZI+O1)
  331.         SHPTOT(4,20,I) =  O25 *(O1-  XI)*(O1-YI*YI)
  332.       ENDDO
  333.  
  334.       IPINTE = MINTE
  335.  
  336. c      RETURN
  337.       END
  338.  
  339.  
  340.  

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