发明名称 引入偏微分方程求解移动机器人路径规划的方法
摘要 本发明公开一种引入偏微分方程求解移动机器人路径规划的方法,其包含:1、根据路径规划问题用建模偏微分方程模型;1.1用二维空间的辐射热传导方程建模;1.2、设置移动机器人路径规划模型的障碍物边界条件;2、采用自适应网格求解路径规划模型;2.1、初始化移动机器人可通行区域为拟正则三角剖分网格;2.2、判断网格精度是否正常;2.3求三角形中心解梯度;2.4、比较三角形中心解的梯度与预置能量最小值;2.5、细化网格;2.6、规划路径。本发明引入偏微分方程路径规划模型和设置障碍物边界条件,采用自适应网格求解此模型,进而得到规划结果;理论研究和计算结果表明该方法是可行、有效的,为移动机器人路径规划研究提供了新的思路和方法。
申请公布号 CN102708423A 申请公布日期 2012.10.03
申请号 CN201210183239.8 申请日期 2012.06.06
申请人 上海第二工业大学 发明人 康亮
分类号 G06Q10/04(2012.01)I 主分类号 G06Q10/04(2012.01)I
代理机构 上海信好专利代理事务所(普通合伙) 31249 代理人 张妍
主权项 1.一种引入偏微分方程求解移动机器人路径规划的方法,其特征在于,该方法包含以下步骤:步骤1、根据路径规划问题用建模偏微分方程模型;步骤1.1用二维空间的辐射热传导方程建模,在二维平面直角坐标系下,有以下方程:<img file="754777DEST_PATH_IMAGE001.GIF" wi="420" he="66" /><img file="6767DEST_PATH_IMAGE002.GIF" wi="441" he="68" />其中<img file="968907DEST_PATH_IMAGE003.GIF" wi="79" he="26" />分别表示电子、离子、光子;<img file="597334DEST_PATH_IMAGE004.GIF" wi="23" he="25" />是介质密度,<img file="462784DEST_PATH_IMAGE005.GIF" wi="73" he="28" />为比容;<img file="885676DEST_PATH_IMAGE006.GIF" wi="27" he="32" />为单位质量的内能,<img file="335111DEST_PATH_IMAGE007.GIF" wi="25" he="32" />为温度,<img file="767230DEST_PATH_IMAGE008.GIF" wi="23" he="25" />为压强,<img file="221607DEST_PATH_IMAGE009.GIF" wi="22" he="29" />为能量项,<img file="18662DEST_PATH_IMAGE010.GIF" wi="25" he="23" />为交换能;步骤1.2、设置移动机器人路径规划模型的障碍物边界条件;针对二维能量方程,做如下假设:(1)不同介质的界面固定不动,且<img file="689815DEST_PATH_IMAGE011.GIF" wi="302" he="63" />;(2)同一材料的介质密度为常数;(3)能源项为零;经过如上假设,对应偏微分定义后,得到如下定解条件:(1)固壁边界条件<img file="925624DEST_PATH_IMAGE012.GIF" wi="69" he="52" />;(2)障碍物边界的绝热条件<img file="998622DEST_PATH_IMAGE013.GIF" wi="68" he="51" />;(3)起始区域的Dirichlet条件<img file="530360DEST_PATH_IMAGE014.GIF" wi="74" he="30" />;(4)终点区域的Neumann条件<img file="688809DEST_PATH_IMAGE015.GIF" wi="94" he="54" />;(5)初始条件<img file="462730DEST_PATH_IMAGE016.GIF" wi="62" he="32" />;取<img file="390234DEST_PATH_IMAGE017.GIF" wi="257" he="34" />形式的控制函数,其中<img file="92873DEST_PATH_IMAGE018.GIF" wi="120" he="29" />;步骤2、采用自适应网格求解路径规划模型;设物理域<img file="738618DEST_PATH_IMAGE019.GIF" wi="22" he="22" />为工作环境中除去障碍物的移动机器人可通行区域;步骤2.1、移动机器人路径中采用的区域创建网格为三角形网格;设参考域<img file="112968DEST_PATH_IMAGE020.GIF" wi="33" he="34" />中一个边长为<img file="98242DEST_PATH_IMAGE021.GIF" wi="18" he="27" />的正三角形<img file="971782DEST_PATH_IMAGE022.GIF" wi="76" he="26" />在变换<img file="104823DEST_PATH_IMAGE023.GIF" wi="20" he="32" />之下映射到物理域<img file="220546DEST_PATH_IMAGE024.GIF" wi="24" he="24" />中的仿射三角形<img file="60327DEST_PATH_IMAGE025.GIF" wi="67" he="29" />;令<img file="603303DEST_PATH_IMAGE026.GIF" wi="82" he="34" />,<img file="447807DEST_PATH_IMAGE027.GIF" wi="87" he="35" />,设<img file="367222DEST_PATH_IMAGE028.GIF" wi="65" he="29" />和<img file="123825DEST_PATH_IMAGE029.GIF" wi="81" he="28" />面积分别为<img file="775386DEST_PATH_IMAGE030.GIF" wi="62" he="39" />和<img file="883020DEST_PATH_IMAGE031.GIF" wi="80" he="43" />,则<img file="842011DEST_PATH_IMAGE032.GIF" wi="113" he="67" />,<img file="453121DEST_PATH_IMAGE033.GIF" wi="114" he="66" />,<img file="603479DEST_PATH_IMAGE034.GIF" wi="442" he="74" />;初始化移动机器人可通行区域为拟正则三角剖分网格,设<img file="870513DEST_PATH_IMAGE035.GIF" wi="17" he="25" />为预置的均匀网格单元边长值,<img file="131730DEST_PATH_IMAGE036.GIF" wi="25" he="31" />-mini设为可通行区域能量预置的最小值;步骤2.2、在<img file="98811DEST_PATH_IMAGE037.GIF" wi="23" he="23" />的移动机器人通行的起始区域三角剖分网格上用有限元差分格式求得<img file="420071DEST_PATH_IMAGE036.GIF" wi="25" he="31" />的近似解;判断在随后的网格生成中三角形网格精度是否正常,<img file="174400DEST_PATH_IMAGE036.GIF" wi="25" he="31" />近似解的误差是否满足预先给定的容许数值;若是,则网格精度正常,跳转到步骤2.3;若否,则跳转到步骤2.2.1;步骤2.2.1、重新确定并生成该部分网格,跳转到步骤2.2,重新判断网格精度是否正常;步骤2.3在三角形结点和中点之间进行线性内插,求得三角形中心解的梯度;步骤2.4、判断三角形中心解的梯度是否小于预置能量最小值<img file="973729DEST_PATH_IMAGE036.GIF" wi="25" he="31" />-mini;若是,则说明移动机器人可通行区域全部网格生成结束,计算完成,线性内插结束,并跳转到步骤2.5;若否,则跳转到步骤2.3,继续进行下一步的梯度和网格参数计算;步骤2.5、细化网格;步骤2.6、集合三角形<img file="559431DEST_PATH_IMAGE036.GIF" wi="25" he="31" />-mini边,规划移动机器人可通行区域的路径。
地址 201206 上海市浦东新区金海路2360号