Télécharger la plaquette Cast3M


fr en

Télécharger geradk.eso

Retour à la liste

Numérotation des lignes : 

geradk

  1. C GERADK    SOURCE    CHAT      05/01/13    00:17:55     5004
  2.       SUBROUTINE GERADK(XEL,GEOM)
  3. C |====================================================================|
  4. C |                                                                    |
  5. C |        CALCULE DES CONSTANTES D UN TRIANGLE POUR L ELEMENT DKT     |
  6. C |                                                                    |
  7. C |==== ENTREES                                                        |
  8. C |        XEL(3,3)= COORDONNEES DE L ELEMENT                          |
  9. C |==== SORTIES                                                        |
  10. C |        GEOM(17)=DES CONSTANTES POUR LE DKT                         |
  11. C |   CODE SUO X.Z.                                                    |
  12. C |====================================================================|
  13.       IMPLICIT INTEGER(I-N)
  14.       IMPLICIT REAL*8(A-H,O-Z)
  15.       DIMENSION XEL(3,*),GEOM(*)
  16.       DATA TRSQUA/0.75D0/,UNQUA/0.25D0/,UNDEMI/0.5D0/
  17.       X12=XEL(1,1)-XEL(1,2)
  18.       X23=XEL(1,2)-XEL(1,3)
  19.       X31=XEL(1,3)-XEL(1,1)
  20.       Y12=XEL(2,1)-XEL(2,2)
  21.       Y23=XEL(2,2)-XEL(2,3)
  22.       Y31=XEL(2,3)-XEL(2,1)
  23.       XSQ12=X12*X12
  24.       XSQ23=X23*X23
  25.       XSQ31=X31*X31
  26.       YSQ12=Y12*Y12
  27.       YSQ23=Y23*Y23
  28.       YSQ31=Y31*Y31
  29.       DIST12=XSQ12+YSQ12
  30.       DIST23=XSQ23+YSQ23
  31.       DIST31=XSQ31+YSQ31
  32.       A4=-X23/DIST23
  33.       A5=-X31/DIST31
  34.       A6=-X12/DIST12
  35.       D4=-Y23/DIST23
  36.       D5=-Y31/DIST31
  37.       D6=-Y12/DIST12
  38.       B4=TRSQUA*X23*Y23/DIST23
  39.       B5=TRSQUA*X31*Y31/DIST31
  40.       B6=TRSQUA*X12*Y12/DIST12
  41.       C4=(UNQUA*XSQ23-UNDEMI*YSQ23)/DIST23
  42.       C5=(UNQUA*XSQ31-UNDEMI*YSQ31)/DIST31
  43.       C6=(UNQUA*XSQ12-UNDEMI*YSQ12)/DIST12
  44.       E4=(UNQUA*YSQ23-UNDEMI*XSQ23)/DIST23
  45.       E5=(UNQUA*YSQ31-UNDEMI*XSQ31)/DIST31
  46.       E6=(UNQUA*YSQ12-UNDEMI*XSQ12)/DIST12
  47.       GEOM(1)=A4
  48.       GEOM(2)=A5
  49.       GEOM(3)=A6
  50.       GEOM(4)=B4
  51.       GEOM(5)=B5
  52.       GEOM(6)=B6
  53.       GEOM(7)=C4
  54.       GEOM(8)=C5
  55.       GEOM(9)=C6
  56.       GEOM(10)=D4
  57.       GEOM(11)=D5
  58.       GEOM(12)=D6
  59.       GEOM(13)=E4
  60.       GEOM(14)=E5
  61.       GEOM(15)=E6
  62.       GEOM(16)=1.D0
  63.       GEOM(17)=1.D0
  64.       RETURN
  65.       END
  66.  
  67.  

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