TÍNH TOÁN NHIỆT CÁC QUÁ TRÌNH
Quá trình nạp
1 ¿ 16,5−1 b.Hệ số khí sót γ r c.Nhiệt độ cuối quá trình nạp
AI Quá trình nén a.Tính n 1
● Tỷ nhiệt mol đẳng tích trung bình
● Tỷ nhiệt mol đẳng tích trung bình của sản phẩm cháy
- Khi α > 1 tíí́nh cho động cơ diesel theo công thức sau mc´ ” v = (19,867+ 1,634 α ) + 1 2 (427,38 + 184,36 α ).10 -5 T c (2)
● Tỷ nhiệt mol đẳng tích trung bình của hỗn hợp khí trong quá trình nén ´ ' mc´ v + γ r mc´ ô v m c vc = a' v = 19,84 b' v = 0,0047
● Xác định chỉ số nén đa biến trung bình n 1 : n 1 =
Thông thường để chọn n 1 ta chọn n 1 trong khoảng từ 1,340 ÷ 1,40
Thay dần các giá tri n 1 vào 2 vế của phương trình đến khi cân bằng, ta được: n 1 = 1,3693 b.Áp suất cuối quá trình cháy Pc
Pc = P a ε n 1 ¿0,085.16,50 1,3693 =3,9490 MPa c.Nhiệt độ cuối quá trình nén Tc
BI Quá trình cháy ” a.Lượng không khí lý thuyết cần để đốt cháy 1kg nhiên liệu Mo:
M0 0,211 ( 12 C+ H 4 − 32 O ) Đối với nhiên liệu của động cơ Diesel ta có:
M 0 =0,4357( kmol kk / kg nl ) b Lượng khí nạp mới M 1
32 4 d.Hệ số thay đổi phân tử lí thuyết 0 β 0 = M 2
=1,0898 M 1 0,6536 e.Hệ số thay đổi phân tử thực tế β β=1+ β
−1 =1+ 1,0898−1 =1,0872 1+ γ r 1+ 0,03 f Hệ số thay đổi phân tử thực tế tại điểm z( z ) β z =1+ β
1+γ r ξ b 1+ 0,03 0,85 g Tỷ nhiệt mol đẳng tích trung bình của môi chất tại điểm Z
2 vz =0,002965h Nhiệt độ cuối quá trình cháy T z
Theo bảng 1.14 sách tíí́nh toán nhiệt và động lực học động cơ đốt trong
Q H = 42530(Kj/kg) vì α > 1 nên ∆ Q H = 0 (Kj/kgnl)
T c 4,1019 giải phương trình trên và loại giá trị âm, ta có:
T z #18,72 K i Áp suất cuối quá trình cháy P z
IV Quá trình giãn nở
1 Tỷ số giản nở đầu β Z động cơ diesel ρ = ❑ b Tỷ số giãn nở sau δ= ε ρ = 1,5323 16,5
,7684 c Xác định chỉ số dãn nở đa biến trung bình (n¿¿2)¿ n 2 1 ¿ (0,9−0,8)
Thay dần các giá tri n 2 vào 2 vế của phương trình đến khi cân bằng, ta được n 2 = 1,1962
T b =¿ 1454 K c.Áp suất cuối quá trình giãn nở
P = P z δ n 2 b d.Kiểm nghiệm nhiệt độ khí sót
V Tính toán các thông số đặc trưng của chu trình a.Áp suất chỉ thị trung bình tính toán p i ' p a ¿ 0,0850.16,5−1 [1,7 (1,5323−1)+
1,0462 [ MN / m 2 ] b Áp suất chỉ thị trung bình thực tế
=1,0462.0,93=0,9730 [ MN / m 2 ] c Áp suất tổn hao cơ khí:
30 30 d Áp suất có ích trung bình
P e = P i −P m = 0,9730-0,1777 = 0,7953 [ MN / m 2 ] e Hiệu suất cơ giới η m =1− Pm
Q H P k η v ¿ 8,314.8 e Hiệu suất có ích η e =η m η i =0,8174.0,4737=0,3872 f Suất tiêu hao nhiên liệu chỉ thị g = i g Suất tiêu hao nhiên liệu có ích g = e
VI Tính thông số kết cấu của động cơ:
Tính thể tích công tác Vh:
Tíí́nh đường kíí́nh piston:
VII Vẽ đồ thị công chỉ thị
- Chọn tọa độ vuông góc:
Biễu diễn áp suất khíí́ thể (p) trên trục tung và thể tíí́ch khíí́ (V) trên trục hoành với tỉ lệ xíớ́ch à v và à p phự hợp với khổ giấy vẽ
- Xác định các điểm đặt của đồ thị công: Điêm a: cuôi qua trinh nap, ap suât P a , thê tích V a
V a = V h + V c = 1,2574 + 0,0811 = 1,3385 (dm 3 ) Điêm c: cuôi qua trinh nen
V c = 0,0811 (dm 3 ) Điêm z: cuôi qua trinh chay
(dm 3 ) Điêm b : điêm cuôi qua trinh giãn nơ
V b = V a = 1,3385 (dm 3 ) Điêm r : cuôi hanh trinh xa
Trong hành trình nén khíí́ trong xylanh bị nén với chỉ số đa biến trung bình n 1 =1.373, từ phương trình :
Pxn,Vxn là áp suất và thể tíí́ch khíí́ tại một điểm bất kì trên đướng cong nén Bằng cách cho các giá trị V xn đi từ V a đến V c ta lần lượt xác định được các giá trị p xn , bước nhảy của V xn là 20cm 3
+Dựng đường cong dãn nở
Trong quá trình giãã̃n nở khíí́ cháy được giãã̃n nở theo chỉ số giãã̃n nở đa biến n 2 , tương tự như trên, ta có: p z V zn 2 p xg V xgn 2 const p xg p z ( V z
Trong đó: p xg , V xg là áp suất và thể tíí́ch khíí́ tại một điểm bất kì trên đường cong giãã̃n nở.
Bằng cách cho các giá trị của V xg đi từ V z đến V b ta lần lượt xác định được xác định được các giá trị p xg bước nhảy của V xg là 20cm 3
Quá trình giãn nở
1 Tỷ số giản nở đầu β Z động cơ diesel ρ = ❑ b Tỷ số giãn nở sau δ= ε ρ = 1,5323 16,5
,7684 c Xác định chỉ số dãn nở đa biến trung bình (n¿¿2)¿ n 2 1 ¿ (0,9−0,8)
Thay dần các giá tri n 2 vào 2 vế của phương trình đến khi cân bằng, ta được n 2 = 1,1962
T b =¿ 1454 K c.Áp suất cuối quá trình giãn nở
P = P z δ n 2 b d.Kiểm nghiệm nhiệt độ khí sót
Tính toán các thông số đặc trưng của chu trình
a.Áp suất chỉ thị trung bình tính toán p i ' p a ¿ 0,0850.16,5−1 [1,7 (1,5323−1)+
1,0462 [ MN / m 2 ] b Áp suất chỉ thị trung bình thực tế
=1,0462.0,93=0,9730 [ MN / m 2 ] c Áp suất tổn hao cơ khí:
30 30 d Áp suất có ích trung bình
P e = P i −P m = 0,9730-0,1777 = 0,7953 [ MN / m 2 ] e Hiệu suất cơ giới η m =1− Pm
Q H P k η v ¿ 8,314.8 e Hiệu suất có ích η e =η m η i =0,8174.0,4737=0,3872 f Suất tiêu hao nhiên liệu chỉ thị g = i g Suất tiêu hao nhiên liệu có ích g = e
Tính thông số kết cấu của động cơ
Tính thể tích công tác Vh:
Tíí́nh đường kíí́nh piston:
Vẽ đồ thị công chỉ thị
- Chọn tọa độ vuông góc:
Biễu diễn áp suất khíí́ thể (p) trên trục tung và thể tíí́ch khíí́ (V) trên trục hoành với tỉ lệ xíớ́ch à v và à p phự hợp với khổ giấy vẽ
- Xác định các điểm đặt của đồ thị công: Điêm a: cuôi qua trinh nap, ap suât P a , thê tích V a
V a = V h + V c = 1,2574 + 0,0811 = 1,3385 (dm 3 ) Điêm c: cuôi qua trinh nen
V c = 0,0811 (dm 3 ) Điêm z: cuôi qua trinh chay
(dm 3 ) Điêm b : điêm cuôi qua trinh giãn nơ
V b = V a = 1,3385 (dm 3 ) Điêm r : cuôi hanh trinh xa
Trong hành trình nén khíí́ trong xylanh bị nén với chỉ số đa biến trung bình n 1 =1.373, từ phương trình :
Pxn,Vxn là áp suất và thể tíí́ch khíí́ tại một điểm bất kì trên đướng cong nén Bằng cách cho các giá trị V xn đi từ V a đến V c ta lần lượt xác định được các giá trị p xn , bước nhảy của V xn là 20cm 3
+Dựng đường cong dãn nở
Trong quá trình giãã̃n nở khíí́ cháy được giãã̃n nở theo chỉ số giãã̃n nở đa biến n 2 , tương tự như trên, ta có: p z V zn 2 p xg V xgn 2 const p xg p z ( V z
Trong đó: p xg , V xg là áp suất và thể tíí́ch khíí́ tại một điểm bất kì trên đường cong giãã̃n nở.
Bằng cách cho các giá trị của V xg đi từ V z đến V b ta lần lượt xác định được xác định được các giá trị p xg bước nhảy của V xg là 20cm 3
TÍNH TOÁN ĐỘNG LỰC HỌC CƠ CẤU TRỤC KHUỶU THANH TRUYỀN DỰNG VÀ HIỆU ĐÍNH ĐỒ THỊ CÔNG
Động học của piston
Khi trục khuỷu quay một góc α thì piston dịch chuyển được một khoảng S p so với vị tríí́ ban đầu (ĐCT) Chuyển vị của piston trong xylanh động cơ được tíí́nh bằng công thức:
2 :bán kíí́nh quay trục khuỷu (cm) b Tốc độ của piston:
Ta xác định được phương trình tốc độ của piston là hàm phụ thuộc theo góc quay của trục khuỷu bằng cách vi phân biểu thức (2.1) theo thời gian, ta được phương trình tốc độ của piston theo α:
(rad/s) c Gia tốc của piston:
Phương trình gia tốc của piston có được bằng cách lấy đạo hàm phương trình (2.2) theo thời gian, ta được:
Động lực học của cơ cấu khuỷu trục- thanh truyền
Lực khíí́ thể được xác định theo công thức sau:
Trong đó : p 0 ≈0,1 MN/m 2 –áp suất khíí́ quyển
Khối lượng chuyển động tịnh tiến của cơ cấu trục khuỷu thanh truyền:
Ta tíí́nh khối lượng này trên đơn vị diện tíí́ch đỉnh piston: m ' j m ' (0, 275 0, 35).m
Trong đó: m ' p , m ' tt ,m ' A lần lượt là khối lượng nhóm piston, tổng thanh truyền, đầu nhỏ thanh truyền trên một đơn vị diện tíí́ch.
Chọn m p’ =(g/cm 2 ); m tt’ =(g/cm 2 ) m A’ = (g/cm 2 ) m j’ = m p’ + m A’ = (g/cm 2 )= (kg/m 2 )
Lực quán tíí́nh trên đơn vị diện tíí́ch đỉnh piston:
B = 0.01180; phi 0.93; QH = 42530; fz 0.7; fb = 0.85; dentaT
= 23; lamda = 1.7; lamda1 = 1.02; lamda2 1; lamdaT = 1.11; ep 16.5;
Nv = (1/(ep-1))*(Tk/(Tk + dentaT))*(Pa/Pk)*(ep*lamda1 - lamdaT*lamda2*(Pr/Pa)^(1/m));
Yr = (lamda2*Pr*Tk)/((ep-1)*Nv*Pk*Tr);
% - Nhi?t ?? cu?i quá trình n?p Ta
Ta = (Tk + dentaT + lamdaT*Yr*Tr*(Pa/Pr)^(m-1/m))/(1+Yr);
%- Ch? s? nén ?a bi?n trung bình n1 syms x eqn = 8.314/(19.806+ (0.00419/2)*Ta*(ep^(x-1)+1)) - x == -1; solx vpasolve(eqn,x); n1 = solx; n1 = double(n1);
%- áp su?t quá trình nén Pc
%- Nhi?t ?? cu?i quá trình nén Tc
%- T? nhi?t mol ??ng tích trung bình
%- T? nhi?t mol ??ng tích trung bình c?a s?n ph?m cháy
MCv2phay = (19.867 + (1.634/alpha))+ 1/2*(427.38+(184.36/alpha))*10^(-5)*Tc;
%- T? nhi?t mol ??ng tích trung bình c?a h?n h?p khí trong quá trình nén
MCvc1phay = (MCv + Yr*MCv2phay)/(1+Yr);
%L??ng khí n?p m?i th?c t? n?p vào xilanh M1
%H? s? bi?n ??i phân t? khí lý thuy?t Beta0
%H? s? bi?n ??i phân t? khí th?c t? Beta
%H? s? bi?n ??i phân t? khí t?i ?i?m Bz xz = fz/fb;
%T?n th?t nhi?t l??ng do cháy không hoàn toàn dentaQH = 0;
%T? nhi?t mol ??ng tính trung bình c?a môi ch?t t?i ?i?m Z
MCvz2phay = (M2*(xz + Yr/Beta0)*MCv2phay+ M1*(1-xz)*MCv)/(M2*(xz + Yr/Beta0)+ M1*(1-xz));
% Nhi?t ?? cu?i quá trình cháy Tz
Tz = ((fz*QH)/(M1*(1 + Yr)) + (MCvc1phay + 8.314*lamda)*Tc)/(BetaZ*(MCvz2phay + 8.314));
%áp su?t cu?i quá trình cháy Pz
%TÍNH TOÁN QUÁ TRÌNH DÃN N?
%T? s? dãn n? ban ??u p = (BetaZ*Tz)/(lamda*Tc);
%Ch? s? dãn n? ?a bi?n trung bình n2, nhi?t ?? cu?i quá trính dãn n? Tb syms Tb n2 sol= vpasolve([Tb == Tz/(delta^(n2 - 1)) , n2 - 1 == 8.314/(((fb-fz)*QH)/(M1*(1 + Yr)*Beta*(Tz - Tb)) + 21.101 + (0.00593/2)*(Tz + Tb))],[Tb,n2]);
%áp su?t cu?i quá trình dãn n? Pb
%Ki?m nghi?m nhi?t ?? khí sót Trtt
%TINH TOÁN THÔNG S? ??C TR?NG CHU K?
%áp su?t ch? th? trung bình tính toán Pi1phay
Pi1phay = Pa*(ep^n1)/(ep - 1)*((lamda*(p-1)) + (lamda*p/(n2 - 1))*(1 -
%áp su?t ch? th? trung bình th?c t? Pi
%Xác ??nh áp su?t có ích trung bình Pe
Ni = 8.314*(M1*Pi*Tk)/(QH*Pk*Nv);
Ne = 8.314*(M1*Pe*Tk)/(QH*Nv*Pk);
%Tiêu hao nhiên li?u ch? th?
%Tiêu hao nhiên li?u có ích
%TINH THONG SO KET CAU CUA ??NG C?
%Th? tích công tác Vh Vh
%Ve tung diem dac biet x = [x0 x1 x2 x3 x4]; y
%DIEM XAC DINH TREN DO THI plot(Vc,Pz,'b.')%Z' plot((Vc+x2)/2,Pz,'b.')%Z'' plot(x1,(Pz-y1)/3+y1,'b.')%C'' plot(x0,y0+(y3-y0)/2,'b.')%B'' lc= Pa*(Va.^n1./0.1126.^n1); plot(0.1126,lc,'b.')%C' lb=(Pz)*(Vz.^n2./1.1583.^n2); plot(1.1583,lb,'b.')%B'
% DUNG DUONG CONG NEN phu1 = (Va - 0.09012)/0.001;
Pxn = Pa*(Va.^n1./Vxn.^n1); plot(Vxn,Pxn,'r-');
% DUNG DUONG DAN NO phu2= (1.1583 -0.1401)/0.001;
Pxg = (Pz)*(Vz.^n2./Vxg.^n2); plot(Vxg,Pxg,'r-'); ylim([-1 7]) xlim([0 1.5]) axis fill;
%%% doan 1 x= 0.08112:0.001:0.09013; y= ( 3.42e-09)*x.^-8.13+2.345; plot (x,y,'r-'); plot ([x1, x1],[ (Pz-y1)/3+y1 , 5.958],'r-'); x= 0.08612:0.001:0.1403; y= 408.2*x.*exp(-9.958*x) -8.364; plot (x,y,'r-'); x= 0.08112:0.001:0.08612; y= (-1.716e-15)*x.^-13.56+7.022 plot (x,y,'r-');
%VE DUONG NUA DUONG TRON BEN DUOI
Ophay = (lamdaphay*(Vh/2))/2; xCenter_2 = (Va+Vc)/2; yCenter_2 0; radius = (Va-Vc)/2 ; plot(xCenter_2,0,'.'); plot(xCenter_2 + Ophay,0,'r.'); fxy =@(x,y) (x - xCenter_2).^2 + (y - yCenter_2).^2 - radius.^2 +(y>=0) ;
% DIEM XAC DINH TREN DUONG TRON plot(0.1126,-
0.1913,'b.'); % nay la thetaS plot(0.1686,-0.3209,'b.'); % nay la theta1& plot(0.1451,-0.2753,'b.'); % nay la theta4" plot(1.2173,-0.3718,'b.'); % nay la theta2H plot(1.1583,-0.4406,'b.'); % nay la theta3V hold off;
%%%%%%%%%%%% title('Do thi cong P - V') ylabel('P'); xlabel('V'); close all; clc;
%TINH TOAN DONG LUC HOC CO CAU TRUC KHUYU THANH TRUYEN
R= S/2*10^-1; goc_alpha = 0 : pi/30 : 2*pi; goc_beta= asin(lamdaphay*sin(goc_alpha));
Sp = R*((1- cos(goc_alpha))+ (lamdaphay/4)*(1 - cos(2*goc_alpha))); subplot ( 3,1,1); plot(goc_alpha*180/pi,Sp,'r-'); title('Chuyen vi cua piston') legend('Sp','Sp1','Sp2') ; ylabel('Sp'); xlabel('alpha'); hold on omega = pi*n/30;
Vp = R*omega*(sin(goc_alpha)+ lamdaphay/2*sin(2*goc_alpha));
Vp1 = R*omega*sin(goc_alpha);
Vp2 = R*omega*(lamdaphay/2)*sin(2*goc_alpha); subplot(3,1,2); plot(goc_alpha*180/pi, Vp, 'r-'); hold on plot(goc_alpha*180/pi, Vp1, 'b-'); plot(goc_alpha*180/pi, Vp2, 'g-'); title('Van toc cua piston') legend('Vp','Vp1','Vp2') ; ylabel('Vp'); xlabel('alpha');
Jp = R*omega^2*(cos(goc_alpha) + lamdaphay*cos(2*goc_alpha));
Jp1 = R*omega^2*cos(goc_alpha);
Jp2 = Jp - Jp1; subplot(3,1,3); plot(goc_alpha*180/pi, Jp , 'r-'); hold on plot(goc_alpha*180/pi, Jp1 , 'b-'); plot(goc_alpha*180/pi, Jp2 , 'g-'); title('Gia toc cua piston') legend('Jp','Jp1','Jp2') ; ylabel('Jp'); xlabel('alpha'); maphay = 0.35*mttphay; mbphay = 0.65*mttphay; mjphay = mpphay + maphay; mrphay = mkphay + mbphay; mj = mjphay*10.^-3*pi*(D*10)^2/4; mr = mrphay*pi*D^2/4; hold on; ylim([-7 7]) xlim([-4 4]) axis fill; i=0:1:180;
Pj = -mjphay*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10.^-5; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi;
Z=Ptt.*cos((i+goc_beta).*pi/180)./cos(goc_beta*pi/180);
T=Ptt.*sin((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(T,Z,'r-'); i0:1:345; j=1:0.0934:16.5;
Pj = -mjphay*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10.^-5; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi;
Z=Ptt.*cos((i+goc_beta).*pi/180)./cos(goc_beta*pi/180);
T=Ptt.*sin((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(T,Z,'r-'); i45:2:361;
Pj = -mjphay*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10.^-5; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi;
Z=Ptt.*cos((i+goc_beta).*pi/180)./cos(goc_beta*pi/180);
T=Ptt.*sin((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(T,Z,'r-'); i60:1:518; length(i); j=1:0.0618:delta; length(j);
Pj = -mjphay*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10.^-5; Ptt= Pkt+Pj;
21 goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi;
Z=Ptt.*cos((i+goc_beta).*pi/180)./cos(goc_beta*pi/180);
T=Ptt.*sin((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(T,Z,'r-'); i= 518:1:530;
Pj = -mjphay*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10.^-5; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi;
T=Ptt.*sin((i+goc_beta).*pi/180)./cos(goc_beta*pi/180);
Z=Ptt.*cos((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(T,Z,'r-'); i= 530:1:570;
Pj = -mjphay*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10.^-5; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi;
Z=Ptt.*cos((i+goc_beta).*pi/180)./cos(goc_beta*pi/180);
T=Ptt.*sin((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(T,Z,'r-'); iW0:1:720;
Pj = -mjphay*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10.^-5; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi;
Z=Ptt.*cos((i+goc_beta).*pi/180)./cos(goc_beta*pi/180);
T=Ptt.*sin((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(T,Z,'r-'); title('Do thi phu tai tac dung len chot khuyu') ylabel('- Z'); xlabel('alpha"'); hold on; ylim([-3 7]) xlim([0 800]) axis fill; i=0:1:180;
Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi;
Z=Ptt.*cos((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(i,Z,'r-'); i0:1:345;
Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi;
Z=Ptt.*cos((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(i,Z,'r-'); i45:2:361;
Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi;
Z=Ptt.*cos((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(i,Z,'r-'); i60:1:518; length(i); j=1:0.0618:delta; length(j);
Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi;
Z=Ptt.*cos((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(i,Z,'r-'); i= 518:1:660;
Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi;
Z=Ptt.*cos((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(i,Z,'r-'); if0:1:720;
Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi;
Z=Ptt.*cos((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(i,Z,'r-'); title('Luc phap tuyen gay uon truc khuyu') ylabel('Z (MN)'); xlabel('GQTK'); hold on; ylim([-3 7]) xlim([0 800]) axis fill; i=0:1:180;
Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj;
23 goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi;
T=Ptt.*sin((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(i,T,'r-'); i0:1:345;
Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi;
T=Ptt.*sin((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(i,T,'r-'); i45:2:361;
Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi;
T=Ptt.*sin((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(i,T,'r-'); i60:1:518; length(i); j=1:0.0618:delta; length(j);
Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi;
T=Ptt.*sin((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(i,T,'r-'); i= 518:1:660;
Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi;
T=Ptt.*sin((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(i,T,'r-'); if0:1:720;
Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi;
T=Ptt.*sin((i+goc_beta).*pi/180)./cos(goc_beta*pi/180); plot(i,T,'r-'); title('Luc tiep tuyen lam quay truc khuyu') ylabel('T (MN)'); xlabel('GQTK');
%plot(359,Pz-0.1,'r.'); hold on; ylim([-3 7]) xlim([0 800]) axis fill; i=0:15:180;
Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; plot(i,Pj,'r-'); plot(i,Ptt,'b-'); plot([0 180],[ (Pa - 0.1) (Pa - 0.1)],'g-'); i0:1:345;
Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; plot(i,Pj,'r-'); plot(i,Ptt,'b-'); plot(i,Pkt,'g-'); i45:2:361;
Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; plot(i,Pj,'r-'); plot(i,Ptt,'b-'); plot(i,Pkt,'g-'); i60:1:518; length(i); j=1:0.0618:delta; length(j);
Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; plot(i,Pj,'r-'); plot(i,Ptt,'b-'); plot(i,Pkt,'g-');
Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; plot(i,Pj,'r-'); plot(i,Ptt,'b-'); plot(i,Pkt,'g-'); if0:1:720;
Pj = -mjphay*10*R*omega.^2*(cos(i*pi/180)+lamdaphay*cos(2*i*pi/180))*10^-6; Ptt= Pkt+Pj; goc_beta= asin(lamdaphay*sin(i*pi/180))*180/pi; plot(i,Pj,'r-'); plot(i,Ptt,'b-'); plot([660 720],[ Pr-0.1 Pr-0.1],'g-'); title('Do thi luc khi the Pkt, luc quan tin Pj, luc tong hop') legend('Pj','Ptt','Pkt') ; ylabel('P'); xlabel('alpha');
