第11章程序设计337 1.3。 则根据上述原理可知,i、j应满足下面条件: 图112國弧滑裂面与边坡相 图11.3具有n个折点的边坡线与圆弧滑裂面相交的关系 交的关系 R1+1≤R<R 找到了上、下交点所在的线段号、j之后,就可以根据该线段两端点的坐标值和圆弧的 圆心坐标值x0、地和半径值R,算得其交点的坐标值。对于点m(m=i,,计算公式为 (x-x0)2+(y-y0)2=R (11.3) (-ym)/(x-xm)= (114) 联立以上方程组,求解得: x(+r>,+sT+(+s)Ro-(, -T)'I (11.6) T (11.7) 此联立方程组有两个根,计算下交点时取“+号,计算上交点时取‘-’号 11.2.3构造任意形状滑裂面的步骤 如果要将若干点用光滑的曲线相连,可以采用样条函数予以实现,在STAB程序中,采 用了第二种边界条件的三次样条函数。 第二种边界条件三次样条函数曲线两端点的二阶导数值为已知值,即y12=y(x) yn2=y(xn)。在构造滑裂面时取y2=y(x1)=0,yn2=y(xn)=0,由此确定的样条称为第 11 章 程序设计 337 11.3 则根据上述原理可知 i j 应满足下面条件 图 11. 3 具有 n 个折点的边坡线与圆弧滑裂面相交的关系 图 11. 2 圆弧滑裂面与边坡相 交的关系 (11.1) Ri+1 ≤ R0 < Ri (11.2) R j ≤ R0 < R j+1 找到了上 下交点所在的线段号 i j 之后 就可以根据该线段两端点的坐标值和圆弧的 圆心坐标值 x0 y0和半径值 R 算得其交点的坐标值 对于点 m (m = i, j) 计算公式为 (11.3) 2 2 0 2 0 (x − x ) + ( y − y ) = R (11.4) ( y − ym )/(x − xm ) = s (11.5) ( )/( ) s = ym+1 − ym xm+1 − xm 联立以上方程组 求解得 [ (1 ) ( ) ] (1 ) 1 2 0 0 2 2 x0 sT s R sx T s x + ± + − − + = (11.6) (11.7) 0 T = sxm − ym + y 此联立方程组有两个根 计算下交点时取 + 号 计算上交点时取 − 号 11. 2. 3 构造任意形状滑裂面的步骤 如果要将若干点用光滑的曲线相连 可以采用样条函数予以实现 在 STAB 程序中 采 用了第二种边界条件的三次样条函数 第二种边界条件三次样条函数曲线两端点的二阶导数值为已知值 即 ( ) 12 1 y = y′′ x ( ) n2 n y = y′′ x 在构造滑裂面时取 y12 = y′′(x1 ) = 0 yn2 = y′′(xn ) = 0 由此确定的样条称为