船体外板水平截面展开方法研究

船体外板水平截面展开方法研究郭培军，李学军，唐兴玲，刘四平(中远船务工程集团有限公司技术中心，辽宁 大连116600)摘要应用空间解析几何坐标转换方法建立船体外板局部坐标系，选取外板展开准线，并由参数样条插值和参数样条曲线拟合方法得到准线方程以及各水平截面与外板的交线，从而求得水平截面展开曲线的拱高和交线的实长，最后根据拱高、交线实长和准线展开线得到外板平面展开形状，并为船体生产加工提供外板展开数据，建立了一种基于水平截面的船体外板展开方法。通过实船设计和生产建造的检验，该方法计算结果可以应用于实际生产。关键词外板展开水平截面展开法参数样条函数 坐标转换1引言各种船舶的外表面大都是由复杂的不可展的空间曲面构成，船体外板根据其弯曲的形状和程度，可以分为平直外板、单曲度外板和双曲度外板三种。在船舶生产设计过程中，需要将船体双曲度外板进行外板展开并根据展开形状数控套料。 传统的船体外板展开的主要依据是放样提供的已作了板缝排列的肋骨型线图(即包括了船体外板在肋骨型线图上的投影 )。展开的方法有十字线法、测地线法、短程线法等，它们都是基于手工作图产生的展开方法 [1卜⑶，其展开形状是一个封闭的平面图形。目前，国内各大船厂普遍采用了三维船舶设计软件如 Tribon软件，Tribon中应用等轴变换等方法对船体外板进行展开 ，但并不输出直接生产用的数据文件⑷。为了给生产加工提供如舵外板、艏艉柱包板等类型的船体外板展开数据，本文应用空间解析几何坐标转换、参数样条插值和参数样条曲线拟合等方法，建立局部坐标系和外板展开准线方程，并求得水平截面展开曲线的拱高以及截面与外板交线的实长，进而得到外板的平面展开形状和外板展开数据。最后将该本文所建立的方法应用于实船生产设计之中。2基本原理2.1坐标转换提取全船坐标系xyz下外板的型值，为了将其展开，建立外板的局部坐标系 uvw。图1为舵外板的局部坐标系，其坐标原点取作挂舵臂外板上端向艏中点， u轴取舵水平剖面中心线向艉，取垂直水平面向下为w轴方向，通过右手定则确定v轴方向。应用空间解析几何坐标转换法将全船坐标系 xyz下的外板型值转换到局部坐标系uvw下：_MU_MUIMv=Mx-OxMy-Oy〕Mw一UiMz-Oz"ujVi wiIvjwjVk Wk(1)式中，(MuMvMw)T表示局部坐标系下舵外板型值阵； (Mx-OxMy-OyMz-Oz)表示位移阵;uivWi示位移阵;ujvjWj表示局部坐标系的 u轴、v轴和w轴在全船坐标系中的方向余弦阵;UkvkWk一(Mx,My,Mz)表示全船坐标系下舵外板型值；(Ox,Oy,Oz)表示局部坐标系原点在全船坐标系中的坐标值。22)斜率连续 P』(0)=dp'⑴仝P'⑴ (5)图1舵的局部坐标系图1舵的局部坐标系图2曲线段端点的点矢和切矢示意图2.2参数三次样条函数插值原理⑸岡2.2.1参数三次曲线段2.2.1参数三次曲线段对于任一空间曲线,其参数三次曲线段的表达式为P(u)=AUT式中：U为参数，u•[0,1]；P(u)二X(u)y(u)z(u)l，称为点矢量，其导数 式中：U为参数，u•[0,1]；P(u)二X(u)y(u)z(u)l，称为点矢量，其导数 P'(u)称为切矢量;A-ax2ax1ax0ay2ay1ay0az2az1az0，为待定系数阵;ax3ay3灵3给定曲线段首、末端的点矢P(0),P(1)和切矢U=u3p'(0),p'(1)(如图2)，则可以求得Ferguson参数三次曲线段的具体表达式为:P(u)=|uP(u)=|u3u2■2-21们_P(0)1u11-3030-2110P(1)|p'(0)J000一-P'(1)一(3)2.2.2合成参数三次样条曲线连续性条件将若干参数三次曲线段拼接成一条完整的曲线，称为合成参数三次样条曲线。合成曲线内相邻曲线段间应满足如下连续性条件：1)位置连续P(1)=P+(0)1)位置连续P(1)=P+(0)(4)其中，Pi'(1)=:iti(l)；pi'(0^:i.iti1(0)；:-i和：•id分别为两曲线段在拼接处切失的模长；t二ti(1)=ti1(0)为两曲线段拼接处的单位切矢。3)曲率连续 pi4i"(0)=巴2p"(1)+HR'(1) (6)上式几何意义可表述为：为使两曲线段 R(Ui)和Pr(Ui41)拼接时达到曲率连续，曲线段 Phh(Uhh)首端的二阶导矢应处在由曲线段 P(uJ末端二阶导矢Pi"(1)和一阶导矢Pi'(1)所张成的平面内。式中口为一常数。2.2.3合成参数三次样条曲线的构造给定型值点P(i=0,1,…，n)，我们可以用型值点处的切矢 Pi'表达合成Ferguson曲线，其连续性方程为：Pj-+4Pj+匕州=3(Pj州—Pj_)，j=1,2,…，n—1 (7)对于空间曲线，上式包含三组方程，分别对应于 x,y和z坐标。2.2.4边界条件取曲线两端为自由端边界条件，即 P0"=Pn"=°，则补充方程为：[4P0'+2R'=6(P-P。)丿(8)2巳』+4巳'=6(巳一巳_4)由式(7)和(8)构成方程组，并用追赶法求解 P'(i二0,1,…，n)。3外板水平截面展开算法3.1水平截面展开曲线拱高本文以挂舵臂外板展开为例说明水平截面展开方法。在 2.1所建立的uvw坐标系下，挂舵臂外板与uow坐标平面的交线为B0Bn，选取该交线为外板展开准线，并应用参数三次样条函数方法得到准线方程。挂舵臂外板对 uow坐标平面的投影为 AoAnB^o(如图3)，将AAn分为n等份，得等分点A1,A2,…，And，分别求取过点A1,A2,,An4且平行于u轴的直线，交准线BoB.于点B「B2,…，，

线段A1B1,A2B2,…,An4Bn^即为垂直于平面坐标系UOW截取挂舵臂外板所得到的一组水平截面线， 令丨0=A0B0,h=AB1，…，In=AnBn。每条截面线展开后将为一曲线， 称为水平截面展开曲线 (如图4)。以截面线li为例说明水平截面展开曲线MiBiNi拱高的求取方法。将截面线AiBi分为..m等份，等分点为E1,E2,…，EmJ，以各等分点及A点为圆心建立圆方程，分别与准线BoBn相切于F1,F2,…,Fm4,Fm点，则BiF^BiF?,…,BiFm4,AiFm的长度即为水平截面展开曲线MiBi'Ni的一组拱高，记作h1,h2,,hm_,hm。

A L*iri [dn图3挂舵臂外板投影图4挂舵臂外板展开图3.2图4挂舵臂外板展开图为了准确地将挂舵臂外板展开，还需要求得各水平截面与挂舵臂外板的交线以及交线的实长。截面与外板的交线数据可以根据挂舵臂外板的型值，应用参数样条函数插值的方法得到。下面来求取水平截面i与挂舵臂外板的交线i的实长。以Bi为原点，BiAi方向为ui轴，左舷方向为Vi轴，建立局部坐标系UiOiVi（如图5）。图中，Ei,E2,…，Em4为截面展开曲线拱高求取时，截面线 AiBi的等分点。求取过点El,E2,…，Em^Ai且垂直于Ui轴垂线，与Ui轴两侧曲线的交点分别为Ci,C2,…,Cm和Di,D2,…，Dm，应用参数样条函数拟合交线 i，进而求得交线上各曲线段OiG,CiC2,…，Cm：Cm的实长，记为Si,S2,…,Sm。同理求得OiDi,DiD2,…，Dm」Dm的实际长度，记为S'i,S‘2,…，"。Cm-■.iHi-4■-J1Ji■i昨-i—Elh-JFjE，i)■)1图5截面与挂舵臂外板交线3.3挂舵臂外板的展开 _建立如图4所示v'o'w'坐标系，在w'轴上求取点B'1，使得o'B；=B0B1，其中，B0B1为图3中准线上Bo点到Bi点之间的曲线段长度。同理，可依次在 W'轴上得到B'2,B'3,…,B'n点，则线段O'B'n即为准线展开线。通过3.1和3.2求得了水平截面展开曲线拱高以及交线上各对应点之间的曲线段实长，下面根据该结果将水平截面 i与挂舵臂外板的交线i展开到平面上。以B'i点为圆心，分别以hi,h2,…，hm4,hm为半径建立圆方程，与w'轴正向分别交于点Fi,F2,…，Fm4,Fm。求得过已，F?,…，Fm_,Fm点且垂直于W'轴的垂线，并建立以B'i点为圆心，Si为半径的圆方程，与过Fi点的垂线交于2点，取v'坐标值大者为C'i点，然后再建立以C'i为圆心，S2为半径的圆方程，与过F2点的垂线也交于2点，取V'坐标值大者为C'2点，以此类推，可以得到一组点C'i,C'2,…,C'm,同理可得到点D'l,D'2，…,D'm。应用参数三次样条函数将点C'm,C'm」,…，C1,B'i,D1,D'2,…，D'm拟合为曲线，即得到交线i的水平截面展开曲线MiBi，Ni，同理求得其它交线的截面展开曲线。 分别将点Mo,Mi,…,Mn和N°,Ni,…,Nn用参数三次样条函数拟合为曲线，最后得到挂舵臂外板的平面展开形状（如图 4）。3.4水平截面展开算法流程水平截面展开算法的实现流程如图 6所示:图6水平截面展开算法流程4水平截面展开算法的应用应用水平截面展开算法对中远船务工程集团有限公司技术中心设计的 8160kw推船的舵外板进行展开，并根据舵外板展开结果进行钢板套料。 由于篇幅所限，本文仅给出了该船挂舵臂外板的型值数据 （如表1）和外板展开数据（如表2）。表1挂舵臂外板型值数据（mm）截面z=5500z=5000z=4170z=3600z=3160数据xyxyxyxyxy1504.0308.7504.0299.8504.0284.6504.0273.8504.0265.32853.3289.7834.3280.0802.6263.7780.8252.1764.0243.131202.6232.51164.6222.81101.2209.11057.6198.91024.0191.041552.101498.801399.701334.401284.0051202.6-232.51164.6-222.81101.2-209.11057.6-198.91024.0-191.06853.3-289.7834.3-280.0802.6-263.7780.8-252.1764.0-243.17504.0-308.7504.0-299.8504.0-284.6504.0-273.8504.0-265.3表2挂舵臂外板展开数据（mm）展开曲线12345数据fvw'fvw'fvw'fvw'fvw'11142.6119.11083.0615.8983.91440.7915.82007.0863.32443.82798.279.4755.0578.3685.61406.7639.21975.5604.42414.23449.539.7418.8540.8384.51372.7359.51944.0340.52384.64000503.301338.701912.402354.95-449.539.7-418.8540.8-384.51372.7-359.51944.0-340.52384.66-798.279.4-755.0578.3-685.61406.7-639.21975.5-604.42414.27-1142.6119.1-1083.0615.8-983.91440.7-915.82007.0-863.32443.8上面的计算结果已应用于该船的设计与生产，效果良好。5结论1） 本文针对船舶舵外板、艏艉柱包板等类型的双曲度外板，建立了可以提供船体外板展开数据的水平截面展开法。2） 将该方法应用于8160kw推船的舵外板展开，其展开精度和提供的数据满足实际生产要求。3） 该方法可以作为三维船舶设计软件外板展开计算的补充和验证。参考文献1戴寅生，于金生.船体外板计算展开[M].北京：国防工业出版社，1985.2王勇毅，董守富，林哲.计算机辅助船体建造[M].北京：人民交通出版社，1995.3徐兆康.船体曲板测地线作图方法分析[J].