README for data in Chapter 2 in Developing parametric finite element models of the residual limb – prosthetic socket system towards clinical application


Author: J. W. Steer
Date: 19 June 2019

This contains instructions for opening the .pkl file containing the data behind the figures. The .pkl file is strutured using dictionaries

Opening the data file
---------------------

The data file is encoded using the pickle package in python 3.6. The data is all stored within numpy arrays. The authors recommend downloading the anaconda distribution of python. 

"""
import pickle
import numpy
with open('data.pkl', 'rb') as f:
	Data = pickle.load(f)
"""

Figure 3.6a
--------

Mesh convergence analysis

data = ['Residuum tip pressure, kPa', 
		'Fibula head pressure, kPa', 
		'Tibial tuberosity pressure, kPa',
		'Distal tibia tissue strain, %']


Example
Data['Figure 3.6a']['Residuum tip pressure, kPa']
Returns a numpy array of the RMSE error vs mesh size 

Figure 3.6b
--------

Mesh convergence analysis

data = ['Residuum tip pressure, kPa', 
		'Fibula head pressure, kPa', 
		'Tibial tuberosity pressure, kPa',
		'Distal tibia tissue strain, %']

Example
Data['Figure 3.6b']['Residuum tip pressure, kPa']
Returns a numpy array of the RMSE error vs solver time  

Figure 3.7
--------

Biomechanical outputs for different anatomical structures

data = ['Residuum tip pressure, kPa',
		'Fibula head pressure, kPa', 
		'Tibial tuberosity pressure, kPa', 
        'Posterior calf pressure, kPa', 
		'Posterior calf shear, kPa', 
		'Distal strain, %']

model = ['B_QuasiStatic', 'B_NoTendon', 'B_wSkin', 'B_wSkinFat']

Example
Data['Figure 3.7']['Residuum tip pressure, kPa']['Baseline']
Returns a numpy array of the biomechanical value for the '-2%', '0%', '+2%' and '+4%' socket designs

Figure 3.8
--------

Biomechanical outputs for different soft tissue models

data = ['Residuum tip pressure, kPa', 'Fibula head pressure, kPa', 
          'Tibial tuberosity pressure, kPa', 'Distal strain, %']
models = ['linear', 'hyperelastic']

Example
Data['Figure 3.8']['Residuum tip pressure, kPa']['linear']
Returns a numpy array of the biomechanical value for the soft tissue stiffness 

Figure 3.9
--------

Pressure distributions for different donning methods 

data = ['Overclosure\nPosition = 0mm', 'Overclosure\nPosition = 6mm',
          'Donnning\nPosition = 6mm', 'Donnning\nPosition = 30mm']


Example
Data['Figure 3.9']['Overclosure\nPosition = 0mm']
Returns a numpy array of the pressure distribution 

Figure 3.10
--------

Residuum tip pressure and tibial tuberosity shear with overclosure and donning for 4 socket designs

sockets = ['0%', '2%', '4%', '3% distal only']
data = ['Overclosure - Residuum Tip Pressure',
        'Overclosure - Tibial Tuberosity Shear',
        'Donning - Residuum Tip Pressure',
        'Donning - Tibial Tuberosity Shear']


Example
Data['Figure 3.10']['0%']['Overclosure - Residuum Tip Pressure']
Returns a numpy array of the interface stress vs initial position 

Figure 3.11
--------

Biomechanical outputs for difference sockets and quasi-static load cases 

data = ['Residuum tip pressure, kPa',
		'Fibula head pressure, kPa',
		'Tibial tuberosity pressure, kPa', 
		'Tibial tuberosity shear, kPa',
		'Posterior calf shear, kPa',
		'Distal strain, %']

socket = ['-2%', '0%', '+2%', '+4%']

Example
Data['Figure 3.11']['Residuum tip pressure, kPa']['0%']
Returns a numpy array of the biomechanical value for the Heel-strike, mid-stance and toe-off

Figure 3.12
--------

Pressure distribution for difference sockets and quasi-static load cases 


socket = ['-2%', '0%', '+2%', '+4%']

Example
Data['Figure 3.12']['-2%']
Returns a numpy array of the pressure for the Heel-strike, mid-stance and toe-off

Figure 3.13
--------

Biomechanical outputs for difference sockets and quasi-static load cases 

data = ['Residuum tip pressure, kPa', 
		'Fibula head pressure, kPa', 
		'Patella tendon pressure, kPa', 
		'Soft tissue strain, %']

socket = ['-2%', '0%', '+2%', '+4%']

Example
Data['Figure 3.13']['Residuum tip pressure, kPa']['0%']
Returns a numpy array of the biomechanical value across the gait cycle 

Figure 3.14
--------

Biomechanical outputs for difference sockets and quasi-static load cases 

data = ['Residuum tip pressure, kPa', 
		'Fibula head pressure, kPa', 
		'Patella tendon pressure, kPa', 
		'Soft tissue strain, %']


Example
Data['Figure 3.14']['Residuum tip pressure, kPa']
Returns a 5 x 5 numpy array of the biomechanical value with changing COF for the '-2%', '0%', '+2%' and '+4%' sockets

Figure 3.15
--------

Pressure profiles for COF = 0.3 and 0.7


Example
Data['Figure 3.15']
Returns a N x 2 numpy array of the pressure for 0.3 and 0.7