Télécharger la plaquette Cast3M


fr en

Télécharger ndkt.eso

Retour à la liste

Numérotation des lignes : 

ndkt

  1. C NDKT      SOURCE    CHAT      05/01/13    01:56:26     5004
  2.       SUBROUTINE NDKT (IGAU,XEL,EXCENT,SHPTOT,SHP,XGENE,DJAC)
  3. ************************************************************************
  4. *
  5. *                             N D K T
  6. *                             -------
  7. *
  8. * FONCTION:
  9. * ---------
  10. *
  11. *     CALCUL DE LA MATRICE "N" EN UN POINT DE GAUSS DE L'ELEMENT "DKT",
  12. *     POUR LE CALCUL DE LA MATRICE-MASSE = SOMME SUR L'ELEMENT DE "NT.N"
  13. *
  14. *     FORME DE LA MATRICE "N":
  15. *
  16. *        A1 00 E1 F1 G1 00  A2 00 E2 F2 G2 00  A3 00 E3 F3 G3 00
  17. *        00 A1 H1 I1 J1 00  00 A2 H2 I2 J2 00  00 A3 H3 I3 J3 00
  18. *        00 00 B1 C1 D1 00  00 00 B2 C2 D2 00  00 00 B3 C3 D3 00
  19. *
  20. *     LES TERMES "E" A "J" SONT NON NULS LORSQUE L'ELEMENT EST EXCENTRE.
  21. *
  22. * MODULES UTILISES:
  23. * -----------------
  24. *
  25.       IMPLICIT INTEGER(I-N)
  26. -INC CCREEL
  27. *
  28. * PARAMETRES:   (E)=ENTREE   (S)=SORTIE   (+ = CONTENU DANS UN COMMUN)
  29. * -----------
  30. *
  31. *     IGAU    (E)  NUMERO DU POINT DE GAUSS.
  32. *     XEL     (E)  COORDONNEES LOCALES DES NOEUDS DE L'ELEMENT
  33. *                  (1 NOEUD PAR COLONNE).
  34. *     EXCENT  (E)  VALEUR D'EXCENTREMENT DE L'ELEMENT.
  35. *     SHPTOT  (E)  FONCTIONS DE FORMES ET DERIVEES (DANS L'ORDRE USUEL).
  36. *     SHP     E/S  TABLEAU DE TRAVAIL.
  37. *     XGENE   (S)  MATRICE "N".
  38. *     DJAC    (S)  JACOBIEN.
  39. *
  40.       INTEGER IGAU
  41.       REAL*8 EXCENT,DJAC
  42.       REAL*8 XEL(3,3),SHPTOT(6,3,*),SHP(6,3),XGENE(3,18)
  43. *
  44. * CONSTANTES:
  45. * -----------
  46. *
  47. *     XX, YY = COORDONNEES DANS LE TRIANGLE DE REFERENCE DES POINTS DE
  48. *              GAUSS.
  49. *
  50.       PARAMETER (NBNO = 3)
  51.       REAL*8 XX(3),YY(3)
  52. *
  53. * VARIABLES:
  54. * ----------
  55. *
  56.       REAL*8 BB(18)
  57. *
  58. * AUTEUR, DATE DE CREATION:
  59. * -------------------------
  60. *
  61. *     PASCAL MANIGOT     15 JUIN 1987
  62. *
  63. * LANGAGE:
  64. * --------
  65. *
  66. *     FORTRAN77 + EXTENSION: DECLARATION "REAL*8".
  67. *
  68. ************************************************************************
  69. *
  70.       DATA XX/.5D0,.0D0,.5D0/
  71.       DATA YY/.0D0,.5D0,.5D0/
  72. *
  73.       CALL ZERO(XGENE,3,18)
  74. *
  75.       DO 127 NPOI=1,NBNO
  76.       SHP(1,NPOI)=SHPTOT(1,NPOI,IGAU)
  77.       SHP(2,NPOI)=SHPTOT(2,NPOI,IGAU)
  78.   127 SHP(3,NPOI)=SHPTOT(3,NPOI,IGAU)
  79.       CALL JACOBI(XEL,SHP,2,NBNO,DJAC)
  80. *     RQ: A CE NIVEAU, SHP(2,*) ET SHP(3,*) CONTIENNENT LES DERIVEES
  81. *         DANS L'ESPACE GEOMETRIQUE.
  82.       K=1
  83.       DO 227 NPOI=1,NBNO
  84.       XGENE(1,K  )=SHP(1,NPOI)
  85.       XGENE(2,K+1)=SHP(1,NPOI)
  86.   227 K=K+6
  87.       CALL MFDKT (XX(IGAU),YY(IGAU),XEL,BB)
  88.       KK=3
  89.       IA=1
  90.       DO 327 NPOI=1,NBNO
  91.       XGENE(3,KK  )=BB(IA  )
  92.       XGENE(3,KK+1)=BB(IA+1)
  93.       XGENE(3,KK+2)=BB(IA+2)
  94.       KK=KK+6
  95.       IA=IA+3
  96.  327  CONTINUE
  97. *
  98.       IF (ABS(EXCENT) .GT. XPETIT) THEN
  99.          CALL MFDKT2 (XX(IGAU),YY(IGAU),XEL,BB)
  100.          KK=3
  101.          IA=1
  102. *
  103. *  CORRECTION P. MANIGOT  : CHANGEMENT DE SIGNE DES EXCENTREMENTS
  104. *
  105.          DO 427 NPOI=1,NBNO
  106.          XGENE(1,KK  )= EXCENT*BB(IA)
  107.          XGENE(1,KK+1)= EXCENT*BB(IA+1)
  108.          XGENE(1,KK+2)= EXCENT*BB(IA+2)
  109.          XGENE(2,KK  )= EXCENT*BB(IA+ 9)
  110.          XGENE(2,KK+1)= EXCENT*BB(IA+10)
  111.          XGENE(2,KK+2)= EXCENT*BB(IA+11)
  112.          KK=KK+6
  113.          IA=IA+3
  114.  427     CONTINUE
  115.       END IF
  116. *
  117.       END
  118.  
  119.  
  120.  
  121.  

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