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 20 : CALL ZERO(SHP,6,9) 21 : C 22 : X21=XEL(1,2)-XEL(1,1) 23 : Y31=XEL(2,3)-XEL(2,1) 24 : DJAC=Y31*X21 25 : C 26 : GOTO (100,200,300) ,IGAU 27 : 100 CONTINUE 28 : DX=XEL(1,2)-XEL(1,1) 29 : DY=XEL(2,2)-XEL(2,1) 30 : IDEB=0 31 : IFIN=3 32 : GOTO 500 33 : 200 CONTINUE 34 : DX=XEL(1,3)-XEL(1,1) 35 : DY=XEL(2,3)-XEL(2,1) 36 : IDEB=0 37 : IFIN=6 38 : GOTO 500 39 : 300 CONTINUE 40 : DX=XEL(1,3)-XEL(1,2) 41 : DY=XEL(2,3)-XEL(2,2) 42 : IDEB=3 43 : IFIN=6 44 : GOTO 500 45 : 500 CONTINUE 46 : XLONG= DX*DX + DY*DY 47 : XLON2= SQRT(XLONG) 48 : CC=DX/XLON2 49 : SS=DY/XLON2 50 : C 51 : C 52 : SHP(1,IDEB+1)= UNDEMI 53 : SHP(2,IDEB+1)=-TRDEMI*CC/XLON2 54 : SHP(3,IDEB+1)=-TRDEMI*SS/XLON2 55 : C 56 : SHP(1,IDEB+2)=-UNHUIT*XLON2*CC 57 : SHP(2,IDEB+2)= -TRQUA*CC*SS 58 : SHP(3,IDEB+2)= - SS*SS*UNQUA+CC*CC*UNDEMI 59 : C 60 : SHP(1,IDEB+3)=-UNHUIT*XLON2*SS 61 : SHP(2,IDEB+3)= CC*CC*UNQUA-SS*SS*UNDEMI 62 : SHP(3,IDEB+3)= TRQUA*CC*SS 63 : C 64 : 65 : DO 510 IA=1,3 66 : DO 510 IB=1,3 67 : SHP(IA,IFIN+IB)=SHP(IA,IDEB+IB) 68 : 510 CONTINUE 69 : C 70 : SHP(2,IFIN+1)=-SHP(2,IDEB+1) 71 : SHP(3,IFIN+1)=-SHP(3,IDEB+1) 72 : C 73 : 74 : SHP(1,IFIN+2)=-SHP(1,IDEB+2) 75 : SHP(1,IFIN+3)=-SHP(1,IDEB+3) 76 : C 77 : RETURN 78 : END 79 :
© Cast3M 2003 - All rights reserved.
Disclaimer