Homec4science
Diffusion MSM (master) Browse

msm/ShoulderModel_HHT/postprocessing_initilize_contact_model_glenoid.m
8406fed2a624master

Last Change
Raw File

postprocessing_initilize_contact_model_glenoid.m
View Options

function PPDATA = postprocessing_initilize_contact_model_glenoid(CMDATA)
%{
Function for initializing the contact on the glenoid.
--------------------------------------------------------------------------
Syntax :
PPDATA = postprocessing_initilize_contact_model_glenoid(CMDATA)
--------------------------------------------------------------------------
Function Description :
This function spans the internal surface of the glenoid to define the
points h_i on the surface. It also defines the constant values of the span
step size for this illustration.
--------------------------------------------------------------------------
%}
% define the position of the hi in the humerus frame {H}:
PPDATA.SpanStepSize = 2*pi/5; % spanning step size
theta = [0:PPDATA.SpanStepSize:2*pi]; % longtitude 0<theta<2pi on the xy from x
phi = [0:PPDATA.SpanStepSize:pi]; % latitude 0<phi<pi from z
% define the top and bottom points of the sphere (and their associated
% surafce areas)
PPDATA.h{1,1} = CMDATA.r_g*[cos(theta(1))*sin(phi(1));
cos(phi(1));
sin(theta(1))*sin(phi(1))];
PPDATA.a(1,1) = size(theta,2)*(0.5*CMDATA.r_g^2*PPDATA.SpanStepSize^2*...
sin(phi(2))*cos(PPDATA.SpanStepSize/2));
PPDATA.h{size(theta,2),size(phi,2)} = CMDATA.r_g*[cos(theta(1))*sin(phi(size(phi,2)));
cos(phi(size(phi,2)));
sin(theta(1))*sin(phi(size(phi,2)))];
PPDATA.a(size(theta,2),size(phi,2)) = PPDATA.a(1,1);
% span the whole surface of the humeral head and define the center points
% of the mesh grids. And rightafter their associated surace areas.
for i = 1:size(theta,2)-1 % theta = 0 and 2pi result in the same points
for j = 2:size(phi,2)-1 % exclude the top and bottom points
PPDATA.h{i,j} = CMDATA.r_g*[cos(theta(i))*sin(phi(j));
cos(phi(j));
sin(theta(i))*sin(phi(j))]; % hi in the {H} frame
PPDATA.a(i,j) = CMDATA.r_g^2*sin(phi(j))*PPDATA.SpanStepSize^2; % surface area of the mesh grid [m^2]
end
end
return;
c4science · Help