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

Discrepancy in Expected Deformation Results with Given Units #290

Open
antonagafonov opened this issue Dec 10, 2023 · 2 comments
Open

Discrepancy in Expected Deformation Results with Given Units #290

antonagafonov opened this issue Dec 10, 2023 · 2 comments
Labels
bug Something isn't working

Comments

@antonagafonov
Copy link

Hello,

I'm encountering an issue with the units in my simulation and could use some assistance in pinpointing the problem. My 3D mesh is scaled in meters, and based on this, I'm using Newtons (N) for forces and Pascals (Pa) for the Young's modulus. The simulation aims to replicate deformation forces on human soft tissue, specifically in the arm region. The expected Young's modulus range for this tissue is between 0.1 to 1 GPa.

However, to achieve deformations within a realistic range (approximately 0.5-3 mm), I find myself having to set the Young's modulus $E$ to an unusually low value of 2 $Pa$, while the total force is set to 1000 $N$. These values seem inconsistent with standard expectations.

Could there be an error in my understanding of the units or in the simulation setup? Any insights or suggestions would be greatly appreciated.

Thank you.

@antonagafonov antonagafonov added the bug Something isn't working label Dec 10, 2023
@zfergus
Copy link
Member

zfergus commented Jan 17, 2024

How are you applying the forces? Are they Neumann boundary conditions? If so, the units of NBCs are not force but force over area (i.e., pressure).

@antonagafonov
Copy link
Author

antonagafonov commented Jan 18, 2024

Hi @zfergus ,
Thanks for your reply.
As I see here, it is Force over Face in [N], which is Force in [N] and not Pressure in [Pa].
https://polyfem.github.io/tutorials/ipc_quick_start_guide/ipc_quick_start_guide/?h=neuman#neumann
Also, if I apply pressure over a triangle surface, should it be a scalar?
I am confused...
Can you please explain?
Appreciate your help,
Thanks

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
None yet
Development

No branches or pull requests

2 participants