Prestrain on local axes

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • birdave
    Junior Member
    • Dec 2014
    • 6

    Prestrain on local axes

    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
    Attached Files
  • maas
    Lead Code Developer
    • Nov 2007
    • 3400

    #2
    Hi David,

    In your file "Prestrain_lig_PS.feb" (which Dr. Ateshian forwarded to me), on line 11879 the value for "d" is NAN. This is the mat_axis definition for element 22:

    Code:
    			<element id="22">
    				<mat_axis>
    					<a>-0.96534,0.15514,-0.20988</a>
    					<d>nan,nan,nan</d>
    				</mat_axis>
    			</element>
    In fact, it looks like there are quite a few nans in the input file. How did you generate this data? Was this PreView, FEBioStudio, or a custom code?

    Steve
    Department of Bioengineering, University of Utah
    Scientific Computing and Imaging institute, University of Utah

    Comment

    • birdave
      Junior Member
      • Dec 2014
      • 6

      #3
      Steve,

      Hello and thanks for the speedy reply. The local axes were generated with a custom code. Somewhere I divided by zero, and failed to send the code without the nans. I have some code without the nans, but it still generates an "error nan detected". Where could I send it?

      David

      Comment

      • maas
        Lead Code Developer
        • Nov 2007
        • 3400

        #4
        Hi David,

        You can send it to steve dot maas at Utah dot edu.

        Cheers,

        Steve
        Department of Bioengineering, University of Utah
        Scientific Computing and Imaging institute, University of Utah

        Comment

        • bracamontej
          Junior Member
          • Mar 2019
          • 14

          #5
          Hello! I've been trying to achieve prestrain convergence on cylindrical geometries. I got really excited reading that you were able to get prestrain to converge on an ideal cylinder.

          Do you have any advice on how to get convergence? Which boundary conditions did you apply? Was the prescribed prestrain field such that the material remained isochoric (incompressible)?

          Any hint or advice is welcomed!

          Thanks!

          Comment

          Working...
          X
          😀
          😂
          🥰
          😘
          🤢
          😎
          😞
          😡
          👍
          👎