CN107512597B - 基于变速率学习的螺杆式多组份物料下料方法 - Google Patents
基于变速率学习的螺杆式多组份物料下料方法 Download PDFInfo
- Publication number
- CN107512597B CN107512597B CN201710905893.8A CN201710905893A CN107512597B CN 107512597 B CN107512597 B CN 107512597B CN 201710905893 A CN201710905893 A CN 201710905893A CN 107512597 B CN107512597 B CN 107512597B
- Authority
- CN
- China
- Prior art keywords
- blanking
- component
- screw
- multiple groups
- groups part
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000000463 material Substances 0.000 title claims abstract description 170
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000002156 mixing Methods 0.000 claims abstract description 21
- 238000009825 accumulation Methods 0.000 claims abstract description 13
- 239000000203 mixture Substances 0.000 claims abstract description 5
- 238000005303 weighing Methods 0.000 claims description 77
- 239000007921 spray Substances 0.000 claims description 10
- 238000005086 pumping Methods 0.000 claims description 9
- 238000002474 experimental method Methods 0.000 claims description 7
- NJPPVKZQTLUDBO-UHFFFAOYSA-N novaluron Chemical compound C1=C(Cl)C(OC(F)(F)C(OC(F)(F)F)F)=CC=C1NC(=O)NC(=O)C1=C(F)C=CC=C1F NJPPVKZQTLUDBO-UHFFFAOYSA-N 0.000 claims description 6
- 238000003756 stirring Methods 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 5
- 238000002360 preparation method Methods 0.000 claims description 3
- 239000004744 fabric Substances 0.000 claims description 2
- 230000004044 response Effects 0.000 claims description 2
- 230000013016 learning Effects 0.000 abstract description 18
- 238000012804 iterative process Methods 0.000 abstract description 3
- 238000005457 optimization Methods 0.000 abstract description 3
- 230000008569 process Effects 0.000 description 13
- 238000010586 diagram Methods 0.000 description 11
- 238000011049 filling Methods 0.000 description 10
- 238000009826 distribution Methods 0.000 description 8
- 238000004519 manufacturing process Methods 0.000 description 8
- 238000006073 displacement reaction Methods 0.000 description 7
- 238000001514 detection method Methods 0.000 description 5
- 239000004615 ingredient Substances 0.000 description 5
- 239000000843 powder Substances 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000004806 packaging method and process Methods 0.000 description 4
- 238000005520 cutting process Methods 0.000 description 3
- 239000002245 particle Substances 0.000 description 3
- 239000000047 product Substances 0.000 description 3
- 239000002994 raw material Substances 0.000 description 3
- 230000006641 stabilisation Effects 0.000 description 3
- 238000011105 stabilization Methods 0.000 description 3
- CDBYLPFSWZWCQE-UHFFFAOYSA-L Sodium Carbonate Chemical compound [Na+].[Na+].[O-]C([O-])=O CDBYLPFSWZWCQE-UHFFFAOYSA-L 0.000 description 2
- 235000013339 cereals Nutrition 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000001276 controlling effect Effects 0.000 description 2
- 239000003814 drug Substances 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000012856 packing Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 239000004593 Epoxy Substances 0.000 description 1
- 244000068988 Glycine max Species 0.000 description 1
- 235000010469 Glycine max Nutrition 0.000 description 1
- 240000002853 Nelumbo nucifera Species 0.000 description 1
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 1
- 244000046052 Phaseolus vulgaris Species 0.000 description 1
- 235000010627 Phaseolus vulgaris Nutrition 0.000 description 1
- 239000004743 Polypropylene Substances 0.000 description 1
- 239000004793 Polystyrene Substances 0.000 description 1
- 239000006004 Quartz sand Substances 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 244000062793 Sorghum vulgare Species 0.000 description 1
- 241000826860 Trapezium Species 0.000 description 1
- 240000008042 Zea mays Species 0.000 description 1
- 235000005824 Zea mays ssp. parviglumis Nutrition 0.000 description 1
- 235000002017 Zea mays subsp mays Nutrition 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000003213 activating effect Effects 0.000 description 1
- 238000012271 agricultural production Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000004566 building material Substances 0.000 description 1
- 210000002421 cell wall Anatomy 0.000 description 1
- 239000004568 cement Substances 0.000 description 1
- 239000011248 coating agent Substances 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 239000002361 compost Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 235000005822 corn Nutrition 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000006052 feed supplement Substances 0.000 description 1
- 239000003337 fertilizer Substances 0.000 description 1
- 235000013305 food Nutrition 0.000 description 1
- 231100001261 hazardous Toxicity 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000003317 industrial substance Substances 0.000 description 1
- 229920000609 methyl cellulose Polymers 0.000 description 1
- 239000001923 methylcellulose Substances 0.000 description 1
- 235000019713 millet Nutrition 0.000 description 1
- 239000000575 pesticide Substances 0.000 description 1
- 231100000614 poison Toxicity 0.000 description 1
- 230000007096 poisonous effect Effects 0.000 description 1
- 229920002239 polyacrylonitrile Polymers 0.000 description 1
- -1 polypropylene Polymers 0.000 description 1
- 229920001155 polypropylene Polymers 0.000 description 1
- 229920002223 polystyrene Polymers 0.000 description 1
- 229920000915 polyvinyl chloride Polymers 0.000 description 1
- 239000004800 polyvinyl chloride Substances 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 235000021067 refined food Nutrition 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 239000002893 slag Substances 0.000 description 1
- 229910000029 sodium carbonate Inorganic materials 0.000 description 1
- 235000017550 sodium carbonate Nutrition 0.000 description 1
- 238000009987 spinning Methods 0.000 description 1
- 238000005406 washing Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
- 230000004584 weight gain Effects 0.000 description 1
- 235000019786 weight gain Nutrition 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B65—CONVEYING; PACKING; STORING; HANDLING THIN OR FILAMENTARY MATERIAL
- B65G—TRANSPORT OR STORAGE DEVICES, e.g. CONVEYORS FOR LOADING OR TIPPING, SHOP CONVEYOR SYSTEMS OR PNEUMATIC TUBE CONVEYORS
- B65G65/00—Loading or unloading
- B65G65/30—Methods or devices for filling or emptying bunkers, hoppers, tanks, or like containers, of interest apart from their use in particular chemical or physical processes or their application in particular machines, e.g. not covered by a single other subclass
- B65G65/34—Emptying devices
- B65G65/40—Devices for emptying otherwise than from the top
- B65G65/46—Devices for emptying otherwise than from the top using screw conveyors
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B65—CONVEYING; PACKING; STORING; HANDLING THIN OR FILAMENTARY MATERIAL
- B65G—TRANSPORT OR STORAGE DEVICES, e.g. CONVEYORS FOR LOADING OR TIPPING, SHOP CONVEYOR SYSTEMS OR PNEUMATIC TUBE CONVEYORS
- B65G65/00—Loading or unloading
- B65G65/30—Methods or devices for filling or emptying bunkers, hoppers, tanks, or like containers, of interest apart from their use in particular chemical or physical processes or their application in particular machines, e.g. not covered by a single other subclass
- B65G65/32—Filling devices
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B65—CONVEYING; PACKING; STORING; HANDLING THIN OR FILAMENTARY MATERIAL
- B65G—TRANSPORT OR STORAGE DEVICES, e.g. CONVEYORS FOR LOADING OR TIPPING, SHOP CONVEYOR SYSTEMS OR PNEUMATIC TUBE CONVEYORS
- B65G2203/00—Indexing code relating to control or detection of the articles or the load carriers during conveying
- B65G2203/04—Detection means
- B65G2203/042—Sensors
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Filling Or Emptying Of Bunkers, Hoppers, And Tanks (AREA)
Abstract
本发明公开了基于变速率学习的螺杆式多组份物料下料方法，基于单次和累积下料误差，控制器通过迭代学习对螺旋输送器进行控制，迭代过程中基于下料误差的变化分别对单次和累积下料误差的学习因子进行动态调整；各螺旋输送器依次动作，在完成一次配方量下料后，打开落料阀，然后在检测到混料斗中的物料达到设定值时，将物料混合均匀后打开推板，将混合料排出。本发明采用距离传感器和搅拌器对下料仓内的物料堆积进行检测和调节，保证物料密实度稳定，学习因子能自动调整优化，与现有技术相比，本装置不需要对参数反复试凑，能快捷地获得收敛性能较好的迭代式，且下料效率高，在迭代学习过程中的落料能被有效利用。
Description
技术领域
本发明涉及定量下料领域，具体涉及一种基于变速率学习的螺杆式多组份物料下料方法。
背景技术
在工农业制造和商品包装中，有大量的粉粒物料，如聚丙烯、聚苯乙烯、聚氯乙烯、轻甲基纤维素、聚丙烯睛、环氧树脂粉末涂料等化工原料，石英砂、水泥等建材原料，洗衣粉等日用化工产品，小米、大豆等谷物豆类农产品，或粉、渣、粒状加工食品，饲料、化肥、农药等农业生产物料，以及粉粒状的保健品、中西药剂、调味品等均需要自动定量包装或者配料制造。
目前我国有很多企业仍然采用手工定量配料或者包装，一方面劳动强度大，速率慢，经济效益差；另一方面，食品、药品等手工定量往往不能满足卫生要求，有毒有害的物料，人工参与定量容易对人体造成伤害。因此对生产企业来说，急需提供成本较低但具有较高速率和准确度的多组份自动定量下料设备或者装置，满足大量的物料定量包装或者配料制造要求。
目前国内外粉粒物料自动定量下料装置常用方法有两种，容积式和称重式。容积式定量依据物料容积进行计量充填或者投料，定量投料迅速，但定量物料质量受到物料密度变化而变化。为提高下料精度，出现了多种调节方法，如申请号为201320001933.3的中国专利，对螺杆采用变频调速，在接近目标值逐渐减慢喂料速度，减少空中落差值；申请号为201310234280.8的中国专利，在纯碱包装机三速变频给料工艺中采用大小螺杆分多阶段下料；申请号为200920248298.2的中国专利考虑到快速下料时难以控制定量而通过先快后慢的方法来减小供料落差的影响；这些非称重式方案的下料终值只能接近期望值，准确度不高。
称重式定量依据所设定的物料一次量进行计量充填或者投料，需要在下料过程中不断称重，根据称重结果反馈控制下料量，由于称重受到下料冲击和空中滞后物料影响较大，组份下料速度和精度都面临很多困难。为了补偿空中物料对计量精度的干扰，很多方案采用提前关闭阀门的技术，如申请号为201410230888.8的中国专利将配料称重过程划分为三个阶段，并在最后一个阶段采用迭代学习控制方式来计算关闭提前控制量，但该方案中迭代学习的学习因子需要通过反复试凑和反馈观察来进行优选，因此需要长时间的实验和调试，而且该方案仅能提高学习完成后的下料精度，在学习过程中的累积下料精度则无法保证。
发明内容
单纯的螺杆式送料器属于容积式定量范畴，容积式定量充填基于容积来计量充填物料的数量，其结构简单，成本低，但定量充填速度稳定性及精度依赖于物料视比重的稳定性，受物料松散程度、颗粒均匀程度、吸湿性等物理化学性质的影响较大。
从原理上来看，目前容积式定量充填可以分为控制充填物料的流量或时间、以及用相同的计量容器量取物料这两种类型。前者常采用控制振动给料机的振动时间或控制螺旋充填机的旋转时间实现定量充填；后者用量杯、量筒或柱塞的定量充填机等。不管采用那一种容积式定量充填，都有一个共同的问题，就是要尽量使物料的视比重稳定。
由于普通容积式本质上是换算式的，无法像称重式一样掌握下料的确切质量，后来虽然出现了结合称重的方案，但由于没有空中量预测而只能依靠下料最后阶段极低的送料速度来保证精度。
为了对空中量进行预测，由于下料过程中空中落料的多少受到输送装置关闭速度、下料口到秤斗料面间落差大小、物料下落形态流率等因素影响，提前关闭输送装置的时间难以通过离线实验一次性确定。在迭代学习控制中，往往需要对学习因子反复试凑，通过观察误差变化来优化调整学习因子的选择。因而，普通的迭代学习控制需要较长时间的反复实验来获取优化的学习因子，这对研发过程中多组份配方实验和多组份原料的快捷制造来说，是无法满足要求的。
为此，本发明对下料装置的下料仓和计量斗进行改进，减小物料的空中落差及形态变化；对迭代学习中下料误差的变化过程进行实时检测并据此自动调整学习因子取值；同时通过在迭代预测中将累积误差作为被控量，从而能快速实现高精度的连续下料，并且在控制中使得螺旋输送器以较高速度运转，避免出现一般螺杆式下料装置在最后阶段运转特别慢的问题。
本发明的技术解决方案是，提供一种基于变速率学习的螺杆式多组份物料下料方法，包括以下步骤：
S1、根据一次量和各配方比例，确定各组份的一次下料量Ws，给各组份的累积下料误差E赋初始值0；将当前组份设定为第一种组份；
S2、对当前组份从下料仓下料，控制器读取称重模块的传感值，记录初始时刻为t0时计量斗的初始重量G0，起动螺旋输送器并以μ·amax加速度加速，当速度达到vH＝λ·vR时保持速度不变；
其中，amax为螺旋输送器的螺杆额定最大加速度，vR为最大速度，μ为(0.5～0.9)之间的加速度系数，λ为(0.85～1.0)之间的速度系数；
S3、当检测到计量斗重量达到(G0+Ws-Wa-0.5·λ·d·ts)时，开始关闭螺旋输送器并以μ·amax加速度开始减速，直至停止；
其中，Wa为空中量预测值，d为螺杆以最大速度运转时螺旋输送器的下料速率，ts为从开始减速时刻至停止时刻之间的时间长度，
S4、等待物料完全下落至计量斗，读取称重模块的传感值，获得当前实际下料量Wr，计算本次下料误差e＝Wr-Ws；
S5、更新累积下料误差E′＝E+e，计算空中量预测值：
Wa′＝αk·Wa+βk·ek+γ·E，
其中，学习因子α、β和γ分别按如下方式进行动态调整：
其中，k大于等于1，sign()为符号函数，α在单次下料误差e大于等于零及小于零两种情形下分别以初始值1.1和0.9为初始值进行迭代，β的初值取为0.7，γ在最初两次取零值且从k等于3开始按上式取值；
S6、迭代，令E＝E′，Wa＝Wa′，ek-2＝ek-1，ek-1＝ek，为下一次下料准备；
S7、更换下料组份，如果全部组份下料完毕，则转下一步，否则，转至步骤2；
S8、打开计量斗底部的落料阀，使得多组份物料组成的一次配方量物料落入混料斗，读取混料斗中料位传感器的状态，若检测到料位超过设定阈值，则控制混料斗中混料器旋转搅拌，将多组份物料混合均匀后，打开混料斗底部的推板，将混合物料从输料管输出；
S9、如果预设下料批次已经完成，则结束下料；否则，将组份设定为第一种组份，转步骤S2。
作为优选，所述步骤S1之前，还有以下步骤：
(T1)通过离线实验，对称重模块和下料仓内的距离传感器进行标定；
(T2)通过控制器的触摸屏进行参数设置，包括一次量、配方表、批量值、下料速率标定的时长Tb与重复次数、稳定称重延时Td；
(T3)对各组份进行下料标定：从0时刻开始按匀加速-恒速-匀减速的调速曲线运转螺旋输送器一定时长Tb，在关闭螺旋输送器Tb时刻及称重稳定后的Tb+Td时刻，分别读取并记录称重模块的重量值Wcb和Wdb；重复多次后，计算本组份的下料速率d＝AVG[Wdb/(Tb-ts)]/λ，空中量初始值Wa＝AVG(Wdb-Wcb)。
作为优选，还包括以下步骤：
步骤S1还包括：按各组份一次下料量Ws的值从大到小排序，确定各组份的下料次序；
步骤S3还包括：当时间t≥t0+0.8Ws/d时，控制器禁止除安全响应之外的所有外部中断；
步骤S8还包括：混料器间断地多次旋转一个小扇角，再次读取料位传感器的状态，若仍然料位超过设定阈值，才进行混料搅拌。
作为优选，还包括以下步骤：将所述下料仓侧壁上安装的两个距离传感器分别指向下料仓内预设的料位高点和料位低点；当两个距离传感器检测到的物料距离超出一设定比例范围，或者发现单位时间下料量波动超过设定阈值如5％后，控制器命令安装于下料仓侧壁的搅拌器动作，改善物料分布。
作为优选，所述下料仓上方设有储料仓和进料泵，所述进料泵为螺杆式进料泵，其后端进料管的出口有一个物料喷头，所述物料喷头为球冠形，其表面分布有圆形小孔，所述进料泵转速按如下式进行控制：
其中，V进0为一设定最大进料速度，l为下料仓当前料位，LM和Lm分别为所预设的最高、最低下料仓料位。
作为优选，在计量斗的上部设一个分料器，所述分料器为上部为圆锥体下部为压扁的锥体结构的沙漏形分料器，其上部为开口形，下部则仅在长度方向的两端有斜坡形喷嘴；所述计量斗面向所述喷嘴的方向上分布有错落的球冠状凸起，下料时物料经分料器喷嘴和凸起落入计量斗中的料堆。
作为优选，所述下料仓底部有一个抽板；所述螺旋输送器包括螺杆箱、输送螺杆、连接器和电机，所述电机外壳通过连接器与螺杆箱相连，位于螺杆箱内的输送螺杆通过轴套与电机轴相连，螺杆箱上表面相对下料仓底部开口处有一进料口，螺杆箱与电机相对的另一端部还连接有一个竖直放置的下料管，下料时打开抽板，物料从从下料仓底部的开口经螺旋输送器到下料管，再经分料器落到计量斗中的料堆上。
作为优选，所述搅拌器包括依次相连的底座、两个支臂、连接两个支臂的支臂转轴、爪手转轴和爪手，搅拌器动作时爪手做螺旋形翻转。
采用本发明下料方法，与现有技术相比，具有以下优点：本发明分别采用距离传感器和机械手形搅拌器对下料仓内的物料堆积形态进行检测和调节，保证物料密实度稳定；又通过在计量斗中设置分料器来减小物料空中落差与冲击量的变化，能帮助减小迭代收敛需要的次数；通过对迭代学习因子的自动优化调整，能减小离线实验量，快速达到下料误差超调小且调节时间短的学习效果。因而，本发明方法能应用于小批量的快速配料，并且通过对下料累积误差的控制，使得迭代预测过程中的物料能有效利用，防止了物料的浪费；同时由于下料过程中螺杆可以保持较高运转速度，因此时提升了下料效率。
附图说明
图1为螺杆式多组份物料下料装置的组成结构图；
图2为螺杆式多组份物料下料装置的外形结构图；
图3为物料下落过程示意图；
图4为储料仓及下料仓局部结构示意图；
图5为下料仓内料位示意图；
图6为下料仓中搅拌器结构示意图；
图7为分料器及计量斗侧壁结构示意图；
图8为计量斗内多组份物料分层示意图；
图9为电机调速曲线；
图10为固定因子迭代学习物料单次下料误差变化图；
图11为学习因子分区示意图。
其中：1、下料仓 2、螺旋输送器 3、计量斗 4、称重模块 5、落料阀 6、混料斗 7、推板 8、输料管 9、控制器 10、储料仓 11、进料泵 12、搅拌器 13、混料器 14、料位传感器15、进料管 16、物料喷头 17、小孔 18、距离传感器 19、料位面 20、分料器 21、喷嘴 22、凸起 23、下料管
30、机架 101、抽板
121、底座 122、支臂 123、支臂转轴 124、爪手转轴 125、爪手
201、螺杆箱 202、输送螺杆 203、连接器 204、电机
具体实施方式
以下结合附图对本发明的优选实施例进行详细描述，但本发明并不仅仅限于这些实施例。本发明涵盖任何在本发明的精神和范围上做的替代、修改、等效方法以及方案。
为了使公众对本发明有彻底的了解，在以下本发明优选实施例中详细说明了具体的细节，而对本领域技术人员来说没有这些细节的描述也可以完全理解本发明。
在下列段落中参照附图以举例方式更具体地描述本发明。需说明的是，附图均采用较为简化的形式且均使用非精准的比例，仅用以方便、明晰地辅助说明本发明实施例的目的。
如图1和图2所示，本发明基于螺杆式多组份物料下料装置，所述装置包括下料仓1、螺旋输送器2、计量斗3、称重模块4、落料阀5、混料斗6和控制器9，其中每种组份的物料都有一组下料仓1和螺旋输送器2对应，常用的组份类别为2～6种，还可以根据需要增加组份类别。作为优选，下料仓1选用直角梯形和矩形组成的料仓形结构，所述下料仓和螺旋输送器为2～6组。
下料仓1底部有一个抽板101，下料时抽板打开，物料从下料仓底部开口处流出。
螺旋输送器2包括螺杆箱201、输送螺杆202、连接器203和电机204，电机204的外壳通过连接器203与螺杆箱201相连，位于螺杆箱201内的输送螺杆202通过轴套与电机204的轴相连；螺杆箱201上表面相对下料仓1底部开口处有一进料口，其与电机相对的另一端部连接到下料管23，下料管23固定在机架30上。
作为优选，下料管23内壁可设螺旋槽，且所述螺旋槽的槽壁在靠近下料管底部的尾段为断续分布。
作为优选，下料管23内可设置一个由两圆锥台组成、其锥顶相对且下锥台底有莲蓬头喷头的部件，使物料均匀下落。
结合图1和图2所示，下料时，控制器9将抽板101打开，物料从下料仓1落入螺旋输送器2的螺杆箱201中，控制器9又命令电机204起动，输送螺杆202随电机一起旋转，将物料输送到端部的下料管23，并从下料管23落入下方的计量斗3中。
机架30作为设备的框架，用来固定和支撑其他各个部件。称重模块4固定在机架30上，计量斗3则活动式扣压在称重模块4上，计量斗3的底部有开口，所述开口的打开与关闭受落料阀5的控制。计量斗3位于下料管23的下部，多个螺旋输送器2相对下料管23及计量斗3的中心呈径向分布。
控制器9采用触摸式操作方式，其触摸屏上有人机界面供进行多组份物料的配方及其他参数的设置，配方包括一次下料的总重量和每个组份占该重量的百分比。控制器9动态读取称重模块4的当前读数，通过控制各动作部件来实现按配方的下料。
混料斗6位于落料阀5下方，且其底部有一个推板7，推板下方连接有一个输料管8，后者将多组份的混合物料输送到包装袋或者生产设备。
作为优选，在混料斗6的侧壁上安装有一个料位传感器14，其内部还有一个混料器15，所述混料器15采用螺旋形桨叶搅拌器。混料斗6的容量是计量斗3的若干如15倍，在完成多个一次量下料后，控制器9读取料位传感器14的状态，若检测到料位超过设定阈值，则控制混料器旋转搅拌，将多种物料混合均匀后，在控制器9的控制下，推板7打开，混合物料从输料管8输出。
图3示意了物料下落过程中料位落差与落料速度对计量斗冲击的变化，物料以初速度v0从螺旋输送器2中落下，螺旋输送器2出口与计量斗3底部的距离为H，随着计量斗中料位h2的增加，空中落差h1将变小。
称重模块检测到的物料质量当量变化可用下式表示：
其中，在t时刻，dm为螺旋输送器2出口的单位时间落料质量(g/s)，v0为物料下落时的初始速度，Δm的物料在落到计量斗时的速度在Δt时间内从速度v1变为0。
从式(1)可以看出，随着空中落差h1的变化，物料对计量斗的冲击也随着改变，因此，计量斗的重量变化是随时间改变的。
另一方面，式(1)中单位时间落料质量当量还受到下料仓1中物料形态分布的影响。
单纯的螺杆式送料器属于容积式定量范畴，本发明采用动态称重对物料下落量进行检测，并通过迭代学习进行物料空中量的预测。结合图4、5、6和7所示，为了减小螺杆输送时物料下落速率的变化，加快学习收敛速度并抑制变量波动范围，本发明通过距离传感器和机械手形搅拌器对下料仓内的物料堆积形态进行检测和调节，使得其下料口上方交替出现动态料拱的形成与坍塌，保证物料密实度和落料形态的稳定。
如图4所示，下料仓1不断出料，当仓内料位降低到一定值时，需要对其进行补料。为此，在下料仓1上方设置一个储料仓10，储料仓10中的物料通过进料泵11和进料管15进入下料仓1。为使得物料颗粒均匀下料，在进料管15的末端出口处设有一个物料喷头16，物料喷头16表面为球冠形，其表面分布有圆形小孔17，小孔孔径根据物料的粒度进行优选。进料泵11采用螺杆式进料泵，其动作由控制器进行控制。
如图5所示，进料泵11在控制器的控制下，使得下料仓内物料顶面的料位保持在预设值LT附近，其转速按下式进行控制：
其中，V进0为一设定最大进料速度，l为下料仓当前料位，LM和Lm分别为所预设的在LT附近的最高、最低料位。
结合图4、6所示，本发明通过距离传感器和搅拌器的检测与动作配合，保证物料的均匀分布。如图4所示，下料仓1内的料位通过安装在其侧壁上的两个小锥角的距离传感器18来进行检测，由于物料喷头16的中心D与下料仓1的下料口B一般不在一条竖直线上，因此，将两个距离传感器18按偏心布置，高低两个距离传感器P、Q分别指向下料仓内预设的料位高点C和料位低点A，其中，C、A分别与D、B处于同一竖直线上。在下料仓1下料过程中，随着料位面19的降低，当距离传感器Q检测到的距离大于QA所对应的距离时，控制器命令进料泵11动作，开始向下料仓1送料；然后，位面升高，当距离传感器P检测到的距离小于PC所对应的距离时，控制器命令进料泵11动作，停止向下料仓1送料。
作为优选，可以在两个距离传感器上增设可调节上下倾角的底座，从而对下料仓内的料层分布进行分辨率更高的检测。
作为优选，对不同下料仓1由于尺寸和外形的区别，两个距离传感器不足以检测到所有物料板结时，可以增加距离传感器的个数并使得它们指向不同料层高度。
如图6所示，本发明通过在下料仓1的侧壁上安装一个搅拌器12来改善物料的分布。搅拌器12包括依次相连的底座121、两个支臂122、连接两个支臂的支臂转轴123、爪手转轴124和爪手125，其中底座121也含有一个旋转轴。
下料过程中，本发明分别通过距离传感器的检测和对单位时间下料率的跟踪来判断下料仓内物料的分布，使得下料仓内的物料面保持近似抛物线面形。结合图4和图6所示，优化距离传感器的布置，使得当物料均匀分布时，距离传感器P与Q检测到的物料距离在一定比例范围内。当物料局部发生板结或稳定的料拱时，距离传感器的读数将超出该范围，如在一定距离内P检测到物料而Q则未检测到物料时，说明此时出现了局部异常。同时，通过称重模块对各下料仓的下料速度进行实时跟踪。当距离传感器检测到上述异常状态或者发现单位时间下料量波动超过设定阈值如5％后，控制器命令搅拌器动作，通过转轴的旋转，其爪手从起点开始经料位高点区域到料位低点区域，做螺旋形翻转，从而破除偶尔形成的板结或料拱，使物料分布恢复均匀。
通过对物料分布的动态检测和控制，使下部仓斗内的物料活化，改善物料的流动，减小物料密实度的波动，从而保证单位时间充填量的稳定。在搅拌器动作动作的同时，暂停下料，且将抽板关上。
结合图3和图7所示，从式(1)可以看出，由于物料空中落差h1的变化，物料对计量斗的冲击也随着改变，造成称重模块单位时间内的重量增加值是变化的。如图7所示，为减小空中落差变化的影响，本发明在计量斗3上部设置一个分料器20，分料器20呈上部为圆锥体下部为压扁的锥体结构的沙漏形分料器；其中上部为开口形，接纳下料仓中的物料；下部则仅在长度方向的两端对称地分布有斜坡形喷嘴21。计量斗3面向喷嘴21的方向上分布有错落的球冠状凸起22，作为优选，凸起直径为0.2～0.6毫米或为所落物料直径的2～3倍。
通过分料器的作用，物料下落分为三阶段，第一阶段从下料仓底部的开口经螺旋输送器到下料管，第二阶段从下料管下落到分料器，第三阶段从分料器喷嘴到计量斗中的料堆。其中，前两个阶段的物料冲击是近似不变的，第三阶段，由于分料器和计量斗壁上错落分布凸起的作用，物料颗粒冲击计量斗中物料面的速度已经大大降低，从分料器喷嘴到达计量斗中不同高度料堆面的冲击力区别非常小，从而为控制器的迭代预测提供了条件。
图8示意了4种组份下料时计量斗内的物料分布。
传统的迭代学习采用固定学习因子，且不考虑累积误差，如申请号为201410230888.8的中国专利中给料机关闭提前量的迭代式为：
uk+1＝uk+q·ek。
采用固定因子的迭代学习，对下料过程的空中量进行预测，图10示意了迭代过程中物料单次下料误差的变化，其中，横坐标为下料次数，纵坐标为每次的下料相对误差。从图中可以看出，图10a对应的下料误差超调大，收敛太慢；而图10b对应的下料超调量合适，前几次收敛趋势较快，但后段稳定慢，过渡时间太长。
由于传统迭代学习需要对学习因子进行试凑，经过反复实验并基于操作经验才能摸索出较好的参数。因此，本发明通过对下料迭代过程的观察和分析，采用变速率的迭代学习对下料的空中量进行预测，通过对相邻两次下料误差的比较和连续三次下料误差的变化范围对比，分别对迭代学习中单次下料误差和累积下料误差的学习因子进行动态调整，同时根据下料过多还是不足分别调整前一次预测量的学习因子。
如图11为根据相邻两次下料误差的比较对单次下料误差的学习因子进行调整的分区示意图。如图所示，横轴为x，曲线a和d对应包络线为曲线b和c则对应包络线为为使得下料误差迅速接近零且尽早收敛，将相邻两次下料误差的相对关系分成四个区域，分别为曲线a的下侧、曲线a和横轴之间、曲线c和横轴之间及曲线c上侧这四个区域。
基于变速率迭代学习，控制器动态预测物料下落的空中量，结合图9所示，本发明采用以下步骤进行下料控制：
(1)根据一次量和各配方比例，确定各组份的一次下料量Ws，给各组份的累积下料误差E赋初始值0；将当前组份设定为第一种组份；
(2)对当前组份下料，控制器读取称重模块的传感值，记录计量斗的初始重量G0，起动螺旋输送器并以μ·amax加速度加速，当速度达到vH＝λ·vR时保持速度不变；
其中，amax为螺旋输送器的螺杆额定最大加速度，vR为最大速度，μ为(0.5～0.9)之间的加速度系数，λ为(0.85～1.0)之间的速度系数；
(3)当检测到计量斗重量达到(G0+Ws-Wa-0.5·λ·d·ts)时，开始关闭螺旋输送器并以μ·amax加速度开始减速，直至停止；
其中，d为螺杆以最大速度运转时螺旋输送器的下料速率，ts为从开始减速时刻t2开始至停止时刻t3之间的时间长度：
(4)等待物料完全下落至计量斗，读取称重模块的传感值，获得当前实际下料量Wr，计算本次下料误差e＝Wr-Ws；
(5)更新累积下料误差E′＝E+e，计算空中量预测值：
Wa′＝αk·Wa+βk·ek+γ·E，
其中，学习因子α、β和γ分别按如下方式进行动态调整：
其中，k大于等于1，sign()为符号函数，α在单次下料误差e大于等于零及小于零两种情形下分别以初始值1.1和0.9为初始值进行迭代，β的初值取为0.7，γ在最初两次取零值且从k等于3开始按上式取值；
(6)迭代，令E＝E′，Wa＝Wa′，ek-2＝ek-1，ek-1＝ek，为下一次下料准备；
(7)更换下料组份，如果全部组份下料完毕，则转下一步，否则，转至步骤2；
(8)打开落料阀，使得多组份物料组成的一次配方量物料落入混料斗，读取料位传感器的状态，若检测到料位超过设定阈值，则控制混料器旋转搅拌，将多组份物料混合均匀后，打开推板，将混合物料从输料管输出；
(9)如果预设下料批次已经完成，则结束下料；否则，将组份设定为第一种组份，转步骤2。
下料期间，控制器还通过对距离传感器和称重模块信号的计算分析，对下料仓内的料堆形态进行实时检测，若发现非正常下料，则及时命令机械手形搅拌器动作，保障下料时的整体均匀分布和密实度稳定。
在连续下料之前，还要进行以下操作：
(i)通过离线实验，对称重模块和距离传感器进行标定；
(ii)通过控制器的触摸屏进行参数设置，包括一次量、配方表、批量值、下料速率标定的时长Tb与重复次数、稳定称重延时Td；
(iii)对各组份进行下料标定：从0时刻开始按或参考图9所示的调速曲线运转螺旋输送器一定时长Tb，在关闭螺旋输送器Tb时刻及称重稳定后的Tb+Td时刻，分别读取并记录称重模块的重量值Wcb和Wdb；重复多次后，计算本组份的下料速率d＝AVG[Wdb/(Tb-ts)]/λ，空中量初始值Wa＝AVG(Wdb-Wcb)。
应用本发明装置进行下料，不需要依赖人工经验来对学习因子进行调整，控制器能够根据下料误差的变化自动对其进行优化，因此，能快捷地获得收敛性能较好的迭代式，适用于需要快速适应的场合，如研发过程中多组份配方实验和多组份原料的快捷制造等。而且，相比于其他迭代学习，本装置不需要将迭代学习过程中的物料废弃，而是可以直接应用于后续生产，所以又适用于小批量的快速配料下料。
以上所述的实施方式，并不构成对该技术方案保护范围的限定。任何在上述实施方式的精神和原则之内所作的修改、等同替换和改进等，均应包含在该技术方案的保护范围之内。
Claims (8)
1.基于变速率学习的螺杆式多组份物料下料方法，包括以下步骤：
S1、根据一次量和各配方比例，确定各组份的一次下料量Ws，给各组份的累积下料误差E赋初始值0；将当前组份设定为第一种组份；
S2、对当前组份从下料仓下料，控制器读取称重模块的传感值，记录初始时刻为t0时计量斗的初始重量G0，起动螺旋输送器并以μ·amax加速度加速，当速度达到vH＝λ·vR时保持速度不变；
其中，amax为螺旋输送器的螺杆额定最大加速度，vR为最大速度，μ为0.5～0.9之间的加速度系数，λ为0.85～1.0之间的速度系数；
S3、当检测到计量斗重量达到(G0+Ws-Wa-0.5·λ·d·ts)时，开始关闭螺旋输送器并以μ·amax加速度开始减速，直至停止；
其中，Wa为空中量预测值，d为螺杆以最大速度运转时螺旋输送器的下料速率，ts为从开始减速时刻至停止时刻之间的时间长度，ts＝λ·vR/μ·amax；
S4、等待物料完全下落至计量斗，读取称重模块的传感值，获得当前实际下料量Wr，计算本次下料误差e＝Wr-Ws；
S5、更新累积下料误差E′＝E+e，计算空中量预测值：
Wa′＝αk·Wa+βk·ek+γ·E，
其中，学习因子α、β和γ分别按如下方式进行动态调整：
αk记为
其中，k大于等于1，sign()为符号函数，α在单次下料误差e大于等于零及小于零两种情形下分别以初始值1.1和0.9为初始值进行迭代，β的初值取为0.7，γ在最初两次取零值且从k等于3开始按上式取值；
S6、迭代，令E＝E′，Wa＝Wa′，ek-2＝ek-1，ek-1＝ek，为下一次下料准备；
S7、更换下料组份，如果全部组份下料完毕，则转下一步，否则，转至步骤2；
S8、打开计量斗底部的落料阀，使得多组份物料组成的一次配方量物料落入混料斗，读取混料斗中料位传感器的状态，若检测到料位超过设定阈值，则控制混料斗中混料器旋转搅拌，将多组份物料混合均匀后，打开混料斗底部的推板，将混合物料从输料管输出；
S9、如果预设下料批次已经完成，则结束下料；否则，将组份设定为第一种组份，转步骤S2。
2.根据权利要求1所述的基于变速率学习的螺杆式多组份物料下料方法，其特征在于，所述步骤S1之前，还有以下步骤：
T1、通过离线实验，对称重模块和下料仓内的距离传感器进行标定；
T2、通过控制器的触摸屏进行参数设置，包括一次量、配方表、批量值、下料速率标定的时长Tb与重复次数、稳定称重延时Td；
T3、对各组份进行下料标定：从0时刻开始按匀加速-恒速-匀减速的调速曲线运转螺旋输送器一定时长Tb，在关闭螺旋输送器Tb时刻及称重稳定后的Tb+Td时刻，分别读取并记录称重模块的重量值Wcb和Wdb；重复多次后，计算本组份的下料速率d＝AVG[Wdb/(Tb-ts)]/λ，空中量初始值Wa＝AVG(Wdb-Wcb)。
3.根据权利要求1所述的基于变速率学习的螺杆式多组份物料下料方法，其特征在于，还包括以下步骤：
步骤S1还包括：按各组份一次下料量Ws的值从大到小排序，确定各组份的下料次序；
步骤S3还包括：当时间t≥t0+0.8Ws/d时，控制器禁止除安全响应之外的所有外部中断；
步骤S8还包括：混料器间断地多次旋转一个小扇角，再次读取料位传感器的状态，若仍然料位超过设定阈值，才进行混料搅拌。
4.根据权利要求1所述的基于变速率学习的螺杆式多组份物料下料方法，其特征在于，还包括以下步骤：
将所述下料仓侧壁上安装的两个距离传感器分别指向下料仓内预设的料位高点和料位低点；当两个距离传感器检测到的物料距离超出一设定比例范围，或者发现单位时间下料量波动超过设定阈值如5％后，控制器命令安装于下料仓侧壁的搅拌器动作，改善物料分布。
5.根据权利要求1所述的基于变速率学习的螺杆式多组份物料下料方法，其特征在于，所述下料仓上方设有储料仓和进料泵，所述进料泵为螺杆式进料泵，其后端进料管的出口有一个物料喷头，所述物料喷头为球冠形，其表面分布有圆形小孔，所述进料泵转速按如下式进行控制：
其中，V进0为一设定最大进料速度，l为下料仓当前料位，LM和Lm分别为所预设的最高、最低下料仓料位。
6.根据权利要求1所述的基于变速率学习的螺杆式多组份物料下料方法，其特征在于：在计量斗的上部设一个分料器，所述分料器为上部为圆锥体下部为压扁的锥体结构的沙漏形分料器，其上部为开口形，下部则仅在长度方向的两端有斜坡形喷嘴；所述计量斗面向所述喷嘴的方向上分布有错落的球冠状凸起，下料时物料经分料器喷嘴和凸起落入计量斗中的料堆。
7.根据权利要求6所述的基于变速率学习的螺杆式多组份物料下料方法，其特征在于：所述下料仓底部有一个抽板；所述螺旋输送器包括螺杆箱、输送螺杆、连接器和电机，所述电机外壳通过连接器与螺杆箱相连，位于螺杆箱内的输送螺杆通过轴套与电机轴相连，螺杆箱上表面相对下料仓底部开口处有一进料口，螺杆箱与电机相对的另一端部还连接有一个竖直放置的下料管，下料时打开抽板，物料从下料仓底部的开口经螺旋输送器到下料管，再经分料器落到计量斗中的料堆上。
8.根据权利要求4所述的基于变速率学习的螺杆式多组份物料下料方法，其特征在于：所述搅拌器包括依次相连的底座、两个支臂、连接两个支臂的支臂转轴、爪手转轴和爪手，搅拌器动作时爪手做螺旋形翻转。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710905893.8A CN107512597B (zh) | 2017-09-19 | 2017-09-19 | 基于变速率学习的螺杆式多组份物料下料方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710905893.8A CN107512597B (zh) | 2017-09-19 | 2017-09-19 | 基于变速率学习的螺杆式多组份物料下料方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107512597A CN107512597A (zh) | 2017-12-26 |
| CN107512597B true CN107512597B (zh) | 2019-04-05 |
Family
ID=60726274
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710905893.8A Active CN107512597B (zh) | 2017-09-19 | 2017-09-19 | 基于变速率学习的螺杆式多组份物料下料方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107512597B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108414040A (zh) * | 2018-03-23 | 2018-08-17 | 四川大学 | 基于步进电机的中药粉末精密计量仪 |
| KR20220049407A (ko) | 2020-10-14 | 2022-04-21 | 삼성전자주식회사 | 디스플레이 장치 및 그 제어 방법 |
| CN112783126B (zh) * | 2020-12-31 | 2022-03-22 | 浙江中控技术股份有限公司 | 一种基于dcs的流体物料输送提前量计算方法及流体物料输送方法 |
| CN113580483A (zh) * | 2021-06-21 | 2021-11-02 | 安徽美博智能科技有限公司 | 一种注塑加工用原料颗粒储存分配装置 |
| CN114879767A (zh) * | 2022-05-31 | 2022-08-09 | 太仓摩丹卡勒多尼塑料机械有限公司 | 基于调速电机和重量反馈的塑料颗粒给料重量的控制方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2271435A1 (en) * | 1996-11-12 | 1998-05-22 | Voest-Alpine Industrieanlagenbau Gmbh | Device for discharging or mixing lump materials |
| JPH11165832A (ja) * | 1997-12-05 | 1999-06-22 | Babcock Hitachi Kk | スクリュー移動式のホッパ |
| CN101019081A (zh) * | 2003-12-05 | 2007-08-15 | 贝朗医疗有限公司 | 药物配制系统和方法 |
| JP2014019087A (ja) * | 2012-07-20 | 2014-02-03 | Handy Techno Kk | 木質合成建材製造用ペレットの製造装置及び製造方法 |
| CN204162108U (zh) * | 2014-10-28 | 2015-02-18 | 三门峡强芯铸造材料有限公司 | 宝珠砂自动配料装置 |
| CN206068976U (zh) * | 2016-08-26 | 2017-04-05 | 连云港东睦新材料有限公司 | 粉末冶金用辅料自动配料机 |
-
2017
- 2017-09-19 CN CN201710905893.8A patent/CN107512597B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2271435A1 (en) * | 1996-11-12 | 1998-05-22 | Voest-Alpine Industrieanlagenbau Gmbh | Device for discharging or mixing lump materials |
| JPH11165832A (ja) * | 1997-12-05 | 1999-06-22 | Babcock Hitachi Kk | スクリュー移動式のホッパ |
| CN101019081A (zh) * | 2003-12-05 | 2007-08-15 | 贝朗医疗有限公司 | 药物配制系统和方法 |
| JP2014019087A (ja) * | 2012-07-20 | 2014-02-03 | Handy Techno Kk | 木質合成建材製造用ペレットの製造装置及び製造方法 |
| CN204162108U (zh) * | 2014-10-28 | 2015-02-18 | 三门峡强芯铸造材料有限公司 | 宝珠砂自动配料装置 |
| CN206068976U (zh) * | 2016-08-26 | 2017-04-05 | 连云港东睦新材料有限公司 | 粉末冶金用辅料自动配料机 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107512597A (zh) | 2017-12-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107673083B (zh) | 基于变速率学习的螺杆式物料下料装置及其控制器 | |
| CN107715727B (zh) | 螺杆式多组份物料配料装置及其控制器 | |
| CN107512597B (zh) | 基于变速率学习的螺杆式多组份物料下料方法 | |
| CN107684846B (zh) | 直落式多组份物料下料方法 | |
| CN107694469B (zh) | 基于变速率学习的直落式多组份物料配料方法 | |
| CN107572016B (zh) | 直落式多组份物料下料装置及其控制器 | |
| CN107661728B (zh) | 基于变速率学习的直落式物料配料装置及其控制器 | |
| CN107720311B (zh) | 基于神经网络的螺杆失重式物料下料机及其控制器 | |
| CN107684847B (zh) | 螺杆式多组份物料配料方法 | |
| CN107601083B (zh) | 基于神经网络的直落失重式物料下料方法 | |
| CN107640609B (zh) | 基于机器学习的螺杆式物料配料机控制器 | |
| CA2327599C (en) | Dispensing measured quantities of materials for mixing into a larger batch | |
| CN107697660B (zh) | 基于机器学习的螺杆式物料配料机控制方法 | |
| CN110697438B (zh) | 基于神经网络的直落失重式物料下料机控制器 | |
| CN105600500B (zh) | 一种自动计重喂料装置 | |
| CN208882221U (zh) | 一种有机质肥料的生产系统 | |
| CN106475001B (zh) | 一种小型配肥机 | |
| CN108002062B (zh) | 基于神经网络的螺杆失重式物料下料方法 | |
| CN209080205U (zh) | 一种用于婴幼儿辅食定量下料装置 | |
| RU131478U1 (ru) | Многокомпонентный микродозатор | |
| JPS6117025A (ja) | 粉粒体の計量装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20240929 Address after: 230000 B-1015, wo Yuan Garden, 81 Ganquan Road, Shushan District, Hefei, Anhui. Patentee after: HEFEI MINGLONG ELECTRONIC TECHNOLOGY Co.,Ltd. Country or region after: China Address before: 310018, No. 258, source street, Xiasha Higher Education Park, Hangzhou, Zhejiang Patentee before: China Jiliang University Country or region before: China |