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

Holzapfel Ogden Modified Anisotropy (HO-ma) Constitutive Model #272

Merged
merged 17 commits into from
Oct 16, 2024

Conversation

divyaadil23
Copy link
Contributor

@divyaadil23 divyaadil23 commented Sep 13, 2024

Same as the previous pull request, but the commits have been squashed into one.

Current situation

#249

Release Notes

  • Added the Holzapfel Ogden Modified Anisotropy (HO-ma) Constitutive Model to material models.
  • Updated Holzapfel Ogden Model.
  • Both HO and HO-ma models are taken from https://github.com/vvedula22/svFSI/blob/master/Code/Source/svFSI/MATMODELS.f and modified. Added additional terms containing the second derivative of smoothed heaviside funciton to the elasticity tensor to make it exact. Changed the approximate derivatives of heaviside function to exact. Removed fds (dot product of fiber directions in 8th invariant).
  • Tested it with the new framework of unit tests - the unit tests now test the order of convergence for finite difference estimates.
  • Added some integrated tests for struct and ustruct for both models.

Documentation

Where can I document how to use the new material model? Also, should I add the derivation?

Testing

Added 4 integration tests for HO and HO-ma models. Added new unit tests for both.

Code of Conduct & Contributing Guidelines

…del to material models.

Updated Holzapfel Ogden Model.
Both HO and HO-ma models are taken from https://github.com/vvedula22/svFSI/blob/master/Code/Source/svFSI/MATMODELS.f and modified. Added additional terms containing the second derivative of smoothed heaviside funciton to the elasticity tensor to make it exact. Changed the approximate derivatives of heaviside function to exact. Removed fds (dot product of fiber directions in 8th invariant).
Tested it with the new framework of unit tests - the unit tests now test the order of convergence for finite difference estimates.
Added some integrated tests for struct and ustruct for both models.
@aabrown100-git
Copy link
Collaborator

Looks like in your most recent commit (merging SimVascular::main), your new tests failed. I think you might just have to rename your new test input files from svFSI.xml to svFSIplus.xml.

Code/Source/svFSI/read_files.cpp Outdated Show resolved Hide resolved
Code/Source/svFSI/set_material_props.h Outdated Show resolved Hide resolved
tests/cases/fluid/pipe_RCR_3d_trilinos_bj/README.md Outdated Show resolved Hide resolved
Copy link
Collaborator

@aabrown100-git aabrown100-git left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

All looks good to me!

…del to material models.

Updated Holzapfel Ogden Model.
Both HO and HO-ma models are taken from https://github.com/vvedula22/svFSI/blob/master/Code/Source/svFSI/MATMODELS.f and modified. Added additional terms containing the second derivative of smoothed heaviside funciton to the elasticity tensor to make it exact. Changed the approximate derivatives of heaviside function to exact. Removed fds (dot product of fiber directions in 8th invariant).
Tested it with the new framework of unit tests - the unit tests now test the order of convergence for finite difference estimates.
Added some integrated tests for struct and ustruct for both models.
@aabrown100-git
Copy link
Collaborator

Looks like your HO and HO-ma integration tests failed. It's probably because the reference solutions were generated with the default solid_viscosity value, which was erroneously set to 0.9 previously. The most recent commit changes the default value to 0.0, so the integration tests now use a solid_viscosity value of 0.0.

You can verify this by setting the solid_viscosity parameter to 0.9 in your test input files
<Solid_viscosity> 0.9 </Solid_viscosity>
The tests should pass. If so, I think you can just rerun your tests with 0 solid viscosity (just delete it from the input file) and update your reference solution with the result.

@ktbolt ktbolt removed the request for review from msbazzi October 3, 2024 17:45
Copy link

codecov bot commented Oct 5, 2024

Codecov Report

Attention: Patch coverage is 85.24845% with 95 lines in your changes missing coverage. Please review.

Project coverage is 66.38%. Comparing base (82c9ea6) to head (9f79d78).

Files with missing lines Patch % Lines
tests/unitTests/test.h 67.49% 92 Missing ⚠️
Code/Source/svFSI/Parameters.cpp 50.00% 2 Missing ⚠️
Code/Source/svFSI/read_files.cpp 0.00% 1 Missing ⚠️
Additional details and impacted files
@@            Coverage Diff             @@
##             main     #272      +/-   ##
==========================================
+ Coverage   65.67%   66.38%   +0.71%     
==========================================
  Files         116      116              
  Lines       29574    30334     +760     
==========================================
+ Hits        19422    20138     +716     
- Misses      10152    10196      +44     

☔ View full report in Codecov by Sentry.
📢 Have feedback on the report? Share it here.

README.md Outdated Show resolved Hide resolved
@divyaadil23
Copy link
Contributor Author

@ktbolt we have resolved the comments on this pull request. Could you review it again? Thanks!

@ktbolt ktbolt merged commit f9fe6fe into SimVascular:main Oct 16, 2024
5 checks passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

6 participants