Hello everyone,
I am working on a 3D linear elasticity problem and developing a lower order stress simulation framework from scratch. For this, I’m using the Finite Cell Method (a fictitious domain approach) with axis-parallel hexahedral cells. These cells are refined into tetrahedra or smaller hexahedra, categorized as either fully inside or fully outside the domain boundary.
The meshing is working well, but I am stuck on implementing adaptive integration for cells cut by the domain surface. I’ve been using 2x2x2 Gauss quadrature for regular hexahedra, but I am unsure how to properly account for the geometry of intersected cells. My current approach involves weighting the integration based on the volume fraction inside the domain (using a very small weight for outside points), but the results don’t seem correct.
Has anyone worked with similar setups or implemented adaptive integration for cut cells? I’ve attached my routine for stiffness calculation of surface-intersected cells where all elements are inside the domain. Any insights or code review would be greatly appreciated!
I am working on a 3D linear elasticity problem and developing a lower order stress simulation framework from scratch. For this, I’m using the Finite Cell Method (a fictitious domain approach) with axis-parallel hexahedral cells. These cells are refined into tetrahedra or smaller hexahedra, categorized as either fully inside or fully outside the domain boundary.
The meshing is working well, but I am stuck on implementing adaptive integration for cells cut by the domain surface. I’ve been using 2x2x2 Gauss quadrature for regular hexahedra, but I am unsure how to properly account for the geometry of intersected cells. My current approach involves weighting the integration based on the volume fraction inside the domain (using a very small weight for outside points), but the results don’t seem correct.
Has anyone worked with similar setups or implemented adaptive integration for cut cells? I’ve attached my routine for stiffness calculation of surface-intersected cells where all elements are inside the domain. Any insights or code review would be greatly appreciated!
