[SolidMechanics] Add component for source term integration based on new FEM framework#6196
Open
th-skam wants to merge 2 commits into
Open
[SolidMechanics] Add component for source term integration based on new FEM framework#6196th-skam wants to merge 2 commits into
th-skam wants to merge 2 commits into
Conversation
Contributor
|
I am not a FEM expert so I don't feel confident to review this PR. Nevertheless I checked to make sure this is compatible with higher order element and it seems ton be the case. So it is fine by me given that @alxbilger and/or @hugtalbot review it. Could you add an example ? |
This file contains hidden or 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.
This is to allow integrating a source term (load) on the rhs of the PDE. Let's say we have
So$f$ is force per unit volume (provided by the user) and is sampled at the nodes. We integrate it with a quadrature rule to compute nodal forces $F_i$ .
I modelled the component after
FEMMassthough I only implementedaddForcefor now. For my current use-case, the source terms are set on the reference configuration and are only explicit contributions. It's been tested through the verification process with MMS.Discussion
There's some similarity with what
FEMMassis doing. It would be worthwhile to discuss a possible code factorization of these components even with this PR.And, on top of that, perhaps we could factorize the computations in the scene: one component handling the matrix assembly and the other one responsible for using it to compute addForce. This way, if there are more than one source terms, we do not visit the elements/quadratures multiple times. That's to be seen.
Fixes #6166