**************************************************************
** Geometry and Property Files
**************************************************************
**
** Bone Geometry and Properties
*INCLUDE, INPUT = BONE_TIBIA.inp
*INCLUDE, INPUT = BONE_FOOT.inp
*INCLUDE, INPUT = PROPS_BONES_TIBIA.inp
*INCLUDE, INPUT = PROPS_BONES_FOOT.inp
**
** Cartilage Geometry and Properties
*INCLUDE, INPUT = CARTILAGE_TIBIA.inp
*INCLUDE, INPUT = CARTILAGE_TALUS.inp
*INCLUDE, INPUT = CARTILAGE_PROPS.inp
**
**************************************************************
** Muscle Geometry and Properties
*INCLUDE, INPUT = MUSC_TIB_ANT.inp
*INCLUDE, INPUT = MUSC_SOLEUS.inp
**
** Joint Alignment/Connectors
*INCLUDE, INPUT = JOINT_CONFIGURATION.inp
**
** Contact Defititions
*INCLUDE, INPUT = CONTACT_INCLUSIONS.inp
**
**************************************************************
** Step 1
**************************************************************
**
*STEP
*DYNAMIC, EXPLICIT, DIRECT USER CONTROL
2.5E-5, 0.5
*OUTPUT,FIELD,VARIABLE=PRESELECT, TIME INTERVAL=0.025
**
**************************************************************
** Write files for Restart
**************************************************************
*RESTART, WRITE, NUMBER INTERVAL=1
**
**************************************************************
** Prescribe Kinematics
**************************************************************
**
*BOUNDARY, OP=NEW
** tibia rbrn
10, 1,6, 0.0
**
*** Ankle Kinematics *****************************************
*AMPLITUDE, VALUE=RELATIVE, SMOOTH=0.3, NAME=ANKLE_SI
0.0,	0.0
0.5,	5.0
*AMPLITUDE, VALUE=RELATIVE, SMOOTH=0.3, NAME=ANKLE_ML
0.0,	0.0
0.5,	2.0
**
*CONNECTOR MOTION, OP=NEW, AMPLITUDE=ANKLE_SI
Link_AnkleSI, 1, 1
*CONNECTOR MOTION, OP=NEW, AMPLITUDE=ANKLE_ML
Link_AnkleML, 1, 1
** Constrain other links in necessary degrees of freedom
*CONNECTOR MOTION, OP=NEW
Link_AnkleML, 4, 0
Link_AnkleAP, 4, 0
Link_AnkleAP, 1, 0
Link_AnkleSI, 4, 0
**
**************************************************************
** Muscle Forces
**************************************************************
**
*AMPLITUDE, VALUE=RELATIVE, SMOOTH=0.1, NAME=amp_tibant1
0.0, 5, 0.5, 5, 1.0, 5
*AMPLITUDE, VALUE=RELATIVE, SMOOTH=0.1, NAME=amp_soleus1
0.0, 5, 0.5, 5, 1.0, 5
**
*CONNECTOR LOAD, OP=NEW, AMPLITUDE=amp_tibant1
musc_tibant, 1, -1
*CONNECTOR LOAD, OP=NEW, AMPLITUDE=amp_soleus1
musc_soleus, 1, -1
**
**************************************************************
** Outputs
**************************************************************
**
*OUTPUT, HISTORY, TIME INTERVAL=0.025
*ELEMENT OUTPUT, ELSET=musc_soleus
CTF, CU
*ELEMENT OUTPUT, ELSET=musc_tibant
CTF, CU
*ELEMENT OUTPUT, ELSET=Link_AnkleML
CTF, CU
*ELEMENT OUTPUT, ELSET=Link_AnkleAP
CTF, CU
*ELEMENT OUTPUT, ELSET=Link_AnkleSI
CTF, CU
**
**************************************************************
*END STEP
**
**************************************************************
** Step 2
**************************************************************
**
*STEP
*DYNAMIC, EXPLICIT, DIRECT USER CONTROL
2.5E-5, 0.5
*OUTPUT,FIELD,VARIABLE=PRESELECT, TIME INTERVAL=0.025
**
***************************************************************
** Write files for Restart 
**
*RESTART, WRITE, NUMBER INTERVAL=1
*OUTPUT,FIELD,VARIABLE=PRESELECT, TIME INTERVAL=0.025
**
**************************************************************
** Contact
**************************************************************
**
*INCLUDE, INPUT = CART_CONTACT.inp
**
**************************************************************
** Prescribe Kinematics
**************************************************************
**
*BOUNDARY, OP=NEW
** tibia rbrn
10, 1,6, 0.0
**
*** Ankle Kinematics *****************************************
**
*CONNECTOR MOTION, OP=NEW
**
*CONNECTOR LOAD, OP=NEW
Link_AnkleML, 4, 0.0
*CONNECTOR LOAD, OP=NEW
Link_AnkleAP, 4, 0.0
*CONNECTOR LOAD, OP=NEW
Link_AnkleSI, 4, 0.0
*CONNECTOR LOAD, OP=NEW
Link_AnkleML, 1, 0.0
*CONNECTOR LOAD, OP=NEW
Link_AnkleAP, 1, 0.0
*CONNECTOR LOAD, OP=NEW
Link_AnkleSI, 1, 0.0
**
**************************************************************
** Muscle Forces
**************************************************************
**
*AMPLITUDE, VALUE=RELATIVE, SMOOTH=0.1, NAME=amp_tibant2
0.0, 5, 0.2, 10, 0.5, 10, 1.0, 10
*AMPLITUDE, VALUE=RELATIVE, SMOOTH=0.1, NAME=amp_soleus2
0.0, 5, 0.2, 5, 0.5, 5, 1.0, 5
**
*CONNECTOR LOAD, OP=NEW, AMPLITUDE=amp_tibant2
musc_tibant, 1, -1
*CONNECTOR LOAD, OP=NEW, AMPLITUDE=amp_soleus2
musc_soleus, 1, -1
**
**************************************************************
** Compressive Load to Bring into Contact
**************************************************************
**
*CLOAD, OP=NEW
4295633, 3, -10
4771823, 3, -10
4295633, 1, -5
4771823, 1, -5
**
**************************************************************
** Outputs
**************************************************************
**
*OUTPUT, HISTORY, TIME INTERVAL=0.025
*ELEMENT OUTPUT, ELSET=musc_soleus
CTF, CU
*ELEMENT OUTPUT, ELSET=musc_tibant
CTF, CU
*ELEMENT OUTPUT, ELSET=Link_AnkleML
CTF, CU
*ELEMENT OUTPUT, ELSET=Link_AnkleAP
CTF, CU
*ELEMENT OUTPUT, ELSET=Link_AnkleSI
CTF, CU
**
**************************************************************
*END STEP
**