32
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
DO 5 I=1,NC-1
DO 5 J=I+1,NC
UM=UM+2.*Y(I)*Y(J)*(UIJ(I,J)**5-1D0)*(EI(I)*EI(J))**2.5
GM=GM+2.*Y(I)*Y(J)*(GIJ(I,J)-1D0)*(GI(I)+GI(J))
5 KM=KM+2.*Y(I)*Y(J)*(KIJ(I,J)**5-1D0)*(KI(I)*KI(J))**2.5
KM=KM**.6
UM=UM**.2
DO 7 N=1,13
B1(N)=0D0
DO 9 I=1,NC
9 В1(N)=B1(N)+Y(I)*Y(I)(GIJM(I,I)+ 1D0-GN(N))**GN(N)*
*(QI(I)*QI(I)+1D0-QN(N))**QN(N)*(FI(I)+1D0-FN(N))*FN(N)*
*EIJM(I,I)"UN(N)*KI(I)*KI(I)*KI(I)
DO 11 I=1,NC-1
DO 11 J=I+1,NC