数据18个 6*6*No，A； ∧∧??∧6*6* 6*6*B，R24，Fi2，R4， Fi，D2，D4B，D4， H0，Qn，E1，E2， Ml，K2，K4，K； ∧∧??∧6*6* 数据18个 6*6*No，A； ∧∧??∧6*6* 6*6*B，R24，Fi2，R4， Fi，D2，D4B，D4， H0，Qn，E，E2， Ml，K2，K4，E； ∧∧??∧6*6* 数据17个 6*6*No，A； ∧∧??∧6*6* 6*6*，Ld，H，D1，D2， D5，H0，Qn，E， E2，Qd，Ml，K， K3，K5，E； ∧∧??∧6*6* 数据18个 6*6*No，A； ∧∧??∧6*6* 6*6*Ld，R24，H，D1， D2，D5，H0，Qn， E1，E2，Qd，Ml， K1，K2，K5，E； ∧∧??＜6*6* 本断面含立轴及横轴长圆
(3)数据穿孔 使用ALGOL程序(TQ—16)机，数据用‘N’格式分两段穿孔，也可用机器允许的简化格式穿孔。使用FORTRAN程序(FELIX机)，数据用卡片穿孔，前两个数据穿一张卡片，其余每8个穿一张卡片，每个数据后穿一逗点。 2.上机操作
按机器要求的操作命令，先输入源程序，待编译结束，进行标准变量赋值(当需要时)，再输入算题数据，机器开始算题。最后输出计算成果，计算结束。 本程序在TQ—16机安排标准变量赋值功能为：
H00000=1多组数据自动输入，若H00000不赋1或H00000≠l，则无此功能。
H00001=q，q为衬砌组成构件(半跨顶拱、边墙、底板等)的等分计算段数，可取10或30，40，50等分段。为提高计算精度，大断面隧洞可取较多的计算段数。当H000001不赋值，或所赋值不等于10或30，40，50时，程序自动按等分20段执行。 H00002=1水头按均匀水头计算。
“本程序的FORTRAN程序，已有多组数据自动输入功能，分段数q值的改变，可用调换q值赋值语句卡片来实现。等水头计算功能，可用改变编译卡，使程序中的注解成为计算语句来实现。” 3.计算成果
程序列表输出各计算点的内力和位移，格式为：
单位及符号：
M——弯矩，以衬砌内壁受拉为正，t·m。 N——轴向力，拉为正，t。
Q——剪力，逆时针转动为正，t。
V——法向位移，以向外法线方向位移为正，m。 U——切向位移，顺外法线向右为正，m。 W——转角位移，逆时针为正，弧度。
输出是由衬砌顶部的中点开始，以半跨顶拱(或顶板)、边墙、侧拱(或底圆角)、半跨底板(或底拱)为序，每段按等分10小段11点输出结果。 4.算例
圆拱直墙形无压隧洞，受顶部垂直山岩压力7t/m，齐墙顶内水压力及衬砌自重等荷载，断面尺寸如附图1。
并已知：
围岩弹性抗力系数 K=5×104t/m3； 混凝土弹性模量 E=23×102t/m2； 混凝土容重 Wl=2.4t/m3；
题号 NO=198306(任选)； 断面型号 A=4。
附图1
采用TQ—16机ALGOL程序原始数据： 6*6*198306，4；∧∧??∧6*6*
6*6*3，6，60，0.6，0.6，0.6，0.6，6，7，0，0，0，2.4，5410，5410，5410，23510；∧∧??∧6*6*输出成果；
原 始 数 据 NO：198306 A=4
3.00 6.00 60.00 0.60 0.60 0.60 6.00 6.00 7.00 0.00 0.00 0.00 2.40 50000 50000 50000 2300000
计 算 结 果 MX M 弯矩 7.125 6.822 5.943 4.522 2.683 0.547 ﹣1.728 ﹣3.964 ﹣5.970 ﹣7.554 ﹣8.526 ﹣8.526 ﹣2.127 1.238 2.530 2.571 1.820 0.353 ﹣2.020 ﹣5.522 ﹣10.027 ﹣14.521 ﹣14.521 ﹣5.658 0.487 4.437 6.717 7.810 8.136 8.037 7.781 7.553 7.465 N 轴向力 ﹣13.210 ﹣13.457 ﹣14.186 ﹣15.363 ﹣16.929 ﹣18.808 ﹣20.907 ﹣23.120 ﹣25.336 ﹣27.439 ﹣29.315 ﹣25.984 ﹣26.848 ﹣27.712 ﹣28.576 ﹣29.440 ﹣30.304 ﹣31.168 ﹣32.032 ﹣32.896 ﹣33.760 ﹣34.624 ﹣6.322 ﹣6.322 ﹣6.322 ﹣6.322 ﹣6.322 ﹣6.322 ﹣6.322 ﹣6.322 ﹣6.322 ﹣6.322 ﹣6.322 Q 剪力 ﹣0.000 ﹣1.660 ﹣3.207 ﹣4.535 ﹣5.546 ﹣6.157 ﹣6.302 ﹣5.936 ﹣5.038 ﹣3.607 ﹣1.192 13.625 7.884 3.613 0.928 ﹣0.662 ﹣1.823 ﹣3.126 ﹣4.856 ﹣6.795 ﹣7.959 ﹣6.322 34.624 24.730 16.528 10.101 5.369 2.151 0.209 ﹣0.716 ﹣0.889 ﹣0.567 0.000 V 法向位移 ﹣0.0014320 ﹣0.0014131 ﹣0.0013576 ﹣0.0012689 ﹣0.0011525 ﹣0.0010155 ﹣0.0008661 ﹣0.0007133 ﹣0.0005659 ﹣0.0004318 ﹣0.0003176 0.0002053 0.0001822 0.0001383 0.0001037 0.0000902 0.0000985 0.0001221 0.0001482 0.0001558 0.0001147 ﹣0.0000137 0.0008591 0.0007541 0.0006363 0.0005192 0.0004114 0.0003180 0.0002414 0.0001825 0.0001410 0.0001164 0.0001083 U 切向位移 0.0000000 ﹣0.0001458 ﹣0.0002876 ﹣0.0004215 ﹣0.0005443 ﹣0.0006532 ﹣0.0007466 ﹣0.0008235 ﹣0.0008840 ﹣0.0009292 ﹣0.0009608 ﹣0.0009908 ﹣0.0009794 ﹣0.0009675 ﹣0.0009553 ﹣0.0009426 ﹣0.0009297 ﹣0.0009163 ﹣0.0009026 ﹣0.0008884 ﹣0.0008739 ﹣0.0008591 ﹣0.0000137 ﹣0.0000124 ﹣0.0000110 ﹣0.0000096 ﹣0.0000082 ﹣0.0000069 ﹣0.0000055 ﹣0.0000041 ﹣0.0000027 ﹣0.0000014 0.0000000 W 转角位移 0.0000000 ﹣0.0000615 ﹣0.0001178 ﹣0.0001640 ﹣0.0001958 ﹣0.0002102 ﹣0.0002050 ﹣0.0001800 ﹣0.0001362 ﹣0.0000766 ﹣0.0000056 ﹣0.0000056 0.0000674 0.0000708 0.0000415 0.0000034 ﹣0.0000293 ﹣0.0000460 ﹣0.0000351 0.0000181 0.0001299 0.0003090 0.0003090 0.0003803 0.0003975 0.0003785 0.0003373 0.0002840 0.0002259 0.0001672 0.0001098 0.0000543 0.0000000 顶拱 点号 0 1 2 3 4 5 6 7 8 9 10 10 11 12 13 14 15 16 17 18 19 20 底板 20 21 22 23 24 25 26 27 28 29 30 由计算成果正理绘制衬砌内力及抗力分布图，见附图2。 附图2 算例内力及抗力分布图
第二节 源程序
6*6* begin
real r1，r24，r4，fi1，fi2，fi，d1，d2，d3，d4，d5，d4b，e，ml，k1，k2，k3，k4，k5，h3，h，ld，ho，qn，e1，e2，qd，fo2，fn2，go3，gn3，eo4，en4，ds1，ds2，ds3，ds4，ds5，ab，ab3，ab4，hf，hf3，hf4，sf2，cf4，ho3，ho4，ho5，en2，eo2；
in teger m1，m2，m3，m4，m5，m，a，b，no，q，t，s；
switch sw1：=la1，la2，la3，la4，la5，la6，la7，la8，la9，la10，la11； m24：#read(0，‘N‘，no，a)；
#print(0，‘80s-，2/，28X，4Hyuan，4X，3Hshi， 4X，3Hshu，4X，2H
u，/，28X，24s*，2/，28X，3Hno=，I8，9X，2Ha=，I2，/‘，no，a，)；
m1∶=m2∶=m3∶=m4∶=m5∶=0；m∶=0.01745329；q∶=20； for s：10，20，30，40，50do if H00001=s then q∶=H00001； goto sw1[a]；
la1：#read(0，‘N‘，r24，d2，ho，qn，e1，e2，qd，mL，k2，k4，e)；
#print(0，‘8F10.2，3I10‘，r24，d2，ho，qn，e1，e2，qd，ml，k2，k4，e)； d4b∶=d4∶=d2；r4∶=r24；fi∶=fi2∶=90；fn2∶=e2； en4∶=fo2∶=(e1+e2)/2；ab∶=#abs(ho)；
ho4∶=if ab＞r24 then #sign(ho)*(ab﹣r24)else 0； goto Loo；
la2：#read(o，‘N‘，ld，fi，d4，d4b，ho，qn，e1，e2，ml，k4，e)；
#print(0，‘9F10.2，2I10‘，ld，fi，d4，d4b，ho，qn，e1，e2，ml，k4，e)； en4∶=e2；ho4∶=ho；r4∶=Ld/#sin(fi*m)； goto Loo
la3：#read(0，‘N‘，b，ld，h，fi，d3，d4b，d4，ho，qn，e1，e2，ml，k3，k4，e)；
#print(0，‘I10，11F10.2，3I10‘，b，ld，h，fi，d3，d4b，d4，ho，qn，e1，e2，ml，k3，k4，e)； m3∶=q；ds3：=h/m3；h3：=h；
la31：r4：=ld/#sin(fi*m)；ho3：=ho；gn3：=e2；
ab3：=#abs(ho3)；hf3：=#sign(ho3)；ab：=#abs(ho)；