在衣片缝合过程中,为了将不同的衣片缝在一起,我们在衣片对应缝合边上加载缝合力。在模型中,缝合力被定义成对应缝合点之间距离的线性函数。对两个缝合点ｐｉ,ｊ和ｑｉ,ｊ间的缝合力,可以按如下公式计算:ｆ缝合力(ｉ,ｊ)=ｃｓ ｄｉｓ(ｐｉ,ｊ,ｑｉ,ｊ) ｎｐｉ,ｊ-ｑｉ,ｊ式中ｃｓ为缝合力系数,该系数与织物的缝合性能有关,通常,较难变形的布料采用较大的缝合力系数;ｄｉｓ(ｐｉ,ｊ,ｑｉ,ｊ)表示两缝合点ｐｉ,ｊ和ｑｉ,ｊ间的距离;ｎｐｉ,ｊ-ｑｉ,ｊ表示从ｐｉ,ｊ点指向ｑｉ,ｊ点的单位方向矢量。为了获得较真实的仿真效果,我们在变形模型中考虑了衣片所受的重力。质点所受的重力可按如下公式计算:ｆ重力(ｉ,ｊ)=ｍｉ,ｊ ｇ式中ｍｉ,ｊ为质点ｐｉ,ｊ的质量。在弹簧质点模型中,唯一考虑的弹性内力是弹簧的弹性变形力,由于采用的是理想的弹簧质点系统,可以利用胡克(ｈｏｏｋｅ)定律来计算弹簧的弹性变形力:ｆ内力(ｉ,ｊ)=-∑(ｋ,ｌ)∈ｒｋ(ｐｉ,ｊｐｋ,ｌ-∥ｐｉ,ｊｐｋ,ｌ∥0ｐｉ,ｊｐｋ,ｌ∥ｐｉ,ｊｐｋ,ｌ∥)其中,ｋ是弹簧的弹性变形系数,ｒ是ｐ(ｉ,ｊ)邻点的集合,∥ｐｉ,ｊｐｋ,ｌ∥0表示质点ｐ(ｉ,ｊ)与质点ｐ(ｋ,ｌ)之间的原始距离,弹簧的弹性变形系数ｋ可以依据所选用织物的材料性能参数曲线确定。

   2.3织物变形模型的求解我们选择显式欧拉方法来求解织物变形模型。求解公式如下:ａｉ,ｊ(ｔ+△ｔ)=1ｍｉ,ｊｆｉ,ｊ(ｔ)ｖｉ,ｊ(ｔ+△ｔ)=ｖｉ,ｊ(ｔ)+△ｔａｉ,ｊ(ｔ+△ｔ)ｐｉ,ｊ(ｔ+△ｔ)=ｐｉ,ｊ(ｔ)+△ｔｖｉ,ｊ(ｔ+△ｔ)其中,ｆｉ,ｊ是质点ｐ(ｉ,ｊ)所受所有力的合力,ｍｉ,ｊ(ｔ)是质点ｐ(ｉ,ｊ)的质量,ａｉ,ｊ(ｔ)、ｖｉ,ｊ(ｔ)和ｐｉ,ｊ(ｔ)分别是质点ｐ(ｉ,ｊ)在时间ｔ的加速度,速度和位置。△ｔ是系统选定的时间步长。

   3基于ａａｂｂ树层次包围盒的碰撞检测

   3.1建立ａａｂｂ树一个碰撞体的ａａｂｂ被定义为包含该碰撞体,且边平行于坐标轴的最小六面体。因此,描述一个ａａｂｂ,仅需六个标量。在构造ａａｂｂ包围盒时,需沿着碰撞体局部坐标系统的轴向(ｘ,ｙ,ｚ)来构造,所以所有的ａａｂｂ包围盒具有一致的方向。

   ａａｂｂ树是基于ａａｂｂ的二叉树,按照由上至下的递归细分方式构造生成的。在每一次递归过程中,要求取最小的ａａｂｂ,需沿所选择的剖分面将碰撞体分为正负两半,并将所对应的原始几何元素(如三角面)分别归属正、负两边,整个递归过程类似于空间二叉剖分,只是每次剖分的对象是ａａｂｂ,而不是空间区域。递归细分一直要进行到每一个叶子节点只包容一个原始几何元素为止,所以具有ｎ个原始几何元素的ａａｂｂ树具有ｎ-1个非叶子节点和ｎ个叶子节点。对于剖分面的选择,在本文中,选择垂直ａａｂｂ的最长轴,且平分该轴的平面。经试验证明,这种方式,在大多数情况下的算法复杂度仅为ｏ(ｎｌｏｇｎ),较其它的剖分面选择方法有了极大的提高。至于原始几何元素的归属则应依据几何元素的重心ｐ在最长轴上的投影坐标。若投影坐标大于剖分面的坐标(ｍｉｄ),则在剖分面的正向,否则在负向,如图2所示。图2三角面归属负区域,因为其质心投影坐标小于剖分面的基准坐标

   3.2ａａｂｂ的相交判断ａａｂｂ间的相交测试比较简单,两个ａａｂｂ相交当且仅当它们在三个坐标轴上的投影区间均相交。通过投影,我们即将三维求交问题转化为一维求交问题。而对一维求交问题,我们则采用ｓａｔ(ｓｅｐａｒａｔｉｎｇａｘｅｓｔｅｓｔ)[2]法。因ｓａｔ无需求交计算,只需比较两个包围盒分别在三个轴向上投影的重叠情况,即可得出相交测试结果,非常简单。现以在一个轴向上的投影情况为例说明:图3ａａｂｂｓ在ｘ轴向相交判断。

   设ａ,ｂ为两包围盒,ｘ为投影轴,ｃａ,ｃｂ分别为ａ,ｂ的中心点,ｐａ,ｐｂ为点ｃａ,ｃｂ在ｘ上的投影。ｒａ,ｒｂ分别为包围盒ａ,ｂ在ｘ上的投影。若ｒａ+ｒｂ

  ｐａｐｂ,(如图3所示)则在轴向ｘ上ａ和ｂ不相交,反之在轴向ｘ上ａ和ｂ邻接或相53第5期高成英等:虚拟穿衣中织物模型的建立和碰撞检测的处理交。当包围盒ａ,ｂ在三条轴向上的投影均相交时,则ａ,ｂ相交。定义ａａｂｂ的六个最大最小值分别确定了它在三个坐标轴上的投影区间,因此ａａｂｂ间的相交测试最多只需六次比较运算,非常简单快速。

   3.3ａａｂｂ树的更新当衣片移动、旋转后,需要对ａａｂｂ进行更新,根据定义ａａｂｂ的6个最大最小值的组合,可以得到ａａｂｂ的8个顶点,对这8个顶点进行相应的旋转和平移变化,并根据变化后的顶点计算新的ａａｂｂ。当衣片发生变形时,需要重新计算ａａｂｂ树中发生变形了的叶结点的ａａｂｂ,再利用变形叶节点的新ａａｂｂ来重新计算它们父节点的ａａｂｂ。这种计算必须严格按照从下到上的方式进行。父节点ａａｂｂ的具体求法为:令(ｘｍａｘ1,ｘｍｉｎ1,ｙｍａｘ1,ｙｍ

上一页  [1] [2] [3] [4] 下一页