CN104156520B - 一种基于直线投影的凸多面体碰撞检测方法 - Google Patents

Abstract

本发明涉及一种基于直线投影的凸多面体碰撞检测方法，步骤1：输入两凸多面体顶点集；步骤2：输入两凸多面每个面的顶点索引；步骤3：利用中心线正投影分离检测方法检测两凸多面体是否分离，如果分离，转步骤10，否则，进入步骤4；步骤4：计算两凸多面体的面摩擦值；步骤5：按面摩擦值降序原则构造相向面集合；步骤6：将相向面集中的棱边取出，构造投影分离线簇，投影分离线簇个数为λ，令k＝0；步骤7：如果k＝λ,则转步骤10，否则令k＝k+1，将凸多面体沿第k条棱在坐标平面做投影；步骤8：提取投影边界获得凸多边形；步骤9：判断凸多边形是否相交，如果分离，转步骤10，否则，转步骤7；步骤10：输出分离或碰撞结果。

Description

一种基于直线投影的凸多面体碰撞检测方法

技术领域

本发明涉及一种基于直线投影的凸多面体碰撞检测方法。

背景技术

碰撞检测在计算机图形学、虚拟装配、系统仿真、机器人路径设计等领域都有广泛的应用。由于任意物体均可表示成凸多面体的组合，因此实时准确的凸多面体碰撞检测方法对于复杂环境下多个物体的碰撞检测具有重要意义。

现有凸多面体碰撞检测方法多采用包围盒法、基于图形硬件的方法、线性规划法或分离面法等。现有凸多面体碰撞检测方法存在的不足之处在于：一是不能满足任意凸多面体碰撞检测的要求，二是碰撞检测过程中不能实时准确的给出检测结果，例如分离面法，通过凸多面体向二维坐标平面投影的方法，将3D空间凸多面体的相交问题转换为2D平面凸多边形的相交问题，但是该方法的前提是凸多面体的某个面恰能作为分离面，但当发生碰撞的点均位于各自凸多面体的某棱边，该算法失效。三是计算复杂，计算量大，例如分离面法中，传统构造分离面法实现凸多面体分离，难度大、算法复杂、计算量大。

发明内容

本发明目的在于提供一种基于直线投影的凸多面体碰撞检测方法，可适用任意凸多面体碰撞检测，并且检测实时、快速、准确。

实现本发明目的技术方案：

一种基于直线投影的凸多面体碰撞检测方法，其特征在于:

步骤1：输入两凸多面体V_A,V_B顶点集:V_A(A₁,A₂,…A_n),V_B(B₁,B₂,…B_m)；

步骤2：输入两凸多面体V_A,V_B每个面的顶点索引,两凸多面体V_A,V_B的面数分别为j,l.；

步骤3：利用中心线正投影分离检测方法检测两凸多面体是否分离，如果分离，转步骤10，否则，进入步骤4；

步骤4：计算两凸多面体的面摩擦值；

步骤5：按面摩擦值降序原则构造相向面集合；

步骤6：将相向面集中的棱边取出，构造投影分离线簇，投影分离线簇个数为λ，令

k＝0；

步骤7：如果k＝λ,则转步骤10，输出碰撞结果，否则令k＝k+1，将凸多面体沿第k条棱在坐标平面做投影；

步骤8：提取投影边界获得凸多边形；

步骤9：判断步骤8获得的凸多边形是否相交，如果分离，转步骤10，否则，转步骤7；

步骤10：输出分离或碰撞结果。

步骤3中，若2个凸多面体V_A,V_B的中心分别为C_A,C_B，直线C_AC_B的参数方程为y＝(1-t)C_A+tC_B.分别过2个凸多面体的顶点P_i作线段C_AC_B的垂线，则垂足的参数为：t_i＝(1-t)C_A+tC_B，i＝1,2,…,n,n+1,…,n+m；

其中，当1≤i≤n时，P_i＝A_i；当n+1≤i≤n+m时，P_i＝B_i，设V_A,V_B在中心线向量C_AC_B上的投影区间为(t_Amin,t_Amax),(t_Bmin,t_Bmax),若t_Amax＜t_Bmin或t_Bmax＜t_Amin时，则判断两个凸多面体分离。

步骤4中，凸多面体V_A,V_B的每个面分别为^AF_i、^BF_i，^AF_i、^BF_i上的外单位法向量分别为^AN_i、^BN_i，两凸多面体间中心向量为C_AC_B或C_BC_A，则面^AF_i的摩擦值^AI_i为， ^BF_i的摩擦值^BI_i为，凸多面体V_A的面摩擦值集合为{^AI_i}＝{^AI_i|i＝1,2,…,j}，凸多面体V_B的面摩擦值集合为{^BI_i}＝{^BI_i|i＝1,2,…,l})。

步骤5中，按摩擦值降序从V_A,V_B面集{^AF_i},{^BF_i}中分别取出j₁,l₁个面，构成了相向面集{^ABF_i}。

步骤6中，将相向面集{^ABF_i}中的棱边取出，构造投影分离线簇{^ABE_i|i＝1,2,…,λ}。

步骤8中，具体包括以下步骤，

步骤8.1：计算凸多面体顶点沿空间直线在坐标面的投影；

步骤8.2：利用快速排序法将所有投影顶点横坐标由小到大排序；

步骤8.3：根据凸包上边界线斜率下降准则，构造上边界；

步骤8.4：根据凸包下边界线斜率上升准则，构造完下边界；

步骤8.5：将上、下边界连接最终形成凸多边形。

本发明具有的有益效果：

本发明以直线投影为基础，针对分离面算法提出改进，首先通过中心线上的正投影剔除不可能发生碰撞的凸多面体，其次通过寻找凸多面体上的棱边构造投影线分离线簇，联系了三维凸多面体和二维平面多边形以及一维线段的相交检测关系，通过有限次投影运算和相交检测即可得出凸多面体碰撞检测结果，将投影线分离法成功运用到二维凸多边形的碰撞检测和最近距离计算。与现有凸多面体碰撞检测方法相比，本发明具有可适用任意凸多面体碰撞检测，并且检测实时、快速、准确的优点。

本发明改进了以往的凸多面体碰撞检测方法，黎自强在《计算机辅助设计与图形学学报》中提出“凸多面体快速碰撞检测的投影分离算法”(简称投影分离面法)，该方法将某一凸多面体表面作为分离面，并且沿着该表面与坐标平面的交线方向投影，通过投影的相交性判断空间凸多面体相交性，但当两凸多面体处于棱相交，棱对点相交时，凸多面体表面无法作为分离面，因此该方法失效。本发明改进投影分离面法，提出投影分离线的概念，沿着投影分离线向坐标平面作投影可将空间凸多面体的碰撞检测降维为平面多边形碰撞检测问题，本发明支持将凸多面体棱线作为投影分离线并且给出构造投影分离线簇的方法，通过投影分离线簇实现降维相交测试，实现凸多面体碰撞检测。

本发明对于任意复杂形状和任意相对位置的凸多面体均能实时准确的给出碰撞检测结果，且成功将基于直线的投影分离算法应用到凸多边形的碰撞检测和最近距离求解。

本发明应用面广，针对任何物体均可由凸多面体组合近似，因此本发明提出的一种基于直线投影的凸多面体碰撞检测方法可以应用到机器人路径规划、卫星仓布局优化设计、虚拟装配等诸多领域。

附图说明

图1本发明基于直线投影的凸多面体碰撞检测方法流程图；

图2凸多边形边界形成示意图；

图3凸多面体投影线分离示意图；

图4凸多边形投影线分离示意图。

具体实施方式

如图1所示，本发明基于直线投影的凸多面体碰撞检测方法包括以下步骤：

步骤1：输入两凸多面体V_A,V_B顶点集:V_A(A₁,A₂,…A_n),V_B(B₁,B₂,…B_m)；

步骤2：输入两凸多面体V_A,V_B每个面的顶点索引,两凸多面体V_A,V_B的面数分别为j,l.；

步骤3：利用中心线正投影分离检测方法检测两凸多面体是否分离，如果分离，转步骤10，否则，进入步骤4；

利用中心线上的正投影判断2个凸多面体分离，加快分离检测。若2个凸多面体V_A,V_B的中心分别为C_A,C_B，直线C_AC_B的参数方程为y＝(1-t)C_A+tC_B.分别过2个凸多面体的顶点P_i作线段C_AC_B的垂线，则垂足的参数为：t_i＝(1-t)C_A+tC_B，i＝1,2,…,n,n+1,…,n+m；

其中，当1≤i≤n时，P_i＝A_i；当n+1≤i≤n+m时，P_i＝B_i，设V_A,V_B在中心线向量C_AC_B上的投影区间为(t_Amin,t_Amax),(t_Bmin,t_Bmax),若t_Amax＜t_Bmin或t_Bmax＜t_Amin时，则判断两个凸多面体分离。

步骤4：计算两凸多面体的面摩擦值；

凸多面体V_A,V_B的每个面分别为^AF_i、^BF_i，^AF_i、^BF_i上的外单位法向量分别为^AN_i、^BN_i，两凸多面体间中心向量为C_AC_B或C_BC_A，则面^AF_i的摩擦值 ^BF_i的摩擦值凸多面体V_A的面摩擦值集合为{^AI_i}＝{^AI_i|i＝1,2,…,j}，凸多面体V_B的面摩擦值集合为{^BI_i}＝{^BI_i|i＝1,2,…,l})。

步骤5：按面摩擦值降序原则构造相向面集合；

按摩擦值降序从V_A,V_B面集{^AF_i},{^BF_i}中分别取出j₁,l₁个面，构成了相向面集{^ABF_i}。

步骤6：将相向面集{^ABF_i}中的棱边取出，构造投影分离线簇{^ABE_i|i＝1,2,…,λ}，投影分离线簇个数为λ，应注意筛选，避免重复投影，减少计算量。令k＝0；

步骤7：如果k＝λ,则转步骤10，输出碰撞结果，否则令k＝k+1，选取第k条投影分离线l_k(a_k,b_k,c_k)，将凸多面体在所选分离线在坐标系上作投影。以凸多面体V_A为例，将V_A沿着投影分离线l_B(a_i,b_i,c_i)在坐标系上作投影，如图3所示。其中点A_i(x_i,y_i,z_i)在坐标系xoz平面上投影顶点为A′_i(x′_i,y′_i,z′_i)，计算公式为：

若a_i≠0，A_i沿l_B向坐标平面yoz作投影；若c_i≠0，A_i沿l_B向坐标平面xoy作投影。

步骤8：提取投影边界获得凸多边形；

由于凸多面体在坐标系上的投影区域边界为凸多边形，因此需要根据投影顶点集提取投影凸多边形，如图2所示，

步骤8.1：按照公式(1)计算凸多面体顶点沿空间直线在坐标面的投影，如图2中(a)图所示；

步骤8.2：利用快速排序法将所有投影顶点横坐标由小到大排序，此时同一横坐标可能对应多个顶点，如图2中(b)图所示；

步骤8.3：根据凸包上边界线斜率下降准则，如图2中(c)图所示,图示顶点构成第一条棱边AB，下一个顶点为C,位于棱AB的左侧，则与准则冲突并从上边集中剔除B，将C加入上边集，同理删除D保留E，如果同一横坐标对应多个顶点，则选取纵坐标值最大的点。最终构造上边界，如图2中(d)图所示。

步骤8.4：根据凸包下边界线斜率上升准则，构造完下边界，如图2中(e)图所示；

步骤8.5：将上、下边界连接最终形成凸多边形,如图2中(f)图所示

步骤9：判断步骤8获得的凸多边形是否相交，如果分离，转步骤10，否则，转步骤7

如图4所示，将凸多面体投影分离法应用到凸多边形，将凸多边形相向面上的棱边作为投影分离线簇，再将凸多边形沿投影分离线向坐标轴作投影(如果投影分离线与所投影坐标轴平行，则改向另一坐标轴作投影)即可由投影区域π_A,π_B的相交性得出凸多边形是否相交，由于凸多边形的投影分离原理与凸多面体类似，因此不再详述。

当凸多边形分离时，凸多边形P_A,P_B间的最短距离h′与投影区域的最短距离h″满足：

h′＝h″|cosβ|

其中，β为投影分离线l_A与坐标轴的夹角，如图4所示。若h″＞0，凸多边形分离，可判定凸多面体分离

步骤10：输出分离或碰撞结果。

Claims (6)

1.一种基于直线投影的凸多面体碰撞检测方法，其特征在于:

步骤1：输入两凸多面体V_A,V_B顶点集:V_A(A₁,A₂,…A_n),V_B(B₁,B₂,…B_m)；

步骤2：输入两凸多面体V_A,V_B每个面的顶点索引,两凸多面体V_A,V_B的面数分别为j,l；

步骤3：利用中心线正投影分离检测方法检测两凸多面体是否分离，如果分离，转步骤10，否则，进入步骤4；

步骤4：计算两凸多面体的面摩擦值；

步骤5：按面摩擦值降序原则构造相向面集合；

步骤6：将相向面集中的棱边取出，构造投影分离线簇，投影分离线簇个数为λ，令

k＝0；

步骤7：如果k＝λ,则转步骤10，输出碰撞结果，否则令k＝k+1，将凸多面体沿第k条棱在坐标平面做投影；

步骤8：提取投影边界获得凸多边形；

步骤9：判断步骤8获得的凸多边形是否相交，如果分离，转步骤10，否则，转步骤7；

步骤10：输出分离或碰撞结果。

2.根据权利要求1所述的基于直线投影的凸多面体碰撞检测方法，其特征在于：步骤3中，若2个凸多面体V_A,V_B的中心分别为C_A,C_B，直线C_AC_B的参数方程为y＝(1-t)C_A+tC_B，分别过2个凸多面体的顶点P_i作线段C_AC_B的垂线，则垂足的参数为：t_i＝(1-t)C_A+tC_B，i＝1,2,…,n,n+1,…,n+m；

其中，当1≤i≤n时，P_i＝A_i；当n+1≤i≤n+m时，P_i＝B_i，设V_A,V_B在中心线向量C_AC_B上的投影区间为(t_Amin,t_Amax),(t_Bmin,t_Bmax),若t_Amax＜t_Bmin或t_Bmax＜t_Amin时，则判断两个凸多面体分离。

3.根据权利要求2所述的基于直线投影的凸多面体碰撞检测方法，其特征在于：步骤4中，凸多面体V_A,V_B的每个面分别为^AF_i、^BF_i，^AF_i、^BF_i上的外单位法向量分别为^AN_i、^BN_i，两凸多面体间中心向量为C_AC_B或C_BC_A，则面^AF_i的摩擦值^AI_i为， ^BF_i的摩擦值^BI_i为，凸多面体V_A的面摩擦值集合为{^AI_i}＝{^AI_i|i＝1,2,…,j}，凸多面体V_B的面摩擦值集合为{^BI_i}＝{^BI_i|i＝1,2,…,l}。

4.根据权利要求3所述的基于直线投影的凸多面体碰撞检测方法，其特征在于：步骤5中，按摩擦值降序从V_A,V_B面集{^AF_i},{^BF_i}中分别取出j₁,l₁个面，构成了相向面集{^ABF_i}。

5.根据权利要求4所述的基于直线投影的凸多面体碰撞检测方法，其特征在于：步骤6中，将相向面集{^ABF_i}中的棱边取出，构造投影分离线簇{^ABE_i|i＝1,2,…,λ}。

6.根据权利要求5所述的基于直线投影的凸多面体碰撞检测方法，其特征在于：步骤8中，具体包括以下步骤，

步骤8.1：计算凸多面体顶点沿空间直线在坐标面的投影；

步骤8.2：利用快速排序法将所有投影顶点横坐标由小到大排序；

步骤8.3：根据凸包上边界线斜率下降准则，构造上边界；

步骤8.4：根据凸包下边界线斜率上升准则，构造完下边界；

步骤8.5：将上、下边界连接最终形成凸多边形。