SUBROUTINE BTE(IE,NJ,NE,LND,X,Y,AE,B)

DIMENSION LND(NE,3),X(NJ),Y(NJ),B(3,6)

I=LND(IE,1)

J=LND(IE,2)

K=LND(IE,3)

DO 10 II=1,3

DO 10 JJ=1,6

10 B(II,JJ)=0.

B(1,1)=Y(J)-Y(K)

B(1,3)=Y(K)-Y(I)

B(1,5)=Y(I)-Y(J)

B(2,2)=X(K)-X(J)

B(2,4)=X(I)-X(K)

B(2,6)=X(J)-X(I)

B(3,1)=B(2,2)

B(3,2)=B(1,1)

B(3,3)=B(2,4)

B(3,4)=B(1,3)

B(3,5)=B(2,6)

B(3,6)=B(1,5)

DO 20 I1=1,3

DO 20 J1=1,6

20 B(I1,J1)=.5/AE*B(I1,J1)

END

*-------------------------------------------------------

C SUBPROGRAM-6

C CALCULATE GLOBAL STIFFNESS MATRIX

SUBROUTINE STIFF(NJ,NE,NJ2,NW,LND,X,Y,E,PR,T,KS) DIMENSION LND(NE,3),X(NJ),Y(NJ)

REAL KS(NJ2,NW),KE(6,6)

DO 5 I=1,NJ2

DO 5 J=1,NW

5 KS(I,J)=0.

DO 10 IE=1,NE

CALL STE(IE,NJ,NE,LND,X,Y,E,PR,T,KE)

DO 10 I=1,3

IZ=LND(IE,I)

DO 10 II=1,2

IH =2*(I -1)+II

IDH=2*(IZ-1)+II

DO 10 J=1,3

JZ=LND(IE,J)

DO 10 JJ=1,2

L =2*(J -1)+JJ

IL=2*(JZ-1)+JJ