[Fortran] 纯文本查看 复制代码
PROGRAM FSTIF
REAL A,B,X0,FK(4,4),XM,DENF
DIMENSION A1(4,4),A2(4,4),A3(4,4)
A=18.0 ！ 测试程序我自己赋的值
B=12.0 ！同上
X0=54.0 ！同上
DENF=0.035885 ！同上
XM=1.0 ！同上
DO 12 I=1,4
DO 12 J=1,4
A1(I,J)=0.
A2(I,J)=0.
12 A3(I,J)=0.
V1=X0*B/A/6.
A1(1,1)=2.*V1
A1(2,2)=2.*V1
A1(3,3)=2.*V1
A1(4,4)=2.*V1
A1(1,2)=-2.0*V1
A1(2,1)=-2.0*V1
A1(3,4)=-2.0*V1
A1(4,3)=-2.0*V1
A1(1,3)=-1.0*V1
A1(3,1)=-1.0*V1
A1(2,4)=-1.0*V1
A1(4,2)=-1.0*V1
A1(2,3)=V1
A1(3,2)=V1
A1(1,4)=V1
A1(4,1)=V1
2 V2=X0*A/B/6.
A2(1,1)=(2.-A/X0)*V2
A2(4,4)=(2.-A/X0)*V2
A2(1,3)=-V2
A2(3,1)=-V2
A2(2,4)=-V2
A2(4,2)=-V2
A2(2,2)=(2.+A/X0)*V2
A2(3,3)=(2.+A/X0)*V2
A2(1,4)=-(2.-A/X0)*V2
A2(4,1)=-(2.-A/X0)*V2
A2(2,3)=-(2.+A/X0)*V2
A2(3,2)=-(2.+A/X0)*V2
A2(1,2)=V2
A2(2,1)=V2
A2(3,4)=V2
A2(4,3)=V2
3 V3=B/A/A/12.
IF(A .EQ. X0) X0=X0+.001
E1=((A+X0)*(A+X0)*ALOG((X0+A)/(X0-A))-2.*A*(2.*A+X0))*V3
E2=((A-X0)*(A-X0)*ALOG((X0+A)/(X0-A))+2.*A*(2.*A-X0))*V3
E3=((A-X0)*(A+X0)*ALOG((X0+A)/(X0-A))+2.*A*X0)*V3
IF(A .EQ. X0) X0=X0-.001
A3(1,1)=2.*E1
A3(4,4)=2.*E1
A3(2,2)=2.*E2
A3(3,3)=2.*E2
A3(1,2)=2.*E3
A3(2,1)=2.*E3
A3(3,4)=2.*E3
A3(4,3)=2.*E3
A3(1,3)=E3
A3(3,1)=E3
A3(2,4)=E3
A3(4,2)=E3
A3(1,4)=E1
A3(4,1)=E1
A3(2,3)=E1
A3(3,2)=E2
DO 10 I=1,4
DO 10 J=1,4
10 FK(I,J)=3.14159*(A1(I,J)+A2(I,J)+A3(I,J)*XM*XM)/DENF
PRINT *, ((FK(I,J),J=1 ,4),I=1,4)
END
[Fortran] 纯文本查看 复制代码
SUBROUTINE FSTIF(A,B,XO,FK,XM,DENF)
!DEC$ attributes dllexport ::FSTIF
REAL A,B,X0,FK(4,4),XM,DENF
DIMENSION A1(4,4),A2(4,4),A3(4,4)
DO 12 I=1,4
DO 12 J=1,4
A1(I,J)=0.
A2(I,J)=0.
[Fortran] 纯文本查看 复制代码
PROGRAM SFTIFss
REAL A,B,X0,XM,DENF,FK(4,4)
DO 10 I=1,4
DO 10 J=1,4
10 FK(I,J)=0.0
A=18.0
B=12.0
X0=54.0
DENF=0.035885
XM=1.0
CALL FSTIF(A,B,X0,FK,XM,DENF)
PRINT *, ((FK(I,J),J=1,4),I=1,4)
END