The Question about Pre-strain

Overall purposes of the project are:

1) To provide free access to a three-dimensional finite element representation of the knee joint

2) Open development of a knee model for collaborative testing and use
User avatar
Ariel Schwartz
Posts: 16
Joined: Tue Feb 20, 2018 7:20 am

Re: The Question about Pre-strain

Post by Ariel Schwartz » Mon Nov 23, 2020 8:03 am

Hi Marloes,

You should check the Outputs section of you febio file, to make sure you are requesting the "fiber stretch" output. It is not automatically added. See the febio users manual for details on requesting outputs.

Best,
Ariel

User avatar
Marloes vdH
Posts: 8
Joined: Fri Nov 06, 2020 1:07 am

Re: The Question about Pre-strain

Post by Marloes vdH » Wed Nov 25, 2020 2:40 am

Hi dr Schwartz,

Thank you for your quick reply. The fiber stretch was included inside the output as <var type="fiber stretch"/>.
I thought this should be the way to request the fiber stretch in the output, but it didn't give any results. It seems like the PreStrain isn't applied to the fibers.

Kind Regards,
Marloes
Attachments
Schermafbeelding 2020-11-20 om 16.38.22.png
Schermafbeelding 2020-11-20 om 16.38.22.png (203.28 KiB) Viewed 1839 times

User avatar
Ahmet Erdemir
Posts: 77
Joined: Sun Sep 10, 2006 1:35 pm

Re: The Question about Pre-strain

Post by Ahmet Erdemir » Mon Dec 14, 2020 6:26 am

Dear Marloes,

Did you try by defining fiber direction outside the elastic material definition? As

<material id="11" name="acl" type="uncoupled prestrain elastic">
<density>1.5e-09</density>
<k>73.2</k>
<elastic type="trans iso Mooney-Rivlin">
<c1>1.95</c1>
<c2>0</c2>
<c3>0.0139</c3>
<c4>116.22</c4>
<c5>535.039</c5>
<lam_max>1.046</lam_max>
</elastic>
<prestrain type="in-situ stretch">
<stretch lc="1">1</stretch>
<isochoric>1</isochoric>
</prestrain>
<fiber type="local">1,4</fiber>
</material>

User avatar
Marloes vdH
Posts: 8
Joined: Fri Nov 06, 2020 1:07 am

Re: The Question about Pre-strain

Post by Marloes vdH » Tue Dec 15, 2020 4:42 am

Dear Ahmet Erdemir,

I did try that, however, in 2.5 files 'fiber' is seen as a synonym for 'mat_axis'. In 3.0 files this is no longer the case. When loading it into FEBio Studio (3.0 files) it shifts into the MRTI part. I also tried to include this axes after loading the 2.5 file with the fiber outside the MRTI however, nothing seems to work... The prestretch seems to be applied on the matrix in stead of the fibers. Furthermore, when I don't fix the femur in all DOF there is a +0.2 translation in z direction (upwards) and the orientation of the fibers changes.

Besides, I had an other small question. Since the assumption of a hyperelastic material is made when using the prestrain, I suppose it will be difficult to go to a viscoelastic material.

In the end I want to simulate several gait cycles and look at the tibial cartilage stresses. Therefore, I want to include the viscoelastic behavior of the ligaments, which is time dependent, hypoelastic. I can imagine that when applying the prestrain on a viscoelastic material there will be problems with the relaxation of the material.

Is there an other possibility to apply prestrain to the material?

Kind regards,
Marloes

User avatar
Ahmet Erdemir
Posts: 77
Joined: Sun Sep 10, 2006 1:35 pm

Re: The Question about Pre-strain

Post by Ahmet Erdemir » Tue Dec 15, 2020 6:49 am


I did try that, however, in 2.5 files 'fiber' is seen as a synonym for 'mat_axis'. In 3.0 files this is no longer the case. When loading it into FEBio Studio (3.0 files) it shifts into the MRTI part.
We are currently using FEBio version 2.9. We also create our models in a scripted fashion, without using PreView nor FEBio Studio. We have seen some compatibility issues importing our files to FEBio Studio, i.e. while they are compliant with FEBio file format 2.5, sometimes things are not fully recognized. I recommend you to start with a simple model. See https://simtk.org/svn/openknee/app/PrestrainExamples/ for one that works with FEBio version 2.9 and PreStrain plug in. You may want to first run this in FEBio version 2.9 and then port it to FEBio version 3 to get accustomed to format changes and make sure that you can visualize fiber stretch and prestrain stretch. You may want to reach out to Steve Maas in regard to nuances of FEBio Studio and the updated file format.


I also tried to include this axes after loading the 2.5 file with the fiber outside the MRTI however, nothing seems to work... The prestretch seems to be applied on the matrix in stead of the fibers.
These material models are essentially solid mixtures. There is a ground substance (matrix) and a fiber component. When you specify in situ stretch it should be applied to fiber. In the example above, it is applied in a weak way: A target overall in situ stretch is specified. However, this is allowed to adjust with the aim of decreasing deviations from geometry due to in situ stretch application. In the end, when in situ stretch is applied, the geometry is deformed (at t=1 of the simulation, fiber stretch (indicating deformation) is different than 1 and prestrain stretch (indicating both applied in situ stretch and deformation)

Furthermore, when I don't fix the femur in all DOF there is a +0.2 translation in z direction (upwards) and the orientation of the fibers changes.
When you apply prestrain, the ligament becomes loaded and all the material points should be in equilibrium. If you want to apply the prestrain exactly and if this is not compatible with the geometry (i.e., equilibrium is not possible at the reference geometry), it will deform the geometry (therefore your femur may move). If you want to keep geometry fixed, you need to apply pre-strain weakly, e.g. allow it to adjust while maintaining equilibrium and minimizing deformations from the reference geometry. The above example tries that. By increasing number of augmentations during in situ application, you try to reduce deviations from the reference geometry.

Besides, I had an other small question. Since the assumption of a hyperelastic material is made when using the prestrain, I suppose it will be difficult to go to a viscoelastic material.
I suspect that this will be possible. We never tried it but the way material properties are defined in FEBio should allow it. I guess this will be an adventure to learn from.
Is there an other possibility to apply prestrain to the material?
I recommend you to read the PreStrain manual and also check https://pubmed.ncbi.nlm.nih.gov/27131609/.

POST REPLY