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 :
© Cast3M 2003 - All rights reserved.
Disclaimer