对数螺旋线型双曲拱坝几何计算程序使用说明书.docxVIP

对数螺旋线型双曲拱坝几何计算程序使用说明书 对数螺旋线型双曲拱坝计算机辅助设计几何计算程序采用 QBASIC 语言编制，在一般微机上运行，该程序可解决对数螺旋线型双曲拱坝平面拱圈、各种横缝和孔口等的施工放样问题。 一、 坐标系及单位 1、三维直角坐标系的 Y 轴就是拱坝的“ 对称” 中心轴线，并指向下游； X 轴指向左岸； Z 轴垂直向下；座标系原点设在坝顶（一般在顶拱拱冠上游 点）。 2、 单位 程序输入、输出所用单位，长度以 m 计；角度以度计。 二、描述体型的主要参数及其函数关系 描述对数螺旋线型拱坝体型的主要几何参数有： 1、Yc Yc 是拱圈中心轴线在拱冠点处的 Y 座标值，或者说是拱冠梁中心轴线 上各点的 Y 座标。 2、Tc Tc 是拱冠梁各高程处的厚度 3、Tal 及 Tar Tal 及 Tar 分别是左、右两半拱拱圈的端部厚度。 4、Rl 及 Rr Rl 及 Rr 分别是左、右两半拱拱圈轴线在拱冠处的曲率半径。 5、θ l 及θ r θl 及θ r 分别是左、右两半拱拱圈轴（对数螺旋线）线方程中的初始角。 6、XDl 及 XDr XDl 及 XDr 分别是左、右两半拱拱圈下游端点 X 座标。 一般地说，上述参数都是 Z 座标的多项（ n+1 项）式： 在作施工放样座标计算时，上述全部参数的函数关系应尽知。 这些参数的函数式，其次数往往是不同的，设其中最高的次数是 n0 次， 用户在使用程序时，应把坝顶高程 H 和 n 的数值，库存在程序的第 21 行， o 0 前述各参数函数式中的系数 [A] 都要按序紧接 n0 库存，中间不允许插入任何别 的内容，而且， Tc 的系数 [A] 应从程序的第 23 行开始库存，每个参数的系数 都必须是 n0+1 个，不能多也不能少， 不足部分或未知者均须用若干个零按位 1 补足。 三、 主要计算公式 如图 1 示，某高程左右水平拱圈中轴线各为某对数螺旋线的一段，其极座标方程为： ek 相应参数方程为： xc 0 [ ek sin( ) sin ] y c Yc 0 [cos ek cos( )] 0 0 / 1 k 2 o e kφ 其中， k=tg θ，ρ =R ， R=R 式中： θ ：对数螺旋线的初始角； ρo：初始极半径； φ ：称为“ 似中心角”（拱中心角）； Ro ：拱轴线在拱冠处的曲率半径； R：轴线上任一点的曲率半径； Yc：拱轴线在拱冠处的 y 坐标； θ、φ 均以左曲线为正，右曲线为负。 拱圈厚度按下列公式变化 α T=T c+（TA—TC ）（Sc/SA） 式中， T：任一点处拱厚； TA：拱端厚度； Tc：拱冠梁厚度； Sc：拱轴线上一点至拱冠的弧长 Sc=Ro/k(e kφ-1)； SA：拱端至拱冠的拱轴线弧长； α：拱圈厚度变化指数，一般取 2.0 ～4.0 ，本工程取 4.0。 拱圈上、下游坝面曲线参数方程为： x = xc T sin 2 y y c T cos 2 曲线上一点 Q（Xq ，Yq ），其对应的曲率中心点是 OO （ XX， YY）， Sc 和 SA 以左半拱为正，右半拱为负。 设拱圈上有一条“ 平行” 于坝面（包括坝面）的曲线 S，若平行于上游 面，它到上游面的径向距离为 Lu；若平行于下游面，它到下游面的径向距离 为 Ld，规定， Lu 及 Ld 以在坝面内为正，反之为负。令： 2 Ls=Lu 或 Ld A=C{ （ 1+DD ）Tc/2-Ls} B=C ·DD·Tc/2 Lo=A-Bcos 1 （当 Ls=Lu 时） 系数： C= -1 （当 Ls=Ld 时） 则 S 的参数方程为 X=Xq+L 0 sin Y=Yq-L 0cos S 上一点（似中心角为 φ 处）到拱冠“ 弧距” SA（左正右负）用下式计 算： 1 Acos Bcoscos 2 2 2 2 SA 0 B cos sin d cos 当 DD=0 时， SA 的计算变得甚为简单： SA=Sc+C （Tc/2-Ls ） 拱冠到似中心角为 φ处的部分拱圈面积为： 0 A0= Tc( Sc DD (Sc 1 k 2 (ek (k cos sin ) k ))) 四、 拱坝横缝类型简介 拱坝横缝一般有四类，下面逐类简单介绍。 （一）一类曲面缝 一类曲面系全径向横缝面。指定某层（控制高程）拱圈中心轴上的一点 Xq， Yq），该层拱轴的通过该点的径向线，将与缝段区间内任意高程拱圈轴相交一点，通过交点的相应层的径向线在坝面内的线段就迭成一类曲面。 （二）二类曲面缝 二类曲面也是全径向横缝面，它是由通过驻点——控制点（ Xp， Yp）的区间内任意高程拱圈的径向线在坝面内的线段迭成。 （三）垂直平面缝 垂直平面缝通常形象地称为“ 一刀切” 面