超高速水下航行体利用超空泡技术具有的减阻效应,实现了超高速航行。在超高速航行体的航行过程中,加速段是空泡形成阶段[1],该阶段空泡由局部空泡发展为超空泡,雷体的流体动力特性发生很大改变[2],航速在0.08 s内从30 kn迅速增长到200 kn左右。
航速跨度大时,常规舵机舵面因处于全沾湿状态,舵机伺服系统无法到达操舵要求,而双自由度舵机伺服系统,不仅可以控制舵面偏转改变舵角,还能通过减小舵面的沾湿面积改变舵效,满足超高速航行体的控制要求。
在航速迅速变化过程中,舵板会承受水流产生的瞬态冲击载荷,由于舵面直接与双自由度舵机传动机构相连,所以校核齿轮箱与传动机构的抗冲击能力是必要的。
1 旋转机构旋转机构是控制舵面偏转改变舵角的执行机构,设计中采用了九级齿轮减速机构,减速比约为803。机构输出轴的轴端通过过渡块与舵板固定连接,连接方式均为螺栓连接,装配图如图1所示。
![]() |
图 1 舵板装配图 Fig. 1 Rudder plate assembly drawing |
超高速航行体在加速阶段,舵板上的瞬态冲击载荷会对旋转结构产生影响。根据经验,对机构最后一级传动与齿轮箱影响最大,通过三维建模软件UG导出实体模型,并导入到Ansys workbench软件中进行仿真分析[3],验证结构在冲击下的瞬态动力学响应。
1.1.1 模型前处理瞬态动力学分析是在随时间变化载荷下分析结构变化并确定动态响应的分析方式。将模型导入软件中,采用45钢材料[4],泊松比
![]() |
图 2 旋转机构载荷与约束加载模型 Fig. 2 Load and constrained loading model of rotating mechanism |
在瞬态动力学仿真计算后,查看结构总变形、等效应力以及等效位移云图,如图3所示。
![]() |
图 3 旋转机构强度校核结果 Fig. 3 Rotating mechanism strength check results |
根据分析结果可知,输出轴总的最大变形为0.0030081 mm,最大等效应力为160.14 MPa,都发生在舵板连接的螺栓孔中,最大等效应力低于材料的屈服极限355 MPa,最小安全系数为1.5611,说明输出轴韧性满足设计要求,可以承受瞬态冲击载荷。
1.2 输出轴模态分析旋转机构输出轴为齿轮轴,为了防止齿轮啮合的激振频率接近轴的固有频率而发生共振[7],需要对输出轴进行模态分析。经分析,得到轴的各阶振型如图4所示,固有频率如表1所示。
![]() |
图 4 前7阶模态振型图 Fig. 4 Figure of the first seven modes |
![]() |
表 1 旋转机构输出轴模态频率与振型 Tab.1 Mode frequency and shape of output shaft of rotating mechanism |
齿轮啮合的激振频率[8]为:
${f} = Z \times {n}/60\text{。}$ | (1) |
式中:Z为齿轮轴上的齿数,值为31,n为轴的转速,频率约为0.58 Hz,激振频率与轴的固有频率差距较大,不会发生共振。
2 展开机构展开机构通过2个短轴带动旋转机构和舵板一起做展开运动,改变舵板的沾湿面积,采用8级齿轮传动,传动比约为235.2。在受到冲击载荷时,整个旋转机构对展开机构最后一级啮合齿轮产生力矩,验证该机构最后一级齿轮在冲击载荷下的响应。
2.1 前处理将模型导入workbench的Transient Structural模块中,同样采用45钢材料,自动生成网格,得到31094个节点,17204个单元。对两齿轮添加Body-Ground的旋转副,在主动轮的旋转副上添加joint load,类型为rotational velocity,速度值为0,被动轮上添加joint load,类型为moment,值为50 N·m,方向与齿轮转向相反,两齿轮齿面上添加摩擦接触,添加结果如图5所示。
![]() |
图 5 展开机构输出级齿轮载荷与约束加载模型 Fig. 5 Load and constrained loading model of output gear of expansion mechanism |
仿真时间设为0.08 s,结构总变形、等效应力以及安全系数如图6所示。
![]() |
图 6 展开机构输出级齿轮校核结果 Fig. 6 Check results of output gear of expansion mechanism |
由仿真结果可知,最大变形为0.053187 mm,最大等效应力为534.25 MPa,出现在两齿面啮合处,高于材料的屈服极限355 MPa,从安全系数中可以看到最小安全系数为0.4679,说明使用该材料时,在载荷冲击工况下展开机构输出级啮合齿轮会发生破坏变形,需要选用屈服极限更高的材料或对其进行热处理提高齿轮韧性。
3 旋转机构齿轮箱校核旋转机构齿轮箱分为上盖板、箱体以及下盖板3部分组成,通过螺栓连接。其中,输出轴上的轴承安装在上盖板的孔中,在舵面受到载荷冲击时,减速箱体同样会受到冲击,对箱体的校核同样重要。
3.1 前处理将箱体实体文件导入workbench的Transient Structural模块中,材料采用Q235结构钢,泊松比
![]() |
图 7 旋转机构减速箱载荷与约束加载模型 Fig. 7 Load and constrained loading model of rotating gear reducer |
仿真时间设为0.08 s,结构总变形和等效应力如图8所示。
![]() |
图 8 旋转机构减速箱校核结果 Fig. 8 Check results of rotating gear reducer |
根据分析结果可知,减速箱总的最大变形为0.01244 mm,发生在与上盖板的轴承孔中,最大等效应力为155.56 MPa,小于Q235的屈服极限235 MPa,说明箱体结构设计与材料选用可以承受该工况下的载荷冲击。
4 结 语本文针对超高速水下航行体在加速过程中,由于速度跨度大导致双自由度舵机承受较大瞬态冲击载荷的工况进行了仿真校核,通过建立旋转机构输出轴、减速箱体和展开机构最后一级齿轮副的多刚体动力学仿真模型,进行有限元分析计算,校核各结构抵抗冲击的韧性,并对旋转机构输出轴进行模态分析,仿真与校核结果为机构优化设计提供了依据。
[1] |
王改娣. 超空泡鱼雷技术特点分析[J]. 鱼雷技术, 2007(5): 1-4. WANG Gai-di. Analysis of technica lfeatures of supercavitating torpedoes[J]. Torpedo Technology, 2007(5): 1-4. |
[2] |
陈超倩, 曹伟, 王聪, 等. 超空泡航行体加速段控制设计[J]. 哈尔滨工业大学学报, 2016, 48(8): 147-153. CHEN Chaoqian, CAO Wei, WANG Cong, et al. Acceleration stage control design for supercavitating vehicle[J]. Journal of Harbin Institute of technology, 2016, 48(8): 147-153. |
[3] |
刘涛. 基于Ansys workbench的汽车高速轴动态分析[J]. 科技风, 2019(7): 137. |
[4] |
舒彪, 喻道远, 王灯, 等. 采用UG、HyperMesh和Ansys的齿轮轴模态分析[J]. 现代制造工程, 2012(2): 71-73+121. DOI:10.3969/j.issn.1671-3133.2012.02.017 |
[5] |
向禹豪. 基于Ansys齿轮轴的模态分析[J]. 内燃机与配件, 2019(11): 37-39. DOI:10.3969/j.issn.1674-957X.2019.11.016 |
[6] |
张进军. 有限元分析及Ansys workbench工程应用[M]. 西安: 西北工业大学出版社, 2018.
|
[7] |
彭飞, 黄勇, 李准. 基于Ansys Workbench振动筛激振器轴强度有限元分析[J]. 煤矿机械, 2020, 41(3): 173-176. PENG Fei, HUANG Yong, LI Zhun. Finite Element Analysis of Strength of Vibrating Screen Vibrator Shaft Based on ANSYS Workbench[J]. Coal Mine Machinery, 2020, 41(3): 173-176. |
[8] |
史金奇, 张秦赓, 王立文, 等. 鱼雷电动伺服舵机动力学仿真[J]. 鱼雷技术, 2016, 24(6): 463-468. SHI Jin-qi, ZHANG Qin-geng, WANG Li-wen. Dynamic Simulation of Torpedo′s Electric Servo Actuator[J]. Torpedo Technology, 2016, 24(6): 463-468. |
[9] |
岗海明. 高速齿轮箱体的有限元分析及结构优化设计[D]. 北京: 机械科学研究总院, 2016.
|