用常应变三角形单元解弹性力学平面问题的程序(5)

IF(IL.GE.IDH) THEN

IDL=IL-IDH+1

KS(IDH,IDL)=KS(IDH,IDL)+KE(IH,L)

ENDIF

10 CONTINUE

END

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

C SUBPROGRAM-7

C CALCULATE NODAL LOAD VECTOR

SUBROUTINE EQUPE(NJ,NE,NPJ,NJ2,T,V,LND,X,Y,PJ,P) DIMENSION LND(NE,3),X(NJ),Y(NJ),PJ(NPJ,3),P(NJ2) DO 10 I=1,NJ2

10 P(I)=0.

DO 20 I=1,NPJ

II=PJ(I,1)

P(2*II-1)=PJ(I,2)

20 P(2*II)=PJ(I,3)

30 IF(V.EQ.0.) GOTO 50

DO 40 IE=1,NE

CALL ATE(IE,NJ,NE,LND,X,Y,AE)

PE=-V*AE*T/3.

DO 40 I=1,3

II=LND(IE,I)

40 P(2*II)=P(2*II)+PE

50 RETURN

END

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

C SUBPROGRAM-8

C INTRODUCE BOUNDARY CONDITION

SUBROUTINE INSCD(NS,NW,NJ2,JR,KS,P)

DIMENSION P(NJ2),JR(NS,3)

REAL KS(NJ2,NW)

DO 30 I=1,NS

IR=JR(I,1)

DO 30 J=2,3

IF(JR(I,J).EQ.0) GOTO 30

II=2*IR+J-3

KS(II,1)=1.

DO 10 JJ=2,NW

10 KS(II,JJ)=0.

IF(II.GT.NW) JO=NW

IF(II.LE.NW) JO=II

DO 20 JJ=2,JO

20 KS(II-JJ+1,JJ)=0.

P(II)=0.