本发明专利技术公开了一种自由基聚合物分子量分布操作条件优化方法。本发明专利技术对于自由基聚合模型,以给定的分子量分布曲线作为产品指标,在满足产品指标的前提下,求取使得反应物的转化率最大的反应温度条件和反应时间。该方法在目标分子量分布曲线上取特定的点,取顶点和浓度大小为顶点一半的点对应的聚合物浓度值作为优化命题中的约束,通过控制变量离散化方法得到优化结果,以结果为初值模拟计算分子量分布,与给定的分子量分布曲线比较误差,若最大误差超出5%,则可以通过再次随机给初值方式进行优化,使得最终得到的优化结果满足最大误差在5%之内。本发明专利技术有益处在于实现聚合物产品操作条件优化,可应用于离线优化和产品设计。
【技术实现步骤摘要】
本专利技术涉及高分子自由基聚合分子量分布的操作条件优化方法,尤其涉及一种自 由基聚合物分子量分布操作条件优化方法。
技术介绍
自由基聚合物产品的使用性能和加工性能由聚合物的分子量分布曲线决定,反应 物的利用率主要由反应物转化率指标决定。现有的自由基聚合机理模型的模拟可以很好 吻合实际的反应过程。在模拟优化中,通过给定分子量分布曲线可以确定最终产品的使用 性能及加工性能,此时根据最大化反应物转化率的目标来得到符合给定的分子量分布曲线 的反应操作条件,可以为实际的工业生产提供指导。由于自由基聚合反应模型维度巨大( IO5 口 IO6 ),模拟困难,通过引入序贯法(见《自由基聚合产品设计及分子量分布大规模计算 研究》冯剑一(硕士毕业论文)),实现全链长范围的分子量分布模拟,这为给定分子量分 布条件下的操作条件优化提供了可能,在工业生产中,对于给定分子量分布的操作条件优 化在产品设计,离线优化方面具备广阔的前景。动态优化一般方法有同时联立法,序贯求解法,变分法。对于超大规模模型的动态 优化求解策略,采取控制变量离散化方法,将时间域上的随时间变化的变量值离散为时间 区间的代数表达式,从而将微分变量离散为时间域上的代数变量,求解变换后的非线性代 数优化问题。需要自由基聚合反应模型提供模拟的结果作为优化命题的目标和约束。自由基聚合反应模型由规模为IO5CHO6的微分代数方程组描述。其中包含了反应 速率系数动态变化模型,采用同时联立求解法例如五步后向差分法求解会出现物理内存过 小从而使得计算的最大链长值过小(Dell 2G机器求解最大链长为7327),无法表征完整的 分子量分布曲线。通过变步长多尺度变量解耦序贯方法(见《自由基聚合产品设计及分子量 分布大规模计算研究》冯剑——(硕士毕业论文)),实现基于自由基聚合反应模型的分子量 分布模拟,链长可以达到IO5 ,描述整个分子量分布状况。从而可以为优化提供分子量分布 的信息。
技术实现思路
本专利技术的目的是针对自由基聚合反应分子量分布操作条件优化的应用场景,提供 一种。是在已知聚合物分子量分布曲 线_条件下,选取曲线中的特征点聚合物分子量浓度最大的链长点,聚合物分子量 浓度为最大聚合物分子量浓度一半的链长点,以这三个特征点的聚合物分子量浓度值 MWDiiI = 1,2,3作为约束,以单体转化率X作为优化目标,求解符合聚合物分子量分布曲线的操作条件反应时间和反应温度;给定一个随机初始温度~和反应时间&作为控制变量,通过控制变量的离散化,调用序列二次规划算法求解优化命题;自由基聚合反应模型通 过解耦方式分解为小规模自由基聚合反应模型和大规模自由基聚合反应模型,通过变步长变 阶次后向差分方法求解小规模自由基聚合反应模型,得到单体转化率X和反应速率系数在时 间离散点的数值,单体转化率X作为优化目标值;将各离散时间点的反应速率系数代入自由 基聚合反应模型,变求解自由基聚合反应模型为大规模自由基聚合反应模型,通过序贯法求 解大规模自由基聚合反应模型得到聚合物分子量分布曲线,特征点的分子量浓度值作为约束方程值;目标函数和约束方程导数值由差动方式提供,设置扰动量?;厶,以差分形式求出 导数值,序列二次规划算法根据目标函数和约束方程函数值及导数值决定下一步寻优方向及步长,并判断下一步迭代点是否满足收敛条件,若收敛,返回最优解(7/々),以最优解(7/4)作为初值计算自由基聚合反应模型聚合物分子量分布曲线权利要求1. 一种,其特征在于在已知聚合物分子 量分布曲线胃条件下,选取曲线中的特征点聚合物分子量浓度最大的链长点,聚合物 分子量浓度为最大聚合物分子量浓度一半的链长点,以这三个特征点的聚合物分子量浓度 值= 1, 3作为约束,以单体转化率X作为优化目标,求解符合聚合物分子量分布曲线的操作条件反应时间和反应温度;给定一个随机初始温度,和反应时间6作为控制变量,通过控制变量的离散化,调用序列二次规划算法求解优化命题;自由基聚合反应模型通 过解耦方式分解为小规模自由基聚合反应模型和大规模自由基聚合反应模型,通过变步长 变阶次后向差分方法求解小规模自由基聚合反应模型,得到单体转化率X和反应速率系数 在时间离散点的数值,单体转化率X作为优化目标值;将各离散时间点的反应速率系数代 入自由基聚合反应模型,变求解自由基聚合反应模型为大规模自由基聚合反应模型,通过 序贯法求解大规模自由基聚合反应模型得到聚合物分子量分布曲线,特征点的分子量浓度值作为约束方程值;目标函数和约束方程导数值由差动方式提供,设置扰动量?IA,以差分形式求出导数值,序列二次规划算法根据目标函数和约束方程函数值及导数值决定下一步寻优方向及步长,并判断下一步迭代点是否满足收敛条件,若收敛,返回最优解而),以最优解(7>而)作为初值计算自由基聚合反应模型聚合物分子量分布曲线MfD ,并计算与给定的聚合物分子量分布曲线_误差分布曲线b I,2,···,Ctsx ,若max^},! =1,2"·,C^ffi在-5%至5%之间,则认为最终的产品符合给定聚合物性能要求,否则产生新的一组初值重新开始优化,直至最后得到的误差曲线均在-5% 至5%之间;若求解不收敛,返回改变初值重新优化,得到优化后的操作条件; 方法具体步骤如下1)读取聚合物分子量分布曲线_数据,找出聚合物分子量浓度最大的链长点‘,得到链长点 对应的聚合物分子量浓度;找出聚合物分子量浓度为链长&对应的聚合物分子量浓度—半的链长点,得到链长点‘4·对应的聚合物分子量浓度MWDiMm.,以这三个链长点的聚合物分子量浓度作为优化命题的不等式约束,令-力 反应温度,t为反应时间,使满足2.根据权利要求1所述的一种, 其特征在于所述的聚合物为由单一单体连锁聚合反应,以自由基为活性中心生成的聚合物。3.根据权利要求1所述的一种,其特征在 于所述的变求解自由基聚合反应模型为大规模自由基聚合反应模型为4.根据权利要求1所述的一种,其特征 在于所述的通过序贯法求解大规模自由基聚合反应模型得到聚合物分子量分布曲线为大规模自由基聚合反应模型未知量为自由基《4 = 2,···,C^ax,和聚合物巧3 = -,C^mx, ,以乓为一组计算单元,先计算茑,马,按链长增长方向依次计算茑,召,《,乓直至全文摘要本专利技术公开了一种。本专利技术对于自由基聚合模型,以给定的分子量分布曲线作为产品指标,在满足产品指标的前提下,求取使得反应物的转化率最大的反应温度条件和反应时间。该方法在目标分子量分布曲线上取特定的点,取顶点和浓度大小为顶点一半的点对应的聚合物浓度值作为优化命题中的约束,通过控制变量离散化方法得到优化结果,以结果为初值模拟计算分子量分布,与给定的分子量分布曲线比较误差,若最大误差超出5%,则可以通过再次随机给初值方式进行优化,使得最终得到的优化结果满足最大误差在5%之内。本专利技术有益处在于实现聚合物产品操作条件优化,可应用于离线优化和产品设计。文档编号G06F19/00GK102142060SQ20111000049公开日2011年8月3日 申请日期2011年1月4日 优先权日2011年1月4日专利技术者姚臻, 范伟伟, 邵之江, 钱积新, 陈曦 申请人:浙江大学本文档来自技高网...
【技术保护点】
1.一种自由基聚合物分子量分布操作条件优化方法,其特征在于:在已知聚合物分子量分布曲线 条件下,选取曲线中的特征点:聚合物分子量浓度最大的链长点,聚合物分子量浓度为最大聚合物分子量浓度一半的链长点,以这三个特征点的聚合物分子量浓度值作为约束,以单体转化率X作为优化目标,求解符合聚合物分子量分布曲线的操作条件:反应时间和反应温度;给定一个随机初始温度和反应时间作为控制变量,通过控制变量的离散化,调用序列二次规划算法求解优化命题;自由基聚合反应模型通过解耦方式分解为小规模自由基聚合反应模型和大规模自由基聚合反应模型,通过变步长变阶次后向差分方法求解小规模自由基聚合反应模型,得到单体转化率X和反应速率系数在时间离散点的数值,单体转化率X作为优化目标值;将各离散时间点的反应速率系数代入自由基聚合反应模型,变求解自由基聚合反应模型为大规模自由基聚合反应模型,通过序贯法求解大规模自由基聚合反应模型得到聚合物分子量分布曲线,特征点的分子量浓度值作为约束方程值;目标函数和约束方程导数值由差动方式提供,设置扰动量,以差分形式求出导数值,序列二次规划算法根据目标函数和约束方程函数值及导数值决定下一步寻优方向及步长,并判断下一步迭代点是否满足收敛条件,若收敛,返回最优解,以最优解作为初值计算自由基聚合反应模型聚合物分子量分布曲线,并计算与给定的聚合物分子量分布曲线误差分布曲线,若在-5%至5%之间,则认为最终的产品符合给定聚合物性能要求,否则产生新的一组初值重新开始优化,直至最后得到的误差曲线均在-5%至5%之间;若求解不收敛,返回改变初值重新优化,得到优化后的操作条件;方法具体步骤如下:1)读取聚合物分子量分布曲线数据,找出聚合物分子量浓度最大的链长点,得到链长点对应的聚合物分子量浓度;找出聚合物分子量浓度为链长对应的聚合物分子量浓度一半的链长点,得到链长点对应的聚合物分子量浓度,以这三个链长点的聚合物分子量浓度作为优化命题的不等式约束,令为反应温度,为反应时间,使满足:;其中为以为反应温度,为反应时间的操作条件下得到的链长点为的聚合物分子量浓度值;为需要优化的目标聚合物分子量浓度值;2)设置自由基聚合模型计算参数,计算最大链长点为约束中使用到的最大链长点,最大迭代次数设置为500,目标函数为最大化单体转化率;设置优化命题:3)随机产生初值:反应温度条件及反应持续时间,使用序列二次规划算法开始优化;4)给定操作条件,使用变步长变阶次后向差分方法求取小规模自由基聚合反应模型得到单体转化率,间隔为0.1分钟的从时间0到的离散时间点上的反应速率系数数值;将反应速率系数代入自由基聚合反应模型,得到大规模自由基聚合反应模型,使用序贯法求取大规模自由基聚合反应模型得到对应的聚合物分子量浓度;5)给定扰动大小,得到两组初值、,同理计算相应的单体转化率及聚合物分子量浓度值,采用有限差分方法,得到单体转化率和聚合物分子量浓度值在时间和温度的导数值:,,,,计算得到目标函数与约束方程的导数值;6)序列二次规划算法处理目标函数及约束方程在的函数值和导数值,寻优得到下一个迭代点,判断约束收敛条件与目标函数收敛条件,若不满足迭代收敛条件且迭代次数未超出则返回4),从继续寻优;(1)若迭代次数超出最大迭代次数且未收敛,则未找到最优解,求解失败;(2)若收敛条件满足,得到最优反应条件,以优化解为自由基聚合模型的反应条件,设置,利用序贯法模拟大规模自由基聚合模型,得到全链长范围的分子量分布数据,与给定的分子量分布曲线数据求相对误差,a)若,则所求的优化结果不满足聚合物分子量分布质量要求,求解失败;b)若,则所求的优化结果满足聚合物分子量分布质量要求,求解成功,返回最优解。...
【技术特征摘要】
【专利技术属性】
技术研发人员:范伟伟,陈曦,姚臻,邵之江,钱积新,
申请(专利权)人:浙江大学,
类型:发明
国别省市:86
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。