三维结构复合材料预制体轮廓约束机理_应哲赟.pdf

第 51 卷，第 6 期2023 年 6 月工程塑料应用Vol.51，No.6Jun.2023ENGINEERING PLASTICS APPLICATION三维结构复合材料预制体轮廓约束机理应哲赟1，单忠德1，2，刘丰1，吴晓川1，杜悟迪1，陈哲1(1.中国机械科学研究总院集团有限公司先进成形技术与装备国家重点实验室，北京 100083；2.南京航空航天大学，南京 210016)摘要：为实现三维结构复合材料预制体的高性能、快速成形织造，研究了基于柔性导向三维多针成形原理的预制体轮廓约束方法，利用图像处理技术获得了预制体纤维环形状特征，建立了工艺参数对纤维环形状特征影响的理论模型，搭建了纤维环成环试验装置，探究了织造针伸出长度、织造层数和织造高度对纤维环形状特征的影响规律。并通过该理论模型对纤维环进行形状特征调控，完成典型件的试织。结果表明，纤维环侧边高度以及最大内接圆直径随伸出长度和织造高度增大而增大，随织造层数增大而减小，该模型预测与实际测得的纤维环侧边高度及最大内接圆直径的最大相对误差分别为6.29%和4.63%，进一步通过调节织造高度来改变纤维环的形状特征，完成层数为90层的预制体织造，说明该理论模型可以准确描述纤维环形状特征，并可通过该模型对纤维环形状特征进行调控，扩大了预制体轮廓约束工艺窗口，为实现预制体轮廓自动、高质量约束提供了理论指导和实践基础。关键词：预制体；轮廓约束；纤维环；形状调控中图分类号：TH16 文献标识码：A 文章编号：1001-3539（2023）06-0054-08Contour Constraint Mechanism of Three-Dimensional Structure Composite PreformYing Zheyun1，Shan Zhongde1，2，Liu Feng1，Wu Xiaochuan1，Du Wudi1，Chen Zhe1(1.State Key Laboratory of Advanced Forming Technology and Equipment，China Academy of Machinery Science and Technology Group Co.，Ltd.，Beijing 100083，China；2.Nanjing University of Aeronautics and Astronautics，Nanjing 210016，China)Abstract：In order to realize the high-performance and rapid forming weaving of three-dimensional structural composite preforms，the contour constraint method of preforms based on the principle of flexible-oriented three-dimensional multi-needle forming was studied.The shape characteristics of preform annulus were obtained by image processing technology.The theoretical model of the influence of process parameters on the shape characteristics of annulus was established，and the annulus forming test device was built.The influence of weaving needle extension length，weaving layers and weaving height on the shape characteristics of annulus was explored.Through the theoretical model，the shape characteristics of the fiber annuluses were regulated，and the weaving of typical parts was completed.The results show that the side height and the maximum inscribed circle diameter of the fiber annulus increase with the increase of the extension length and the weaving height，and decrease with the increase of the weaving layers.The maximum relative error between the model prediction and the actual measurement of the side height and the maximum inscribed circle diameter is 6.29%and 4.63%respectively.The shape characteristics of the fiber annulus are further changed by adjusting the weaving height，and the preform weaving with 90 layers is completed，which shows that the theoretical model can accurately describe the shape characteristics of the fiber annulus，and the shape characteristics of the fiber annulus can be regulated by the model.It expands the process window of the preform contour constraint，and provides theoretical guidance and practical basis for realizing the automatic and high-quality constraint of the preform contour.Keywords：preforms；contour constraint；fiber annulus；regulation of shape三维增强复合材料由于其厚度方向引入了Z向纤维束，克服了二维复合材料抗分层能力和抗冲击性差的缺陷，受到了越来越多的国内外学者的关注1。由于三维增强复合材料具有轻质、高比强度、doi:10.3969/j.issn.1001-3539.2023.06.009基金项目：国防科技基础加强计划项目通信作者：单忠德，研究员，研究方向为数字化机械装备与先进成形制造收稿日期：2023-03-17引用格式：应哲赟，单忠德，刘丰，等.三维结构复合材料预制体轮廓约束机理J.工程塑料应用，2023，51(6)：5461.Ying Zheyun，Shan Zhongde，Liu Feng，et al.Contour constraint mechanism of three-dimensional structure composite preformJ.Engineering Plastics Application，2023，51(6)：5461.54网络首发时间：2023-06-21 14:28:37网络首发地址：https:/ 纤维成环工艺分析与建模在柔性导向三维多针织造过程中，织造针携带纤维束伸出导向阵列一定长度后，织造针、导向阵列和纤维束形成的封闭图形称之为纤维环，如图1所示。纤维环形状特征决定了置入辅助边棒的难易程度。当纤维束在一定张力条件下，纤维环近似于直角三角形。基于纤维环直角三角形假设和三维织造工艺特点，建立如图2所示的纤维环成形理论分析模型。在图2中，建立以O1点为原点、垂直向下为x轴正方向、水平向右为y轴正方向的坐标系。图2中模型各参数物理意义见表1。织造针伸出导向阵列后，为保证纤维束处于保3EIDDUr3434(n图1纤维环成环过程示意图l1l2l3HSlhDDb/DO1xy3EI34K图2纤维环成环结构示意图表1纤维成环模型参数物理意义表参数名称HLSlhbl1l2l3物理意义织造高度纤维预留长度织造针尖至导向阵列距离织造针伸出长度纤维环侧边高度进针方向预制体幅宽织造针与纤维束夹角辅助边棒切点至织造针尖距离辅助边棒纤维束包覆段辅助边棒切点至导向阵列距离55工程塑料应用2023 年，第 51 卷，第 6 期持绷直状态，S应满足以下关系：S=b+l+4r(1)式中：r辅助边棒半径。在纤维束送纱过程中，纤维预留长度保持不变，送纱完成后纤维预留长度由三部分构成，分别为l1，l2，l3。纤维预留长度及三部分长度分别为：L=H-2(n-1)r 2+S2(2)l1=rsin+r+b+ll l2+h2(3)l2=(-)r(4)l3=r(5)式中：L纤维预留长度；n织造层数。辅助边棒沿导向阵列边缘置入，纤维环侧边高度可以决定辅助边棒沿x正方向的运动空间，联立式(2)、式(3)、式(4)、式(5)可得侧边高度与其它各参数的关系：rsin()arctanhl+r+b+ll l2+h2+(-arctanhl)r+r=H-2(n-1)r 2+S2(6)为进一步描述纤维环形状特征，提出其最大内接圆直径和圆心位置参数，最大内接圆直径可以决定可置入辅助边棒的直径大小，圆心位置可以为轮廓约束装置置入位置提供参考，最大内接圆直径确定如式(7)。d=2hlh+l+l2+h2(7)式中：d最大内接圆直径。圆心坐标值为：(d/2，d/2)。根据式(6)和式(7)可知，当参数一定时，可通过纤维环侧边高度以及最大内接圆直径来表征纤维环形状特征。2 试验材料及方法选取直径为1 mm的304不锈钢棒，布置276的矩形导向阵列，导向棒中心距为2.4 mm，即进针方向导向阵列等效幅宽为180 mm。采用T300-3K碳纤维进行纤维束铺放，使用直径为3 mm的辅助边棒进行轮廓约束，通过MV-CA060-10GC工业相机进行纤维成环图像采集，搭建纤维环数据采集试验台，如图3所示。试验环境保持相对湿度为30%50%，温度为2025，设定织造排针运动速度为20 mm/s，开展纤维成环理论模型验证试验以及探究织造针伸出长度、预制体织造层数以及织造排针初始高度对纤维环形状特征的影响规律试验。对纤维环图像依次进行灰度化、Otsu二值化、ROI区域提取、最大连通域提取、8邻域边缘提取处理，可得纤维环侧边高度、最大内接圆直径以及其圆心坐标值，其图像处理流程及分析过程如图4和图5所示。3 结果与讨论在预制体轮廓约束过程中，理想情况下辅助边棒在圆心位置以一定初始速度进入多束纤维环中最大内接圆形成的圆柱形通道，对纤维束进行固定，由于织造针刚性较差，在纵向存在一定的挠度，经过多次试验表明，其轮廓可约束的纤维环最大内接圆直径需大于10 mm。在三维多针织造成形过程中，对纤维环形状特征规律变化的把握至关重要，因此针对各工艺参数对纤维环形状特征的影响规律展开探究。图3纤维环数据采集试验台%Otsu08EEED*0h0dEE)8EEED*D*2%K/0)%KAAK%!d,KBKE3%4(a)(b)a图像分析流程；b最大内接圆直径(d)求解流程图4图像分析流程图56应哲赟，等：三维结构复合材料预制体轮廓约束机理3.1工艺参数对纤维环形状特征的影响规律(1)伸出长度对纤维环形状特征的影响规律。选取织造高度为100 mm，织造层数为1层。根据织造针的有效长度和布置导向阵列的幅宽，可以知道织造针伸出长度的最大值为95 mm；由织造经验可知，织造针伸出长度最小值为30 mm。获得织造针伸出长度对纤维环侧边高度、最大内接圆直径及圆心坐标的影响规律，如图6所示，经20余次试验得到的最大内接圆圆心坐标斜率统计如图7所示。由图6可以看出，在织造高度、织造层数一定时，织造针伸出长度越大，纤维环侧边高度越大，纤维最大内接圆直径也越大；纤维环侧边高度试验值和理论值最大相对误差为6.29%，纤维环最大内接圆直径试验值和理论值最大相对误差为4.63%。圆心x，y坐标值相对于理论x，y坐标值最大相对误差分别为7.25%，7.63%，圆心x，y坐标值随着织造针伸出长度的增大，均逐渐增大；实际圆心坐标斜率均匀分布于理论圆心坐标斜率两侧(图7)，最大斜率为051015200.000.250.500.751.001.250.9430405060708090(AOl/mm11.18图7不同伸出长度(l)下最大内接圆圆心坐标斜率统计分布图34(DP/mmJ/mm3040506070809051015202530/mmJ/mmAO)304050607080905101520253035AO)30405060708090468101214x/mmJ/mmAO)30405060708090468101214y/mmJ/mmAO)(a)(b)(c)(d)a纤维环侧边高度；b最大内接圆直径；c圆心x坐标；d圆心y坐标图中坐标值为圆心位置到设定原点的距离，下同图6伸出长度对纤维环形状特征的影响规律(a)(b)(c)(d)10020030040001002003004005006007008002%y2%x10020030040001002003004005006007008002%y2%x(e)(f)a原图；b灰度图；c最大连通域图；d边缘轮廓图；e抽稀处理图；f最大内接圆图图5图像分析过程图57工程塑料应用2023 年，第 51 卷，第 6 期1.18，最小的斜率为0.94。在该工艺参数下，考虑到进针侧预制体幅宽和织造针有效长度，为了保证织造效率(较长的织造针伸出长度需要较长的运动时间)和轮廓顺利约束，织造针伸出长度以4060 mm为佳。(2)织造层数对纤维环形状特征的影响规律。选取织造针伸出长度为 40 mm，织造高度为100 mm，探究织造层数对纤维环形状特征的影响规律。通过对纤维环的图像分析，获得织造层数对纤维环侧边高度、最大内接圆直径及圆心坐标的影响规律和圆心坐标斜率分布(共试验50次)，如图8和图9所示。由图8和图9可知，纤维环侧边高度和最大内接圆直径与织造层数均呈负相关，纤维环侧边高度试验值和理论值最大相对误差为5.59%，最大内接圆直径试验值和理论值之间的最大相对误差为3.31%，最大相对误差均小于10%。圆心x坐标值的最大相对误差为5.08%，圆心y坐标值的最大相对误差为3.83%。随着织造层数的增加，圆心坐标靠近设定原点，当织造层数大于8层时，圆心坐标斜率的试验值基本分布于理论值上方(图9圆圈之处)。这是由于随着织造层数的增加，辅助边棒增多，辅助边棒与纤维的摩擦力增大，从而导致纤维张紧力增大，织造针在纤维张紧力作用下产生沿x方向的挠度变形，使得实际圆心坐标沿x轴负方向移动，导致试验值大于理论值。为便于将辅助边棒置入纤维环中，所形成的纤维环最大内接圆直径应不小于10 mm，则在该选定织造参数下，理论最大可织造层数为10层。当确定织造工艺参数时，可根据纤维成环理论模型进行最大可织造层数的预测。(3)织造高度对纤维环形状特征的影响规律。进一步地，探究织造高度对纤维环形状特征的影响规律，选取织造针伸出长度为40 mm，织造层数为1层，获得织造高度对纤维环侧边高度、最大内接圆直径及圆心坐标的影响规律和圆心坐标斜率分布，如图10和图11所示。由图10和图11可知，随着织造高度的增加，纤维环侧边高度和最大内接圆直径增加，最大内接圆圆心坐标x，y值也随之变大。纤维环侧边高度试验值和理论值最大相对误差为4.30%，最大内接圆直径试验值和理论值之间的最大相对误差为2.78%；010203040500.000.250.500.751.001.251.500.781.311 1 2 3 4 5 6 7 8 9 10(AOn/图9不同织造层数(n)下圆心坐标斜率统计分布图024681068101214161834(DP/mm3E/AO)02468106810121416/mm3E/AO)02468104.55.05.56.06.57.07.5x/mm3E/AO)02468104.55.05.56.06.57.07.5y/mm3E/AO)(a)(b)(c)(d)a纤维环侧边高度；b最大内接圆直径；c圆心x坐标；d圆心y坐标图8织造层数对纤维环形状特征的影响规律58应哲赟，等：三维结构复合材料预制体轮廓约束机理圆心斜率统计分布图中，实际最大斜率为1.19，最小斜率为0.83，均匀分布于理论值两侧。在实际织造过程中，选择较大的织造高度，使得纤维环形状特征便于轮廓约束。3.2典型预制体试织基于工艺参数对纤维环形状特征影响的理论计算模型，选取典型复合材料预制体进行织造验证，其织造技术规格见表2。导向阵列布置为4343，如图12所示，根据纤维密排堆积可计算所需层数为90层。在织造过程中选定织造针伸出长度为50 mm，初始织造高度为100 mm，通过式(6)、式(7)计算该织造参数下的纤维环形状特征，如图13所示。由图13可知，在该织造参数下，当织造层数为21层时，最大内接圆直径为9.44 mm，不满足边棒置入要求，即其最大可织造层数为20层。根据纤维环形状特征理论计算模型可知，可通过改变织造针伸出长度或织造高度来调节纤维环的最大内接圆直径，从而获得更多的织造层数以保证预制体的可织造性。以最大内接圆直径10 mm为目标，通过改变织造针伸出长度来调节纤维环的形状特征，其计算结图12导向阵列布置图1001201401601802001216202428323634(DP/mm3EP/mm100120140160180200101214161820222426/mm3EP/mmAO)1001201401601802006789101112x/mm3EP/mmAO)1001201401601802006789101112y/mm3EP/mmAO)(a)(b)(c)(d)a纤维环侧边高度；b最大内接圆直径；c圆心x坐标；d圆心y坐标图10织造高度对纤维环形状特征的影响规律0510152010152025303534(DP,3E/34(DP/mm81012141618202224,/mm9.44图13各层纤维环理论形状特征表2预制体织造技术规格纤维规格x，y向6K1z向6K4外形尺寸长/mm100宽/mm100层厚/mm20010203040500.000.250.500.751.001.250.831.191 100 110 120 130 140 150 160 170 180 190(AOH/mm图11不同织造高度(H)下圆心坐标斜率统计分布图59工程塑料应用2023 年，第 51 卷，第 6 期果见表3。由表3可知，当织造层数为24层时，织造针伸出导向阵列的距离为257.31 mm，受织造装备限制，织造排针末端将与导向阵列发生碰撞。所以通过调节织造针伸出导向阵列的长度来调整纤维环形状特征时，其最大织造层数为23层。同样的，以最大内接圆直径10 mm为目标，通过改变织造高度来调节纤维环的形状特征，其计算结果如图14所示。由图14可知，为保证每一层纤维环最大内接圆直径为10 mm，织造层数与织造高度呈线性相关，即每一层调整高度为3 mm，由于受设备允许最大织造高度及导向棒有效长度限制，在103层以下可以任意织造。在典型预制体样件的试织过程中采用通过织造高度来调节纤维环形状特征，完成织造。预制体样件试织如图15所示。4 结论根据复合材料柔性导向三维多针织造成形工艺特点，基于预制体轮廓约束方法，建立了织造工艺参数对纤维环形状特征影响的理论计算模型，搭建了纤维环测量试验装置，并进行了模型验证和规律探究，通过实验结果得到以下结论。(1)不同工艺参数下，纤维环侧边高度试验值与理论值的最大相对误差为6.29%，最大内接圆直径试验值与理论值的最大相对误差为4.63%，证明该理论模型可以准确描述纤维环形状特征。(2)通过探究纤维环成形规律，发现纤维环形状随着织造针伸出长度和织造高度增大而增大，随织造层数的增大而减小，为保证所形成的纤维环便于辅助边棒的置入，根据织造预制体幅宽和织造针有效长度，在合理范围内，可选取较大的织造针伸出长度和织造高度。(3)进一步依据提出的理论计算模型，通过改变织造高度调整纤维环形状特征，完成典型预制体样件的试织。(4)笔者所建立的纤维环理论计算模型，为后续实现预制体轮廓自动约束技术提供理论指导和技术支撑，进一步地提高三维多针织造成形设备的自动化程度，缩短预制体的成形周期。参 考 文 献1 荆云娟，韦鑫，张元，等.三维编织复合材料的发展及应用现状J.棉纺织技术，2019，47(11):7984.Jing Yunjuan，Wei Xin，Zhang Yuan，et al.Development and application status of three-dimensional knitting composite materialJ.Cotton Textile Technology，2019，47(11):7984.2 刘澜涛，陈冠宇，易黎明，等.固体火箭发动机炭/炭复合材料喉衬烧蚀行为研究J.炭素技术，2021，40(3):3135.Liu Lantao，Chen Guanyu，Yi Liming，et al.Study on the ablation mechanism of C/C composite throat for solid rocket motorJ.Carbon Techniques，2021，40(3):3135.3 张天昊，查柏林，李志农，等.轴棒法编织C/C复合材料等离子烧蚀性能分析J.兵器装备工程学报，2021，42(9):237243.Zhang Tianhao，Zha Bailin，Li Zhinong，et al.Analysis on plasma ablation properties of axial carbon rod woven 4D C/C compositesJ.Journal of Ordnance Equipment Engineering，2021，42(9):237243.4 许建伟，李春林，高凌，等.三维织物复合材料天线罩的应用研究J.电子机械工程，2018，34(2):3639.Xu Jianwei，Li Chunlin，Gao Ling，et al.Study on application of 3D fabric composite radomeJ.Electro-Mechanical Engineering，2018，34(2):3639.5 孟祥福.三维编织复合材料叶轮制件的RTM工艺设计D.西安：西安工程大学，2019.Meng Xiangfu.RTM process design of 3D braided composite impeller partsD.Xi an:Xi an Polytechnic University，2019.6 Emonts C，Grigat N，Merkord F，et al.Innovation in 3D braiding 20 30 40 50 60 70 80 901001101201001502002503003504003EP,3E/3EP/mm103.739.09.510.010.511.011.512.0,/mm图14织造高度调节纤维环形状特征图15典型预制体件试织表3伸出长度调节纤维环形状特征织造层数/层21222324伸出长度/mm54.9467.53119.64257.31最大内接圆直径/mm1010101060应哲赟，等：三维结构复合材料预制体轮廓约束机理technology and its applicationsJ.Textiles，2021，1(2):185205.7 Wei Qingsong，Yang Dan，Gao Bo.Ballistic penetration simulation of a 3D woven fabric using HigH strain-rate dependent yarn modelJ.The Journal of The Textile Institute，2022，113(7):1 4011 410.8 翟军军，王露晨，孔祥霞.三维多向编织复合材料温度效应综述：热传导、热膨胀性质和力学响应J.复合材料学报，2021，38(8):2 4592 478.Zhai Junjun，Wang Luchen，Kong Xiangxia.A review on temperature effect of 3D multi-directional braided composites：Heat conduction，thermal expansion properties and mechanical responseJ.Acta Materiae Compositae Sinica，2021，38(8):2 4592 478.9 单忠德，刘丰.复合材料预制体数字化三维织造成形M.北京：机械工艺出版社，2019.Shan Zhongde，Liu Feng.Digital precision forming technology and equipment for composite materialsM.Beijing：China Machine Press，2019.10 李文启.剑杆织机空间四连杆引纬机构动力学分析与仿真D.天津：天津工业大学，2018.Li Wenqi.Dynamic analysis and simulation of space four-bar linkage weft insertion mechanism of rapier loomD.Tianjin：Tianjin Polytechnic University，2018.11 汪文煜，周小红.混纬织物喷气引纬工艺探索J.现代纺织技术，2015，23(3):13，6.Wang Wenyu，Zhou Xiaohong.Exploration of jetting weft insertion technology of blended-weft fabricJ.Advanced Textile Technology，2015，23(3):13，6.12 Li Dandan，Yang Jiancheng，Wang Huaqing，et al.3D design and performance analysis of weft insertion mechanism of carbon fiber multi-layer loomC/Proceedings of 2013 3rd International Conference on Mechanical Materials and Manufacturing Engineering(ICMMME 2013).S.l.:Trans Tech Publications Ltd，2013:605-610.13 蔡志城，房家惠，马颜雪，等.段纬自动引纬技术的研发J.产业用纺织品，2021，39(1):2023，43.Cai Zhicheng，Fang Jiahui，Ma Yanxue，et al.Research and development of automatic segment filling insertion technologyJ.Technical Textiles，2021，39(01):2023，43.14 冯伟，王盛，张守伟.挠性剑杆织机打纬机构高速化分析J.机械制造与自动化，2011，40(6):5152，85.Feng Wei，Wang Sheng，Zhang Shouwei.High speed analysis of flexible rapier looms beating-up mechanismJ.Machine Building&Automation，2011，40(6):5152，85.15 蔡锡松，胥光申.基于多眼综丝开口方式的辅助开口系统J.轻工机械，2022，40(6):8187.Cai Xisong，Xu Guangshen.Auxiliary opening system based on multi-eye heald opening modeJ.Light Industry Machinery，2022，40(6):8187.16 董红坤，贺辛亥，钟鹏，等.三维机织物织边造型工艺设计J.纺织学报，2017，38(4):5054.Dong Hongkun，He Xinhai，Zhong Peng，et al.Design of 3-D woven fabric selvage modeling processJ.Journal of Textile Research，2017，38(4):5054.17 蒋宇，李志雄，汤双清，等.三维五向编织复合材料飞轮预制件的工艺设计J.现代制造工程，2009(10):8387.Jiang Yu，Li Zhixiong，Tang Shuangqing，et al.Process design of three dimensional and five directional braided composite flywheelJ.Modern Manufacturing Engineering，2009(10):8387.14 余春平，袁明.肉桂醛的研究与应用进展J.山东化工，2022，51(19):9799.Yu Chunping，Yuan Ming.Recent advances of research and application on cinnamaldehydeJ.Shandong Chemical Industry，2022，51(19):9799.15 于文喜，段华伟，曾少甫，等.山梨酸钾抗菌包装材料研究进展J.包装学报，2019，11(6):716.Yu Wenxi，Duan Huawei，Zeng Shaofu，et al.Advances in antimicrobial packaging materials of potassium sorbateJ.Packaging Journal，2019，11(6):716.16 李凤红，师岩，陈延明，等.一种可降解复合增塑剂改性的PVA及其制备方法:CN111303570AP.2020-06-19.Li Fenghong，Shi Yan，Chen Yanming，et al.A degradable composite plasticizer modified PVA and its preparation method:CN111303570AP.2020-06-19.17 Gao H F，Yang H.Preparation and characterization of cinnamaldehyde/polyvinyl alcohol/silver nanoparticles ternary composite filmsJ.International Journal of Polymer Analysis and Characterization，2021，26(1):2436.18 王帅阳.添加剂对聚乙烯醇/木聚糖复合膜性能的影响D.广州：华南理工大学，2014.Wang Shuaiyang.Influence of additives on properties of poly vinyl alcohol/xylan composite filmsD.Guang Zhou：South China University of Technology，2014.19 Fasihnia S H，Peighambardoust S H，Peighambardoust S J，et al.Development of novel active polypropylene based packaging films containing different concentrations of sorbic acidJ.Food Packaging and Shelf Life，2018，18:8794.20 Ardekani-Zadeh A H，Hosseini S F.Electrospun essential oil-doped chitosan/poly(-caprolactone)hybrid nanofibrous mats for antimicrobial food biopackaging exploitsJ.Carbohydrate Polymers，2019，223.DOI:10.1016/j.carbpol.2019.115108.21 韩笑.抑菌性聚乙烯醇膜对冷鲜猪肉的保鲜研究D.南宁：广西大学，2020.Han Xiao.Study on preservation of cold fresh pork by polyvinyl alcohol filmD.Nanning：Guangxi University，2020.22 Ke J X，Xiao L Y，Yu G X，et al.The study of diffusion kinetics of cinnamaldehyde from corn starch-based film into food simulant and physical properties of antibacterial polymer filmJ.International Journal of Biological Macromolecules，2019，125:642650.(上接第26页)61