版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認(rèn)領(lǐng)
文檔簡介
1、汽車系統(tǒng)動力學(xué)matlab作業(yè)報告小組成員:一、組內(nèi)任務(wù)分配二、matlab程序與圖形1、不同轉(zhuǎn)向特性車輛在不同車速下的系統(tǒng)特征根m=1000;i=1500;a1=1.15;b1=1.35;caf=53000;car=53000;i=1;r=;for uc=10:5:100; d=(i*(caf+car)+m*(a12*caf+b12*car)/(m*i*uc); s=(a1+b1)2*caf*car/(m*i*uc2)+(b1*car-a1*caf)/i; p=1 d s; r=roots(p); r(i,1)=r(1,1);r(i,2)=r(2,1);i=i+1;endplot(real(
2、r(:,1),imag(r(:,1),bo);holda2=1.25;b2=1.25;t=1;s=; for uc=10:5:100 p=m 0;0 i;q=(caf+car)/uc,m*uc+(a2*caf-b2*car)/uc;(a2*caf-b2*car)/uc,(a22*caf+b22*car)/uc; r=caf;a2*caf; a=-p(-1)*q; d=eig(a); i=imag(d); r=real(d); s(t,1)=r(1); s(t,2)=i(1); t=t+1;endplot(s(:,1),s(:,2),*)a3=1.35;b3=1.15; for uc=10:5:
3、100 p=m 0;0 i; q=(caf+car)/uc,m*uc+(a3*caf-b3*car)/uc;(a3*caf-b3*car)/uc,(a32*caf+b32*car)/uc; r=caf;a3*caf; a=-p(-1)*q; d=eig(a); i=imag(d); r=real(d); s(t,1)=r(1); s(t,2)=i(1); t=t+1;endgrid onplot(s(:,1),s(:,2),d);axis(-14 2 0 3);xlabel(實軸(re);ylabel(虛軸(im);text(-8,2.8,不足轉(zhuǎn)向);text(0,0.2,過多轉(zhuǎn)向);text
4、(-3,0.2,中性轉(zhuǎn)向)set(gca,fontname,helvetica,fontsize,10)title(不同轉(zhuǎn)向特性車輛在不同車速下的系統(tǒng)特征根,fontsize,12); 2.1、具有不同轉(zhuǎn)向特性車輛的橫擺角速度幅頻和相頻響應(yīng)m=1000;i=1500;a1=1.15;a2=1.25;a3=1.35;b1=1.35;b2=1.25;b3=1.15;caf=53000;car=53000;d=;c=;m=;uc=20;i=1;for f=0.1:0.1:10 w=2*pi*f; e=1 0;0 1; p=m 0;0 i; r1=caf;a1*caf; r2=caf;a2*caf;
5、r3=caf;a3*caf; q1=(caf+car)/uc m*uc+(a1*caf-b1*car)/uc;(a1*caf-b1*car)/uc (a12*caf+b12*car)/uc; q2=(caf+car)/uc m*uc+(a2*caf-b2*car)/uc;(a2*caf-b2*car)/uc (a22*caf+b22*car)/uc; q3=(caf+car)/uc m*uc+(a3*caf-b3*car)/uc;(a3*caf-b3*car)/uc (a32*caf+b32*car)/uc; a1=-inv(p)*q1; a2=-inv(p)*q2; a3=-inv(p)*q
6、3; b1=inv(p)*r1; b2=inv(p)*r2; b3=inv(p)*r3; hw1=-inv(a1-1i*w*e)*b1; hw2=-inv(a2-1i*w*e)*b2; hw3=-inv(a3-1i*w*e)*b3; d(i,1)=angle(hw1(2)*180/pi;d(i,2)=abs(hw1(2); c(i,1)=angle(hw2(2)*180/pi;c(i,2)=abs(hw2(2); m(i,1)=angle(hw3(2)*180/pi;m(i,2)=abs(hw3(2); i=i+1;endsubplot(2,1,1)f=0.1:0.1:10;semilogx(
7、f,d(:,2),k,f,c(:,2),r-.,f,m(:,2)grid onset(gca,xtick,0.1,0.3,1.5,3,10)set(gca,fontname,helvetica,fontsize,10)legend(不足轉(zhuǎn)向,中性轉(zhuǎn)向,過多轉(zhuǎn)向)title(具有不同轉(zhuǎn)向特性車輛的橫擺角速度幅頻響應(yīng),fontsize,12);xlabel(頻率/hz)ylabel(橫擺角速度增益/(/s)/()subplot(2,1,2)f=0.1:0.1:10;semilogx(f,d(:,1),k,f,c(:,1),r-.,f,m(:,1)m=1000;i=1500;a1=1.15;a2=1
8、.25;a3=1.35;b1=1.35;b2=1.25;b3=1.15;caf=53000;car=53000;d=;c=;m=;uc=20;i=1;for f=0.1:0.1:10 w=2*pi*f; e=1 0;0 1; p=m 0;0 i; r1=caf;a1*caf; r2=caf;a2*caf; r3=caf;a3*caf; q1=(caf+car)/uc m*uc+(a1*caf-b1*car)/uc;(a1*caf-b1*car)/uc (a12*caf+b12*car)/uc; q2=(caf+car)/uc m*uc+(a2*caf-b2*car)/uc;(a2*caf-b2
9、*car)/uc (a22*caf+b22*car)/uc; q3=(caf+car)/uc m*uc+(a3*caf-b3*car)/uc;(a3*caf-b3*car)/uc (a32*caf+b32*car)/uc; a1=-inv(p)*q1; a2=-inv(p)*q2; a3=-inv(p)*q3; b1=inv(p)*r1; b2=inv(p)*r2; b3=inv(p)*r3; hw1=-inv(a1-1i*w*e)*b1; hw2=-inv(a2-1i*w*e)*b2; hw3=-inv(a3-1i*w*e)*b3; d(i,1)=angle(hw1(2)*180/pi;d(
10、i,2)=abs(hw1(2); c(i,1)=angle(hw2(2)*180/pi;c(i,2)=abs(hw2(2); m(i,1)=angle(hw3(2)*180/pi;m(i,2)=abs(hw3(2); i=i+1;endsubplot(2,1,1)f=0.1:0.1:10;semilogx(f,d(:,2),k,f,c(:,2),r-.,f,m(:,2)grid onset(gca,xtick,0.1,0.3,1.5,3,10)set(gca,fontname,helvetica,fontsize,10)legend(不足轉(zhuǎn)向,中性轉(zhuǎn)向,過多轉(zhuǎn)向)title(具有不同轉(zhuǎn)向特性車
11、輛的橫擺角速度幅頻響應(yīng),fontsize,12);xlabel(頻率/hz)ylabel(橫擺角速度增益/(/s)/()subplot(2,1,2)f=0.1:0.1:10;semilogx(f,d(:,1),k,f,c(:,1),r-.,f,m(:,1)grid onset(gca,xtick,0.1,0.3,1.5,3,10) set(gca,fontname,helvetica,fontsize,10)legend(不足轉(zhuǎn)向,中性轉(zhuǎn)向,過多轉(zhuǎn)向)title(具有不同轉(zhuǎn)向特性車輛的橫擺角速度相頻響應(yīng),fontsize,12);xlabel(頻率/hz)ylabel(橫擺角速度相位())gr
12、id onset(gca,xtick,0.1,0.3,1.5,3,10) set(gca,fontname,helvetica,fontsize,10)legend(不足轉(zhuǎn)向,中性轉(zhuǎn)向,過多轉(zhuǎn)向)title(具有不同轉(zhuǎn)向特性車輛的橫擺角速度相頻響應(yīng),fontsize,12);xlabel(頻率/hz)ylabel(橫擺角速度相位())2.2、具有不同轉(zhuǎn)向特性車輛的側(cè)向加速度和相頻響應(yīng)m=1000;i=1500;a1=1.15;a2=1.25;a3=1.35;b1=1.35;b2=1.25;b3=1.15;caf=53000;car=53000;d=;c=;m=;uc=20;l=a1+b1;i=
13、1;for f=0.1:0.1:10 w=2*pi*f; e=1 0;0 1; vi=w*i*caf; rr=l*caf*car/uc; vr1=(l*b1*caf*car/uc-m*a1*caf*uc); vr2=(l*b2*caf*car/uc-m*a2*caf*uc); vr3=(l*b3*caf*car/uc-m*a3*caf*uc); ri1=w*m*a1*caf; ri2=w*m*a2*caf; ri3=w*m*a3*caf; dr1=-w2*m*i+l2*caf*car/(uc2)+m*(b1*car-a1*caf); dr2=-w2*m*i+l2*caf*car/(uc2)+m
14、*(b2*car-a2*caf); dr3=-w2*m*i+l2*caf*car/(uc2)+m*(b3*car-a3*caf); di1=w*(i*(caf+car)+m*(a12*caf+b12*car)/uc; di2=w*(i*(caf+car)+m*(a22*caf+b22*car)/uc; di3=w*(i*(caf+car)+m*(a32*caf+b32*car)/uc; hay1=1i*w*(vr1+1i*vi)/(dr1+1i*di1)+uc*(rr+1i*ri1)/(dr1+1i*di1); hay2=1i*w*(vr2+1i*vi)/(dr2+1i*di2)+uc*(rr
15、+1i*ri2)/(dr2+1i*di2); hay3=1i*w*(vr3+1i*vi)/(dr3+1i*di3)+uc*(rr+1i*ri3)/(dr3+1i*di3); d(i,1)=angle(hay1)*180/pi;d(i,2)=abs(hay1); c(i,1)=angle(hay2)*180/pi;c(i,2)=abs(hay2); m(i,1)=angle(hay3)*180/pi;m(i,2)=abs(hay3); i=i+1;endsubplot(2,1,1)f=0.1:0.1:10;semilogx(f,d(:,2)*(pi/(180*9.8),k,f,c(:,2)*(p
16、i/(180*9.8),r-.,f,m(:,2)*(pi/(180*9.8)grid onset(gca,xtick,0.1 0.3 1 3 10)set(gca,fontname,helvetica,fontsize,10)legend(不足轉(zhuǎn)向,中性轉(zhuǎn)向,過多轉(zhuǎn)向)title(具有不同轉(zhuǎn)向特性車輛的側(cè)向加速度幅頻響應(yīng),fontsize,12);xlabel(頻率/hz)ylabel(側(cè)向加速度增益/g/()subplot(2,1,2)f=0.1:0.1:10;semilogx(f,d(:,1),k,f,c(:,1),r-.,f,m(:,1)grid onlegend(不足轉(zhuǎn)向,中性轉(zhuǎn)向,過
17、多轉(zhuǎn)向)set(gca,xtick,0.1 0.3 1 3 10)set(gca,fontname,helvetica,fontsize,10)title(具有不同轉(zhuǎn)向特性車輛的側(cè)向加速度幅頻響應(yīng),fontsize,12);axis(0.1 10 -120 120)xlabel(頻率/hz)ylabel(側(cè)向加速度相位())2.3、 abs控制器設(shè)計w0=120;v0=30;tb=600;ki=4500;kd=5000;ts=0.05;iw=12;rd=0.25;uh=0.8;ug=0.6;s0=0.2;m=300;g=9.8;w=;v=;s=;w(1)=120;v(1)=30;s(1)=0;
18、k=0;i=2;while(v00) k=k+1; sb=(v0-rd*w0)/v0; s(i)=sb; if sb0.22 tb=tb-kd*ts; elseif sb0.18 tb=tb+ki*ts; else tb=tb; end w0=w0+(fxb*rd-tb)/iw*ts; w(i)=w0; v0=v0-fxb/m*ts; v(i)=v0; i=i+1;endt=0:ts:ts*k;subplot(2,1,1)plot(t,v,t,w*rd);grid onlegend(車輪前進速度,車輪線速度);set(gca,fontname,helvetica,fontsize,10)tit
19、le(車輪前進速度與車輪線速度關(guān)系曲線,fontsize,12);xlabel(時間/(s));ylabel(速度/(m/s));axis(0,ts*k,0,32);subplot(2,1,2)plot(t,s);axis(0,ts*k,0,1);grid onset(gca,fontname,helvetica,fontsize,10)title(abs控制的滑移率時域仿真結(jié)果,fontsize,12);xlabel(時間/(s));ylabel(滑移率);set(gca,ytick,0:0.2:1)3、扭振系統(tǒng)振型圖j1=1.986*10(-3);1.910*10(-3);1.931*10
20、(-3);1.931*10(-3);1.910*10(-3);1.924*10(-3);7.8426*10(-2);2.258*10(-3);2.641*10(-2);2.91*10(-3);2.51*10(-3);1.77*10(-3);7.836*10(-2);3.238;k1=7.95*104;7.95*104;6.95*104;7.95*104;7.95*104;6.90*104; 8.93*103;1.41*104;1.02*104;4.4*103;1.38*104;1.616*102;2.51*102;j=diag(j1);k(1,1)=k1(1);k(1,2)=-k1(1);i=
21、2;for i=2:13; k(i,i-1)=-k1(i-1); k(i,i)=k1(i-1)+k1(i); k(i,i+1)=-k1(i);endk(14,13)=-k1(13);k(14,14)=k1(13);a=inv(j)*k;g,d=eig(a);f=sqrt(d)/(2*pi); m=1; while m=14 g(:,m)=g(:,m)/g(1,m) m=m+1endsubplot(3,2,1)plot(g(:,13)xlabel(質(zhì)點號);set(gca,fontname,helvetica,fontsize,10)grid on title(單節(jié)點振型圖(4.26hz),fo
22、ntsize,12)axis(1 15 -0.2 1.1) subplot(3,2,2) plot(g(:,12)xlabel(質(zhì)點號);set(gca,fontname,helvetica,fontsize,10)grid onaxis(1 15 -5 1.2)title(雙節(jié)點振型圖(11.99hz),fontsize,12)subplot(3,2,3)plot(g(:,11)xlabel(質(zhì)點號);set(gca,fontname,helvetica,fontsize,10)grid onaxis(1 15 -7 1.2)title(三節(jié)點振型圖 (73.9hz),fontsize,12
23、)subplot(3,2,4)plot(g(:,10)xlabel(質(zhì)點號);set(gca,fontname,helvetica,fontsize,10)grid ontitle(四節(jié)點振型圖 (147.7hz),fontsize,12)axis(1 15 -9 37)subplot(3,2,5)plot(g(:,9)xlabel(質(zhì)點號);set(gca,fontname,helvetica,fontsize,10)grid onaxis(1 15 -0.2 1.1)title(五節(jié)點振型圖 (252hz),fontsize,12)5、福特granada轎車后懸架單輪模型頻率響應(yīng)函數(shù)和系統(tǒng)
24、相應(yīng)輸入功率譜密度mb=317.5;mw=45.4;ks=22000;kt=192000;cs=1500;go=5*10(-6);uc=20;n=1;b=;d=;for f=0:0.01:15 w=2*pi*f; sf=4.47*10(-4)*power(f,-2.5); a=i*cs*w+(ks+kt-w2*mw),-i*cs*w-ks;-i*cs*w-ks,i*cs*w+(ks-mb*w2); c=kt;0; d=-a(-1)*c; b(1,n)=abs(d(2)-d(1); b(2,n)=abs(-d(2)*w2); b(3,n)=abs(kt*(d(1)-1); b(4,n)=sf;
25、b(5,n)=b(2,n)2*sf; b(6,n)=b(1,n)2*sf; b(7,n)=b(3,n)2*sf; n=n+1;endf=0:0.01:15subplot(3,2,1)plot(f,b(1,:);xlabel(頻率/hz)ylabel(懸架動行程增益)grid onsubplot(3,2,2);plot(f,b(2,:);grid onxlabel(頻率/hz)ylabel(不舒適性參數(shù)增益/(m/s2)/m)subplot(3,2,3);plot(f,b(3,:);grid onxlabel(頻率/hz)ylabel(輪胎動載荷增益/n/m)subplot(3,2,4)plot
26、(f,b(6,:);xlabel(頻率/hz)ylabel(懸架動行程功率譜密度/m2/hz)grid onsubplot(3,2,5);plot(f,b(5,:);grid onxlabel(頻率/hz)ylabel(不舒適性功率譜密度/(m/s2)2/hz)subplot(3,2,6);plot(f,b(7,:);grid onxlabel(頻率/hz)ylabel(輪胎動載荷功率譜密度/n2/hz)6、魔術(shù)公式xdata = 0 0.78 1.88 2.79 3.80 4.82 6.29 7.82 9.31 11.80;ydata = 0 7.99 16.05 20.65 23.82 26.00 27.
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 2025年寶雞貨運從業(yè)資格證試題庫及答案
- 2025年玉溪貨運考試題目
- 2025年丹東c1貨運從業(yè)資格證考試題
- 2024年市政工程建設(shè)項目借款合同3篇
- 2024年木材供需合同3篇
- 2024年度高品質(zhì)天然石材荒料采購與銷售合同3篇
- 2025存量房購買居間合同
- 2025種薯買賣合同范本
- 2024年版智慧醫(yī)療系統(tǒng)研發(fā)合同
- 2024年標(biāo)準(zhǔn)化汽車維修服務(wù)合同版B版
- 2023-2024學(xué)年廣州大附屬中學(xué)中考一模物理試題含解析
- 綠化養(yǎng)護工作日記錄表
- 2024美的在線測評題庫答案
- 2024版高考數(shù)學(xué)二輪復(fù)習(xí):解析幾何問題的方法技巧
- 輿情監(jiān)測服務(wù)方案
- 北京市海淀區(qū)2023-2024學(xué)年八年級上學(xué)期期末英語試卷
- 果品類原料的烹調(diào)應(yīng)用課件
- 地彈簧行業(yè)分析
- 如何發(fā)揮采購在公司高質(zhì)量發(fā)展中作用
- 民事糾紛及其解決機制課件
- 美術(shù)高考總結(jié)匯報
評論
0/150
提交評論