弹性地基梁计算程序 Basic 源码

1 REM "G-2-2"
2 CLEAR : LPRINT "** ";
3 LPRINT "弹性地基梁计算程序(G-2-2)";
4 LPRINT " **"
5 OPEN "I",#1,"G-2-2.DAT"
40 INPUT #1, R,N,LM,L,E0,E,U,B
42 INPUT #1, NH: DIM X(NH),H(NH)
43 FOR I1=1 TO NH: INPUT #1, X(I1),H(I1): NEXT I1
45 DIM Z(R,4),AK(N+2),Y(N),D(N+2),B1(N+2)
47 DIM F((N+2)*(N+3)/2),X2(N*(N+1)/2)
50 FOR I=1 TO R: FOR J=1 TO 4: INPUT #1, Z(I,J): NEXT J: NEXT I
55 INPUT #1, RL: IF RL=0 THEN 65
60 DIM PL(RL),AL(RL): FOR I=1 TO RL: INPUT #1, PL(I),AL(I): NEXT I
65 INPUT #1, RR: IF RR=0 THEN 75
70 DIM PR(RR),AR(RR): FOR I=1 TO RR: INPUT #1, PR(I),AR(I): NEXT I
75 INPUT #1, U0
76 INPUT "需要打印原始数据吗 (Y/N)";G$
80 IF G$="N" OR G$="n" THEN 152
85 LPRINT "原始数据:"
90 LPRINT "地基的压缩模量 E0=";E0;"(KN/m)"
95 LPRINT "梁的弹性模量 E=";E;"(KN/m)"
100 LPRINT "梁的泊松比 U=";U; "梁宽 B=";B;"(m)"
105 IF NH=1 THEN 113
107 LPRINT : LPRINT TAB(2);"X(m)";TAB(11);"梁高H(m)"
108 X(0)=0 : FOR I=1 TO NH: LPRINT TAB(1);X(I-1);"-";X(I);TAB(14);H(I)
110 NEXT I: LPRINT : GOTO 115
113 LPRINT "梁高 H=";H;"(m)"
115 LPRINT TAB(1);"梁长L=";L;"(m)"
116 LPRINT TAB(1);"分段数N=";N;TAB(13);"荷载个数R=";R
120 LPRINT TAB(1);"问题类型(空间/平面/文克勒尔(1/2/3))LM=";LM
122 LPRINT "梁上荷载要素:"
125 LPRINT TAB(2);"Z(I,1)";TAB(11);"Z(I,2)";TAB(20);"Z(I,3)";TAB(29);"Z(I,4)"
127 FOR I=1 TO R: LPRINT TAB(1);Z(I,1);TAB(10);Z(I,2);TAB(19);Z(I,3);TAB(28);Z(I,4)
129 NEXT I: LPRINT
130 IF RL=0 THEN 135
131 LPRINT "边荷载要素:"
132 LPRINT TAB(2);"左PL(I)";TAB(11);"AL(I)"
133 FOR I=1 TO RL: LPRINT TAB(1);PL(I);TAB(10);AL(I):NEXT I: LPRINT
135 IF RR=0 THEN 140
136 LPRINT TAB(2);"右PR(I)";TAB(11);"AR(I)"
138 FOR I=1 TO RR: LPRINT TAB(1);PR(I);TAB(10);AR(I):NEXT I: LPRINT
140 LPRINT "U0=";U0;"(地基泊松比,文氏问题为梁下可压缩层厚度)
145 LPRINT "- - - - - - - - -":
152 IF LM<>2 THEN 165
154 INPUT "按平面应力/平面应变计算 (L/B)";A$
156 IF A$="L" OR A$="l" THEN LPRINT "按平面应力计算":GOTO 165
158 IF A$="B" OR A$="b" THEN LPRINT "按平面应变计算":GOTO 160
159 GOTO 154