通过对芳烃分子3D结构的有限元分析，建立了分子结构的刚度矩阵和振动方程，用MATLAB软件计算得到相应分子结构的固有频率和矩阵特征值，用多元回归分析方法建立了液体芳烃分子导热率的QSPR模型，用线性规划法对模型指数进行优化、均值化和差值化处理. 结果表明，301个液体芳烃导热率实验数据的训练集实验值与计算值的相关系数高达0.979，平均绝对误差小于0.0024 W/(m?K)，平均相对绝对误差小于2.02%；22个液体芳烃导热率实验数据的测试集实验值和预测值的相关系数0.978，平均绝对偏差小于0.0028 W/(m?K)，平均相对绝对偏差小于2.99%.
The 3D molecular structure of aromatic hydrocarbon was analyzed by finite element analysis to establish the stiffness matrix and vibration equation, from which nature frequency and stiffness matrix eigenvalues are calculated with MATLAB software, then QSPR model for the conductivity of liquid aromatic hydrocarbons was developed by multivariate regression method, and further performing optimization, mean value and D-value treatment on the model. The results showed that for training set (301 liquid aromatic hydrocarbons), the correlation coefficient between the results calculated by QSPR model and the experimented data is 0.979, the mean absolute error is less than 0.0024 W/(m?K), the relative error is less than 2.02%, for test set (22 liquid aromatic hydrocarbons), the correlation coefficients is 0.978, the mean absolute error is less than 0.0028 W/(m?K), the relative error is less than 2.99%.
[1] NEINDRE B. L.，GARRABOS Y. and NIKRAVECH M.；Measurements of the Thermal Conductivity of Propane in the Supercritical Region；J. Chem. Eng. Data，2014, 59, 3422?3433.
[2] HIDEO W. , HIDEYUKI K.；Thermal Conductivity and Thermal Diffusivity of Twenty-Nine
Liquids: Alkenes, Cyclic (Alkanes, Alkenes, Alkadienes, Aromatics),and Deuterated Hydrocarbons；J. Chem. Eng. Data，2004, 49, 809-825.
[3] ALI Z. H., SONA R., and MOSTATA L.；Estimation of Thermal Conductivity of Ionic Liquids Using a Perceptron Neural Network；Ind. Eng. Chem. Res. 2012, 51, 9886?9893.
[4] CHEN Q. L., WU K. J., and HE C. H.；Thermal Conductivity of Ionic Liquids at Atmospheric Pressure：Database, Analysis, and Prediction Using a Topological Index Method；Ind. Eng. Chem. Res. 2014, 53, 72247232.
[5] 高 硕，曹晨忠；拓扑量子方法估算液相链烷烃导热率；物理化学学报，2006，22（12）：1478～1483.
GAO, S.，CAO，C.-Z.；A Topological-Quantum Method for the Estimation of the Thermal Conductivity of Liquid Alkanes；Acta Phys. Chim. Sin.； 2006，22（12）：1478～1483.
[6] JES U. C., TRINIDADM E.-M., MANUEL G., etc.；Thermal Conductivity of Ionic Liquids: A Pseudolattice Approach；J. Phys. Chem. C 2012, 116, 1265–1273.
[7] 张克武，张玉英，液体烃的导热率[J]，化工学报，1999，50（2）：247-252.
ZHANG，K. W.；ZHANG，Y. Y.； Extended Application of Tensg on thermal Conductivity and Molecular Structure of Liquid Hydrocarbons[J]. J. Chem. Ind. and Eng. 1999，50(2)：247-252.
[8] 万丽华,梁德青,关进安. 烃类水合物导热特性的分子动力学模拟[J]. 化工学报, 2014, 65(3): 792-796. WAN L. H.，LIANG D. Q., GUAN J. A.,Characteristic of thermal conduction in hydrocarbon hydrates using molecular dynamics method[J]. CIESC Journal, 2014, 65(3): 792-796.
[9] GREGORY W. K. , PETER C. J.；Prediction of Surface Tension, Viscosity, and Thermal Conductivity for Common Organic Solvents Using Quantitative Structure Property Relationships；J. Chem. Inf. Comput. Sci.；2001, 41： 408-418.
[10] PAUL M. M.， VIPUL S. P.，and EDWIN J. M.；Prediction and Correlation of the Thermal Conductivity of Pure Fluids and Mixtures, Including the Critical Region；Ind. Eng. Chem. Res. 2002， 41：989-999.
[11] RICHARD A. P.， MARIAL. V. R.，CARLOS A. N. etc.；Measurement and Correlation of the Thermal Conductivity of Butane from 135 K to 600 K at Pressures to 70 MPa；J. Chem. Eng. Data； 2002， 47：1263-1271.
[12] LI X.-J.，WU J.-T.，and LIU Z.-G；Thermal Conductivity of Liquid 2?Methoxyethyl Acetate, 2?Ethylhexyl Acetate, and Diethyl Succinate；J. Chem. Eng. Data；2012，57，2863?2868.
[13] 仇明华，刘万强，岳 明，刘凤萍，陈冠凡；用振动力学有限元分析方法估算多类烃的液体导热率[J]；化工学报, 2016，67(7)：2672-2678.
QIU M.H，LIU W.Q， YUE M.，LIU F. P，CHEN G. F；Estimation Conductivity of Liquid Hydrocarbons by Finite Element Analysis in Vibratuon Mechanics[J]；CIESC Journal，2016，67（7）：2672-2678.
[14] 仇明华，谢文林，刘凤萍，陈东初，一种磺酰胺类碳酸酐酶II抑制剂3D-QSAR分析新方法[J].化学学报，2010, 68（24）：2581-2589.
QIU M.H., XIE, W. L., LIU F. P., CHEN, D.C., A New Method of 3D-QSAR Analysis for Carbonic Anhydrase II Inhibitors[J]. Acta. Chim. Sinic. 2010, 68(24): 2581-2589.
[15] 刘延柱，弹性细杆的非线性力学——DNA力学模型的理论基础[M]，清华大学出版社, 北京，2006. LIU, Y.-Z. Nonlinear Mechanics of Thin Elastic Rod: Theoretical Basis of Mechanical Model of DNA, Tsinghua University Press, Beijing, 2006.
[16] 方 同，薛 璞，振动理论及应用[M]，西北工业大学出版社, 西安，p99-125，2002.
FANG T., XUE, P., Vibration Mechanics and Application, NWPU press, Xi’an, 2002.
[17] GUDZINOWICZ B. J.，CAMPBELL R. H.，ADAMS J. S., Thermal Conductivity Measurements of Complex Saturated Hydrocarbons[J]. J. Chem. Eng. Data，1964，（9）: 79-81.
[18] 张玉英，张克武；分子热力学性质手册——计算方法与最新实验数据[M]；化学工业出版社，2009，205-238.
ZHANG, Y. Y., ZHANG, K.-W. , Handbook for Thermodynamics Properties of Molecules：Calculation Method and Latest Experimental Data；Chemical Industry press, Beijing, 2009, 205-238.
[19] 马沛生；有机化合物实验物性数据手册——含碳、氢、氧、卤部分；化学工业出版社，2006,594-600.
MA，P. S.； Handbook for Experimental properties Data of organic compounds：Section of Carbon, Hydrogen, Oxygen, and Halogens[M]；Chemical Industry press, Beijing, 2006, 594-600.
[20] POLING B. E. , PRAUSNITZ J. M., O’CONNELL J. P.，The Properties of Gases and Liquids, Fifth Edition [M]. Boston, The McGraw-Hill Companies，Inc.，2001:420-436.
