利用MacCormack两部差分格式求解一维激波管问题Fortran程序

! MacCormack1D.for

!------------------------------------------------------------------------------------------------------------

!利用MacCormack两部差分格式求解一维激波管问题

!

!-----------------------------------------------------------------------------------------------------------

program MacCormack1D

implicit double precision (a-h,o-z)

parameter (M=1000)

common /G_def/ GAMMA,PI,J,JJ,dL,TT,Sf

dimension U(0:M+1,0:2),Uf(0:M+1,0:2)

dimension Ef(0:M+1,0:2)

GAMMA=1.4 !气体常数

PI=3.1415926

J=M !网格数

dL=2.0 !计算区域

TT=0.4 !总时间

Sf=0.8 !时间步长因子

call Init(U,dx)

T=0

1 dt=CFL(U,dx)

T=T+dt

write(*,*)'T=',T,'dt=',dt

call MacCormack_1D_Solver(U,Uf,Ef,dx,dt)

call bound(U,dx)

if(T.lt.TT)goto 1

call Output(U,dx)

end

!-------------------------------------------------------

!计算时间步长

!入口: U，当前物理量，dx, 网格宽度；

!返回: 时间步长。

!-------------------------------------------------------

double precision function CFL(U,dx)

implicit double precision (a-h,o-z)

common /G_def/ GAMMA,PI,J,JJ,dL,TT,Sf

dimension U(0:J+1,0:2)

dmaxvel=1e-10

do 10 i=1,J

uu=U(i,1)/U(i,0)

p=(GAMMA-1)*(U(i,2)-0.5*U(i,0)*uu*uu)

vel=dsqrt(GAMMA*p/U(i,0))+dabs(uu)

if(vel.gt.dmaxvel)dmaxvel=vel

10 continue

CFL=Sf*dx/dmaxvel

end

!-------------------------------------------------------

!初始化

!入口: 无；

!出口: U, 已经给定的初始值，dx，网格宽度。

!-------------------------------------------------------

subroutine Init(U,dx)

implicit double precision (a-h,o-z)

common /G_def/ GAMMA,PI,J,JJ,dL,TT,Sf

dimension U(0:J+1,0:2)

!初始条件

rou1=1.0

u1=0

v1=0

p1=1.0

rou2=0.125

u2=0

v2=0

p2=0.1

dx=dL/J

do 20 i=0,J/2

U(i,0)=rou1

U(i,1)=rou1*u1

U(i,2)=p1/(GAMMA-1)+0.5*rou1*u1*u1 20 continue

do 21 i=J/2+1,J+1

U(i,0)=rou2

U(i,1)=rou2*u2

U(i,2)=p2/(GAMMA-1)+0.5*rou2*u2*u2 21 continue

end

!-------------------------------------------------------

!边界条件

!入口: dx，网格宽度；

!出口: U，已经给定边界。

!-------------------------------------------------------

subroutine bound(U,dx)

implicit double precision (a-h,o-z)

common /G_def/ GAMMA,PI,J,JJ,dL,TT,Sf

dimension U(0:J+1,0:2)

!左边界

do 30 k=0,2

U(0,k)=U(1,k)

30 continue

!右边界

do 31 k=0,2

U(J+1,k)=U(J,k)

31 continue

end

!-------------------------------------------------------

!根据U计算E

!入口: U，当前U矢量；

!出口: E，计算得到的E矢量，

! U、E定义见Euler方程组。

!-------------------------------------------------------

subroutine U2E(U,E,is,in)

implicit double precision (a-h,o-z)

common /G_def/ GAMMA,PI,J,JJ,dL,TT,Sf

dimension U(0:J+1,0:2),E(0:J+1,0:2)

do 40 i=is,in

uu=U(i,1)/U(i,0)

p=(GAMMA-1)*(U(i,2)-0.5*U(i,1)*U(i,1)/U(i,0))

E(i,0)=U(i,1)

E(i,1)=U(i,0)*uu*uu+p

E(i,2)=(U(i,2)+p)*uu

40 continue

end

!-------------------------------------------------------

!一维差分格式求解器

!入口: U，上一时刻U矢量，

! Uf、Ef，临时变量，

! dx，网格宽度，dt,，时间步长；

!出口: U，计算得到得当前时刻U矢量。

!-------------------------------------------------------

subroutine MacCormack_1D_Solver(U,Uf,Ef,dx,dt)

implicit double precision (a-h,o-z)

common /G_def/ GAMMA,PI,J,JJ,dL,TT,Sf

dimension U(0:J+1,0:2),Uf(0:J+1,0:2)

dimension Ef(0:J+1,0:2)

r=dt/dx

dnu=0.25

do 60 i=1,J

do 60 k=0,2

!开关函数

q=dabs(dabs(U(i+1,0)-U(i,0))-dabs(U(i,0)-U(i-1,0)))/(dabs(U(i+1,0)-U(i,0))+dabs( U(i,0)-U(i-1,0))+1e-10)

!人工黏性项

Ef(i,k)=U(i,k)+0.5*dnu*q*(U(i+1,k)-2*U(i,k)+U(i-1,k))

60 continue

do 61 k=0,2

do 61 i=1,J

U(i,k)=Ef(i,k)

61 continue

call U2E(U,Ef,0,J+1)

do 63 i=0,J

do 63 k=0,2

!U(n+1/2)(i+1/2)

Uf(i,k)=U(i,k)-r*(Ef(i+1,k)-Ef(i,k))

63 continue

!E(n+1/2)(i+1/2)

call U2E(Uf,Ef,0,J)

do 64 i=1,J

do 64 k=0,2

!U(n+1)(i)

U(i,k)=0.5*(U(i,k)+Uf(i,k))-0.5*r*(Ef(i,k)-Ef(i-1,k))

64 continue

end

!-------------------------------------------------------

!输出结果, 用数据格式画图

!入口: U，当前时刻U矢量，

! dx，网格宽度；

!出口: 无。

!-------------------------------------------------------

subroutine Output(U,dx)

implicit double precision (a-h,o-z)

common /G_def/ GAMMA,PI,J,JJ,dL,TT,Sf

dimension U(0:J+1,0:2)

open(1,file='MacCormack_1D result.txt',status='unknown')

do 80 i=0,J+1

rou=U(i,0)

uu=U(i,1)/rou

p=(GAMMA-1)*(U(i,2)-0.5*U(i,0)*uu*uu)

write(1,81)i*dx,rou,uu,p,U(i,2)

80 continue

close(1)

81 format(D20.10,D20.10,D20.10,D20.10,D20.10)

end

ABAQUS中Fortran子程序调用方法

第一种方法： / o/ J5 @6 U/ ^- o$ 1. 建立工作目录/ ]" 2. 将Abaqus安装目录\6.4-pr11\site下的aba_param_dp.inc或aba_param_sp.inc拷贝到工作目录，并改名为aba_param.inc； # ~/ |0 I0 E6 {, @4 X3 q: W3. 将编译的fortran程序拷贝到工作目录； 4. 将.obj文件拷贝到工作目录； 5. 建立好输入文件.inp； 6. 运行abaqusjob=inp_name user=fortran name即可。 第二种方法： 在Job模块里，创建工作，在EditJob对话框中选择General选项卡，在Usersubroutine file中点击Select 按钮，从弹出对话框中选择你要调用的子程序文件（后缀为.for或.f）。 , D8 i7 d/r c6 @" | 以下是网上摘录的资料，供参考：. |$ t/ }$W7 Y6 m4 h6 D6 j 用户进行二次开发时，要在命令行窗口执行下面的命令： 4 O. R+ ^,@( ? abaqus job=job_name user=sub_name ABAQUS会把用户的源程序编译成obj文件，然后临时生成一个静态库standardU.lib和动态库standardU.dll，还有其它一些临时文件，而它的主程序（如standard.exe和explicit.exe等）则没有任何改变，由此看来ABAQUS是通过加载上述2个库文件来实现对用户程序的连接，而一旦运行结束则删除所有的临时文件。这种运行机制与ANSYS、LS-DYNA、marc等都不同。 : j6 g' R-o( {0 [* N2 J3 X这些生成的临时文件要到文件夹C:\Documentsand Settings\Administrator\Local Settings\Temp\中才能找到，这也是6楼所说的藏了一些工作吧，大家不妨试一下。 1子程序格式(程序后缀是.f; .f90; .for;.obj??) 答：我试过，.for格是应该是不可以的，至少6.2和6.3版本应该是不行，其他的没用过，没有发言权。在Abaqus中，运行abaqusj=jobname user=username时，默认的用户子程序后缀名是.for（.f,.f90应该都不行的，手册上也有讲过），只有在username.for文件没有找到的情况下，才会去搜索username.obj，如果两者都没有，就会报错误信息。 如果username包括扩展名for或obj，那么就根据各自的扩展名ABAQUS会自动选择进行操作。 2CAE中如何调用？Command下如何调用？ 答：CAE中在creat job的jobmanager中的general中可以指定子程序； Command下用命令:abaqus j=jobnameuser=userfilename （无后缀）； 3若有多个子程序同时存在，如何处理 答：将其写在一个文件中即可，然后用一个总的子程序调用（具体参见手册） 4我对VF不是很熟，是否可以用VC，C＋＋编写子程序？ A: 若要在vf中调试，那么应该根据需要把SITE文件夹中的ABA_PARAM_DP.INC（双精度）或ABA_PARAM_SP.INC（单精度）拷到相应的位置，并改名为ABA_PARAM.INC即可。 据说6.4的将可以，6.3的你可以尝试着将VC，C＋＋程序编译为obj文件，没试过。在你的工作目录下应该已经存在ufield.obj和uvarm.obj这两个文件（这两个文件应该是你分别单独调试ufield.FOR和uvarm.FOR时自动编译生成的，你可以将他们删掉试试看），但是由于你的FOR文件中已经有了UV ARM 和UFIELD这两个subroutine，显然会造成重复定义，请查实。 用户子程序的使用 假设你的输入文件为：a.inp b.for 那么在ABAQUS Command 中的命令应该是这样的： abaqusjob=a user=b

FORTRAN 90 程序编程规范

FORTRAN 90 程序编程规范 Fortran 90 编程规范，使程序代码高度组织化，更加易读、易懂、易于维护，程序更加高效。使编出的程序更易懂、易于维护。 1 语言选择 数值预报创新系统软件开发应避免使用Fortran77 的某些过时特征以Fortran 90不一致的特征。选择Fortran 90 作为开发语言，并采用Fortran 90 的新功能，如动态内存的分配(dynamic memory allocation)、递归（recursion ）, 模块(modules)、POINTER 、长变量名、自由格式等。 Fortran 77其中某些只是一些冗余的功能，这些功能已经过时，另外，还有一些在Fortran90 中被证明是不好的用法，建议不要使用。 2 Fortran 90 的新特性 2.1.1 建议使用的Fortran 90 新特性 建议使用Fortran 90 提供的模块（module ），并用Use ONLY 指定module 中哪些变量或派生类型定义可用于调用程序。 尽量使用数组下标三元组，这样可优化并减少所需的代码行数。为提高可读性，要在括号内表明数组的维数，例如： 1dArrayA(:) = 1dArrayB(:) + 1dArrayC(:) 2dArray(: , :) = scalar * Another2dArray(: , :) 当访问数组的子集时，例如在有限差分等式中，可以通过使用下标三元组实现。例如：2dArray(: , 2:len2) = scalar *（ & Another2dArray(:, 1:len2 -1) & - Another2dArray(:, 2:len2) & ) 对程序单元（program units ）命名，并使用End program ，End subroutine ，End interface ，End module 等结构再次指定“program unit ”的名称。 在逻辑表达式中使用>、 >=、 ==、 <、 <=、 /=，它们分别代 替.gt.、.ge.、.eq.、.lt.、.le.、.ne. 。新的表示方法更接近标准的数学符号 在变量定义中始终使用“：：”；始终用“DIMENSION ”定义数组形状；始终用（len=）的语法格式声明字符变量的长度。

FORTRAN程序设计复习题及答案

FORTRAN程序设计复习题 一、选择题 B (1)下列各FORTRAN表达式中合法的是 A) S+T*2P >= B) .NOT. (A*B+C) C) A2+B2/(C+D) <= D) (A+B).NOT．A*B．GT．(.NOT.只跟一个表达式) C (2)数学式(3/5)ex+y的FORTRAN表达式是 A) 3*EXP(X+Y)/5 B) 3*E* *(X+Y)/ C) (3/5)*EXP(X+Y)D) EXP(X+Y) D (3)下列FORTRAN77表达式中不合法的是 A) A．GT．B．EQV．C．GT．D B) A．AND．B．AND．C．AND．D C) .NOT．(X．LE．D) A．LT．B．LT．C．LT．D D(4)下列叙述中不正确的是 A) FORTRAN子程序可以单独编译 B) 对一个FORTRAN源程序进行编译和连接无误后可生成可执行文件 C) 即使编译和连接都正确无误，FORTRAN程序运行时仍可能出错 D) FORTRAN连接的主要任务是把函数库中的函数翻译成机器指令(正确描述：主要任务为连接目标文件) B (5)在下列FORTRAN77运算符中，优先级最高的是 A) .AND. B) .NOT. C) .OR. D) .EQ. B (6)FORTRAN表达式"6/5+9/2**3/2"的值为 A) 33 B) 1 C) 5 D) 3 A (7)下列FORTRAN77表达式中，合法的是： A) .AND.. B) 10.0 C) D) 提示：A）相当于 .AND.（.NOT.（）） D (8)关于编译一个FORTRAN源程序文件，下列说法中错误的是 A) 允许编译只有一个主程序而没有子程序的源文件 B) 允许编译有多个子程序的源文件 C) 允许编译只有一个子程序而没有主程序的源文件 D) 允许编译有多个主程序的源文件 C (9)在FORTRAN77源程序中，续行标志符必须放在 A) 第1列 B) 第1-6列C) 第6列D) 第5列 D (10)下列关于"SUBROUTIN E MAP(X，Y)"语句行的叙述中，不正确的是 A) 这是子程序的第一个语句 B) 字符串"MAP"是子程序名 C) 变量X是子程序的形参D) 子程序执行后，MAP将返回整型数据 提示：子程序无返回值，自定义函数才有） A (11)FORTRAN表达式"2/4+"的值是 A) B) 1 C) D) 0 提示：2/4默认等于整型，=》 D (12)FORTRAN表达式"MOD，"的值是 A) B)0.0 C) D) A (13下列FORTRAN运算符中，优先级最低的是 A）逻辑运算符.AND. B）算术运算符*

中国海洋大学Fortran程序考试典型题

1、找出100-150之间和400-450之间能被9整除的数 2、找出100-999之间的回文数 3、找出水仙花数

4、输出“*”图形 5、输入n个数，找出大于平均值的数和最小数

6、把一个数值型数组的相同数删除到只剩一个 7、形成一个5*5矩阵，对角线元素为“i”，其余为“j”program juzhen implicit none integer I,J integer,parameter::size=5 integer::a(size,size) forall(I=1:size,J=1:size,I>J) a(I,J)=j forall(I=1:size,J=1:size,I==J) a(I,J)=i forall(I=1:size,J=1:size,I

integer,allocatable::c(:) write(*,*) "输入数组A的数据个数" read(*,*) n write(*,*) "输入数组B的数据个数" read(*,*) m l=m+n allocate (a(n)) allocate (b(m)) allocate (c(l)) write(*,*) "从小到大输入A的元素" do i=1,n read(*,*) a(i) end do write(*,*) "从小到大输入B的元素" do i=1,m read(*,*) b(i) end do do i=1,n c(i)=a(i) end do i=1 do while(i=c(n))then n=n+1 c(n)=b(i) goto 10 else if(b(i)>c(j).and.b(i)

fortran程序实例

1)实例3—求多个半径下的圆周长 ! z3.f90 --Fortran95 ! FUNCTIONS: ! z3 - Entry point of console application. !************************************************************************* ! PROGRAM: z3 ! PURPOSE: Entry point for the console application. !************************************************************************ program z3 ! 求多个半径下的圆周长 ! 主程序 ! PROGRAM Z3 PRINT *, 'R=',1.2,'C=',C(1.2) PRINT *, 'R=',3.4,'C=',C(3.4) PRINT *, 'R=',15.6,'C=',C(15.6) PRINT *, 'R=',567.3,'C=',C(567.3) END program z3 !子程序 FUNCTION C(R) PI=3.1415926 C=2*PI*R RETURN ! Body of z3 end

2)实例4—键盘与显示器输入/输出 a)Fortran 基本操作 b)程序指令 ! ZXZ_I_O.f90 ! FUNCTIONS: ! ZXZ_I_O - Entry point of console application. ! PROGRAM: ZXZ_I_O ! PURPOSE: Entry point for the console application. !***************输入、输出样式种种************************** program ZXZ_I_O implicit none !变量声明的位置 INTEGER(2) i; INTEGER(4) j; INTEGER(4) m; REAL n INTEGER A,B ! Variables PRINT*,'输入整数A'; READ*, A PRINT*,'输入整数B'; READ*, B B=A+B PRINT*,'B=A+B=',B WRITE(*,*) 'A*B=',A*B PRINT* ,'以上为计算机的计算结果，注意B的值'

aspenfortran子程序的编译、连接、运行

Aspen plus fortran 子程序的编译、连接与运行 Aspen Plus可以用来模拟反应精馏，但提供的动力学方程形式比较死板，很多类型的动力学方程都不能直接输入，不得不外接一个fortran 子程序。本来用到fortran 的机会就不多，所以会的人也不多。何况编写好fortran 子程序后，还要用fortran 编译器来编译生成obj 文件，而且每次运行的时候都得连接一次，比较麻烦。我的电脑以前是装了编译器这些软件 的，但最近出了问题，只能编译fortran 子程序，运行aspen外接fortran 子程序的时候就会出错。恰好这几天要帮师兄写个fortran 子程序，写完后需要调试，但又运行不了，重装 呢也没把握，想起网上有兄弟说过可以用asplink命令生成一个dll文件，这样就不用运行 aspen的时候每次都连接编译器了。于是找来资料学习学习，居然成功了，欣喜之余，分享给大家。第一步当然是编写动力学子程序了（此处省略若干字） 第二步编译写好的动力学子程序（以网友soaty的动力学子程序为例） 打开Aspen Plus Simulation Engine JS M... 打开后是这样的 接下来改变文件夹路径，转到存放动力学子程序的那个文件夹 比如我的文件路径是C:\ProgramData\Aspe nTech\Aspen Plus \IPP 于是输入cd C:\ProgramData\AspenTech\Aspen Plus \IPP

回车后就转到这个文件夹了，这时候就可以进行编译了：编译命令：aspcomp 输入命令aspcomp ,回车 到文件夹一看你会发现已经生成了obj文件了,

几个fortran程序

计算圆周率 REAL R,R1,R2,PI ISEED=RTC() N0=0 N=300000 DO I=1,N R1=RAN(ISEED) R2=RAN(ISEED) R=SQRT(R1*R1+R2*R2) IF(R<1.0)N0=N0+1 END DO PI=4.0*N0/N WRITE(*,*)PI END 一）蒙特卡洛计算生日问题 假设有N个人在一起，各自的生日为365天之一，根据概率理论，与很多人的直觉相反，只需23个人便有大于50％的几率人群中至少有2个人生日相同。 INTEGER M(1:10000), NUMBER1(0:364), NUMBER2 REAL X,Y ISEED=RTC() DO J=1, 10000 NUMBER1=0 X=RAN(ISEED) NUMBER1(0)=INT(365*X+1) JJJ=1 DO I=1,364 Y=RAN(ISEED) NUMBER2=INT(365*Y+1) ETR=COUNT(NUMBER1.EQ.NUMBER2) IF (ETR= =1) THEN EXIT ELSE JJJ=JJJ+1 M(J)=JJJ NUMBER1(I)=NUMBER2 END IF END DO END DO DO I=1,10000 IF(M(I).LE.23) SUM=SUM+1 END DO PRINT *,SUM/10000 END 二）MONTE CARLO SIMULATION OF ONE DIMENSIONAL DIFFUSION 蒙特卡罗计算一维扩散问题

INTEGER X,XX(1:1000,1:1000) REAL XXM(1:1000) ! X:INSTANTANEOUS POSITION OF ATOM ! XX(J,I)：X*X ,J:第几天实验,I:第几步跳跃 ! XXM(I): THE MEAN OF XX WRITE(*,*) "实验天数JMAX，实验次数IMAX" READ(*,*) JMAX,IMAX ISEED=RTC() DO J=1,JMAX !第几天实验 X=0 !!! DO I=1,IMAX !第几步跳跃 RN=RAN(ISEED) IF(RN<0.5)THEN X=X+1 ELSE X=X-1 END IF XX(J,I)=X*X END DO END DO OPEN(1,FILE="C:\DIF1.DAT") DO I=1,IMAX XXM=0.0 XXM(I)=1.0*SUM(XX(1:JMAX,I))/JMAX !! WRITE(1,*) I, XXM(I) END DO CLOSE(1) END 三维的！ 三）通过该程序了解FORTRAN语言如何画图（通过像素画图）USE MSFLIB INTEGER XR，YR !在的区域中画一个圆 PARAMETER XR=400,YR=400 INTEGER R, S(1:XR,1:YR) X0=XR/2 ! 圆心位置X0，YO Y0=YR/2 R=MIN(X0-10,Y0-10) !圆半径 S=0 !像素的初始状态（颜色） DO I=1,XR DO J=1,YR IF((I-X0)**2+(J-Y0)**2<=R**2)S(I,J)=10 IER=SETCOLOR(S(I,J)) IER=SETPIXEL(I,J) END DO END DO END

fortran程序30个

程序1 PROGRAM chaper1272 INTEGER :: A(10)=(/5,7,4,8,12,2,10,3,9,11/) INTEGER :: sum=0,ave !打开一数据文件,设置一个有格式直接存取文件,将10个数分2个记录写入文件。 !数据文件生成2个记录,每个记录行长度相同,记录长度为25 OPEN(1,FILE='input21.dat',FORM='FORMATTED',ACCESS='D IRECT',RECL=25) WRITE(1,"(5I5)",REC=1)(A(I)+10,I=1,5) !按格式说明将头5个数写入第1个记录 WRITE(1,"(5I5)",REC=2)(A(I)+10,I=6,10)!按格式说明将后5个数写入第2个记录 READ(1,"(5I5)",REC=2)(A(I),I=6,10) !按格式说明从第2个记录中读取后5个数 READ(1,"(5I5)",REC=1)(A(I),I=1,5) !按格式说明从第1个记录中读取头5个数 DO I=1,10 sum=sum+A(I) ENDDO ave=sum/10 !打开一个最大记录长度为22的有格式顺序存取文件

OPEN(2,FILE='input22.dat',FORM='FORMATTED',ACCESS='D IRECT',RECL=22) WRITE(2,"('10个数之和为：',I5)",REC=1) sum !输出1记录行,记录长度为22 WRITE(2,"('10个数平均值为：',I5)",REC=2)ave !输出1记录行,记录长度为22 WRITE(2,"(A)",REC=3) '程序运行正常结束。' !输出一个记录行,记录长度为22 END 程序2 PROGRAM average REAL sum,ave INTEGER n OPEN(1,file='score.dat') PRINT*,'正在统计平均成绩,请等待。' sum=0.0;n=0 DO READ(1,*,END=100) S sum=sum+s n=n+1 ENDDO 100 ave=sum/n

完整word版ABAQUS中Fortran子程序调用方法—自己总结

第一种方法：在Job模块里，创建工作，在Edit Job对话框中选择General选项卡，在User subroutine file中点击Select按钮，从弹出对话框中选择你要调用的子程序文件（后缀为.for或.f）。 第二种方法： 1. 建立工作目录?崠 2. 将Abaqus安装目录\6.4-pr11\site下的aba_param_dp.inc 或aba_param_sp.inc 拷贝到工作目录，并改名为aba_param.inc； 3. 将编译的fortran程序拷贝到工作目录； 4.将.obj文件拷贝到工作目录； 5. 建立好输入文件.inp； 6. 运行abaqus job=inp_name user=fortran name即可。 以下是网上摘录的资料，供参考： 用户进行二次开发时，要在命令行窗口执行下面的命令： abaqus job=job_name user=sub_name ABAQUS会把用户的源程序编译成obj文件，然后临时生成一个静态库standardU.lib和动态库standardU.dll，还有其它一些临时文件，而它的主程序（如standard.exe和explicit.exe等）则没有任何改变，由此看来ABAQUS是通过加载上述2个库文件来实现对用户程序的连接，而一旦运行结束则删除所有的临时文件。这种运行机制与ANSYS、LS-DYNA、marc等都不同。 这些生成的临时文件要到文件夹C:\Documents and Settings\Administrator\Local Settings\Temp\中才能找到，这也是6楼所说的藏了一些工作吧，大家不妨试一下。1子程序格式(程序后缀是.f; .f90; .for;.obj??) 答：我试过，.for格是应该是不可以的，至少6.2和6.3版本应该是不行，其他的没用过，没有发言权。 在Abaqus中，运行abaqus j=jobname user=username时，默认的用户子程序后缀名是.for（.f,.f90应该都不行的，手册上也有讲过），只有在username.for文件没有找到的情况下，才会去搜索username.obj，如果两者都没有，就会报错误信息。如果username包括扩展名for或obj，那么就根据各自的扩展名ABAQUS会自动选择进行操作。 2CAE中如何调用？Command下如何调用？ 答：CAE中在creat job的job manager中的general中可以指定子程序；Command下用命令:abaqus j=jobname user=userfilename （无后缀）； 3若有多个子程序同时存在，如何处理 答：将其写在一个文件中即可，然后用一个总的子程序调用（具体参见手册）＋＋编写子程序？C，VC不是很熟，是否可以用VF我对4． A: 若要在vf中调试，那么应该根据需要把SITE文件夹中的 ABA_PARAM_DP.INC（双精度）或ABA_PARAM_SP.INC（单精度）拷到相应

Fortran95程序设计课后习题答案(word版方便)

第四章 1.program main implicit none write(*,*) "Have a good time." write(*,*) "That's not bad." write(*,*) '"Mary" isn''t my name.' end program 2.program main real, parameter :: PI=3 implicit none.14159 real radius write(*,*) "请输入半径长" read(*,*) radius write(*,"(' 面积='f8. 3)") radius*radius*PI end program 3.program main implicit none real grades write(*,*) "请输入成绩" read(*,*) grades write(*,"(' 调整后成绩为'f8.3)") SQRT(grades)*10.0 end program 4.integer a,b real ra,rb a=2 b=3 ra=2.0 rb=3.0 write(*,*) b/a ! 输出1, 因为使用整数计算, 小数部分会无条件舍去write(*,*) rb/ra ! 输出1.5 5.p rogram main implicit none type distance real meter, inch, cm end type type(distance) :: d write(*,*) "请输入长度:" read(*,*) d%meter d%cm = d%meter*100 d%inch = d%cm/2.54 write(*,"(f8.3'米='f8.3'厘米='f8.3'英寸')") d%meter, d%cm, d%inch end program 第五章 1.program main implicit none integer money real tax write(*,*) "请输入月收入" read(*,*) money if ( money<1000 ) then tax = 0.03 else if ( money<5000) then tax = 0.1 else tax = 0.15 end if write(*,"(' 税金为'I8)") nint(money*tax) end program 2.program main implicit none integer day character(len=20) :: tv write(*,*) "请输入星期几" read(*,*) day select case(day) case(1,4) tv = "新闻" case(2,5) tv = "电视剧" case(3,6) tv = "卡通" case(7) tv = "电影" case default write(*,*) "错误的输入" stop end select write(*,*) tv end program 3.program main implicit none integer age, money real tax write(*,*) "请输入年龄" read(*,*) age write(*,*) "请输入月收入" read(*,*) money if ( age<50 ) then if ( money<1000 ) then tax = 0.03 else if ( money<5000 )then tax = 0.10 else tax = 0.15 end if else if ( money<1000 ) then tax = 0.5 else if ( money<5000 )then tax = 0.7 else tax = 0.10 end if end if write(*,"(' 税金为'I8)") nint(money*tax) end program 4.program main implicit none integer year, days logical mod_4, mod_100, mod_400 write(*,*) "请输入年份" read(*,*) year mod_4 = ( MOD(year,4) == 0 ) mod_100 = ( MOD(year,100) == 0 ) mod_400 = ( MOD(year,400) == 0 ) if ( (mod_4 .NEQV. mod_100) .or. mod_400 ) then days = 366 else days = 365 end if write(*,"('这一年有'I3'天')") days stop end program 第六章 1.program main implicit none integer i do i=1,5 write(*,*) "Fortran" end do stop end program 2.program main implicit none integer i,sum sum = 0 do i=1,99,2 sum = sum+i end do write(*,*) sum stop end program 3.program main implicit none integer, parameter :: answer = 45 integer, parameter :: max = 5 integer weight, i do i=1,max write(*,*) "请输入体重" read(*,*) weight if ( weight==answer ) exit end do if ( i<=max ) then write(*,*) "猜对了" else write(*,*) "猜错了" end if stop end program 4.program main implicit none integer, parameter :: max=10 integer i real item real ans ans = 1.0 item = 1.0 do i=2,max item = item/real(i) ans = ans+item

邓肯张模型FORTRAN子程序源代码

邓肯张模型FORTRAN子程序源代码 SUBROUTINE UMA T(STRESS,STA TEV,DDSDDE,SSE,SPD,SCD, 1 RPL,DDSDDT,DRPLDE,DRPLDT,STRAN,DSTRAN, 2 TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,MA TERL,NDI,NSHR,NTENS, 3 NSTA TV,PROPS,NPROPS,COORDS,DROT,PNEWDT,CELENT, 4 DFGRD0,DFGRD1,NOEL,NPT,KSLAY,KSPT,KSTEP,KINC) C INCLUDE 'ABA_PARAM.INC' C CHARACTER*80 MA TERL DIMENSION STRESS(NTENS),STA TEV(NSTA TV), 1 DDSDDE(NTENS,NTENS),DDSDDT(NTENS),DRPLDE(NTENS), 2 STRAN(NTENS),DSTRAN(NTENS),TIME(2),PREDEF(1),DPRED(1), 3 PROPS(NPROPS),COORDS(3),DROT(3,3), 4 DFGRD0(3,3),DFGRD1(3,3) C DIMENSION PS(3),DSTRESS(NTENS),TDSTRESS(NTENS),TSTRESS(NTENS) PARAMETER (ONE=1.0D0,TWO=2.0D0,THREE=3.0D0,SIX=6.0D0) K=PROPS(1) N=PROPS(2) RF=PROPS(3) C=PROPS(4) FAI=PROPS(5)/180.0*3.1415926 G=PROPS(6) D=PROPS(7) F=PROPS(8) KUR=PROPS(9) PA=PROPS(10) DFAI=PROPS(11)/180.0*3.1415926 S1S3O=STA TEV(1) S3O=STA TEV(2) SSS=STA TEV(3) CALL GETPS(STRESS,PS,NTENS) FAI=FAI-DFAI*LOG10(S3O/PA) CALL GETEMOD(PS,K,N,RF,C,FAI,ENU,PA,KUR,EMOD,S,S3O,G,D,F 1 ,SSS,S1S3O) EBULK3=EMOD/(ONE-TWO*ENU) EG2=EMOD/(ONE+ENU) EG=EG2/TWO EG3=THREE*EG ELAM=(EBULK3-EG2)/THREE CALL GETDDSDDE(DDSDDE,NTENS,NDI,ELAM,EG2,EG)

Fortran求逆矩阵子程序

! 求逆矩阵 subroutine inverse(A,IA) implicit none real :: A(:,:), IA(:,:) real, allocatable :: B(:,:) integer :: i,j,N N = size(A,1) allocate(B(N,N)) ! 先把IA设定成单位矩阵 forall(i=1:N,j=1:N,i==j) IA(i,j)=1.0 forall(i=1:N,j=1:N,i/=j) IA(i,j)=0.0 ! 保存原先的矩阵A, 使用B来计算 B=A ! 把B化成对角线矩阵(除了对角线外,都为0) call Upper(B,IA,N) ! 先把B化成上三角矩阵call Lower(B,IA,N) ! 再把B化成下三角矩阵 ! 求解 forall(i=1:N) IA(i,:)=IA(i,:)/B(i,i) return end subroutine ! 求上三角矩阵的子程序 subroutine Upper(M,S,N) implicit none integer :: N

real :: M(N,N) real :: S(N,N) integer :: I,J real :: E do I=1,N-1 do J=I+1,N E=M(J,I)/M(I,I) M(J,I:N)=M(J,I:N)-M(I,I:N)*E S(J,:)=S(J,:)-S(I,:)*E end do end do return end subroutine Upper ! 求下三角矩阵的子程序subroutine Lower(M,S,N) implicit none integer :: N real :: M(N,N) real :: S(N,N) integer :: I,J real :: E do I=N,2,-1 do J=I-1,1,-1 E=M(J,I)/M(I,I)

Fortran程序总结

1.行的书写(行的长度、分行、续行) 一行可以是0~132个字符，空格有意义， 语句最长不超过2640个字符 一行可以有多个语句，用“;”分隔 一个语句可分行写，读行标记为&（放在尾部），但如为关键字，首尾均加&。 最多可有511个续行。 2.说明语句必须出现在可执行语句之前，格式说明语句（FORMAT语句）除外。 3.注释标志符： 自由格式：！固定格式：C * 语句分隔符：分号；（仅自由格式可以使用） 续行符：自由格式：& 申明标号：1到5位无符号整数 空格：关键字、变量、常量内部不能用空格，但相邻两者之间须用空格 4.信息处理的分类：数值处理和信息处理 现代计算机工作原理：程序存储和程序控制（冯·诺依曼） 1、运算器——算术运算、逻辑运算 2、控制器——根据指令控制计算机工作 运算器、控制器和寄存器称为中央处理器CPU 3、I/O设备——提供数据传输服务 4、总线——数据传输的公共通道 1.机器语言：二进制代码形式，可以被计算机直接执行，不可移植 2.汇编语言：用助记符来代替机器指令，容易记忆，不可移植 3.高级语言：接近自然语言（英语）的程序设计语言，又称算法语言，易学、易用、易 维护，可移植性好例：FORTRAN、BASIC、PASCAL、C、LISP、PROLOG等 5.FORTRAN90程序是一种分块结构，由若干个程序单元块组成： 主程序、外部子程序、模块、块数据单元无论是主程序单元，还是子程序单元，都是独立的程序单位，应该独立编写，它们的形式相似。 非语句行：注释语句：！后的所有字符都被编译器忽略。可独占一行，也可在其它语句之后，空行为注释行（固定格式用C和*） 6.常量的定义：常量是程序执行过程中不能变化的量。 基本数据类型有五种：整型、实型、复型、字符型和逻辑型 前三种属于数值类型，后两种为非数值类型，主要用于文字处理和控制。 FORTRAN 90通过KIND值确定整数的存储开销、最大位数和取值范围，如表所示。缺省KIND 值则取默认值4。 KIND值字节数二进位数取值范围 1 1 8 -128~127 2 2 16 -32768~32767 4 4 32 -2147483648~2147483647 格式：SELECTEND_INT_KIND(N) ： 功能：返回范围在-10N~10N的整数的kind-Value 在FORTRAN90中，也可以表示二进制、八进制、十六进制形式的无符号整数。其形式如下：二进制数：B’101101’ 或B”101101” 八进制数：O’76210’ 十六进制数：Z’1FA2’

平面四边形四节点等参单元Fortran源程序

C ************************************************ C * FINITE ELEMENT PROGRAM * C * FOR Two DIMENSIONAL ELASticity PROBLEM * C * WITH 4 NODE * C ************************************************ PROGRAM ELASTICITY character*32 dat,cch DIMENSION SK(80000),COOR(2,300),AE(4,11),MEL(5,200), & WG(4),JR(2,300),MA(600),R(600),iew(30),STRE(3,200) COMMON /CMN1/ NP,NE,NM,NR COMMON /CMN2/ N,MX,NH COMMON /CMN3/ RF(8),SKE(8,8),NN(8) WRITE(*,*)'PLEASE ENTER INPUT FILE NAME' READ(*,'(A)')DAT OPEN(4,FILE=dat,STATUS='OLD') OPEN(7,FILE='OUT',STATUS='UNKNOWN') READ(4,*)NP,NE,NM,NR WRITE(7,'(A,I6)')'NUMBER OF NODE---------------------NP=',np WRITE(7,'(A,I6)')'NUMBER OF ELEMENT------------------NE=',ne WRITE(7,'(A,I6)')'NUMBER OF MATERIAL-----------------NM=',nm WRITE(7,'(A,I6)')'NUMBER OF surporting---------------NC=',Nr CALL INPUT (JR,COOR,AE,MEL) CALL CBAND (MA,JR,MEL) DO I=1,NH SK(I)= enddo CALL SK0(SK,MEL,COOR,JR,MA,AE) do I=1,N R(I)= enddo pause 'aaa' stop READ(4,*)NCP,NBE,iz WRITE(*,'(5i8)')NCP,NBE,iz WRITE(7,'(5i8)')NCP,NBE,iz IF CONCR(NCP,R,JR) IF CALL BODYR(NBE,R,MEL,COOR,JR,AE) IF do jj=1,iz READ (4,*)Js,nse,(WG(I),I=1,4) read(4,*)(iew(m),m=1,nse) CALL FACER(iew,NSE,R,MEL,COOR,JR,WG) enddo endif CALL DECOP (SK,MA)

Fortran95实验指导 实验八 子程序

实验八 函数与子例行程序 子程序是FORTRAN95提供的重要功能之一，FORTRAN95允许用户自己定义子程序，允许用户通过子程序设计和编写程序。在设计和编写一个复杂的程序时，通常需要使用科学的程序设计方法，如：结构化程序设计方法、模块化程序设计方法、面向对象程序设计方法等，子程序功能完全支持结构化和模块化程序设计方法，结构化和模块化程序设计方法可通过子程序得到应用。使用子程序可提高程序的可读性、可理解性、可维护性和可修改性，是编写高质量、高水平、高效率程序的有效手段。 子程序分内部子程序和外部子程序。内部子程序是包含在程序单元CONTAINS 结构中的子程序，只有定义它们的程序单元才能调用执行它们。内部子程序包括函数子程序和子例行程序两类，通常将函数子程序简称为函数，子例行程序简称为子程序。 本次实验是学习和掌握子程序的一次系统实验活动。通过本次实验，能够熟练使用语句函数、内部子程序、外部子程序设计和编写程序，求解比较复杂的实际问题。 一、实验目的 1. 了解函数的有关基本概念，如：函数名称、函数自变量、函数调用、函数表达式等。 2. 掌握语句函数定义、调用、参数的基本内容和使用规则。 3. 掌握使用语句函数设计和编写程序的方法。 4. 了解内部子程序的有关基本概念，如：定义、调用、形式参数、实在参数、参数传递关系等。 5. 掌握内部子程序定义、调用、形式参数、实在参数、参数传递的基本内容和使用规则。 6. 掌握结构化程序设计方法。 7. 掌握使用内部子程序设计和编写程序的方法。 8. 了解递归概念，掌握递归问题的求解方法。 9. 掌握递归子程序定义、调用、参数传递关系的基本内容和使用规则。 10. 了解外部子程序概念，掌握并行、协调、合作编写大型复杂程序的方法。 11. 掌握外部子程序定义、调用、参数传递关系的基本内容和使用规则。 12. 掌握使用递归和外部子程序设计和编写程序的方法 二、实验内容1 1、问题描述 在解决某工程计算问题时，遇到下列函数： 2()12y x x x =++ 2 ()f x = ()p x = ()q x =计算x=0.1,0.2,0.3时，上述函数的值，并输出。使用语句函数编写程序实现之。