gibbonCode / GIBBON

The Geometry and Image-Based Bioengineering add-On for MATLAB
http://gibboncode.org/
GNU Affero General Public License v3.0
191 stars 60 forks source link

Apparent surface inversion in contact problem #93

Closed TGadzella closed 3 years ago

TGadzella commented 3 years ago

Using the MATLAB PDE toolbox to do the meshing, but otherwise closely following the contact demos to define contact surfaces in a multibody simulation of a highly-simplified periodontal complex (represented in a quarter section). Simulation runs in FEBio 3 seem to show inversion of some of the contact surfaces which cause negative Jacobians. Tried flipping element orders but the surface normals appear to be correct. I'm sure I may have made a mistake somewhere but cannot for the life of me figure out what it is. Thanks!

Project files

InvertedSurface Bodies Body1_Normals Body2_Normals

TJG_PDL_Builder.zip

TGadzella commented 3 years ago

This was my mistake - there was a typo in line 247-248 that was calling the wrong surface!