feat: correct solid mechanics bc#3986
Merged
castelletto1 merged 8 commits intodevelopfrom Feb 27, 2026
Merged
Conversation
rrsettgast
approved these changes
Feb 27, 2026
src/coreComponents/physicsSolvers/solidMechanics/SolidMechanicsLagrangianFEM.cpp
Outdated
Show resolved
Hide resolved
src/coreComponents/physicsSolvers/solidMechanics/SolidMechanicsLagrangianFEM.cpp
Outdated
Show resolved
Hide resolved
src/coreComponents/fieldSpecification/TractionBoundaryCondition.cpp
Outdated
Show resolved
Hide resolved
castelletto1
approved these changes
Feb 27, 2026
jafranc
reviewed
Feb 27, 2026
src/coreComponents/physicsSolvers/solidMechanics/SolidMechanicsLagrangianFEM.cpp
Outdated
Show resolved
Hide resolved
…sLagrangianFEM.cpp Co-authored-by: Nicola Castelletto <38361926+castelletto1@users.noreply.github.com>
…sLagrangianFEM.cpp Co-authored-by: Nicola Castelletto <38361926+castelletto1@users.noreply.github.com>
…n.cpp Co-authored-by: Nicola Castelletto <38361926+castelletto1@users.noreply.github.com>
Contributor
Author
|
The following are the justifications for a update the baseline ID:
|
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 PR introduces exact integration of traction terms on 2D boundaries with highly distorted cells.
The setup employed for verification defines a unit cube domain,$\Omega = [0,1]^3$ , subjected to the following boundary conditions:
The “negative” faces ($x=0$ , $y=0$ , $z=0$ ) are assigned zero normal displacement, $u \cdot n = 0$ .
The “positive” faces ($x=1$ , $y=1$ , $z=1$ ) are subjected to a uniform normal compressive stress, $\sigma \cdot n = -10 \text{MPa}$
The resulting configuration is in a state of pure hydrostatic stress throughout the body (no shear stress components).
Before (note the nonplanar displacement contours and the unphysical shear components):

After (note the planar displacement contours and the numerically zero shear components):

The ideas of this implementation and its coverage are provided in testFEMConsistency.cpp from PR-3957.
Credits to @castelletto1 for providing the test mesh.