您好,欢迎来到华佗健康网。
搜索
您的当前位置:首页基于模腔压力的注塑过程注射速度设定曲线优化

基于模腔压力的注塑过程注射速度设定曲线优化

来源:华佗健康网
第61卷第5期 化工 学 报 Vo1.61 NO.5 May 2O10 2010年5月 CIESC Journa1 基于模腔压力的注塑过程注射速度 设定曲线优化 李 闯 ,王福利 ,刘 阳 ,常玉清 ( 东北大学信息科学与工程学院; 东北大学流程工业综合自动化教育部重点实验室,辽宁沈阳110004) 摘要:注塑成型中,把熔前速度恒定(熔体均匀填充模具)作为设定注射速度参数的重要依据。针对熔前速度 无法测量问题,提出对于一类注塑制品,用模腔压力随时问线性增加来近似熔前速度恒定。其次,利用实验数 据,设计并建立了模腔压力递归神经网络(RNN)模型。在此基础上,提出一种注射速度设定曲线的分级优化 策略,其中曲线分级是利用Douglas—Peucker方法自动确定速度分级数目和位置,优化过程则采取基于滑动窗口 的迭代优化方法来实现。实验结果表明,该优化方法可以根据模腔复杂程度自动选择合适的注射曲线分级数目, 并且其优化精度较常规等间隔分级优化方法更高。 关键词:模腔压力;注射速度曲线;Douglas—Peucker方法;迭代优化 中图分类号:TQ 320.662 文献标识码:A 文章编号:O438—1157(2010)05—1163—06 Optimization of inj ection velocity profile in inj ection molding based on cavity pressure LI Chuang ,WANG Fuli ,LIU Yang ,CHANG Yuqing ・ ( School of Information Science 8L Engineering,Northeastern University,Shenyang 1 1 0004,Liaoning,China Key Laboratory of Process Industry Automation of Ministry of Education,Northeastern University, Shenyang 110004,Liaoning,China) Abstract:Uniform melt front(MFV)velocity throughout the filling of the mold cavity is recommended as the design criterion for inj ection profile.Due to the difficulty of measuring MFV,a criterion that cavity pressure followed a ramp with time was proposed as an approximate substitute to constant MFV in this study.Then a recurrent neural network was developed for estimating the cavity pressure from online measurable variables,based on which the profile problem was transformed to an optimization problem for the predicted cavity pressure to track a given ramp.Subsequently the velocity profile was treated as a piecewise ramp one with respect to the increment of screw stroke,and its turning points and the slopes of each section were the optimization variables.The Douglas—Peucker method was introduced to determine the number and positions of separation points of the profile and an iterative optimization method with sliding window was presented to solve the optimization problem.Experimental results showed that this optimization method could automatically give the number of pieces of injection velocity profile and its precision of optimization was higher than the regular optimization method with equal fixed intervals. Key words:cavity pressure;injection velocity profile;Douglas—Peucker method;iterative optimization 2009—09—23收到初稿,2010 O1—25收到修改稿。 Received date:2009一O9—23. 联系人:王福利。第一作者: 李闯(1982),男,博士研 究生。 Corresponding author:Prof. WANG Full, wangfuli @ise.neu.edu.cn 基金项目:国家重点基础研究发展计划项(2009CB32O6O1);国家自然科学基金项目(60774068)。 目 Foundation itern: supported by the Nationa1 Basic Research Program of China( 2009CB320601)and the National Natural Science Foundation of China(60774068). ・j164・ 化工学报 第61卷 引 言 一个典型的注塑成型过程可分为注射、保压和 冷却三个阶段。在注射阶段,料筒内已塑化的高温 塑料熔体在螺杆的推动下,经流道、浇口注入封闭 的型腔,直至充满型腔约95 为止。整个注射过 程通常在几秒内完成,约占整个注塑生产循环周期 的5 ,甚至更少。尽管注射过程非常短暂,但其 对注塑产品的质量性能有着重要的影响。注射速度 设定的不合理将导致熔体前沿流动速度(以下简称 熔前速度)不稳,产生收缩不均、翘曲、甚至局部 烧焦等缺陷口 。 研究表明,理想的填充模式为,整个充填过程 中熔体以固定的熔前速度填充制品,并能同时到达 模具型腔的各个角落_2]。其中,后者主要与浇口的 位置有关,这在模具的设计阶段就已决定。而前者 则可通过注射速度设定进行在线调节,从而保证固 定的熔前速度,这便是本文研究的问题。 然而,由于熔体在封闭的模腔中行进,在模腔 形状未知的情况下,熔前速度通常是无法检测的。 Chen等口 研制出了熔料在模腔内填充位移的电容 传感器,它能有效地在线检测出熔体在模腔内的填 充位移和速度。但是,由于电容传感器安装较复 杂、增加模具生产成本,且检测精度易受干扰,文 献所提方法在生产实际中很难推广。郁滨等_4]通过 保持熔体沿充模方向的压力梯度保持恒定,来设定 最优注射速度曲线。然而,实际测量梯度时,螺杆 位置被用于替代不可测的熔体充模长度,这使得优 化结果在理论上偏离真实值。此外,还有通过数值 模拟技术来获取熔前速度的方法 。],但这类方法 都需要预先知道模腔的形状,而且优化结果准确与 否取决于所用机理模型的精确程度,然而建立一个 精确的机理模型是非常困难的。 对于速度曲线的优化问题,Chen等 伽先建立 熔前速度软测量模型,之后将注射曲线划分为等间 隔的固定级数,最后利用序贯二次规划方法进行优 化。但是对于形状多样的模具,等间隔划分方法在 一定程度上了速度曲线的变化,当然通过增加 分级段数可以更精确地描述速度曲线,但这样会增 加优化的计算负担,并且实际注塑机对分级数都有 一定的,不可能无增加级数。Zhang 等|1叩通过一次恒定注射速度实验得到的填充位移 曲线,自动确定与实验模具相对应的注射速度曲线 分级位置。该方法属于不等问隔分级法,能根据模 腔形状更合理地确定分级位置。但最优速度曲线的 分级特征点往往较填充位移曲线的分级特征点更 多,利用后者来简单地代替前者,往往不能充分描 述最优注射速度曲线,从而降低优化结果的精度。 本文针对一类注塑制品,首先根据机理分析, 提出以模腔压力是否随时间线性上升,作为模腔是 否均匀填充判定依据。其次,通过实验数据,建立 模腔压力的RNN模型来近似注射填充过程。之 后,基于此模型,提出一种基于Douglas—Peucker 曲线特征提取方法和滑动窗口的分级迭代优化策 略,该方法先优化出近似最优速度曲线,再利用 Douglas—Peucker方法对该曲线进行特征提取,确 定曲线分级数和分级点近似位置,随后通过迭代优 化方法获得精确最优注射速度曲线。最后,通过实 验说明了该方法较常规等间隔分级优化方法的 优势。 1 熔体的均匀填充 注射阶段熔体在模腔内的运动属于压力流动。 为简化问题,对熔体填充如图1所示的矩形模腔, 作如下假设: (1)模具是薄壁的,即模腔厚度2H远小于宽 度W; (2)熔体不可压缩,且其流动为等温流动; (3)质量力和惯性力忽略不计。 图1矩形模腔 Fig.1 Rectangular mold cavity 取模腔长度L方向上dz长的微元,其两侧压 力梯度为 等一吉『dZ H  I 2 WH ’lQT  ㈩ … 这里 代表非牛顿流变指数;k代表熔体黏度 系数,它是熔体温度和剪切速率的函数;Q代表熔 体的体积流速 川。 令 代表熔体前沿的平均流速,即熔前速度, 那么显然 Q一2WHuf (2) 第5期 dl一 李闯等:基于模腔压力的注塑过程注射速度设定曲线优化 dt (3) 1165。 必须对其降维才能实施优化。由于速度曲线的差 异,每次注射时间会不同,但此时对应的熔体填充 量是相同的,因此可按螺杆行程(填充量)把速度 将式(2)和式(3)代入式(1)中可得 一 一 (\ )  + (4) 由式(4)可知,如果认为 与足近似不变, 对于截面积为矩形的薄壁制品,在厚度不变的条件 下,只要保证 恒定,就可以保证 恒定。换言 之,对于薄壁制品,只要其薄厚不变,通过调整注 射速度设定,使模腔压力随时问线性增加,就能保 证薄壁制品的均匀填充。 2 模腔压力RNN模型 2.1建模变量的选择 在模具和注射材料固定的条件下,模腔压力P 是熔体在模具内熔体流速的反映。那么与熔体流速 最相关的输入变量是螺杆位置SD和注射速度n , 前者可反映模具形状变化的位置,后者决定熔体流 速的快慢。熔体温度和模具温度对熔前速度也有影 响,但通过实验数据观察,其值在注射过程中变化 不大,所以未被选为输入。此外,由于模腔压力的 累积作用,建模时不仅要选择其当前时刻 的采样 值CP ,还要选择前是时刻的采样CP 作为输入, 才能得到合理的预测值,这里尼称为CP的阶次。 因此本文要建立的是一个ARMA模型,它是一个 非线性动态模型。 2.2 RNN模型 递归神经网络(RNN)是拟合非线性动态函 数关系的有力工具,设计模腔压力RNN的关键是 确定输入阶次、网络结构和网络权值_1 。在训练 和测试数据已知的情况下,这里采用凑试法,寻找 使测试数据泛化误差达到最小的各输入变量阶次和 网络结构,而网络权值的训练则采用Levenberg— Marquart优化方法。 离线优化注射速度时,还需要给出与注射速度 n 对应的螺杆位置SD,来预测模腔压力CP。由 于n,与SD之间的关系相对简单,这里使用线性 模型对其预测,用凑试法确定模型阶次,用最小二 乘方法对模型参数进行拟合 。 3 注射速度的优化方法 3.1注射速度曲线优化问题 由于优化对象是一条无穷维的速度曲线,所以 曲线分为m级。令初始注射速度为 ,自模腔压 力上升时刻 。开始进入第1级,速度在第1级随 时问线性变化,变化斜率为R ,到达螺杆行程s 后,转为2级,以此类推,直至填充完成。令该时 刻为£ ,那么连续的注射速度曲线C可简化为 级的分段折线,且该折线可表示为2m维向量 C一[ R R …R S S2…S ]or C一 s] 其中,R一[ Rl R2…R ],Js—Es1 S2…s 1]。 这样速度曲线优化问题就可以描述为 minr ̄RNN(R,s, )一slope(f )]zdf (5) R,S Jtn 3.2迭代优化 仿真经验表明,序贯二次规划(SQP)方法更 适合于变化斜率R的优化,而单纯形法则适于s 的优化。因此这里把问题(5)拆分为两阶段进行, 先固定S对R优化,再固定R对S优化,如此循 环,直至满足收敛准则。 此外,在进行固定s对R优化时,由于模腔 压力的累积作用,那么第i级模腔压力曲线P ,只 跟第i级速度曲线c 及其之前级的速度曲线有关。 并且一旦前i级速度曲线优化好,即前i级模腔压 力曲线已达到设定,那么在优化P汁 时,就不需 要再改变之前的速度优化结果,只要优化C汁 就可 以了。但实际上由于动态模型为是阶,还需要改变 甚至更前面的速度曲线来调整P汁 。 这样就可以把原优化问题拆分成若干子优化问 题求解。具体来说,就是从前向后,每次只优化连 续的M级速度曲线,完成后再一级级向后滑动, 直到所有级都优化完。这里将这种优化方法称为滑 动窗口法,采用该方法可以降低优化问题的维数, 大大提高优化计算速度和优化结果精度。 3.3速度曲线的自动分级 速度曲线分级多少应与被充填模腔形状复杂程 度有关,模腔形状越复杂分级点越多。恒定注射速 度注射时,模腔形状的变化可以反映在熔体充填长 度位置或喷嘴(模腔)压力的拐点上,这些拐点信 息可用于指导分级位置点选取l4 。但通过这种方 法选取特征点常出现遗漏,从而降低优化精度,因 为真实注射速度曲线的特征点往往多于上述的充填 长度曲线拐点。 ・ 】】66 ・ 化工学报 第6l卷 为了更全面地定位速度最优曲线特征点,这里 的数据采集输出板卡,采样周期为10 ms。控制系 先利用RNN模型初步优化出速度曲线,然后用这 个解来选取特征点,显然该结果更能代表最优曲线 的特征。这里选取Douglas—Peucker矢量曲线特征 统的软件部分是在嵌入式实时操作系统QNX平台 下开发,采用C语言编程。为便于比较研究,这 里采用如图2所示的工形和梯形两种模具,它们具 有相同尺寸的主流道和扇形浇口,工件厚度均为2 提取方法I】引,用于对已知矢量曲线分级点的提取 定位。通过调整该方法的分级精度参数 ,得到特 征点由少到多的一系列原始曲线的插值曲线,再从 mm,图中P点为压力传感器位置。实验用塑料材 料为HDPE。 中挑选出能足够表达最优速度曲线的一组特征点, 作为速度曲线分级点。 3.4算法流程 注射速度曲线离线智能分级优化算法可分为 4步。 (1)给出设定模腔压力曲线P =slope( 一 )。 (2)获得初步优化睦线C。。先把速度曲线等 间隔分为m级得到s。, 应超过注塑机注射速度 最大分级数。在s。下使用滑动窗口SOP方法优化 得 ,窗口包含优化级数M应与模型输入阶次有 关。至此得到曲线C。一[R。,s。]。 (3)提取特征插值曲线c 。使用Douglas— Peucker方法,通过调整分级精度参数 ,获得c。 的一组插值曲线 ,i表示C 里含有分级点的数 目。设定阈值 ,从c 中挑选出集合 ( … —C。IJ< ) (6) 其中I1.1l代表两条曲线上所有对应采样点速度值 的均方根偏差(RMSE)。再从集合式(6)中挑选 出最小i对应的一组插值曲线C 一[R ,Js ],s 即为最优速度曲线的特征分级点。 (4)迭代优化获得最优速度设定曲线c 。利 用迭代优化方法,先固定s ,以R 为初始值,使 用滑动窗口sQP优化得R。;再固定R ,以s 为 初始值,使用单纯形法优化得s ,这样就得到新 曲线C。一[R。,s。]。如此迭代,直至满足如下任 一条件迭代终止:①lI P 一P。JI≤e;②迭代次数达 到最大次数,z,从而获得Cr。这里P。为第q一 1次迭代的速度 线对应的模腔压力输出,P。为 设定模腔压力曲线。 4 优化结果 实验所用注塑机为宁波海天公司生产的往复螺 杆式注塑机,其型号为HTF58X2一B。模腔压力是 由安装在模腔内的压力传感器(Green Sensor TPT一250)测量的,螺杆位置是由电子尺(NOVO Technik type—LMN一130)测量的。系统采用研华 L }_t——————————————————————————————————————————,I  (b)B 图2实验用的两种模具 Fig.2 Two molds used in experiment(mm) 4.1模型建立 首先通过实验数据获取训练和测试数据。针对 模具大小,将注射速度的范围限定在10~25 mm・s 之间,选择3类19条速度曲线进行实验, 并记录其对应模腔压力输出,其中训练数据14组, 测试数据5组,具体情况如表1所示。所有数据截 取从模腔压力开始上升(熔料进入模腔),到注射 保压切换时刻(基本注满模腔),作为有效填充 量。两个时刻对应的螺杆位置差即为该模具的充 填量。注射保压切换方法采取预设的螺杆位置转 换法 。 这里仅采用模具A的实验数据来确定神经网 络结构,因为模具A较模具B形状复杂。通过凑 试法,确定单隐层网络,隐层节点数为7,各输入 阶次如图3所示。采用该结构,经网络权值训练 后,模具A的训练RMSE=0.1038,测试RMSE= 第5期 李闯等:基于模腔压力的注塑过程注射速度设定曲线优化 表1输入数据的实验设计 Table 1 Design of injection velocity profile 2 景 至 △ ①If step profile jumps from high speed to low speed,turning point is situated on 2/3 of filling amount of mold,otherwise it is on 1/3 of that amount. 星 墨 图3模腔压力递归神经网络模型 0 图 Fig.3 RNN of cavity pressure 0.5008;模具B的训练RMSE一0.0861,测试 ∞ ㈣融 形 制 工 RMSE一0.3987。 此外,再通过表1的训练和测试数据,确定螺 杆位置线性预测模型的阶次为2,表示如下 比 O 结 3。 SD 一a1SD—l+12zⅣ +a3Ⅳ +a4Ⅳ,卜2 果 f ~ J 并用最小二乘法拟合出模型的参数,模具A的训 m 练RMSE=0.2304,测试RMSE=0.4564;模具B  :一一一一一一vd e、l,o。 ci 。ty 。p oif。1。 的训练RMSE一0.1032,测试RMSE=0.1571。 喜 II.I  . .4.2优化结果 。 50 100 step 根据3.4节所列的优化步骤,首先确定目标模 (b) 腔压力曲线斜率,对于模具A和B,设定slope分 图5梯形制品优化结果 别为0.12和0.08。随后,再等间隔分级,分级数 Fig.5 Optimization results of mold B m一15,滑动窗口包含2级速度曲线的条件下,优 模具A,其分级数为8,而对于模具B,其分级数 化出初步最优速度曲线C。。再通过C。找出一系列 为4。最后利用迭代优化方法优化C ,得到最终的 特征点插值曲线,如表2所示。取阈值 一2,根 速度优化曲线C 。 据式(6)可得到简化特征曲线C ,对于较复杂的 工形制品和梯形制品的Cr分别如图4(b)和 表2速度曲线与其各个插值曲线的RMSE Table 2 RMSE between velocity profile and its each interpolation curve ・j J68・ 化工学报 第61卷 图5(b)所示,其对应的模腔压力曲线如图4(a) 和图5(a)中虚线所示。可见对于工形制品的C 在4O~100采样之间先降后升,而梯形制品c 在 50个采样后逐渐增加,这完全符合模具截面积变 化的定性分析。C 前段的上升是由于初始速度 。 较小造成的。 表3给出C。、C 与C 对应的模腔压力输出, 与设定模腔压力所有对应采样点的均方根偏差。由 比较可知,C 较C。虽然在分级数方面少了,但其 对应输出却较C。更加靠近设定模腔压力了。这说 明利用智能分级方法得到的分级位置更合理,达到 相同输出精度所需的分级数更少。此外,C 较C 的RMSE显著减小,说明了迭代优化方法的有 效性。 表3优化结果比较 Table 3 Comparison of optimization results 5 结 论 本文针对一类厚度不变的注塑制品,提出将模 腔压力是否随时间线性上升,作为制品是否均匀填 充判定依据。而后,在利用实验数据建立的模腔压 力RNN模型的基础上,提出了一种注射速度的自 动分级和迭代优化方法,该方法可自动按照模腔复 杂度,给出适合该模型的速度分级方式,并能准确 地获得与模具相对应的最优速度曲线。实验结果表 明,使用智能分级优化方法得到的注射速度曲线分 级数更少,优化精度更高。此外,采用何种手段 (比如选取某种质量指标)来证明文中优化依据的 有效性,以及如何在RNN模型失配时,自适应地 调整速度设定优化曲线,将是未来进一步的研究 工作。 References [1] Li Dequn(李德群).Principle,Method and Application of Modern Plastic Injection Molding(现代塑料注射成型的原 理、方法与应用).Shanghai:Shanghai Jiao Tong University Press,2005 [2] Liu Laiying(刘来英).Injection Molding Technics(注塑 成型工艺).Beijing:China Machine Press,2004 Chen Xi, Chen Guohua, Gao Furong. Capacitive transducer for in—mold monitoring of injection molding. Polymer Engineering and Science,2004,44(8): 1571 1578 [4] Yu Bin(郁滨),Yuen M M F.Research for automatic setting velocity profiles based on nozzle pressure curves. China Plastics(中国塑料),2002,16(2):63—71 E5] Zhai Ming(翟明),Gu Yuanxian(顾元宪),Shen Changyu(申长雨).Optimization on injection rate of filling phase in injection molding.Acta Polymerica Sinica(高分 子学报),2003(1):35 38 [6] Cui Kun(崔醌),Ouyang Jie(欧阳洁),Zheng Supei(郑 素佩),Zhao Zhifeng(赵智峰),Ruan Chunlei(阮春蕾). Level Set method for 3D meh front surface.Journal of Chemical Industry and Engineering(China)(化工学报), 2008,59(12):3020—3026 [7] Liu Yongzhi(刘永志),Zhao Zhenfeng(赵振峰),Ma Lan(马兰),Shen Changyu(申长雨).Numerical simulation for velocity and pressure in filling stage of injection molded parts with uneven thickness.Journal of Chemical Industry and Engineering(China)(化工学报), 2009,60(7):1818 1822 [8] Chen X,Gao F R,Chen G H.A soft sensor development for melt—-flow--length measurement during injection mold filling.Materials Science and Engineering A,2004,384: 245—254 [9] Chen X,Gao F R.Profiling of injection velocity for uniform mold filling.Advances in Polymer Technology, 2006,25(1):13 21 [10] Zhang L,Chen X,Gao F R,Lu J n An intelligent profiling strategy of injection velocity for uniform cavity filling based on recurrent neural network model//18th International Congress of Chemical and Process Engineering.Praha:Czech Republic,2008:24—28 [11] Zhou Chixing(周持兴),Yu wei(俞炜).Polymer Processing Theory(聚合物加工理论).Beijing:Science Press,2004 [12] Faur C,Cougnard A,Drefus G,Le Cloirec P.Modeling the breakthrough of activated carbon filters by pesticides in surface waters with static and recurrent neural networks. Chemical Engineering Journal,2008,145:7-15 [13] Wang Jing(王净),Jiang Gangwu(江刚武).Researching and realization of the quick compression method aimed at the non—topology vector data.Acta Geodaetica et Cartographica Sinica(测绘学报),2003,32(2):173—177 [14] Li Chuang(李闯),Wang Fuli(王福利),Chen Juan(陈 娟),Dong Zhaoran(董召然).Wavelet detection of filling— to—packing switchover point in injection molding.Journal of Chemical Industry and Engineering(Chian)(化工学报), 2009,60(5):1213—1218 

因篇幅问题不能全部显示,请点此查看更多更全内容

Copyright © 2019- huatuo0.com 版权所有 湘ICP备2023021991号-1

违法及侵权请联系:TEL:199 1889 7713 E-MAIL:2724546146@qq.com

本站由北京市万商天勤律师事务所王兴未律师提供法律服务