Mô phỏng các bước tính toán chi tiết cơ khí bằng phần mềm Matlab

disp' Tinh toan kiem nghiem ben'disp' Tinh suc ben dau nho thanh truyen' disp' Tinh suc ben dau nho khi chiu keo' disp'coi dau nho la mot dam cong duoc ngam 2 dau tiet dien I-i' disp'ung

Trang 1

BẢNG PHỤ LỤC BẢNG PHỤ LỤC 1

A Mô phỏng các bước tính toán bằng phần mềm

1 Phần tính toán Thanh truyền sử dụng Matlap.

%%CHUONG TRINH MATLAP TINH TOAN THANH TRUYEN

%%CAC BUOC THUC HIEN

disp('Nhap thong so cho truoc')

n=input('n=')

Dxl=input('Duong kinh xylanh (m): Dxl=')

S=input('Hanh trinh piston (m): S=')

epsilon=input('epsilon=')

l=input('Chieu dai thanh truyen (m): l=')

mtt=input('Khoi luong thanh truyen (kg): mtt=')

Mnp=input('Khoi luong nhom piston (kg): Mnp=')

pa=input('pa=')

disp('Nhap thong so tinh toan')

Dpt=Dxl

dcp=input('duong kinh chot dcp=')

db=input('duong kinh be chot db=')

d0=input('duong kinh lo tren chot d0=')

denta=input('chieu dai bac lot denta=')

dentaphay=input('Khe ho huong kinh bac lot va chot piston dentaphay=')ld1=input('ld1=')

elseif denta<0.08*dcp | denta>0.085*dcp

disp('chieu dai bac lot sai roi')

return

elseif dentaphay<0.0001*dcp | dentaphay>0.00015*dcp

disp('Khe ho huong kinh bac lot va chot piston sai roi')

Trang 2

disp(' Tinh toan kiem nghiem ben')

disp(' Tinh suc ben dau nho thanh truyen')

disp(' Tinh suc ben dau nho khi chiu keo')

disp('coi dau nho la mot dam cong duoc ngam 2 dau (tiet dien I-i)')

disp('ung voi goc gama la')

ro1=input('ban kinh goc luon noi dau nho voi than ro1=')

H=2*r1

gama=90+acos((H/2+ro1)/(r2+ro1))

disp('do tinh chat doi xung cua ngam nen khi tinh toan, ta cat bo mot nua

va thay the bang cac luc phap tuyen va momen uon NA, MA')

disp('dam cong bao gom hai cung ung voi goc gamax nho hon hoac bang 90 do

va lon hon 90 do')

Trang 3

disp('Tren mat trong')

otj=((-2*Mj*(6*ro-s))/(s*(2*ro-s))+xima*Nj)*1/(ld1*s)

disp('b, tinh suc ben dau nho khi chiu nen')

disp('Luc nen tac dung len dau nho thanh truyen la hop luc cua luc khi the

va luc quan tinh cua khoi luong piston')

disp('Momen uon va luc phap tuyen tren cung BC ung voi gamax >= 90')

gamax=input('Nhap gamax tinh theo radian, gamax=')

cos(gamax)/pi)

Mz2=MA+NA*ro*(1-cos(gamax))-Pxichma*ro*(sin(gamax)/2-gamax*sin(gamax)/pi-Nz2=NA*cos(gamax)-Pxichma*ro*(sin(gamax)/2-gamax*sin(gamax)/pi-cos(gamax)/pi)

disp('Nhu da phan tich o tren, do lap ghep cang bac lot tren dau nho, nen luc phap tuyen tac dung tren dau nho ma khong phai toan bo Nz chi la mot phan cua Nz tuc la xima*Nz')

disp('Ung suat tong gay ra trong dau nho khi chiu nen la')

disp('Tren mat ngoai')

disp('Ban nhap sai roi')

t=input('Nhap t o khoang(100 den 140) t=')

P=(denta+dentat)/(d1*(((d2^2+d1^2)/(d2^2-d1^2)+nuy)/En)+((d2^2+db^2)/(d2^2-disp('ung suat bien dang theo cong thuc Lame')

disp('ung suat ben ngoai mat dau nho')

ondenta=P*(2*d1^2)/(d2^2-d1^2)

disp('ung suat tren mat trong')

otdenta=P*(d2^2+d1^2)/(d2^2-d1^2)

disp('d he so an toan cua dau nho thanh truyen')

disp('tinh toan cho tiet dien nguy hiem va tren mat ngoai')

onz=onz2*1

omax=onj+onz

omin=ondenta+onz

Trang 4

disp(' Tinh suc ben than thanh truyen')

disp('a Tinh tiet dien nho nhat(tiet dien I-I)')

disp('b Tinh o tiet dien trung binh (tiet dien II-II)')

disp('ung suat keo')

Trang 5

dck=input('duong kinh chot khuyu dck=')

tblm=input('chieu day bac lot mong tblm=')

tbld=input('chieu day bac lot day tbld=')

C=input('Khoang cach tam bulong C=')

Trang 6

disp('luc xiet bulong khi lap ghep')

disp('chon he so kinh nghiem a nam trong khoang 2-4 a=')

a=input('he so kinh nghiem a=')

while a<2 | a>4

disp('ban nhap sai roi')

a=input('a=')

end

Px=a*Pb

disp('luc tac dung len bulong thanh truyen')

disp('he so bien dang thanh truyen nam trong khoang 0,15-0,3')

x1=input('chon he so bien dang x1=')

while x1<0.15 | x1>0.3

disp('ban chon sai roi')

x1=input('he so bien dang x1=')

Trang 7

2 Phần tính toán Trục Khuỷu sử dụng Matlap.

%%CHUONG TRINH MATLAP TINH TOAN TRUC KHUYU

%%CAC BUOC THUC HIEN

disp('chuong trinh tinh toan truc khuyu cua loai dong co xang khong tang

ap kieu dong co mot hang')

disp('Nhap thong so dau vao')

Nemax=input('cong suat dong co Nemax=')

n=input('toc do dong co n=')

Dxl=input('Duong kinh xilanh (don vi la m) Dxl=')

s=input('Hanh trinh piston (don vi la m) s=')

l=input('Chieu dai thanh truyen (don vi la m)l=')

mnp=input('Khoi luong nhom piston (don vi la kg) mnp=')

mtt=input('Khoi luong nhom thanh truyen (don vi la kg) mtt=')

pkt=input('Ap suat khi the (don vi la MN) pkt=')

disp('Nhap thong so chon')

dct=input('Duong kinh co truc khuyu (m) dct=')

lct=input('Chieu dai co truc khuyu(m) lct=')

dch=input('Duong kinh chot khuyu(m) dch=')

lch=input('Chieu dai chot khuyu(m) lch=')

h=input('chieu rong ma khuyu(m) h=')

b=input('do day ma khuyu(m) b=')

r=input('Ban kinh goc luon(m) r=')

disp('do trung diep la')

disp('Do so vong quay cua dong co khi khoi dong nho nen ta co the bo qua

Trang 8

anh huong cuar luc quan tinh')

disp('Do vay luc nay luc tac dung chi con lai luc do ap suat khi ther lonnhat trong xilanh Pzmax')

disp('khoang cach tu tam ma khuyu den tam chot khuyu')

disp('a, tinh suc ben cua chot khuyu')

disp('Momen uon chot khuyu la')

disp('b, tinh suc ben cho ma khuyu')

disp('ung suat uon ma khuyu la')

disp('c, tinh suc ben cho co truc khuyu')

disp('ung suat uon co truc khuyu')

ou=Zphay*bphay/Wu

disp('truong hop truc khuyu chiu luc Zmax')

m1=input('Khoi luong thanh truyen quy ve dau nho m1=')

disp('Tinh toan kiem nghiem ben')

disp('Luc tac dung len truc khuyu luc nay la Zmax va vi tri truc khuyu la

Trang 9

disp('Luc quan tinh ly tam cua chot khuyu')

disp('Ta co van toc goc cua truc khuyu w=')

disp('+ tinh suc ben cua chot khuyu')

disp('ung suat uon cua chot khuyu')

disp('Tinh suc ben cua co truc khuyu')

disp('Ung suat uon cua co truc khuyu')

disp('Tinh suc ben cua ma khuyu')

disp('Trong qua trinh lam viec ma khuyu chiu nen va uon theo truc x-x')disp('Ung suat nen cua ma khuyu')

3 Phần tính toán Phanh sử dụng Matlap.

%TINH TOAN THIET KE CO CAU PHANH XE Ma_500

Trang 10

%CAC BUOC THUC HIEN

disp('Nhap cac thong so co ban')

G=input('nhap trong luong toan bo cua xe khi day tai(N): G= ');

G_1=input('nhap trong luong toan bo cua xe phan ra cau truoc(N): G_1= ');G_2=input('nhap trong luong toan bo cua xe phan ra cau sau(N): G_2= ');g=input('nhap gia tri gia toc trong truong (m/s^2):g= ');

lamda=input('nhap he so bien dang cua lop :lamda= ');

L=input('nhap chieu dai co so cua xe (m): L= ');

hg=input('nhap gia tri do cao trong tam cua xe(m): hg= ');

j_max=input('nhap gia toc cham dan cuc dai cua oto khi phanh(m/s^2): j_max= ');

phi=input('nhap gia tri he so bam cua lop voi mat duong phi= ');

B_bx=input('nhap be rong cua lop banh xe(m):B_bx= ');

d_bx=input('nhap duong kinh cua vanh banh xe(m):d_bx= ');

r_t=input('nhap ban kinh tang trong(m):r_t= ');

disp('1.Xac dinh mo men sinh ra o cac co cau phanh ')

disp(' Ban kinh banh xe theo thiet ke la (m)')

disp('2 Xac dinh diem dat luc')

beta_0t=input('nhap goc om tam ma sat cua co cau phanh

Trang 11

disp('3 Xac dinh luc can thiet tac dung len guoc phanh')

disp('Do su dung phuong phap hoa do luc Nen ta dua vao ket qua tinh toan o ngoai ta duoc ket qua ')

% Luc N1 la luc huong tam ,luc T1 la luc tiep tuyen , luc R1 la luc tong

% hop cua N1 va T1.Luc R1 tao voi N1 goc phi1

% cac luc tac dung len guoc phanh la: luc P do cam quay sinh ra, luc R phan

% luc cua trong tac dung len ma phanh , luc U phan luc tu chot phanh tac

% dung len guoc phanh

disp('co cau phanh truoc ');

R_1t=input('nhap gia tri R_1t=');

P_1t=input('nhap gia tri P_1t=');

U_1t=input('nhap gia tri U_1t=');

disp('co cau phanh sau ');

R_1s=input('nhap gia tri R_1s=');

P_1s=input('nhap gia tri P_1s=');

U_1s=input('nhap gia tri U_1s=');

disp('4 xac dinh kich thuoc guoc phanh ')

% a guoc phanh

V_0=input('nhap gia tri van toc cua oto khi bat dau phanh(km/h): V_0=');

%dien tich toan bo cua cac ma phanh o tat ca cac co cau phanh oto

m=input('nhap so luong ma phanh m= ');

b_mp=input('nhap chieu rong ma phanh(cm)= ');

xicma_tongF=2*(beta_0t*r_t*b_mp)+2*(beta_0s*r_t*b_mp)

('cong ma sat rieng la (J/cm2)')

L=(G*V_0^2)/(2*g*xicma_tongF)

N_1=('khong thoa man cong cho phep nhap them so luong ma phanh ');

D_l=('thoa man cong ma sat rieng cho phep ');

disp(' 5 Ap suat tren be mat ma sat ')

disp('ap suat tren ma phanh la (MN/m2):')

q=M_pt*180/(muy*b_mp*r_t^2*beta_0s*pi)

N_q=('khong thoa man ap suat tren be mat ma phanh ');

D_q=('thoa man ap suat tren be mat ma phanh ');

disp('6 Tinh toan nhiet phat ra khi phanh')

disp(' do day tang trong va de tang trong(cm)')

denta=1.2

denta_d=1.5

b_tt=input('nhap chieu rong tang trong(cm) = ');

v1=input('nhap toc do ban dau khi phanh(m/s)=');

v2=input('nhap toc do cuoi khi phanh(m/s)=');

disp('the tich vat lieu lam tang trong (cm^3)')

V_t=2*pi*(r_t+0.6)*b_tt*denta+pi*denta_d

disp(' khoi luong rieng cua vat lieu lam tang trong(g/cm^3)')

gama=7

Trang 12

disp(' nhiet dung cua vat lieu lam tang trong (gang)la(J/(kg*do))')

c=500

G_bt=input('Nhap trong luong ban than cua xe = ')

disp('khoi luong cac trong phanh m_t(kg)')

m_t=4*gama*V_t/1000

disp('su tang nhiet do trong phanh la:')

t_0=G_bt*((v1/3.6)^2-(v2/3.6)^2)/(2*g*m_t*c)

N_t0=('khong thoa man nhiet do khi phanh');

D_t0=('thoa man nhiet do khi phanh');

disp('7 Tinh ben trong phanh ')

disp(' ban kinh trong va ngoai cua trong phanh(cm)')

N_x=('khong thoa ung suat cua trong phanh ');

D_x=('thoa man ung suat cua trong phanh ');

disp('8 tinh ben chot phanh ')

disp('ung suat cat')

d_c=input('Nhap duong kinh chot phanh(cm)= ');

l=input('Nhap chieu dai tiep xuc cua chot voi guoc phanh=');

disp(' 9 tinh ben loxo hoi vi')

disp(' so vong loxo lam viec (vong)la chon:')

n_0=20

disp(' mo dun dan hoi dich chuyen cua loxo (MN/m2)la ')

G_lx=8*10^4

d_lx=input('nhap duong kinh day loxo (mm)d_lx=');

disp('duong kinh vong loxo (mm) chon ')

D=4*d_lx

disp(' do cung cua loxo la :')

C=(d_lx^4*G_lx)/(8*D^3*n_0)

Trang 13

disp('luc tac dung len lo xo ')

N_lx=('khong thoa ung suat cua loxo Can tang duong kinh loxo ');

D_lx=('thoa man ung suat cua loxo ');

disp('10 tinh ben guoc phanh')

a_gp=input('Nhap chieu dai phan tren chu T(m):=');

b_gp=input('Nhap chieu rong phan tren chu T(m):=');

c_gp=input('Nhap chieu rong phan duoi chu T(m):=');

d_gp=input('Nhap chieu dai phan duoi chu (m): =');

disp(' kich thuoc che tao guoc phanh (m): ')

R1=input('ban kinh tam cua banh xe den dinh chu T(m)= ');

disp('ban kinh trong tam cua phan dien tich tren ,tinh den tam tang trong')

R_1=R1-b_gp/2

R2=input('ban kinh tam cua banh xe den day phan tren chu T (m)= ');

disp('ban kinh trong tam cua phan dien tich duoi ,tinh den tam tang trong')

disp('11 kiem tra ben guoc phanh')

% ta dat cac luc nay vao guoc truoc cua co cau phanh sau , o tai diemdat

% luc tong hop R1 , ta phan ra hai thanh phan luc huong kinh N1 , va luc

% tiep tuyen T1 , o tai chot quay cua guoc phanh tacung phan luc tong hop

% U1 ra hai thanh phan U1x va U1y sau do guoc phanh nay o tai luc R1 ta

% cat huong kinh U1y va luc U1x ,sau do guoc phanh nay o tai luc R1 ta cat

% ra thanh hai thayvao mat cat do luc huong tam Nz1 va Qy1 va Mu1 O nua

% duoi la Nz2, Qy2 va Mu2 nguoc voi cac thanh phan luc va momen cua phan

Trang 14

disp('Xet nua tren guoc phanh ')

disp('dau am chi chieu cua luc nguoc chieu gia thiet')

disp('Xet su can bang cho nua duoi guoc phanh ')

disp('xet tai diem B')

%dua vao xe tham khao ta chon

disp('xet su can bang tai diem C :beta_0=0')

disp('F la dien tich cua tiet dien tinh toan')

disp('Rth la ban kinh duong trung hoa')

disp('Ri la ban kinh tai diem dang xet')

disp('Ung suat tai diem A co ')

disp('Sx mo men tinh phan bi cat doi voi truc quan tinh trung tam')

disp('Jx mo men quan tinh tiet dien dang xet')

disp('Nx luc cat ')

disp('b chieu day phan bi cat')

disp('Ung suat tai diem A ')

Trang 15

N_th=('khong thoa ung suat tong hop Can tinh lai ');

D_th=('thoa man ung suat tong hop ');

%TINH TOAN THIET KE CO CAU PHANH XE Ma_500

%CAC BUOC THUC HIEN

disp('Nhap cac thong so co ban')

G=input('nhap trong luong toan bo cua xe khi day tai(N): G= ');

G_1=input('nhap trong luong toan bo cua xe phan ra cau truoc(N): G_1= ');G_2=input('nhap trong luong toan bo cua xe phan ra cau sau(N): G_2= ');g=input('nhap gia tri gia toc trong truong (m/s^2):g= ');

lamda=input('nhap he so bien dang cua lop :lamda= ');

L=input('nhap chieu dai co so cua xe (m): L= ');

hg=input('nhap gia tri do cao trong tam cua xe(m): hg= ');

j_max=input('nhap gia toc cham dan cuc dai cua oto khi phanh(m/s^2): j_max= ');

phi=input('nhap gia tri he so bam cua lop voi mat duong phi= ');

B_bx=input('nhap be rong cua lop banh xe(m):B_bx= ');

d_bx=input('nhap duong kinh cua vanh banh xe(m):d_bx= ');

r_t=input('nhap ban kinh tang trong(m):r_t= ');

disp('1.Xac dinh mo men sinh ra o cac co cau phanh ')

disp(' Ban kinh banh xe theo thiet ke la (m)')

Trang 16

disp('2 Xac dinh diem dat luc')

beta_0t=input('nhap goc om tam ma sat cua co cau phanh

disp('3 Xac dinh luc can thiet tac dung len guoc phanh')

disp('Do su dung phuong phap hoa do luc Nen ta dua vao ket qua tinh toan o ngoai ta duoc ket qua ')

% Luc N1 la luc huong tam ,luc T1 la luc tiep tuyen , luc R1 la luc tong

% hop cua N1 va T1.Luc R1 tao voi N1 goc phi1

% cac luc tac dung len guoc phanh la: luc P do cam quay sinh ra, luc R phan

% luc cua trong tac dung len ma phanh , luc U phan luc tu chot phanh tac

% dung len guoc phanh

disp('co cau phanh truoc ');

disp('co cau phanh sau ');

disp('4 xac dinh kich thuoc guoc phanh ')

% a guoc phanh

V_0=input('nhap gia tri van toc cua oto khi bat dau phanh(km/h): V_0=');

Trang 17

%dien tich toan bo cua cac ma phanh o tat ca cac co cau phanh oto

m=input('nhap so luong ma phanh m= ');

b_mp=input('nhap chieu rong ma phanh(cm)= ');

xicma_tongF=2*(beta_0t*r_t*b_mp)+2*(beta_0s*r_t*b_mp)

('cong ma sat rieng la (J/cm2)')

L=(G*V_0^2)/(2*g*xicma_tongF)

N_1=('khong thoa man cong cho phep nhap them so luong ma phanh ');

D_l=('thoa man cong ma sat rieng cho phep ');

disp(' 5 Ap suat tren be mat ma sat ')

disp('ap suat tren ma phanh la (MN/m2):')

q=M_pt*180/(muy*b_mp*r_t^2*beta_0s*pi)

N_q=('khong thoa man ap suat tren be mat ma phanh ');

D_q=('thoa man ap suat tren be mat ma phanh ');

disp('6 Tinh toan nhiet phat ra khi phanh')

disp(' do day tang trong va de tang trong(cm)')

denta=1.2

denta_d=1.5

b_tt=input('nhap chieu rong tang trong(cm) = ');

v1=input('nhap toc do ban dau khi phanh(m/s)=');

v2=input('nhap toc do cuoi khi phanh(m/s)=');

disp('the tich vat lieu lam tang trong (cm^3)')

G_bt=input('Nhap trong luong ban than cua xe = ')

disp('khoi luong cac trong phanh m_t(kg)')

m_t=4*gama*V_t/1000

disp('su tang nhiet do trong phanh la:')

t_0=G_bt*((v1/3.6)^2-(v2/3.6)^2)/(2*g*m_t*c)

N_t0=('khong thoa man nhiet do khi phanh');

D_t0=('thoa man nhiet do khi phanh');

disp('7 Tinh ben trong phanh ')

disp(' ban kinh trong va ngoai cua trong phanh(cm)')

Trang 18

disp('ung suat tong hop tac dung len trong phanh co them he so an toan 1.5 (kg/cm2)')

xicma_tong=1.5*sqrt(xicma_n^2+xicma_t^2)

N_x=('khong thoa ung suat cua trong phanh ');

D_x=('thoa man ung suat cua trong phanh ');

disp('8 tinh ben chot phanh ')

disp('ung suat cat')

d_c=input('Nhap duong kinh chot phanh(cm)= ');

l=input('Nhap chieu dai tiep xuc cua chot voi guoc phanh=');

disp(' 9 tinh ben loxo hoi vi')

disp(' so vong loxo lam viec (vong)la chon:')

n_0=20

disp(' mo dun dan hoi dich chuyen cua loxo (MN/m2)la ')

G_lx=8*10^4

d_lx=input('nhap duong kinh day loxo (mm)d_lx=');

disp('duong kinh vong loxo (mm) chon ')

N_lx=('khong thoa ung suat cua loxo Can tang duong kinh loxo ');

D_lx=('thoa man ung suat cua loxo ');

disp('10 tinh ben guoc phanh')

a_gp=input('Nhap chieu dai phan tren chu T(m):=');

b_gp=input('Nhap chieu rong phan tren chu T(m):=');

c_gp=input('Nhap chieu rong phan duoi chu T(m):=');

d_gp=input('Nhap chieu dai phan duoi chu (m): =');

disp(' kich thuoc che tao guoc phanh (m): ')

Trang 19

R1=input('ban kinh tam cua banh xe den dinh chu T(m)= ');

disp('ban kinh trong tam cua phan dien tich tren ,tinh den tam tang trong')

R_1=R1-b_gp/2

R2=input('ban kinh tam cua banh xe den day phan tren chu T (m)= ');

disp('ban kinh trong tam cua phan dien tich duoi ,tinh den tam tang trong')

disp('11 kiem tra ben guoc phanh')

% ta dat cac luc nay vao guoc truoc cua co cau phanh sau , o tai diemdat

% luc tong hop R1 , ta phan ra hai thanh phan luc huong kinh N1 , va luc

% tiep tuyen T1 , o tai chot quay cua guoc phanh tacung phan luc tong hop

% U1 ra hai thanh phan U1x va U1y sau do guoc phanh nay o tai luc R1 ta

% cat huong kinh U1y va luc U1x ,sau do guoc phanh nay o tai luc R1 ta cat

% ra thanh hai thayvao mat cat do luc huong tam Nz1 va Qy1 va Mu1 O nua

% duoi la Nz2, Qy2 va Mu2 nguoc voi cac thanh phan luc va momen cua phan

disp('dau am chi chieu cua luc nguoc chieu gia thiet')

disp('Xet su can bang cho nua duoi guoc phanh ')

disp('xet tai diem B')

disp('xet su can bang tai diem C :beta_0=0')

beta=0

Nz2=-U1x*cos(beta)-U1y*sin(beta)

Qy2=U1y*sin(beta)-U1x*cos(beta)

Trang 20

disp('Ung suat do Qy2 va Mu1 gay ra duoc tinh theo cong thuc :xicma=Qy2/F+ Mu1/F*(1-R_th/Ri)')

disp('F la dien tich cua tiet dien tinh toan')

disp('Rth la ban kinh duong trung hoa')

disp('Ri la ban kinh tai diem dang xet')

disp('Ung suat tai diem A co ')

disp('Sx mo men tinh phan bi cat doi voi truc quan tinh trung tam')

disp('Jx mo men quan tinh tiet dien dang xet')

disp('Nx luc cat ')

disp('b chieu day phan bi cat')

disp('Ung suat tai diem A ')

N_th=('khong thoa ung suat tong hop Can tinh lai ');

D_th=('thoa man ung suat tong hop ');

if xicma_tha>=4*10^7|xicma_thb>=4*10^7|xicma_thc>=4*10^7

N_th

Trang 21

%%CHUONG TRINH MATLAB TINH TOAN THIET KE HOP SO CO KHI

%%CAC BUOC TINH TOAN HOP SO CO KHI

disp('chon ty so truyen')

Memax=input('Memax=')

rota=45

while rota <40 | rota >50

rota=input('hay phimax khoang tren,rota=')

disp('1.che do tai trong thiet ke')

disp('a tai trong tinh tu dong co den chi tiet dang tinh cua hop so')

disp('ap dung cong thuc Mt=Me.i')

disp('trong do')

disp('+Mt momen tinh toan chi tiet can tinh')

disp('+Me momen cua dong co')

disp('Me=Memax')

disp('+i ty so truyen tu dong co den chi tiet can tinh toan')Mt=Memax*ih1

phimax=0.8

while phimax <0.7 | phimax >0.9

phimax=input('hay phimax khoang tren,phimax=')

end

disp('ty so truyen tu banh xe den hop so')

Mphimax=(phimax*Gphi*rbx)/i0*If*icc

disp('2.xac dinh khoang cach giua cac truc')

disp('Vi hop so co dinh nen khoang cach truc so bo duoc xac dinh theo cong thuc')

disp('A=a.(Memax)^(1/3)')

disp('+Memax momen cuc dai cua dong co')

Trang 22

disp('+a la he so kinh nghiem')

a=15.88

while a<14.5 | a>16

a=input('hay nhap a khoang tren,a=')

end

A=a*(Memax)^(1/3)

round(A)

disp(' 3 chon modun banh rang');

disp('+ chon medun theo do thi');

disp(' doi voi o to, modun phap tuyen(mn) cua banh rang duoc chon theo do thi o hinh(I-3a) sach thiet ke hop so co chinh oto-may keo');

disp(' modun phan tuyrn phu thuoc vao');

disp('mo men quay cuc dai tren truc thu capcua hop so');

disp('chon modun cua cap banh rang tru rang thang la m=3,5');disp('chon modun cau cap banh rang nghieng la mn=3');

disp('goc beta chon theo kinh nghiem')

disp('doi voi o to con')

mn=3

m=3.5

beta=35

while beta<30 | beta>50

beta=input('ban nhap beta khoang tren,beta=')

end

disp('xac dinh so rang cua cac banh rang')

disp('sau khi da xac dinh duoc khoang cah truc so bo(A)giua cac truc,modun cua cac banh rang(m), goc nghieng cau cac banh rang (beta) thi so luong rang duoc xac dinh nhu sau')

disp('so luong rang cua banh rang chu dong cua cap banh rang lon luon an khop co thr chon theo dieu kien khong cat chan rang, nghia

disp('xac dinh ty so truyen cua cac cap banh rang duoc gai o cac

so truyen khac nhau cua hop so')

disp(' ig1 ty so truyen cua cap banh rang duoc gai o so 1')

disp(' ih1 ty so truyen cua cua hop so o so 1')

disp('ig1=ih1/ia')

disp('ig2=ih2/ia')

Trang 23

disp('tinh chinh xac khoang cach truc A')

disp('Sau khi da xac dinh chinh xac so luong cac banh rang va da lam tron so rang thi co the xay ra khoang cah truc khong con nhu nhau nua Vi vay ta co the tinh lai khoang cach truc A cua cac capbanh rang an khop voi nhau trong hop so')

disp('cap banh rang luon an khop')

goxi0=0.06015

disp('Neu Ac>A1 thi lamda0 va goxi0 se duong va anpha<20')

disp('Neu Ac<A1 thi lamda0 va goxi0 se duong va anpha>20')

disp('xac dinh he so dich chinh tong cong epchilont ')

goxit=0.5*goxi0*(zg1+zg11);

disp('goxit=')

disp(goxit)

Trang 24

disp('XAC DINH CAC THONG SO HINH HOC CUA RANG')

disp('doi voi cap banh rang nghieng luon an khop')

disp('ty so truyen')

iankhop=round(za1)/za

disp('mo dun phap tuyen mn=3.5')

disp('buoc phap tuyen')

Trang 25

disp('tong dich chinh tuong doichon theo phuc luc 4')

disp('he so dich chinh tong cong')

Trang 26

Trang 27

disp('+modun mat dau')

disp('goc nhieng cua rang betasl=0')

disp('duong kinh vong tron chia')

disp('Tinh Ben hop so')

disp('+ che do tai trong de tinh toan hop so')

disp(' a.momen truyen den cac truc hop so')

disp('truc trung gian')

disp('tu dong co truyen den')

Mtg=Memax*ia

disp('tu banh xe truyen den')

Mphimaxtg=(phimax*Gphi*10*rbx)/(icc*i0*If*ig1)

disp('truc thu cap')

disp('tu dong co truyen den')

Mtc1=Memax*ih1phay

Mtc2=Memax*ih2phay

Trang 28

Mtclui=Memax*igl

disp('+tu banh xe truyen den')

Mphimaxtc=(phimax*Gphi*rbx)/(ih1phay)

disp('b.tinh luc tac dung len cap banh rang')

disp('+cap rang luon an khop')

disp(' 2 tinh ben rang')

disp('a tinh ben suc uon')

kd=1.8

while kd<1.5 | kd>2.5

disp('ban nhap khoang tren,kd=')

disp('ban nhap khoang tren,ktp=')

end

kc=1

kgc=1

disp('doi voi rang tru rang thang')

disp('doi voi banh rang so 1')

kbeta1=1

disp('doi voi banh rang tru rang nghieng')

disp('doi voi banh rang an khop')

epchilonphay=(1.88-3.2*(1/za+1/round(za1)))*betaak

kbeta=1/epchilonphay

epchilon2phay=(1.88-3.2*(1/round(zg2)+1/round(zg22)))*betaso2kbeta2=1/epchilon2phay

epchilon3phay=(1.88-3.2*(1/round(zg3)+1/round(zg33)))*betaso3kbeta3=1/epchilon3phay

disp('y: he so dang rang phu thoc vao do thi he so dang rang hinhI-5 sach thiet ke tinh toan hop so chinh o to may keo')

ya=0.135

ya1=0.139

y1=0.125

y11=0.138

Trang 29

disp('doi cap banh rang luon an khop')

disp('banh rang chu dong')

disp('doi cap banh rang so 1')

disp('doi cap banh rang so 2')

Trang 30

while E<1.9 | E>2.2

disp('ban nhap khoang tren,E=')

end

disp('doi voi banh rang luon an khop')

3*sin(pi/9)*cos(pi/9)))*(10^3/ra+10^3/ra1))

Trang 31

disp('Tinh toan truc hop so')

disp('chon so kich thuc truc')

disp('duong kinh truc so cap')

disp('truc trung gian')

disp('tinh phan luc tai cac goi do')

Trang 32

disp('do vong cua truc so cap')

disp('trong mat phang zox')

jsc=(pi*dsc^4)/64

Q1=0

f1=(((R1+Ra1)*b2*(a2+b2))/(3*E*jsc)-(Q1*d01banhrang1*(2*a2+3*b2))/(6*E*jsc))*10^-6

disp('trong mat phang yox')

f2=(((P1+Pa1)*b2*(a2+b2))/(3*E*jsc)-(Q1*d01banhrang1*(2*a2+3*b2))/(6*E*jsc))*10^-6

Trang 33

disp('do vong truc thu cap')

disp('goc xoay cua truc so cap')

disp('goc xoay truc thu cap')

Trang 34

disp('chon o do theo bang 17p sach thiet ke chi tiet may')

disp('doi voi truc thu cap')

disp('momen quan tinh khoi luong cua cac banh rang')

disp('cap banh rang chung')

Trang 35

disp('thoi gian chuyen so')

disp('tu so thap len so cao')

Trang 36

disp('goc nghieng cua be mat ham')

k3=(nguy*Rms)/(sin((7*pi)/180)*Rbetac)

tan(beta)=k3

disp('tinh toan kiem tra cac thong so co ban cua dong toc')

disp('momen ma sat thuc te')