三维动画形变算法(Mixed Finite Elements)

Posted shushen

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了三维动画形变算法(Mixed Finite Elements)相关的知识,希望对你有一定的参考价值。

  混合有限元方法通入引入辅助变量后可以将高阶微分问题变成一系列低阶微分问题的组合。在三维网格形变问题中,我们考虑如下泛函极值问题:

技术图片

  其中u: Ω0 → R3是变形体的空间坐标,上述泛函极值问题对应的欧拉拉格朗日方程就是双调和方程?2u = 0。

  通过引入额外变量v,我们可以将上述无约束高阶优化问题转变为带约束的低阶优化问题:

技术图片

  引入拉格朗日函数,并利用格林公式可以得到:

技术图片

  将变量vu,λ写成基函数的线性组合形式,并且基函数选择为分段线性帽函数:

技术图片

  将上式分别对变量viuiλi求偏导可以得到:

技术图片

  这里讨论两种边界条件情况:

  第一种是固定区域边界条件,就是说在变形体上存在区域Ωf,其对应的空间位置为uf

技术图片

  求解形式:

技术图片

  消除额外变量后得到:

技术图片

  

  第二种是固定曲线边界条件,就是说在变形体上存在曲线C,其满足:

技术图片

  求解形式:

技术图片

  消除额外变量后得到:

技术图片

 

技术图片

效果:

技术图片

技术图片

技术图片

本文为原创,转载请注明出处:http://www.cnblogs.com/shushen

 

 

参考文献:

[1] Jacobson A, Tosun E, Sorkine O, et al. Mixed Finite Elements for Variational Surface Modeling. symposium on geometry processing, 2010, 29(5): 1565-1574.

 

附录

基于三维网格的微分问题求解

  对于求解三维网格上带Dirichlet边界条件的泊松方程:

技术图片

  将上式写成等价的弱形式:

技术图片

其中v是任意的测试函数,但是其满足在边界处值为0。

  利用分部积分,可以将上式等式左边进行变化:

技术图片

  由于测试函数在边界处值为0,所以上式中等式右边第一项为0:

技术图片

  然后将待求解函数u和测试函数v表示成基函数的线性组合形式,并且基函数选择为分段线性帽函数:

技术图片

其中I是网格上所有的顶点集合。

  三维网格上分段线性帽函数如下图所示:

技术图片

  那么对于网格内部任意顶点j,我们都取测试函数v = Φj,这样可以形成一个方程:

技术图片

  将积分里的求和提取到积分外,得到:

技术图片

  由于边界处u的数值已知,那么可以将已知项移到等式右侧:

技术图片

  这样就形成了n个方程组,n为三维网格除边界外的顶点个数。

 

  对于求解三维网格上带Neumann边界条件的泊松方程:

技术图片

  同样转换成弱形式,并利用分部积分:

技术图片

  在Neumann边界条件下,我们选择边界处值不为0的测试函数,那么可以得到:

技术图片

 

Cotangent权重矩阵L

  将Φi定义为分段线性帽函数,那么矩阵L的具体数值就可以根据三维网格模型唯一确定了,另外矩阵L之所以称为Cotangent权重,是因为其具体表达式是一个三角函数形式。

  下面具体推导下矩阵L的形式,根据其原始定义表达式:

技术图片

  在每个三角片内∇Φi是恒定的,并且只有在顶点i相邻的周边三角片上是非零值的。对于一个三角片,∇Φi指向的方向与顶点i相对的底边ei垂直,数值大小为三角片在底边ei上高h的倒数,即:

技术图片

上式中||ei||是边ei的长度,A是三角片的面积。

技术图片

  对于相邻顶点ij,∇Φi与∇Φj的方向分别垂直于各自的底边eiej,它们之间的夹角记为θ,于是我们可以得到:

技术图片

  由于ΦiΦj只有在相邻三角片TαTβ内才同时为非零值,那么:

技术图片

 

Mass质量矩阵M

  矩阵M的原始定义表达式为:

技术图片

  由于Φi定义为分段线性帽函数,所以Mij只有当ij为三维网格上相邻顶点或者ij相同时其值才为非零。

技术图片

  根据积分规则:

技术图片

其中A是三角片T的面积,x1x2x3是三角片三条边的中点。

  由于Φi在每个三角片内是个简单的线性函数,那么上式可以进一步表示为:

技术图片

  于是可以得到矩阵M的表达式:

技术图片

  通常矩阵M可以用对角矩阵Md来近似代替,近似过程为:

技术图片

  矩阵Md的对角线元素之和等于三维网格Ω的表面积,一般有两种方式来计算对角线元素,分别为重心质量矩阵和Voronoi质量矩阵,如下图所示。

技术图片

以上是关于三维动画形变算法(Mixed Finite Elements)的主要内容,如果未能解决你的问题,请参考以下文章

第六十二天 形变 动画 与 表格

2018.12.15 2d形变,动画,表格,多行文本垂直居中

iOS形变之CATransform3D

iOS形变之CGAffineTransform

PPT里面怎样做变形动画

jquery动画函数里面可以跟一个回调函数,表示动画结束后执行的代码