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ộ - 22

end

Li = L1*u /(u - h/L1); H = 0.5*q*Li /u;

Ryc = H*u;

Ryt = H*u*abs(1-2*L1/Li); Tc=H*sqrt(1+u^2);

Tt= H*sqrt(1+(1-2*L1/Li)^2); f0 = (Li*u/4);

KetQuaf(i+1,j+1)= f0; KetQuaRyc(i+1,j+1)= Ryc; KetQuaRyt(i+1,j+1)= Ryt; KetQuaH(i+1,j+1)= H;

KetQuaLi(i+1,j+1)= Li;

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

end

if L <= LKT

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

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

bay0= ' Khong du do dai cua L , L > ' ; son0= num2str(LKT);

bay00=[bay0 son0];

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ộ - 22

%bay=[bay LKT] disp(bay00 ); end

end % if h > 0

%$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ end

end

Bay = floor(KetQuaf) Son =floor(KetQuaLi)

for j=2:5

Bay1(1, (j-2)*3+2)= KetQuaf(1,j); end

for i=2: 4

Bay1(i+1, 1)= KetQuaf(i,1); for j=2:5

Bay1(i+1, (j-2)*3+2)= KetQuaf(i,j) -Mangh(j-1); Bay1(i+1, (j-2)*3+3)= KetQuaRyc(i,j);

Bay1(i+1, (j-2)*3+4)= KetQuaRyt(i,j); end

end

Do_vong_Max = max(Bay1(2:4, 1))

xlswrite(Ten_ghi, {'Tinh cap lech'}, 'Sheet9', [ 'A' num2str(1)]); xlswrite(Ten_ghi, Bay1, 'Sheet3', [ 'A' num2str(3)]); xlswrite(Ten_ghi, {'Do dai day (cm)'}, 'Sheet3', [ 'A'

num2str(4)]);

xlswrite(Ten_ghi, {'h (cm)'}, 'Sheet3', [ 'B' num2str(2)]); for j=1:4

xlswrite(Ten_ghi, {'f(cm)'}, 'Sheet3', [ setstr((j-1)*3+66) num2str(4)]);

xlswrite(Ten_ghi, {'Rc (N)'}, 'Sheet3', [ setstr((j-1)*3+67) num2str(4)]);

xlswrite(Ten_ghi, {'Rt (N)'}, 'Sheet3', [ setstr((j-1)*3+68) num2str(4)]);

end

xlswrite(Ten_ghi, {''}, 'Sheet3', [ 'A' num2str(3)]);

%+++++++++++++++++++++++++++++++++++ VE DO THI++++++++++++++++++++

Duong_day = 1; Cao_h =1;

f0 = Bay(1+Duong_day,1+ Cao_h); Li = Son(1+Duong_day,1+ Cao_h);

Nve = 100; Y0 = f0; X0 = 20;

X = zeros() ; Y = zeros(); delX = L1/Nve;

X(1) = 0;

% Y(1) = Y0;

Y(1) =0;

for i =1: Nve

X(i+1) = X(i) + delX;

Y(i+1) = - (4*f0*X(i+1)*(Li-X(i+1))/Li^2);

end plot(X,Y);

axis([0 2600 -(Y0+5) 0 ]);

axis on

xlabel(["Chieu dai nhip cap (cm)" ; "( h= 10 cm , L = 2501 cm) "]); ylabel('Do vong nhip cap (cm)');

grid on

%+++++++++++++++++++++++++++++++++++ HET VE DO THI++++++++++++++++++++

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

function hamf =ham(x) tgf=(2*h/L1 - x)/x ;

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

end

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

++++++++

function hamf1 =ham1(x)

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

end

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

++++++++

function hamU1 = hamU_lech(a0,b0) 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

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

++++++++

end

3) Tính độ dãn dây cáp

function Tinh_do_dan_day_chuan


clc;

Ten_ghi = [ 'SoLieu_EXCEL_MATLAB.xlsx']; E = 1980*10^4 ; % N/cm2

F = 0.28; % cm2

L1 = 2400 ; % L0 = 2510 ; cm

q = 0.05 ; % N/cm sai_so = 10^(-8); nq=4; nL=3; Lap =5;

qk=zeros(); Lk =zeros();

qk(1) =0.08 ; qk(2)=0.1; qk(3)=0.15; qk(4)=0.2; % Tinh voi 4

gia tri q

for k=1:3

Lk(k)=L1+ 0.4 +(k-1)*0.3;

end

Bang=zeros(3,2*nq+1); Bang_sai=zeros(3,2*nq+1); Tit1(1,1)={'L0(cm)'}; Titq(1,1)={'q(kg/cm)'};

for i=1:nq

Tit1(1,2*i)={'L (cm)'}; Tit1(1,2*i+1)={'f (cm)'};

ch=num2str(qk(i)); Titq(1,2*i)= { ch};

end

for i=1:nL L0=Lk(i);

Bang(i,1)= L0; Bang_sai(i,1)=L0;

for j=1:nq

q=qk(j);

L=L0;

for k=1:Lap

A0= 2*L/L1;

a= sqrt(A0^2-2*A0); b = sqrt(A0^2-1); kt = abs(a-b);

while kt > sai_so c= 0.5*(a+b);

fb = sqrt(1+b^2) + log(b+sqrt(1+b^2))/b - A0; fc= sqrt(1+c^2) + log(c+sqrt(1+c^2))/c - A0; if fb*fc < 0

a=c;

end

if fb*fc>=0 b =c;

end

kt = abs(a-b); end

u = (b+a)/2;

del_L = 0.5*q*L1^2*(1/u+u/3)/(E*F); L=del_L + L0;

% disp(L*1000 );

s = sprintf( '%6f', L); f = u*L1/4;

f0 = sprintf( '%6f', f); if k>1

s_L = abs(L_s-L); s_f = abs(f_s-f); end


L_s=L; f_s=f; end

Bang(i,2*j)= L; Bang(i,2*j+1)=f; Bang_sai(i,2*j)= s_L; Bang_sai(i,2*j+1)=s_f;


end


end


cong = ['Tinh do dan day ung voi do dai nhip ' num2str(L1) ' cm']; xlswrite(Ten_ghi, {cong}, 'Sheet2', ['A' num2str(1)]); xlswrite(Ten_ghi, Titq, 'Sheet2', ['A' num2str(3)]); xlswrite(Ten_ghi, Tit1, 'Sheet2', ['A' num2str(4)]);

xlswrite(Ten_ghi, Bang, 'Sheet2', ['A' num2str(5)]); xlswrite(Ten_ghi, Bang_sai, 'Sheet2', ['A' num2str(11)]);


xlswrite(Ten_ghi, {'Sai so tuong ung '}, 'Sheet2', ['B' num2str(9)]);

end

4) Tính công suấ t

function tinh_cong_suat clc;

%============= Doc tu file EXCEL

=========================================

Ten_ghi = 'SoLieu_EXCEL_MATLAB.xlsx'; [num txt]=xlsread(Ten_ghi);

So_cua=zeros(); Ton_hao_cua = zeros();

So_cua(1) = 4; So_cua(2) = 5; So_cua(3) = 6; So_cua(4) = 7;

So_cua(5) = 8;

Tong_S=num(4,1); f1=zeros(); N=num(2,1)

A=zeros(N); VP=zeros(); Nhap = zeros(N,4); Nhap0 = zeros(N,4); U=zeros(); q=zeros();d=zeros(); S=zeros(); d0=zeros();

B = zeros(); B1 = zeros(); f = zeros(); fw = zeros(); Li = zeros(); L1 =zeros();

DEL_U=zeros(N,1); VP0 = zeros(N,1); tgL=zeros(); tgS=0;

Gia_do =zeros(N,3);

R =zeros(); R1=zeros(); R2=zeros(); S0=zeros(); h=zeros(); L=zeros(); q0=zeros(); saiso = 10^(-6);

Nghiem01=zeros(N,3); for j=1:N

Nhap0(j,1)=num(4,j+2);

Nhap0(j,2)= num(6,j+2);

Nhap0(j,4)=num(5,j+2);

tgL(j) =sqrt(Nhap0(j,1)^2+Nhap0(j,2)^2); tgS=tgS+ tgL(j);

Nhap(j,1)=Nhap0(j,1);

Nhap(j,2)= abs(Nhap0(j,2));

Nhap(j,4)=Nhap0(j,4);

d(j)= Nhap(j,2); q(j)= Nhap(j,4); S(j) = Nhap(j,1); d0(j)= Nhap0(j,2);

end


Tong_S=tgS; Tong_L0=num(3,1) Bang_ghi = zeros();

%&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&( 1 ) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&

%for K_lap =1:1

for K_lap =1:15


Tong_L = (Tong_L0 + 5*K_lap) ;

% Tong_L = Tong_L0 ; Bang_ghi(K_lap,1) = Tong_L;

T = zeros(); H = zeros();


%

$$$$$$$$$$$$$$$$&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&========================= if (Tong_L > Tong_S) & (Tong_S > 0)

hs = Tong_L/Tong_S; for j=1:N

Nhap(j,3)= hs*tgL(j) ;

end

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

=======


%$$$$$$$$$$$ Tinh cac gia tri U ban dau 4444444444444444 for k=1:N

Nhap_vao = zeros(); for j=1: 4

Nhap_vao(j) = Nhap(k,j); end

if Nhap_vao(2)>0 Nghiem0=Tinh_vong_lech(Nhap_vao);

Nghiem01(k,1)= Nghiem0(1); Nghiem01(k,2)= Nghiem0(2); Nghiem01(k,3)= Nghiem0(3); U(k)= Nghiem0(1);

end % if h > 0


%$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$44444444444444 444444

if Nhap_vao(2)==0 Nghiem0=Tinh_vong(Nhap_vao);

Nghiem01(k,1)= Nghiem0(1); Nghiem01(k,2)= Nghiem0(2); Nghiem01(k,3)= Nghiem0(3);

U(k)= Nghiem0(1); end


end % for k=1:N va Tinh xong Nghiem U ban dau

%============ Nhap cac gia tri cua ma tran he so tinh theo Newton- raphson

for k=1:N-1

A(k,k)= 1/(q(k)*S(k));

A(k,k+1)= - 1/(q(k+1)*S(k+1));

end

%++++++++++=================== j=0; KTr =1;

while KTr > saiso j=j+1;

for k=1:N-1

tgd = d(k)/(q(k)*S(k)^2) -d(k+1)/(q(k+1)*S(k+1)^2);

VP(k) = U(k+1)/(q(k+1)*S(k+1))- U(k)/(q(k)*S(k)) +tgd ;

end

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

for k=1:N if d(k)==0

B(k)= -1; B1(k) =0;

end

if d(k)>0

B(k) = 2* d(k)/(U(k) *S(k)) -1 ;

B1(k) = -2* d(k)/(U(k)^2 *S(k)) ; end

end

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

for k=1:N

tg1 = sqrt(1+U(k)^2) + log(U(k)+sqrt(1+U(k)^2))/U(k);

tg2= - B(k)*sqrt(1+(B(k)*U(k))^2) -

log(U(k)*B(k)+sqrt(1+(B(k)*U(k))^2))/U(k); V(k) = tg1+tg2;

tg3 = sqrt(1+U(k)^2)/U(k) - log(U(k)+sqrt(1+U(k)^2))/U(k)^2;

tg4 = -(B(k)+2*U(k)*B1(k))*sqrt(1+(B(k)*U(k))^2)/U(k)

+log(U(k)*B(k)+sqrt(1+(B(k)*U(k))^2))/U(k)^2; V1(k) = tg3 + tg4;

end

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


tg1=0; tg2=0; for k=1:N

tg2 = 1/(S(k)*U(k)- d(k));

A(N,k) = 0.25*S(k)^2*( -d(k)* tg2^2*V(k) + U(k)*tg2*V1(k));

L1(k)=0.5*S(k)* V(k)/(1-B(k));

tg1=tg1+ L1(k); end

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

VP(N)= Tong_L-tg1;

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

Thu = det(A) G = inv(A);

DEl_U = G*VP';

% DEl_U = AVP';

U = U + DEl_U';

VP0= A*DEl_U;

tg3=0; tg4=0; for k=1:N tg3=tg3+VP0(k)^2; tg4 =tg4+ VP(k)^2; end

KTr = sqrt(tg3) end

%++++++++++========================= while KTr > saiso

for k=1:N

Li(k) = S(k)^2*U(k) /(S(k)*U(k) - d(k)); end

f= (Li.*U)/4;

f1=f-abs(d); for k=1:N

H(k)= 0.5* q(k)*Li(k)/U(k); end


for k=1:N-1


if d0(k) > 0

Gia_do(k,2)= H(k)*U(k)*abs(1-2*S(k)/Li(k)); Gia_do(k+1,1)= H(k)*U(k);


fw(k)= f(k)- d(k); end

if d0(k) == 0 Gia_do(k,2)= H(k)*U(k);

Gia_do(k+1,1)= H(k)*U(k); fw(k)= f(k);

end


if d0(k) < 0

Gia_do(k,2)= H(k)*U(k);

Gia_do(k+1,1)= H(k)*U(k)*abs(1-2*S(k)/Li(k)); fw(k)= f(k);

end


end

if d0(N) > 0

Gia_do(N,2)= H(N)*U(N)*abs(1-2*S(N)/Li(N)); Gia_do(1,1)= H(N)*U(N);

fw(N)= f(N)- d(N);

end

if d0(N) == 0 Gia_do(N,2)= H(N)*U(N);

Gia_do(1,1)= H(N)*U(N); fw(N)= f(N);

end

if d0(N) < 0

Gia_do(N,2)= H(N)*U(N);

Gia_do(1,1)= H(N)*U(N)*abs(1-2*S(N)/Li(N)); fw(N)= f(N);

end for k=1:N

Gia_do(k,3) = Gia_do(k,1) + Gia_do(k,2);

% T(k)=H(k)*sqrt(1+ q(k)^2*Li(k)^2/(4*H(k)^2));

T(k)= (0.5*q(k)*Li(k)/U(k))* sqrt(1+ U(k)^2);


R(k)= Gia_do(k,3); R1(k)= Gia_do(k,1);R2(k)= Gia_do(k,2); end

% Gia_do

% T

%+++++++================================

TongMasat =0; TongLuckeo=0; TaiTrong =0; for k=1: N

TongMasat = TongMasat+ R(k); % don vi N

TongLuckeo = TongLuckeo + q(k)*fw(k); % don vi N

end

VanToc = 2; % 2 km/h


VanToc = VanToc*1000/3600; % m/s TongMasat = TongMasat * 0.01;

Congsuat = TongLuckeo* VanToc + TongMasat*VanToc;


for i =1: 5

Ton_hao_cua(i) = H(1)*sqrt(2)*So_cua(i)*0.01*VanToc; Mang_conguat(i)= Ton_hao_cua(i)+Congsuat;

end


Bang_ghi(K_lap,2)= H(1); Bang_ghi(K_lap,3)= Congsuat;

Bang_ghi(K_lap,4)= Ton_hao_cua(1) ; for i=1:5

Bang_ghi(K_lap,4+i)= Congsuat+Ton_hao_cua(i) ; end

end % if (Tong_L > Tong_S) & (Tong_S > 0)


%

$$$$$$$$$$$$$$$$&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&========================= if (~((Tong_L > Tong_S) & (Tong_S > 0)))

bay= ' Do dai L cua day Cap phai > ' ; son=num2str(Tong_S);

bay=[bay son]; disp(bay ); end

end % end K_lap

%&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& &&&&&&&&&&&&&&&&&&&&&&&&& k_ve =1;

if k_ve ==1


k_ghi=0 % Ghi bang lien he H và Cong suat voi so chuyen huong if k_ghi==1

Tieu_de = {'L (cm)' ' H (N)' ' W_k ' 'W_c1' 'W_01' 'W02' 'W03' 'W04' 'W05'};

Tieu_de_cua = { num2str(So_cua(1)) num2str(So_cua(2)) num2str(So_cua(3)) num2str(So_cua(4)) num2str(So_cua(5))}


xlswrite(Ten_ghi, Tieu_de_cua, 1, [ 'G' num2str(30)]);


xlswrite(Ten_ghi, Tieu_de, 1, [ 'C' num2str(31)]); xlswrite(Ten_ghi, Bang_ghi, 1, [ 'C' num2str(32)]);


xlswrite(Ten_ghi, VanToc, 1, [ 'G' num2str(29)]); end

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

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