![Page 1: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/1.jpg)
Robust Treatment of Degenerate Elements inInteractive Corotational FEM Simulations
O. Civit-Flores and A. Susın
UPC-BarcelonaTech
June 11, 2014
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 1 / 36
![Page 2: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/2.jpg)
Outline
1 Introduction
2 Corotational FEM
3 Rotation extraction
4 Degeneration-Aware Polar Decomposition
5 Results
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 2 / 36
![Page 3: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/3.jpg)
Interactive FEM
Interactive simulation of deformable solids using FEMApplications:
Virtual reality, surgery, training...Videogames
Requirements:User interactionEfficiencyRobustnessRealism
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 3 / 36
![Page 4: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/4.jpg)
Contribution
Element degeneration threatens robustness and realism:We identify issues with existing degenerate element treatmentschemesWe propose a new method that avoids them
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 4 / 36
![Page 5: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/5.jpg)
Outline
1 Introduction
2 Corotational FEM
3 Rotation extraction
4 Degeneration-Aware Polar Decomposition
5 Results
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 5 / 36
![Page 6: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/6.jpg)
Finite Element Method
Partition the computational domain Ω into sub-domains Ωi with Nshared nodes
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 6 / 36
![Page 7: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/7.jpg)
Linear FEM
Tetrahedral elements:
re =
r1...
r4
, xe =
x1...
x4
, f e =
f 1...
f 4
The elastic forces on the nodes are:
f e = −Keue, ue = xe − re
Properties:Constant stiffness matrix Ke
Invariant to translation, but not to rotation
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 7 / 36
![Page 8: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/8.jpg)
Linear FEMlinearization error
From [1]
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 8 / 36
![Page 9: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/9.jpg)
Corotational Linear FEM (I)
Idea: Apply LFEM in a reference system local to each element
R
r1
r2 r3
x2
x1
x3
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 9 / 36
![Page 10: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/10.jpg)
Corotational Linear FEM (II)
1 Compute deformation matrix F :
F = D(xe)D(re)−1, D(ve) =[v2 − v1 v3 − v1 v4 − v1
]2 Factorize into rotation and scaling:
F = RS
3 Apply linear elasticity in local coordinates:
f e = −ReKe(RTe xe − re)
Properties:Geometrically non-linearInvariant to translation and rotation
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 10 / 36
![Page 11: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/11.jpg)
Corotational Linear FEM (II)
1 Compute deformation matrix F :
F = D(xe)D(re)−1, D(ve) =[v2 − v1 v3 − v1 v4 − v1
]2 Factorize into rotation and scaling:
F = RS
3 Apply linear elasticity in local coordinates:
f e = −ReKe(RTe xe − re)
Properties:Geometrically non-linearInvariant to translation and rotation
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 10 / 36
![Page 12: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/12.jpg)
Corotational Linear FEM (II)
1 Compute deformation matrix F :
F = D(xe)D(re)−1, D(ve) =[v2 − v1 v3 − v1 v4 − v1
]2 Factorize into rotation and scaling:
F = RS
3 Apply linear elasticity in local coordinates:
f e = −ReKe(RTe xe − re)
Properties:Geometrically non-linearInvariant to translation and rotation
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 10 / 36
![Page 13: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/13.jpg)
Corotational Linear FEM (II)
1 Compute deformation matrix F :
F = D(xe)D(re)−1, D(ve) =[v2 − v1 v3 − v1 v4 − v1
]2 Factorize into rotation and scaling:
F = RS
3 Apply linear elasticity in local coordinates:
f e = −ReKe(RTe xe − re)
Properties:Geometrically non-linearInvariant to translation and rotation
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 10 / 36
![Page 14: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/14.jpg)
Corotational Linear FEM (II)
1 Compute deformation matrix F :
F = D(xe)D(re)−1, D(ve) =[v2 − v1 v3 − v1 v4 − v1
]2 Factorize into rotation and scaling:
F = RS
3 Apply linear elasticity in local coordinates:
f e = −ReKe(RTe xe − re)
Properties:Geometrically non-linearInvariant to translation and rotation
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 10 / 36
![Page 15: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/15.jpg)
Dynamics and Quasistatics
Node positions x ∈ R3N are the DOF:Dynamics:
Mx = f s(x) + f d (x , x) + f ext
Quasistatics:f s(x) = −f ext
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 11 / 36
![Page 16: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/16.jpg)
Element Degeneration
Collapse with |det(F)| < ε or inversion with det(F) < 0UnphysicalUnavoidable with (finite) linear forcesUnavoidable due to discretizationUnavoidable due to user interactiondet(F) < ε affects F = RS factorization
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 12 / 36
![Page 17: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/17.jpg)
Element Degeneration
Collapse with |det(F)| < ε or inversion with det(F) < 0UnphysicalUnavoidable with (finite) linear forcesUnavoidable due to discretizationUnavoidable due to user interactiondet(F) < ε affects F = RS factorization
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 12 / 36
![Page 18: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/18.jpg)
Element Degeneration
Collapse with |det(F)| < ε or inversion with det(F) < 0UnphysicalUnavoidable with (finite) linear forcesUnavoidable due to discretizationUnavoidable due to user interactiondet(F) < ε affects F = RS factorization
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 12 / 36
![Page 19: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/19.jpg)
Element Degeneration
Collapse with |det(F)| < ε or inversion with det(F) < 0UnphysicalUnavoidable with (finite) linear forcesUnavoidable due to discretizationUnavoidable due to user interactiondet(F) < ε affects F = RS factorization
From [2]
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 12 / 36
![Page 20: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/20.jpg)
Element Degeneration
Collapse with |det(F)| < ε or inversion with det(F) < 0UnphysicalUnavoidable with (finite) linear forcesUnavoidable due to discretizationUnavoidable due to user interactiondet(F) < ε affects F = RS factorization
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 12 / 36
![Page 21: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/21.jpg)
Element Degeneration
Collapse with |det(F)| < ε or inversion with det(F) < 0UnphysicalUnavoidable with (finite) linear forcesUnavoidable due to discretizationUnavoidable due to user interactiondet(F) < ε affects F = RS factorization
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 12 / 36
![Page 22: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/22.jpg)
Outline
1 Introduction
2 Corotational FEM
3 Rotation extraction
4 Degeneration-Aware Polar Decomposition
5 Results
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 13 / 36
![Page 23: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/23.jpg)
Methods
Several methods to extract R from F are possible:Polar Decomposition [1]QR Factorization [3]Hybrid PD-QR [5]Modified Singular Value Decomposition (SVD1) [2]Coherent Singular Value Decomposition (SVD2) [4]Project/Reflect [6]Degeneration-Aware Polar Decomposition [?]
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 14 / 36
![Page 24: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/24.jpg)
Polar Decomposition
Factorizes F = RS, where R is orthonormal and S is symmetricBest matching, minimizes ‖F −R‖2FFails if |det(F)| ≤ ε (collapsed)Reflected R with det(R) = −1 if det(F) < 0 (inverted)
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 15 / 36
![Page 25: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/25.jpg)
Polar Decomposition
Factorizes F = RS, where R is orthonormal and S is symmetricBest matching, minimizes ‖F −R‖2FFails if |det(F)| ≤ ε (collapsed)Reflected R with det(R) = −1 if det(F) < 0 (inverted)
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 15 / 36
![Page 26: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/26.jpg)
Polar Decomposition
Factorizes F = RS, where R is orthonormal and S is symmetricBest matching, minimizes ‖F −R‖2FFails if |det(F)| ≤ ε (collapsed)Reflected R with det(R) = −1 if det(F) < 0 (inverted)
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 15 / 36
![Page 27: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/27.jpg)
QR Factorization
Factorizes F = RE using Gram-Schmidt orthonormalization, where Ris orthonormal and E is upper-triangular
Fast and RobustHandles collapsed and inverted elements seamlesslyInduces AnisotropyCritical point on collapse plane
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 16 / 36
![Page 28: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/28.jpg)
QR Factorization
Factorizes F = RE using Gram-Schmidt orthonormalization, where Ris orthonormal and E is upper-triangular
Fast and RobustHandles collapsed and inverted elements seamlesslyInduces AnisotropyCritical point on collapse plane
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 16 / 36
![Page 29: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/29.jpg)
QR Factorization
Factorizes F = RE using Gram-Schmidt orthonormalization, where Ris orthonormal and E is upper-triangular
Fast and RobustHandles collapsed and inverted elements seamlesslyInduces AnisotropyCritical point on collapse plane
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 16 / 36
![Page 30: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/30.jpg)
Hybrid PD-QR
Use Polar Decomposition for undegenerate elements and QR fordegenerate ones below a threshold det(F) < α
Inherits good properties of PD and QR......but also the drawbacks of QR......and adds discontinuity across transition
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36
![Page 31: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/31.jpg)
Hybrid PD-QR
Use Polar Decomposition for undegenerate elements and QR fordegenerate ones below a threshold det(F) < α
Inherits good properties of PD and QR......but also the drawbacks of QR......and adds discontinuity across transition
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36
![Page 32: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/32.jpg)
Hybrid PD-QR
Use Polar Decomposition for undegenerate elements and QR fordegenerate ones below a threshold det(F) < α
Inherits good properties of PD and QR......but also the drawbacks of QR......and adds discontinuity across transition
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36
![Page 33: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/33.jpg)
Modified Singular Value Decomposition (SVD1)
SVD factorizes F = UFVT , where U and VT are orthonormal and F isdiagonal, with singular values σi ≥ σi+1 ≥ 0
Handles collapsed elementsHandles inverted elements negating the smallest σi
R = U VT equivalent to “invertible” PDComputationally expensiveCritical point at repeated singular values
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 18 / 36
![Page 34: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/34.jpg)
Modified Singular Value Decomposition (SVD1)
SVD factorizes F = UFVT , where U and VT are orthonormal and F isdiagonal, with singular values σi ≥ σi+1 ≥ 0
Handles collapsed elementsHandles inverted elements negating the smallest σi
R = U VT equivalent to “invertible” PDComputationally expensiveCritical point at repeated singular values
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 18 / 36
![Page 35: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/35.jpg)
Modified Singular Value Decomposition (SVD1)
SVD factorizes F = UFVT , where U and VT are orthonormal and F isdiagonal, with singular values σi ≥ σi+1 ≥ 0
Handles collapsed elementsHandles inverted elements negating the smallest σi
R = U VT equivalent to “invertible” PDComputationally expensiveCritical point at repeated singular values
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 18 / 36
![Page 36: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/36.jpg)
Coherent Singular Value Decomposition (SVD2)
Tries to improve SVD1 by enforcing temporal coherence in the inferredinversion direction (smallest σi )
Better aligned forcesComputationally expensiveDiscontinuous
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 19 / 36
![Page 37: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/37.jpg)
Coherent Singular Value Decomposition (SVD2)
Tries to improve SVD1 by enforcing temporal coherence in the inferredinversion direction (smallest σi )
Better aligned forcesComputationally expensiveDiscontinuous
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 19 / 36
![Page 38: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/38.jpg)
Coherent Singular Value Decomposition (SVD2)
Tries to improve SVD1 by enforcing temporal coherence in the inferredinversion direction (smallest σi )
Better aligned forcesComputationally expensiveDiscontinuous
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 19 / 36
![Page 39: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/39.jpg)
Comparison: Force Fields
Identity and PD QR
SVD1 and SVD2 Project and Reflect
Forcefields experienced by x1 (red) by in a [−3,3]2 region when x2,x3(green, blue) are fixed
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 20 / 36
![Page 40: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/40.jpg)
Comparison: Trajectories
DAPDQR MSVD
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 21 / 36
![Page 41: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/41.jpg)
Discussion
The main problems are:Discrete-time heuristic ignores trajectories: QR, SVD1, SVD2Critical points induce counter-intuitive rotations: QR, SVD1Discontinuity induces jitter: Hybrid PD-QR, SVD2
Our solution: Degeneration-Aware Polar Decomposition
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 22 / 36
![Page 42: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/42.jpg)
Discussion
The main problems are:Discrete-time heuristic ignores trajectories: QR, SVD1, SVD2Critical points induce counter-intuitive rotations: QR, SVD1Discontinuity induces jitter: Hybrid PD-QR, SVD2
Our solution: Degeneration-Aware Polar Decomposition
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 22 / 36
![Page 43: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/43.jpg)
Outline
1 Introduction
2 Corotational FEM
3 Rotation extraction
4 Degeneration-Aware Polar Decomposition
5 Results
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 23 / 36
![Page 44: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/44.jpg)
Idea
Use continuous time to detect collapseCompute R avoiding rotation-past-collapse
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 24 / 36
![Page 45: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/45.jpg)
Detecting collapse (I)
Given a timestep [t0, t1]:1 New degeneration if det(F(t0)) > α and det(F(t1)) ≤ α2 Compute tc ∈ [t0, t1] solving det(F(tc)) = α
3 Compute the collapse feature pair (Ac ,Bc) at tcDetails:
In 2D only V-E case, in 3D V-F and E-E casesDegeneration threshold α > 0 improves numerical robustness
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 25 / 36
![Page 46: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/46.jpg)
Detecting collapse (I)
Given a timestep [t0, t1]:1 New degeneration if det(F(t0)) > α and det(F(t1)) ≤ α2 Compute tc ∈ [t0, t1] solving det(F(tc)) = α
3 Compute the collapse feature pair (Ac ,Bc) at tcDetails:
In 2D only V-E case, in 3D V-F and E-E casesDegeneration threshold α > 0 improves numerical robustness
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 25 / 36
![Page 47: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/47.jpg)
Detecting collapse (II)
x1
x2 x3
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 26 / 36
![Page 48: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/48.jpg)
Computing R (I)
We compute R for degenerate elements as follows:1 Compute degeneration direction dc from (Ac ,Bc)
2 Displace nodes of feature Ac as x i = x i + λ(α, β)dc
3 Compute F from displaced nodes x i
4 Compute R using Polar Decomposition on FDetails:
λ(α, β) computes a displacement length that reverts elementdegeneration and guarantees det(F) > ε
Parameter β controls force alignment
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 27 / 36
![Page 49: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/49.jpg)
Computing R (I)
We compute R for degenerate elements as follows:1 Compute degeneration direction dc from (Ac ,Bc)
2 Displace nodes of feature Ac as x i = x i + λ(α, β)dc
3 Compute F from displaced nodes x i
4 Compute R using Polar Decomposition on FDetails:
λ(α, β) computes a displacement length that reverts elementdegeneration and guarantees det(F) > ε
Parameter β controls force alignment
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 27 / 36
![Page 50: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/50.jpg)
Computing R (II)
∆h
dc
Reflect dc
Project dc
xc
Project (β = 1) and Reflect (β = 2). Height ∆h depends on α
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 28 / 36
![Page 51: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/51.jpg)
Computing R (III)
The effect of parameter β for values 1, 2 and 10
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 29 / 36
![Page 52: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/52.jpg)
Computing R (VI)
dc
Reflect dc
x (b)c (t)
det(F) = 0
x (a)c (t)
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 30 / 36
![Page 53: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/53.jpg)
Discussion
DAPD:Cached pair (Ac ,Bc) ensures temporal coherencyCoherently aligned recovery forcesContinuous and differentiable everywhere......except when Ac or Bc collapse during inversion
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 31 / 36
![Page 54: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/54.jpg)
Outline
1 Introduction
2 Corotational FEM
3 Rotation extraction
4 Degeneration-Aware Polar Decomposition
5 Results
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 32 / 36
![Page 55: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/55.jpg)
Results
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 33 / 36
![Page 56: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/56.jpg)
Results
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 33 / 36
![Page 57: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/57.jpg)
Results
(a) QR, (b) SVD1, (c) DAPD β = 1, (d) DAPD β = 5, (e) Nonlinear
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 33 / 36
![Page 58: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/58.jpg)
Results
http://www.lsi.upc.edu/˜ocivit/videos/DCFEM_Full.ogghttp://www.lsi.upc.edu/˜ocivit/videos/DCNLFEM-Draft.avi
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 33 / 36
![Page 59: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/59.jpg)
Conclusions
Advantages of DAPD:Increased realismShorter recovery timeFaster than SVD1
Limitations:Cannot handle initially degenerate elementsLimited to triangular/tetrahedral elementsSlower than QR
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 34 / 36
![Page 60: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/60.jpg)
References
M.Muller, M.GrossInteractive virtual materialsProceedings of Graphics Interface 2004
G.Irving, J.Teran, R.FedkiwInvertible finite elements for robust simulation of large deformation,Symposium on Computer Animation 2004
M.Nesme, Y.Payan, F.FaureEfficient, Physically Plausible Finite ElementsEurographics 2005
R.Schmedding, M.TeschnerInversion handling for stable deformable modelingThe Visual Computer 2008
E.Parker, J.O’BrienReal-time deformation and fracture in a game environmentSymposium on Computer Animation 2009
O. Civit-Flores, A. SusınRobust treatment of degenerate elements in interactive corotational FEM simulations(to appear) Computer Graphics Forum 2014
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 35 / 36
![Page 61: Robust Treatment of Degenerate Elements in Interactive ...ocivit/slides/SlidesDAPD.pdf · O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 17 / 36 Modified Singular Value Decomposition](https://reader030.vdocument.in/reader030/viewer/2022040216/5f36c1fc75658820a140d340/html5/thumbnails/61.jpg)
Q&A
?
O. Civit-Flores & A. Susin (UPC) DAPD June 11, 2014 36 / 36