发明名称 一种水下滑翔机能耗最优的运动参数优化方法
摘要 本发明涉及一种水下滑翔机能耗最优的运动参数优化方法,输入水下滑翔机初始化参数;将运动参数约束条件数值中的最小值作为待优化运动参数的初始值;采用枚举计算方法分别计算出运动参数约束条件数值中所有滑翔角度条件下的G函数值;确定最大G函数值对应的滑翔角度;将滑翔速度增加一个迭代步长;如果超过最大滑翔速度约束条件,则运动参数约束条件数值中的滑翔深度增加一个迭代步长;如果超过最大滑翔深度约束条件,确定G函数值序列中的最大值,该G函数值序列中的最大值所对应的滑翔速度、滑翔深度和滑翔角度即为最优的滑翔运动参数。本发明综合优化了水下滑翔机的系统能耗特征和运动性能,提高水下滑翔机的能源利用率,提高水下滑翔机系统的续航能力。
申请公布号 CN102486633B 申请公布日期 2014.05.07
申请号 CN201010573970.2 申请日期 2010.12.04
申请人 中国科学院沈阳自动化研究所 发明人 俞建成
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 沈阳科苑专利商标代理有限公司 21002 代理人 李晓光
主权项 1.一种水下滑翔机能耗最优的运动参数优化方法,其特征在于包括以下步骤:输入水下滑翔机初始化参数,包括水下滑翔机基本特征参数、运动参数约束条件数值、以及运动参数迭代步长;将运动参数约束条件数值中的最小值作为待优化运动参数的初始值;根据水下滑翔机系统特征函数,采用枚举计算方法分别计算出运动参数约束条件数值中所有滑翔角度条件下的G函数值;根据计算的G函数值确定最大G函数值对应的滑翔角度,即为给定滑翔深度和滑翔速度条件下的最优滑翔角度;将滑翔速度增加一个迭代步长;判断增加一个迭代步长后的滑翔速度是否超过运动参数约束条件数值中的最大滑翔速度约束条件;如果超过最大滑翔速度约束条件,则运动参数约束条件数值中的滑翔深度增加一个迭代步长;判断增加一个迭代步长后的滑翔深度是否超过运动参数约束条件数值中的最大滑翔深度约束条件;如果超过最大滑翔深度约束条件,则根据上述最优滑翔角度对应的G函数值序列,确定G函数值序列中的最大值,该G函数值序列中的最大值所对应的滑翔速度、滑翔深度和滑翔角度即为最优的滑翔运动参数,结束本次优化过程;所述水下滑翔机系统特征函数为G函数:<maths num="0001"><![CDATA[<math><mrow><mi>G</mi><mrow><mo>(</mo><mi>U</mi><mo>,</mo><mi>h</mi><mo>,</mo><mi>&gamma;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mn>2</mn><mi>h</mi></mrow><mrow><mi>tan</mi><mo>|</mo><mi>&gamma;</mi><mo>|</mo><msub><mi>E</mi><mi>G</mi></msub><mrow><mo>(</mo><mi>U</mi><mo>,</mo><mi>h</mi><mo>,</mo><mi>&gamma;</mi><mo>)</mo></mrow></mrow></mfrac></mrow></math>]]></maths>式中,γ为滑翔角度;h为滑翔深度;U为滑翔速度;E<sub>G</sub>为水下滑翔机完成一个滑翔周期消耗的总能量;水下滑翔机完成一个滑翔周期消耗的总能量E<sub>G</sub>通过以下公式得到:<maths num="0002"><![CDATA[<math><mrow><mfenced open='' close=''><mtable><mtr><mtd><msub><mi>E</mi><mi>G</mi></msub><mrow><mo>(</mo><mi>U</mi><mo>,</mo><mi>h</mi><mo>,</mo><mi>&gamma;</mi><mo>)</mo></mrow><mo>=</mo><mo>|</mo><mfrac><mrow><msup><mrow><mn>2</mn><mi>U</mi></mrow><mn>2</mn></msup><mrow><mo>(</mo><msub><mi>K</mi><mrow><mi>L</mi><mn>0</mn></mrow></msub><mo>+</mo><msub><mi>K</mi><mi>L</mi></msub><mi>&alpha;</mi><mrow><mo>(</mo><mi>&gamma;</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow><mrow><mi>&rho;</mi><mi>g</mi><mi>cos</mi><mi>&gamma;</mi></mrow></mfrac><mo>|</mo><mrow><mo>(</mo><mfrac><msub><mi>P</mi><mi>v</mi></msub><msub><mi>q</mi><mi>v</mi></msub></mfrac><mo>+</mo><msub><mi>&eta;</mi><mi>bd</mi></msub><mrow><mo>(</mo><mi>h</mi><mo>)</mo></mrow><mi>h</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>+</mo><mn>4</mn><mfrac><msub><mi>mP</mi><mi>p</mi></msub><mrow><msub><mi>m</mi><mi>p</mi></msub><msub><mi>v</mi><mi>p</mi></msub></mrow></mfrac><mo>|</mo><msub><mi>z</mi><mi>G</mi></msub><mi>tan</mi><mrow><mo>(</mo><mi>&gamma;</mi><mo>+</mo><mi>&alpha;</mi><mrow><mo>(</mo><mi>&gamma;</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>|</mo><mo>+</mo><mfrac><mrow><msub><mrow><mn>2</mn><mi>P</mi></mrow><mi>c</mi></msub><mi>h</mi></mrow><mrow><mi>U</mi><mi>sin</mi><mo>|</mo><mi>&gamma;</mi><mo>|</mo></mrow></mfrac></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,α(γ)为滑翔运动攻角;K<sub>L0</sub>、K<sub>L</sub>为升力系数;P<sub>v</sub>为浮力调节装置在水面以流量q<sub>v</sub>工作时的功率;η<sub>bd</sub>(h)为浮力调节装置在滑翔深度h处工作时的效率;z<sub>G</sub>为水下滑翔机的稳心高;m为水下滑翔机的总重量;m<sub>p</sub>为俯仰调节装置移动质量;P<sub>p</sub>为俯仰调节装置以速度v<sub>p</sub>移动质量块时的功率;P<sub>c</sub>为控制系统作业过程中的平均功率。
地址 110016 辽宁省沈阳市东陵区南塔街114号