wheel

(./wheel.m)

Help text

%%%
 Script file to generate the descripiton of a bicycle wheel.

Cross-Reference Information

plot_truss     ./plot_truss.m
solve_truss    ./solve_truss.m

Listing of ./wheel.m

%%%%
% Script file to generate the descripiton of a bicycle wheel.

dof = 6;    % Frame model
r = 20;        % Radius of the rim.
spokes = 32;
% nodes = [x,y,z;x,y,z;...]
     theta = ((1:spokes)*2*pi/spokes)';
    nodes = [r*sin(theta),r*cos(theta),zeros(spokes,1)];
end;
nodes = [0,0,1;0,0,-1;nodes]

[num_nodes,nc] = size(nodes)
active_nodes = [3:num_nodes];


%members = [:,member] where member == [pi,pj,a,e,beta,iss,itt,j,g]

arim = .1;
erim = 11000000 ; %psi
Iss_r = .01;
Itt_r = .005; 
beta = 0;
j_r = .01;
g_r = 4000000;%

members = [[3:num_nodes]',[4:num_nodes,3]',arim*ones(spokes,1),...
                erim*ones(spokes,1),zeros(spokes,1),...
                Iss_r*ones(spokes,1),Itt_r*ones(spokes,1),...
                j_r*ones(spokes,1),g_r*ones(spokes,1)];


mem2 = zeros(spokes,9);

a_spoke = .0025;
e_spoke = 29000000;

mem2(:,1) = ones(spokes,1);
mem2(:,2) = (3:num_nodes)';
mem2(:,3) = a_spoke*ones(spokes,1);
mem2(:,4) = e_spoke*ones(spokes,1);

mem2(2:2:spokes,1) = 2*mem2(2:2:spokes,1);

members = [members;mem2];

%loads = [node, fx,fy,fz,mx,my,mz;...
loads = [4,0,0,1,0,0,0];

displacements = solve_truss(nodes,members,active_nodes,loads,dof);
plot_truss(nodes,members,displacements);