当前位置：文档之家 › FORTRAN程序设计大作业

FORTRAN程序设计大作业

  • 格式：doc
  • 大小：128.50 KB
  • 文档页数：6

1 上 海 师 范 大 学

FORTRAN 大作业

题目： 核检钢筋砼矩形截面梁的安全性能及

求算施工托架在荷载作用下的各杆轴力

专 业 土木工程（ 中英合作）

学 号 姓名 120445055 王璐祎

120445056 张双

指导老师 郭剑虹

完成日期 2014年11月

一．问题描述1： 2

现核检一钢筋混凝土矩形截面梁，已知该混凝土保护层为25mm（二a类环境），宽度为b，高度为h，所能承受弯矩设计值M=160kN·m，采用C20级混凝土，HRB400级钢筋，箍筋直径为8mm，截面配筋如下图1-1所示，请核查该核矩形截面的安全性。

图1-1

附表1：混凝土强度设计值（N/mm2） 3

附表2 普通钢筋强度设计值（N/mm2）

二．程序代码1：

PROGRAM EXAM

IMPLICIT NONE

INTEGER B,H

REAL FY,AS,H0,X,A1,FC(8),FT(8),CB,M,PMIN,P

READ *,B,H

A1=1.0 强度种类 符号 混凝土强度等级

C15 C20 C25 C30 C35 C40 C45 C50

轴心抗压 fc 7.2 9.6 11.9 14.3 16.7 19.1 21.1

23.1

轴心抗拉 ft 0.91 1.10 1.27 1.43 1.57 1.71 1.80 1.89

牌号 抗拉强度设计值fy 抗拉强度设计值fy’

HPB300 270 270

HRB335,HRBF335 300 300

HRB400,HRBF400,RRB400 360 360

HRB500,HRBF500 435 410 4 FC=(/7.2,9.6,11.9,14.3,16.7,19.1,21.1,23.1/)

FT=(/0.91,1.10,1.27,1.43,1.57,1.71,1.80,1.89/)

AS=1256

FY=360

CB=0.518

H0=H-25-8-10

PMIN=0.45*FT(8)/FY

P=AS/(B*H)

DO WHILE (P

选择语句当P

DO WHILE (P<0.002)

选择语句当P<0.002时

PRINT *,"FAIL"

输出 FAIL

END DO

结束选择语句

X=FY*AS/(A1*FC(8)*B)

那么X=FY乘以AS除以(A1乘以FC(8)乘以B)

IF(X<(CB*H)) THEN

选择语句当X小于CB乘以H时

M=FY*AS*(H0-X/2)

END IF

结束选择语句

IF (M>160) THEN

选择语句当M大于160时 5 PRINT *,"PASS"

输出“PASS”

ELSE

否则

PRINT *,"FAIL"

输出"FAIL"

END IF

结束选择

WRITE(*,10)FC

10 FORMAT(1X,2X4)

PAUSE

END

结束

三.问题描述2：

下所示为一施工托架的计算简图。在所示荷载作用下，试求各杆的轴力。

四．程序代码2：

Fya=Fyb=(8*4+6)/2=19

USE section A

19-8-Fyad=0

PRINT *,”Fyad=”,Fyad

READ *,Fyad 6 Fxad=Fyad*1.5/0.5

PRINT*,”Fxad=”,Fxad

READ *,Fxad

Fnac=-Fxad

PRINT*,”Fnac=”,Fnac

DO section C,D

PRINT*,”Fnce=”,Fnce;”Fncd=”,Fncd;”Fyde=”,Fyde;”Fnde=”,Fnde;”

Fxde=”,Fxde;”Fndf=”,Fndf

END DO

USE section E

IF Fye+Fyd+Fyf=0 THEN

PRINT*,Fnad,Fnac,Fnce,Fncd,Fnde,Fxde,Fndf

ELSE

DO section A,C,D

END DO

END IF

PAUSE

END

五．组员分配

张双 负责：提出结构力学课题，计算整体解题方案

王璐祎 负责：主程序的大致编程

王刚 负责：子程序、数组编程，电子打印，word制作

合集下载

fortran程序案例题汇编(14道)

fortran程序案例题汇编(14道)

1 / 18

1．Fibonacci数列定义如下：

F1=1F2=1Fn=F

n-1+F

n-2(n>2)

求Fibonacci数列的前30项。

integer f

(30),if(1)=1f(2)=2

do i=3,30

f(i)=f(i-1)+f(i-2)

enddo

print*,f

end

2．输入10个学生的总分，求每个学生的名次

integer s

(10),a

(10),i,j

do i=1,10

read*,s(i)

enddo

do i=1,10

a(i)=1

2 / 18

do j=1,10

if(s(i)

enddo

enddo

do i=1,10

print*,s(i),a(i)

enddo

end

3．给定一组数，按照从小到大的顺序输出。

integer a

(10)

integer p

do i=1,10

read *,a(i)

enddo

do j=1,9

p=j

do i=j+1,10

if (a(i)

enddo

if (p/=j) then

3 / 18

t=a(p);a(p)=a(j);a(j)=t

endif

enddo

print *,(a(i),i=1,10)

end

4．输入若干名学生的学号和三门课程（语数英）的成绩，要求从键盘输入一个学生的学号，能打印出该学生的三门功课成绩和总分。

character*6,dimension(:

),allocatable:

:xue

integer,dimension(:

,:),allocatable:

:g

integer,dimension(:

),allocatable:

:zong

integer i,j,n

character*6,xh

print *,"请输入学生的个数"

read *,n

allocate(xue(n))

allocate(g(n,3))

4 / 18

allocate(zong(n))

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

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.program 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

Fortran作业、期末复习、西南交大

Fortran作业、期末复习、西南交大

编写两个子程序，分别求两个整数的最大公约数和最小公倍数，并编写主程序调用这两个子程序，输出结果

program x

integer m,n,s

write(*,*)"请输入两个数:"

read(*,*)m,n

s=m*n

call acb(n,m)

call asd(n,m,s)

end

subroutine acb(n,m)

integer temp

if(m

temp=m;m=n;n=temp

end if

do while(mod(m,n)/=0)

temp=n;n=mod(m,n);m=temp

end do

write(*,"('最大公约数为：'i4)") n

return

end

subroutine asd(n,m,s)

if(m

temp=m;m=n;n=temp

end if

do while(mod(m,n)/=0)

temp=n;n=mod(m,n);m=temp end do write(*,"('最小公倍数为：'i4)") s/n

return

end

编写一个递归子程序，求Fibonacci数列某项的值。并编写主程序调用该递归子程序，输出结果。

program x

integer n,temp

write(*,*)"输入Fibonacci数列某项数,n大于三"

read(*,*)n

call abc(n)

end

subroutine abc(n)

integer i,j,temp,a(n)

a(1)=1

a(2)=1

do i=3,n

a(i)=a(i-1)+a(i-2)

end do

temp=a(n)

write(*,"('Fibonacci数列n项的值为：',i8)")temp

return

end

1：输入10名学生的学号和一门课程的考试成绩，分别存入顺序文件和直接文件中。每名学生的数据占一个记录。

分别对顺序文件和直接文件进行以下各项操作：

fortran作业

fortran作业

1、反射波走时计算：

program fanshebo

implicit none

real t,sum,t1,p,x

integer i,n

real,allocatable,dimension(:)::v,h !! Define arrays type

write(*,*)'imput Model strata n'

read(*,*) n

allocate(v(n),h(n)) !!For dynamic array allocate memory

write(*,102)'please imput',n,' level of velocity value'

read(*,*)(v(i),i=1,n)

write(*,102)'please imput ',n,' level of height value'

read(*,*)(h(i),i=1,n)

write(*,*)'please imput shot interval x'

read(*,*) x

sum=0.0

call tp(x,p,n,v,h) !!carry out intrinsic procedure

!! calculate the time T

do i=1,n

t1=h(i)/(v(i)*sqrt(1-(p**2)*(v(i)**2)))

sum=sum+t1

end do

t=2*sum

!!output the time T

write(*,101) 'Light through', n ,'layer needs time:'

write(*,100)t,'s'

100 format(1x,f10.4,1x,A)

101 format(1x,A,2x,i2,2x,A)

102 format(1x,A,1x,i2,1x,A)

deallocate(v,h)

contains