%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%PROGRAMA PRINCIPAL SIMULACIÓN MOTOCICLETA%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%desde aqui se ejecutan todos los demás%
global MBody
global joy
joy = vrjoystick(1);
MBody = Motocicleta; %Carga parámetros del sistema en la estructura pcpal
Parametros; %los llama
nind = length(MBody.Ind);
n=nind;
y0ind=[0 0 0.97 0 0 0 0 0 0, 10 0 0 0 0 0 0 10/0.34 10/0.34]';%al simular contacto normal con disip y darle t, se ve la posic eq, algo mas profunda la trasera por cdg global
MBody.i=1; %SE ACTUALIZA EN ECMOV, COMO EL .t
%extrae del vector y0
q0=y0ind(1:n); %las n primeras soluc son q
v0 =y0ind(nind+1:2*nind); %las ultimas son v
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %%% SIMULACIÓN COORDENADAS %%%
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n = nind; %nº coordenadas
tic
At = 0.01; %increm de tiempo para RK4
MBody.At=At;
tspan = [0 30]; %intervalo de simulación
%ahora integra las ecs mov 2orden a partir de
%q0,v0, ec mov...para "todo" instante de t
[t y] = RungeKutta2(@EcMov_Baumgarte, tspan, [q0' v0']', At);
toc
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
