[DamApplication] Use reference geometry in SmallDisplacementInterfaceElement - #14686
Open
joaquinirazabal wants to merge 1 commit into
Open
[DamApplication] Use reference geometry in SmallDisplacementInterfaceElement#14686joaquinirazabal wants to merge 1 commit into
joaquinirazabal wants to merge 1 commit into
Conversation
…ement - Introduced methods for calculating current and reference rotation matrices, reference Jacobian determinants, and creating reference geometry in small_displacement_interface_element.hpp. - Updated Python bindings to include the new DamPScheme class. - Cleaned up whitespace in dam_analysis.py. - Added new test cases for linear solver and global tangent consistency in the DamApplication. - Implemented tests for interface tangent consistency, ensuring the tangent matrix is consistent across different loading branches and that the reference geometry remains invariant during opening.
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.
📝 Description
This PR makes the quasi-static formulation of
SmallDisplacementInterfaceElementconsistently use the reference configuration.Previously, the local interface frame and integration measure were computed from the current nodal coordinates. Since DamApplication normally uses
move_mesh_flag = true, these quantities changed during the Newton iterations, while the corresponding geometric derivatives were not included in the tangent matrix. This resulted in an inconsistent assembled tangent, especially noticeable for compressed interfaces with large penalty stiffness.The element originates from the Poromechanics small-strain interface-element family, where this geometry dependence is normally masked by the standard fixed-mesh configuration (
move_mesh_flag = false). In DamApplication, the moved-mesh configuration exposed this implicit dependency.Key changes
jump_local = R0 * Nu * u.Validation
move_mesh_flag = truein both opening and mixed opening/compression configurations.1e-5to1e-9.1e3to1e10.small,nightly, andalltest suites pass.🆕 Changelog
SmallDisplacementInterfaceElementto consistently use reference geometry in its quasi-static formulation.move_mesh_flag.