汽車?yán)碚撜n程設(shè)計

1、汽車?yán)碚撈嚴(yán)碚撛O(shè)計報告設(shè)計報告 汽車動力性、經(jīng)濟(jì)性的計算機(jī)模擬汽車動力性、經(jīng)濟(jì)性的計算機(jī)模擬 學(xué)號 姓名張少波張少波 學(xué)院交通學(xué)院交通學(xué)院 班級車輛本車輛本 15011501 班班 2017 年 12 月 18 日 目 錄 1汽車驅(qū)動力圖汽車驅(qū)動力圖.1 1.1汽車驅(qū)動力圖簡介.1 1.2汽車驅(qū)動力圖.2 2汽車驅(qū)動力汽車驅(qū)動力- -行駛阻力平衡圖行駛阻力平衡圖.2 2.1汽車驅(qū)動力-行駛阻力平衡圖簡介.2 2.1.1滾動阻力.3 2.1.2空氣阻力.3 2.2汽車驅(qū)動力-行駛阻力平衡圖.3 3汽車爬坡度圖汽車爬坡度圖.3 3.1汽車爬坡度圖簡介.3 3.2汽車爬坡度圖.4 4汽車加速度圖汽

2、車加速度圖.5 4.1汽車加速度圖簡介.5 4.2汽車加速度圖.5 5汽車動力特性圖汽車動力特性圖.6 5.1汽車動力特性圖簡介.6 5.2汽車動力特性圖.6 6汽車功率平衡圖汽車功率平衡圖.6 6.1汽車功率平衡圖簡介.6 6.2汽車功率平衡圖.7 7汽車百公里油耗圖汽車百公里油耗圖.7 7.1汽車百公里油耗圖簡介.7 7.2汽車百公里油耗圖.8 參考文獻(xiàn)參考文獻(xiàn) .8 附錄附錄 .8 汽車動力性、經(jīng)濟(jì)性的計算機(jī)模擬 張少波 （交通學(xué)院，車輛工程，1501 班， ） 摘要：摘要：通過 MATLAB 計算機(jī)軟件進(jìn)行汽車動力性、經(jīng)濟(jì)性的計算機(jī)模擬，模擬 得出汽車的各種特性曲線。包括：驅(qū)動力圖，驅(qū)

3、動力-行駛阻力平衡圖，爬坡度圖， 加速度圖，動力特性圖，功率平衡圖，百公里油耗圖。同時對汽車特性曲線的計算 公式加以統(tǒng)計匯總學(xué)習(xí)。 關(guān)鍵詞：關(guān)鍵詞：驅(qū)動力圖，驅(qū)動力-行駛阻力平衡圖，爬坡度圖，加速度圖，動力特 性圖，功率平衡圖，百公里油耗圖，汽車特性曲線的計算公式 1汽車驅(qū)動力圖 1.1汽車驅(qū)動力圖簡介 驅(qū)動力是由發(fā)動機(jī)的轉(zhuǎn)矩經(jīng)傳動系統(tǒng)傳至驅(qū)動輪上得到的。汽車發(fā)動機(jī)產(chǎn) 生的轉(zhuǎn)矩，經(jīng)傳動系傳至驅(qū)動輪上。此時作用于驅(qū)動輪上的轉(zhuǎn)矩 Tt產(chǎn)生一對地 面的圓周力 F0,地面對驅(qū)動輪的反作用力 Ft（方向與 F0相反）既是驅(qū)動汽車的 外力，此外力稱為汽車的驅(qū)動力。其數(shù)值為 （1.1） 式中，Tt為作用于驅(qū)

4、動輪上的轉(zhuǎn)矩；r 為車輪半徑。 作用于驅(qū)動輪上的轉(zhuǎn)矩 Tt是由發(fā)動機(jī)產(chǎn)生的轉(zhuǎn)矩經(jīng)傳動系傳至車輪上的。 若令 Ttq表示發(fā)動機(jī)轉(zhuǎn)矩，ig 表示變速器的傳動比，i0表示主減速器的傳動比， 表示傳動系的機(jī)械效率，則有 T （1.2） 對于裝有分動器、輪邊減速器、液力傳動等裝置的汽車，上式應(yīng)計入相應(yīng) 的傳動比和機(jī)械效率。 因此驅(qū)動力為 （1.3） 汽油發(fā)動機(jī)使用外特性的 Tq-n 曲線的擬合公式為 （1.4） 一般用根據(jù)發(fā)動機(jī)外特性確定的驅(qū)動力與車速之間的函數(shù)關(guān)系曲線 Ft-ua來 全面表示汽車的驅(qū)動力，稱為汽車的驅(qū)動力圖。設(shè)計中的汽車有了發(fā)動機(jī)的外 r T 0 t Tgtq iiT Tgtq iiT

5、 0t T r Tt t F 432 ) 1000 (8445 . 3 ) 1000 (874.40) 1000 (44.165) 1000 (27.25913.19Tq nnnn 特性曲線、傳動系的傳動比、傳動效率、車輪半徑等參數(shù)后，即可用式（1.1） 求出各個檔位的 Ft值，再根據(jù)發(fā)動機(jī)轉(zhuǎn)速與汽車行駛速度之間轉(zhuǎn)換關(guān)系求出 ua，即可求得各個檔位的 Ft于 ua曲線。發(fā)動機(jī)轉(zhuǎn)速與汽車行駛速度之間的關(guān)系 式為 （1.5） 式中，ua為汽車行駛速度（km/h）;n 為發(fā)動機(jī)轉(zhuǎn)速（r/min）;r 為車輪半徑 （m）; ig為變速器傳動比；i0為主減速器傳動比。 另：實際行駛中，發(fā)動機(jī)常在節(jié)氣門部

6、分開啟下工作，相應(yīng)的驅(qū)動力要比它小。 1.2 汽車驅(qū)動力圖 其車輛基本參數(shù)見附錄 1，其 Matlab 程序見附錄 2。 附圖： 2汽車驅(qū)動力-行駛阻力平衡圖 2.1汽車驅(qū)動力-行駛阻力平衡圖簡介 汽車在水平道路上等速行駛時，必須克服來自地面的滾動阻力和來自空氣 0 a 377. 0u ii rn g 的空氣阻力。滾動阻力以符號 Ff表示，空氣阻力以符號 Fw表示。當(dāng)汽車在坡道 上上坡行駛時，還必須克服重力沿坡道的分力，稱為坡道阻力，以符號 Fi表示。 汽車加速行駛時還需要克服加速阻力，以符號 Fj表示。因此，汽車行駛的總阻 力為 jiwf FFFFF （2.1） 在汽車驅(qū)動力圖上把汽車行駛中

7、經(jīng)常遇到的滾動阻力和空氣阻力一并畫出， 做出汽車驅(qū)動力-行駛阻力平衡圖。 2.1.1 滾動阻力 滾動阻力等于滾動阻力系數(shù)與車輪負(fù)載之乘積。即 Wf f F （2.2） 其中滾動阻力系數(shù)由試驗確定。滾動阻力系數(shù)與路面的種類、行駛車速以 及輪胎的構(gòu)造、材料、氣壓等有關(guān)。這里選取滾動阻力系數(shù)為良好的瀝青或混 凝土路面滾動阻力系數(shù)。 （詳見附錄 1） 。 2.1.2 空氣阻力 汽車直線行駛時受到的空氣作用力在行駛方向上的分力稱為空氣阻力。此 處討論無風(fēng)條件下汽車的運(yùn)動，空氣阻力計算公式為 （2.3） 式中，A 為迎風(fēng)面積，即汽車行駛方向的投影面積（m2） ；CD為空氣阻力系 數(shù)；ua2為車速，km/h

8、。 （詳見附錄 1） 。 2.2 汽車驅(qū)動力-行駛阻力平衡圖 其車輛基本參數(shù)見附錄 1，其 Matlab 程序見附錄 3。 15.21 F 2 w aDAu C 3汽車爬坡度圖 3.1汽車爬坡度圖簡介 汽車的上坡能力是用滿載（或某一載質(zhì)量）時汽車在良好路面上的最大爬 坡度 imax表示的。另外 imax代表了汽車的極限爬坡能力，它應(yīng)比實際行駛中遇 到的道路最大坡度超出很多，這是因為應(yīng)考慮到在實際坡道行駛時，在坡道上 停車后順利起步加速、克服松軟坡道路面的大阻力、克服坡道上崎嶇不平路面 的局部大阻力等要求的緣故。 也可以理解汽車的爬坡能力，指汽車在良好路面上克服 Ff+Fw后的余力全部 用來（即

9、等速）克服坡度阻力時能爬上的坡度。 因此 Fi=Ft-（Ff+Fw） （3.1） 一般汽車最大爬坡度達(dá) 30%左右，因此利用汽車行駛方程式確定 1 檔及低 檔爬坡能力時，應(yīng)采用 Gsin 作為坡道阻力，即上式應(yīng)為 （3.2） ) 15.21 A cos(sin 2 D 0 a Tgtq uC fG r iiT G G FFF wft )( sinarc 即 （3.3） 根據(jù)驅(qū)動力-行駛阻力圖求出汽車能爬上之坡道角，相應(yīng)地根據(jù) tan=i 可求 出坡度值。 3.2汽車爬坡度圖 其車輛基本參數(shù)見附錄 1，其 Matlab 程序見附錄 4。 4汽車加速度圖 4.1汽車加速度圖簡介 汽車的加速能力可用

10、它在水平良好路面上行駛時能產(chǎn)生的加速度來評價。 而汽車加速行駛時，需克服的其質(zhì)量加速運(yùn)動時的慣性力，就是加速阻力 Fj 。汽車的質(zhì)量分為平移質(zhì)量和旋轉(zhuǎn)質(zhì)量兩部分。為了簡易計算旋轉(zhuǎn)質(zhì)量，常用 系數(shù) 作為計算計入旋轉(zhuǎn)質(zhì)量慣性力偶矩后的汽車旋轉(zhuǎn)質(zhì)量換算系數(shù)。 主要 與飛輪的轉(zhuǎn)動慣量、車輪的轉(zhuǎn)動慣量以及傳動系的傳動比有關(guān)。 其公式如下： （4.1） 2 2 0 2 2 1 m 1 1 r iiI mr I Tgf w 式中，Iw為車輪的轉(zhuǎn)動慣量（kgm2） ；If為飛輪的轉(zhuǎn)動慣量（kgm2） ； i0為主傳動比；ig為變速器的速比。 汽車加速時其驅(qū)動力需全力加速，即坡道阻力為 0 。 4.2汽車加速度

11、圖 其車輛基本參數(shù)見附錄 1，其 Matlab 程序見附錄 5。 5汽車動力特性圖 5.1汽車動力特性圖簡介 用動力特性圖來分析汽車動力性。 將汽車行駛方程兩邊除以汽車重力并整理如下： （5.1） (5.2) 令 為汽車的動力因數(shù)并以符號 D 表示，則 dt duF FFFF t jwif G F- F w t G Fw t F dt du D （5.3） 汽車在各檔下的動力因數(shù)與車速的關(guān)系曲線稱為動力特性圖。 5.2汽車動力特性圖 其車輛基本參數(shù)見附錄 1，其 Matlab 程序見附錄 6。 6汽車功率平衡圖 6.1汽車功率平衡圖簡介 正文部分用宋體小四號 在汽車行駛的每一瞬間，發(fā)動機(jī)發(fā)出的

12、功率始終 等于機(jī)械傳動損失功率與全部運(yùn)動阻力多消耗的功率。 汽車運(yùn)動阻力所消耗的功率有滾動阻力功率 Pf、空氣阻力功率 Pw、坡度阻 力功率 Pi及加速阻力功率 Pj。 將汽車行駛方程式兩邊乘以行駛車速 ua，并經(jīng)單位換算整理出汽車功率平 衡方程式（式中功率單位為 kW） ，如下 （6.1） ) 3600761403600 i 3600 ( 1 P 3 e dt dumuAuCuGGfu aaDaa T 將汽車經(jīng)常遇到的阻力功率 （Pf+Pw） 。發(fā)動機(jī)功率對車速的曲線繪制。 得到汽車功率平衡圖。 6.2汽車功率平衡圖 其車輛基本參數(shù)見附錄 1，其 Matlab 程序見附錄 7。 7汽車百公里

13、油耗圖 7.1汽車百公里油耗圖簡介 正文部分用宋體小四號 根據(jù)參考文獻(xiàn)1，等速行駛工況燃油消耗量的 計算。在萬有特性圖上有等燃油消耗率曲線。根據(jù)曲線可以確定發(fā)動機(jī)在一定 轉(zhuǎn)速 n,發(fā)出一定功率 Pe時的燃油消耗率 b。 （7.1） （7.2） T 1 4 4 3 3 2 2101 Bb eeee PBPBPBPB gu bP a e 02.1 Qs （7.3） 7.2汽車百公里油耗圖 其車輛基本參數(shù)見附錄 1，其 Matlab 程序見附錄 8。 參考文獻(xiàn) 1 余志生.汽車?yán)碚揗.北京：機(jī)械工業(yè)出版社，2009.3. 附錄 附錄 1 車輛基本參數(shù)為發(fā)動機(jī)最低轉(zhuǎn)速 600r/min ， 最高轉(zhuǎn)速

14、4000r/min ， 裝載質(zhì) 量 2000kg ，整車整備質(zhì)量 1800kg ,車輪半徑 0.367m ，傳動系效率 0.85 , 滾動阻力系數(shù) 0.013 ，空氣阻力系數(shù)與迎風(fēng)面積之積為 2.77m2 ， 主減速器 傳動比 5.83 ,飛輪轉(zhuǎn)動慣量 0.218kg m2 ,前輪總轉(zhuǎn)動慣量為 1.798kg m2 ,后 輪總的轉(zhuǎn)動慣量為 3.598kg m2,軸距為 3.2m 質(zhì)心至前軸距離 1.947m ，質(zhì)心 高 0.9m ， 四擋變速器的變速比分別為 6.09 ， 3.09 ， 1.71 ， 1 。 LQs17.12 附錄 2 close all; n=linspace(600,400

15、0,100);%設(shè)定轉(zhuǎn)速范圍(linspace(x1,x2,N)在 x1x2 范圍內(nèi)選取 N 個點(diǎn)，均分。) ua1=0.377*0.367*n/(5.83*6.09);%公式計算各檔車速范圍 ua2=0.377*0.367*n/(5.83*3.09); ua3=0.377*0.367*n/(5.83*1.71); ua4=0.377*0.367*n/(5.83*1); Ttq=(-19.13+259.27*(n/1000)-165.44*(n/1000).2+40.874*(n/1000).3- 3.8445*(n/1000).4);%?Ttq-n 曲線擬合公式 x1=6.09; x2=3.

16、09; x3=1.71; x4=1; Ft1=x1*Ttq*5.83*0.85/0.367;%公式計算各檔對應(yīng)轉(zhuǎn)速下的驅(qū)動力 Ft2=x2*Ttq*5.83*0.85/0.367; Ft3=x3*Ttq*5.83*0.85/0.367; Ft4=x4*Ttq*5.83*0.85/0.367; plot(ua1,Ft1,g,ua2,Ft2,r,ua3,Ft3,y,ua4,Ft4,b);%汽車驅(qū)動力-行駛 阻力平衡圖 xlabel(ua/(km/h);%汽車驅(qū)動力-行駛阻力平衡圖 ylabel(F/N); legend(1,2,3,4); grid on; title(汽車驅(qū)動力圖); 附錄 3

17、close all; n=linspace(600,4000);%設(shè)定轉(zhuǎn)速范圍 ua1=0.377*0.367*n/(5.83*6.09);%公式計算各檔車速范圍 ua2=0.377*0.367*n/(5.83*3.09); ua3=0.377*0.367*n/(5.83*1.71); ua4=0.377*0.367*n/(5.83*1); Ttq=(-19.13+259.27*(n/1000)-165.44*(n/1000).2+40.874*(n/1000).3- 3.8445*(n/1000).4);%?Ttq-n 曲線擬合公式 x1=6.09; x2=3.09; x3=1.71; x4

18、=1; Ft1=x1*Ttq*5.83*0.85/0.367;%公式計算各檔對應(yīng)轉(zhuǎn)速下的驅(qū)動力 Ft2=x2*Ttq*5.83*0.85/0.367; Ft3=x3*Ttq*5.83*0.85/0.367; Ft4=x4*Ttq*5.83*0.85/0.367; F1=3800*9.8*0.013+(2.77/21.15)*ua1.2;%公式計算各檔對應(yīng)的各個車速下的行駛阻 力 F2=3800*9.8*0.013+(2.77/21.15)*ua2.2; F3=3800*9.8*0.013+(2.77/21.15)*ua3.2; F4=3800*9.8*0.013+(2.77/21.15)*ua

19、4.2; plot(ua1,Ft1,g,ua2,Ft2,b,ua3,Ft3,y,ua4,Ft4,m,ua1,F1,r,ua2,F 2,r,ua3,F3,r,ua4,F4,r);%汽車驅(qū)動力-行駛阻力平衡圖? xlabel(ua/(km/h);%汽車驅(qū)動力-行駛阻力平衡圖 ylabel(F/N); legend(1,2,3,4,行駛阻力); grid on; title(汽車驅(qū)動力-行駛阻力平衡圖); 附錄 4 close all; n=600:10:4000; Tq=-19.13+259.27*(n/1000)-165.44*(n/1000).2+40.874*(n/1000).3- 3.84

20、45*(n/1000).4; m=3800;g=9.8; nmin=600;nmax=4000;G=m*g; ig=6.09 3.09 1.71 1 ; nT=0.85;%傳動系效率 r=0.367;%車輪半徑 f=0.013;%滾動阻力系數(shù) CDA=2.77;%空氣阻力系數(shù)與迎風(fēng)面積之積 i0=5.83;%主減速器傳動比 L=3.2;%軸距 a=1.947;%質(zhì)心至前軸距離 hg=0.9;%質(zhì)心高 If=0.218;%飛輪轉(zhuǎn)動慣量 Iw1=1.798;%前輪總轉(zhuǎn)動慣量 Iw2=3.598;%后輪總轉(zhuǎn)動慣量 Ft1=Tq*ig(1)*i0*nT/r; Ft2=Tq*ig(2)*i0*nT/r;

21、 Ft3=Tq*ig(3)*i0*nT/r; Ft4=Tq*ig(4)*i0*nT/r; ua1=0.377*r*n/ig(1)/i0; ua2=0.377*r*n/ig(2)/i0; ua3=0.377*r*n/ig(3)/i0; ua4=0.377*r*n/ig(4)/i0; Fw1=CDA*ua1.2/21.15; Fw2=CDA*ua2.2/21.15; Fw3=CDA*ua3.2/21.15; Fw4=CDA*ua4.2/21.15; Ff=G*f; i1=asin(Ft1-Ff-Fw1)/G); i2=asin(Ft2-Ff-Fw2)/G); i3=asin(Ft3-Ff-Fw3)

22、/G); i4=asin(Ft4-Ff-Fw4)/G); plot(ua1,i1,ua2,i2,ua3,i3,ua4,i4); title(汽車的爬坡度圖); xlabel(ua/(km*h-1); ylabel(i/%); legend(1,2,3,4); grid on; 附錄 5 close all; n=600:10:4000; Tq=-19.13+259.27*(n/1000)-165.44*(n/1000).2+40.874*(n/1000).3- 3.8445*(n/1000).4; m=3880;g=9.8; nmin=600;nmax=4000; G=m*g; ig=6.09

23、 3.09 1.71 1.00; nT=0.85;r=0.367; f=0.013;CDA=2.77;i0=5.83; L=3.2;a=1.947;hg=0.9; If=0.218;Iw1=1.798;Iw2=3.598; Ft1=Tq*ig(1)*i0*nT/r; Ft2=Tq*ig(2)*i0*nT/r; Ft3=Tq*ig(3)*i0*nT/r; Ft4=Tq*ig(4)*i0*nT/r; ua1=0.377*r*n/ig(1)/i0; ua2=0.377*r*n/ig(2)/i0; ua3=0.377*r*n/ig(3)/i0; ua4=0.377*r*n/ig(4)/i0; Fw1=

24、CDA*ua1.2/21.15; Fw2=CDA*ua2.2/21.15; Fw3=CDA*ua3.2/21.15; Fw4=CDA*ua4.2/21.15; Ff=G*f; deta1=1+(Iw1+Iw2)/(m*r2)+(If*ig(1)2*i02*nT)/(m*r2);%汽車旋轉(zhuǎn)質(zhì)量換 算系數(shù) deta2=1+(Iw1+Iw2)/(m*r2)+(If*ig(2)2*i02*nT)/(m*r2); deta3=1+(Iw1+Iw2)/(m*r2)+(If*ig(3)2*i02*nT)/(m*r2); deta4=1+(Iw1+Iw2)/(m*r2)+(If*ig(4)2*i02*nT)/

25、(m*r2); a1=(Ft1-Ff-Fw1)/(deta1*m); a2=(Ft2-Ff-Fw2)/(deta2*m); a3=(Ft3-Ff-Fw3)/(deta3*m); a4=(Ft4-Ff-Fw4)/(deta4*m); plot(ua1,a1,ua2,a2,ua3,a3,ua4,a4); title(汽車的行駛加速度曲線); xlabel(ua/(km*h-1); ylabel(a/(m*s-2); legend(1,2,3,4); grid on; 附錄 6 close all; n=linspace(600,4000);%設(shè)定轉(zhuǎn)速范圍 ua1=0.377*0.367*n/(5.

26、83*6.09);%公式計算各檔車速范圍 ua2=0.377*0.367*n/(5.83*3.09); ua3=0.377*0.367*n/(5.83*1.71); ua4=0.377*0.367*n/(5.83*1.00); Ttq=(-19.13+259.27*(n/1000)-165.44*(n/1000).2+40.874*(n/1000).3- 3.8445*(n/1000).4);%Ttq-n 曲線擬合公式 m=3800;g=9.8; G=m*g; x1=6.69; x2=3.09; x3=1.71; x4=1; Ft1=x1*Ttq*5.83*0.85/0.367;%公式計算各檔

27、對應(yīng)轉(zhuǎn)速下的驅(qū)動力 Ft2=x2*Ttq*5.83*0.85/0.367; Ft3=x3*Ttq*5.83*0.85/0.367; Ft4=x4*Ttq*5.83*0.85/0.367; F1=(2.77/21.15)*ua1.2;%公式計算各檔對應(yīng)的各個車速下的空氣阻力 F2=(2.77/21.15)*ua2.2; F3=(2.77/21.15)*ua3.2; F4=(2.77/21.15)*ua4.2; D1=(Ft1-F1)/G; D2=(Ft2-F2)/G; D3=(Ft2-F2)/G; D4=(Ft3-F3)/G; plot(ua1,D1,ua2,D2,ua3,D3,ua4,D4);

28、 title(汽車的動力特性圖); xlabel(ua/(km*h-1); ylabel(D); legend(1,2,3,4) grid on; 附錄 7 close all; n=linspace(600,4000);%設(shè)定轉(zhuǎn)速范圍? ua1=0.377*0.367*n/(5.83*6.09);%公式計算各檔車速范圍? ua2=0.377*0.367*n/(5.83*3.09); ua3=0.377*0.367*n/(5.83*1.71); ua4=0.377*0.367*n/(5.83*1.00); Ttq=(-19.13+259.27*(n/1000)-165.44*(n/1000).

29、2+40.874*(n/1000).3- 3.8445*(n/1000).4);%?Ttq-n?曲線擬合公式? x1=6.09; x2=3.09; x3=1.71; x4=1; Ft1=x1*Ttq*5.83*0.85/0.367;%公式計算各檔對應(yīng)轉(zhuǎn)速下的驅(qū)動力? Ft2=x2*Ttq*5.83*0.85/0.367; Ft3=x3*Ttq*5.83*0.85/0.367; Ft4=x4*Ttq*5.83*0.85/0.367; F1=3800*9.8*0.013+(2.77/21.15)*ua1.2;%公式計算各檔對應(yīng)的各個車速下的行駛阻 力? F2=3800*9.8*0.013+(2.7

30、7/21.15)*ua2.2; F3=3800*9.8*0.013+(2.77/21.15)*ua3.2; F4=3800*9.8*0.013+(2.77/21.15)*ua4.2; Pe1=Ft1.*ua1./(3600*0.85);%計算各檔對應(yīng)轉(zhuǎn)速下的功率? Pe2=Ft2.*ua2./(3600*0.85); Pe3=Ft3.*ua3./(3600*0.85); Pe4=Ft4.*ua4./(3600*0.85); P1=F1.*ua1./(3600*0.85);%計算各檔對應(yīng)的各個車速下的行駛阻力功率? P2=F2.*ua2./(3600*0.85); P3=F3.*ua3./(3600*0.85); P4=F4.*ua4./(3600*0.85); plot(ua1,Pe1,b,ua2,Pe2,y,ua3,Pe3,m,ua4,Pe4,y,ua1,P1