活塞式壓縮機(jī)設(shè)計(jì)說(shuō)明

1、活塞式壓縮機(jī)的設(shè)計(jì)說(shuō)明姓名：班級(jí)：學(xué)號(hào)：指導(dǎo)老師：1.題目：復(fù)算19WY-9/150型氫氦氣壓縮機(jī)在目前操作條件下的各級(jí)壓力、排氣溫度、排氣量、功率，作出計(jì)算示功圖、切向力圖、活塞力圖、標(biāo)明最大活塞力與切向力，核算配用電機(jī)功率是否適當(dāng)？2.19WY-9/150型氫氦氣壓縮機(jī)簡(jiǎn)介19WY-9/150型氫氦氣壓縮機(jī)是我省投產(chǎn)3000噸小型化肥廠的氮?dú)錃鈮嚎s機(jī)，二列之間為飛輪，由電機(jī)通過(guò)三角皮帶拖動(dòng)。壓縮機(jī)為臥式、兩列、門(mén)型、四線壓縮。原料（半水煤氣）經(jīng)脫硫后進(jìn)入I級(jí)，經(jīng)I級(jí)壓縮后送去變換、水洗、碳化，碳化后為碳化氣。碳化氣返回II級(jí)、III、IV級(jí)壓縮后去洗銅、合成。3.目前操作條件與有關(guān)數(shù)據(jù)（1

2、）操作條件：吸氣壓力：0.15MPa（絕）排氣壓力：16.0MPa（絕）I級(jí)出口與II級(jí)進(jìn)口壓力差為DP=0.09MPa吸氣溫度：級(jí)IIIIIIIVTs()30354040I級(jí)進(jìn)口相對(duì)濕度F=1（2）氣體組成組分N2H2CO2COO2ArCH4半水煤氣0.210.390.120.260.0050.00450.0105碳化氣0.2370.7090.0050.0290.0030.0050.012（3）有關(guān)數(shù)據(jù)：活塞行程：S=310mm，活塞桿直徑d=60mm轉(zhuǎn)速：n=209rpm，連桿長(zhǎng)度l=700mm；I、IV列超前II、III列90度往復(fù)運(yùn)動(dòng)件重量：IIV列210.9kg；IIIII列193.

3、7kg級(jí)IIIIIIIV直徑D34021013565相對(duì)余隙a0.0650.150.160.1抽氣系數(shù)m10.890.890.89飛輪矩GD2為471.0kgm2，配用電機(jī)額定功率：155kw。設(shè)計(jì)計(jì)算一 . 計(jì)算各級(jí)的行程容積。I級(jí)：II級(jí)：III級(jí)：IV級(jí)：二 計(jì)算各級(jí)名義壓力和名義壓力比已知 三 排氣溫度 1 .絕熱指數(shù)計(jì)算首先需算出各氣體的絕熱指數(shù)，查活塞式壓縮機(jī)附表1常用氣體物理性質(zhì)得：各氣體絕熱指數(shù)組分絕熱指數(shù)1.401.411.311.401.401.681.32注：表中所列的氣體絕熱指數(shù)為時(shí)的參數(shù)，但對(duì)于理想氣體絕熱指數(shù)的大小與溫度的關(guān)系不大，顧可不考慮溫度的影響，可把k看做是

4、常數(shù)?；旌蠚怏w的k值計(jì)算式為：。對(duì)于半水煤氣來(lái)說(shuō)有1k-1=0.211.40-1+0.391.41-1+0.121.31-1+0.261.40-1+0.0051.40-1+0.00451.68-1+0.01051.32-1 =2.5652k=1.3898同理對(duì)碳化氣有：1k-1=0.2371.40-1+0.7091.41-1+0.0051.31-1+0.0291.40-1+0.0031.40-1+0.0051.68-1+0.0121.32-1 =2.4627k=1.4061所k1=1.3898 k2=k3=k4=1.4061。 2 計(jì)算排氣溫度 由式可得級(jí)次吸氣溫度排氣溫度壓力比I303442

5、.073.845II308486.374.864III313447.813.456IV313379.871.955四 第一次精度計(jì)算 1 各級(jí)容積系數(shù) 根據(jù)過(guò)程流體機(jī)械表2-1計(jì)算膨脹指數(shù)m得： 對(duì)于理想氣體，容積系數(shù) 2 計(jì)算工況系數(shù)C和精度B 由 精度不達(dá)標(biāo) 進(jìn)行第一次復(fù)算五 第一次復(fù)算 1. 計(jì)算各級(jí)進(jìn)出口壓力 由式計(jì)算各級(jí)壓力 2 計(jì)算各級(jí)壓力比 3 確定膨脹系數(shù) 根據(jù)過(guò)程流體機(jī)械表2-1計(jì)算膨脹指數(shù)得： 4 計(jì)算容積系數(shù) 5 計(jì)算工況系數(shù)和精度 精度不達(dá)標(biāo) 進(jìn)行第二次復(fù)算六 第二次復(fù)算 重復(fù)第一次復(fù)算步驟 ，列出計(jì)算結(jié)果如下表 級(jí)IIIIIIIV進(jìn)氣壓力0.150.53842.325

6、57.0617排氣壓力0.62842.32557.016716壓力比4.18934.31933.03662.2657容積系數(shù)0.84940.68960.79730.9211工況常數(shù)C''0.046200.046620.047860.04858精度B''0.95<0.97七 第三次復(fù)算 重復(fù)第一次復(fù)算步驟 ，列出計(jì)算結(jié)果如下表 級(jí)IIIIIIIV進(jìn)氣壓力Psi'''0.150.53832.28526.889排氣壓力Pdi'''0.62832.28526.88916壓力比i'''4.188

7、74.24523.01462.3225容積系數(shù)vi'''0.84950.69570.79930.9179工況常數(shù)C'''0.046200.047020.047150.04723精度B'''0.98>0.97八 修正排氣溫度 由式，可得 級(jí)次吸氣溫度排氣溫度壓力比I303452.814.1887II308467.624.2452III313430.483.0146IV313399.242.3225九 排氣量計(jì)算 1 容積系數(shù) 2 確定壓力系數(shù)對(duì)于第一級(jí)，壓力系數(shù)可取為0.950.98，故??；對(duì)于其余各級(jí)，壓力系數(shù)可取為

8、0.981，故??； 3 確定溫度系數(shù) 根據(jù)各級(jí)壓力比，查過(guò)程流體機(jī)械圖2-12可得： ； ； ； 4 確定泄露系數(shù)對(duì)于泄露系數(shù)的選取，有如下規(guī)則：對(duì)于不嚴(yán)密，延遲關(guān)閉的閥門(mén)：；對(duì)于單作用式氣缸的活塞環(huán)：；對(duì)于雙作用式氣缸的活塞環(huán)：；對(duì)于填料，j為級(jí)次。泄露系數(shù)計(jì)算表級(jí)次IIIIIIIV氣閥I0.02II0.02III0.02IV0.02活塞環(huán)I0.009II0.009III0.03IV0.03填料I0.0005II0.00100.0010III0.00150.00150.0015IV0.00200.00200.00200.0020總相對(duì)泄漏量0.03400.03350.05350.052泄露系

9、數(shù)0.96710.96760.94920.9506 5 綜上并且將結(jié)果列表如下級(jí)次行程容積容積系數(shù)壓力系數(shù)溫度系數(shù)泄露系數(shù)轉(zhuǎn)速n（r/min)I0.054390.84950.970.940.9671 209II0.017040.69570.990.940.9676III0.003560.79930.990.950.9492IV0.001030.91790.990.950.9506 6 計(jì)算排氣量 由式，并代入上表數(shù)據(jù)可得 十 計(jì)算功率并核算電機(jī)功率是否適當(dāng) 1 計(jì)算各級(jí)指示功率和總指示功率 計(jì)算相對(duì)壓力損失 根據(jù)各級(jí)進(jìn)排氣壓力，查過(guò)程流體機(jī)械圖2-14可得： 計(jì)算實(shí)際循環(huán)指示功率 由公式 取n

10、1=0.97k1 總指示功率 2 計(jì)算軸功率 取機(jī)械效率 則軸功率 3 計(jì)算驅(qū)動(dòng)機(jī)功率 取傳動(dòng)效率 則驅(qū)動(dòng)機(jī)功率 考慮到實(shí)際運(yùn)行當(dāng)中壓縮機(jī)會(huì)受到載荷波動(dòng)，氣體泄漏等因素影響，會(huì)引起功率的增加，驅(qū)動(dòng)機(jī)因留有的儲(chǔ)備功率，故所選用壓縮機(jī)的功率為,取壓縮機(jī)功率為1100kW，而配用電機(jī)額定功率為155kW ，故不適當(dāng)。十 示功圖繪制示功圖及壓力指示圖是壓縮機(jī)循環(huán)過(guò)程在P-V圖上的表示，從圖中可反映出循環(huán)的功耗。在作示功圖時(shí)，可將膨脹變壓和壓縮變壓過(guò)程可看做是絕熱膨脹和壓縮過(guò)程，其體積和壓力之間存在著關(guān)系式：，C為常數(shù)。對(duì)于I,II,III,IV級(jí)汽缸，其相關(guān)參數(shù)列表如下：級(jí)次IIIIIIIV行程容積0

11、.054390.017040.003560.00103相對(duì)余隙容積0.0650.150.160.1余隙容積0.0035350.0025560.00056960.000103容積系數(shù)0.84950.69570.79930.91790.046200.011850.0028460.0009454絕熱指數(shù)1.38981.40611.40611.40610.057930.019600.004130.001130.02070.007010.00190.000620.0035350.0025560.00056960.0001030.011730.0077460.0012840.000188 利用Matlab

12、畫(huà)出示功圖，圖及程序如下：第一級(jí)：>> x1=0.01173:0.00001:0.05793;>> y1=0.15;>> x2=0.0207:0.00001:0.05793; >> y2=0.15*0.057931.3898.*x2.(-1.3898); >> x3=0.003535:0.00001:0.0207;>> y3=0.6283;>> x4=0.003535:0.00001:0.01173; >> y4=0.6283*0.0035351.3898.*x4.(-1.3898);>>

13、; plot(x1,y1,x2,y2,x3,y3,x4,y4)第二級(jí)：>> x1=0.007646:0.00001:0.01960;>> y1=0.5383;>> x2=0.00701:0.00001:0.01960; >> y2=0.5383*0.019601.4061.*x2.(-1.4061);>> x3=0.002556:0.00001:0.00701;>> y3=2.2852;>> x4=0.002556:0.00001:0.007746; >> y4=2.2852*0.0025561.4

14、061.*x4.(-1.4061);>> plot(x1,y1,x2,y2,x3,y3,x4,y4)第三級(jí)：>> x1=0.001284:0.00001:0.00413;>>y1=2.2852>>x2=0.0019:0.00001:0.00413; >> y2=2.2852*0.004131.4061.*x2.(-1.4061);>>x3=0.0005696:0.00001:0.0019;>> y3=6.889;>> x4=0.0005696:0.00001:0.001284; >> y

15、4=6.889*0.00056961.4061.*x4.(-1.4061);>> plot(x1,y1,x2,y2,x3,y3,x4,y4)第四級(jí)：>> x1=0.000188:0.00001:0.00113;>>y1=6.889;>>x2=0.00062:0.00000001:0.00113; >> y2=6.889*0.001131.4061.*x2.(-1.4061);>>x3=0.000103:0.00001:0.00062;>>y3=16;>>x4=0.000103:0.00000001:

16、0.000188; >> y4=16*0.0001031.4061.*x4.(-1.4061);>>plot(x1,y1,x2,y2,x3,y3,x4,y4)十一 綜合活塞力的計(jì)算 綜合活塞力為活塞中的氣體力、往復(fù)慣性力和往復(fù)摩擦力的代數(shù)和 1 氣體力計(jì)算 確定蓋側(cè)面積和軸側(cè)面積 活塞運(yùn)動(dòng)分析活塞的位移公式為： 活塞的加速度為：式中：其中I-IV列運(yùn)動(dòng)的角度為，II-III列為，以 II-III列為例，則蓋側(cè)對(duì)應(yīng)于氣缸容積的活塞當(dāng)量位移：軸側(cè)對(duì)應(yīng)于氣缸容積的活塞當(dāng)量位移： 計(jì)算任意過(guò)程中的某處壓力 膨脹過(guò)程： 吸氣過(guò)程： 壓縮過(guò)程： 排氣過(guò)程： 其中余隙容積的當(dāng)量行程S

17、C=S* I-IV列 氣體力計(jì)算 蓋側(cè)氣體力計(jì)算 膨脹過(guò)程: 吸氣過(guò)程： 壓縮過(guò)程： 排氣過(guò)程 ：軸側(cè)氣體力計(jì)算 壓縮過(guò)程： 排氣過(guò)程：膨脹過(guò)程：吸氣過(guò)程 II-III列 氣體力計(jì)算 蓋側(cè)氣體力計(jì)算 膨脹過(guò)程:吸氣過(guò)程 壓縮過(guò)程 排氣過(guò)程 ： 軸側(cè)氣體力計(jì)算 壓縮過(guò)程： 排氣過(guò)程 ：膨脹過(guò)程: 吸氣過(guò)程：2 往復(fù)慣性力的計(jì)算 對(duì)I-IV列 其中一階往復(fù)慣性力：=二階往復(fù)慣性力：=3.47 對(duì)II-III列一階往復(fù)慣性力：=二階往復(fù)慣性力：=3.18因?yàn)橐凰牧谐岸?0°，所以=+90 一階往復(fù)慣性力合力： = + = 14.38cos15.66sin 二階往復(fù)慣性力合力：+=3.

18、183.473 往復(fù)摩擦力的計(jì)算 取，式中為第i列指示功率；為壓縮機(jī)的機(jī)械效率；s為活塞行程，n為壓縮機(jī)轉(zhuǎn)速 對(duì)I-IV列 對(duì)II-III列 十二 活塞力圖的繪制1）由以上的計(jì)算結(jié)果，利用matlab繪制I-IV列活塞力圖及程序如下：蓋側(cè)活塞力程序及Matlab圖：>> a1=0:0.0001:0.5793;b1=-(3.2254./(5.6461-cos(a1)-4.5161.*(1-0.04903.*(sin(a1).2).(1./2).1.3898+5.5263./(5.7161-cos(a1)-4.5161.*(1-0.04903.*(sin(a1).2).(1./2).1

19、.4061);a2=0.5793:0.0001:pi;b2=-35.61;a3=pi:0.0001:5.01;b3=-(37.5274./(5.6461-cos(a3)-4.5161.*(1-0.04903.*(sin(a3).2).(1./2).1.3898+69.3051./(5.7161-cos(a3)-4.5161.*(1-0.04903.*(sin(a3).2).(1./2).1.4061);a4=5.01:0.0001:2.*pi;b4=-108.0774;plot(a1,b1,a2,b2,a3,b3,a4,b4)蓋側(cè)活塞力軸測(cè)活塞力Matlab程序及Matlab圖：>>

20、; a1=0:0.0001:1.84;b1=37.7362./(-3.3861+cos(a1)+4.5161.*(1-0.04903.*(sin(a1).2).0.5).1.3898;a2=1.84:0.0001:pi;b2=55.27;a3=pi:0.0001:3.915;b3=3.2438./(-3.3861+cos(a3)+4.5161.*(1-0.04903.*(sin(a3).2).0.5).1.3898;a4=3.915:0.0001:2.*pi;b4=13.19;plot(a1,b1,a2,b2,a3,b3,a4,b4)軸測(cè)氣體力圖摩擦力和往復(fù)慣性力Matlab程序及Matlab

21、圖： 對(duì)I-IV列>> a1=0:0.0001:2.*pi;b1=16.71;a2=0:0.0001:2.*pi;b2=15.15.*(cos(a2)+0.2214.*cos(2.*a2);plot(a1,b1,a2,b2)摩擦力和往復(fù)慣性力圖綜合活塞力的matlab程序及matlab圖：>> a1=0:0.0001:0.5793;b1=-(3.2254./(5.6461-cos(a1)-4.5161.*(1-0.04903.*(sin(a1).2).(1./2).1.3898+5.5263./(5.7161-cos(a1)-4.5161.*(1-0.04903.*(s

22、in(a1).2).(1./2).1.4061)+37.7362./(-3.3861+cos(a1)+4.5161.*(1-0.04903.*(sin(a1).2).0.5).1.3898+16.34;a2=0.5793:0.0001:1.8392;b2=-35.61+37.7362./(-3.3861+cos(a2)+4.5161.*(1-0.04903.*(sin(a2).2).0.5).1.3898+15.15.*(cos(a2)+0.2214.*cos(2.*a2)+16.34;a3=1.8392:0.0001:pi;b3=-35.61+55.27+15.15.*(cos(a3)+0.

23、2214.*cos(2.*a3)+16.34;a4=pi:0.0001:3.9151;b4=-(37.5274./(5.6461-cos(a4)-4.5161.*(1-0.04903.*(sin(a4).2).(1./2).1.3898+69.3051./(5.7161-cos(a4)-4.5161.*(1-0.04903.*(sin(a4).2).(1./2).1.4061)+ 3.2438./(-3.3861+cos(a4)+4.5161.*(1-0.04903.*(sin(a4).2).0.5).1.3898+15.15.*(cos(a4)+0.2214.*cos(2.*a4)+16.3

24、4;a5=3.9151:0.0001:5.01;b5=-(37.5274./(5.6461-cos(a5)-4.5161.*(1-0.04903.*(sin(a5).2).(1./2).1.3898+69.3051./(5.7161-cos(a5)-4.5161.*(1-0.04903.*(sin(a5).2).(1./2).1.4061)+13.19+15.15.*(cos(a5)+0.2214.*cos(2.*a5)+16.34;a6=5.01:0.0001:2.*pi;b6=-108.0774+13.19+15.15.*(cos(a6)+0.2214.*cos(2.*a6)+16.34;

25、plot(a1,b1,a2,b2,a3,b3,a4,b4,a5,b5,a6,b6)綜合活塞力圖綜合活塞力的最大值約在1.8弧度的位置，其大小約為25kN2） 由以上的計(jì)算結(jié)果，利用matlab繪制II-III列活塞力圖及程序如下：蓋側(cè)活塞力程序及Matlab圖：>> a1=0:0.0001:0.9934;b1=-79.16.*(0.3./(5.8161-cos(a1)-4.5161.*(1-0.04903.*(sin(a1).2).0.5).1.4061;a2=0.9934:0.0001:pi;b2=-18.55;a3=pi:0.0001:5.299;b3=-18.65.*(2.3

26、./(5.8161-cos(a3)-4.5161.*(1-0.04903.*(sin(a3).2).0.5).1.4061;a4=5.299:0.0001:2.*pi;b4=-79.16;plot(a1,b1,a2,b2,a3,b3,a4,b4)蓋側(cè)活塞力軸測(cè)活塞力Matlab程序及Matlab圖：>> a1=0:0.0001:1.7926;b1=10.94.*(2.3./(-3.2161+cos(a1)+4.5161.*(1-0.04903.*(sin(a1).2).0.5).1.4061+26.26.*(2.32./(-3.1961+cos(a1)+4.5161.*(1-0.0

27、4903.*(sin(a1).2).0.5).1.4061;a2=1.7926:0.0001:pi;b2=125.58;a3=pi:0.0001:4.211;b3=46.44.*(0.3./(-3.2161+cos(a3)+4.5161.*(1-0.04903.*(sin(a3).2).0.5).1.4061+79.14.*(0.32./(-3.1961+cos(a3)+4.5161.*(1-0.04903.*(sin(a3).2).0.5).1.4061;a4=4.211:0.0001:2.*pi;b4=37.20;plot(a1,b1,a2,b2,a3,b3,a4,b4)軸測(cè)氣體力圖摩擦力

28、和往復(fù)慣性力Matlab程序及Matlab圖：>> a1=0:0.0001:2.*pi;b1=16.71;a2=0:0.0001:2.*pi;b2=14.38.*(sin(a2)+0.2214.*cos(2.*a2); plot(a1,b1,a2,b2) 摩擦力和往復(fù)慣性力圖綜合活塞力的matlab程序及matlab圖：>> a1=0:0.0001:0.9934;b1=-79.16.*(0.3./(5.8161-cos(a1)-4.5161.*(1-0.04903.*(sin(a1).2).0.5).1.4061+10.94.*(2.3./(-3.2161+cos(a1

29、)+4.5161.*(1-0.04903.*(sin(a1).2).0.5).1.4061+26.26.*(2.32./(-3.1961+cos(a1)+4.5161.*(1-0.04903.*(sin(a1).2).0.5).1.4061+16.71+14.38.*(sin(a2)+0.2214.*cos(2.*a2);a2=0.9934:0.0001:1.7926;b2=-18.55+10.94.*(2.3./(-3.2161+cos(a1)+4.5161.*(1-0.04903.*(sin(a1).2).0.5).1.4061+26.26.*(2.32./(-3.1961+cos(a1)

30、+4.5161.*(1-0.04903.*(sin(a1).2).0.5).1.4061+16.71+14.38.*(sin(a2)+0.2214.*cos(2.*a2);a3=1.7926:0.0001:pi;b3=-18.55+125.58+16.71+14.38.*(sin(a3)+0.2214.*cos(2.*a3);a4=pi:0.0001:4.211;b4=-18.65.*(2.3./(5.8161-cos(a4)-4.5161.*(1-0.04903.*(sin(a4).2).0.5).1.4061+46.44.*(0.3./(-3.2161+cos(a4)+4.5161.*(1

31、-0.04903.*(sin(a4).2).0.5).1.4061+79.14.*(0.32./(-3.1961+cos(a4)+4.5161.*(1-0.04903.*(sin(a4).2).0.5).1.4061+16.71+14.38.*(sin(a4)+0.2214.*cos(2.*a4);a5=4.211:0.0001:5.299;b5=-18.65.*(2.3./(5.8161-cos(a5)-4.5161.*(1-0.04903.*(sin(a5).2).0.5).1.4061+37.20+16.71+14.38.*(sin(a5)+0.2214.*cos(2.*a5);a6=5

32、.299:0.0001:2.*pi;b6=-78.14+37.20+16.71+14.38.*(sin(a5)+0.2214.*cos(2.*a5); plot(a1,b1,a2,b2,a3,b3,a4,b4,a5,b5,a6,b6)綜合活塞力圖綜合活塞力的最大值約在1.8弧度的位置，其大小約為135kN。十三、切向力計(jì)算及切向力圖的繪制由化工機(jī)器P72公式得切向力的公式為：式中：為列的綜合活塞力，T為切向力由以上公式和已算得的綜合活塞力，使用Matlab容易得、列的切向力圖：>> a1=0:0.0001:0.5793;b1=(-(3.2254./(5.6461-cos(a1)-4

33、.5161.*(1-0.04903.*(sin(a1).2).(1./2).1.3898+5.5263./(5.7161-cos(a1)-4.5161.*(1-0.04903.*(sin(a1).2).(1./2).1.4061)+37.7362./(-3.3861+cos(a1)+4.5161.*(1-0.04903.*(sin(a1).2).0.5).1.3898+16.34).*(sin(a1)+0.2214.*sin(2.*a1)./(2.*(1-0.04903.*(sin(a1).2).0.5);a2=0.5793:0.0001:1.8392;b2=(-35.61+37.7362./

34、(-3.3861+cos(a2)+4.5161.*(1-0.04903.*(sin(a2).2).0.5).1.3898+15.15.*(cos(a2)+0.2214.*cos(2.*a2)+16.34).*(sin(a2)+0.2214.*sin(2.*a2)./(2.*(1-0.04903.*(sin(a2).2).0.5);a3=1.8392:0.0001:pi;b3=(-35.61+55.27+15.15.*(cos(a3)+0.2214.*cos(2.*a3)+16.34).*(sin(a3)+0.2214.*sin(2.*a3)./(2.*(1-0.04903.*(sin(a3).

35、2).0.5);a4=pi:0.0001:3.9151;b4=(-(37.5274./(5.6461-cos(a4)-4.5161.*(1-0.04903.*(sin(a4).2).(1./2).1.3898+69.3051./(5.7161-cos(a4)-4.5161.*(1-0.04903.*(sin(a4).2).(1./2).1.4061)+ 3.2438./(-3.3861+cos(a4)+4.5161.*(1-0.04903.*(sin(a4).2).0.5).1.3898+15.15.*(cos(a4)+0.2214.*cos(2.*a4)+16.34).*(sin(a4)+0

36、.2214.*sin(2.*a4)./(2.*(1-0.04903.*(sin(a4).2).0.5);a5=3.9151:0.0001:5.01;b5=(-(37.5274./(5.6461-cos(a5)-4.5161.*(1-0.04903.*(sin(a5).2).(1./2).1.3898+69.3051./(5.7161-cos(a5)-4.5161.*(1-0.04903.*(sin(a5).2).(1./2).1.4061)+13.19+15.15.*(cos(a5)+0.2214.*cos(2.*a5)+16.34).*(sin(a5)+0.2214.*sin(2.*a5).

37、/(2.*(1-0.049035.*(sin(a5).2).0.5);a6=5.01:0.0001:2.*pi;b6=(-108.0774+13.19+15.15.*(cos(a6)+0.2214.*cos(2.*a6)+16.34).*(sin(a6)+0.2214.*sin(2.*a6)./(2.*(1-0.04903.*(sin(a6).2).0.5);plot(a1,b1,a2,b2,a3,b3,a4,b4,a5,b5,a6,b6)切向力圖I-IV列最大切向力出現(xiàn)在約為5.0弧度位置處，其大小約為78kN。由以上公式和已算得的綜合活塞力，使用Matlab容易得II-III列的切向力圖：

38、>> a1=0:0.0001:0.9934;b1=(-79.16.*(0.3./(5.8161-cos(a1)-4.5161.*(1-0.04903.*(sin(a1).2).0.5).1.4061+10.94.*(2.3./(-3.2161+cos(a1)+4.5161.*(1-0.04903.*(sin(a1).2).0.5).1.4061+26.26.*(2.32./(-3.1961+cos(a1)+4.5161.*(1-0.04903.*(sin(a1).2).0.5).1.4061+16.71+14.38.*(sin(a2)+0.2214.*cos(2.*a2).*(sin(a1)+0.2214.*sin(2.*a1)./(2.*(1-0.04903.*(sin(a1).2).0.5);a2=0.9934:0.00