Télécharger la plaquette Cast3M


fr en

Télécharger elgaus.eso

Retour à la liste

Numérotation des lignes : 

elgaus

  1. C ELGAUS    SOURCE    CHAT      05/01/12    23:33:00     5004
  2.       SUBROUTINE ELGAUS(MINTE,MMAT1,SGAUSS,IFOIS,IFOI2)
  3.       IMPLICIT INTEGER(I-N)
  4.       IMPLICIT REAL*8(A-H,O-Z)
  5. -INC CCREEL
  6. -INC SMINTE
  7.  
  8.       SEGMENT MMAT1
  9.          REAL*8 VALMAT(NMATR)
  10.          REAL*8 XE(3,NBNN),XE1(3,NBNN)
  11.          REAL*8 SHP(6,NBNN),GRAD(NDIM,NBNN,NBPGAU)
  12.          REAL*8 COSD1(3),COSD2(3)
  13.       ENDSEGMENT
  14.       SEGMENT SGAUSS
  15.          REAL*8 XGAUSS(3,NBPGAU)
  16.          REAL*8 DX(NBPGAU)
  17.       ENDSEGMENT
  18. *
  19.       NBPGAU=DX(/1)
  20.       NBNN=XE(/2)
  21.       NDIM=GRAD(/1)
  22.       IDIM=NDIM
  23.       NDIM2=IDIM-1
  24.       IFOIS=0
  25.       IFOI2=0
  26.       DO 20 IGAU=1,NBPGAU
  27. *
  28. *      CALCUL DE LA MATRICE GRADIENT DES FONCTIONS DE FORME ET
  29. *      DU JACOBIEN(DANS LE PLAN), EN UN POINT DE GAUSS
  30. *
  31.        DO 90 NP=1,NBNN
  32.         DO 90  I=1,IDIM
  33.           SHP(I,NP)=SHPTOT(I,NP,IGAU)
  34.   90   CONTINUE
  35. *      COORDONNEES GLOBALES DU POINT DE GAUSS DANS
  36. *      L'ELEMENT COURANT
  37.        DO 70 I=1,IDIM
  38.          XGAUSS(I,IGAU)=0.
  39.          DO 70 NP=1,NBNN
  40.            XGAUSS(I,IGAU)=XGAUSS(I,IGAU)+XE(I,NP)*SHP(1,NP)
  41.   70   CONTINUE
  42. *
  43. *      DERIVES DES FONCTIONS DE FORME DANS LA GEOMETRIE REELLE
  44. *      LOCALE (2D) ET LE JACOBIEN
  45.        CALL JACOBI(XE1,SHP,NDIM2,NBNN,DJAC)
  46. *
  47. *      ON FAIT TOURNER LE GRADIENT DE -PI/2 DANS LE REPERE LOCAL
  48. *      POUR ETRE PARALLELE AU COTE OPPOSE AU SOMMET
  49.        DO 60 NP=1,NBNN
  50.          XG=SHP(2,NP)
  51.          YG=SHP(3,NP)
  52.          SHP(2,NP)=YG
  53.          SHP(3,NP)=-XG
  54.    60  CONTINUE
  55. *
  56. *      ON REVIENT DANS LE REPERE GLOBAL 3D
  57.        DO 100 NP=1,NBNN
  58.          DO 100 I= 1,IDIM
  59.            GRAD(I,NP,IGAU)=SHP(2,NP)*COSD1(I)+SHP(3,NP)*COSD2(I)
  60.   100  CONTINUE
  61. *
  62.        IF(DJAC.LT.XZERO)IFOIS=IFOIS+1
  63.        IF(ABS(DJAC).LT.XPETIT)IFOI2=IFOI2 +1
  64. *
  65. *      ON MULTIPLIE LE JACOBIEN PAR LE POIDS D'INTEGRATION,POUR LE
  66. *      POINT DE GAUSS CONSIDERE
  67. *
  68.        DJAC=ABS(DJAC)*POIGAU(IGAU)
  69.        DX(IGAU)=DJAC
  70. *
  71.    20 CONTINUE
  72.       RETURN
  73.       END
  74.  
  75.  
  76.  
  77.  
  78.  

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