Реферат Моделювання руху голки в замку в'язального механізму

Тема: Курсовые проекты

діють в зазорі голка-клин і прискорення деформації для

% побудови графіковk=1: length (t)

L (k)=t (k) * v;

[y (k), vy (k), wy (k)]=ksi (L (k));

[F (k), Ft (k), P (k)]=sila (t (k), dlt (k), vdlt (k)); (k)=(F ( k) -Ft (k) -P (k) -m * gm * v * v * wy (k))/m;

end;

% Побудова графіків законів руху голки

figure (1) (311); plot (L, y, -k raquo ;, LineWidth raquo ;, 2); on, xlabel ( L, m ), ylabel ( xi, m ); (312); plot (L, vy * v, -k raquo ;, LineWidth raquo ;, 2); on, xlabel ( L, m ), ylabel ( xi prime, m/c ); ( 313); plot (L, wy * v ^ 2, -k raquo ;, LineWidth raquo ;, 2); on, xlabel ( L, m ), ylabel ( xi prime prime, m/c ^ 2 );

% Побудова графіків залежностей деформації, швидкості і ускоернія

% деформації в зазорі голка-клин (2)

subplot (311), plot (t, dlt, -k raquo ;, LineWidth raquo ;, 2), grid on ( t, c ), ylabel ( Delta , m ) (312), plot (t, vdlt, -k raquo ;, LineWidth raquo ;, 2), grid on ( t, c ), ylabel ( Delta prime, m/ c ) (313), plot (t, wdlt, -k raquo ;, LineWidth raquo ;, 2), grid on ( t, c ), ylabel ( Delta prime prime, m/c ^ 2 )

% Побудова графіка залежності сили в зазорі голка-клин

figure (3); (311), plot (t, Ft, -k raquo ;, LineWidth raquo ;, 2), grid on ( t, c ), ylabel ( Ft, H ) (312), plot (t, P, -k raquo ;, LineWidth raquo ;, 2), grid on ( t, c ), ylabel ( P, H ) (313), plot (t, F, -k raquo ;, LineWidth raquo ;, 2), grid on ( t, c ), ylabel ( F, H )

% ========================================= ================================

% підпрограми

% ========================================= ================================

% Підпрограма обчислення залежності ksi (L) і її першої та другої

% передавальних функцій [y, vy, wy]=ksi (L)

% Опис глобальних параметровL1 L2 L3 R1 R2 b1 z1 z2 z3

% Перебір участковL lt;=L1

% Перша ділянка

y=b1 ^ 2 * L ^ 2/(4 * z1);=b1 ^ 2 * L/(2 * z1);=b1 ^ 2/(2 * z1); L gt;=L2

% Якщо рух на третьому ділянці

y=z3 + (z3-z2-z1) * ((L-L3)/(L3-L2) - 1/pi * sin (pi * (L-L3)/(L3-L2 )));=(z3-z2-z1)/(L3-L2) * (1-cos (pi * (L-L3)/(L3-L2)));=(z3-z2-z1) * ( pi/(L3-L2) ^ 2) * sin (pi * (L-L3)/(L3-L2));

else

% Якщо рух на другій ділянці

y=z1 + b1 * (L-L1);=b1;=0;

% -----------------------------------------------------

% Підпрограма обчислення сил, що діють в моделі

function [F, Ft, P]=sila (t, dlt, vdlt)

% Опис глобальних параметровc b dlt1 m Ft0 P0 gv mu

% Обчислення сили F (dlt lt; 0)

F=-c * dlt-b * vdlt; (dlt gt;=0 amp; dlt lt;=dlt1)=0.0;=- c * (dlt-dlt1) -b * vdlt;

end

% Обчислення зовнішнього кінематичного впливу=v * t;

[y, vy, wy]=ksi (L);=v * vy + vdlt;