活塞式压缩机设计说明

1、活塞式压缩机的设计说明班级：学号：指导老师:1.题目复算19WY-9/150型氢氦气压缩机在目前操作条件下的各级压力、排气温度、排气量、功率，作出计算示功图、切向力图、活塞力图、标明最大活塞力与切向 力，核算配用电机功率是否适当？2.19WY-9/150型氢氦气压缩机简介19WY-9/150型氢氦气压缩机是我省投产 3000吨小型化肥厂的氮氢气压缩 机，二列之间为飞轮，由电机通过三角皮带拖动。压缩机为卧式、两列、门型、 四线压缩。原料（半水煤气）经脱硫后进入I级，经I级压缩后送去变换、水洗、 碳化，碳化后为碳化气。碳化气返回II级、山、IV级压缩后去洗铜、合成。019WY制刘爪绑机几很闊3.目

2、前操作条件与有关数据（1）操作条件：吸气压力：0.15MPa（绝）排气压力：16.0MPa（绝）I级出口与II级进口压力差为P=0.09MPa吸气温度:级IIIIIIIVTs( C)30354040I级进口相对湿度=1(2)气体组成组分N2H2COCOQArCH半水煤气0.210.390.120.260.0050.00450.0105碳化气0.2370.7090.0050.0290.0030.0050.012(3)有关数据：活塞行程：S=310mm活塞杆直径d=60mm转速：n=209rpm,连杆长度l=700mmI、IV列超前II、III列90度往复运动件重量：I IV列210.9kg ；

3、II III列193.7kg级IIIIIIIV直径D34021013565相对余隙0.0650.150.160.1抽气系数10.890.890.89飞轮矩GD为471.0kgm2，配用电机额定功率：155kw。设计计算一.计算各级的行程容积I级：n(D2d4d2d2)Sn(°.34 22°.°652倔2)°310.05439m3II级:VS2- (2D240.0 仃04m 3d2)s冗4(0.21 2 20.135 2)III 级：VS3n(D2d2)sn4(0.135 20.062)0.310.00356m3IV级：、，冗fVS4D2s冗0. 065

4、20.31S440.00103m3计算各级名义压力和名义压力比已知 PSi0.15MPa Py 16MPa乂1艮2乂2耳1人2PV1120.15 0.05439 308o.4867MPaTS1VS20.01704 303丘1FS20.090.5767MPa用2Pd3£ 1£ 2PS3Ps4Pd1Ps1Pd2Pd3Ps3Ps1Vs1Ts3站3艮1 %1 Ts4Ts1 Vs40.57670.152.36730.48678.18222.36730.150.054383130.003563030.150.054393130.001033033.8454. 8643. 4562. 3

5、673MPa8. 1822 MPaPs48.18221 . 955三排气温度1. 绝热指数计算首先需算出各气体的绝热指数，查活塞式压缩机附表 1常用气体物理性质得:各气体绝热指数组分N2H2CO2COO2ArCH4绝热指数1.401.411.311.401.401.681.32注：表中所列的气体绝热指数为t 0：C P 1.013 105pa时的参数，但对于理想气体绝热指数的大小与温度的关系不大，顾可不考虑温度的影响，可把k看做是常数10.210.390.120.260.0050.00451).0105-+4+ -k 11.40-1L41 -:1.31-11.40-1140-11.68 - 1

6、L32-1混合气体的k值计算式为:ki1。对于半水煤气来说有=2.5652k=1.3898同理对碳化气有:10.2370.7090.0050.0290.0030.0050.012+ +卜+ -11.40 - 1L41 - I131-11.40-11.40-11.68-11*32-1=2.4627k=1.40612计算排气温度Td1F k 1 由式 TdTs(=)kFsk 1k 1Ts£k可得1.3898 13.845 1-3898442.07KTs1£1 k1303k2 1Ts2£2 k21.4061 1Td23084. 864 匸4061486 37Kk3 11

7、.40611Td3Ts3£3 kg3133. 456 14061447. 81Kk4 1Ts4£4 k41.4061 1Td43131. 955 1.4061379 87K级次吸气温度排气温度压力比Ts(K )Td (K )I303442.073.845II308486.374.864III313447.813.456IV313379.871.955四第一次精度计算1 各级容积系数根据过程流体机械表2-1计算膨胀指数m得:FSi1.5barm110.5(k 11)10.5(1.38981)1.1949艮24. 867bar10.75(k 21)10.75(1.40611)1

8、. 3046Fs323. 673barm310.88(k 31)10.88(1.406 11)1. 3574Ps481.822barm4k41.40612 计算工况系数C和精度BC1入 V1VS1Cpsi由BCminPs1CmaxC1Ts1入 Vi0. 86440. 054390. 04701C2Ps2入v2*20. 4867Ps120. 15C3Ps3“32. 3673Ps1T.3 30. 15CPs4G入v4%48. 1822C4Ps1Ts440. 15v1Vs13030. 64570. 017040. 039463080. 893030. 76110.003560. 046513130.

9、 893030. 93890.001030. 057383130. 89BCminCmax0. 039460. 057380.690.971对于理想气体，容积系数入V 1 a ( £ m 1)11入v11衍（m111)1 - 0.065(3.8441.194911)0.8644v21a2(m221 m; 311)1 - 0.15(4.8641.304611)0.6457入v31a3(1)1 - 0.16(3.4561.357411)0.7611入v41a4(41)1 - 0.1(1.9551.40611)0.9389精度不达标进行第一次复算五第一次复算1.计算各级进出口压力由式Psi

10、 - Psi (1C1)计算各级压力2CiPsiPs2Psi2 Ps2(i0.15MPa0.50.4867(10.04701 )0.03946)0.5333MPaPs3 2 Ps3(1 g)0.52.3673(12. 3800 MPa0.047010.046511Ps4Ps4(10.58. 1822(10.047010.057387. 4428 MPaPd1Ps20.090.6233MPaP d2Ps32 . 3800 MPap d3p s47 . 4428 MPaP d4P d41 6 MPa2计算各级压力比£1用10.62334 155%0.15£2Pd22.38004

11、.4628艮20.5333£3Pd37.44283.1272艮32.3800£4Pd4162.1497Ps47. 4428FS11.5barm10.5(k 15. 333bar1) 10.5(1.38981)1.1949m10.75(k 21) 10.75(1.40611)1.3046Pj323. 800 barms10.88(k 31) 10.88(1.40611)1.3574m3艮474. 428barm4k41.40614计算容积系数入v11a1(£1 m11)1 -0.065(4.1551.194911)0.8509入v21a2(£2讥11)1

12、-0.15(4.46281.304611)0.6779入v31a3(£3 m1£4 m41)1 -0.16(3.12721.357411)0.7894入v41a4(1)1-0.1(2.14971.40611)0.9276115计算工况系数C和精度BC1MS1，C血 汕 *，B 护3确定膨胀系数根据过程流体机械表2-1计算膨胀指数m得:Ps1siMiCmaxC2Ps2Tsi入v20. 5333Psi玉2口 20. i5C3Ps3Tsi“3s32.3800PsiTs3 3S30. i5C4Ps4Tsi入v4- 7.4428PsiTs4 40. i5BCminC0.045400.

13、88CmaxC40.05i56精度不达标进行第二次复算3030. 67790. 0i7040. 045403080.893030. 78940. 003560. 048503i30. 893030. 92760. 00i030. 05i563i30. 890. 97六第二次复算重复第一次复算步骤，列出计算结果如下表级IIIIIIIV进气压力Psi0.i50.53842.32557.06i7排气压力Pdi0.62842.32557.0i67i6压力比i4.i8934.3i933.03662.2657容积系数vi0.84940.68960.79730.92ii工况常数Q0.046200.04662

14、0.047860.04858精度贞0.95<0.97七第三次复算重复第一次复算步骤，列出计算结果如下表级IIIIIIIV1 l-t p进气压力si0.i50.53832.28526.889R VV排气压力p由0.62832.28526.889i6压力比片4.18874.24523.01462.3225i容积系数1r 10.84950.69570.79930.9179工况常数E m0.046200.047020.047150.04723精度0.98>0.97八修正排气温度由式TdTd1k 1 丁場）kk11环1 k1k 1kTs，可得3031.389814. 1887 匸389845

15、2.81Kk2 1 忘£2 k21.40611羽23084. 2452 1.4061467. 62Kk311.40611Td3人3匕3 k33133 . 0146 1.4061430. 48Kk4 11.40611Td4人4£4 k43132. 3225 j4061399. 24K级次吸气温度排气温度压力比Ts（K ）Td（K ）I303452.814.1887II308467.624.2452III313430.483.0146IV313399.242.3225九排气量计算1 容积系数v10.8495v20.6957v30.7993v40.9179确定压力系数对于第一级，

16、压力系数可取为 0.950.98，故取p1 0.97 ;对于其余各级，压力系数可取为 0.981，故取p2 p3P40.99 ;确定温度系数根据各级压力比，查过程流体机械图 2-12可得:£14.1887£24.2452心 0.94；入 t2 0.94&33.0146入 t30.95&42.3225入 t40.95确定泄露系数对于泄露系数的选取，有如下规则:对于不严密，延迟关闭的阀门：vv0.01 0.04 ；对于单作用式气缸的活塞环：vr0.01 0.05 ；对于双作用式气缸的活塞环：vr 0.003 0.015 ;对于填料vpi (0.00050.001

17、0)j , j为级次。泄露系数计算表级次IIIIIIIV气阀I0.02II0.02III0.02IV0.02活塞环I0.009II0.009III0.03IV0.03填料I0.0005II0.00100.0010III0.00150.00150.0015IV0.00200.00200.00200.0020总相对泄漏量0.03400.03350.05350.052泄露系数0.96710.96760.94920.95065综上并且将结果列表如下级次行程容积容积系数人i压力系数温度系数入i泄露系数转速n(r/min)I0.054390.84950.970.940.9671II0.017040.695

18、70.990.940.9676III0.003560.79930.990.950.9492209IV0.001030.91790.990.950.95066 计算排气量由式qvi V人i入入入小，并代入上表数据可得qvi8. 5153m?/minqv22.2310m3/minqv30.5309m3/minqv40.1767m3/min十计算功率并核算电机功率是否适当1 计算各级指示功率和总指示功率计算相对压力损失根据各级进排气压力，查过程流体机械图 2-14可得：Ps10.15MPa, 5s16.4%; pd10.6283MPa, 5d19.6%Ps20.5383MPa, 5s24.4%; P

19、d22.2852MPa, 5d26.6%Ps32.2852MPa, 5s32.6%; pd36.889MPa, 5d34.8%Ps46.889MPa, 5s42%; Pd416MPa, 5d4 3.8%S01"s1"d116.0%5 025s2 5d211.0%S035s3 5d3 7.4%5 045s45d45.8%计算实际循环指示功率ni 1由公式Nidn(1605 )psi ziV,ni1w 5)ni1sini取rtj 037匕2091.38980.97N1d(10.064)0.150.84950.054391d 601.3898097 - 11.38980.97-

20、14.1887(10.16)1.38980. 971103293. 36kWN3d209(i3.0146(10.026)2.28521.406110.074) 1.40610.79930.003561.40611. 406111 103158.69kWN4d209(1601.40611. 406111.406110.058) 1.40611103104. 74kW0.02)6. 8890.91790.001032.3225(1209N2d(1600.044)0.53830.69570 017041.40611.406111.4061 14.2452(10.11)1.406111Q3248. 3

21、4kW总指示功率 NdN1dN2dN3dN4d805. 13kW2计算轴功率取机械效率m 0.88则轴功率Nd 805.13 n m0.88914. 92kW3 计算驱动机功率取传动效率c 0.97则驱动机功率Nnc914.920.97943. 22kW考虑到实际运行当中压缩机会受到载荷波动，气体泄漏等因素影响，会引起功率的增加，驱动机因留有5 15%的储备功率，故所选用压缩机的功率为 990.38kW - 1084.70kW,取压缩机功率为1100kW而配用电机额定功率为 155kW，故不适当。十示功图绘制示功图及压力指示图是压缩机循环过程在 P-V图上的表示，从图中可反映出循环的功耗在作示

22、功图时，可将膨胀变压和压缩变压过程可看做是绝热膨胀和压缩过对于l,llll,IV级汽缸，其相关参数列表如下:级次IIIIllIV行程容积Vsi0.054390.017040.003560.00103相对余隙容积i0.0650.150.160.1余隙容积VciiVsi0.0035350.0025560.00056960.000103容积系数vi0.84950.69570.79930.9179V s1vi Vsi0.046200.011850.0028460.0009454绝热指数kTi1.38981.40611.40611.4061VVcVs0.057930.019600.004130.0011

23、3V1V2斫0.02070.007010.00190.00062V3 Vc0.0035350.0025560.00056960.000103V4VcVsVs10.011730.0077460.0012840.000188利用Matlab画出示功图，图及程序如下: 第一级：>> x仁0.01173:0.00001:0.05793;>> y1=0.15;>> x2=0.0207:0.00001:0.05793;>> y2=0.15*0.05793M.3898*x2.A(-1.3898);>> x3=0.003535:0.00001:0.0

24、207;>> y3=0.6283;>> x4=0.003535:0.00001:0.01173;>> y4=0.6283*0.003535M.3898*x4.A(-1.3898);>> Plot(x1,y1,x2,y2,x3,y3,x4,y4)第二级:>> x1=0.007646:0.00001:0.01960;>> y1=0.5383;>> x2=0.00701:0.00001:0.01960;>> y2=0.5383*0.01960A1.4061.*x2.A(-1.4061);>>

25、x3=0.002556:0.00001:0.00701;>> y3=2.2852;>> x4=0.002556:0.00001:0.007746;>> y4=2.2852*0.002556A1.4061 *x4.A(-l.4061);>> Plot(x1,y1,x2,y2,x3,y3,x4,y4)第三级：>> x仁0.001284:0.00001:0.00413;>>y仁2.2852>>x2=0.0019:0.00001:0.00413;>> y2=2.2852*0.00413A1.4061.*x2

26、.A(-1.4061);>>x3=0.0005696:0.00001:0.0019;>> y3=6.889;>> x4=0.0005696:0.00001:0.001284;>> y4=6.889*0.0005696A1.4061 *x4.A(-l.4061);>> Plot(x1,y1,x2,y2,x3,y3,x4,y4)n 10'3第四级:>> x1=0.000188:0.00001:0.00113;>>y1=6.889;>>x2=0.00062:0.00000001:0.00113;&

27、gt;> y2=6.889*0.00113A1.4061.*x2.A(-1.4061);>>x3=0.000103:0.00001:0.00062;>>y3=16;>>x4=0.000103:0.00000001:0.000188;>> y4=16*0.000103A1.4061.*x4.A(-1.4061);>>plot(x1,y1,x2,y2,x3,y3,x4,y4)卜一综合活塞力的计算综合活塞力为活塞中的气体力、往复惯性力和往复摩擦力的代数和1气体力计算确定盖侧面积AGi和轴侧面积AziA G1D：)n 22評34 20.

28、。65 2)A Z1九24d2)n22-(0.34 20.06 2)40.08796m 2A G2-D40.21 20.03464m 24A Z2n(D：4d3)n22-(0.21 20.135 2)40.02032m 2A G3A Z3n(Dl4d2)n (0.135 20.05 2)40.01149m 2A G4n D24 D400.00332m 2活塞的位移公式为：x r(1 cos ) -(1 '.12sin2)2活塞的加速度为：a r (coscos 2 ),rr1552 n2209式中：0.2214,21.89radl7006060其中l-IV列运动的角度为，11-111列

29、为，且=+90°活塞运动分析/ sA Z4以II-III 列为例，贝U盖侧对应于气缸容积的活塞当量位移：xGir(1cos )丄(112 sin2)轴侧对应于气缸容积的活塞当量位移：xZir(1cos )丄(1i12 sin2)膨胀过程：PpiPd (sc)kScX计算任意过程中的某处压力/ s S k压缩过程：PyiPs（）Sc X排气过程：Ppi Pd其中余隙容积的当量行程''I-IV列气体力计算盖侧气体力计算 膨胀过程：AgiSc1Pd1(Sd-kAG4 Pd4( 二)4 Sc4x 4Sc40.087470.62831033100.065()1.38983100

30、.0651551 cose 4. 5161(1,10.04903 si n 2e)0.0033216 1033100.1(3100.11551 cos e 4.5161(1)1.4061 10.04903 sin 2e)54.9570.13(5.6461 cos e 4.516110.04903sin)1.38982e53.12(5.7161 cos0.2e 4.516110.04903sin-)14061 kN2e吸气过程:(A G1 ps1A G4 Ps4)(0.087470.151030.003326.7728103)35.61kN压缩过程:Ag1Sc1-)k1Ag4X1Ps4(学Sc4

31、x4)k40.087470.151033100.0653100.003326. 889103(-_ - -) 1.389823100.065155 1 cose 4.5161(110.04903 sin 2B)3100.131013.121(3100.1155 1 cos e 4.5161(12.13)1.4061 10.04903 sin 2B)22. 871()1.38985.6461 cos e 4.5161 10.04903sin 2e2.2()1.4061 kN5.7161 cos e 4.5161 10.04903sin 2eFAZ1Ps1( J $ )kSc1x10.087960

32、.15103(3100.0651551ps4(Sc4Sc4s)k4x43100.0653101.3898cos 04.5161(110.04903sin 20)13.194(-3.3861cos 02.134.5161120.04903sin 20)1.3898kN排气过程：FAZ1 pd1AZ4pd40.087960.628310355.27kN膨胀过程：FAZ1Pd1(&1)k1AZ4 Pd4(Sc4)k4Sc1x1Sd4x40.087960.62831033100.065(3100.0651551cos 04.5161(1”)1.3898 10.04903 sin 20)55.

33、27(0.13)1.3898 kN/ 2-3.3861 cos04.5161 d0.04903sin 20吸气过程F AZ1 ps1 AZ4 ps4 0.08796 0.15 10313.19kNII-III 列气体力计算排气过程：F(AG1 pd1AG4 Pd4)(0.087470.62831030.0033216103)108. 0774kN轴侧气体力计算压缩过程:盖侧气体力计算膨胀过程:AG2Pd2(Sc2Sc2x；)k2AG3陶衆旳0.034642.28521033100.15)1.4061 3100.151551 cos a 4.5161(1.10.04903 sin 2a)79.1

34、6(0.32)1.4061 kN5.8161 cos a 4.5161.10.04903sin 2a吸气过程F (Ag2 Ps2 Ag3 Ps3)0.03464 0.5354 10318.55kN压缩过程Sr2s k2Sc3s k3F Ag2 Ps2(笑)k2Ag3 Ps3( ®)心Sc2x2Sc3x30.034640.5383103(3100.15310)1.40613100.151551 cos a 4.5161(1 门 0.04903 sin 2a)18.65(5.8161cos a2.34.5161 <1)1.4061 kN0.04903sin 2a排气过程：F (A

35、G2 pd2 AG3 Pd3)0.034642.2852 10379. 16kN轴侧气体力计算压缩过程:Sd2s )k2X2Sc3S k,AZ3Ps3( C3f3Sc3X30.020320.53831033100.153100.011492 2852103(v v. 一v)1.4061310 0.15155 1 cos 4.5161 (1, 10.04903 sin2 )310 0.16310(3100.16155 1 cos 4.5161(1-.1)1.40610.04903 sin 2 )10. 942.30(-3.2161 cos 4.5161 10.04903sin2)1.406126

36、.26(-3.19612.32)14061kN2cos 4.5161 10.04903si n2排气过程FAZ2Pd2AZ3pd30.020322. 2852 1030.011496. 889103125. 58kN膨胀过程:FAZ2 Pd2(02ST-X2Azs Pd3(Sc)k3Sc3X30.02032 2.2852103310 0.150.01149 6.889103(亠二一)1.4061310 0.15155 1 cos 4.5161 (1. 10.04903 sin2 )310 0.16(310 0.16155 1 cos 4.5161(1)1.406110.04903 sin2 )

37、0.30.3279.14(-3.1961 cos1.406146.44( )-3.2161 cos 4.5161.10.04903sin2)1.4061 kN4.5161 10.04903s in2吸气过程：FAZ2 ps2AZ3 ps30.020320.53831030.011492. 2852 10337. 20kN2往复惯性力的计算ms1210.9kg,ms2193.7kg对l-IV列F|s1ns1rcos6 入 cos29)2109 0155 (209 n、cos60.2214cos29)60 )151$ cosB 0.2214cos29 )kN其中一阶往复惯性力:I 22092冗 2

38、Flms1r w2 cos 9210.9 0.155 () cos 9 = 15. 66 cos060二阶往复惯性力：II 2F1ms1r w cos 22092冗 2=0. 2214210.90.155()2 cos 060=3.47 cos 0对II-III 列2209 n 2F|s2ms2r wf cos 入 cos2 )1937 0.155 ( )f cos 02214cos2 )6014.38 cos02214cos2 )kN一阶往复惯性力：1 22092冗 2F2ms1r w2cos0193. 7 0.155 ()2 cos =14.38 cos2 s 1 60二阶往复惯性力：f2

39、Ims1r w 2 cos 22092冗 2=0. 2214 193. 7 0.155 ()2 cos 260=3.18 cos 2因为一四列超前二三列90°,所以&二"十90一阶往复惯性力合力：斤! = F + f2 = 14.38 C0S-i5T66sind二阶往复惯性力合力：F|!F1II +F2I =3.18 cos 2 3.47 cos 2 3往复摩擦力的计算1 P(1)取Ffs0.65 ns/30 ，式中P为第i列指示功率；m为压缩机的机械效率；s为活塞行程，n为压缩机转速对I-IV列Ffsl(Ng0.65 -1N4d)(1)nmns/30(293.36

40、0.65104 74)( 0.881)2090.31/3016. 34kN对II-III 列(N 2d0.65 -1N3d)(1)n mns/300.65(248.34158.69)(0.881)2090.31/3016.71kN十二活塞力图的绘制1)由以上的计算结果，利用 matlab绘制I-IV列活塞力图及程序如下: 盖侧活塞力程序及Matlab图：>> a仁0:0.0001:0.5793; b1=-(3.2254./(5.6461-cos(a1)-4.5161.*(1-0.04903.*(si n( a1).A2).A(1./2).M.3898+5.5263./(5.7161

41、-cos(a1)-4.5161.*(1-0.04903.*(si n(a1). A2).A(1./2).A1.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.*(si n(&3).八2).八 (1./2).A1.3898+69.3051./(5.7161-cos(a3)-4.5161.*(1-0.04903.*(s in (a3) ).A2).A(1./2).A1.4061);a4=5.01:0.0001:2.*pi;b4=-10

42、8.0774;Plot(a1,b1,a2,b2,a3,b3,a4,b4)盖侧活塞力轴测活塞力Matlab程序及Matlab图：>> a1=0:0.0001:1.84;b1=37.7362./(-3.3861+cos(a1)+4.5161.*(1-0.04903.*(si n( a1).A2).A0 .5)".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.*(s in 3)八2)八0. 5).A1.3898;a4=3.915:

43、0.0001:2.*pi;b4=13.19;Plot(a1,b1,a2,b2,a3,b3,a4,b4)轴测气体力图摩擦力和往复惯性力 Matlab程序及Matlab图: 对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)摩擦力和往复惯性力图综合活塞力的matlab程序及matlab图：>> a仁0:0.0001:0.5793;b1=-(3.2254./(5.6461-cos(a1)-4.5161.*(1-0

44、.04903.*(si n( a1).A2).A(1. /2).M.3898+5.5263./(5.7161-cos(a1)-4.5161.*(1-0.04903.*(si n(a1). A2)A(l./2)Ai.4061)+37.7362./(-3.3861+cos(a1)+4.5161.*(1-0.04903.*( sin( a1).A2).A0.5).A1.3898+16.34;a2=0.5793:0.0001:1.8392; b2=-35.61+37.7362./(-3.3861+cos(a2)+4.5161.*(1-0.04903.*(s in (a2) .A2).A0.5).A1.

45、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.2214.*cos(2.*a3)+16.34;a4=pi:0.0001:3.9151;b4=-(37.5274./(5.6461-cos(a4)-4.5161.*(1-0.04903.*(s in 4)八2)八 (1./2).A1.3898+69.3051./(5.7161-cos(a4)-4.5161.*(1-0.04903.*(s in (a4) ).A2).A(1./2).A1.4061)

46、+ 3.2438./(-3.3861+cos(a4)+4.5161.*(1-0.04903.*(s in 4)八2)八0.5)八1.3 898+15.15.*(cos(a4)+0.2214.*cos(2.*a4)+16.34;a5=3.9151:0.0001:5.01;b5=-(37.5274./(5.6461-cos(a5)-4.5161.*(1-0.04903.*(si n(&5)八2)八 (1./2).A1.3898+69.3051./(5.7161-cos(a5)-4.5161.*(1-0.04903.*(s in (a5) ).A2).A(1./2).A1.4061)+13.

47、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;Plot(a1,b1,a2,b2,a3,b3,a4,b4,a5,b5,a6,b6)综合活塞力图综合活塞力的最大值约在1.8弧度的位置，其大小约为25kN2)由以上的计算结果，利用 matlab绘制II-III列活塞力图及程序如下:盖侧活塞力程序及Matlab图：>> a1=0:0.0001:0.9934;b1=-79.16.*(0.3./(5.8161-cos(a1)-4.5161.*(1-0.04903.*(si n( a1).A2 ).A0.5).A1.40