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

[GeoMechanicsApplication] Modified pressure points location for different order element. #11558

Open
wants to merge 1 commit into
base: master
Choose a base branch
from

Conversation

mnabideltares
Copy link
Contributor

Done for 10 and 15 noded geometries to match the 6 and 10 noded reduced element for pressure.

📝 Description
The pressure for 15 noded element is solved on 10 noded element because for different order elements the pressure must be solved on an element with one order lower.
for 6 noded element, the pressure is solved on 3 noded, which is geometrically compitable, because the 3 noded elemet is a subset of 6 noded element.
However, it is not the case for 15 and 10 noded elements. The location of the points on 10 noded element does not match the 15 noded ones.
For 10 and 6 noded elements, the same issue exist as 6 noded element is not a subset of 10 node element.
Here an interpolation precedure is added to generate new set of points

🆕 Changelog

  • Use interpolation for the location of points for pressure.

Done for 10 and 15 noded geometries to match the 6 and 10 noded reduced element for pressure.
@mnabideltares mnabideltares self-assigned this Sep 13, 2023
@mnabideltares mnabideltares added GeoMechanics Issues related to the GeoMechanicsApplication Bug labels Sep 13, 2023
@roigcarlo
Copy link
Member

Seems this Pr is also affected by #11557. Merge it once its in the master pls

@mnabideltares
Copy link
Contributor Author

Modified pressure points with interpolations for pressure, as they are not at the same locations for 10 and 15 nodes elements.
However, we still see incorrect visualization in Plaxis interface.
Then we found the visualization is well correct when we used 15 nodes elements (with proper numbering for integration points), but not as different order element.
Whether the Plaxis interface visualization for different order element, needs to be investigated.
It is more probably that the interpolation of pressure in Kratos (for difforder) needs to be adjusted, because the pressure finally gets interpolated back to the 15 noded element,

@mnabideltares mnabideltares changed the title Modified pressure points location for different order element. [GeoMechanicsApplication] Modified pressure points location for different order element. Feb 28, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Bug Bugfix GeoMechanics Issues related to the GeoMechanicsApplication
Projects
Development

Successfully merging this pull request may close these issues.

[GeoMechanicsApplication] 15 Noded element leads to irregular stresses
2 participants