344
Appendix B
C
CI=X(MK)-X(MJ)
CJ=X(MI)-X(MK)
TRS=ABS(CJ*BI-BJ*CI)
C Calculating PPR and PAD
C
C
C
C
C
TRS1=TRS/6.
PSE(MI)=PSE(MI)+ TRS1*SE(N)
PSE(MJ)=PSE(MJ)+ TRS1*SE(N)
PSE(MK)=PSE(MK)+ TRS1*SE(N)
PPR(MI)=PPR(MI)+ TRS1*PR(N)
PPR(MJ)=PPR(MJ)+ TRS1*PR(N)
PPR(MK)=PPR(MK)+ TRS1*PR(N)
PAD(MI)=PAD(MI) + TRS1*AD(N)
PAD(MJ)=PAD(MJ)+ TRS1*AD(N)
PAD(MK)=PAD(MK)+ TRS1*AD(N)
IF (IHC .EQ. 1) THEN
AVHO=(HO(MI)+HO(MJ)+HO(MK))/3.
AVTR=(TR(MI)+TR(MJ)+TR(MK))/3.
AVBR=(BR(MI)+BR(MJ)+BR(MK))/3.
ELSE
IF (AVHO .GT. AVTR) THEN
AVTH=AVTR-AVBR
SS=SE(N)
ELSE
AVTH=AVHO-AVBR
SS=PR(N)
ENDIF
TT=TK(N)*AVTH
IF (ISO .EQ. 0) THEN
D11=TT
D12=0.
D22=TT
ELSE
Dll=DXX(N) *AVTH
D12=DXY(N) *AVTH
D22=DYY(N) *AVTH
ENDIF
SS=l.
011=OXX(N)
012=OXY(N)
022=OYY(N)
ENDIF
00 J=l, 3
ME(J)=IJM(J,N)
END 00
DO L=l, 3
MI=ME(l)
MJ=ME(2)
MK=ME(3)
XM=(X(MI)+X(MJ)+X(MK))/3
YM=(Y(MI)+Y(MJ)+Y(MK))/3
BI=Y(MJ)-Y(MK)
BJ=Y(MK) -Y(MI)
BK=Y(MI)-Y(MJ)
CI=X(MK)-X(MJ)
CJ=X(MI)-X(MK)