CN113806987A - 一种推进系统多目标优化设计方法及电子设备 - Google Patents
一种推进系统多目标优化设计方法及电子设备 Download PDFInfo
- Publication number
- CN113806987A CN113806987A CN202111130307.XA CN202111130307A CN113806987A CN 113806987 A CN113806987 A CN 113806987A CN 202111130307 A CN202111130307 A CN 202111130307A CN 113806987 A CN113806987 A CN 113806987A
- Authority
- CN
- China
- Prior art keywords
- power consumption
- shafting
- bearing
- friction power
- calculation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Sliding-Contact Bearings (AREA)
Abstract
本发明涉及船舶动力系统设计技术领域，尤其涉及一种推进系统多目标优化设计方法及电子设备。通过以轴系传动效率为重点，兼顾提高轴系校中质量、改善振动特性为目标开展多目标综合优化分析，能够降低轴系功耗，同时兼顾提高轴系校中质量、改善振动特性；通过建立推力轴承摩擦功耗、径向轴承摩擦功耗、密封摩擦功耗以及转轴摩擦功耗的计算模型，能够准确高效地计算推进系统的传动效率；通过采用有限差法和数学归纳，使得摩擦功耗的计算结果更加准确，且能有效降低摩擦功耗的计算时间。
Description
技术领域
本发明涉及船舶动力系统设计技术领域，尤其涉及一种推进系统多目标优化设计方法及电子设备。
背景技术
船舶在航行过程中，主机输出的功率通过减速齿轮箱经长、短轴系传递给螺旋桨，螺旋桨产生的推力又通过轴系上的推力轴承传递给船体，实现推进船舶的使命。在船舶主机齿轮箱输出到螺旋桨之间的功率传递过程中，推力轴承、中间轴承、艉轴承等接触应力较大的部位，以及在水中运转的部分轴系，会消耗掉主机的一部分输出功率，使螺旋桨接收到的功率降低。
针对大型船舶的长轴系，由于轴系长度较长，且功率和扭矩较高，轴径也较大，在轴系的功率传递过程中经过主推力轴承摩擦损失、中间轴承摩擦损失、密封装置摩擦损失、水润滑轴承摩擦损失、水中轴段水阻力损失等，对轴系传动效率有较大的影响。
目前，常规的推进系统设计方法，仅对轴系的校中和振动进行分析，优化支撑负荷和振动特性，但无法兼顾轴系的传动效率和轴系功耗，因此行业内急缺一种能够同时兼顾传动效率、校中质量和振动情况的推进系统多目标优化设计方法。
发明内容
本发明要解决的技术问题是：提供一种推进系统多目标优化设计方法及电子设备，能够提高推进系统传动效率，降低轴系功耗，同时兼顾提高轴系校中质量、改善振动特性。
为解决上述技术问题，本发明采用的技术方案是：
一种推进系统多目标优化设计方法，采用多目标综合分析策略，在轴系校中分析的基础上进行摩擦功耗分析，优化轴系的传动效率，同时结合轴系振动分析，优化轴系的振动特性，所述多目标综合分析策略具体包括：
A、结合静态校中结果进行轴承的润滑特性计算，获得轴承的摩擦功耗，同时进行动态校中，获得最终的摩擦功耗以及轴承动刚度和阻尼系数；
B、以动态校中获得的动刚度、阻尼系数为输入，获得各支撑处的扰动信息，进行考虑振动的轴承摩擦功耗计算，将振动分析与传动效率计算进行耦合；
C、构建轴系综合分析模型，包括轴系校中模块、轴系功耗模块、轴系振动模块，通过所述摩擦功耗、动刚度、阻尼系数、扰动信息将上述模块相关联，以轴承的径向变位及轴系转速作为优化变量，以轴系校中的质量作为约束条件，以轴系传动效率、支撑负荷、扰度曲线、临界转速、振动响应作为优化目标，进行推进系统多目标优化设计。
进一步的，所述一种推进系统多目标优化设计方法，具体流程如下：
S1、建立输入参数库，包括轴系几何、工况参数；
S2、轴系校中模块中的静态校中模块计算得出各轴承负荷，并传递给轴系功耗模块中的轴承模块；
S3、所述轴承模块进行轴承摩擦功耗、动刚度、阻尼系数的计算，并将动刚度、阻尼系数传递至轴系校中模块中的动态校中模块，进行动态校中，获取最终轴承摩擦功耗、动刚度、阻尼系数、轴段扰度；
S4、根据最终的轴承动刚度和阻尼系数进行轴系回旋振动计算，得到轴颈响应值，传递给所述轴承模块，计算考虑振动的轴承摩擦功耗；
S5、所述轴系动态校中模块计算求得的轴段扰度，传递给轴系功耗模块中的转轴模块，计算求得转轴摩擦功耗；
S6、所述轴系功耗模块在完成轴承摩擦功耗、密封摩擦功耗、转轴摩擦功耗以及其它功耗计算的基础上，计算求得轴系传动效率；
S7、以轴承的径向变位及轴系转速作为优化变量，以轴系校中的质量作为约束条件，将轴系传动效率、支撑负荷、扰度曲线、临界转速、振动响应作为优化目标，进行推进系统多目标优化设计。
进一步的，所述轴承摩擦功耗包括推力轴承摩擦功耗、中间轴承摩擦功耗和艉轴承摩擦功耗，所述中间轴承摩擦功耗和艉轴承摩擦功均为径向轴承摩擦功耗。
进一步的，所述推力轴承摩擦功耗的计算包括如下步骤：
S1、根据推力轴承结构、工况参数以及轴承初始支点膜厚，求解膜厚方程；
S2、根据膜厚方程计算结果，求解稳态雷诺方程；
S3、判断计算出的轴承压力是否收敛，若是，则计算液膜水平力和垂直力，若否，则返回执行步骤S2；
S4、计算液膜合力；
S5、判断轴承载荷与液膜合力是否相等，若是，则执行步骤S6，若否，则修正支点膜厚并返回执行步骤S1；
S6、求解静特性方程，并输出推力轴承摩擦功耗计算结果。
进一步的，所述径向轴承摩擦功耗计算包括如下步骤：
S1、根据径向轴承结构、工况参数以及初始偏心率和偏位角，求解膜厚方程；
S2、根据膜厚方程计算结果，求解稳态雷诺方程；
S3、判断计算出的轴承压力是否收敛，若是，则计算液膜水平力和垂直力，若否，则返回执行步骤S2；
S4、判断偏位角是否收敛，若是，则计算液膜合力，若否，则修正偏位角并返回执行步骤S1；
S5、判断轴承载荷与液膜合力是否相等，若是，则执行步骤S6，若否，则修正偏心率并返回执行步骤S1；
S6、求解静特性方程和瞬态雷诺方程；
S7、求解径向轴承动刚度和阻尼系数；
S8、输出径向轴承摩擦功耗计算结果。
进一步的，所述密封摩擦功耗计算包括如下步骤：
S1、根据密封件的结构和转轴工况参数建立密封结构模型；
S2、对密封结构模型进行网格划分；
S3、进行有限元计算；
S4、判断有限元计算结果是否收敛，若是，则计算摩擦力矩并输出计算结果，若否，则返回执行步骤S3；
S5、对多工况点的摩擦力矩计算结果进行数学归纳；
S6、判断数学归纳结果误差是否在允许范围内，若是，则输出密封摩擦功耗计算结果，若否，则返回执行步骤S5。
进一步的，所述转轴摩擦功耗计算包括如下步骤：
S1、根据敞水中轴段轴系结构和工况参数建立敞水中轴段流场模型；
S2、对敞水中轴段流场模型进行网格划分；
S3、进行Fluent流体有限元计算；
S4、判断流体有限元计算结果残差值是否收敛，若是，则计算摩擦力矩并输出计算结果，若否，则返回执行步骤S2；
S5、对多工况点的摩擦力矩计算结果进行数学归纳；
S6、判断数学归纳结果误差是否在允许范围内，若是，则输出转轴摩擦功耗计算结果，若否，则返回执行步骤S5。
进一步的，所述其他功耗包括隔音联轴器、隔舱密封装置、遥测功率仪的功率损失，通过加入与所述其他功耗相匹配的修正系数对轴系功率损失进行修正。
进一步的，所述轴系几何、工况参数具体包括轴系几何参数、船体变形、螺旋桨激励以及轴系转速。
一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于：所述处理器执行所述程序时实现如上述任一项所述的一种推进系统多目标优化设计方法。
本发明与现有技术相比具有以下主要的优点：
1、通过以轴系传动效率为重点，兼顾提高轴系校中质量、改善振动特性为目标开展多目标综合优化分析，能够降低轴系功耗，同时兼顾提高轴系校中质量、改善振动特性；
2、通过建立推力轴承摩擦功耗、径向轴承摩擦功耗、密封摩擦功耗以及转轴摩擦功耗的计算模型，能够准确高效地计算推进系统的传动效率；
3、通过采用有限差法和数学归纳，使得摩擦功耗的计算结果更加准确，且能有效降低摩擦功耗的计算时间。
附图说明
图1为本发明技术路线示意图；
图2为船舶推进系统能量耗散图；
图3为本发明推力轴承摩擦功耗计算流程图；
图4为本发明径向轴承摩擦功耗计算流程图；
图5为本发明密封摩擦功耗计算流程图；
图6为本发明转轴摩擦功耗计算流程图；
图7为本发明推进系统多目标优化设计方法架构示意图；
图8为本发明推进系统多目标优化设计方法流程示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
需要指出，根据实施的需要，可将本申请中描述的各个步骤/部件拆分为更多步骤/部件，也可将两个或多个步骤/部件或者步骤/部件的部分操作组合成新的步骤/部件，以实现本发明的目的。
本发明实施的一种推进系统多目标优化设计方法，如图1所示，在常规推进系统设计方法优化支撑负荷的基础上，基于轴系校中分析结果进行摩擦功耗分析，优化轴系的传动效率；同时进行轴系振动分析，优化轴系的振动特性。主要包括如下步骤：
A、结合静态校中结果进行各轴承的润滑特性计算，获得各轴承的摩擦功耗。同时，在静态校中的基础上，进一步考虑动态校中，获得动态校中下的动刚度、阻尼信息与对应的摩擦功耗，将动态校中分析与传动效率计算进行耦合。
B、以动态校中下的动刚度、阻尼信息为输入，获得各支撑处的扰动信息，进一步进行考虑振动的轴承摩擦功耗计算，将振动分析与传动效率计算进行耦合。
C、构建轴系综合分析模型，包括轴系校中模型、轴系功耗模型、轴系振动模型，通过步骤A、B中所述关联参数将上述模型相关联，结合智能优化算法，以轴系传动效率、校中质量和振动性能为目标参数，进行轴系设计多目标优化，提出优化方案。
基于上述的主要步骤，本发明实施的一种推进系统多目标优化设计方法具体包括：
一、推进系统功耗计算
船舶推进系统的任务是连接主机与螺旋桨，将主机发出的功率传递给螺旋桨，同时将螺旋桨所产生的推力通过推力轴承传递给船体，推进船舶。
如图2所示，在船舶推进系统中，主要的能耗因素包括轴承摩擦功耗、密封摩擦功耗、转轴摩擦功耗以及其它功耗。
其中，所述轴承摩擦功耗包括推力轴承摩擦功耗、中间轴承摩擦功耗和艉轴承摩擦功耗，所述中间轴承和艉轴承同为径向轴承，它们的润滑模型除了黏度方程不同，其它方程基本相同。
推力轴承与径向轴承的功耗计算模型方程类型相同，主要包括雷诺方程、膜厚方程和润滑性能方程等，只是因为结构和坐标系不同，方程形式有所区别。由于推力轴承与径向轴承的润滑模型的方程相似，它们的模型求解算法基本相同，首先对模型进行离散化，然后采用有限差分法进行数值计算。
如图3所示，推力轴承摩擦功耗计算包括如下步骤：
S1、根据推力轴承结构、工况参数以及轴承初始支点膜厚，求解膜厚方程；
S2、根据膜厚方程计算结果，求解稳态雷诺方程；
S3、判断计算出的轴承压力是否收敛，若是，则计算液膜水平力和垂直力，若否，则返回执行步骤S2；
S4、计算液膜合力；
S5、判断轴承载荷与液膜合力是否相等，若是，则执行步骤S6，若否，则修正支点膜厚并返回执行步骤S1；
S6、求解静特性方程，并输出推力轴承摩擦功耗计算结果。
如图4所示，径向轴承摩擦功耗计算包括如下步骤：
S1、根据径向轴承结构、工况参数以及初始偏心率和偏位角，求解膜厚方程；
S2、根据膜厚方程计算结果，求解稳态雷诺方程；
S3、判断计算出的轴承压力是否收敛，若是，则计算液膜水平力和垂直力，若否，则返回执行步骤S2；
S4、判断偏位角是否收敛，若是，则计算液膜合力，若否，则修正偏位角并返回执行步骤S1；
S5、判断轴承载荷与液膜合力是否相等，若是，则执行步骤S6，若否，则修正偏心率并返回执行步骤S1；
S6、求解静特性方程和瞬态雷诺方程；
S7、求解径向轴承动刚度和阻尼系数；
S8、输出径向轴承摩擦功耗计算结果。
进一步的，密封摩擦功耗主要包括各轴承密封结构的摩擦功耗，根据密封件的结构和转轴工况参数通过有限元方法完成密封结构的功耗计算，为提高轴系功耗计算效率降低计算时间，通过多工况点的仿真计算结果进行数学归纳，形成密封摩擦功耗计算的拟合公式。
如图5所示，密封摩擦功耗计算包括如下步骤：
S1、根据密封件的结构和转轴工况参数建立密封结构模型；
S2、对密封结构模型进行网格划分；
S3、进行有限元计算；
S4、判断有限元计算结果是否收敛，若是，则计算摩擦力矩并输出计算结果，若否，则返回执行步骤S3；
S5、对多工况点的摩擦力矩计算结果进行数学归纳；
S6、判断数学归纳结果误差是否在允许范围内，若是，则输出密封摩擦功耗计算结果，若否，则返回执行步骤S5。
更进一步的，转轴摩擦功耗主要针对敞水中轴段轴系摩擦功耗，所述敞水中轴段轴系摩擦功耗定义为轴段或轴承、轴颈的四周处于全液体状况下的摩擦功耗。根据结构和工况参数构建水中轴段流场模型，通过有限元方法完成摩擦功耗计算，为提高轴系功耗计算效率降低计算时间，通过多工况点的仿真计算结果进行数学归纳，形成转轴摩擦功耗计算的拟合公式。
如图6所示，转轴摩擦功耗计算包括如下步骤：
S1、根据敞水中轴段轴系结构和工况参数建立敞水中轴段流场模型；
S2、对敞水中轴段流场模型进行网格划分；
S3、进行Fluent流体有限元计算；
S4、判断流体有限元计算结果残差值是否收敛，若是，则计算摩擦力矩并输出计算结果，若否，则返回执行步骤S2；
S5、对多工况点的摩擦力矩计算结果进行数学归纳；
S6、判断数学归纳结果误差是否在允许范围内，若是，则输出转轴摩擦功耗计算结果，若否，则返回执行步骤S5。
更进一步的，其他功耗包括隔音联轴器、隔舱密封装置、遥测功率仪等轴系附属设备功率损失，可通过修正系数对轴系功率损失进行修正。
二、推进系统多目标优化设计方法架构
如图7所示，构建轴系综合分析模型，包括轴系校中模型、轴系功耗模型、轴系振动模型，三个模型构成循环链路架构，由基于响应面分析的优化内核进行驱动。
分析系统首先由轴系校中模型发起，将轴承负荷、扰度等关联参数传递给轴系功耗模型，轴系功耗模型在计算得到传动效率的同时，也将动刚度、阻尼等信息传递回轴系校中模块进行动态校中，也传递给轴系振动模型中的回旋振动计算得到临界转速、振动响应等关联参数，这些关联参数传回轴系功耗模型进行计入振动的轴承摩擦功耗计算。
在整个链路中，分析系统同构关联参数(例如载荷、位移、刚度等轴系静、动态数据)将模型进行关联。此外，轴系校中模型是上述子模型的数据来源，为轴系功耗模型和轴系振动模型计算提供轴系运行状态数据。
三、推进系统多目标优化设计方法流程
如图8所示，推进系统多目标优化设计方法具体流程如下：
S1、建立输入参数库，包括轴系几何、工况参数，如轴径、坐标、弹性模量、刚度、阻尼信息、载荷等；
S2、轴系静态校中模块(也称直线校中模块)计算得出各轴承负荷，并传递给轴系功耗模块中的轴承模块；
S3、轴系功耗模块中的轴承模块进行轴承摩擦功耗、动刚度、阻尼系数的计算，并将动刚度、阻尼系数传递至轴系动态校中模块，进行动态校中，获取最终的动特性(动刚度、阻尼系数)、轴承摩擦功耗、轴段的扰度；
S4、根据最终的轴系动特性进行轴系回旋振动计算，得到轴颈响应值，传递给轴承模块，计算考虑振动的轴承摩擦功耗；
S5、轴系动态校中模块计算求得的轴段扰度，传递给轴系功耗模块中的转轴模块，计算求得转轴摩擦功耗；
S6、轴系功耗模块在完成轴承摩擦功耗、密封摩擦功耗、转轴摩擦功耗以及其它功耗计算的基础上，计算求得轴系传动效率；
S7、以轴承的径向变位及轴系转速作为优化变量，以轴系校中的质量作为约束条件，将轴系传动效率、支撑负荷、扰度曲线、临界转速、振动响应作为优化目标，进行推进系统多目标优化设计。
通过以某轴系试验平台为对象对本实例进行了仿真验证，以推进系统中两个轴承的径向变位及转速作为优化变量，以轴系功耗最小、后轴架轴承的支反力最小、同类型支承或相邻支承间负荷的偏差最小、各支撑位置振动响应幅值最小等作为优化目标，进行本实例所述的推进系统多目标优化设计获得优化结果，结果显示优化目标均得到了不同程度的优化。
采用上述的一种推进系统多目标优化设计方法，通过以轴系传动效率为重点，兼顾提高轴系校中质量、改善振动特性为目标开展多目标综合优化分析，能够降低轴系功耗，同时兼顾提高轴系校中质量、改善振动特性；通过建立推力轴承摩擦功耗、径向轴承摩擦功耗、密封摩擦功耗以及转轴摩擦功耗的计算模型，能够准确高效地计算推进系统的传动效率；通过采用有限差法和数学归纳，使得摩擦功耗的计算结果更加准确，且能有效降低摩擦功耗的计算时间。
基于同一发明构思，本申请实施例还提供了一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现上述方法的所有方法步骤或部分方法步骤。
本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种推进系统多目标优化设计方法，其特征在于，采用多目标综合分析策略，在轴系校中分析的基础上进行摩擦功耗分析，优化轴系的传动效率，同时结合轴系振动分析，优化轴系的振动特性，所述多目标综合分析策略具体包括：
A、结合静态校中结果进行轴承的润滑特性计算，获得轴承的摩擦功耗，同时进行动态校中，获得最终的摩擦功耗以及轴承动刚度和阻尼系数；
B、以动态校中获得的动刚度、阻尼系数为输入，获得各支撑处的扰动信息，进行考虑振动的轴承摩擦功耗计算，将振动分析与传动效率计算进行耦合；
C、构建轴系综合分析模型，包括轴系校中模块、轴系功耗模块、轴系振动模块，通过所述摩擦功耗、动刚度、阻尼系数、扰动信息将上述模块相关联，以轴承的径向变位及轴系转速作为优化变量，以轴系校中的质量作为约束条件，以轴系传动效率、支撑负荷、扰度曲线、临界转速、振动响应作为优化目标，进行推进系统多目标优化设计。
2.根据权利要求1所述的一种推进系统多目标优化设计方法，其特征在于，具体流程如下：
S1、建立输入参数库，包括轴系几何、工况参数；
S2、轴系校中模块中的静态校中模块计算得出各轴承负荷，并传递给轴系功耗模块中的轴承模块；
S3、所述轴承模块进行轴承摩擦功耗、动刚度、阻尼系数的计算，并将动刚度、阻尼系数传递至轴系校中模块中的动态校中模块，进行动态校中，获取最终轴承摩擦功耗、动刚度、阻尼系数、轴段扰度；
S4、根据最终的轴承动刚度和阻尼系数进行轴系回旋振动计算，得到轴颈响应值，传递给所述轴承模块，计算考虑振动的轴承摩擦功耗；
S5、所述轴系动态校中模块计算求得的轴段扰度，传递给轴系功耗模块中的转轴模块，计算求得转轴摩擦功耗；
S6、所述轴系功耗模块在完成轴承摩擦功耗、密封摩擦功耗、转轴摩擦功耗以及其它功耗计算的基础上，计算求得轴系传动效率；
S7、以轴承的径向变位及轴系转速作为优化变量，以轴系校中的质量作为约束条件，将轴系传动效率、支撑负荷、扰度曲线、临界转速、振动响应作为优化目标，进行推进系统多目标优化设计。
3.根据权利要求2所述的一种推进系统多目标优化设计方法，其特征在于，所述轴承摩擦功耗包括推力轴承摩擦功耗、中间轴承摩擦功耗和艉轴承摩擦功耗，所述中间轴承摩擦功耗和艉轴承摩擦功均为径向轴承摩擦功耗。
4.根据权利要求3所述的一种推进系统多目标优化设计方法，其特征在于，所述推力轴承摩擦功耗的计算包括如下步骤：
S1、根据推力轴承结构、工况参数以及轴承初始支点膜厚，求解膜厚方程；
S2、根据膜厚方程计算结果，求解稳态雷诺方程；
S3、判断计算出的轴承压力是否收敛，若是，则计算液膜水平力和垂直力，若否，则返回执行步骤S2；
S4、计算液膜合力；
S5、判断轴承载荷与液膜合力是否相等，若是，则执行步骤S6，若否，则修正支点膜厚并返回执行步骤S1；
S6、求解静特性方程，并输出推力轴承摩擦功耗计算结果。
5.根据权利要求3所述的一种推进系统多目标优化设计方法，其特征在于，所述径向轴承摩擦功耗计算包括如下步骤：
S1、根据径向轴承结构、工况参数以及初始偏心率和偏位角，求解膜厚方程；
S2、根据膜厚方程计算结果，求解稳态雷诺方程；
S3、判断计算出的轴承压力是否收敛，若是，则计算液膜水平力和垂直力，若否，则返回执行步骤S2；
S4、判断偏位角是否收敛，若是，则计算液膜合力，若否，则修正偏位角并返回执行步骤S1；
S5、判断轴承载荷与液膜合力是否相等，若是，则执行步骤S6，若否，则修正偏心率并返回执行步骤S1；
S6、求解静特性方程和瞬态雷诺方程；
S7、求解径向轴承动刚度和阻尼系数；
S8、输出径向轴承摩擦功耗计算结果。
6.根据权利要求2所述的一种推进系统多目标优化设计方法，其特征在于，所述密封摩擦功耗计算包括如下步骤：
S1、根据密封件的结构和转轴工况参数建立密封结构模型；
S2、对密封结构模型进行网格划分；
S3、进行有限元计算；
S4、判断有限元计算结果是否收敛，若是，则计算摩擦力矩并输出计算结果，若否，则返回执行步骤S3；
S5、对多工况点的摩擦力矩计算结果进行数学归纳；
S6、判断数学归纳结果误差是否在允许范围内，若是，则输出密封摩擦功耗计算结果，若否，则返回执行步骤S5。
7.根据权利要求2所述的一种推进系统多目标优化设计方法，其特征在于，所述转轴摩擦功耗计算包括如下步骤：
S1、根据敞水中轴段轴系结构和工况参数建立敞水中轴段流场模型；
S2、对敞水中轴段流场模型进行网格划分；
S3、进行Fluent流体有限元计算；
S4、判断流体有限元计算结果残差值是否收敛，若是，则计算摩擦力矩并输出计算结果，若否，则返回执行步骤S2；
S5、对多工况点的摩擦力矩计算结果进行数学归纳；
S6、判断数学归纳结果误差是否在允许范围内，若是，则输出转轴摩擦功耗计算结果，若否，则返回执行步骤S5。
8.根据权利要求2所述的一种推进系统多目标优化设计方法，其特征在于，所述其他功耗包括隔音联轴器、隔舱密封装置、遥测功率仪的功率损失，通过加入与所述其他功耗相匹配的修正系数对轴系功率损失进行修正。
9.根据权利要求2所述的一种推进系统多目标优化设计方法，其特征在于，所述轴系几何、工况参数具体包括轴系几何参数、船体变形、螺旋桨激励以及轴系转速。
10.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于：所述处理器执行所述程序时实现如权利要求1至9中任一项所述的一种推进系统多目标优化设计方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111130307.XA CN113806987B (zh) | 2021-09-26 | 2021-09-26 | 一种推进系统多目标优化设计方法及电子设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111130307.XA CN113806987B (zh) | 2021-09-26 | 2021-09-26 | 一种推进系统多目标优化设计方法及电子设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113806987A true CN113806987A (zh) | 2021-12-17 |
| CN113806987B CN113806987B (zh) | 2023-09-26 |
Family
ID=78938674
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111130307.XA Active CN113806987B (zh) | 2021-09-26 | 2021-09-26 | 一种推进系统多目标优化设计方法及电子设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113806987B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115876460A (zh) * | 2023-03-08 | 2023-03-31 | 山东双力现代农业装备有限公司 | 一种拖拉机传动性能测试方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108073755A (zh) * | 2017-05-25 | 2018-05-25 | 烟台仙崴机电有限公司 | 电动车开关磁阻电机系统多目标优化设计方法 |
| CN108416159A (zh) * | 2018-03-22 | 2018-08-17 | 中国人民解放军海军工程大学 | 一种舰船轴系优化方法及其优化平台 |
| CN110712731A (zh) * | 2019-10-22 | 2020-01-21 | 中国人民解放军海军工程大学 | 一种舰船轴系校中及回旋振动多学科优化方法 |
-
2021
- 2021-09-26 CN CN202111130307.XA patent/CN113806987B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108073755A (zh) * | 2017-05-25 | 2018-05-25 | 烟台仙崴机电有限公司 | 电动车开关磁阻电机系统多目标优化设计方法 |
| CN108416159A (zh) * | 2018-03-22 | 2018-08-17 | 中国人民解放军海军工程大学 | 一种舰船轴系优化方法及其优化平台 |
| CN110712731A (zh) * | 2019-10-22 | 2020-01-21 | 中国人民解放军海军工程大学 | 一种舰船轴系校中及回旋振动多学科优化方法 |
Non-Patent Citations (2)
| Title |
|---|
| 段明德;姬晓晓;胡伟奇;刘康华;张峰涛;: "高速精密数控车床主轴箱多目标优化设计", 河南科技大学学报(自然科学版), no. 02, pages 4 - 5 * |
| 袁海良;邱俊;王智;王德伦;: "风电传动系统轴承多目标优化设计", 机械设计与制造工程, no. 10, pages 25 - 30 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115876460A (zh) * | 2023-03-08 | 2023-03-31 | 山东双力现代农业装备有限公司 | 一种拖拉机传动性能测试方法 |
| CN115876460B (zh) * | 2023-03-08 | 2023-05-09 | 山东双力现代农业装备有限公司 | 一种拖拉机传动性能测试方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113806987B (zh) | 2023-09-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112765748B (zh) | 机理-数据异构信息融合的旋转机械数字孪生建模方法 | |
| Zhang et al. | Propeller excitation of longitudinal vibration characteristics of marine propulsion shafting system | |
| EP2211160A1 (en) | Method for analyzing torsional vibration of power transmission system, analyzing device, analyzing program, and shafting device between engine and driven device | |
| CN102741847B (zh) | 用于经受氢致应力裂纹的海底装备的设计方法 | |
| Young et al. | Reliability-based design and optimization of adaptive marine structures | |
| Wang et al. | Effects of floating ring bearing manufacturing tolerance clearances on the dynamic characteristics for turbocharger | |
| CN105653781A (zh) | 一种复合材料螺旋桨空泡性能的计算方法 | |
| CN113806987A (zh) | 一种推进系统多目标优化设计方法及电子设备 | |
| Zou et al. | The hydroelastic analysis of marine propellers with consideration of the effect of the shaft | |
| Murawski | Thermal interaction between main engine body and ship hull | |
| Zeng et al. | Modeling of assembly deviation with considering the actual working conditions | |
| Huang et al. | Dynamical response of the shaft-bearing system of marine propeller shaft with velocity-dependent friction | |
| Liu et al. | A model accounting for stiffness weakening of curvic couplings under various loading conditions | |
| CN116702538A (zh) | 一种基于有限元理论的用于阻尼涂层叶片的动力学建模计算方法 | |
| Liu et al. | Parametric modelling of vibration response for high-speed gear transmission system | |
| Martelli et al. | A design framework for combined marine propulsion control systems: From conceptualisation to sea trials validation | |
| Chen et al. | A new method of micro-slip contact normal pressure modeling and its application in the dynamic analysis of the dovetail–tenon–mortise jointed blades | |
| CN110287570B (zh) | 一种船舶艏侧推器激励载荷转换方法 | |
| CN109388907B (zh) | 一种具有预设纵向振动动柔度的轴系的设计方法 | |
| CN116738589A (zh) | 适用于非连续承力结构系统的振动能量传递分析方法 | |
| Luu et al. | Automatic calculation of torsional vibrations on marine propulsion plant using marine two-stroke diesel engine: Algorithms and Software | |
| CN110909420B (zh) | 考虑轴承因素的船舶推进轴系校中迭代计算方法 | |
| Yu | General influence coefficient algorithm in balancing of rotating machinery | |
| Yang et al. | Numerical investigation of after stern tube bearing during ship turning maneuver | |
| CN106503375A (zh) | 一种基于cn群理论确定汽轮机转子临界转速的方法及系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |