[Mapping] MPM-FEM mortar coupling #7886
Open
+2,018
−23
Conversation
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
Description
This PR introduces MPM-FEM mortar coupling made possible by the great work and help of @tteschemacher and @philbucher 🎉 .
This PR mostly introduces a
structure_mpm_modelerthat provides the necessary input for thecoupling_geoms_mapper. Coupling quadrature points are spread along the FEM edge, according togauss_integration_orderin the parameters, which are them mapped to the MPM background grid to later construct the mortar mappers. A current usage limitation (not astructure_mpm_modelerlimitation but acoupling_geom_mapperlimitation) is that the FEM edge must be the slave of the mortar mapper (however, the FEM can be the CoSim origin).Two tests added - conforming and nonconforming. These tests and some other small mortar tests were moved from the nightsuite to the smallsuite.
This PR relies heavily on other pending quadrature point PRs (#7552 ) - it would be great if those can be expedited too!
An example static cantilever beam (solved through CoSim Gauss Seidel) is below (reference answer = -0.01):