作为一种新型水下航行器, 水下滑翔机因其航行时间长、航行距离远、监测范围广的特点成为了探测海洋环境的首选。水下滑翔机发展至今, 按照其外形可以分为常规布局水下滑翔机和非常规布局水下滑翔机。翼身融合水下滑翔机作为一种非常规的水下滑翔机, 采用机翼与机体平滑连接的方式, 增大了升力面积, 减小了表面阻力, 因此具有更好的升阻比特性[1]。但是这种独特的外形也带来了设计上的困难。随着计算流体力学(CFD)、最优化等学科的发展, 基于CFD的外形设计在工业上取得了极大进步[2]。
利用CFD开展外形优化设计首要的就是参数化方法。现有的参数化方法包括基于CAD的参数化、PARSEC参数化、CST参数化、B样条参数化、FFD参数化等[3]。对于翼身融合水下滑翔机而言, 基于CAD的参数化方法是将三维设计空间分割为一些二维截面, 通过对二维截面的优化, 得到优化后的三维几何外形[4]。PARSEC参数化、CST参数化[5]、B样条参数化等作为成熟的翼型参数化方法, 是将翼身融合水下滑翔机分割为几个翼型截面, 通过对翼型的优化实现对整机的优化。而FFD参数化方法通过对控制体的变形来操作所映射外形的变形, 很好地保证了所研究外形的光滑性和连续性, 且保留了相当好的局部变形能力, 适用于复杂的三维外形。近年来, 越来越多的学者在工程实际中选用这种参数化方法, 包括气动外形优化[6], 车辆外形优化[7]等, 均取得了较好的效果。
本文旨在利用FFD参数化方法, 实现对翼身融合水下滑翔机精细化变形; 配合CFD求解器, 优化算法及网格变形算法, 形成一套外形优化设计框架, 并对一种翼身融合水下滑翔机开展外形优化设计工作, 以提升其水动力性能。
1 外形优化设计框架本文搭建的外形优化设计框架如图 1所示[8], 运行过程如下:
![]() |
图 1 外形优化设计框架 |
·绘制几何外形初始网格及FFD控制体;
·使用FFD参数化方法建立FFD控制体与几何目标的映射关系;
·将网格文件导入CFD求解器中求解目标函数及其梯度;
·将计算结果传递至SQP算法, 确定FFD控制体的变形量, 间接变形几何外形;
·根据几何表面网格变形更新整个计算域网格;
·重复迭代输出最优外形;
下面主要介绍框架中的FFD参数化方法、网格变形算法及SQP优化算法。
1.1 FFD参数化方法FFD参数化方法通过建立FFD控制点与目标外形表面点之间的数学映射关系, 实现通过控制点的移动变形目标外形的目的, 从而实现目标外形的参数化[9]。其基本原理如下:
首先在FFD控制体上添加一个局部坐标系O′-STU, 任何点X在这个系统中具有(s, t, u)坐标, 所以在全局坐标系O-XYZ里, 任意点X具有如下形式
![]() |
(1) |
式中,X0为局部坐标系的原点在全局坐标系下的向量表示。
任意点X的局部坐标坐标(s, t, u)很容易用线性代数方式表达, 写成向量形式为
![]() |
(2) |
由(2)式可知, X的局部坐标(s, t, u)应满足
![]() |
接下来在FFD控制体上建立控制点Pi, j, k, Pi, j, k为FFD控制体上第(i, j, k)个控制点(i, j, k表示FFD控制体上分别沿着坐标轴U, V, W方向上的下标, 每个方向被分为l, m, n份)。则Pi, j, k可表示为
![]() |
(3) |
故而, FFD控制体空间内任一点的笛卡尔坐标X(s, t, u)可表示为
![]() |
(4) |
式中,Ni, p(s), Nj, q(t), Nk, r(u)为以s, t, u为变量的p, q, r次基函数。
在按照(4)式建立了初始模型与FFD控制体上控制点的映射关系后, 就可以通过改变FFD控制体上控制点Pi, j, k的位置来改变内部几何模型的形状。
1.2 网格变形算法在1.1节中, FFD参数化方法建立了控制点与几何外形表面网格节点之间的映射关系, 当控制点移动时, 表面网格节点也会发生变形, 为了对变形后的网格进行CFD求解, 需要对空间网格更新, 这就离不开网格变形。网格变形算法发展至今分成3大类:代数法、物理模型法和混合方法[10]。其中代数法运算速度快, 求解效率高, 易于编程实现。本文采用的网格变形算法为基于弧长的代数型方法, 其基本原理如下
![]() |
(6) |
式中:Xs是几何外形表面网格点的位移; Xv是空间网格节点的位移。二者通过弧长的归一化数F进行求解。F的求解以i方向为例, 如(7)式所示。
![]() |
(7) |
这样就可以实现根据几何表面网格节点的位移实现更新空间网格的目的。
1.3 SQP优化算法序列二次规划法(SQP)是目前公认的求解非线性优化问题的最有效方法之一, 其收敛性好、计算效率高[11]。SQP的算法分为两部分, 主迭代和子迭代。主迭代产生一个序列, 并求其收敛。在主迭代的每步迭代上, 新产生一个二次规划子问题, 用来指定主迭代的迭代方向, 而子迭代本身的求解也是一个迭代过程。
SQP算法迭代如下:
① 给定初始点x0、收敛精度ε、令二次导数矩阵H0=E, 置k=0;
② 在点xk简化原问题为二次规划问题;
③ 求解二次规划问题, 并令Sk=S*;
④ 在方向Sk上对原问题目标函数进行一维搜索, 得点xk+1;
⑤ 判断xk+1是否满足给定精度, 符合则令k=k+1、f*=f(xk+1), 输出最优解, 终止计算, 否则转⑥;
⑥ 修正二次导数矩阵Hk+1, 令k=k+1, 转到步骤②继续迭代。
2 外形优化仿真实例 2.1 前处理利用外形优化设计框架对翼身融合水下滑翔机开展外形优化设计研究。初始水下滑翔机由5个NACA翼型放样得到, 如图 2所示。考虑到数值模拟的计算量及直接优化整机可能造成左右不对称问题, 因此选用一半的滑翔机外形进行优化。
![]() |
图 2 翼身融合式水下滑翔机 |
使用ICEM软件绘制结构化网格, 计算域为长方体域, 大小为20 m×20 m×8 m, 水下滑翔机位于计算域中心。设置入口边界条件为速度入口, 速度大小为滑翔机航速0.514 m/s, 出口压力设置为0。边界层采用O-Block来模拟, 取y+值为1, 第一层网格为0.000 2 m, 最终网格数共计932 432。数值模拟采用RANS方程作为控制方程, k-ω湍流模型作为封闭方程, 计算稳态条件下的滑翔机流体动力参数。数值离散采用二阶迎风格式, 设置收敛条件为迭代2 000步或残差收敛小于1×10-6。此外, 求解器耦合有Adjoint方法[12], 可以快速高效地求得目标函数对设计变量的梯度。
绘制其外形FFD控制体如图 3所示。控制体共计120个控制点, 内部嵌有滑翔机表面网格, 利用FFD参数化实现建立控制点与几何表面网格节点的映射关系。
![]() |
图 3 FFD控制体 |
在开展外形优化设计工作之前, 需要对网格无关性进行验证。对4种不同大小的网格进行数值仿真, 网格量分别为385 235(y+=1), 485 938(y+=1), 932 432(y+=1), 205 058 9(y+=0.5)。分别计算在0, 1°, 2°3°, 4°, 5°攻角下的升力因数和阻力因数, 计算结果如图 4所示。结果表明, 随着网格量的增加, 仿真结果越趋于一致。考虑到数值模拟和优化过程的计算量, 最终选用932 432的网格。
![]() |
图 4 网格无关性验证 |
浮力驱动式的水下滑翔机通过重浮心调节, 实现水下的周期性、锯齿形运动。考虑水下滑翔机的周期性上浮与下沉运动, 以3°攻角下沉和-3°攻角上浮, 0°侧滑角, 1节流速, 环境温度20℃, 海水密度998.2 kg/m3, 一个标准大气压作为工作条件。外形优化设计应满足:
·一个运动周期内实现总的阻力因数最小。具体为:分别计算在3°和-3°攻角下的阻力系数, 取其均值为总的阻力因数, 即令Cd, total=(Cd下沉+Cd上浮)/2达到最小;
·体积应约束在合理范围内, 确保内部空间足够安装设备及电路。具体为优化后外形与初始外形体积比约束在0.99~1之间, 即0.99≤V/Vbase≤1。
按照以上要求提出外形优化问题:
目标函数为最小化一个周期运动的阻力因数。
设计变量为FFD控制点的位移, 控制点可以沿x, y, z 3个方向运动(命名为x, y, z)。
约束条件为体积不等式约束。
完整的外形优化设计问题如表 1所示。
优化迭代收敛曲线如图 5所示, 共迭代43次, 相比于初始外形, 一个运动周期的总阻力系数降低了16.45%, 同时, 最优外形与初始外形的体积比为0.9971, 满足0.99~1的体积比约束条件。
![]() |
图 5 迭代过程 |
对比优化前后的FFD控制体外形如图 6所示, 其中红色方点为初始控制点位置, 绿色圆点为优化后控制点位置。在优化过程中, 控制点产生位移进行移动, 由于控制体和几何外形之间存在的映射关系, 初始外形形状产生相应的变形。以滑翔机机梢、机翼与机体连接处、纵向对称面处的3个纵向截面和1个横向对称处截面为例, 对比优化前后的外形变形如图 7所示, 绿色圆点代表优化后截面形状, 红色方点代表优化前的截面形状。优化后滑翔机后缘机翼与机体连接处附近弦长减少; 机翼处的厚度减少, 后掠角略微增大; 其余位置也进行了调整。
![]() |
图 6 优化前后FFD控制体 |
对比优化前后的相关流体参数如表 2所示。
参数 | 优化前 | 优化后 | 变化率/% |
下沉阻力因数 | 0.024 87 | 0.020 91 | 15.95 |
上浮阻力因数 | 0.024 87 | 0.020 65 | 17.00 |
总阻力因数 | 0.024 87 | 0.020 78 | 16.45 |
滑翔机下沉和上浮的阻力因数均得到降低, 一个周期的总阻力因数大大降低。分析阻力因数降低的原因, 水下滑翔机的阻力主要为摩擦阻力和压差阻力, 其中摩擦阻力主要与表面积大小有关, 压差阻力与外形表面压力分布有关。图 7d)明显可以看出优化后外形的表面积减少。对比优化前后在3°攻角下的压力云图, 如图 8所示。可以看出在保持原有流线平行分布的同时, 头部低压区减小, 压力分布得到调整。
![]() |
图 7 优化前后截面形状对比 |
![]() |
图 8 压力云图 |
本文对一种翼身融合水下滑翔机开展了外形优化设计研究, 提出了一套外形优化框架, 包括FFD参数化方法用于参数化, CFD求解器用于数值模拟, 网格变形算法用于计算网格更新, SQP算法用于优化迭代。外形优化设计以最小化一个周期运动的平均阻力系数为目标函数, 在保证体积约束的情况下, 翼身融合水下滑翔机一个周期的平均阻力系数降低了16.45%, 同时在下沉与上浮阶段的阻力系数均得到了降低。
[1] | DHANAK M R, XIROS N I. Springer Handbook of Ocean Engineering[M]. Cham, Switzerland: Springer, 2016: 301-322. |
[2] | CARMICHAEL R L, ERICKSON L L. Pan Air-a Higher Order Panel Method for Predicting Subsonic or Supersonic Linear Potential Flows about Arbitrary Configurations[C]//AIAA 14th Fluid & Plasma Dynamics Conference, 1981 |
[3] |
王丹.飞行器气动外形优化设计方法研究与应用[D].西安: 西北工业大学, 2015 WANG Dan. Approach and Application Research on Aerodynamic Shape Optimization Design[D]. Xi'an: Northwestern Polytechnical University, 2015(in Chinese) |
[4] |
孙春亚.翼身融合水下滑翔机外形设计与运动分析[D].西安: 西北工业大学, 2017 SUN Chunya. Shape Design and Motion Analysis of Blended-Wing-Body Underwater Glider[D]. Xi'an: Northwestern Polytechnical University, 2017(in Chinese) |
[5] | CHENGSHAN L, PENG W, HUACHAO D, et al. A Simplified Shape Optimization Strategy for Blended-Wing-Body Underwater Gliders[J]. Structural and Multidisciplinary Optimization, 2018, 58: 2189-2202. DOI:10.1007/s00158-018-2005-4 |
[6] |
陈颂, 白俊强, 史亚云, 等. 民用客机翼/机身/平尾构型气动外形优化设计方法研究[J]. 航空学报, 2015, 36(10): 3195-3207.
CHEN Song, BAI Junqiang, SHI Yayun, et al. Aerodynamic Optimization Design of Civil Jet Wing-Body-Tail Configuration[J]. Acta Aeronauticaet Astronautica Sinica, 2015, 36(10): 3195-3207. (in Chinese) |
[7] |
汪怡平, 王涛, 黎帅. 基于自由变形技术的汽车气动减阻优化[J]. 机械工程学报, 2017, 53(9): 135-143.
WANG Yiping, WANG Tao, LI Shuai. Aerodynamic Drag Reduction of Vehicle Based on Free Form Deformation[J]. Journal of Mechanical Engineering, 2017, 53(9): 135-143. (in Chinese) |
[8] | LI J, WANG P, DONG H, et al. Shape Optimisation of Blended-Wing-Body Underwater Gliders Based on Free-Form Deformation[J]. Ships and Offshore Structures, 2020, 15(3): 227-235. DOI:10.1080/17445302.2019.1611989 |
[9] | WANG X, SONG B, PENG W, et al. Hydrofoil Optimization of Underwater Glider Using Free-Form Deformation and Surrogate-Based Optimization[J]. International Journal of Naval Architecture & Ocean Engineering, 2018, 10(6): 730-740. |
[10] |
张来平, 邓小刚, 张涵信. 动网格生成技术及非定常计算方法进展综述[J]. 力学进展, 2010, 40(4): 424-447.
ZHANG Laiping, DENG Xiaogang, ZHANG Hanxin. Reviews of Moving Grid Generation Techniques and Numerical Methods for Unsteady Flow[J]. Advances in Mechanics, 2010, 40(4): 424-447. (in Chinese) |
[11] | GILL P E, MURRAY W, SAUNDERS M A. SNOPT:an SQP Algorithm for Large-Scale Constrained Optimization[J]. Society for Industrial and Applied Mathematics, 2005, 47(1): 99-131. |
[12] | WU X, WANG P, LI J, et al. Adjoint-Based Optimization for Blended-Wing-Body Underwater Gliders' Shape Design[C]//Oceans-MTS/IEEE Kobe Techno-Oceans Conference, 2018 |