I'm using the PreStrain Plugin in FEBio 2.8.4
My model is a disk held by the edge and pressurized. If I try to run it without the prestrain plugin it converges fine. However, if I try to use it with the prestrain material, even if I set the prestrain to the identity, the problem does not converge. I'm wandering if I'm doing anything wrong. Attached are the two feb files for demonstrating this behavior:
Without prestrain: minidisk_tilt_Equation_Membrane_In_NOPreStrain.feb
With prestrain (but F0=Identity): minidisk_tilt_Equation_Membrane_In_PreStrain.feb
The two feb files are identical except for the following differences:
No Prestrain:
Prestrain:
My model is a disk held by the edge and pressurized. If I try to run it without the prestrain plugin it converges fine. However, if I try to use it with the prestrain material, even if I set the prestrain to the identity, the problem does not converge. I'm wandering if I'm doing anything wrong. Attached are the two feb files for demonstrating this behavior:
Without prestrain: minidisk_tilt_Equation_Membrane_In_NOPreStrain.feb
With prestrain (but F0=Identity): minidisk_tilt_Equation_Membrane_In_PreStrain.feb
The two feb files are identical except for the following differences:
No Prestrain:
Code:
<Material> <material id="1" name="Material1" type="uncoupled solid mixture"> <solid type="Mooney-Rivlin"> ... </solid> <solid type="fiber-exp-pow-uncoupled"> ... </solid> </material> </Material>
Prestrain:
Code:
<Material> <material id="1" name="Material1" type="uncoupled prestrain elastic"> <prestrain type="prestrain gradient"/> <elastic type="uncoupled solid mixture"> <solid type="Mooney-Rivlin"> ... </solid> <solid type="fiber-exp-pow-uncoupled"> ... </solid> </elastic> </material> </Material> <MeshData> <ElementData elem_set="Part1" var="F0"> <elem lid="1">1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0</elem> ... <elem lid="149">1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0</elem> </ElementData> </MeshData> <Constraints> <constraint type="prestrain"> <tolerance>0.1</tolerance> <min_iters>3</min_iters> </constraint> </Constraints>
Comment