read-subway.给学生.txt

地铁框架计算书

			README - SUBWAY                    
   地铁箱形框架计算程序：subfram.for，适用于明挖框架结构 。
   算例数据文件：subfram.dat--划分为67个单元。
   结构坐标：以结构最左上角节点为坐标原点，向右为x轴，向下为y轴。

   一、数据输入
     read(8,*) m,node,nfix,ig,ld
     m-total number of subfram elements.
     node-total number of subfram nodes.
    nfix-total number of subfram node restraince(=0).
    IG-total number of nodal points probably supported by springs.
    LD-number of loading groups.
       read(8,*) (ckkr(i),i=1,ig)
       read(8,*)(e(i),i=1,m)
   CKKR(M+1)-elastic resistance moduli of rock.
   E(M)-elastic moduli of lining elements.
       read(8,*) mblock
  mblock- 结构分段数,用于输出每根衬砌单元的i端、j端的节点编号。注意：
          分段总的原则--前后节点号差为等值者，且号差递增段为单独一段；号差递减段另为一段。
          具体解释为：前后节点号差为等值者，且凡偶数节点号集中的区域为一段；
                      前后节点号差为等值者，凡有奇、偶数节点相邻的区域为一段；
                      前后节点号差为等值者，凡奇数节点号集中的区域为一段；
                      一般来说，最小mblock＝6，且为偶数。
  mblock-number of the blocks divided (mblock=6,8,10...),can not single number.
       do 10 n=1,mblock
        read(8,*) kelem,noi,noj
  kelem-number of the elems of each block,noi,noj are node number of
        i point & j point.
        ...
10     continue
       read(8,*)  narea
 narea－结构分段数，与mblock分段原则不同，用途也不同:
        分段原则:
        ①截面厚度相同且轴线水平者为一段，截面厚度相同且轴线垂直者为另一段;
        ②角域变截面段应单独分段，亦按水平者为一段，垂直者为另一段。例如，将每个
          角域在水平方向分为两段，第一段为1个等厚单元，则第二段为2个变截面单元
          (有的角域正好相反，即先是2个变截面单元，然后才是1个等厚单元)。
          垂直方向划分方法一样。这样每个角域都分别按水平和垂直两个方向分开划分，共计6个单元。
        用途: 用于计算节点坐标值、单元截面积a(i)、单元轴线长度clg(i)、单元转角al(i)。
 narea-total number of the top table, down table,left,rihgt side wall(narea=5).
 xi,yi,xj,yj,ai,aj-X,Y coordinate and section area .
       do 20 i=1,narea
        read(8,*)bwxx,kelem,xi,yi,xj,yj,ai,aj
         do 30 ifix=1,nfix
         read(8,*) nofix(ifix),(ipres(nofix(ifix),j),j=1,3)
 nofix(ifix),(ipres(nofix(ifix),j),j=1,3)-example:16   0  1  0 ,
     in the 16th node  level displacement =0, rotate =0.
         ...
 30      continue
       read(8,*) cg,cgl,rg1,rg2,b0,h1,h2,h3,dd,ql,pt,pn,fai,posa,rsui
 cg-concrete's GAMA.
 cgL-materal's GAMA the road surface (generaly ,concrete).
 rg2-soil's GAMA up the subwater.
 rg3-soil's GAMA under the subwater.
 b0-the width of subfram.
 h1-the thickness of road surface.
 h2-the thickness of soil up subwater.
 h3-the thickness of under the subwater.
 dd-the thickness of top table of subfram.
 ql-road vehicle load(=10kpa).
 pt-the special load on top table(=150kpa).
    the special load on down table(120kpa=0.8*pt) and
    the special load on side wall (80kpa=0.533*pt) are treated in 
    program,need not reading.
 pn-train axial weight(=250kpa).
 fai-inter friction angle of soil.
 posa-折减系数,砂土=1.0,粘土=0.7,见教材p64
 rsui-水容重,1t/m3=10kpa=10kN/m3,在已有*.dat中均取为10kN/m3.

  算例文件：subfram.dat:

    67  66  2  21  1

    10
    25  10  60
     1  60  59
    29  59   1
     1   1   2
     4   2  10
     1  31  61
     2  61  65
     1  65  66
     2  66  62
     1  62  30


    33
     1     10.0   5.0   9.8  5.0  0.9  0.9
     2      9.8   5.0   9.3  5.0  0.9  0.4
     4      9.3   5.0   5.7  5.0  0.4  0.4
     2      5.7   5.0   5.2  5.0  0.4  0.9
     2      5.2   5.0   4.8  5.0  0.9  0.9
     2      4.8   5.0   4.3  5.0  0.9  0.4
     4      4.3   5.0   0.7  5.0  0.4  0.4
     2      0.7   5.0   0.2  5.0  0.4  0.9
     1      0.2   5.0   0.0  5.0  0.9  0.9

     1      0.0   5.0   0.0  4.8  0.9  0.9
     2      0.0   4.8   0.0  4.3  0.9  0.4
     4      0.0   4.3   0.0  0.7  0.4  0.4
     2      0.0   0.7   0.0  0.2  0.4  0.9
     1      0.0   0.2   0.0  0.0  0.9  0.9

     1     0.0    0.0   0.2  0.0  0.9  0.9
     2     0.2    0.0   0.7  0.0  0.9  0.4
     4     0.7    0.0   4.3  0.0  0.4  0.4
     2     4.3    0.0   4.8  0.0  0.4  0.9
     2     4.8    0.0   5.2  0.0  0.9  0.9
     2     5.2    0.0   5.7  0.0  0.9  0.4
     4     5.7    0.0   9.3  0.0  0.4  0.4
     2     9.3    0.0   9.8  0.0  0.4  0.9
     1     9.8    0.0  10.0  0.0  0.9  0.9

     1     10.0   0.0   10.0  0.2  0.9  0.9
     2     10.0   0.2   10.0  0.7  0.9  0.4
     4     10.0   0.7   10.0  4.3  0.4  0.4
     2     10.0   4.3   10.0  4.8  0.4  0.9
     1     10.0   4.8   10.0  5.0  0.9  0.9

     1     5.0    0.0   5.0  0.2  1.4  1.4
     2     5.0    0.2   5.0  0.7  1.4  0.4
     1     5.0    0.7   5.0  4.3  0.4  0.4
     2     5.0    4.3   5.0  4.8  0.4  1.4
     1     5.0    4.8   5.0  5.0  1.4  1.4

    21*0.9e5

    67*0.27e8

    65  0  1  0
    66  0  1  0

    25.0  25.0  20.0  16.0  10.0  0.5  2.0  8.0  0.5
    10.0  150.0  250.0  25.0


   二、 关于画图程序 TUPLOTG.FOR
      适用于FORTRAN POWER STATION(加油站)的版本。
 	 其原始数据与TUPLOT.FOR	完全相同。
	 注意：如果图形中座标轴X，Y字体不隹，那是因为原用的ROMAN字体不合
	 适，此时只要将语句：
	    fontnum = SETFONT ('t''roaman''h20w10B')
	 中的ROMAN 换成courier 就可以了，即改为：
	    fontnum = SETFONT ('t''courier''h20w10B')

   三、打印说明(亦可看FORTRAN 加油站的“帮助”的“PRINT”条款)
    ①用键盘的ALT+PRINT SCREEN 键(或PRINT SCREEN键)复制当前屏幕图形
    ②粘贴到WINDOWS的画图(在附件中)中
    ③右键用反色将兰色反为黄色
    ④在菜单的“图象”中点“属性”，用“黑白”将图象变为成黑白
    ⑤打印