Trang chủ

Nghiên cứu động lực học đường cáp vận chuyển trái thanh long ở vùng Tây Nam Bộ - 23


[M_i,M_j] = min(Bang_ghi)


Diem_thap = zeros() for i=1:5

Diem_thap(i,1)= Bang_ghi(M_j(4+i),2); Diem_thap(i,2)= Bang_ghi(M_j(4+i),4+i);


end


goc = 1; toi= size(Bang_ghi,1)


Có thể bạn quan tâm!

Xem toàn bộ 196 trang tài liệu này.

Tieu_de_cua = { [ ' ' num2str(So_cua(1)) ' Chuyen huong'] [' ' num2str(So_cua(2)) ' Chuyen huong'] [' ' num2str(So_cua(3)) ' Chuyen huong'] [' ' num2str(So_cua(4)) ' Chuyen huong'] [' ' num2str(So_cua(5)) ' Chuyen huong'] [ ' Cong suat min'] };


Nghiên cứu động lực học đường cáp vận chuyển trái thanh long ở vùng Tây Nam Bộ - 23

plot(Bang_ghi(goc:toi,2),Bang_ghi(goc:toi,5),'R- ',Bang_ghi(goc:toi,2),Bang_ghi(goc:toi,6),'B--',...

Bang_ghi(goc:toi,2),Bang_ghi(goc:toi,7),'R-.'); hold on

plot(Bang_ghi(goc:toi,2),Bang_ghi(goc:toi,8),'k:','linewidth',2); hold on

plot(Bang_ghi(goc:toi,2),Bang_ghi(goc:toi,9),'r: +'); hold on

plot(Diem_thap(:,1),Diem_thap(:,2), 'k*') ;

%text(2300,370,'* - Cong suat min');

xlabel('Luc cang ngang H (N)'); ylabel('Cong suat tieu thu (W)'); legend(Tieu_de_cua);

grid on hold off

end % end if k_ve =1

%&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& Cac ham con &&&&&&&&&&&&&&&&&&&&&&&&&


%++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

function Nghiem = Tinh_vong(Nhap0) S0=Nhap0(1) ; h=Nhap0(2); L=Nhap0(3); q0=Nhap0(4); tgA= 2*L/S0;

a=sqrt(tgA^2-2*tgA) ; b=sqrt(tgA^2-1);


if L > S0

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


tgu0=hamU(a,b);


% Tinh luc H :


tgH = 0.5*q0*S0 /tgu0; Nghiem(1)=tgu0; Nghiem(2)=S0; Nghiem(3)= tgH;


% ==========

end end

%+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

++++++++


function Nghiem_Lech = Tinh_vong_lech(Nhap0) Nhap01=Nhap0;


S0=Nhap01(1) ; h=Nhap01(2); L = Nhap01(3); q0=Nhap01(4);


u0=2*h/S0; u1= u0/2; u2 = sqrt(16*L^2 -S0^2)/S0;


L0=0.5*S0*(sqrt(1+u0^2)+log(u0+sqrt(1+u0^2))/u0); LKT = sqrt(S0^2+h^2);


if L > LKT

%==========


if L < L0

a=u1; b=u0 ; end


if L > L0

a=u0; b=u2;

end


if L==L0

tgu0=u0; else

tgu0=hamU_lech(a,b);

end


% Tinh luc H :


tgLi = S0*tgu0 /(tgu0 - h/S0); tgH = 0.5*q0*tgLi /tgu0; Nghiem_Lech(1)=tgu0; Nghiem_Lech(2)=tgLi; Nghiem_Lech(3)= tgH;


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

end end


%+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

++++++++

function hamf =ham(x)


tgf=(2*h/S0 - x)/x ;

tg0 = tgf*sqrt(1+(x*tgf)^2)-log(x*tgf+sqrt(1+(x*tgf)^2))/x;

hamf =0.5*S0*(sqrt(1+x^2)+log(x+sqrt(1+x^2))/x - tg0)-(1-tgf)*L; end


%+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

++++++++

function hamf1 =ham1(x)


hamf1 =sqrt(1+x^2)+log(x+sqrt(1+x^2))/x-2*L/S0;

end


%++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

+++

function hamU1 = hamU_lech(a01,b01) a0=a01; b0=b01;

KT = abs(b0-a0); while KT > saiso c=(a0+b0)/2;

f1= ham(c); f2= ham(b0);

B01=f1*f2; if B01<0

a0 =c;

end

if B01 >=0

b0=c;


end

end

KT = abs(b0-a0);


hamU1=(a0+b0)/2;


end


%++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

+++


function U1 = hamU(a01,b01) a0=a01; b0=b01;

KT = abs(b0-a0); while KT > saiso c=(a0+b0)/2;


f1= ham1(c);

f2= ham1(b0);

B01=f1*f2; if B01<0

a0 =c;

end

if B01 >=0

b0=c;


end

end

KT = abs(b0-a0);

U1=(a0+b0)/2;

end


%++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

+++

end


5) Tính toán dao độ ng khi chuyể n hướng

function Dao_dong_khi_chuyen_huong_chuan clc;

Ten_ghi = 'SoLieu_EXCEL_MATLAB.xlsx';

% Hình th?

H=3500; % N Luc cang ngang L=0.30; % m do dài dây treo gio

S= 0.80; % m khoang cách giua hai gio M =2; % kg khoi luong gio hàng g=9.80; % m/s^2

% gia_toc= 5.*rand; van_toc = 0.40; % m/s R = 0.60;

t_start = 0;

t_end = 57; %final time in seconds. time_span =t_start:0.001:t_end; [t,x]=ode45(@rhs,time_span,[0 0]); [dong,cot]= size(x)

%************************************** subplot(1,1,1);

cla;

toi=round(dong/3) plot(t(:),L*x(:,1),'B-','linewidth',2);

hold on

xlabel("Thoi gian (s)")

ylabel("Bien do Dao dong chuyen huong (m)") grid on

%************************************** function xdot=rhs(t,x)

dxdt_1 = x(2);

dxdt_2 =(-g/L)*x(1)+((van_toc^2/(L*R+L^2*x(1)))+ forceZ(t))*(1- 0.5*x(1)^2);


xdot=[dxdt_1; dxdt_2];

end

%********************

% The forcing function, edit to change as needed.

%******************** function fz=forceZ(t)

% P_gio = 5.*rand -2;


end

a=5.6 ;

Fz = 2*sin(a*t); fz= Fz/(M*L);

end

6) Tính biên độ dao động cực đại ngang và dọc của

giỏ


function Dao_dong_ngang_doc_chuan

clc; % Tinh dao dong Ngang - Doc và Bien do Max Ten_ghi = 'SoLieu_EXCEL_MATLAB.xlsx';


n=31; % so luong gio trên mot nhip cáp H=3500; % N Luc cang ngang

L=0.2; % m do dài dây treo gio

S= 0.80; % m khoang cách giua hai gio M =2; % kg khoi luong gio hàng g=9.80; % m/s^2

gia_toc0= 1.*rand Bien_Fz=3.*rand;

m=n-1;


Dao_dong_Ngang=1;

Dao_dong_Doc=1;

if Dao_dong_Ngang==1

Dao_dong_Doc=0;

end Bien_do_Max=1 ;

if Bien_do_Max==1

Dao_dong_Ngang=0; Dao_dong_Doc=0;

end

Gio = round(m/2) Gio2= round(m/3) Gio1= round(m/4)

kx=0;


L_day=zeros(); H_ngang = zeros();


for k=1:5

L_day(k) = 0.2+(k-1)*0.05;

end

for k=1:5

H_ngang(k)= 1900+ (k-1)*400; end

Mang_ve=zeros();

%&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& for k_H=1:1

b= (H/(M*g))*(L/S) A=zeros(m);

K_ve = rem(m,2) ; ve = m;

if K_ve == 0

ve = m+1;

end

for i=1:m

A(i,i)= 2*b+1; end

for i=1: m-1 A(i,i+1)=-b;

A(i+1,i)=-b;

end Son=inv(A);

K=zeros(4*m); C = zeros(4*m,1); B = zeros(4*m,1); x0=zeros(1,4*m);

for i=1:m

K(2*i-1,2*i)=1;

end

for i=1:m

B(2*i,1)=1;

end

x1 = zeros(); Z=zeros(); M_1=zeros(); M1=zeros(); D=zeros(); t_start = 0;

t_end = 100; %final time in seconds. time_span =t_start:0.001:t_end;

[t,x]=ode45(@rhs,time_span,x0); cot = size(x,2)

dong = (size(x,1)-1) Goc=1; Toi = dong;

so_tinh = Toi-Goc+1; end % for k_H

%&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&

% TTTTTTTTTTTTTT het thay doi H và L


% #########################################################

% Xet cac bien do Max tai cac gio


if Bien_do_Max==1

i=1 ; k =0; chua=zeros();


while i < dong-1 i=i+1;

a1 = x(i-1,ve); a2=x(i,ve); % a3 =x(i+1,3);

%if (a2>0) & (a2>a1) & (a2>a3) if (a2 * a1) < 0

k=k+1; chua(k,1) = a2;

end

end

so_dao_dong = size(chua,1) [M_i,M_j] =max(abs(x)); MaxA = zeros();

for k=1:m

MaxA(k,1) = k; MaxA(k,2) = M_i(2*k-1); end

plot(MaxA(:,1),100*MaxA(:,2),'+');

xlabel(" Thu tu gio ")

ylabel("Bien do Max cua Dao dong Ngang (cm)") grid on


%************************************** ghi_bang_Bien_do_Max = 1

if ghi_bang_Bien_do_Max==1 Bang_ghi=zeros();


Tieu_de1 ={'Bien do Max cua cac gio hang'}; xlswrite(Ten_ghi, Tieu_de1, 4, [ 'B' num2str(1)]);

for i=1:m-1


end

Bang_ghi(1,i) = i; Bang_ghi(2,i) = MaxA(i,2);

xlswrite(Ten_ghi,MaxA, 4, [ 'B' num2str(3)]);

end

%************************************** end % if Bien_do_Max==1


% #########################################################

% Ve_do_thi_Doc_Ngang = 1;

if (Dao_dong_Ngang==1) | (Dao_dong_Doc==1)

so_gio(1) = 1; so_gio(2) = 8; so_gio(3) =round(m/2); Thu_tu_gio = { ['Gio ' num2str(so_gio(3))]}

if Dao_dong_Ngang==1

cot_x(1)= 2*so_gio(1)-1; cot_x(2)= 2*so_gio(2)-1; cot_x(3)= 2*so_gio(3)-1;


Goc = 1;

plot(t(Goc:Toi),x(Goc:Toi,cot_x(3)),'B-','linewidth',2); hold on

legend(Thu_tu_gio); xlabel("Thoi gian (s)")

ylabel("Bien do Dao dong Ngang (m)") grid on

end

if Dao_dong_Doc==1 Toi= round(dong/3);

cot_x(1)= 2*m+2*so_gio(1)-1; cot_x(2)= 2*m+2*so_gio(2)-1; cot_x(3)=2*m+ 2*so_gio(3)-1;


plot(t(Goc:Toi),x(Goc:Toi,cot_x(3)),'B-','linewidth',2); legend(Thu_tu_gio);

hold off

xlabel("Thoi gian (s)")

ylabel("Bien do Dao dong Doc (radian)") grid on

end end

%************************************** function xdot=rhs(t,x)


for i=1:4*m

C(i,1)=x(i);

end

for i=1:m

x1(i,1)=C(2*i-1,1);

end


for i=1:2*m

x2(i,1)= C(2*m+ i,1);

end Z=Son*x1;

..... Xem trang tiếp theo?
⇦ Trang trước - Trang tiếp theo ⇨

Ngày đăng: 21/07/2022