► 在CAD软件中将船体的模型以stl格式导出,本案例中即为deck面、hull和transform面。如果从CAD模型中导出的是一个整体,可以划分为deck面和hull面。
►运行PumpLinx软件,新建一个工程文件,界面如下:
; y' T- h3 C& y
https://pic1.zhimg.com/v2-0f856eb20616187294c96714acf3e968_b.jpg►选择界面左边的Mesh窗口命令(一共4个窗口选项,分别是Mesh,Model,Simulation和Result,分别代表各个步骤)。
►选择Import/Export Geometry or Grid命令,点击Import Surface From STLTriangulation File,选择deck.stl、transom.stl和hull.stl文件,如图所示:
' b5 b$ v' K! p9 _* I/ L: W* W4 Z
https://pic3.zhimg.com/v2-243d8a0555c64966480960de3643df5a_b.jpg2. 生成网格
► 在Mesh窗口下,选择Marine Template Mesher命令。
► 设置Hull Type为Displacement,其中该选项下Planing指快艇类的船型,Dispalcement指集装箱类的大船;
► Analysis Type保持Resistance/Powering设置不变;
► Mesh Size可以保持Normal不变,或设置为Fine,也可设置User Defined进行自定义设置,本案例可设置为Fine;
► 点开Domain Size左边的小三角,设置Size Reference为Maximum Dimension,设置船体模拟的外流域参数如下:
Front 0.5 表示船体外域在船头方向往前延伸0.5倍船长
Back 1.5 表示船体外域在船尾方向往后延伸1.5倍船长
Side 1 表示船体外域在船侧方向往船身外侧延伸1倍船长
Depth 0.8表示船体外域在船体下沉方向延伸0.8倍船长
Up 0.2 表示船体外域在船体上方方向延伸0.2倍船长
Half Domain yes 表示半船模型
► 设置Max Pitch Angle即最大俯仰角为2度(大船);(小船10度)
► 设置Max Heave及最大下沉量为1m;
► 在Geometric Entities窗口下,选择hull面和transform面;在Properties窗口下,点击Hull右侧的 + ,同理选择deck面,将deck面添加;
: o5 ]$ q8 N2 }( {; xhttps://pic3.zhimg.com/v2-65a0e6c88083f32d88bdfbf588db8072_b.jpg► 设置船体的上浮方向为(0,0,1);
► 设置船体前进方向为(1,0,0);
► 设置船体重心坐标为(2.985,0.002,-0.111);
► 设置水线位置为0.111。注意水线位置是相对于船体重心获得的,不是相对于原点获得的;0.111表示水线在重心上方0.111m处。
► 最终参数设置如下:
5 C% _' R, ?8 A$ U
https://pic3.zhimg.com/v2-34ec9fc978a320a2c74104b7cab2f11a_b.png► 在Properties窗口下,点击Build Marine Template Mesher命令。
7 n0 F, a) `7 v1 |+ H2 vhttps://pic3.zhimg.com/v2-3be293684811f0278412ea40db532bfa_b.jpg
. i4 O+ V% d4 o►在Geometric Entities窗口会出现Volumns/marine选项,即新生成的网格模型。总网格数量约为135万(如硬件条件允许,可以适当加密网格),部分网格截图如下:
9 p% R2 G% O+ V) h( A
https://pic4.zhimg.com/v2-837060d30071274b6246aa316a266453_b.png8 V3 L! j8 F1 M
https://pic4.zhimg.com/v2-fe92b34060f5c3dc7dc5affcb54d925b_b.jpg►点开marine左边的小三角,在boundaries下面合并hull面和subfeatures面,随后点击save命令选择保存路径进行模型保存。注意船头部分应尽量避免出现subfeatures选项。
3.模型设置及求解
►点击Model窗口,由于采用marine模板进行的网格划分,此时PumpLinx在Model窗口下已经自动调用:
Marine 船舶模板
Flow/Multiphase and Turbulence VOF多相流和湍流模型
Translation(1DOF)marineHeave 1平动自由度,计算船体的下沉量变化模板
Rotation(1DOF)marinePitch 1转动自由度,计算船体俯仰角度变化的模板。
# ]0 ~# [8 K; L+ v9 [https://pic1.zhimg.com/v2-d35f0451aabac81a8c36befaa48f11b0_b.png►在Model窗口下,点击Marine,在Properties窗口下设置Dynamic Option为Pitch/Heave,表示既要计算下沉方向的运动,又要计算俯仰方向的运动情况。如果只需计算其中一种自由度情况或只考虑静水情况,则可根据分析需求分别选择Pitch Only/Heave Only/No Dynamics;
►设置Propulsion Option为Prescribed Profile,设置Target Velocity为2m/s,即目标速度;设置Ramp Up Time为2s,即加速时间为2s;
注释:另外两种设置
►方法一:也可设置Propulsion Option为Veclocity Profile,设置Velocity Profile为一定值或是随时间变化的曲线;
►方法二:也可设置Propulsion Option为Propulsion Source即添加动力源模型,设置Propulsion Source为Force Vector/Propeller Model;
►以Force Vector为例,Propulsion ID为动力源数量,Position表示动力源位置坐标,Direction表示作用方向,Force Expression表示力的分布;
►以Propeller Model为例,Propulsion ID为动力源数量,Position为螺旋桨中心位置,Direction为作用方向,Propeller Model可设置为Hough-OrdwayModel/Uniform Thrust Model;
4 m3 H, I9 H! w$ D8 f! h( |
! s- T% g! L/ F7 ^. \8 ^; a─ 以Hough-Ordway Model为例,Propeller Diameter为螺旋桨直径,Hub Diameter为螺旋桨轴径,Hub Thickness为螺旋桨轴向高度,PropellerDirection设置为Right-handed即右手法则,Propeller RPM为螺旋桨转速,Kt Curve为推力系数曲线,Kq为转矩系数曲线。
─ 以Uniform Thrust Model(均匀推力模型)为例,Propeller Diameter为螺旋桨直径,Hub Diameter为螺旋桨轴径,Hub Thickness为螺旋桨轴向高度,Propeller RPM为螺旋桨转速,Kt Curve为推力系数曲线。
2 @9 H8 i6 ?3 R
►设置Body Mass为1648.42370284735kg,注意此处应为整船的质量;
►设置Moment of Inertia即转动惯量为6134.9871kg.m2;
►设置Fluid Property如下:
: B6 p% w+ V8 `2 k5 Y H$ F
https://pic2.zhimg.com/v2-874890e58948b5fbedba38b777d366f9_b.png►完成上述设置后,在Simulation窗口下,设置CFD模拟参数如下:
) m1 P, j1 N7 J9 E+ H- O) J$ khttps://pic1.zhimg.com/v2-b68d9462617f404cb6cf1e03e1cda204_b.png►其中模拟时间为35s,时间步数为3325,迭代步数为25,保存频率为5,点击start即可启动计算。
► 注意:如果要查看具体的模型设置,可将Marine设置为Advance mode查看相关的参数设置。
4.计算结果查看
►云图显示。在Geometric Entities窗口的Volumns部分,选择需要查看的流体域,在Results窗口下,选择对应的变量即可显示云图(如压力、速度、水体积分数、涡量等)。
& B, ~8 v6 U* A9 q' z9 ]https://pic3.zhimg.com/v2-f683b0cc94079a61d2fcc2b06abe3632_b.jpg$ R0 P L. l2 k7 d: P
https://pic4.zhimg.com/v2-48e4cb5ad1edcb2bdfaa59ac4d336fbf_b.jpg►查看整船云图。
* [/ Y% \& c# L" G W, z2 z7 t
https://pic3.zhimg.com/v2-cad389e01c2cb402a8734a9760ed30b6_b.jpg►船体阻力、下沉量、俯仰角查看。
/ z, p7 t" q! V0 Lhttps://pic4.zhimg.com/v2-fbbd4f231416bfef0ad2b5b197fed23b_b.jpg►如需查看具体的阻力、下沉量、俯仰角数据或需要对结果数据进行再处理,可以通过copy date命令将数据直接粘贴至Excel或其他数据处理软件。
$ r' q5 ?( @% }( X( G8 Y* @4 ~
https://pic3.zhimg.com/v2-d168c2031570f9f8c1acb54a52b6bdf6_b.png►创建截面云图。在Geometric Entities窗口下,点击Create a Cross-Section命令,在窗口下方新增Derived Surfaces选项,选择新创建的截面,可通过拖动滑块或设置截面位置的方式确定某一截面位置。
& F9 T l T$ J' B' Y+ Shttps://pic1.zhimg.com/v2-0791ec047214da222c6772fccd806450_b.png6 L9 w0 s' V$ Q( E; m
https://pic4.zhimg.com/v2-8c36f57e9b50fb3c3cc19b439427b293_b.jpg►动画制作。在瞬态计算开始前,设置保存频率为5,即每隔5时间步保存一个结果,当计算完成后,在当前目录下则会存储多个结果文件。在模型显示区域,将三维显示效果调至最佳位置,然后在Flie下面选择Save animation,弹出对话框,选择对应的所有结果文件,然后点击“打开”,在给定该动画的名称,点击保存即可生成后缀为.gif的动画文件。
3 @. r8 P" Y, S2 [https://pic3.zhimg.com/v2-0678cad9246d7d76ce028b4105a8208a_b.jpg5 x7 j# f! f" R( M2 y4 _5 Y# d( Z, O