imcs-compsim / MIRCO

A shared-memory parallel BEM code for the contact of rough surfaces
MIT License
3 stars 3 forks source link

Add correction factors for pressure-based Green function #90

Closed RShaw026 closed 10 months ago

RShaw026 commented 10 months ago

Description and Context

In PR #87, The pressure-based Green function was added. In this PR, the correction factors corresponding to the pressure-based Green function are added. These values are calculated by using MIRCO to solve a contact between a linear elastic half-plane and a rigid flat-indenter.

Related Issues and Pull Requests

How Has This Been Tested?

ctest is passing

Checklist

Additional Information

Interested Parties / Possible Reviewers

@isteinbrecher @mayrmt

RShaw026 commented 10 months ago

@isteinbrecher @mayrmt FYI Since MIRCO currently consumes a lot of memory and time to solve high-resolution problems with large contact area percentages (in this case: flat indenters), our MATLAB-based bem-code is used to calculate these factors.

RShaw026 commented 10 months ago

@RShaw026 These factors are constant, so you could make them constant in code as well.

Am I right, that each vector entry corresponds to a certain resolution? Where is it documented, which resolutions these vectors cover?

Yes, you are right. It is not mentioned here, but currently, MIRCO works only from the resolution of 1 to 8. (2^n+1 times 2^n+1 nodes grid)