![比较沸点法测定乙腈的蒸汽压(共11页)_第1页](http://file4.renrendoc.com/view/03c1a631f71d2a86fe8fe5f08bbfe926/03c1a631f71d2a86fe8fe5f08bbfe9261.gif)
![比较沸点法测定乙腈的蒸汽压(共11页)_第2页](http://file4.renrendoc.com/view/03c1a631f71d2a86fe8fe5f08bbfe926/03c1a631f71d2a86fe8fe5f08bbfe9262.gif)
![比较沸点法测定乙腈的蒸汽压(共11页)_第3页](http://file4.renrendoc.com/view/03c1a631f71d2a86fe8fe5f08bbfe926/03c1a631f71d2a86fe8fe5f08bbfe9263.gif)
![比较沸点法测定乙腈的蒸汽压(共11页)_第4页](http://file4.renrendoc.com/view/03c1a631f71d2a86fe8fe5f08bbfe926/03c1a631f71d2a86fe8fe5f08bbfe9264.gif)
![比较沸点法测定乙腈的蒸汽压(共11页)_第5页](http://file4.renrendoc.com/view/03c1a631f71d2a86fe8fe5f08bbfe926/03c1a631f71d2a86fe8fe5f08bbfe9265.gif)
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、 沸点(fidin)比较法测定乙腈的蒸汽压Michael B. Ewing* and Jesus C. Sanchez OchoaDepartment of Chemistry, University College London, 20 Gordon Street, London, WC1H 0AJ, United Kingdom摘 要 目前已测量的相应温度下乙腈的蒸汽压从278 K、4.3 kPa 到540 K、4455 kPa。约在536 K时乙腈开始(kish)分解,所以我们可以评估热解对蒸汽压力的影响。乙腈易潮解,常做为“干”与“湿”的评价样本,我们可以用类似的方法调查当水作为杂质时
2、的影响。这个结果(ji gu)主要由广义瓦格纳方程得来,,其中一个关键特性是使用反向对比温度=1 - T/TC,TC是临界温度。在、1.5、2.5、5标准下,将产生系统误差,对于方程ln(P/PC)=(TC/T)(c1 + c1.51.5+ c22 + c2.52.5 + c5.55.5),PC是临界压力,从291 K到535 K都适用于我们的结果,其中ln P的标准偏差为6310-6,与以前的文献相比,在三相点和临界点处显著增加了相关性的范围。通过外推到TC=545.46 K,我们可以获得临界压力PC为4835 kPa和在Tslg =229.35 K.下的三相压力167 Pa。安托万方程的结
3、果是在压力低于125 kPa且冷凝温度标准偏差为1.5 mK下获得的。前 言 纯物质的蒸汽压力可以使用沸点比较法准确、快速的测量。1-5在不同的沸点测定计中,样品和参考流体在同一气体压力下沸腾,如氦或氮,就可以测量样品的冷凝温度Ts和参考流体的冷凝温度Tr。通过已知参考流体的P、T的变化关系类推到样品的一个状态点(p,Ts),就可以得到这个压力。这种方法有许多有点:可以避免直接测量压力、通过沸腾脱出液体、在热管中沸点测定计可作为高性能的恒温器。尽管这些和静态装置的渐变相比很容易计算且变化总是较小,但相对的缺点是需要大量的测定温度、在高压下缓冲气体的溶解度、由于压头导致的热梯度。但最大的优势是测
4、量的速度,通常静态实验需要一天才能测量露点压力,而用这种方法只需一个小时就能确定一个状态点。 乙腈是一种重要的极性有机溶剂材料,其蒸汽压已被测量多次,尽管如此,公布的数据依然存在很大的差异。测量部分原材料:乙腈是由B. D. H.化学部提供的高效色谱等级材料,其摩尔百分比最小为99.9 %,最大含水量为0.02%。我们没有通过蒸馏提纯原料,因为相对于所认为的纯度影响,这部分水对实验的影响很小。如下所示,水确实是一个严重的杂质。 Burfield 等6,7研究各种干燥剂和干燥情况对有机溶剂的有效性。已知大量的超重水被添加到严格干溶剂中,通过从氚活动得到的残留水来评估干燥过程。他们发现对乙腈大多数
5、化学干燥剂在很大程度上是无效的,只有用分子筛可以得到极干燥的乙腈,其中水的摩尔分数少于10-6。因此,我们遵从Burfield等的意见使用3A分子筛,它在525 K温度下经7天(至少1天)的两级干燥得以活化。实验(shyn)装置和过程 在装载乙腈前,用高于375 K的温度加热样品沸点测定计,使其排空10小时,并用干燥的氩气冲洗。从低压(dy)至175 kPa,用体积为40 cm3的样品来做一些列的测试(csh)(在表1设置为A),并需要用205 cm3的样品在更高的对比温度下测试。第二个系列的测量(在表1设置为B)从大气压力开始,提供了和A组重叠的数据并且扩大了热分解区域的数据。乙腈是来自相同
6、的原料,只是样品是单独干燥,由于乙腈的强吸湿作用,一些污染物随着水转移到沸点测定计是不可避免的。6,7此外,在354 K、99 kPa和431 K、703 kPa之间用样品测量一组数据来观察水对实验的影响,样品由B. D. H.化学部提供。 在设备已经详细描述之前,2在这里只有一个简要描述。样品和参考流体沸点测定计通过金属冷凝器(用水冷却)连接,并放置(用固体CO2冷却)在同一压力线下,配有2 dm3的镇流器和适当的仪表和阀门,可以用来使装置排空或用干燥的氦气加压。沸点测定计本身是由不锈钢管构成,并有很好的凹角温度计,配有成对的辐射防护屏,使32cm长的铂电阻温度计和45cm长的铂电阻温度计能
7、够分别深入样品和参考流体中。实验证实2,当只浸入到15cm时,温度计的温度降低了106倍。环状加热器(最大功率为140W)被夹紧放置在圆桶式锅炉下4cm处,沸点测定计被绝热能够观察到200K的温差。但是,没有二次加热,我们依靠作为热管的具有高导热性的沸点测定计来消除温度梯度。 在沸点比较法中温度测量非常重要,因为温度是唯一被测量的实验数据。温度计是由N. P. L.在ITS-90上校准的,这个过程中它们的电阻在三相点温度下被定期测量检查。从水的冷凝温度开始,用Wagner8提供的水的蒸汽压的相关性来计算压力,在氦和水蒸气或乙腈之前做一些修正。由于会减压,所以测量时应加压来避免水和样品的交叉污染
8、。同样,在测量间歇,装置也常常放置在加压条件下。当装置在室温下4小时后获得第一个状态点,但此后,随着压力增加,平衡态在一小时后获得。结果与讨论Wagner9已研究出一组方程来说明蒸汽压的相关性。方程的一般式为 (1)其中TC、PC为临界温度和临界压力。 (2)这是反向对比温度变量。式1中的数据总和需用到的条件,其可由相关方程式确定。事实上,、1.5总是存在,标准形式由下式给出 (3)该式符合大多数结果,另一个标准形式是用3、6,其运用也很广泛。但是,当要得到大量且精确的数据时,就需要额外的条件。2,8-11临界温度很容易得到(例如,通过密封管直接观察),但临界压力不易得到,特别是对于有机液体可
9、靠的实验太少。方程式1中的ln p很易修改,而ln(p/pc)不易,在回归分析和处理数据时,ln(p/pc)作为目标函数,ln pc是一个可调变量。我们已经证明1,运用Wagner方程可以从很低的对比温度推导出可信度很高的临界压力。因此,现在我们首选的方法是用蒸汽压和临界温度通过回归分析得到临界压力PC。对于乙腈,我们也通过外推法得到了三相点的压力。我们运用Wagner方程局部和整体(zhngt)来观察两组数据的内在联系。Trejo and McLure12讨论了乙腈的临界性质,我们(w men)就运用了他们的实验数据并用进行调整(tiozhng)以适应于ITS-90。表1列出了乙腈的两种样品
10、的蒸汽压,两中样品分别干燥。A组的25个结果包含了从278K、4.3kPa到373K、175kPa,B组的结果包含了从355K、101kPa到535K、4175kPa。下面分析时用到了A组的前五个点,其温度为277.9 T/K 290.5 K;、A组T 290.5 K;、B组;、A组T 290.5 K;、湿样品(yngpn)图4由方程6计算的乙腈蒸汽压的分数偏差lnp=p/p:、在回归谈条件下的结果;、不在回归条件下的结果;、Trejo12等提供的临界点;、Kratzke等13;、Mousa14;、Francesconi等15;、Putnam等16;、Dojcansky17;、Meyer18;
11、、Hall等19;、Heim20;、Warowny22;、DiElsi等25;、Wilson26;、Smith等17-31;*、Owens32;实线、Dykyj等21;点虚线、Yaws23;虚线、DIPPR Project24图5由方程(fngchng)6计算的乙腈蒸汽压的分数偏差lnp=p/p:、在回归谈条件(tiojin)下的结果;、不在回归条件(tiojin)下的结果;、Trejo12等提供的临界点;、Kratzke等13;、Putnam等16;、Dojcansky17;、Meyer18;、Heim20;、Warowny22;、DiElsi等25;、Wilson26;、Smith等17-
12、31;*、Owens32;实线、Dykyj等21;点虚线、Yaws23;虚线、DIPPR Project24图6由方程6计算的乙腈蒸汽压的偏差p:、在回归谈条件下的结果;、不在回归条件下的结果;、Putnam等16;、Dojcansky17;、Meyer18;、Heim20;、DiElsi等25;、Wilson26;、Smith等17-31;*、Owens32;实线、Dykyj等21这个(zh ge)结果是在由样品A在低压下得到的,样品A的品质没有B高,不过在温度高于290K、在ln pc、和2的条件下,Wagner方程适用(shyng)于A组的数据且标准偏差为、。这个(zh ge)方程比一般
13、的Wanger方程少一个条件,并且计算出来的临界压力准确性很高,但是回归性不明显。在低压下,Antoine方程的计算很方便 (5)这个方程的计算结果在温度290K到362K(或压力8kPa到125kPa)范围内符合A组数据,其标准偏差为(lnp)=6110-6、(T)=1.5mK。图2显示了在相同温度范围内,用Antoine方程计算A组和B组结果的偏差。在低压下(4.3到7.2kPa)很难得到水在稳定沸腾时的那五个点的结果,因为沸点测定计的最初设计是用在高压下的,它并不适用于泡罩33。回归分析时并不取用这五个数据,因为与其他数据相比,它的可信度不高,绘制在图2-6是为了这个平滑趋势的偏差,可与
14、其他文献数据形成一个小的比较。尽管样品A和B的实验结果内部具有高度一致性,沸腾温度只相差3.7mK。用同一种方法合理处理两组数据的优点并整合两组数据是很困难的。但是Wagner方程 (6)其中TC=545.46K,有五个,它适用于数据结果,标准偏差(lnp)=6310-6、(T)=2.0mK。方程6的回归系数为PC=4835.3 kPa,、C1=-7.542771、C1.5= -0.549679、C2=4.958924、C2.5= -5.030732、C5.5= -1.699257。图3对比了方程6计算的冷凝温度偏差T,这种绘制偏差点的图对沸点测量很有用,因为沸点的不确定性应该和P和lnP大致
15、相似。在低温下,出现正弦线的系统偏差,极差大约为3.5mK,但是在高于425K、600kPa时方程6能给数据一个典型描述,并且能够推断出临界压力(58350.7)kPa。PC的不确定性是由方程6得到的标准偏差,如果将0.3K的不确定性归咎于Trejo和McLure测量的TC,那么PC将增加18kPa。我们的计算值和Trejo和McLure直接观察的结果(4.830.02) MPa能很好的吻合。虽然正弦偏差表明方程6并不完全实用,但是在相同条件下两种样品的结果能够很好的吻合,事实上,和以前的一些文献相比,2.0mK的标准偏差和4mK的最大偏差是很好的结果。我们还用没经过分子筛干燥的样品做了实验来
16、观察水对实验结果的影响。如图3所示,在大气压附近,冷凝温度大约提高了16mK,但是在700kPa时又减少到-3mK。乙腈的蒸汽压已被多次测量,结果用式6计算的ln p =p/p表示,如图4所示。有机化合物的蒸汽压不一致并不少见,但是乙腈的蒸汽压极其分散。在低于正常沸点温度下,大多数结果会低1%,Hall和 Baldt19的实验结果低了7%,Owens等32的实验结果高了3%。在室温左右,Francesconi等15的实验结果高达了大约6%(在文献中超过12%)但是比350K时高1-4%。 Mousa14的实验结果在临界状态下低了4%但在440K下达到了6%。图5显示了在一个灵敏的范围实验结果和
17、以往文献的对比,但是对于没用(mi yn)的那些数据Mousa、14 Francesconi等15、Hall和Baldt19都没有做进一步的讨论。对于剩下(shn xi)的数据都是在200kPa以上的,且是Kratze和Mller13 、Warowny22和临界点。Trejo和McLure12中非常吻合的临界点的数据都标记(bioj)在上面。由Kratze和Mller13报道的温度在353-475K的14个蒸汽压的数据能很好的我们的实验数据吻合,并由方程6得到方差为0.024%比估计精度0.1%更少。Warowny22的实验扩展到463 K 或1.3 MPa,精度低于我们(除了一个点),但是其
18、平均偏差为0.7%。在低压下分数偏差lnP=P/P相当大,因此图6给出了压力条件下的偏差。Putnam等16的实验结果一直高于我们的测量,但是其偏差只有19Pa,Meyer等18的实验结果一直低于我们的实验结果达到了96Pa。Dojcansky等17获得的蒸汽压比由式6的计算结果低,只有一个点例外,Heim等20的实验结果横跨我们的结果且方差达到了60Pa。我们有五个低于8kPa的结果也标于图6,这几个点没有用于回归分析,因为实验的冷凝温度不稳定。和文献数据相比,偏差P在0.9-2.5之间,表明我们的测量达到了严格的标准。我们已经提出了几个乙腈蒸汽压的关联性,且将他们的偏差标注在图4-6中。直
19、到正常沸点,Antoine方程的结果由温度记录器,并根据Landolt-Brnstein21的提议得到的,这是在Putnam16和Meyer18实验的基础上得到的,并表示在图5、6中。但是,他们提到,在正常沸点到530K温度之间Antoine方程并不适用,在正常沸点时它得到了高达700Pa的压力(见图5),但随后计算的压力慢慢的变低,和Francesconi等15 、Mousa14 、Warowny22的最高压一样,在539K是达到了-3.9%。虽然这些偏差的可靠性比方程要求的低,但是方程不是真的适用于对比温度大约0.75的蒸汽压。Yaws23和DIPPR24提出了相关性同样运用于方程并用相同
20、的参数。因此,计算出来的压力很相似,在图4中的比例几乎一致。如图4、5显示,在大气压下两个方程根据Kratzke和Mller13提出的建议所计算,但是在低于350K下还是比Owens等32高。在图6中由方程6计算的偏差太大,但是也出现在图4、5。蒸汽压方程向下的外推法是很严格的,可能Wanger方程尤其容易收影响因为反向温度变量在较低位温度下会增加。Frenkel等34调查了Wanger方程标准形式的外推法,并显示可能得到dp/dT的负值,就会导致得到负的汽化焓。Putnam等16报道了三相点的温度Ts+l+g =229.349 K但是没有相应的压力。DIPPR24用Putnam等提供的Ts+
21、l+g计算三相点的压力为187Pa。由方程6得到的三相点的温度为229.349K、压力为167Pa,dp/dT缓慢减少,外推法在超过60K的所有阶段中都很准确。虽然两种样品的整合结果在正常沸点时因不同的含水量的差异导致了3.7K的温差和方程6得到的系统偏差,这样的缺陷与文献中数据的分散和我们观察到的未干燥的样品的16K的温差相比是很轻微的。在高于425K的温度下,测量提出可以解决文献中存在差异的问题,并且我们推断对于乙腈蒸汽压的计算方程6是最可靠。特别地,它能平稳的外推到三相点和临界点,并且由它计算的液体蒸汽压和文献估计值相比能较好的吻合参考文献(1) Ewing, M. B.; Sanche
22、z Ochoa, J. C. An Ebulliometer for Measurements of Vapor Pressure at Low Temperatures: The Vapor Pressures and the Critical State of Perflouromethylcyclopentane. J. Chem. Thermodyn. 1998, 30, 189-198.(2) Ewing, M. B.; Sanchez Ochoa, J. C. The Vapor Pressure of Cyclohexane over the Whole Fluid Range
23、Determined Using Comparative Ebulliometry. J. Chem. Thermodyn. 2000, 32, 1157-1167.(3) Ewing, M. B.; Sanchez Ochoa, J. C. The Vapor Pressures of n-Ocatane Determined Using Comparative Ebulliometry. Fluid Phase Equilib. 2003, 210, 277-285.(4) Ambrose, D.; Ewing, M. B.; Ghiassee, N. B.; Sanchez Ochoa,
24、 J.C. The Ebulliometric Method of Vapor-Pressure Measurement:Vapor Pressures of Benzene, Hexafluorobenzene, and Naphthalene. J. Chem. Thermodyn. 1990, 22, 589-605.(5) Weber, L. A. Vapor Pressure of Heptane from the Triple Point to the Critical Point. J. Chem. Eng. Data 2000, 45, 173-176.(6) Burfield
25、, D. R.; Lee, K. H.; Smithers, R. H. Desiccant Efficiency in Solvent Drying. A Reappaisal by Application of a Novel Method for Solvent Water Assay. J. Org. Chem. 1977, 42, 3060-3065.(7) Burfield, D. R.; Gan, G.-H.; Smithers, R. H. Molecular Sieves,Desiccants of Choice. J. Appl. Chem. Biotechnol. 197
26、8, 28, 23-30.(8) Wagner, W.; Pruss, A. International Equations for the Saturation Properties of Ordinary Water Substance. Revised According tothe International Temperature Scale of 1990. J. Phys. Chem. Ref.Data 1993, 22, 783-787.(9) Wagner, W. Habilitationschrift, TU, Braunschweig, 1973; Wagner, W.
27、Forschr.sBer VDI-Z. Reihe 3, 39. 1974. Translated as a shortened version as A New Correlation Method for Thermodynamic Data Applied to the Vapour Pressures of Argon, Nitrogen and Water. IUPAC Thermodynamic Tables Centre, Imperial College, London 1977.(10) Nowak, R.; Kleinram, R.; Wagner, W. Measurem
28、ent and Correlation of the (p,F,T) Relation of Ethylene. Saturated-Liquid and Saturated-Vapor Densities and Vapor Pressures along the Entire Coexistence Curve. J. Chem. Thermodyn. 1996, 28, 1441-1460.(11) Duschek, W.; Kleinram, R.; Wagner, W. Measurement and Correlation of the (p,F,T) Relation of Ca
29、rbon Dioxide.Saturated-Liquid and Saturated-Vapor Densities and Vapour Pressures along the Entire Coexistence Curve J. Chem. Thermodyn. 1990,22, 841-864.(12) Trejo Rodriguez, A.; McLure, I. A. Critical Temperature and Pressure of Ethanenitrile. J. Chem. Thermodyn. 1979, 11, 1113-1114.(13) Kratzke, H
30、.; Muller, S. Thermodynamic Properties of Acetonitrile. 2. (p,F, T) of Saturated and Compressed Liquid Acetonitrile. J. Chem. Thermodyn. 1985, 17, 151-158.(14) Mousa, A. H. N. Vapor Pressure and Saturated-Vapor Volume of Acetonitrile. J. Chem. Thermodyn. 1981, 13, 201-202.(15) Francesconi, A. Z.; Fr
31、anck, E. U.; Lentz, H. Die PVT Daten des Acetonitrils bis 450 C und 2500 bar. Ber. BunsensGes. 1975,79, 897-901.(16) Putnam, W. E.; McEachern, D. M.; Kilpatrick, J. E. Entropy and Related Thermodynamic Properties of Acetonitrile. J. Chem. Phys.1965, 42, 749-755.(17) Dojcansky, J.; Heinrich, J. Quote
32、d by Boublik, T.; Fried, V.; Hala, E. The Vapour Pressures of Pure Substances, 2nd ed.; Elsevier: Netherlands, 1984; p 17.(18) Meyer, E. F.; Renner, T. A.; Stec, K. S. Cohesive Energies in Polar Organic Liquids. II. The n-Alkyl Nitriles and the 1-Chloroalkanes.J. Phys. Chem. 1971, 75, 642-648.(19) H
33、all, H. K.; Baldt, J. H. Thermochemistry of Strained-Ring Bridgehead Nitriles and Esters. J. Am. Chem. Soc. 1971, 93, 140-147.(20) Heim, G. Quoted by Timmermans, J. Physico-Chemical Constants of Organic Compounds; Elsevier: Netherlands, 1950; p 528.(21) Dykyj, J.; Svoboda, J.; Wilhoit, R. C.; Frenke
34、l, M.; Hall, K. R.Vapor Pressures of Chemicals, Vapor Pressure and Antoine Constants for Nitrogen Containing Organic Compounds. Landolt- Bornstein; Springer-Verlag: Berlin, 2001; Group 4, Vol. 20,Subvol. C.(22) Warowny, W. Volumetric and Phase Behavior of Acetonitrile at Temperatures from 363 to 463
35、 K. J. Chem. Eng. Data 1994, 39, 275-280.(23) Yaws, C. L. Handbook of Vapor Pressures; Gulf Publishing Company: USA, 1994; Vol. 1, p 343. .(24) The DIPPR Project 801, Physical and Thermodynamic Properties of Pure Chemicals, Evaluated Process Design Data, Part 1 (Revision date for acetonitrile, 1987)
36、: New York, 1990.(25) DiElsi, D. P.; Patel, R. B.; Abbott, M. M.; Van Ness, H. C. Excess Thermodynamic Functions for Ternary Systems. 3. Total Pressure Data and GE for Acetone-Acetonitrile-Methyl Acetate at 50 C.J. Chem. Eng. Data 1978, 23, 242-245.(26) Wilson, S. R.; Patel, R. B.; Abbott, M. M.; Va
37、n Ness, H. C. Excess Thermodynamic Functions for Ternary Systems. 4. Total Pressure Data and GE for Acetonitrile-Ethanol-Water at 50 C. J. Chem.Eng. Data 1979, 24, 130-132.(27) Maher, P. J.; Smith, B. D. Vapor-Liquid Equilibrium Data for Binary Systems of Chlorobenzene with Acetone, Acetonitrile, Et
38、hyl Acetate, Ethylbenzene, Methanol, and 1-Pentane. J. Chem. Eng. Data 1979, 24, 363-377.(28) Maher, P. J.; Smith, B. D. Vapor-Liquid Equilibrium Data for Binary Systems of Aniline with Acetone, Acetonitrile, Chlorobenzene, Methanol, and 1-Pentane. J. Chem. Eng. Data 1980, 25,61-68.(29) Khurma, J. R.; Muthu, O.;
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 电子产品物流合同要点分析
- 2025年度办公室绿植养护与室内环境美化合同
- 房屋租赁合同公文
- 企业人才测评及职业发展规划支持方案设计
- 云计算服务配置与管理手册
- 解决方案设计与实施指南
- 设计服务合同书
- 企业信息化解决方案操作手册
- 建设工程施工分包委托协议书
- 车床购买合同样本
- 农村宅基地和建房(规划许可)申请表
- (完整版)袱子的书写格式和称呼
- 供应商新增或变更申请表
- 2023年中国农业银行应急预案大全
- 【新版本】华为 H12-711 V4.0 HCIA-Security 认证华为安全题库(含答案)
- 村卫生室2023年度绩效考核评分细则(基本公共卫生服务)
- 关联公司合作合同
- 【建模教程】-地质统计学矿体建模简明教材
- 重源煤矿 矿业权价款计算书
- PSM工艺安全管理
- 7天减肥餐食谱给你最能瘦的一周减肥食谱
评论
0/150
提交评论