%% Nass height
stewart = struct();
stewart.h = 90; %mm
stewart.jacobian = 174.5; %mm

%% Bottom Plate
BP = struct();
BP.rad.int = 0; %mm
BP.rad.ext = 150; %mm
BP.thickness = 10; % mm
BP.leg.rad = 100; %mm
BP.leg.ang = 5; %deg
BP.density = 8000; %kg/m^3
BP.color = [0.5 0.5 0.5]; %rgb

%% TOP Plate
TP = struct();
TP.rad.int = 0;%mm
TP.rad.ext = 100; %mm
TP.thickness = 10; % mm
TP.leg.rad = 90; %mm
TP.leg.ang = 5; %deg
TP.density = 8000; %kg/m^3
TP.color = [0.5 0.5 0.5]; %rgb

%% Leg
Leg = struct();
Leg.stroke = 80e-6; % m
Leg.rad.bottom = 8; %mm
Leg.rad.top = 5; %mm
Leg.sphere.bottom = 10; % mm
Leg.sphere.top = 8; % mm
Leg.density = 8000; %kg/m^3
Leg.lenght = stewart.h; % mm (approximate)
Leg.m = Leg.density*2*pi*((Leg.rad.bottom*1e-3)^2)*(Leg.lenght*1e-3); %kg
Leg.color.bottom = [0.5 0.5 0.5]; %rgb
Leg.color.top = [0.5 0.5 0.5]; %rgb
Leg.k.ax = 5e7; % N/m
Leg.ksi.ax = 10;
Leg = updateDamping(Leg);


%% Sphere
SP = struct();
SP.thickness.bottom = 1; %mm
SP.thickness.top = 1; %mm
SP.rad.bottom = Leg.sphere.bottom; %mm
SP.rad.top = Leg.sphere.top; %mm
SP.height.bottom = 5; %mm
SP.height.top = 5; %mm
SP.density.bottom = 8000; %kg/m^3
SP.density.top = 8000; %kg/m^3
SP.m = SP.density.bottom*2*pi*((SP.rad.bottom*1e-3)^2)*(SP.height.bottom*1e-3); %kg
SP.color.bottom = [0.5 0.5 0.5]; %rgb
SP.color.top = [0.5 0.5 0.5]; %rgb
SP.k.ax = 1e2; % N*m/deg
SP.ksi.ax = 10;
SP = updateDamping(SP);

%%
function element = updateDamping(element)
    field = fieldnames(element.k);
    for i = 1:length(field)
        element.c.(field{i}) = 1/element.ksi.(field{i})*sqrt(element.k.(field{i})/element.m);
    end
end