Télécharger la plaquette Cast3M


fr en

Télécharger dktshp.eso

Retour à la liste

Numérotation des lignes : 

dktshp

  1. C DKTSHP    SOURCE    CHAT      05/02/25    21:15:12     5035
  2.       SUBROUTINE DKTSHP(IGAU,XEL,SHP,DJAC)
  3. C=======================================================================
  4. C   EBERSOLT OCT 85      FONCTIONS DE FORME POUR W DU DKT
  5. C  ENTREES
  6. C    IGAU = NUMERO DU POINT DE GAUSS
  7. C    XEL(3,3) = COORDONNEES LOCALES DE L ELEMENT
  8. C  SORTIES
  9. C    SHP(6,9)  =  W   W,X   W,Y EN CE POINT DE  GAUSS
  10. C    DJAC  =   JACOBIEN AU POINT DE GAUSS
  11. C    NE PAS CONFONDRE AVEC BFDKT DONNANT     THET X,X THET Y,X
  12. C                                            THET Y,X THET Y,Y
  13. C
  14. C======================================================================
  15.       IMPLICIT INTEGER(I-N)
  16.       IMPLICIT REAL*8(A-H,O-Z)
  17.       DIMENSION SHP(6,*),XEL(3,*)
  18.       DATA UNDEMI,UNQUA,TRDEMI/.5D0,.25D0,1.5D0/
  19.       DATA UNHUIT,TRQUA/.125D0,.75D0/
  20. C
  21.       CALL ZERO(SHP,6,9)
  22. C
  23.       X21=XEL(1,2)-XEL(1,1)
  24.       Y31=XEL(2,3)-XEL(2,1)
  25.       DJAC=Y31*X21
  26. C
  27.       GOTO (100,200,300) ,IGAU
  28.   100 CONTINUE
  29.       DX=XEL(1,2)-XEL(1,1)
  30.       DY=XEL(2,2)-XEL(2,1)
  31.       IDEB=0
  32.       IFIN=3
  33.       GOTO 500
  34.   200 CONTINUE
  35.       DX=XEL(1,3)-XEL(1,1)
  36.       DY=XEL(2,3)-XEL(2,1)
  37.       IDEB=0
  38.       IFIN=6
  39.       GOTO 500
  40.   300 CONTINUE
  41.       DX=XEL(1,3)-XEL(1,2)
  42.       DY=XEL(2,3)-XEL(2,2)
  43.       IDEB=3
  44.       IFIN=6
  45.       GOTO 500
  46.   500 CONTINUE
  47.       XLONG= DX*DX + DY*DY
  48.       XLON2= SQRT(XLONG)
  49.       CC=DX/XLON2
  50.       SS=DY/XLON2
  51. C
  52. C
  53.       SHP(1,IDEB+1)= UNDEMI
  54.       SHP(2,IDEB+1)=-TRDEMI*CC/XLON2
  55.       SHP(3,IDEB+1)=-TRDEMI*SS/XLON2
  56. C
  57.       SHP(1,IDEB+2)=-UNHUIT*XLON2*CC
  58.       SHP(2,IDEB+2)= -TRQUA*CC*SS
  59.       SHP(3,IDEB+2)= - SS*SS*UNQUA+CC*CC*UNDEMI
  60. C
  61.       SHP(1,IDEB+3)=-UNHUIT*XLON2*SS
  62.       SHP(2,IDEB+3)= CC*CC*UNQUA-SS*SS*UNDEMI
  63.       SHP(3,IDEB+3)= TRQUA*CC*SS
  64. C
  65.  
  66.       DO 510 IA=1,3
  67.       DO 510 IB=1,3
  68.       SHP(IA,IFIN+IB)=SHP(IA,IDEB+IB)
  69.   510 CONTINUE
  70. C
  71.       SHP(2,IFIN+1)=-SHP(2,IDEB+1)
  72.       SHP(3,IFIN+1)=-SHP(3,IDEB+1)
  73. C
  74.  
  75.       SHP(1,IFIN+2)=-SHP(1,IDEB+2)
  76.       SHP(1,IFIN+3)=-SHP(1,IDEB+3)
  77. C
  78.       RETURN
  79.       END
  80.  
  81.  

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