发明名称 基于二分法进行动力定位能力分析的方法
摘要 本发明提出了一种基于二分法进行动力定位能力分析的方法,包括步骤:制定二分判断函数I;确定风速作用区间[a,b];设定求解容许误差ε为1m/s;计算并实时更新风速作用区间;总体寻优步数估计。本发明能够得到海洋结构物在不同艏向角下的所抵抗的最大环境力(一般由最大风速表示),简化求解程序,提高求解效率,并且相较于传统方法,具有更好的稳定性。
申请公布号 CN104050370A 申请公布日期 2014.09.17
申请号 CN201410269891.0 申请日期 2014.06.17
申请人 上海交通大学 发明人 王磊;徐胜文;金鑫;汪学锋;邱荷珍
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 上海新天专利代理有限公司 31213 代理人 张宁展
主权项 一种基于二分法进行动力定位能力分析的方法,其特征在于,该方法包括如下步骤:步骤1,制定二分判断函数I:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>I</mi><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn><mo>,</mo></mtd><mtd><mi>x</mi><mo>&Element;</mo><mi>A</mi></mtd></mtr><mtr><mtd><mn>0</mn><mo>,</mo></mtd><mtd><mi>x</mi><mo>&NotElement;</mo><mi>A</mi></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000522363960000011.GIF" wi="446" he="154" /></maths>其中,x∈A表示在所求解的区域内存在可行解,此时,I=1,表示最终可承受的风速会大于现在的计算风速;<img file="FDA0000522363960000015.GIF" wi="152" he="66" />表示在所求解的区域内不存在可行解,此时I=0,表示最终可承受的风速会不大于现在的计算风速;步骤2,确定风速作用区间[a,b],设风速作用区间的初始选值为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><mi>a</mi><mo>=</mo><mi>&alpha;</mi><mo>&CenterDot;</mo><msubsup><mi>V</mi><mi>max</mi><mn>90</mn></msubsup></mtd></mtr><mtr><mtd><mi>b</mi><mo>=</mo><mi>&beta;</mi><mo>&CenterDot;</mo><msubsup><mi>V</mi><mi>max</mi><mn>180</mn></msubsup></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000522363960000012.GIF" wi="368" he="198" /></maths>其中,系数α和β的取值根据海洋结构物的特点进行选取;步骤3,设定容许误差ε为1m/s;步骤4,计算并实时更新风速作用区间,具体过程如下:每次计算,对风速作用区间[a,b]进行二分,并根据二分判断函数决定更新后的风速作用区间[a,b],进行后续计算,直到风力作用区间的模|b‑a|小于或等于容许误差ε,得出海洋结构物在该艏向角下所能承受的最大风速<img file="FDA0000522363960000013.GIF" wi="245" he="78" />依次计算其他艏向角度;步骤5,总体寻优步数估计:假设计算中的最初风速区间为[a<sub>0</sub>,b<sub>0</sub>],容许误差为ε<sub>0</sub>,则区间[a<sub>0</sub>,b<sub>0</sub>]内的计算点的总数为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>P</mi><mi>total</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>b</mi><mn>0</mn></msub><mo>-</mo><msub><mi>a</mi><mn>0</mn></msub></mrow><msub><mi>&epsiv;</mi><mn>0</mn></msub></mfrac></mrow>]]></math><img file="FDA0000522363960000014.GIF" wi="397" he="151" /></maths>在总体的计算点P<sub>total</sub>中进行二分,不断的接近最终的计算结果;实际上风力作用区间内的任意一个点表示为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>x</mi><mo>=</mo><mrow><mo>(</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>P</mi><mi>total</mi></msub><mo>&PlusMinus;</mo><mfrac><mn>1</mn><mn>4</mn></mfrac><msub><mi>P</mi><mi>total</mi></msub><mo>&PlusMinus;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&PlusMinus;</mo><mfrac><mn>1</mn><msup><mn>2</mn><mi>N</mi></msup></mfrac><msub><mi>P</mi><mi>total</mi></msub><mo></mo><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>&epsiv;</mi><mn>0</mn></msub></mrow>]]></math><img file="FDA0000522363960000021.GIF" wi="1067" he="141" /></maths>其中,寻优步数N满足2<sup>N‑1</sup><P<sub>total</sub>≤2<sup>N</sup>,寻优N步,计算点数为2<sup>N</sup>‑1最大计算步数可以通过下式进行估计:N=[log<sub>2</sub>P<sub>total</sub>]其中[.]表示取整。
地址 200240 上海市闵行区东川路800号