Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Unphysical oscillations in tangential forces for contact problem #24632

Closed
recuero opened this issue Jun 8, 2023 Discussed in #24562 · 0 comments · Fixed by #24643
Closed

Unphysical oscillations in tangential forces for contact problem #24632

recuero opened this issue Jun 8, 2023 Discussed in #24562 · 0 comments · Fixed by #24643
Assignees

Comments

@recuero
Copy link
Contributor

recuero commented Jun 8, 2023

Look into formulation for penalty node-face contact, particularly the definition of tangential traction from the normal pressure and its Jacobian. Look at one benchmark problem.

See discussion below.

Discussed in #24562

@recuero recuero self-assigned this Jun 8, 2023
recuero added a commit to recuero/moose that referenced this issue Jun 9, 2023
…b#24632)

This commit changes the way the normal contact force is used in the
node-face penalty friction algorithm. The test that failed with this
change only tests a point contact for a couple of time steps. That
failure was related to a too-large penalty factor. With the previous
approach, a too-large penalty factor may not initially have a negative
influence because the algorithm was using a cached normal force value.
Slightly decreasing that value, however, makes a more extended approach
to penalty contact work for the single point test modified in this commit.

Another example of a frictional cylinder is added to avoid regressions.
That test checks for tangential slips in opposite directions. The previous
formulation is unable to run this example.
recuero added a commit to recuero/moose that referenced this issue Jun 9, 2023
…b#24632)

This commit changes the way the normal contact force is used in the
node-face penalty friction algorithm. The test that failed with this
change only tests a point contact for a couple of time steps. That
failure was related to a too-large penalty factor. With the previous
approach, a too-large penalty factor may not initially have a negative
influence because the algorithm was using a cached normal force value.
Slightly decreasing that value, however, makes a more extended approach
to penalty contact work for the single point test modified in this commit.

Another example of a frictional cylinder is added to avoid regressions.
That test checks for tangential slips in opposite directions. The previous
formulation is unable to run this example.
recuero added a commit to recuero/moose that referenced this issue Jun 13, 2023
…b#24632)

This commit changes the way the normal contact force is used in the
node-face penalty friction algorithm. The test that failed with this
change only tests a point contact for a couple of time steps. That
failure was related to a too-large penalty factor. With the previous
approach, a too-large penalty factor may not initially have a negative
influence because the algorithm was using a cached normal force value.
Slightly decreasing that value, however, makes a more extended approach
to penalty contact work for the single point test modified in this commit.

Another example of a frictional cylinder is added to avoid regressions.
That test checks for tangential slips in opposite directions. The previous
formulation is unable to run this example.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

Successfully merging a pull request may close this issue.

1 participant