基于Monte Carlo法的苯乙烯-丙烯腈单链模拟

doi:10.16085/j.issn.1000-6613.2023-1889

化工进展 ›› 2024, Vol. 43 ›› Issue (11): 6068-6076.DOI: 10.16085/j.issn.1000-6613.2023-1889

• 化工过程与装备 • 上一篇

基于Monte Carlo法的苯乙烯-丙烯腈单链模拟

徐燕红¹(), 陈钰¹, 袁向前¹, 刘纪昌¹^,², 赵基钢¹^,²()

^1.华东理工大学化工学院绿色能源化工国际联合研究中心，上海 200237
^2.石河子大学化学化工学院新疆兵团化工绿色过程重点实验室，新疆石河子 832003

收稿日期:2023-10-26 修回日期:2023-11-26 出版日期:2024-11-15 发布日期:2024-12-07
通讯作者: 赵基钢
作者简介:徐燕红（1999－），女，硕士研究生，研究方向为聚合物构效关系。E-mail：xu_y_h@163.com。
基金资助:
第七师石河子大学科技创新专项(QS2022001)

Monte Carlo simulation of styrene-acrylonitrile single chain

XU Yanhong¹(), CHEN Yu¹, YUAN Xiangqian¹, LIU Jichang¹^,², ZHAO Jigang¹^,²()

^1.International Joint Research Center for Green Energy Chemical Engineering, School of Chemical Engineering, East China University of Science and Technology, Shanghai 200237, China
^2.Key Laboratory for Green Processing of Chemical Engineering of Xinjiang Bingtuan, School of Chemistry and Chemical Engineering, Shihezi University, Shihezi 832003, Xinjiang, China

Received:2023-10-26 Revised:2023-11-26 Online:2024-11-15 Published:2024-12-07
Contact: ZHAO Jigang

摘要/Abstract

摘要：

聚合物链组成和分布对聚合物的性能具有十分重要的意义。基于此，本文参考某企业的单体配方，首先建立了二元共聚链增长模型，并采用Monte Carlo法对苯乙烯-丙烯腈共聚单链的链增长过程进行了模拟，分析对比了单体比对转化率及丙烯腈结合率的影响，考察了分批加料对链结构均匀性的影响。结果表明：当转化率高于91%时，丙烯腈结合率的变化率大于0.5%，造成共聚物分布不均；在考察的单体比范围内，原料中丙烯腈含量越低（大于25%），展现出来的共聚产物越均匀，可调节的转化率范围就越宽；当原料中丙烯腈为27%~35%，转化率在85.0%以上时，丙烯腈瞬时结合率的变化率仅在±0.5%的范围内波动，且丙烯腈平均结合率的变化小于1%；分批进料能有效提高共聚物的链结构均匀性，当单体比为67/33的转化率由91%提升至98%时，丙烯腈瞬时结合率的变化率小于±0.5%，且丙烯腈平均结合率的变化小于1%。上述研究结果为其他类似体系研究及工业生产稳定产品提供了一定的数据基础和理论指导。

关键词: 聚合物多元醇, Monte Carlo法, 单体比, 二元共聚, 单链模拟

Abstract:

The composition and distribution of polymer chains are of paramount importance for the performance of polymers. To address this, this paper, based on a monomer formulation from a certain company, firstly established a binary copolymer chain growth model. Employing the Monte Carlo method, we simulated the chain growth process of styrene-acrylonitrile copolymer single chains and analyzed how the monomer ratio affected the conversion rate and acrylonitrile binding ratio. Additionally, we investigated the impact of staged feedstock addition on product structure uniformity. The results reveal that, when the conversion rate exceeded 91%, the variation in acrylonitrile binding rate exceeded 0.5%, leading to an uneven distribution of the copolymer. Within the simulately studied range, a lower acrylonitrile content in the feedstock (greater than 25%) results in a more uniform copolymer product, widening the adjustable conversion rate range. When the acrylonitrile content in the feedstock ranged from 27%—35% and the conversion rate exceeded 85.0%, the instantaneous binding rate of acrylonitrile fluctuated within a range of only ±0.5%, and the average binding rate of acrylonitrile exhibited variations of less than 1%. Furthermore, a staged feedstock approach proved effective in enhancing the sequence uniformity of the copolymer. When the monomer ratio shifted from 67/33 with a conversion rate of 91% to 98%, the instantaneous acrylonitrile binding rate fluctuated by less than ±0.5%, and the average binding rate of acrylonitrile exhibited variations of less than 1%. These findings provided a foundational dataset and theoretical guidance for similar system studies and the stable production of industrial products.

Key words: polymer polyether polyols, Monte Carlo method, monomer ratio, binary copolymerization, single-chain simulation

中图分类号:

TQ316

徐燕红, 陈钰, 袁向前, 刘纪昌, 赵基钢. 基于Monte Carlo法的苯乙烯-丙烯腈单链模拟[J]. 化工进展, 2024, 43(11): 6068-6076.

XU Yanhong, CHEN Yu, YUAN Xiangqian, LIU Jichang, ZHAO Jigang. Monte Carlo simulation of styrene-acrylonitrile single chain[J]. Chemical Industry and Engineering Progress, 2024, 43(11): 6068-6076.

图/表 13

图1 苯乙烯-丙烯腈共聚计算机模拟流程图

图2 R2与C的关系

图3 F2与C的关系

图4 dF2dC与C的关系

图5 f2与C的关系

图6 单体比对R2与C的影响

图7 单体比对F2与C的影响

表1 单体比对C分别为5%、65%和95%时的F2及其差值的影响

原料中丙烯腈质量分数/%	F₂/%			F₂的差值Δ/%
原料中丙烯腈质量分数/%	转化率 5% (C₅)	转化率 65% (C₆₅)	转化率 95% (C₉₅)	Δ₁(C₆₅-C₅)	Δ₂(C₉₅-C₆₅)
45	31.8	33.2	100.0	1.4	66.8
40	30.9	32.4	99.2	1.5	66.8
35	29.9	31.1	74.4	1.2	43.3
33	29.7	31.1	47.8	1.4	16.7
31	29.9	28.7	32.1	-0.4	3.4
30	29.6	29.1	31.4	-0.5	1.3
27	28.1	26.9	25.4	-1.2	-1.5
25	27.3	26.0	18.0	-1.3	-8.0

图8 单体比对dF2dC与C的影响

表2 单体比对dF2dC的峰值和C的影响

原料中丙烯腈质量分数/%	开始突变时的C/%	达到峰值时的C/%	达到峰值时的 $d F 2 d C$ /%
45	69	85	4.5
40	79	91	6.2
35	87	95	16.8
33	91	96	22.1
31	95	98	33.1
30	95.5	99	52.7
27	96.5	100	-29.9
25	91	98.5	-5.7

表2 单体比对dF2dC的峰值和C的影响

原料中丙烯腈质量分数/%	开始突变时的C/%	达到峰值时的C/%	达到峰值时的 $d F 2 d C$ /%
45	69	85	4.5
40	79	91	6.2
35	87	95	16.8
33	91	96	22.1
31	95	98	33.1
30	95.5	99	52.7
27	96.5	100	-29.9
25	91	98.5	-5.7

图9 单体比对f2与C的影响

表3 单体比对f2变化1%时的C及其差值的影响

原料中丙烯腈质量分数/%	转化率为 3.5%时的f₂/%	f₂变化 1%后/%	f₂变化1%时的转化率/%	C差值Δ/%
45	31.2	32.2	69.4	65.9
40	30.5	31.5	77.0	73.5
35	29.6	30.6	87.7	84.2
33	29.2	30.2	94.0	90.5
31	28.6	29.6	97.4	93.9
30	28.3	29.3	98.9	95.4
27	27.5	26.5	100.0	96.5
25	26.9	25.9	93.8	90.3

图10 分批加料对丙烯腈结合率和转化率的影响

参考文献 20

1	安泽胜, 陈昶乐, 何军坡, 等. 中国高分子合成化学的研究与发展动态[J]. 高分子学报, 2019, 50(10): 1083-1132.
	AN Zesheng, CHEN Changle, HE Junpo, et al. Research and development of polymer synthetic chemistry in China[J]. Acta Polymerica Sinica, 2019, 50(10): 1083-1132.
2	李伯耿, 罗英武, 刘平伟. 聚合物产品工程研究内容与方法的思考[J]. 化工进展, 2023, 42(8): 3905-3909.
	LI Bogeng, LUO Yingwu, LIU Pingwei. Consideration on research content and method of polymer product engineering[J]. Chemical Industry and Engineering Progress, 2023, 42(8): 3905-3909.
3	赵继永, 王志鹏, 程世婧, 等. 高性能聚合物泡沫材料的制备性能与应用研究进展[J]. 高分子材料科学与工程, 2020, 36(6): 136-144.
	ZHAO Jiyong, WANG Zhipeng, CHENG Shijing, et al. Preparation, properties and applications of high-performance polymeric foam materials[J]. Polymer Materials Science & Engineering, 2020, 36(6): 136-144.
4	胡方圆, 刘冬明, 李佳乐, 等. 高性能聚合物在新型储能领域的应用进展[J]. 中国材料进展, 2019, 38(10): 990-998, 1016.
	HU Fangyuan, LIU Dongming, LI Jiale, et al. Advances in application of high-performance polymers in the new energy storage field[J]. Materials China, 2019, 38(10): 990-998, 1016.
5	余胜泉, 徐刘杰. 大数据时代的教育计算实验研究[J]. 电化教育研究, 2019, 40(1): 17-24.
	YU Shengquan, XU Liujie. Education computational experiment research in the era of big data[J]. e-Education research, 2019, 40(1): 17-24.
6	陆书来, 张欣颖, 东升魁, 等. 提高转化率对乳聚丁苯橡胶性能的影响[J]. 合成橡胶工业, 2006, 29(3): 170-173.
	LU Shulai, ZHANG Xinying, DONG Shengkui, et al. Effect of increasing conversion on performance of emulsion polymerized styrene butadiene rubber[J]. China Synthetic Rubber Industry, 2006, 29(3): 170-173.
7	孟庆民. 苯乙烯-丙烯腈梯度共聚物序列结构的¹³C-NMR表征及动态Monte Carlo模拟[D]. 上海: 华东师范大学, 2006.
	MENG Qingmin. Characterization by ¹³C-NMR and Monte Carlo simulation on sequence structure of styrene-acrylonitrile gradient copolymers[D]. Shanghai: East China Normal University, 2006.
8	刘春, 张国杰, 王衍, 等. 碳纳米管表面接枝聚合物动力学的粗粒化模拟研究[J]. 中国科学: 化学, 2023, 53(4): 778-788.
	LIU Chun, ZHANG Guojie, WANG Yan, et al. Coarse-grained simulation study of polymer chain-grafted carbon nanotubes prepared via grafting-from strategy[J]. Scientia Sinica Chimica, 2023, 53(4): 778-788.
9	刘睿, 陈曦. 聚合反应过程中的蒙特卡洛方法[J]. 高校化学工程学报, 2021, 35(3): 389-399.
	LIU Rui, CHEN Xi. Review on the Monte Carlo method in polymerization processes[J]. Journal of Chemical Engineering of Chinese Universities, 2021, 35(3): 389-399.
10	NASRESFAHANI Amin, IDOWU Loretta A, HUTCHINSON Robin A. Extractable content of functional acrylic resins produced by radical copolymerization: A comparison of experiment and stochastic simulation[J]. Chemical Engineering Journal, 2019, 378: 122087-122099.
11	PEIKERT Paul, PFLUG Kristina M, BUSCH Markus. Modeling of high-pressure ethene homo- and copolymerization[J]. Chemie Ingenieur Technik, 2019, 91(5): 673-677.
12	REGO Artur S C, BRANDÃO Amanda L T. Parameter estimation and kinetic Monte Carlo simulation of styrene and n-butyl acrylate copolymerization through ATRP[J]. Industrial & Engineering Chemistry Research, 2021, 60(23): 8396-8408.
13	CAVALLI Federica, DE KEER Lies, HUBER Birgit, et al. A kinetic study on the para-fluoro-thiol reaction in view of its use in materials design[J]. Polymer Chemistry, 2019, 10(22): 2781-2791.
14	AGARWAL Umang, DEN DOELDER Jaap. Integrated kinetic Monte Carlo-structure-rheology model for solution copolymerization of ethylene and α-olefins[J]. Macromolecules, 2019, 52(1): 292-310.
15	李黎. 苯乙烯-丙烯腈在聚醚介质中的接枝共聚合[D]. 杭州: 浙江大学, 2012.
	LI Li. Graft copolymerization of stvrene and acrylonitrile in the presence of polyether polyol[D]. Hangzhou: Zhejiang University, 2012.
16	金彦江, 沈本贤, 赵基钢, 等. 乳聚丁苯橡胶中苯乙烯结合量与单体转化率的Monte Carlo模拟[J]. 合成橡胶工业, 2012, 35(4): 268-271.
	JIN Yanjiang, SHEN Benxian, ZHAO Jigang, et al. Monte Carlo simulation for relation between bound styrene content and monomer conversion in synthesis of styrene-butadiene rubber by emulsion polymerization[J]. China Synthetic Rubber Industry, 2012, 35(4): 268-271.
17	赵晓燕. 混合同余法产生随机噪声的FPGA实现[J]. 电子科技, 2016, 29(5): 49-50, 54.
	ZHAO Xiaoyan. FPGA implementation of random noise generation by mix congruence method[J]. Electronic Science and Technology, 2016, 29(5): 49-50, 54.
18	兰胜坤. 基于MATLAB的直接序列扩频通信系统的仿真模型[J]. 信息系统工程, 2023(10): 24-27.
	LAN Shengkun. Simulation model of direct sequence spread spectrum communication system based on MATLAB[J]. China CIO News, 2023(10): 24-27.
19	张大伟, 邵英海, 左垒. 基于线性同余法的伪随机数产生算法[J]. 辽东学院学报(自然科学版), 2018, 25(3): 201-206.
	ZHANG Dawei, SHAO Yinghai, ZUO Lei. Pseudo random number generation algorithm based on linear congruence method[J]. Journal of Eastern Liaoning University (Natural Science Edition), 2018, 25(3): 201-206.
20	龚纯, 王正林. MATLAB语言常用算法程序集[M]. 2版. 北京: 电子工业出版社, 2011: 367.
	GONG Chun, WANG Zhenglin. MATLAB language commonly used algorithm assembly[M]. 2nd ed. Beijing: Publishing House of Electronics Industry, 2011: 367.

基于Monte Carlo法的苯乙烯-丙烯腈单链模拟

Monte Carlo simulation of styrene-acrylonitrile single chain

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

图/表 13

参考文献 20

相关文章 1

编辑推荐

Metrics

本文评价