土木工程basic.docx

土木工程basic.docx

《土木工程basic.docx》由会员分享，可在线阅读，更多相关《土木工程basic.docx（29页珍藏版）》请在冰豆网上搜索。

土木工程basic

《计算机在土木工程中的应用》

大作业

学院:

工学院

班级：

11级土木工程

（1）班

姓名：

学号：

时间：

2013年秋

《计算机在土木工程中的应用》课程考查题

一、数学相关编程

1、编写求解1!

+2!

+………+n!

的程序，并运行程序，输出结果。

2、结合习题2.6，编写一个辛普森法积分程序，并运行程序，输出结果。

二、测量学相关编程：

在角度平差计算中，由于测量产生的误差为：

，如果

时，各内角修正值为：

实测角度进过修正后为：

，试编关于角度平差计算程序，要求输出实测角编号I及角度

，修正值

，修正后角度

。

并以所学测量书上相关例题运行程序，并输出结果。

三、混凝土相关编程：

设检查混凝土强度测得n个强度值

，已知统计公式：

平均值

均方差

若规定某强度小于（

）为不合格，试设计一程序，当输入

后，可求出

、

及打印出不合格

的个数。

四、水力学相关编程：

明渠渠道断面计算课本计算程序为深度h为已知，求底宽b，试编写已知b，求深度h的计算程序。

五、平面桁架结构程序应用：

试编制程序计算如图所示桁架的内力，杆件截面均为A，弹性模量E=2.06×105N/mm2，设A=1.0

六、平面刚架结构程序应用：

已知如刚架结构，根据平面刚架内力分析程序写出（图二）的数据语句（data语句），横梁和柱的I=1.5E-2m4，A=0.3m2，E=1E6kN/m2。

一、1解题过程

10READN

20J=1

30SUM=0

40FORI=1TON

50J=I*J

60SUM=J+SUM

70NEXTI

80PRINT"DATA=";N

90PRINT"1!

+2!

+......+",I-1;"!

";"=";SUM

100DATA10

110END

运行结果

DATA=10

1!

+2!

+......+10!

=4037913

一、2解题过程

10READA,B,N,M

20AX=.0001

30S=0

40Y1=A*N+B

50YN=A*M+B

60JP=（M-N）/AX

70FORJ=1TOJP

80YJ=A*（N+J*AX）+B

82YJ1=A*（N+（J+1）*AX）+B

85SA=（YJ+YJ1）*AX/2

87S=SA+S

90NEXTJ

100PRINT

110PRINT"ALTERNATIONY=";A;"X+";B;",";"WIDEX=（";N;",";M;"）"

112PRINT

115PRINT"S=";S

120DATA3,4,0,9

130END

运行结果

ALTERNATIONY=3X+4,WIDEX=（0,9）

S=157.5036

2、解题过程

10READN

20SUM=0

25PRINT"DFM"

26PRINT"-------------------------------"

30FORI=1TON

40READA（I,1）,A（I,2）,A（I,3）

45PRINTA（I,1）,A（I,2）,A（I,3）

50A（I,1）=（A（I,1）+A（I,2）/60+A（I,3）/3600）

60SUM=SUM+A（I,1）

70NEXTI

80FL=（N-2）*180

90FB=SUM-FL

91F=（60/3600）*SQR（N）

100IFABS（FB）

110BI=-FB/N

120FORI=1TON

130AX（I）=A（I,1）+BI

140NEXTI

150PRINT"BDUBIZSZ（I）"

155PRINT"--------------------------------------------------------------------------------"

160FORI=1TON

170PRINTI,A（I,1）,BI,AX（I）

175NEXTI

180END

185DATA6

190DATA77,51,15

200DATA110,20,18

210DATA125,06,42

220DATA67,29,09

230DATA88,13,11

240DATA97,70,48

运行结果

DFM

--------------------------------------------------

775115

1102018

125642

67299

881311

977048

BDUBIZSZ（I）

-----------------------------------------------------------------------------------

177.8541625.46837103.3225

2110.338325.46837135.8067

3125.111725.46837150.58

467.4858325.4683792.95421

588.2197225.46837113.6881

698.1825.46837123.6484

3、解题过程

10DIMFCI（5）

20FCI

（1）=8.4:

FCI

（2）=9.5:

FCI（3）=10.3:

FCI（4）=11.5:

FCI（5）=12.8

30SUM=FCI

（1）+FCI

（2）+FCI（3）+FCI（4）+FCI（5）

40FCM=SUM/5

50PRINT"FCM=";FCM

60A=（FCI

（1）-FCM）^2+（FCI

（2）-FCM）^2+（FCI（3）-FCM）^2+（FCI（4）-FCM）^2+（FCI（5）-FCM）^2

70QC=SQR（A/4）

80PRINT"QC=";QC

90READN

100NUM=0

110FORI=1TON

120IF5

130NUM=NUM+1

140NEXTI

150PRINT"NUM=";NUM

160END

170DATA5

运行结果

FCM=10.5

QC=1.713184

NUM=1

四、解题过程

100READQ,M,N,J,B

105DATA40,1.5,0.025,0.0003,18.92

110D=.5:

H2=0:

E=.05

120H1=H2

130H=H1

135GOSUB350

140F1=F

150H2=H1+D

160GOSUB350

170F2=F

180IFF>0THEN120

200H=.5*（H1+H2）

210GOSUB350

220IFF>0THEN260

230F2=F

240H2=H

250GOTO300

260F1=F

270H1=H

300IFABS（H2-H1）

310H=.5*（H1+H2）

320GOSUB350

330IFF<=0THEN260

340GOTO230

350A=（B+M*H）*H

360K1=A^1.66667

370S=B+2*H*SQR（1+M^2）

380K2=N*S^.66667

390I2=SQR（J）

400K=Q/I2

410F=K-K1/K2

420RETURN

430H=.5*（H1+H2）

440PRINT"H=";H;"m"

800END

运行结果H=2.234375m

五、解题过程

30READNN,NE,NF

35N2=NN*2

40DIMSK（4,4）,DC（4,4）

45DIMU（N2）,SG（4）,SL（4）,X（NN）,Y（NN）

50DIMXA（NE）,XE（NE）,XI（NE）,YI（NE）

55DIMNS（NF,2）

60PRINT"NNNENF"

65PRINTNN;NE;NF

70PRINT"SUPPRESSEDDISPLACIMENTPOINT"

75FORI=1TONF

80READNS（I,1）,NS（I,2）

85PRINTNS（I,1）,NS（I,2）

90NEXTI

95FORI=1TONN

100READX（I）,Y（I）

110NEXTI

115PRINT"NO.X-COORD.Y-COORD"

120FORI=1TONN

125PRINTI,X（I）,Y（I）

130NEXTI

135PRINT"I","J","A","E"

140MX=0

145FOREL=1TONE

150READXI（EL）,YI（EL）,XA（EL）,XE（EL）

155PRINTXI（EL）,YI（EL）,XA（EL）,XE（EL）

156IDF=ABS（YI（EL）-XI（EL））

160IFIDF>MXTHENMX=IDF

165NEXTEL

170NW=（MX+1）*2

175NT=N2+NW

180DIMA（NT,NW）,Q（NT）,C（NT）

185FORI=1TONT

190FORJ=1TONW

195A（I,J）=0

200NEXTJ

205Q（I）=0

210NEXTI

215FOREL=1TONE

220I=XI（EL）:

J=YI（EL）:

AR=XA（EL）:

E=XE（EL）

225PRINT

226DXX=X（J）-X（I）

227DYY=Y（J）-Y（I）

228L=SQR（DXX*DXX+DYY*DYY）

235KS=（X（J）-X（I））/L

240SN=（Y（J）-Y（I））/L

245SK（1,1）=KS*KS

250SK（3,3）=KS*KS

255SK（2,1）=KS*SN

260SK（1,2）=KS*SN

265SK（3,4）=KS*SN

270SK（4,3）=KS*SN

275SK（1,3）=-KS^2

280SK（3,1）=SK（1,3）

285SK（1,4）=-KS*SN

290SK（4,1）=-KS*SN

295SK（2,3）=-KS*SN

300SK（3,2）=SK（2,3）

305SK（2,2）=SN*SN

310SK（4,4）=SN*SN

315SK（2,4）=-SN*SN

320SK（4,2）=-SN*SN

325CN=AR*E/L

330FORII=1TO4

335FORJJ=1TO4

340SK（II,JJ）=SK（II,JJ）*CN

345NEXTJJ

350NEXTII

355I1=2*I-2

360J1=2*J-2

365FORJJ=1TO2

370IFJJ=1THENNR=I1

375IFJJ=2THENNR=J1

380FORJ9=1TO2

385NR=NR+1:

II=（JJ-1）*2+J9

390FORKK=1TO2

395IFKK=1THENN9=I1

400IFKK=2THENN9=J1

405FORK=1TO2

410LL=（KK-1）*2+K

415NK=N9+K+1-NR

420IFNK<=0THEN430

425A（NR,NK）=A（NR,NK）+SK（II,LL）

430NEXTK

435NEXTKK

440NEXTJ9

445NEXTJJ

450NEXTEL

455FORII=1TONF

460N9=（NS（II,1）-1）*2+NS（II,2）

465A（N9,1）=A（N9,1）*10^11+10^11

470NEXTII

475FORII=1TON2

480READQ（II）

485C（II）=Q（II）

490NEXTII

495PRINT"VECTOROFNODALLOAD"

500PRINT"NO.X-LOAD.Y-LOAD"

510FORII=1TONN

515J2=2*II

520J1=J2-1

525PRINT"（";II;"）";Q（J1）,Q（J2）

530NEXTII

540PRINT:

PRINT"THENODALDISPLACEMENTARE"

545PRINT"NODEU:

V:

"

550GOSUB1000

560FORII=1TONN

565J2=2*II

570J1=J2-1

575PRINT"（";II;"）";TAB（15）;C（J1）;TAB（30）;C（J2）

580NEXTII

585FORII=1TON2

590U（II）=C（II）

595NEXTII

600FORII=1TO4

605FORJJ=1TO4

610DC（II,JJ）=0

615NEXTJJ

620NEXTII

625PRINT"FORCESINEACHMEMBER"

630PRINT"MEMBERFORCE"

635FOREL=1TONE

640AR=XA（EL）

645E=XE（EL）

650I=XI（EL）

655J=YI（EL）

670L=SQR（（X（J）-X（I））^2+（Y（I）-Y（J））^2）

675KS=（X（J）-X（I））/L

680SN=（Y（J）-Y（I））/L

685DC（1,1）=KS

690DC（2,2）=KS

695DC（3,3）=KS

700DC（4,4）=KS

705DC（1,2）=SN

710DC（3,4）=SN

715DC（2,1）=-SN

720DC（4,3）=-SN

725I1=2*I-2

730J1=2*J-2

735FORII=1TO2

740MM=I1+II

745NM=J1+II

750SG（II）=U（MM）

755SG（II+2）=U（NM）

760NEXTII

765FORII=1TO4

770SL（II）=0

775FORJJ=1TO4

780SL（II）=SL（II）+DC（II,JJ）*SG（JJ）

785NEXTJJ

790NEXTII

795FC=AR*E*（SL（3）-SL

（1））/L

800PRINTI;"-";J,FC

805NEXTEL

810GOTO1500

1000FORII=1TON2

1050IK=II

1055FORJJ=2TONW

1060IK=IK+1

1065CN=A（II,JJ）/A（II,1）

1070JK=0

1075FORKK=JJTONW

1080JK=JK+1

1085A（IK,JK）=A（IK,JK）-CN*A（II,KK）

1090NEXTKK

1095A（II,JJ）=CN

1110C（IK）=C（IK）-CN*C（II）

1115NEXTJJ

1120C（II）=C（II）/A（II,1）

1125NEXTII

1130FORIZ=2TON2

1135II=N2-IZ+1

1140FORKK=2TONW

1145JJ=II+KK-1

1150C（II）=C（II）-A（II,KK）*C（JJ）

1152NEXTKK

1155NEXTIZ

1160RETURN

1165DATA6,9,4

1170DATA3,1,3,2,4,1,4,2

1175DATA2,0,4,0,6,0,0,2,2,2,4,2

1180DATA1,2,1,1,1,4,1,1,1,5,1,1,2,3,1,1,2,5,1,1,2,6,1,1,4,5,1,1,5,6,1,1,3,6,1,1

1185DATA0,0,0,0,0,0,0,0,0,-10,0,-10

1500END

运行结果

（4）00

（5）0-10

（6）0-10

THENODALDISPLACEMENTARE

NODEU:

V:

（1）-37.0711-90.87827

（2）-18.04739-103.2353

（3）8.733809E-12-8.912619E-11

（4）-8.733813E-12-8.082971E-11

（5）-21.95263-109.902

（6）-42.92896-102.259

FORCESINEACHMEMBER

MEMBERFORCE

1-29.511852

1-413.45179

1-5-9.511848

2-39.023697

2-5-.6903573

2-6.4881592

4-5-10.97632

5-6-10.48816

3-6-14.8325

六、解题过程

80READNJ,MS,NF,NC

85N3=NJ*3

90PRINT"NJMSNFNC"

100PRINTTAB

（2）;NJ;TAB（6）;MS;TAB（11）;NF;TAB（16）;NC

110DIMU（N3）,K（6,6）,NS（NF,2）,I9（MS）,J9（MS）,RNS（NF）

120DIMDC（6,6）,SP（6）,SQ（6）,SG（6）,SL（6）,X（NJ）,Y（NJ）

130DIMS9（MS）,A9（MS）,U9（MS）,UP（MS）,BS（MS）,BX（2,6）,S2

（2）

140PRINT"NO.X-CORD.Y-CORD"

150FORI=1TONJ

160READX（I）,Y（I）

170PRINTTAB

（2）;I;TAB（6）;X（I）;TAB（16）;Y（I）

180NEXTI

190PRINT"POINTOFSUPPRESSEDDISPLACEMENT"

200PRINT"NO.DIRECTIONVALUES"

210FORI=1TONE

220READNS（I,1）,NS（I,2）,RNS（I）

230PRINTTAB

（2）;NS（I,1）;TAB（8）;NS（I,2）;TAB（18）;RNS（I）

240NEXTI

250READE,HY

260MX=0

270PRINT"MEMBERDETAILS"

280PRINT"NEI-J"

282PRINT"ASAWAWBBS"

290IFHY=0THEN380

300FORI=1TOMS

310READI9（I）,J9（I）,A9（I）,S9（I）,U9（I）,UP（I）,BS（I）

320II=I9（I）

330JJ=J9（I）

340IFABS（JJ-II）>MXTHENMX=ABS（JJ-II）

350PRINT"（";I;"）";I9（I）;"-";J9（I）

352PRINTA9（I）,S9（I）,U9（I）,UP（I）,BS（I）

360NEXTI

370GOTO430

380FORI=1TOMS

390READI9（I）,J9（I）,A9（I）,S9（I）,U9（I）

400PRINTI9（I）,J9（I）,A9（I）,S9（I）,U9（I）

402II=I9（I）:

JJ=J9（I）

410IFABS（JJ-II）>MXTHENMX=ABS（JJ-II）

420NEXTI

430NW=（MX+1）*3:

NT=N3*NW

440DIMA（NT,NW）,Q（NT）,C（NT）,JOD（NC,2）,ROAD（NC）

450FORI=1TONT

451FORJ=1TONW

452A（I,J）=0

453NEXTJ

454Q（