发明名称 带死区输入的海洋柔性立管边界控制方法
摘要 本发明公开了一种带死区输入的海洋柔性立管边界控制方法,对海洋柔性系统进行建模,得到柔性立管系统的控制方程,构建边界控制率计算公式,对于边界控制率未确定的控制增益等参数,通过预先训练得到,然后实际测得柔性立管顶端的横向位移、倾角和剪切力,并通过历史数据得到这三个参数在时刻t的导数,最后将实际参数代入边界控制率计算公式,得到时刻t的边界控制率,驱动装置根据边界控制率向柔性立管顶端施加作用力。本发明可以实现对海洋振动下柔性立管振动的有效抑制,并且测量参数少,易于实现。
申请公布号 CN104537144A 申请公布日期 2015.04.22
申请号 CN201410706873.4 申请日期 2014.11.27
申请人 电子科技大学 发明人 贺威;何修宇;聂双喜;秦晖
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 成都行之专利代理事务所(普通合伙) 51220 代理人 温利平
主权项 一种带死区输入的海洋柔性立管边界控制方法,其特征在于,包括以下步骤:S1:对海洋柔性立管系统进行建模,得到柔性立管系统的控制方程为:<img file="FDA0000618521140000011.GIF" wi="1267" he="91" />其中,以立管与油田的连接处为原点,x表示立管当前位置距离原点的垂直距离,x的取值范围为0≤x≤L,L表示海平面到原点的垂直距离;t表示时间,取值范围为t≥0;<img file="FDA0000618521140000012.GIF" wi="158" he="88" />表示时间t时立管x处的横向位移,<img file="FDA0000618521140000013.GIF" wi="166" he="92" />和<img file="FDA0000618521140000014.GIF" wi="156" he="92" />分别表示函数<img file="FDA0000618521140000015.GIF" wi="160" he="89" />对时间t的一阶导数和二阶导数,<img file="FDA0000618521140000016.GIF" wi="191" he="90" />和<img file="FDA0000618521140000017.GIF" wi="210" he="90" />分别表示函数<img file="FDA0000618521140000018.GIF" wi="164" he="90" />对距离x的二阶导数和四阶导数;ρ表示立管的密度;EI表示立管的弯曲刚度;Q表示张力;f(x,t)表示时变洋流对在时间t对立管x处的干扰作用力;c表示立管的阻尼系数;S2:构建边界控制率v计算公式为:<img file="FDA0000618521140000019.GIF" wi="1524" he="264" />其中,sgn()表示符号函数,<img file="FDA00006185211400000110.GIF" wi="182" he="86" />和<img file="FDA00006185211400000111.GIF" wi="196" he="85" />分别表示<img file="FDA00006185211400000112.GIF" wi="156" he="92" />在立管顶端L处的一阶导数和三阶导数,<img file="FDA00006185211400000113.GIF" wi="179" he="86" />表示顶端的倾角,<img file="FDA00006185211400000114.GIF" wi="202" he="90" />表示顶端的剪切力,<img file="FDA00006185211400000115.GIF" wi="372" he="86" />和<img file="FDA00006185211400000116.GIF" wi="200" he="87" />分别表示<img file="FDA00006185211400000117.GIF" wi="366" he="87" />和<img file="FDA00006185211400000118.GIF" wi="194" he="87" />对时间t的一阶导数;M<sub>s</sub>表示位于海平面处、与立管顶端连接的驱动装置的重量;d<sub>s</sub>表示驱动装置的阻尼系数;k<sub>1</sub>、k<sub>2</sub>和k表示控制增益参数;φ(t)=[φ<sub>l</sub>(t) φ<sub>r</sub>(t)]<sup>T</sup>,上标T表示矩阵转置,其中:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>&phi;</mi><mi>l</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mi>v</mi><mo>&lt;</mo><msub><mi>b</mi><mi>r</mi></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mi>b</mi><mi>r</mi></msub><mo>&le;</mo><mi>v</mi></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA00006185211400000119.GIF" wi="378" he="162" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>&phi;</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn></mtd><mtd><msub><mi>b</mi><mi>l</mi></msub><mo>&lt;</mo><mi>v</mi></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>v</mi><mo>&le;</mo><msub><mi>b</mi><mi>l</mi></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA00006185211400000120.GIF" wi="379" he="157" /></maths>F(t)=[F<sub>l</sub>(v) F<sub>r</sub>(v)]<sup>T</sup>,其中:<img file="FDA00006185211400000121.GIF" wi="683" he="165" /><img file="FDA0000618521140000021.GIF" wi="674" he="179" />其中,b<sub>l</sub>和b<sub>r</sub>为死区边界,上标“'”表示一阶求导,<img file="FDA00006185211400000215.GIF" wi="131" he="79" />表示以v为变量的函数,需要满足的条件为:<img file="FDA0000618521140000022.GIF" wi="585" he="170" /><img file="FDA0000618521140000023.GIF" wi="616" he="164" />u<sub>a</sub>(t)为辅助信号,定义为:<img file="FDA0000618521140000024.GIF" wi="875" he="94" /><img file="FDA0000618521140000025.GIF" wi="114" he="86" />根据<img file="FDA0000618521140000026.GIF" wi="434" he="106" />公式求得,其中<img file="FDA0000618521140000027.GIF" wi="105" he="92" />表示<img file="FDA0000618521140000028.GIF" wi="110" he="86" />对时间t的一阶导数,γ表示待定参数;参数k<sub>1</sub>、k<sub>2</sub>、k和γ通过预先训练得到,训练方法为:预先仿真一个柔性立管系统,设置边界控制率v计算公式中其他各个参数,根据边界控制率v计算公式和步骤S1中的控制方程进行仿真训练,其训练目标是令柔性立管上各处的横向位移均小于等于横向位移<img file="FDA0000618521140000029.GIF" wi="162" he="91" />的M%,0≤M<100;S3:在时刻t,采用位移传感器测量得到柔性立管顶端的横向位移<img file="FDA00006185211400000210.GIF" wi="192" he="90" />采用倾角计测量得到倾角<img file="FDA00006185211400000211.GIF" wi="202" he="92" />采用剪切力传感器测量得到剪切力<img file="FDA00006185211400000212.GIF" wi="225" he="92" />再通过历史数据得到三个参数在时刻t的导数<img file="FDA00006185211400000213.GIF" wi="684" he="93" />S4:将步骤S3得到的实际参数代入步骤S2得到的边界控制率计算公式,得到时刻t的边界控制率v(t),驱动装置根据边界控制率v(t)向柔性立管顶端施加作用力。
地址 611731 四川省成都市高新区(西区)西源大道2006号