发明名称 一种新的起重机吊装动作规划方法
摘要 本发明属于工程机械、计算及控制自动化领域,涉及到利用快速扩展随机树和启发式策略以及基于全新的现代吊装方案设计理论方法的CAD技术——计算机辅助吊装方案设计技术,对履带式起重机进行吊装动作序列规划的方法;有效地解决了现阶段吊装动作规划方法所得路径迂回严重、路径代价大、规划时间长、路径质量差的问题,在障碍物随机分布的环境中,给出了起重机与被吊物在自动规避碰撞的前提下,从起始位姿自动转换到就位位姿所需的具体的动作序列,起重机依据这些动作序列能安全、高效地完成吊装作业。通过对吊装过程的仿真预演,说明了所提方法的有效性和高效性。
申请公布号 CN104484510A 申请公布日期 2015.04.01
申请号 CN201410713436.5 申请日期 2014.11.28
申请人 大连理工大学 发明人 吴迪;王晓丽;王欣
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 大连理工大学专利中心 21200 代理人 梅洪玉
主权项 一种起重机吊装动作规划方法,其特征在于以下步骤,(1)吊装系统位形建模由起重机和被吊物组成的系统为吊装系统,吊装系统的7个自由度为:①起重机下车坐标(x,z),取值范围分别为<img file="FDA0000619231420000011.GIF" wi="260" he="91" />和<img file="FDA0000619231420000012.GIF" wi="272" he="93" />②下车方向角α([‑π,π));③转台回转角β([‑π,π));④臂架仰角γ([0,π/2]);⑤起升绳长h,范围由臂长和臂架仰角决定,表示为[h<sub>lmt</sub>,h<sub>max</sub>(γ)],其中h<sub>lmt</sub>为起重机的限位高度,h<sub>max</sub>(γ)为起升滑轮组到地面的高度;⑥吊钩旋转角ω([‑π,π]);吊装系统每一个稳定的位姿都表示为一个向量X,如下式(1)所示;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>X</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>x</mi><mi>min</mi><mi>site</mi></msubsup></mtd><mtd><msubsup><mi>x</mi><mi>max</mi><mi>site</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>z</mi><mi>min</mi><mi>site</mi></msubsup></mtd><mtd><msubsup><mi>z</mi><mi>max</mi><mi>site</mi></msubsup></mtd></mtr><mtr><mtd><mo>-</mo><mi>&pi;</mi></mtd><mtd><mi>&pi;</mi></mtd></mtr><mtr><mtd><mo>-</mo><mi>&pi;</mi></mtd><mtd><mi>&pi;</mi></mtd></mtr><mtr><mtd><mi>O</mi></mtd><mtd><mi>&pi;</mi><mo>/</mo><mn>2</mn></mtd></mtr><mtr><mtd><msub><mi>h</mi><mi>lmt</mi></msub></mtd><mtd><msub><mi>h</mi><mi>max</mi></msub><mrow><mo>(</mo><mi>&gamma;</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>-</mo><mi>&pi;</mi></mtd><mtd><mi>&pi;</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000619231420000013.GIF" wi="1579" he="547" /></maths>两位姿X<sub>i</sub>和X<sub>i+1</sub>的距离是指被吊物从状态X<sub>i</sub>变换到X<sub>i+1</sub>所经历的轨迹长度,见下式(2):L(X<sub>i</sub>,X<sub>i+1</sub>)=|x<sub>i+1</sub>‑x<sub>i</sub>|+|z<sub>i+1</sub>‑z<sub>i</sub>|+|(r(a<sub>i+1</sub>‑a<sub>i</sub>)|+|(r(b<sub>i+1</sub>‑b<sub>i</sub>)|                                                           (2)+|(l<sub>b</sub>(g<sub>i+1</sub>‑g<sub>i</sub>)|+|(h<sub>i+1</sub>‑h<sub>i</sub>)|+|l<sub>w</sub>(w<sub>i+1</sub>‑w<sub>i</sub>)|其中,X<sub>i</sub>=(x<sub>i</sub>,z<sub>i</sub>,a<sub>i</sub>,b<sub>i</sub>,g<sub>i</sub>,h<sub>i</sub>,w<sub>i</sub>)和X<sub>i+1</sub>=(x<sub>i+1</sub>,z<sub>i+1</sub>,a<sub>i+1</sub>,b<sub>i+1</sub>,g<sub>i+1</sub>,h<sub>i+1</sub>,w<sub>i+1</sub>)分别为随机空间中2个独立的位姿,r为作业半径,l<sub>b</sub>为臂长,l<sub>w</sub>为被吊物长;设搜索到的动作序列Path包含的各节点依次为X<sub>0</sub>,X<sub>1</sub>,……,X<sub>n‑1</sub>,X<sub>n</sub>,则该条动作序列的代价为<img file="FDA0000619231420000014.GIF" wi="575" he="134" />其中,L(X<sub>i</sub>,X<sub>i+1</sub>)为点X<sub>i</sub>与点X<sub>i+1</sub>间的距离;(2)分析履带起重机的非完整性运动以XOZ为世界坐标系,O<sub>r</sub>为履带起重机的转弯中心,点P(x,z)为起重机的位置坐标,α为履带的方向角,v<sub>L</sub>、v<sub>R</sub>分别是左右履带的线速度,B为左右履带的轨距;起重机的行走和转弯所满足的运动学约束表示为式(3)和式(4),其中R为转弯半径:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><mover><mi>x</mi><mo>.</mo></mover><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mrow><mo>(</mo><msub><mi>v</mi><mi>R</mi></msub><mo>+</mo><msub><mi>v</mi><mi>L</mi></msub><mo>)</mo></mrow><mi>cos</mi><mi>&alpha;</mi></mtd></mtr><mtr><mtd><mover><mi>z</mi><mo>.</mo></mover><mo>=</mo><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mrow><mo>(</mo><msub><mi>v</mi><mi>R</mi></msub><mo>+</mo><msub><mi>v</mi><mi>L</mi></msub><mo>)</mo></mrow><mi>sin</mi><mi>&alpha;</mi></mtd></mtr><mtr><mtd><mrow><mover><mi>&alpha;</mi><mo>.</mo></mover><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>R</mi></mrow></mfrac></mrow><mrow><mo>(</mo><msub><mi>v</mi><mi>R</mi></msub><mo>+</mo><msub><mi>v</mi><mi>L</mi></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000619231420000021.GIF" wi="1496" he="533" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>R</mi><mo>=</mo><mfrac><mrow><mi>B</mi><mrow><mo>(</mo><msub><mi>v</mi><mi>R</mi></msub><mo>+</mo><msub><mi>v</mi><mi>L</mi></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>v</mi><mi>R</mi></msub><mo>-</mo><msub><mi>v</mi><mi>L</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000619231420000022.GIF" wi="1430" he="166" /></maths>当v<sub>L</sub>、v<sub>R</sub>同向并大小相等时,履带起重机直线行走;当v<sub>L</sub>、v<sub>R</sub>反向并大小相等时,履带起重机原地转向;当v<sub>L</sub>、v<sub>R</sub>同向并大小不等时,履带起重机以R为转弯半径为进行转向;起重机的动作由式(5)描述:u=[v<sub>L</sub>,v<sub>R</sub>,w<sub>β</sub>,w<sub>γ</sub>,v<sub>h</sub>,w<sub>ω</sub>]<sup>T</sup>    (5)其中w<sub>β</sub>、w<sub>γ</sub>、w<sub>ω</sub>分别为回转、变幅、吊钩旋转的角速度,v<sub>h</sub>为起升的线速度;起重机非完整性运动学约束的状态转移方程表达为式(6),通过积分得到一个新的位形q′,如式(7):<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mover><mi>q</mi><mo>.</mo></mover><mo>=</mo><mi>f</mi><mrow><mo>(</mo><mi>q</mi><mo>,</mo><mi>u</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>0.5</mn><mrow><mo>(</mo><msub><mi>v</mi><mi>R</mi></msub><mo>+</mo><msub><mi>v</mi><mi>L</mi></msub><mo>)</mo></mrow><mi>cos</mi><mi>&alpha;</mi></mtd></mtr><mtr><mtd><mo>-</mo><mn>0.5</mn><mrow><mo>(</mo><msub><mi>v</mi><mi>R</mi></msub><mo>+</mo><msub><mi>v</mi><mi>L</mi></msub><mo>)</mo></mrow><mi>sin</mi><mi>&alpha;</mi></mtd></mtr><mtr><mtd><mn>0.5</mn><mrow><mo>(</mo><msub><mi>v</mi><mi>R</mi></msub><mo>+</mo><msub><mi>v</mi><mi>L</mi></msub><mo>)</mo></mrow><mo>/</mo><mi>R</mi></mtd></mtr><mtr><mtd><msub><mi>w</mi><mi>&beta;</mi></msub></mtd></mtr><mtr><mtd><msub><mi>w</mi><mi>&gamma;</mi></msub></mtd></mtr><mtr><mtd><msub><mi>v</mi><mi>h</mi></msub></mtd></mtr><mtr><mtd><msub><mi>w</mi><mi>&omega;</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000619231420000031.GIF" wi="1574" he="686" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msup><mi>q</mi><mo>&prime;</mo></msup><mo>=</mo><mi>g</mi><mrow><mo>(</mo><mi>q</mi><mo>,</mo><mi>u</mi><mo>)</mo></mrow><mo>=</mo><mi>q</mi><mo>+</mo><mover><mi>q</mi><mo>.</mo></mover><mi>&Delta;t</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000619231420000032.GIF" wi="1418" he="96" /></maths>(3)起重机吊装动作规划为每个随机点X<sub>rand</sub>设置一个平滑度的约束,定义随机目标点X<sub>rand</sub>的平滑度为:R(X<sub>rand</sub>)=L(X<sub>init</sub>,X<sub>rand</sub>)+L(X<sub>rand</sub>,X<sub>goal</sub>)    (8)只有平滑度满足约束:R(X<sub>rand</sub>)≤ε·L(T<sub>init</sub>,T<sub>goal</sub>)(ε&gt;1)的随机点X<sub>rand</sub>才被采纳作为扩展的阶段目标点;同时,设置尝试次数Δ,即在尝试次数内选择平滑度满足约束的第一个点作为树扩展中的阶段目标点,若尝试次数超过Δ,则选择R(X<sub>rand</sub>)最小的点作为树扩展中的阶段目标点,保证树扩展的随机性。
地址 124221 辽宁省盘锦市辽东湾新区大工路2号