1. 引言
水面与水下航行器的减阻一直以来都是相关领域研究的焦点与难点,柔性壁减阻技术由于不需要对航行器结构进行更改,能适用于大部分壳体形状等优点,逐渐成为研究热点之一。
柔性壁减阻最早来源于1960年Kramer [1] 的研究。他模拟海豚表皮,提出了使用柔性结构可延迟层流向湍流的转捩。Riley等人 [2] 在1988年提出了一种均匀分布的弹性支撑蒙皮结构,其由柔性蒙皮和均匀布置的离散弹性支撑组成。他们通过实验测量得出,所设计的结构可以有效抑制流场中扰动。
由于湍流流动的复杂性,目前对于湍流中的减阻研究大多是以数值计算、实验分析的方式进行的。1993年Lee等人 [3] 通过观察粘弹性柔性平板上湍流边界层,发现柔性壁和边界层在一定雷诺数范围内(1700~3000)相互作用,使湍流边界层对数区发生了明显的上移现象。1997年,Choi等人 [4] 再一次通过实验证实了减阻效果与速度型对数区变化的相关性。2002年,赵汉中 [5] 通过数值分析,在理论上证实了上述湍流边界层速度型变化和减阻效果的关联具有普适性。2010年,陈冰鹤 [6] 使用优化后的Falkner-Skan对湍流求解,其结果与试验数据相匹配。2012年,罗敏和胡建成 [7] 进一步验证了Falkner-Skan方程求解的准确性。2018年,郝路亭 [8] 对蒙皮结构进行仿真和试验验证,证明了蒙皮结构具有减阻能力。
![](//html.hanspub.org/file/1-2650260x7_hanspub.png?20220823082305792)
Figure 1. Physics model of non-uniform spring array supported skin
图1. 非均匀弹簧支撑蒙皮结构物理模型
本文首先设计了一种非均匀弹簧阵列支撑蒙皮结构,并根据湍流边界层速度型对数区变化与减阻效果存在的相关性,验证该结构具有的减阻特性,以及各结构参数变化对减阻效果的影响。
2. 非均匀弹簧阵列支撑蒙皮结构物理模型
非均匀弹簧阵列支撑蒙皮结构主要由下层的非均匀布置的弹簧整列和上层的蒙皮材料两部分组成,如图1所示。弹簧阵列固定在刚性表面上,沿流体流动方向按照1:2的距离间隔布置,而垂直于流体流动方向均匀布置;蒙皮固定在弹簧阵列上。
3. 选取湍流边界层速度型
根据Lee、赵汉中等人的研究结果 [3] [5],蒙皮结构减阻效果与其表面湍流边界层速度型存在关系。所以,可以通过分析湍流边界层速度型来研究非均匀弹簧阵列支撑蒙皮结构具备的减阻效果。
在分析边界层速度型之前,需要首先得到边界层速度型的控制方程。设边界层无量纲速度:
经过Cebeci和Smith坐标变换的湍流Falknar-Skan方程为:
(1)
其中
为关于
的函数,
为涡粘系数的无量纲形式。
根据混合长度理论,
表达式为:
(2)
其中
为排挤厚度,
为内层湍流与外层湍流的分界点,U为边界层外流速,
为动力粘度。
为了方便计算,减少变量数量,设
对式(1)做进一步变换如下:
(3)
变换完毕的式子,即式(3),为非均匀弹簧阵列支撑蒙皮结构在湍流中的表面边界层速度型的控制方程。
得到边界层速度型的控制方程后,需要将式(3)导入Matlab编程,使用有限差分方法对式子进行变换,并选用流速
、密度为
、运动黏度为
、阻尼长度
的流体流过结构,迭代求解 [9] 并生成结构在流向不同位置处的边界层速度型如图2所示。由图可知,随着位置不断向后,也就x值不断增大,边界层速度型形状开始逐渐发生变化,意味着流体随着x值增大而发生转捩,由层流转变为层流湍流混合,最终完全转变为湍流状态。观察可得,在x = 2.4处,速度型已经明显体现为湍流状态下的速度型形状,可以认为此时流体已完全转捩;因此,选取x = 2.4处的曲线作为结构在湍流下的减阻特性的研究对象。
![](//html.hanspub.org/file/1-2650260x24_hanspub.png?20220823082305792)
Figure 2. Velocity characteristics at different location
图2. 不同位置下的速度型
蒙皮结构减阻效果与其表面湍流边界层速度型存在关系,速度型对数区发上移越多,结构的减阻效果越好。因此,下文将通过观察速度型对数区上移的程度对非均匀弹簧阵列支撑蒙皮结构的减阻效果进行分析。
4. 非均匀弹簧阵列支撑蒙皮结构与其它结构减阻效果对比
目前广泛研究的蒙皮减阻结构主要有两种,一种为单层蒙皮结构,另一种为均匀弹簧阵列支撑蒙皮结构。使用Matlab求得上述两种结构与非均匀弹簧阵列支撑蒙皮结构的湍流边界层速度型如图3所示;为了控制变量,保证对比结果有效,三种结构使用相同的蒙皮参数,非均匀与均匀弹簧阵列使用相同的弹簧参数。求解。观察图3可知,非均匀弹簧阵列支撑蒙皮结构的速度型对数区上移的程度高于单层蒙皮结构与均匀弹簧阵列支撑蒙皮结构,说明其相较于后两者具有更好的减阻效果。
![](//html.hanspub.org/file/1-2650260x25_hanspub.png?20220823082305792)
Figure 3. Velocity characteristics comparison of three different structures
图3. 三种结构的速度型对比
5. 结构参数对速度型影响
蒙皮结构中各参数的变化将对整体结构的减阻特性产生影响,这些参数包括弹簧间距、弹簧刚度、蒙皮材料阻尼、蒙皮材料刚度。
1) 弹簧间距对速度型影响
研究非均匀弹簧阵列支撑蒙皮结构的弹簧间距对速度型的影响,可以通过研究无量纲单元长度L来实现。各长度L下的速度型如图4所示。观察图4可得,单元长度L增大,各L参数下得到的速度型的对数区出现上移趋势,说明增大单元长度L时,非均匀弹簧阵列支撑蒙皮结构的减阻效果提高。同时注意到,速度型在层流次层出现下移趋势,这是由于粘性力增大而使得层流次层的区域得到扩大。
![](//html.hanspub.org/file/1-2650260x26_hanspub.png?20220823082305792)
Figure 4. Velocity characteristics under different spring spacing
图4. 不同弹簧间距下的速度型
2) 弹簧刚度对速度型影响
通过改变弹簧刚度k来研究其对非均匀弹簧阵列支撑蒙皮结构减阻效果的影响。不同弹簧刚度k下的速度型如图5所示。由图5可知,当弹簧刚度k增大时,各k参数下得到的速度型的对数区出现下移,说明增大弹簧刚度k,非均匀弹簧阵列支撑蒙皮结构的减阻效果降低。
![](//html.hanspub.org/file/1-2650260x27_hanspub.png?20220823082305792)
Figure 5. Velocity characteristics under different spring stiffness
图5. 不同弹簧刚度下的速度型
3) 蒙皮材料阻尼对速度型影响
通过改变蒙皮材料阻尼系数d来研究材料阻尼对非均匀弹簧阵列支撑蒙皮结构的减阻效果的影响。不同蒙皮材料阻尼系数d的速度型如图6所示。由图6可以看出,当蒙皮材料阻尼系数d增大,各阻尼系数d下得到的速度型的对数区逐渐下移,说明增大蒙皮材料阻尼系数d时,非均匀弹簧阵列支撑蒙皮结构减阻效果降低。
![](//html.hanspub.org/file/1-2650260x28_hanspub.png?20220823082305792)
Figure 6. Velocity characteristics under different damping coefficient
图6. 不同阻尼系数下的速度型
4) 蒙皮材料刚度对速度型影响
通过改变蒙皮材料剪切波速
来研究蒙皮刚度对非均匀弹簧阵列支撑蒙皮结构减阻效果的影响。不同蒙皮材料剪切波速
下的速度型如图7所示。由图7可知,当蒙皮材料剪切波速
增大,各蒙皮材料剪切波速
下得到的速度型的对数区发生下移,说明增大蒙皮材料剪切波速
时,非均匀弹簧阵列支撑蒙皮结构的减阻效果降低。
![](//html.hanspub.org/file/1-2650260x34_hanspub.png?20220823082305792)
Figure 7. Velocity characteristics under different shear wave velocity
图7. 不同蒙皮材料剪切波速下的速度型
6. 结论
本文提出了一种非均匀弹簧阵列支撑蒙皮结构,并建立了结构的物理模型,根据Falknar-Skan方程和混合长度理论,建立了结构的速度型控制方程;使用Matlab,对速度型控制方程进行迭代求解,得到结构在湍流边界层下的速度型;通过观察速度型对数区上移的情况对比分析了非均匀弹簧阵列支撑蒙皮结构相对其他结构的减阻效果,并分析了各结构参数对减阻效果的影响。
结果表明,本文提出的非均匀弹簧阵列支撑蒙皮结构的湍流边界层速度型的对数区上移的程度高于同等条件下的单层蒙皮结构与均匀弹簧阵列支撑蒙皮结构,说明此结构在减阻性能上相较于后两者具有优势。同时,根据分析得到,可通过增大弹簧间距,减小弹簧刚度、减小蒙皮材料的阻尼和刚度等方式提高结构得到减阻效果。
参考文献