一种船舶减摇陀螺参数优化方法技术

一种船舶减摇陀螺参数优化方法，步骤是：一是根据船舶结构基本参数与横摇运动基本参数，建立横摇运动数学模型，根据单框架控制力矩陀螺力矩平衡方程，建立陀螺与船舶联合运动数学模型，分别计算减摇陀螺不工作与工作时船舶横摇角的标准偏差值，建立船舶减摇陀螺减摇率的表达式；二是综合考虑减摇陀螺减摇率、占地面积、整体质量以及电机驱动能耗选定目标函数，根据船舶舱室的安装空间大小、电力系统供应功率以及阻尼器加工工艺难度给定参数的约束条件，构建船舶减摇陀螺参数优化模型；三是将原始减摇陀螺参数作为初始值，通过最值求解优化算法对船舶减摇陀螺参数优化模型进行求解，获得减摇陀螺最优参数值。本方法能有效提高减摇陀螺的减摇效果。

A Method for Optimizing the Parameters of Ship Rolling Stabilization Gyroscope

A method for optimizing the parameters of ship anti-rolling gyroscope is presented. Firstly, according to the basic parameters of ship structure and rolling motion, a mathematical model of rolling motion is established. According to the moment balance equation of single-frame control moment gyroscope, a mathematical model of the combined motion of gyroscope and ship is established. The standard deviation of ship roll angle is calculated when the anti-rolling gyroscope is not working or working, and the ship roll angle is established. The expression of anti-rolling rate of ship anti-rolling gyroscope; the second is to select the objective function by considering the anti-rolling rate of ship anti-rolling gyroscope, the area occupied, the overall mass and the energy consumption of motor drive, and to construct the optimization model of ship anti-rolling gyroscope parameters according to the constraints given by the installation space of ship cabin, the power supply of power system and the processing difficulty of damper. The parameters of the gyroscope are taken as initial values, and the optimal parameters of the ship anti-rolling gyroscope are obtained by solving the optimization model of the ship anti-rolling gyroscope parameters through the optimization algorithm of the maximum solution. This method can effectively improve the anti-rolling effect of anti-rolling gyroscope.

【技术实现步骤摘要】
一种船舶减摇陀螺参数优化方法
本专利技术涉及船舶减摇
，更具体地说，是涉及一种船舶减摇陀螺参数优化方法。
技术介绍
船舶减摇陀螺是船舶减摇装置的一种，近年来，随着现代材料科学、制造技术、控制理论等学科的发展，船舶减摇陀螺重新回到了人们的视野，引发了许多学者的关注与研究。TristanPerez(2009)建立了减摇陀螺与船舶联合数学模型，设计进动状态控制器，通过反馈的进动角大小改变进动状态控制器中的阻尼系数，实现减摇。尤方俊(2012)建立了减摇陀螺与船舶联合数学模型，并针对某一船型进行仿真分析。中国专利申请号为201510724324.4的专利公开了“一种船用陀螺减摇器及其减摇陀螺转子系统”，该专利技术介绍了减摇陀螺的基本框架，并在进动力矩方向和减摇力矩方向上通过磁悬浮式轴承施加一定的与进动力矩、减摇力矩方向相反的磁作用力，减少轴承磨损，提高寿命。美国专利申请号为US6973847B2的专利公开了“抑制船舶横摇的陀螺稳定器”，该专利技术改进了船舶减摇陀螺的框架结构，并提出几种可行性控制方法。现阶段，国内外关于减摇陀螺的专利及论文多为减摇陀螺整体结构框架设计和控制方法研究，且现有研究表明，船舶结构与运动参数、陀螺转子与进动阻尼参数的变化均会导致减摇效果大幅度改变，但根据这些影响因素设计其最优参数的方法并未见提及。
技术实现思路
本专利技术的目的是针对现有技术中存在的问题和不足，提供一种船舶减摇陀螺参数优化方法。本专利技术的方法是根据不同船型的舱室可用空间、供电功率、结构强度等各不相同，进一步优化减摇陀螺转子、进动阻尼参数、电机转速等关键影响因素的参数值，进而在船舶各约束条件下实现更好的减摇效果。为达到上述目的，本专利技术采用如下技术方案予以实现。一种船舶减摇陀螺参数优化方法，包括如下步骤：(1)根据船舶结构基本参数与船舶横摇运动基本参数，建立船舶横摇运动数学模型，再根据单框架控制力矩陀螺力矩平衡方程，建立陀螺与船舶联合运动数学模型，分别计算减摇陀螺不工作与工作时船舶横摇角的标准偏差值，进而建立船舶减摇陀螺减摇率的表达式。(2)综合考虑减摇陀螺减摇率、占地面积、整体质量以及电机驱动能耗选定目标函数，根据船舶舱室的可安装空间大小、电力系统可供应功率以及阻尼器加工工艺难度给定参数的约束条件，进而构建船舶减摇陀螺参数优化模型。(3)将原始减摇陀螺参数作为初始值，通过最值求解优化算法对船舶减摇陀螺参数优化模型进行求解，进而获得减摇陀螺最优参数值。进一步优选的，所述步骤(1)的具体内容和方法步骤如下：(a)计算减摇陀螺不工作时船舶横摇角φb的标准偏差值。设船体自身转动惯量为Iφφ，附加质量转动惯量为Jφφ，船舶自身阻尼恢复系数为c1、c3，船舶自身恢复力矩系数为k1、k3、k5，船舶的排水量为D，船舶的横稳性高为h，波浪模型方程为α(t)，则建立船舶横摇数学模型为：式中φb为船舶横摇角。令则将上式转化为一阶微分方程组初值问题：式中：对该微分方程在时间t上进行微分方程求解，得如下数值格式：其中h表示为时间t的迭代步长，k表示为迭代次数，则得减摇陀螺不工作时船舶横摇角随φb时间t变换的表达式：φbk＝φb(k-1)+hxk减摇陀螺不工作时船舶横摇角φb的标准偏差S(φb)为：式中N为在时间t内的取样次数，i∈N，为船舶横摇角φb的平均值。(b)计算减摇陀螺工作时船舶横摇角φa的标准偏差值。船体自身转动惯量为Iφφ，附加质量转动惯量为Jφφ，船舶自身阻尼恢复系数为c1、c3，船舶自身恢复力矩系数为k1、k3、k5，船舶的排水量为D，船舶的横稳性高为h，波浪模型方程为α(t)，则建立减摇陀螺与船舶联合数学模型为：式中：，φb为船舶横摇角，β为减摇陀螺进动角，R为减摇陀螺转子半径，L为减摇陀螺转子厚度，ω为转子转速，M为减摇陀螺进动方向的阻尼系数。令则将上式转化为一阶微分方程组初值问题：式中：f1(t,φa(t),x(t),β(t),y(t))＝-c1x(t)-c3x(t)3-k1φa(t)-k3φa(t)3-k5φa(t)5-h0y(t)cos(β(t))/(Iφφ+Jφφ)+Dhα(t)/(Iφφ+Jφφ)f2(t,x(t),β(t),y(t))＝h0x(t)cos(β(t))/J-My(t)/J对该微分方程在时间t上进行微分方程求解，可得如下数值格式：其中h表示为时间t的迭代步长，k表示为迭代次数，则得减摇陀螺工作时船舶横摇角随φa时间t变换的表达式：φak＝φa(k-1)+hxk减摇陀螺工作时船舶横摇角φa的标准偏差S(φa)为；式中N为在时间t内的取样次数，i∈N，为船舶横摇角φa的平均值。(c)建立船舶减摇陀螺减摇率的表达式。船舶减摇陀螺减摇率为：减摇率：进一步优选的，所述步骤(2)的具体内容和方法步骤如下：(a)选定目标函数根据船舶自身结构与运动参数、舱室可安装空间、动力驱动系统供电功率、材料结构强度等关键因素提出评价指标，包含减摇率TT、减摇陀螺占地面积πR2、减摇陀螺整体质量ρπR2L、减摇陀螺驱动电机转动能耗ω、减摇陀螺进动阻尼系数M等，综合表示为：minZ(x1,x2,...,xi,...,xn)其中xi∈(R,L,M,ω,TT，…)(b)给定约束条件根据船舶舱室的可安装空间大小，给出减摇陀螺转子的尺寸约束条件。假设船舶舱室的可安装空间尺寸为h长×h宽×h高，则陀螺转子的半径R以及厚度L需满足以下约束条件：根据船舶电力系统可供应功率，给出电机转速的约束条件。假设可供应功率为P船，由于其中P为额定功率，T为额定转矩，则电机转速需满足以下约束条件：根据进动阻尼器加工工艺难度，考虑阻尼器工作时散热及强度等问题，给出阻尼器阻尼系数的约束条件，表示为：Z低≤M≤Z高，其中Z低与Z高需根据具体阻尼器大小给出具体值。对上述步骤(1)求出的减摇陀螺减摇率TT同样提出约束要求：TT≥a(0≤a≤1)则约束条件总结为：(c)构建船舶减摇陀螺参数优化模型根据上述的目标函数以及约束条件，建立船舶减摇陀螺参数优化模型：进一步的，构建减摇率最高的船舶减摇陀螺参数优化模型。要求减摇率最高，则取目标函数为：minZ(1-TT)考虑到船舶舱室可用安装空间、船舶电力系统可供电功率、阻尼器设计要求限制以及减摇率要求等问题给出约束条件。则约束条件表示为：建立减摇率最高的船舶减摇陀螺参数优化模型：进一步的，构建转子质量最小的船舶减摇陀螺参数优化模型。通常在优化设计时都会涉及到轻量化设计，即质量最小。故设计：当减摇率满足设定要求时，陀螺转子的质量最小，则目标函数为：minZ(ρπR2L)式ρ中表示密度，依据转子材料而定。考虑到船舶舱室可用安装空间、船舶电力系统可供电功率、阻尼器设计要求限制以及减摇率要求等问题给出约束条件。则约束条件表示为：则建立转子质量最小的船舶减摇陀螺参数优化模型：进一步的，构建转子直径最小的船舶减摇陀螺参数优化模型。考虑到船舶的空间有限，希望减摇陀螺的占地面积较小，对比减摇陀螺的基本结构可知转子的直径对减摇陀螺整体尺寸影响最大，故目标函数设为减摇陀螺的直径最小：minZ(2R)考虑到船舶舱室可用安装空间、船舶电力系统可供电功率、阻尼器设计要求限制、减摇率要求以及转子直径与厚度关系等问题给出约束条件。则约束条件表示为：则建立转子直径本文档来自技高网...

【技术保护点】
1.一种船舶减摇陀螺参数优化方法，其特征是，包括如下步骤：(1)根据船舶结构基本参数与船舶横摇运动基本参数，建立船舶横摇运动数学模型，再根据单框架控制力矩陀螺力矩平衡方程，建立陀螺与船舶联合运动数学模型，分别计算减摇陀螺不工作与工作时船舶横摇角的标准偏差值，进而建立船舶减摇陀螺减摇率的表达式；(2)综合考虑减摇陀螺减摇率、占地面积、整体质量以及电机驱动能耗选定目标函数，根据船舶舱室的可安装空间大小、电力系统可供应功率以及阻尼器加工工艺难度给定参数的约束条件，进而构建船舶减摇陀螺参数优化模型；(3)将原始减摇陀螺参数作为初始值，通过最值求解优化算法对船舶减摇陀螺参数优化模型进行求解，进而获得减摇陀螺最优参数值。

【技术特征摘要】
1.一种船舶减摇陀螺参数优化方法，其特征是，包括如下步骤：(1)根据船舶结构基本参数与船舶横摇运动基本参数，建立船舶横摇运动数学模型，再根据单框架控制力矩陀螺力矩平衡方程，建立陀螺与船舶联合运动数学模型，分别计算减摇陀螺不工作与工作时船舶横摇角的标准偏差值，进而建立船舶减摇陀螺减摇率的表达式；(2)综合考虑减摇陀螺减摇率、占地面积、整体质量以及电机驱动能耗选定目标函数，根据船舶舱室的可安装空间大小、电力系统可供应功率以及阻尼器加工工艺难度给定参数的约束条件，进而构建船舶减摇陀螺参数优化模型；(3)将原始减摇陀螺参数作为初始值，通过最值求解优化算法对船舶减摇陀螺参数优化模型进行求解，进而获得减摇陀螺最优参数值。2.根据权利要求1所述的船舶减摇陀螺参数优化方法，其特征是，步骤(1)的具体内容和方法步骤是：(a)计算减摇陀螺不工作时船舶横摇角φb的标准偏差值设船体自身转动惯量为Iφφ，附加质量转动惯量为Jφφ，船舶自身阻尼恢复系数为c1、c3，船舶自身恢复力矩系数为k1、k3、k5，船舶的排水量为D，船舶的横稳性高为h，波浪模型方程为α(t)，则建立船舶横摇数学模型为：式中φb为船舶横摇角；令则将上式转化为一阶微分方程组初值问题：式中：对该微分方程在时间t上进行微分方程求解，得如下数值格式：其中h表示为时间t的迭代步长，k表示为迭代次数，则得减摇陀螺不工作时船舶横摇角随φb时间t变换的表达式：φbk＝φb(k-1)+hxk减摇陀螺不工作时船舶横摇角φb的标准偏差S(φb)为：式中N为在时间t内的取样次数，i∈N，为船舶横摇角φb的平均值；(b)计算减摇陀螺工作时船舶横摇角φa的标准偏差值船体自身转动惯量为Iφφ，附加质量转动惯量为Jφφ，船舶自身阻尼恢复系数为c1、c3，船舶自身恢复力矩系数为k1、k3、k5，船舶的排水量为D，船舶的横稳性高为h，波浪模型方程为α(t)，则建立减摇陀螺与船舶联合数学模型为：式中：，φb为船舶横摇角，β为减摇陀螺进动角，R为减摇陀螺转子半径，L为减摇陀螺转子厚度，ω为转子转速，M为减摇陀螺进动方向的阻尼系数；令则将上式转化为一阶微分方程组初值问题：式中：f1(t,φa(t),x(t),β(t),y(t))＝-c1x(t)-c3x(t)3-k1φa(t)-k3φa(t)3-k5φa(t)5-h0y(t)cos(β(t))/(Iφφ+Jφφ)+Dhα(t)/(Iφφ+Jφφ)f2(t,x(t),β(t),y(t))＝h0x(t)cos(β(t))/J-My(t)/J对该微分方程在时间t上进行微分方程求解，可得如下数值格式：其中h表示为时间t的迭代步长，k表示为迭代次数，则得减摇陀螺工作时船舶横摇角随φa时间t变换的表达式：φak＝φa(k-1)+hxk减摇陀螺工作时船舶横摇角φa的标准偏差S(φa)为；式中N为在时间t内的取样次数，i∈N，为船舶横摇角φa的平均值；(c)建立船舶减摇陀螺减摇率的表达式船舶减摇陀螺减摇率为：减摇率：3.根据权利要求1所述的船舶减摇陀螺参数优化方法，其特征是，步骤(2)的具体内容和方法步骤是：(a)选定目标函数根据船舶自身结构与运动参数、舱室可安装空间、动力驱动系统供电功率、材料结构强度等关键因素提出评价指标，包含减摇率TT、...

【专利技术属性】
技术研发人员：苏世杰，钱育辰，丁晓阳，唐文献，丁郁华，祝圆圆，陈赟，张建，
申请(专利权)人：江苏科技大学，
类型：发明
国别省市：江苏,32