Hello oh wise and beneficent FEBio developers,
I have successfully implemented the prestrain plugin on a simple cylindrical geometry, but when I apply the same formulation to a complex geometry (tet4), I get the following error message:
************************************************** ***********************
* ERROR *
* *
* NAN Detected. *
************************************************** ***********************
The formulation I am using sets the local axes using "mat_axis" like this:
<ElementData>
<element id="1">
<mat_axis>
<a>-0.96593,0.25882,1.6273e-19</a>
<d>6.6617e-16,-1.6663e-16,1</d>
</mat_axis>
</element>
And the in situ prestrain is defined thus:
<material id="1" name="Material1" type="prestrain elastic">
<density>1</density>
<prestrain type="in-situ stretch">
<stretch>3</stretch>
<isochoric>1</isochoric>
</prestrain>
<elastic type="neo-Hookean">
<density>1</density>
<E>1</E>
<v>0.3</v>
</elastic>
</material>
Why am I getting NAN only for the complex geometry and not the simple one? I have attached the simple and complex geometries. Oops, the complex is 2.4 MB and can't be attached to the forum. -David
I have successfully implemented the prestrain plugin on a simple cylindrical geometry, but when I apply the same formulation to a complex geometry (tet4), I get the following error message:
************************************************** ***********************
* ERROR *
* *
* NAN Detected. *
************************************************** ***********************
The formulation I am using sets the local axes using "mat_axis" like this:
<ElementData>
<element id="1">
<mat_axis>
<a>-0.96593,0.25882,1.6273e-19</a>
<d>6.6617e-16,-1.6663e-16,1</d>
</mat_axis>
</element>
And the in situ prestrain is defined thus:
<material id="1" name="Material1" type="prestrain elastic">
<density>1</density>
<prestrain type="in-situ stretch">
<stretch>3</stretch>
<isochoric>1</isochoric>
</prestrain>
<elastic type="neo-Hookean">
<density>1</density>
<E>1</E>
<v>0.3</v>
</elastic>
</material>
Why am I getting NAN only for the complex geometry and not the simple one? I have attached the simple and complex geometries. Oops, the complex is 2.4 MB and can't be attached to the forum. -David
Comment