<?xml version="1.0" encoding="ISO-8859-1"?> <febio_spec version="2.5"> <!--This document was written by C:\oks_simulations\oks007\model\MedToFebio.py on Mar 27 2019 at 11:03:54 AM--> <Module type="solid" /> <Material> <material id="1" name="PTB" type="rigid body"> <density> 1e-9</density> </material> <material id="2" name="FBB" type="rigid body"> <density> 1e-9</density> </material> <material id="3" name="TBB" type="rigid body"> <density> 1e-9</density> </material> <material id="4" name="FMB" type="rigid body"> <density> 1e-9</density> </material> <material id="5" name="QAT" type="Mooney-Rivlin"> <c1> 0.1</c1> <c2> 0</c2> <k> 100</k> <density> 1e-9</density> </material> <material id="6" name="TBC-L" type="Mooney-Rivlin"> <c1> 0.1</c1> <c2> 0</c2> <k> 100</k> <density> 1e-9</density> </material> <material id="7" name="PCL" type="Mooney-Rivlin"> <c1> 0.1</c1> <c2> 0</c2> <k> 100</k> <density> 1e-9</density> </material> <material id="8" name="PTC" type="Mooney-Rivlin"> <c1> 0.1</c1> <c2> 0</c2> <k> 100</k> <density> 1e-9</density> </material> <material id="9" name="ACL" type="Mooney-Rivlin"> <c1> 0.1</c1> <c2> 0</c2> <k> 100</k> <density> 1e-9</density> </material> <material id="10" name="MCL" type="Mooney-Rivlin"> <c1> 0.1</c1> <c2> 0</c2> <k> 100</k> <density> 1e-9</density> </material> <material id="11" name="PTL" type="Mooney-Rivlin"> <c1> 0.1</c1> <c2> 0</c2> <k> 100</k> <density> 1e-9</density> </material> <material id="12" name="MNS-L" type="Mooney-Rivlin"> <c1> 0.1</c1> <c2> 0</c2> <k> 100</k> <density> 1e-9</density> </material> <material id="13" name="MNS-M" type="Mooney-Rivlin"> <c1> 0.1</c1> <c2> 0</c2> <k> 100</k> <density> 1e-9</density> </material> <material id="14" name="LCL" type="Mooney-Rivlin"> <c1> 0.1</c1> <c2> 0</c2> <k> 100</k> <density> 1e-9</density> </material> <material id="15" name="TBC-M" type="Mooney-Rivlin"> <c1> 0.1</c1> <c2> 0</c2> <k> 100</k> <density> 1e-9</density> </material> <material id="16" name="FMC" type="Mooney-Rivlin"> <c1> 0.1</c1> <c2> 0</c2> <k> 100</k> <density> 1e-9</density> </material> </Material> <Geometry from="Geometry.feb" /> <Boundary> <rigid name="QAT_With_PTB" node_set="QAT_@_PTB_TiesNodes" rb="1" /> <rigid name="TBC-L_With_TBB" node_set="TBC-L_@_TBB_TiesNodes" rb="3" /> <rigid name="PCL_With_TBB" node_set="PCL_@_TBB_TiesNodes" rb="3" /> <rigid name="PCL_With_FMB" node_set="PCL_@_FMB_TiesNodes" rb="4" /> <rigid name="PTC_With_PTB" node_set="PTC_@_PTB_TiesNodes" rb="1" /> <rigid name="PTL_With_PTB" node_set="PTL_@_PTB_TiesNodes" rb="1" /> <rigid name="ACL_With_TBB" node_set="ACL_@_TBB_TiesNodes" rb="3" /> <rigid name="ACL_With_FMB" node_set="ACL_@_FMB_TiesNodes" rb="4" /> <rigid name="LCL_With_FBB" node_set="LCL_@_FBB_TiesNodes" rb="2" /> <rigid name="MCL_With_TBB" node_set="MCL_@_TBB_TiesNodes" rb="3" /> <rigid name="MCL_With_FMB" node_set="MCL_@_FMB_TiesNodes" rb="4" /> <rigid name="PTL_With_TBB" node_set="PTL_@_TBB_TiesNodes" rb="3" /> <rigid name="MNS-L_With_TBB" node_set="MNS-L_@_TBB_TiesNodes" rb="3" /> <rigid name="MNS-M_With_TBB" node_set="MNS-M_@_TBB_TiesNodes" rb="3" /> <rigid name="LCL_With_FMB" node_set="LCL_@_FMB_TiesNodes" rb="4" /> <rigid name="TBC-M_With_TBB" node_set="TBC-M_@_TBB_TiesNodes" rb="3" /> <rigid name="FMC_With_FMB" node_set="FMC_@_FMB_TiesNodes" rb="4" /> <rigid_body mat="1"> <fixed bc="Rx" /> <fixed bc="Ry" /> <fixed bc="Rz" /> <fixed bc="y" /> <fixed bc="x" /> <fixed bc="z" /> </rigid_body> <rigid_body mat="2"> <fixed bc="Rx" /> <fixed bc="Ry" /> <fixed bc="Rz" /> <fixed bc="y" /> <fixed bc="x" /> <fixed bc="z" /> </rigid_body> <rigid_body mat="3"> <fixed bc="Rx" /> <fixed bc="Ry" /> <fixed bc="Rz" /> <fixed bc="y" /> <fixed bc="x" /> <fixed bc="z" /> </rigid_body> <rigid_body mat="4"> <fixed bc="Rx" /> <fixed bc="Ry" /> <fixed bc="Rz" /> <fixed bc="y" /> <fixed bc="x" /> <fixed bc="z" /> </rigid_body> <!--<prescribed bc="z" lc="1"> 1.0</prescribed> --> <!--<force bc="z" lc="1"> 10.0</force> --> </Boundary> <Contact> <contact surface_pair="QAT_To_FMB" type="facet-to-facet sliding"> <penalty> 1</penalty> <maxaug> 10</maxaug> <auto_penalty> 1</auto_penalty> <laugon> 1</laugon> </contact> <contact surface_pair="QAT_To_FMC" type="facet-to-facet sliding"> <penalty> 1</penalty> <maxaug> 10</maxaug> <auto_penalty> 1</auto_penalty> <laugon> 1</laugon> </contact> <contact surface_pair="TBC-L_To_FMC" type="facet-to-facet sliding"> <penalty> 1</penalty> <maxaug> 10</maxaug> <auto_penalty> 1</auto_penalty> <laugon> 1</laugon> </contact> <contact surface_pair="TBC-L_To_MNS-L" type="facet-to-facet sliding"> <penalty> 1</penalty> <maxaug> 10</maxaug> <auto_penalty> 1</auto_penalty> <laugon> 1</laugon> </contact> <contact surface_pair="PCL_To_TBB" type="facet-to-facet sliding"> <penalty> 1</penalty> <maxaug> 10</maxaug> <auto_penalty> 1</auto_penalty> <laugon> 1</laugon> </contact> <contact surface_pair="PCL_To_FMB" type="facet-to-facet sliding"> <penalty> 1</penalty> <maxaug> 10</maxaug> <auto_penalty> 1</auto_penalty> <laugon> 1</laugon> </contact> <contact surface_pair="PCL_To_MNS-M" type="facet-to-facet sliding"> <penalty> 1</penalty> <maxaug> 10</maxaug> <auto_penalty> 1</auto_penalty> <laugon> 1</laugon> </contact> <contact surface_pair="PCL_To_ACL" type="facet-to-facet sliding"> <penalty> 1</penalty> <maxaug> 10</maxaug> <auto_penalty> 1</auto_penalty> <laugon> 1</laugon> </contact> <contact surface_pair="PTC_To_FMC" type="facet-to-facet sliding"> <penalty> 1</penalty> <maxaug> 10</maxaug> <auto_penalty> 1</auto_penalty> <laugon> 1</laugon> </contact> <contact surface_pair="ACL_To_TBB" type="facet-to-facet sliding"> <penalty> 1</penalty> <maxaug> 10</maxaug> <auto_penalty> 1</auto_penalty> <laugon> 1</laugon> </contact> <contact surface_pair="ACL_To_FMB" type="facet-to-facet sliding"> <penalty> 1</penalty> <maxaug> 10</maxaug> <auto_penalty> 1</auto_penalty> <laugon> 1</laugon> </contact> <contact surface_pair="MCL_To_TBB" type="facet-to-facet sliding"> <penalty> 1</penalty> <maxaug> 10</maxaug> <auto_penalty> 1</auto_penalty> <laugon> 1</laugon> </contact> <contact surface_pair="MCL_To_FMB" type="facet-to-facet sliding"> <penalty> 1</penalty> <maxaug> 10</maxaug> <auto_penalty> 1</auto_penalty> <laugon> 1</laugon> </contact> <contact surface_pair="MNS-L_To_FMC" type="facet-to-facet sliding"> <penalty> 1</penalty> <maxaug> 10</maxaug> <auto_penalty> 1</auto_penalty> <laugon> 1</laugon> </contact> <contact surface_pair="MNS-M_To_TBC-M" type="facet-to-facet sliding"> <penalty> 1</penalty> <maxaug> 10</maxaug> <auto_penalty> 1</auto_penalty> <laugon> 1</laugon> </contact> <contact surface_pair="MNS-M_To_FMC" type="facet-to-facet sliding"> <penalty> 1</penalty> <maxaug> 10</maxaug> <auto_penalty> 1</auto_penalty> <laugon> 1</laugon> </contact> <contact surface_pair="LCL_To_FMB" type="facet-to-facet sliding"> <penalty> 1</penalty> <maxaug> 10</maxaug> <auto_penalty> 1</auto_penalty> <laugon> 1</laugon> </contact> <contact surface_pair="TBC-M_To_FMC" type="facet-to-facet sliding"> <penalty> 1</penalty> <maxaug> 10</maxaug> <auto_penalty> 1</auto_penalty> <laugon> 1</laugon> </contact> <contact surface_pair="MCL_To_MNS-M" type="tied"> <penalty> 10</penalty> </contact> </Contact> <LoadData> <loadcurve id="1" type="linear"> <point> 0,0</point> <point> 1,1</point> </loadcurve> <loadcurve id="2" type="step"> <point> 0.00,0.1</point> <point> 0.10,0.1</point> <point> 0.20,0.1</point> <point> 0.30,0.1</point> <point> 0.40,0.1</point> <point> 0.50,0.1</point> <point> 0.60,0.1</point> <point> 0.70,0.1</point> <point> 0.80,0.1</point> <point> 0.90,0.1</point> <point> 1.00,0.1</point> </loadcurve> </LoadData> <Output> <plotfile type="febio"> <var type="contact gap" /> <var type="contact pressure" /> <var type="contact traction" /> <var type="displacement" /> <var type="reaction forces" /> <var type="stress" /> <!--<var type="prestrain stretch"/> --> </plotfile> </Output> <Step name="Initialize"> <Control type="febio"> <time_steps> 10</time_steps> <step_size> 0.1</step_size> <time_stepper> <dtmin> 1e-10</dtmin> <dtmax lc="2"> 0.1</dtmax> <max_retries> 20</max_retries> <opt_iter> 10</opt_iter> <aggressiveness> 1</aggressiveness> </time_stepper> <optimize_bw> 1</optimize_bw> <plot_level> PLOT_MUST_POINTS</plot_level> <analysis type="static" /> <min_residual> 1e-10</min_residual> </Control> <!--<Constraints> <constraint type="prestrain"> <tolerance> 0.01</tolerance> <update> 1</update> </constraint> </Constraints> --> </Step> <Step name="LoadingStep"> <Control type="febio"> <time_steps> 10</time_steps> <step_size> 0.1</step_size> <time_stepper> <dtmin> 1e-10</dtmin> <dtmax lc="2"> 0.1</dtmax> <max_retries> 20</max_retries> <opt_iter> 10</opt_iter> <aggressiveness> 1</aggressiveness> </time_stepper> <optimize_bw> 1</optimize_bw> <plot_level> PLOT_MUST_POINTS</plot_level> <analysis type="static" /> <min_residual> 1e-10</min_residual> </Control> </Step> </febio_spec>