第11章程序设计343 000-25000 70.000-35.000 -40.000-27.000 54.000-31.000 70.000-31.000 52.000-24.000 10-70.000-24.000 计算成果 CIRCLEC: THE ABSCISSA VALUES OF THE UPPER AND LOWEI NO 51.522050 1.094647 50.955680 3-49822940-31.978910 -48.690200 30491210 863838 -47.557460 29278430 770531 46.424730 28265840 685123 7-45.29199027411070 605365 41.893780 25.566700 387709 11-40.761050 -25.147810 319849 1239628310-24.813490 13-38495570 24.558860 188260 14-37.362830-24.380470 123833 3509736024244130 003743 33.964630-24.284500 067425 1832831890 -131382 4.584810 -.19588 30.566420 4848610 261233 30.000050 5.000020 261233 11.3.3构造任意形状滑裂面的程序 SPLINE FOR 1.说明 本程序(刘德贵,1980)用来计算已知第二种边界条件(曲线两端点处的二阶导数值) 三次样条函数的一阶导数值,并利用样条函数对一元函数进行插值,进而构筑包含直线段的 复合滑裂面。对于滑裂面上控制点区间为直线段时,则直接线性插值,而对于非直线段区间 则运用上述介绍的方法进行三次样条插值。滑裂面被NS1个点分为NS1-1段。从上交点向 下交点编号为1,2,…NS1-1。程序用LS代表此NS1-1段中为直线段的线段总数。在数组 LSA(I)(=1,2,…LS)中存入这些直线段的编号 LS,滑裂面上NS1-1条线段中为直线段的总数 SA(LS),滑裂面上线段为直线段的编号;第 11 章 程序设计 343 7 8 2 5 9 3 9 10 3 3 1 2 3 10 1 .000 -25.000 2 -30.000 -25.000 3 -50.000 -35.000 4 -70.000 -35.000 5 -40.000 -27.000 6 -50.000 -29.000 7 -54.000 -31.000 8 -70.000 -31.000 9 -52.000 -24.000 10 -70.000 -24.000 计算成果 CIRCLE.C THE ABSCISSA VALUES OF THE UPPER AND LOWER NO. X Y ALF 1 -51.522050 -35.000000 1.094647 2 -50.955680 -33.901810 1.094647 3 -49.822940 -31.978910 .969015 4 -48.690200 -30.491210 .863838 5 -47.557460 -29.278430 .770531 6 -46.424730 -28.265840 .685123 7 -45.291990 -27.411070 .605365 8 -44.159260 -26.687410 .529825 9 -43.026520 -26.076840 .457513 10 -41.893780 -25.566700 .387709 11 -40.761050 -25.147810 .319849 12 -39.628310 -24.813490 .253490 13 -38.495570 -24.558860 .188260 14 -37.362830 -24.380470 .123833 15 -36.230100 -24.275990 .059923 16 -35.097360 -24.244130 -.003743 17 -33.964630 -24.284500 -.067425 18 -32.831890 -24.397580 -.131382 19 -31.699150 -24.584810 -.195887 20 -30.566420 -24.848610 -.261233 21 -30.000050 -25.000020 -.261233 11. 3. 3 构造任意形状滑裂面的程序SPLINE.FOR 1. 说明 本程序 刘德贵 1980 用来计算已知第二种边界条件(曲线两端点处的二阶导数值) 三次样条函数的一阶导数值 并利用样条函数对一元函数进行插值 进而构筑包含直线段的 复合滑裂面 对于滑裂面上控制点区间为直线段时 则直接线性插值 而对于非直线段区间 则运用上述介绍的方法进行三次样条插值 滑裂面被 NS1 个点分为 NS1−1 段 从上交点向 下交点编号为 1, 2, , NS1−1 程序用 LS 代表此 NS1−1 段中为直线段的线段总数 在数组 LSA(I)(I=1, 2, , LS)中存入这些直线段的编号 LS, 滑裂面上 NS1−1 条线段中为直线段的总数 LSA(LS) 滑裂面上线段为直线段的编号