土木工程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
120IF5130NUM=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(