32
EIJ(7,13)=0.805823D0
EIJ(2,3)=1.03502D0
UIJ(2,3)=1.0805D0
KIJ(2,3)=1.00046D0
EIJ(2,17)=1.16446D0
UIJ(2,17)=1.61666D0
KIJ(2,17)=1.02034D0
UIJ(2,5)=1.25D0
EIJ(2,4)=1.01306D0
UIJ(2,4)=1.25D0
UIJ(2,10)=1.25D0
EIJ(2,9)=1.00532D0
UIJ(2,9)=1.25D0
EIJ(3,17)=1.034787D0
EIJ(3,4)=1.0049D0
EIJ(17,18)=1.1D0
EIJ(5,17)=1.3D0
EIJ(4,17)=1.3D0
RETURN
END
SUBROUTINE PARMI1
IMPLICIT REAL*8(A-H,O-Z)
REAL*8 KI,KIJ,KM
INTEGER GN,QN,FN
DIMENSION EIJM(19,19),GIJM(19,19)
COMMON/Y1/Y(19)/NC1/NC/EFI/EI(19),KI(19),GI(19),QI(19),FI(19)
*/INTER1/EIJ(19,19),UIJ(19,19),KIJ(19,19),GIJ(19,19)
*/KM/KM/COEF1/B1(13),C1(53)/AN/AN(53)
*/GQFN/GN(53),QN(53),FN(53)/UN/UN(53)
DO 1 I=1,NC
EIJM(I,I)=EI(I)
GIJM(I,I)=GI(I)
DO 1 J=1,NC
IF(I.GE.J) GO TO 1
EIJM(I,J)=EIJ(I,J)*(EI(I)*EI(J))**.5
GIJM(I,J)=GIJ(I,J)*(GI(I)+GI(J))/2.
1 CONTINUE
KM=0D0
UM=0D0
KM=0D0
UM=0D0
GM=0D0
QM=0D0
FM=0D0
DO 3 I=1,NC
KM=KM+Y(I)*KI(I)**2.5
UM=UM+Y(I)*EI(I)**2.5
GM=GM+Y(I)*GI(I)
QM=QM+Y(I)*QI(I)
3 FM=FM+Y(I)**2*FI(I)
KM=KM*KM
UM=UM*UM