%% simulation parameters
% same as data aquisition
step_size = 1.25*10^-6;

% gravity
g = 9.81;

%% model parameters
% impactor
mass_impactor = 0.5078;
stiffness_impactor = 85490581.8233;

% force sensor
mass_9333a_half = 0.137 / 2;
stiffness_9333a = 2.3 * 10^9;

% impactor tip
mass_impactor_tip = 0.0838;
stiffness_impactor_tip = 13075858.7233;
damping_impactor_tip = 244.1288;

% hard stop
transition_region_hard_stop = 5e-05;

% femoral head
mass_head = 0.0621;

% friction element at taper junction
mu_kinetic_taper_junction = 0.49276;
mu_static_taper_junction = 0.78597;
normal_force_constant_taper_junciton = 5;
normal_force_factor_taper_junction = 144686080497.672;

% stem
mass_stem = 0.0434;
stiffness_stem = 134264491.3752;
mass_stem_adapter = 0.12065;

% soft tissue voigt model
mass_soft_tissue = 0.26306;
stiffness_soft_tissue = 4001;
damping_soft_tissue = 400.1;

%% input force
load('data');

%% simulation
model_name = 'model_force_transmission';
load_system(model_name);
in = Simulink.SimulationInput(model_name);

out = sim(in);