三体船型是目前船舶领域的一种热门船型,相对于常规船型,三体船型的兴波阻力更小,在恶劣的风浪条件下具有更高的稳定性,这些特点使得三体船的设计、制造成为船舶行业研究的热点,并在护卫舰、巡航舰、军民两用船等船舶类型上展开应用。我国在三体船型的研发和设计上处于世界先进水平,在三体船水动力理论、试验、船型设计等方面投入了大量的时间和物力。
本文的研究方向是基于模块化设计理念的三体船总体设计和优化,首先介绍模块化设计理论的现状,以某高速三体船为研究对象,采用NURBS曲线进行型线的拟合设计,基于型线设计进行三体船的特征曲面设计研究,最后结合流体动力学特性分析对基于模块化设计的三体船特性进行了建模和仿真测试。
1 模块化设计理论基本研究模块化设计理论最早是20世纪50年代欧美一些国家提出的,在产品大型化和快速化生产、装配过程中,基于模块化设计理论的产品不论是在生产效率、通用性和设计合理性方面,还是在产品的成本方面都具有显著的优势[1]。
模块化设计理论体系的基本原理如图1所示。
![]() |
图 1 模块化设计理论体系的基本原理 Fig. 1 Basic principles of modular design theory system |
模块化设计理论可以根据产品开发过程大致分为以下4个阶段:
1)模块化总体设计
总体设计是指在模块化设计的过程中对系统进行总体的策划,模块化总体设计的第1步是进行外部分析,整理来自外部的输入条件和需求,确定模块化总体设计的目标;第2步是进行模块化方案设计,在承接内部和外部的需求前提下,对方案进行决策。
2)模块化划分设计
模块化划分设计环节是模块化设计的核心环节,基于系统的总体功能特征,将系统进行不同功能和特性的分类和组合,结合系统的分析形成系统的宏模型,将其划分为若干个功能单元模块。在划分阶段,不能局限于功能本身,要考虑设计、制造不同阶段可能面临的问题,需要权衡利弊得到最优模块的划分解决方案。
3)模块化重组设计
模块化重组设计是指模块化创建与组合设计,模块化重组设计是模块化设计的第3个关键阶段,其目标在于将总体设计过程中产生的方案、子系统等内容进一步细化与重组,形成一个可以实际应用的方案。但产品的重组设计在设计过程中不强调每个具体的产品特性,而是强调大量产品所构成的产品系列特性。重组设计模块化设计的关键环节,主要流程包括确定子系统性能参数、几何尺寸、结构布局、模块接口、可靠性分析等,如图2所示。
![]() |
图 2 模块化重组设计的流程图 Fig. 2 Flowchart of modular restructuring design |
4)模块化产品设计
产品设计是模块化设计理论的落脚点和最终目标,也是模块化设计的最后一个环节。
2 基于NURBS的模块化三体船总体设计与优化 2.1 三体船的特征曲面设计由于三体船的外型结构和流体力学特性相对于传统船舶更加复杂,因此在三体船的特征曲面优化时需要综合考虑大量的特征参数和曲面建模,其中,型线参数如轮廓线、设计水线等,是决定三体船流体动力学特性的关键。
基于模块化设计思路,将三体船划分为不同的区域,每个区域的特征曲面按照由点到线、由线到面的设计思路,在型线设计时采用NURBS曲线进行型线的拟合设计
NURBS曲线模型用下式表示:
ˉp(m)=n∑i=0widiNi,j(u)n∑i=0wiNi,j(u)。 |
式中:
NURBS曲线的拟合过程用图3表示。
![]() |
图 3 NURBS曲线的拟合过程示意图 Fig. 3 Schematic diagram of NURBS curve fitting process |
基于模块化设计思路的三体船曲面设计分为下述4步:
1)确定特征点元素
根据三体舰船的区域特性,建立特征点控制方程如下式:
f1=minψ(x,y,z),f2=maxL(t),f3=minC(t),f4=Fx,y,z(t)。 |
其中:
2)由特征点建立三体船的特征型线
根据三体船的长度系数
fs(x)=1kxtgC(t)α0ψ(x,y,z)。 |
其中:k为曲线的曲率参数;
fs(x)=1Fx,y,z(t)∗kxtgL(t)α0。 |
3)三体船的型线优化
使用NURBS曲线函数对三体船的特征曲线进行优化,首先建立海水流场的物面条件:
∂ψ∂τ=0,∂ψ∂τ−Vn=0。 |
式中:
Fn(s)=k1vx+k2vy+k3vz√Hx2+Hy2+Hz2。 |
式中:
4)由三体船型线建立三体船的特征曲面
由线到面的设计过程借助三维软件SolidWorks完成,以特征点为坐标原点,在选定的型线上建立参考基准面,然后通过型线与基准面的平移、旋转、缩放、扫略等变换,生成三体船的特征曲面。
图4为基于SolidWorks产生的三体船特征曲面。
![]() |
图 4 基于SolidWorks产生的三体船特征曲面 Fig. 4 Three body ship feature surface based on SolidWorks |
在型线设计的基础上,对三体船的阻力和流体动力学性能进行数学建模和仿真分析,建立三体船的运动坐标系如图5所示。
![]() |
图 5 三体船的运动坐标系 Fig. 5 Trimaran motion coordinate system |
如图,运动坐标系选取速度方向为OX正向,OY与左舷方向重合,OZ垂向向上,在该坐标系下建立三体船的速度势
Δ2ψ = 0, |
建立三体船的运动模型为:
m∂(ρux)∂t+m∂(ρuxuy)∂t=m∂∂t[μ(∂ux∂y+∂uy∂x)]+Fi。 |
式中:
在流体特性分析时,海水和空气两种流体视为不相容,假设三体船附近的海水和空气流体计算域分别为
α(x,t)={0x∈V1,1x∈V2。 |
式中:
∂α∂t+u∂α∂x+v∂α∂y+w∂α∂z=0。 |
使用流体的VOF函数
∫S∂∂t(ρ0φ)dV+∫div(ρ0ˉφ)dV=S(V)。 |
式中:
使用仿真平台Matlab对三体船的流体动力学特性进行仿真,仿真的输入信号包括三体船的型线方程、运动模型和海浪波形函数。其中,输入的海浪波形函数如下:
ζ(t)=∑ni=1ξcos(kψ±ωt+εi)。 |
式中:
三体船仿真模型参数如表1所示。
![]() |
表 1 三体船仿真模型参数表 Tab.1 Parameter table of trimaran imulation model |
图6为三体船优化时设计前后的模型兴波阻力曲线对比。
![]() |
图 6 三体船优化时设计前后的模型兴波阻力曲线对比 Fig. 6 Comparison of wave making resistance curves of trimaran models before and after design |
三体船作为一种水动力特性更好的新型船舶,目前已经成为研究热点。本文结合模块化设计理念,针对三体船的总体重点利用NURBS曲线完成型线优化设计,基于型线设计完成了三体船的特征曲面设计,并对三体船进行流体动力学特性仿真分析。
[1] |
张业成, 詹成胜, 尚浩东. 基于四自由度的三体船操纵运动数值模拟研究[J/OL]. 武汉理工大学学报(交通科学与工程版): 1−13. ZHANG Ye-cheng, ZHAN Cheng-sheng, SHANG Hao-dong. Numerical simulation research on maneuvering motion of trimaran based on four degrees of freedom [J / OL]. Journal of Wuhan University of Technology (Traffic Science and Engineering Edition): 1−13 |
[2] |
田雷, 陈小星, 朱张立. 低噪声拖曳试验平台数值计算分析[J]. 机电工程技术, 2022, 51(6): 185-190. TIAN Lei, CHEN Xiao-xing, ZHU Zhang-li. Numerical calculation and analysis of low noise towing test platform[J]. Electromechanical Engineering Technology, 2022, 51(6): 185-190. DOI:10.3969/j.issn.1009-9492.2022.06.044 |
[3] |
孙小帅, 马骋, 钱正芳. 水面泵喷推进自航因子数值与试验研究[J]. 船舶力学, 2022, 26(5): 609-616. SUN Xiao-shuai, MA Cheng, QIAN Zheng-fang. Numerical and experimental study on self navigation factor of surface pump jet propulsion[J]. Ship Mechanics, 2022, 26(5): 609-616. DOI:10.3969/j.issn.1007-7294.2022.05.001 |
[4] |
石晓, 曹爱霞, 李丹. 三体船型阻力和运动性能分析[J]. 舰船科学技术, 2022, 44(9): 53-56. SHI Xiao, CAO Ai-xia, LI Dan. Analysis of resistance and motion performance of trimaran[J]. Ship Science and Technology, 2022, 44(9): 53-56. DOI:10.3404/j.issn.1672-7649.2022.09.011 |