1. 前言
在当今能源危机和经济危机的状况下,船东们迫切希望对现有船型进行优化来改善船舶阻力性能,降低日常的营运成本 [1]。在船体线型优化的领域,CFD手段的应用已经越来越广泛,多采用势流理论与经验公式相结合的方法,既节省了人力物力,同时又对优化过程中各参数之于阻力性能的影响有了更为深刻地理解和认识。而基于母型船的参数化方法 [2] 则是目前较为先进的方法之一,通过应变能最小原理生成满足诸如面积、形心、曲线的起点或终点的坐标、斜率、曲率等特定几何要求的、光顺的特征曲线,进而基于NURBS理论,利用蒙皮方法生成光顺的船型曲面。船型参数间的耦合对船舶阻力性能的影响非常复杂 [3],寻找合适的优化参数将能够快速得到较优的方案,提高设计效率。其中,浮心纵向位置作为重要的考量船型参数,决定了船体艏部的肥瘦程度,从而对兴波阻力产生较大的影响。Lackenby变换方法是船型优化中较为常用的一种方法,能够实现对于船体浮心纵向位置的变换,并且使设计船较好地保留母型船一些优秀的特征,从而得到性能较优的新船型。
本文即是在对现有母型船进行半参数化建模的基础上,利用FRIENDSHIP-Framework系统中Lackenby变换方法,并应用势流理论,探讨某艏部轮廓较为特殊的破冰调查船浮心纵向位置与兴波阻力的关系。较为直观地给出该船兴波阻力系数相对浮心纵向位置的变化规律曲线,发现此船型存在一个最佳的浮心纵向位置使得兴波阻力系数最优,为该船型的设计提供可以参考的技术手段和方法。
2. Lackenby变换法
Lackenby变换方法 [4] 是由Lackenby H于1950年提出的,它是基于母型船的横剖面面积曲线(Sectional Area Curve,SAC曲线)基础上的船型变换,通过形变函数(
)的表达,使母型船变换工作有了更进一步的拓展。通过给定设计船需要变化的船型参数(棱形系数Cp,浮心纵向位置Xcb等),应用变形函数对母型船进行插值计算,得到设计船的SAC曲线,进而得到满足设计条件的设计船。变换函数在船体的前部和后部均可应用,亦可控制平行中体长度以及进流段、去流段的棱形系数来达到设计要求。通俗地讲,该方法可以在给定Cp、Xcb等参数的情况下,通过改变排水量沿船长方向的分布状态,较好地计算产生出新的船型,改变原有船型的阻力性能,具有较高的实用性和实践性 [5]。
3. 数值模拟工况及目标
选取某艏部轮廓较为特殊的破冰调查船作为研究对象,模型如图1所示,船型有关参数如表1所示。数值模拟工况:温度为20℃,淡水,傅氏数为0.1947。优化目标:在设计吃水,服务航速下,尽可能地减小兴波阻力。
(a) 艏部(b) 艉部
Figure 1. The sketch of the ship
图1. 船体模型示意图
4. 计算过程
4.1. 原型方案
首先,应用FRIENDSHIP-Framework软件 [6] [7] 进行静水力计算,得到设计吃水下船体的SAC曲线,如图2所示。
Figure 2. The SAC curve of original type
图2. 原型SAC曲线
4.2. 网格收敛性分析
在CFD数值模拟中,网格的收敛性分析必不可少 [8] [9] [10]。分别采用粗网格、中网格和细网格对原型的兴波阻力系数进行计算,计算结果如表2所示。由于粗网格尺度较大,计算结果与中网格和细网格的计算结果相差较为明显,而中网格和细网格的计算误差较小。综合考虑计算效率与精度,对改型方案均采用中网格进行计算。
Table 2. Comparison of wave resistance coefficient with different mesh scale
表2. 不同网格尺度下兴波阻力系数对比
4.3. 改型方案
在得到原型的SAC曲线之后,应用FRIENDSHIP-Framework系统中的Lackenby变化方法,对该船的前部进行调整,使后部线型基本不变,来确保兴波阻力为此次优化的主要目标。为保证SAC曲线能够光顺变化,对曲线变化段两端的切角进行限制,从而保证新的船体线型能够平顺过渡,提高生成的船体表面的质量。本次改型限制Cp不变,从而保证船型在变化前后排水量和湿表面积基本不变,Xcb的变化范围为±0.3% Lpp。以兴波阻力优化为目标,采用Sobol算法对40个改型方案进行优化计算。
4.4. 结果分析
图3给出了兴波阻力系数与Xcb之间的变化规律曲线,可以看出Xcb略微向后移动会改善此船的兴波阻力性能,且此变化范围存在一个临界值,超过此临界值继续后移则会对兴波产生不利影响,可能会导致前肩波影响不利。图4给出了静水力计算得出的最优SAC曲线与原型SAC曲线的对比,可看出优化后的船型水线面进水角减小,船体进流段变得尖瘦。新的SAC曲线变化平顺,再次说明了Lackenby变化方法的可靠性。图5给出了优化前后波形的对比情况,可看出优化后的船型艏部纵波与艉部横波均
Figure 3. The variation of wave resistance coefficient versus Delta_Xcb
图3. 兴波阻力系数相对浮心纵向位置变化趋势
有较大改善,说明船体前部变瘦之后产生的艏部波系与艉部波系产生了有利的干扰。图6中的优化前后波切高度对比图亦可看出船体前部产生的波峰和波谷在定量幅值上均有所降低,且后部波形相位与幅值均有所改变,带来了有利的影响。图7给出了优化前后的线型对比,从纵剖线图可看出优化后的线型S型更加明显,对耐波性同样会带来较好的收益。
Figure 4. The SAC curve of the optimized type
图4. 优化后的SAC曲线
Figure 5. The comparison of the wave pattern (up: opt, down: ori)
图5. 优化前后波形对比(上:opt,下:ori)
Figure 6. The comparison of the wave cut (red: opt, black: ori)
图6. 优化前后波切高度对比(红色:opt,黑色:ori)
(a) 站线(b) 纵剖线
Figure 7. The comparison of the shiplines (red: opt, black: ori)
图7. 优化前后船体线型对比(红色:opt,黑色:ori)
5. 结论
本文选取某艏部轮廓较为特殊的破冰调查船为研究对象,首先选取三种尺度网格对原型的兴波阻力进行计算,来验证网格的收敛性。然后应用FRIENDSHIP-Framework系统中的Lackenby变换方法对母型船进行变化,利用势流理论,选取中网格对改型方案进行计算,探讨Xcb的变化对兴波阻力的影响。
(1) Lackenby变化方法对于此种艏部较为特殊的破冰船适用性强,变化后的新船体曲面光顺程度高;
(2) 优化过程中,发现此种破冰船存在临界浮心纵向位置,使得在设计吃水、服务航速下兴波阻力最小;
(3) 浮心纵向位置得到优化后的船体,水线面进水角变小,艏部变得更加尖瘦;
(4) 优化后的船体艏部纵波和艉部横波都有较大的改善,艏部产生的波峰和波谷在幅值上均有降低;
综上所述,Lackenby变化方法对于此种破冰船具有较强的适用性,可以通过改种方法探讨此类船型浮心纵向位置对于船舶快速性能的影响。