Nghiên cứu ảnh hưởng của một số tham số kết cấu đến tham số động lực học, độ giảm thanh, ổn định chuyển động và độ chụm của đạn cối giảm thanh nguyên lý pít tông ngược

Bộ giáo dục Đào tạo VIệN KHOA HọC Và công nghệ quân Đỗ ĐìNH LàO NGHIấN CU NH HƢỞNG CỦA MỘT SỐ THAM SỐ KẾT CẤU ĐẾN THAM SỐ ĐỘNG LỰC HỌC, ĐỘ GIẢM THANH, ỔN ĐỊNH CHUYỂN ĐỘNG VÀ ĐỘ CHỤM CỦA ĐẠN CỐI GIẢM THANH NGUYÊN LÝ PÍT TƠNG NGƢỢC Ln ¸n tiÕn SÜ kü tht Hà Nội - 2019 Bộ giáo dục Đào tạo Bộ quốc phòng VIệN KHOA học công nghệ quân Đỗ đình lào NGHIấN CU NH HNG CA MT SỐ THAM SỐ KẾT CẤU ĐẾN THAM SỐ ĐỘNG LỰC HỌC, ĐỘ GIẢM THANH, ỔN ĐỊNH CHUYỂN ĐỘNG VÀ ĐỘ CHỤM CỦA ĐẠN CỐI GIẢM THANH NGUN LÝ PÍT TƠNG NGC Chuyên ngành: Cơ kỹ thuật MÃ số: 52 01 01 Ln ¸n tiÕn sÜ kü tht Ng-êi h-íng dÉn khoa häc: PGS TS Bïi Ngäc Håi TS Đặng Hồng Triển H Ni - 2019 i LI CAM ĐOAN Tơi xin cam đoan cơng trình nghiên cứu riêng Các số liệu, kết trình bày luận án trung thực chƣa đƣợc cơng bố cơng trình khác, liệu tham khảo đƣợc trích dẫn đầy đủ Tác giả luận án Đỗ Đình Lào ii LỜI CẢM ƠN Tác giả luận án xin bày tỏ lòng biết ơn chân thành, sâu sắc tới thầy hƣớng dẫn: - PGS TS Bùi Ngọc Hồi, Viện Khoa học Công nghệ quân sự; - TS Đặng Hồng Triển, Trung tâm Nhiệt đới Việt Nga Đã trực tiếp hƣớng dẫn, tận tình bảo, tạo điều kiện tốt để tác giả hoàn thành luận án Tác giả xin chân thành cảm ơn Thủ trƣởng Viện Khoa học Công nghệ quân sự, Viện Tên lửa Phòng Đào tạo tạo điều kiện thuận lợi cho tác giả hoàn thành nhiệm vụ học tập nghiên cứu suốt trình làm luận án Tác giả xin chân thành cảm ơn Đảng ủy, huy Viện Vũ khí, Trung tâm Đo lƣờng Thử nghiệm vũ khí, Phịng Ngịi, Phịng Đạn nơi tơi cơng tác động viên, giúp đỡ tạo điều kiện thuận lợi suốt trình học tập, nghiên cứu thực luận án Tác giả xin chân thành cảm ơn thầy giáo, Nhà khoa học ngồi Qn đội, đồng nghiệp, gia đình ngƣời thân động viên, khích lệ, tận tình giúp đỡ tác giả suốt trình làm luận án Tác giả luận án iii MỤC LỤC Trang DANH MỤC CÁC KÝ HIỆU, CÁC CHỮ VIẾT TẮT vi DANH MỤC CÁC BẢNG ix DANH MỤC HÌNH VẼ x MỞ ĐẦU Chƣơng 1: TỔNG QUAN VỀ VŨ KHÍ GIẢM THANH VÀ ĐẠN CỐI GIẢM THANH .5 1.1 Nguyên nhân sinh âm bắn giải pháp đƣợc ứng dụng để giảm cho hệ vũ khí 1.1.1 Nguyên nhân sinh âm bắn 1.1.2 Các giải pháp đƣợc ứng dụng để giảm 1.2 Yêu cầu chất lƣợng sở đánh giá độ giảm thanh, ổn định chuyển động, độ chụm đạn cối giảm 12 1.2.1 Yêu cầu sở đánh giá độ giảm 12 1.2.2 Yêu cầu sở đánh giá ổn định chuyển động đạn 13 1.2.3 Yêu cầu sở đánh giá độ chụm đạn………………………… 14 1.2.4 Các yêu cầu chất lƣợng khác đạn………………………………….15 1.3 Tổng quan tình hình nghiên cứu đạn giảm ngun lý pít tơng thuận ngun lý pít tơng ngƣợc 16 1.3.1 Tình hình nghiên cứu nƣớc 16 1.3.2 Tình hình nghiên cứu nƣớc 17 1.4 Những vấn đề tồn hƣớng nghiên cứu luận án 21 1.4.1 Những vấn đề tồn 21 1.4.2 Hƣớng nghiên cứu luận án 21 1.5 Kết luận chƣơng 22 Chƣơng 2: XÂY DỰNG MÔ HÌNH KẾT CẤU VÀ MƠ HÌNH BÀI TỐN ĐỘNG LỰC HỌC CỦA ĐẠN CỐI GIẢM THANH 23 2.1 Mơ hình kết cấu q trình động lực học đạn cối giảm 23 2.1.1 Mơ hình kết cấu ngun lý làm việc đạn cối giảm 23 2.1.2 Quá trình động lực học đặc điểm thuật phóng đạn cối giảm 24 2.2 Xây dựng hệ phƣơng trình vi phân TPT đạn cối giảm 29 2.2.1 Các giả thiết xây dựng hệ PTVP TPT .29 2.2.2 Phƣơng trình biểu diễn quy luật sinh khí, tốc độ sinh khí thay đổi bề mặt cháy thuốc phóng….……………………………………….30 iv 2.2.3 Phƣơng trình trạng thái khí thuốc lòng xi lanh 31 2.2.4 Phƣơng trình cân lƣợng .33 2.2.5 Phƣơng trình khí qua khe hở pít tơng xi lanh .34 2.2.6 Phƣơng trình mơ tả chuyển động tịnh tiến xi lanh mang đầu đạn .35 2.2.7 Hệ phƣơng trình vi phân TPT đạn cối giảm 39 2.2.8 Phƣơng pháp giải hệ phƣơng trình vi phân TPT 41 2.2.9 Tính tốn TPT đạn cối giảm cỡ 50 mm lý thuyết 43 2.3 Xây dựng hệ PTVP chuyển động đạn có tính đến dao động trục đạn đánh giá ổn định chuyển động đạn cối giảm đƣờng bay 45 2.3.1 Xây dựng hệ PTVP chuyển động đạn có tính đến dao động trục đạn… 45 2.3.2 Phƣơng pháp giải hệ PTVP chuyển động đạn cối giảm có tính đến dao động trục đạn… …………………… ……… …… 51 2.3.3 Tính tốn kiểm chứng mơ hình tốn chuyển động đạn có tính đến dao động trục đạn đƣợc luận án xây dựng 55 2.4 Kết luận chƣơng 59 Chƣơng 3: KHẢO SÁT ẢNH HƢỞNG CỦA MỘT SỐ THAM SỐ KẾT CẤU ĐẾN THAM SỐ ĐỘNG LỰC HỌC, ĐỘ GIẢM THANH, ỔN ĐỊNH CHUYỂN ĐỘNG VÀ ĐỘ CHỤM CỦA ĐẠN CỐI GIẢM THANH 60 3.1 Khảo sát ảnh hƣởng số TSKC đến TSĐLH đạn cối giảm .60 3.1.1 Ảnh hƣởng khối lƣợng đạn khối lƣợng pít tơng 61 3.1.2 Ảnh hƣởng áp suất kéo đứt vành tai pít tông 63 3.1.3 Ảnh hƣởng hành trình chuyển động xi lanh 64 3.1.4 Ảnh hƣởng đƣờng kính xi lanh 65 3.1.5 Ảnh hƣởng thể tích buồng đốt 66 3.1.6 Ảnh hƣởng khe hở lắp ghép pít tông xi lanh 67 3.1.7 Ảnh hƣởng tổng hợp tham số kết cấu đạn 68 3.2 Khảo sát ảnh hƣởng số tham số kết cấu đến ổn định chuyển động nòng súng khả giảm đạn cối giảm .72 3.2.1 Trạng thái ứng suất, biến dạng pít tơng, xi lanh 73 3.2.2 Ảnh hƣởng số tham số kết cấu đến ổn định chuyển động nòng súng khả giảm đạn cối giảm thanh…….83 3.3 Khảo sát ảnh hƣởng tham số kết cấu đến ổn định chuyển động đƣờng bay độ chụm đạn cối giảm 88 v 3.3.1 Phƣơng pháp khảo sát dao động ổn định chuyển động đạn cối giảm đƣờng bay 88 3.3.2 Ảnh hƣởng số tham số kết cấu đến dao động ổn định chuyển động đạn cối giảm cỡ 50 mm đƣờng bay…… 89 3.3.3 Ảnh hƣởng tham số kết cấu đến độ chụm đạn cối giảm 97 3.4 Kết luận chƣơng 98 Chƣơng 4: NGHIÊN CỨU THỰC NGHIỆM 99 4.1 Mục đích nghiên cứu thực nghiệm 99 4.2 Chế tạo đạn cối giảm cỡ 50 mm cho thực nghiệm chọn súng thử nghiệm 99 4.2.1 Chế tạo đạn cối giảm cỡ 50 mm chọn súng thử nghiệm 100 4.2.2 Chế tạo đạn cối giảm cỡ 50 mm chuyên dùng lựa chọn súng thử nghiệm chuyên dùng 105 4.3 Thực nghiệm kiểm chứng mơ hình tính tốn lý thuyết 107 4.3.1 Thực nghiệm đo sơ tốc đạn, biến thiên áp suất khí thuốc theo thời gian lịng xi lanh ứng với số miền khe hở cặp pít tơng-xi lanh… 107 4.3.2 Thực nghiệm đo mức âm đạn cối giảm cỡ 50 mm….117 4.3.3 Thực nghiệm đo tầm bắn độ chụm đạn cối giảm cỡ 50 mm… .124 4.4 Kết luận chƣơng 128 KẾT LUẬN 129 DANH MỤC CÁC CƠNG TRÌNH KHOA HỌC ĐÃ CƠNG BỐ .131 TÀI LIỆU THAM KHẢO 132 PHỤ LỤC vi DANH MỤC CÁC KÝ HIỆU, CÁC CHỮ VIẾT TẮT 2e1 - Bề dày cháy thuốc phóng [m] B - Mơ men qn tính xích đạo, [N.m.s ] c′ cT - Hệ số khí đ - Hệ số khí đ cx - Hệ số lực c cy - Hệ số khí đ m D - Đƣờng kính xi lanh, [m] f - Lực thuốc phóng, [J/kg] fmoi - Lực thuốc mồi, [J/kg] Ik - Xung lƣợng tồn phần áp suất khí thuốc thời gian thuốc phóng cháy, [MPa.s] k - Chỉ số mũ đoạn nhiệt k1, k2 - Các hệ số hiệu chỉnh hệ phƣơng trình thuật phóng L - Quãng đƣờng chuyển động xi lanh thời điểm xét, [m] LD - Sai số trung tâm tản mát tầm, [m LH - Sai số trung tâm tản mát hƣớng LD / X - Độ chụm tầm LH / X - Độ chụm hƣớng Lđ - Hành trình chuyển động xi lan LTC - Chiều dài từ mũi đạn đến tâm cản LTT - Chiều dài từ mũi đạn đến trọng tâm Lvđ - Chiều dài viên đạn, [m] m - Khối lƣợng đạn (trừ pít tơng), [kg] mpt - Khối lƣợng pít tơng, [kg] p - Áp suất trung bình thuật phóng, [MPa] pdt - Áp suất kéo đứt vành tai pít tơng, [MPa] vii pđ - Áp suất khí thuốc thời điểm pít tơng đóng chặt vào xi lanh, [MPa] pk pmax pmoi - Áp suất thời điểm thuốc phóng cháy hết, [MPa] - Áp suất lớn khí thuốc, [MPa] - Áp suất mồi, [MPa] SM St R - Hằng số khí tđ S - Diện tích tiết diện mặt cắt ngang lịng xi lanh, [m ] 2 T1 - Diện tích tiết diện ngang lớn đạn (tiết diện Mi đen), [m ] Tck - Diện tích tiết diện ngang khe hở pít tơng-xi lanh, [m ] t - Thời gian chuyển động đạn thời điểm xét, [s] v0 vđ - Thời gian chuyển động xi lanh hành trình, [s] T - Nhiệt độ khí thuốc buồng đốt thời điểm xét, [ K] - Nhiệt độ cháy thuốc phóng [ K] - Chu kỳ dao động, [s] v - Vận tốc đầu đạn thời điểm xét, [m/s] - Vận tốc đạn rời khỏi cán truyền lực (sơ tốc đạn), [m/s] - Vận tốc đạn pít tơng va chạm đóng chặt vào xi lanh, [m/s] W - Thể tích riêng khí thuốc, [dm ] X - Tầm bắn, [m] z - Bề dày cháy tƣơng đối thuốc phóng α- Lƣợng cộng tích khí thuốc, [m3/kg]  - Góc chƣơng động khơng gian, [độ] δ0 - Góc lắc ban đầu mặ δd - Góc chƣơng động δn - Góc chƣơng động m δtp - Trọng lƣợng riêng δrxl - Biến dạng hƣớng kính δrpt - Biến dạng hƣớng kính δ0 - Tốc độ lắc ban đầu P2 Si=1+2*a; lamda=-(2*a+a^2)/Si; muy=a^2/Si; phi=1.001; % He so cong thu yeu phi2=0.242; % He so ton that luu luong phut %% Xac dinh cac gia tri ban dau del=(omega/Wo); % Mat nhoi Sd=(S-St); % Tiet dien dan c=(k-1)*phi*m/(f*omega*g*100/9.81); K0 = (2/(1 + k))^((k + 1)/2/(k - 1))*sqrt(g*k*10); p0=pmoi*10^4/9.81; xi0=(1/del-1/deltatp)/(fmoi*10^-6/pmoi+alpha1/deltatp); xichma0=sqrt(1+4*lamda*xi0/Si); z0=2*xi0/Si/(xichma0+1); W0=(Wo-(1-xi0)*omega/deltatp-alpha*omega*xi0)*10^3; to=1; neta0=0; v0=0; l0=0; t0=[0:0.000001: 0.01]; y0=[z0 xi0 neta0 v0 l0 W0 to p0]; options=odeset('events',@hevents); [T,Y]=ode45(@HPT DẠN COI GIAM THAHH,t0,y0,options); Y(:,8)=Y(:,8)*9.81/10000; %P= Mpa Y(:,4)=Y(:,4)/10; % v0 m/s Y(:,5)=Y(:,5)/10; %ld=m Tai thoi diem chay het [zk i1]=max(Y(Y(:,1)0 s2=1; else s2=0; end % He phuong trinh vi phan % Be day chay tuong doi cua thuốc phóng dy(1)=s1*(y(8)^Nuy)/(Ik*10^4/9.81); % Toc sinh dy(2)=Si*(1+2*lamda*y(1)+3*muy*y(1)^2)*dy(1); % Hieu suat phut qua khe ho dy(3)=phi2*K0*St*10^2*y(8)/(omega*sqrt(f*10/9.81*y(7))); % Van toc chuyen dong cua dan dy(4)=s2*Sd*100*y(8)*g*10/(phi*m); % Quang duong chuyen dong cua dan dy(5)=y(4); P6 % The tich khoang khong gian sau day dan dy(6)=omega*dy(2)/(deltatp*10^-3)-alpha*10^3*omega*(dy(2)dy(3))+S*100*dy(5); % Toc bien thien nhiet dy(7)=((1-y(7))*dy(2)-(k-1)*y(7)*dy(3)-c*y(4)*dy(4))/ (y(2)-y(3)); % Toc bien thien ap suat dy(8)=(f*10/9.81*omega*(y(7)*(dy(2)-dy(3))+(y(2)y(3))*dy(7))-y(8)*dy(6))/y(6); end function [qdg,df,dg]=hevents(t,y) % Dieu kien dung tai dau nong global ld qdg=y(5)-ld*10; df=1; dg=0; end P7 Phụ lục CHƢƠNG TRÌNH TÍNH QUỸ ĐẠO CHUYỂN ĐỘNG VÀ DAO ĐỘNG CỦA ĐẠN CỐI GIẢM THANH CỠ 50 mm TRONG MẶT PHẲNG ĐỨNG function DAO DONG_DUNG clc; clear all; close all; % Cac thong so dau vao global g ar an CX CYP bc br rokk B Lvd ct m Sm teta0=pi*45/180; % goc phong[radian] v0=77.1; delta0=0.55*57.3; % Toc lac ban dau, 1/s d=50/1000; % co dan [m] m=0.958+0.022;% kG g=9.81;% Gia toc truong [m/s2] td=10*pi/180; rokk=1.226; % Mat khong N.s2/m4 Sm=pi*d^2/4; % Dien tich miden cua dan, m2 Lvd=0.3195; % Chieu dai toan bo vien dan, m B=5837.27.10^-6; % Mo men quan tinh xich dao, N.m.s2 lk=146.5/d/1000; % Chieu dai than dan, lan co lp=65/d/1000; % Chieu dai canh on dinh, lan co lt=108/d/1000; % Chieu dai ong on dinh h1=1.75/d/1000; % Be rong canh, lan co fnd1=0.06; % He so tra theo thi n=8; % so canh LKK=[2.5 3.0 4.0 5.0 6.0]; M=[0.62 0.68 0.78 0.84 0.87]; if lk>3 cx=0.008*lk+0.017*(lp-1.1)+33*fnd1*h1+0.0045*(lk-3); else if lp>=1.1 cx=0.008*lk+0.017*(lp-1.1)+33*fnd1*h1; else cx=0.008*lk+33*fnd1*h1; P8 end end cy0=0.2; cy=cy0+0.011*(n-12)+0.16*(lp-0.5); k=(0.11/(1+0.8*(1+4/n)*lp))*lp; m1=interp1(LKK,M,lk); LTC=((lk+lt+0.25*lp)-0.02*lk*m1/(k))*d; % Vi tri tam can, m CX=cx*(1+30*td^2); % He so luc can chinh dien CYP=cy/td; LTT=0.132; % Vi tri tam, tinh tu mui, m CYDELTA=CX+CYP; ct=CYDELTA*(LTC-LTT)/Lvd; an=rokk/9.81*Sm*CYP*g/(2*m); ar=rokk/9.81*Sm*CX*g/(2*m); bc=rokk*Sm*Lvd^2*ct/B; br=rokk*Sm*Lvd*ct/(2*B); disp('chu ky dao dong dao tien'); T1=2*pi/(v0*sqrt(br)) disp('Chieu dai buoc song'); lamdabs=2*pi/ (sqrt(br)) disp('So mu tat dan') exilon=pi*(bc+an)/(2*sqrt(br)) % dieu kien bien v=[];teta=[];x=[];y=[]; delta=[]; T=[]; t=[]; v(1)=v0;teta(1)=teta0;x(1)=0;y(1)=0; delta(1)=0; T(1)=delta0; % Thuat toan runge_kutta h=0.005; i=1; t(1)=0; vi=[1 2 1]; while (y(i)>=0) k1(1)=h*dT(v(i),teta(i),T(i),delta(i)); k2(1)=h*ddelta(T(i)); k3(1)=h*dteta(v(i),teta(i),delta(i)); P9 k4(1)=h*dv(v(i),teta(i)); k5(1)=h*dx(v(i),teta(i)); k6(1)=h*dy(v(i),teta(i)); for j=2:4 k1(j)=h*dT(v(i)+k4(j-1)/vi(j),teta(i)+k3(j1)/vi(j),T(i)+k1(j-1)/vi(j),delta(i)+k2(j-1)/vi(j)); k2(j)=h*ddelta(T(i)+k1(j-1)/vi(j)); k3(j)=h*dteta(v(i)+k4(j-1)/vi(j),teta(i)+k3(j1)/vi(j),delta(i)+k2(j-1)/vi(j)); k4(j)=h*dv(v(i)+k4(j1)/vi(j),teta(i)+k3(j-1)/vi(j)); k5(j)=h*dx(v(i) +k4(j-1)/vi(j),teta(i)+k3(j-1)/vi(j)); k6(j)=h*dy(v(i)+k4(j-1)/vi(j),teta(i)+k3(j-1)/vi(j)); end T(i+1)=T(i) + (k1(1)+2*k1(2)+2*k1(3)+k1(4))/6; delta(i+1)=delta(i) + (k2(1)+2*k2(2)+2*k2(3)+k2(4))/6; teta(i+1)=teta(i) v(i+1)=v(i) x(i+1)=x(i) y(i+1)=y(i) i=i+1; t(i)=t(i-1)+h; end plot(x,delta,'b','linewidth',1.5); grid on ylabel('Goc Ch.dong (do)','FontSize',14); xlabel('Tam ban (m)','FontSize',13); set(gca,'YTick',[min(delta) -1 max(delta)]); ketqua=[delta']; xlswrite('ketqua_dung',ketqua,'Temperatures','B2' ); end % Cac ham phu %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function dv=dv(v,teta) global ar g dv=-ar*v^2g*sin(teta); end P10 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function dteta =dteta(v,teta,delta) global g an dteta=-g*cos(teta)/v+ an*v*delta; end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function dx=dx(v,teta) dx=v*cos(teta); end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function dy=dy(v,teta) dy=v*sin(teta); end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function dT=dT(v,teta,T,delta) global ar an g bc br dT=g*cos(teta)*(ar+bc+2*g*sin(teta)/v^2)-(bc+an)*v*T- br*v^2*delta; end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function ddelta=ddelta(T) ddelta=T; end P11 Phụ lục CHƢƠNG TRÌNH TÍNH DAO ĐỘNG CỦA ĐẠN CỐI GIẢM THANH CỠ 50 mm TRONG MẶT PHẲNG NGHIÊNG function DAO DONG_NGHIENG clc; clear all; close all; % Cac thong so dau vao global g ar an CX CYP bc br rokk B Lvd ct m Sm teta0=pi*45/180; % goc phong[radian] v0=77.1; delta0=0.55*57.3; % Toc lac ban dau, 1/s d=50/1000; % co dan [m] m=0.958+0.022;% kG g=0;% Gia toc truong tren mat phang nghieng[m/s2] td=10*pi/180; rokk=1.226; % Mat khong N.s2/m4 Sm=pi*d^2/4; % Dien tich miden cua dan, m2 Lvd=0.3195; % Chieu dai toan bo vien dan, m B=5837.27.10^-6; % Mo men quan tinh xich dao, N.m.s2 lk=146.5/d/1000; % Chieu dai than dan, lan co lp=65/d/1000; % Chieu dai canh on dinh, lan co lt=108/d/1000; % Chieu dai ong on dinh h1=1.75/d/1000; % Be rong canh, lan co fnd1=0.06; % He so tra theo thi n=8; % so canh LKK=[2.5 3.0 4.0 5.0 6.0]; M=[0.62 0.68 0.78 0.84 0.87]; if lk>3 cx=0.008*lk+0.017*(lp-1.1)+33*fnd1*h1+0.0045*(lk-3); else if lp>=1.1 cx=0.008*lk+0.017*(lp-1.1)+33*fnd1*h1; else P12 cx=0.008*lk+33*fnd1*h1; end end cy0=0.2; cy=cy0+0.011*(n-12)+0.16*(lp-0.5); k=(0.11/(1+0.8*(1+4/n)*lp))*lp; m1=interp1(LKK,M,lk); LTC=((lk+lt+0.25*lp)-0.02*lk*m1/(k))*d; % Vi tri tam can, m CX=cx*(1+30*td^2); % He so luc can chinh dien CYP=cy/td; LTT=0.132; % Vi tri tam, tinh tu mui, m CYDELTA=CX+CYP; ct=CYDELTA*(LTC-LTT)/Lvd; an=rokk/9.81*Sm*CYP*g/(2*m); ar=rokk/9.81*Sm*CX*g/(2*m); bc=rokk*Sm*Lvd^2*ct/B; br=rokk*Sm*Lvd*ct/(2*B); disp('chu ky dao dong dao tien'); T1=2*pi/(v0*sqrt(br)) disp('Chieu dai buoc song'); lamdabs=2*pi/ (sqrt(br)) disp('So mu tat dan') exilon=pi*(bc+an)/(2*sqrt(br)) % dieu kien bien v=[];teta=[];x=[];y=[]; delta=[]; T=[]; t=[]; v(1)=v0;teta(1)=teta0;x(1)=0;y(1)=0;delta(1)=0; T(1)=delta0; % Thuat toan runge_kutta h=0.005; i=1; t(1)=0; vi=[1 2 1]; while (y(i)>=0) k1(1)=h*dT(v(i),teta(i),T(i),delta(i)); k2(1)=h*ddelta(T(i)); P13 k3(1)=h*dteta(v(i),teta(i),delta(i)); k4(1)=h*dv(v(i),teta(i)); k5(1)=h*dx(v(i),teta(i)); k6(1)=h*dy(v(i),teta(i)); for j=2:4 k1(j)=h*dT(v(i)+k4(j-1)/vi(j),teta(i)+k3(j1)/vi(j),T(i)+k1(j-1)/vi(j),delta(i)+k2(j-1)/vi(j)); k2(j)=h*ddelta(T(i)+k1(j-1)/vi(j)); k3(j)=h*dteta(v(i)+k4(j-1)/vi(j),teta(i)+k3(j1)/vi(j),delta(i)+k2(j-1)/vi(j)); k4(j)=h*dv(v(i)+k4(j1)/vi(j),teta(i)+k3(j-1)/vi(j)); k5(j)=h*dx(v(i) +k4(j-1)/vi(j),teta(i)+k3(j-1)/vi(j)); k6(j)=h*dy(v(i)+k4(j-1)/vi(j),teta(i)+k3(j-1)/vi(j)); end T(i+1)=T(i) + (k1(1)+2*k1(2)+2*k1(3)+k1(4))/6; delta(i+1)=delta(i) + (k2(1)+2*k2(2)+2*k2(3)+k2(4))/6; teta(i+1)=teta(i) + (k3(1)+2*k3(2)+2*k3(3)+k3(4))/6; v(i+1)=v(i) + (k4(1)+2*k4(2)+2*k4(3)+k4(4))/6; x(i+1)=x(i) + (k5(1)+2*k5(2)+2*k5(3)+k5(4))/6; y(i+1)=y(i) + (k6(1)+2*k6(2)+2*k6(3)+k6(4))/6; i=i+1; t(i)=t(i-1)+h; end plot(x,delta,'b','linewidth',1.5); grid on ylabel('Goc Ch.dong (do)','FontSize',14); xlabel('Tam ban (m)','FontSize',13); set(gca,'YTick',[min(delta) -1 max(delta)]); ketqua=[delta']; xlswrite('ketqua_dung',ketqua,'Temperatures','B2' ); end % Cac ham phu %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function dv=dv(v,teta) global ar g dv=-ar*v^2g*sin(teta); P14 end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function dteta =dteta(v,teta,delta) global g an dteta=-g*cos(teta)/v+ an*v*delta; end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function dx=dx(v,teta) dx=v*cos(teta); end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function dy=dy(v,teta) dy=v*sin(teta); end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function dT=dT(v,teta,T,delta) global ar an g bc br dT=g*cos(teta)*(ar+bc+2*g*sin(teta)/v^2)-(bc+an)*v*T- br*v^2*delta; end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function ddelta=ddelta(T) ddelta=T; end P15 Phụ lục CHƢƠNG TRÌNH KHẢO SÁT ẢNH HƢỞNG TỔNG HỢP CÁC THAM SỐ KẾT CẤU ĐẾN THAM SỐ ĐỘNG LỰC HỌC CỦA ĐẠN CỐI GIẢM THANH function KHAO SAT_TONG HOP clear all; clc; close all; q=[0.98 0.984 0.988 0.992 0.996 1]; % Trong luong dan m=[0.02 0.021 0.022 0.023 0.024]; % Trong luong pitong ld=[1.5 1.52 1.54 1.56 1.58 1.6]; %Quang duong chuyen dong D=[0.2000 0.2005 0.201]; % Duong kinh xi lanh Sd=pi*d.^2/4; % pduttai=[14 15 16 17]*10000/9.81;% Ap suat keo dut vanh tai Skh=[1.5 1.7 1.9 2.1 2.3 2.5]*10^-4; % Dien tich tiet dien khe ho giua pitong va xi lanh W=[3.1*10^-3 3.15*10^-3 3.2*10^-3]; %The tich buong dot omega=[2.049*10^-3 2.050*10^-3 2.051*10^-3]; %Trong luong thuoc phong vi=1; for j1=1:length(q) for k1=1:length(m) for i1=1:length(ld) for m1=1:length(Sd) for o1=1:length(Skh) for p1=1:length(W) for n1=1:length(pduttai) for r1=1:length(omega) [PM(vi) VD(vi)] =KS_TPT(ld(i1),Sd(m1),Skh(o1),W(p1),q(j1),m(k1),pduttai(n1),o mega(r1)); D(vi)=d(m1)/10; %m LD(vi)=ld(i1)/10; %m S(vi)=Sd(m1); P16 SKH(vi)=Skh(o1)/100; W0(vi)=W(p1)/1000; Q(vi)=q(j1); M(vi)=m(k1); Pduttai(vi)=pduttai(n1)/10000*9.81; OMEGA(vi)=omega(r1); vi=vi+1; end end end end end end end end [Pmax, [Pmin, [Vmax, [Vmin, i1]= i2]= i3]= i4]= max(PM) min(PM) max(VD) min(VD) KETQUA=[Q' M' LD' D' SKH' W0' Pduttai' OMEGA' PM' VD']; xlswrite('KETQUA',KETQUA,'Temperatures','B2' ); end function [Pm,v0d]=KS_TPT(ld1,S1,St1,Wo1,q1,m1,pdutdai1,omega1) %Cac thong so ban dau global Ik Si lamda muy phi2 K0 St pduttai omega f S q delta Sd k c g phi ld Nuy ld=ld1; S=S1; St=St1; Wo=Wo1; q=q1-m1; % % % % % Chieu dai dan chay dm Tiet dien nong (dm2) Tiet dien khe ho(dm2) The tich ban dau cua buong (dm3) Trong luong dan tru pittong (kG) ... tài luận án: ? ?Nghiên cứu ảnh hưởng số tham số kết cấu đến tham số động lực học, độ giảm thanh, ổn định chuyển động độ chụm đạn cối giảm nguyên lý pít tơng ngược? ?? 2 Mục tiêu nghiên cứu - Xây dựng... công nghệ quân Đỗ đình lào NGHIÊN CỨU ẢNH HƢỞNG CỦA MỘT SỐ THAM SỐ KẾT CẤU ĐẾN THAM SỐ ĐỘNG LỰC HỌC, ĐỘ GIẢM THANH, ỔN ĐỊNH CHUYỂN ĐỘNG VÀ ĐỘ CHỤM CỦA ĐẠN CỐI GIM THANH NGUYấN Lí PT TễNG NGC... cối giảm thanh; Chƣơng Khảo sát ảnh hƣởng số tham số kết cấu đến tham số động lực học, độ giảm thanh, ổn định chuyển động độ chụm đạn cối giảm thanh; Chƣơng Nghiên cứu thực nghiệm Toàn nội dung