Skip to content
Snippets Groups Projects
Commit e3115e4d authored by Dmitri Naumov's avatar Dmitri Naumov
Browse files

[T/M] Fix corner case in Hertz contact bisection

parent 9b307f3a
No related branches found
No related tags found
No related merge requests found
...@@ -45,6 +45,10 @@ class HertzContactBC(OpenGeoSys.BoundaryCondition): ...@@ -45,6 +45,10 @@ class HertzContactBC(OpenGeoSys.BoundaryCondition):
# BCs for the current nonlinear solver iteration, i.e., the # BCs for the current nonlinear solver iteration, i.e., the
# applied BCs and the information based on which the BCs are # applied BCs and the information based on which the BCs are
# applied are off by one iteration. # applied are off by one iteration.
# Reset left side to x_with_y_excess value, which in some cases
# might be lower
self._a_range[0] = min(self._a_range[0], self._max_x_with_y_excess)
self._a_range[1] = self._max_x_with_y_excess self._a_range[1] = self._max_x_with_y_excess
assert self._a_range[0] <= self._a_range[1] assert self._a_range[0] <= self._a_range[1]
self._tendency_is_up = False self._tendency_is_up = False
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment