feat: Mortar Method #3879
Open
feat: Mortar Method #3879
Conversation
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
danielemoretto44
added
type: feature
danielemoretto44
added
the
flag: requires rebaseline
danielemoretto44
requested review from
Guotong-Ren,
OmarDuran,
andrea-borio,
castelletto1,
cssherman,
dkachuma,
jhuang2601 and
matteofrigo5
as code owners
Closed
paveltomin
reviewed
Oct 23, 2025
paveltomin
reviewed
Oct 23, 2025
inputFiles/mortar/VTKplot.pvd
Outdated
Collaborator
There was a problem hiding this comment.
results should not be here
Contributor
|
@danielemoretto44 Could you please provide the reference of the implemented method, if available? A documentation for the work is highly expected. |
paveltomin
reviewed
Oct 27, 2025
| protected: | ||
|
|
||
| /// The array containing the nodal position array. | ||
| //arrayView2d< real64 const, nodes::REFERENCE_POSITION_USD > const m_X; |
paveltomin
reviewed
Oct 27, 2025
| DofManager::Connector::Elem, | ||
| meshTargets ); | ||
|
|
||
| //dofManager.addCoupling( solidMechanics::totalDisplacement::key(), |
paveltomin
reviewed
Oct 27, 2025
| solution.create( dofManager.numLocalDofs(), MPI_COMM_GEOS ); | ||
|
|
||
|
|
||
| //Write the matrix pattern to a file for debugging |
paveltomin
reviewed
Oct 27, 2025
|
|
||
| // ParallelMatrix parallel_matrix_3; | ||
| // parallel_matrix_3.create( localMatrix.toViewConst(), dofManager.numLocalDofs(), MPI_COMM_GEOS ); | ||
| // parallel_matrix_3.write("systemMatrix.mtx"); |
Collaborator
There was a problem hiding this comment.
clean up this function
paveltomin
reviewed
Oct 27, 2025
| /* | ||
| * ------------------------------------------------------------------------------------------------------------ | ||
| * SPDX-License-Identifier: LGPL-2 std::map< ElementShape, array1d< localIndex > > const & | ||
| faceTypesToFaceElementsSlave = m_faceTypeToElementList.at( MortarSide::Slave ); |
Collaborator
There was a problem hiding this comment.
something happened here
Author
@jhuang2601 I updated the description of the PR with some references. Please let me know if more material is needed. |
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 implements the mortar method for solid mechanics to enforce continuity between the displacements of non-conforming mesh bodies.
The formulation uses piecewise constant Lagrange multipliers with one-sided bubble stabilization, as described in
Hauret, P., & Le Tallec, P. (2007). A discontinuous stabilized mortar method for general 3d elastic problems. Computer Methods in Applied Mechanics and Engineering, 196(49-52), 4881-4900.
Useful references for implementation details:
contact searching:
Yang, B., & Laursen, T. A. (2008). A contact searching algorithm including bounding volume trees applied to finite sliding mortar formulations. Computational mechanics, 41(2), 189-205.
segment based quadrature scheme:
Laursen, T. A., Puso, M. A., & Sanders, J. (2012). Mortar contact formulations for deformable–deformable contact: past contributions and new extensions for enriched and embedded interface formulations. Computer methods in applied mechanics and engineering, 205, 3-15.