Hello Osama, I currently do this using Abaqus. Im not sure how it is in ANSYS but I dont think it will be much different. Bone remodels in order to maintain homeostasis. This can be defined as a certain target strain or stress state. If the stimulus is above that target, bone resorbs, if it is above, bone apposition occurs. For more information on this, please refer to Frost HM. Bone ‘mass’ and the ‘mechanostat’: a proposal. Anatomical Record 1987; 219: 1–9. DOI: 10.1002/ar.1092190104.
I have implemented a bone adaptation algorithm in Abaqus. What I do is I load the femur under physiological boundary conditions, extract the stimulus I want from the results of the model (in this case, strain). Then, using a sub-routine implemented in MATLAB, I update the material properties, making bone less or more stiff, less or more dense, as required in order to reach that stimulus. I re-load the model, and continue with the algorithm iteratively until convergence is achieved. I am sure you can do something similar with ANSYS. If you are struggling with it, I can help you implementing it in Abaqus.
For more information on my remodelling system, please refer to: Geraldes D. M., and Phillips A. T. M. (2014), A comparative study of orthotropic and isotropic bone adaptation in the femur, International Journal for Numerical Methods in Biomedical Engineering, 30, pages 873–889, doi: 10.1002/cnm.2633
Dear Dr. Diogo Geraldes,
Thank you very much for your detailed answer. your information and paper are very interesting and useful to me. I am currently trying to implement the approach with some modifications in material assignment using ANSYS. I intend to assign an anisotropic material for the femur based on the ct data.
Best regards
- Badges
- Science topic