[发明专利]一种布料仿真方法、系统、电子设备及存储介质

权利要求书

1.一种基于弹簧质点模型的各项异性材质的布料仿真方法，其特征在于，包括：

对初始状态下的布料进行离散化的三角形网格划分，获得所述布料上每个三角形顶点坐标；

将初始状态下的所述布料在二维纹理坐标系中展开铺平，获得每个所述三角形顶点的二维纹理坐标；

按照牛顿力学以及胡克定律，对每个所述三角形顶点进行弹簧质点形变计算，获得每个时间节点上的每个所述三角形顶点的新坐标；

根据所述三角形顶点的新坐标以及拓扑连接关系，计算出表面的法向量并进行渲染，实时仿真地绘制出动态形变可视化的效果；

其中，对每个所述三角形顶点进行质点弹簧迭代形变计算包括每个所述三角形顶点中任意两个顶点之间弹力，所述弹力的强度系数包括位于所述二维纹理坐标系中X轴和Y轴方向；

对每个所述三角形顶点进行弹簧质点形变计算，获得每个时间节点上的每个所述三角形顶点的新坐标步骤包括：

根据弹性材质方程公式：

其中，p为所有所述三角形顶点坐标；t为当前时刻；t^-为上一时刻；p^[t]为所述布料在仿真的时刻t所有三角形顶点位置；为所述布料在上一时刻所有三角形顶点位置；h＝t-t^-为仿真时间；为所述布料在上一时刻三角形顶点的速度；M为仿真中标准的质量矩阵；W(p)为所述布料三角形顶点间的内力；

所述W(p)包括弹簧拉力E_s(p)和弯曲力E_b(p)；

所述弹簧拉力E_s(p)的拉力强度系数包括所述二维纹理坐标系中X轴和Y轴方向的拉力强度系数为所述弹簧拉力E_s(p)的弯曲强度系数包括所述二维纹理坐标系中X轴和Y轴方向的弯曲强度系数为

任意一个所述三角形中的两个顶点i、j形成边e_ij，所述两个顶点i、j在所述二维纹理坐标系形成边u_ij；

所述边u_ij与所述二维纹理坐标系中X轴的夹角为θ，且弹簧沿不同方向的拉伸量落在一个椭圆，可得所述边e_ij的拉力强度系数为和弯曲强度系数为

另外，所述边u_ij与所述二维纹理坐标系中X轴和Y轴的曲率分别为H^(x)，H^(y)；

由欧拉公式可知，沿所述边u_ij方向的曲率H^(ij)：

H^(ij)＝H^(x)(cosθ)²+H^(y)(sinθ)²；

由胡克定律可知，则

其中，Q′为所述三角形顶点与其对应所述二维纹理坐标轴位置累加得到的全局矩阵；p^T为所有所述三角形顶点组成的列向量；

则W(p)＝E_s(p)+E_b(p)。

2.如权利要求1所述的基于弹簧质点模型的各项异性材质的布料仿真方法，其特征在于，对初始状态下的布料进行离散化的三角形网格划分，获得所述布料上每个三角形顶点坐标的步骤包括：

对所述布料的不同面分别离散化并缝合，形成由三角形面片组成的离散化表示；

采用改进误差扩散算法，在离散化曲面上按一定的密度进行初始布点；

采用基于空间距离的粒子动力松弛算法对点云进行初步均匀化；

采用基于曲面距离的k均值算法进行再次均匀化；

对均匀的点云求基于曲面距离的Voronoi图，并获得三角网格；

对三角网格进行拓扑优化；

采用弹簧质点法对三角网格进行光顺优化。