Télécharger la plaquette Cast3M


fr en

Télécharger bfdkt.eso

Retour à la liste

Numérotation des lignes : 

bfdkt

  1. C BFDKT     SOURCE    CHAT      05/01/12    21:39:09     5004
  2.       SUBROUTINE BFDKT(QSI,ETA,GEOM,B)
  3. C=======================================================================
  4. C
  5. C     EVALUATION DE LA MATRICE B(3,9) AU POINT QSI ETA
  6. C     POUR L'ELEMENT DKT
  7. C
  8. C     CODE EBERSOLT JANVIER 84
  9. C=======================================================================
  10.       IMPLICIT INTEGER(I-N)
  11.       IMPLICIT REAL*8 (A-H,O-Z)
  12.       DIMENSION BXP(9),BYP(9),BXE(9),BYE(9)
  13.       DIMENSION GEOM(*),B(3,*)
  14. C
  15.       DATA UN/1.0D0/,DEUX/2.0D0/,QUATRE/4.0D0/,SIX/6.0D0/
  16. C
  17.       CL4=GEOM(1)
  18.       CL5=GEOM(2)
  19.       CL6=GEOM(3)
  20.       SL4=GEOM(4)
  21.       SL5=GEOM(5)
  22.       SL6=GEOM(6)
  23.       D4= GEOM(7)
  24.       D5= GEOM(8)
  25.       D6= GEOM(9)
  26.       SU4=GEOM(10)
  27.       SU5=GEOM(11)
  28.       SU6=GEOM(12)
  29.       VJ11=GEOM(13)
  30.       VJ12=GEOM(14)
  31.       VJ21=GEOM(15)
  32.       VJ22=GEOM(16)
  33.       DETJ=GEOM(17)
  34.       PS2=UN-DEUX*QSI
  35.       ET2=UN-DEUX*ETA
  36.       SU5E=SU5*ET2
  37. C---- DERIVEE DE BETAX(I) /QSI
  38.       BXP(1)=CL6*PS2+ETA*(CL5-CL6)
  39.       BXP(2)=-ETA*D5
  40.       BXP(3)=-QUATRE+SIX*(QSI+ETA)-ETA*SU5
  41.       BXP(4)=-CL6*PS2+ETA*(CL4+CL6)
  42.       BXP(5)=ETA*D4
  43.       BXP(6)=-DEUX+SIX*QSI+ETA*SU4
  44.       BXP(7)=-ETA*(CL5+CL4)
  45.       BXP(8)=ETA*(D4-D5)
  46.       BXP(9)=-ETA*(SU5-SU4)
  47. C---- DERIVEE DE BETAY(I) /QSI
  48.       BYP(1)=ETA*SL5
  49.       BYP(2)=UN-ETA*SU5
  50.       BYP(3)=-BXP(2)
  51.       BYP(4)=ETA*SL4
  52.       BYP(5)=-UN+ETA*SU4
  53.       BYP(6)=-BXP(5)
  54.       BYP(7)=-ETA*(SL5+SL4)
  55.       BYP(8)=BXP(9)
  56.       BYP(9)=-BXP(8)
  57. C---- DERIVEE DE BETAX(I) /ETA
  58.       BXE(1)=-CL5*ET2-QSI*(CL6-CL5)
  59.       BXE(2)=D5*ET2-QSI*D5
  60.       BXE(3)=-QUATRE+SIX*(QSI+ETA)+SU5E-QSI*SU5
  61.       BXE(4)=QSI*(CL4+CL6)
  62.       BXE(5)=QSI*D4
  63.       BXE(6)=QSI*SU4
  64.       BXE(7)=CL5*ET2-QSI*(CL4+CL5)
  65.       BXE(8)=D5*ET2+QSI*(D4-D5)
  66.       BXE(9)=-DEUX+SIX*ETA+SU5E+QSI*(SU4-SU5)
  67. C---- DERIVEE DE BETAY(I) /ETA
  68.       BYE(1)=-SL5*ET2+QSI*SL5
  69.       BYE(2)=UN+SU5E-QSI*SU5
  70.       BYE(3)=-BXE(2)
  71.       BYE(4)=QSI*SL4
  72.       BYE(5)=BXE(6)
  73.       BYE(6)=-BXE(5)
  74.       BYE(7)=SL5*ET2-QSI*(SL4+SL5)
  75.       BYE(8)=-UN+SU5E+QSI*(SU4-SU5)
  76.       BYE(9)=-BXE(8)
  77. C---- CALCUL DE B
  78.       DO 10 I=1,9
  79.       B(1,I)=(VJ11*BXP(I)+VJ12*BXE(I))
  80.       B(2,I)=(VJ21*BYP(I)+VJ22*BYE(I))
  81.       B(3,I)=(VJ11*BYP(I)+VJ12*BYE(I)+VJ21*BXP(I)+VJ22*BXE(I))
  82. 10    CONTINUE
  83. C
  84. C
  85.       RETURN
  86.       END
  87.  
  88.  

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