CN109165439A - 一种电流体动力学均匀性喷印图案的工艺参数优化方法 - Google Patents
一种电流体动力学均匀性喷印图案的工艺参数优化方法 Download PDFInfo
- Publication number
- CN109165439A CN109165439A CN201810952905.7A CN201810952905A CN109165439A CN 109165439 A CN109165439 A CN 109165439A CN 201810952905 A CN201810952905 A CN 201810952905A CN 109165439 A CN109165439 A CN 109165439A
- Authority
- CN
- China
- Prior art keywords
- representing
- sequence
- jet printing
- electrohydrodynamic
- coefficient
- 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.)
- Granted
Links
- 238000007639 printing Methods 0.000 title claims abstract description 156
- 238000000034 method Methods 0.000 title claims abstract description 81
- 239000007921 spray Substances 0.000 title claims abstract description 76
- 230000008569 process Effects 0.000 title claims abstract description 54
- 238000002474 experimental method Methods 0.000 claims abstract description 72
- 238000005457 optimization Methods 0.000 claims abstract description 13
- 239000000463 material Substances 0.000 claims abstract description 11
- 238000004364 calculation method Methods 0.000 claims abstract description 6
- 238000009825 accumulation Methods 0.000 claims description 58
- 238000002347 injection Methods 0.000 claims description 35
- 239000007924 injection Substances 0.000 claims description 35
- 239000000243 solution Substances 0.000 claims description 34
- 238000005316 response function Methods 0.000 claims description 10
- 230000001186 cumulative effect Effects 0.000 claims description 9
- 238000001802 infusion Methods 0.000 claims description 9
- 230000009471 action Effects 0.000 claims description 6
- 239000000654 additive Substances 0.000 claims description 6
- 230000000996 additive effect Effects 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 3
- 238000010219 correlation analysis Methods 0.000 claims description 3
- 238000011161 development Methods 0.000 claims description 3
- 230000002068 genetic effect Effects 0.000 claims description 3
- 239000002245 particle Substances 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000002922 simulated annealing Methods 0.000 claims description 3
- 230000007246 mechanism Effects 0.000 abstract description 2
- 239000012530 fluid Substances 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 6
- 239000007788 liquid Substances 0.000 description 5
- 230000005684 electric field Effects 0.000 description 4
- 230000005499 meniscus Effects 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000024703 flight behavior Effects 0.000 description 1
- 238000007641 inkjet printing Methods 0.000 description 1
- 239000002082 metal nanoparticle Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002086 nanomaterial Substances 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 229920000642 polymer Polymers 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 239000010409 thin film Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Application Of Or Painting With Fluid Materials (AREA)
- Spray Control Apparatus (AREA)
Abstract
本发明公开了一种电流体动力学均匀性喷印图案的工艺参数优化方法，将计算机优化方法和电流体动力学喷印技术机理相结合，根据少量的实验得到不同工艺参数的电流体动力学喷印图案，采用计算机优化方法进行优化计算得到不同喷印图案的工艺参数，提高图案喷印效率，通过优化计算方法得到符合实验结果的电流体动力学喷印图案厚度模型，并且根据图案喷印厚度给出最佳厚度的工艺参数范围，使得电流体动力学喷印图案厚度更加均匀，不但节省了人力物力，还节省了时间成本，同时提高了流体动力学喷印效率和图案喷印质量。
Description
技术领域
本发明涉及电流体动力学喷印技术领域，具体涉及一种电流体动力学喷印图案的工艺参数优化方法。
背景技术
电流体动力学喷印技术具有设备简单、成本低、效率高、可使用材料广、无需掩模、直接成形等优点，尤其适用于聚合物和金属纳米颗粒等溶液材料的图案喷印，能够制备可穿戴传感器、柔性电子显示器、有机发光二极管、薄膜晶体管、射频识别器件、太阳能电池、电子皮肤、电子报纸等器件。电流体动力学喷印技术是制造微纳米结构与器件的一种新型喷印技术，在微纳制造方面具有巨大的潜能和突出的优势。
电流体动力学喷印技术使用高压电场替代传统喷墨技术中的压电、热气饱或气动等作用，其工作原理为：在基板和喷嘴之间施加电压，在诱导电场力作用下，溶液从喷嘴流出，在喷嘴处形成弯液面，随着电压逐渐升高，电荷在弯液面聚集，电荷间的库仑力导致液体表面产生切向应力，在剪切力的作用下，弯液面在喷嘴顶端形成泰勒锥，随着电场强度增加，库仑力克服液体表面张力，液滴从泰勒锥的顶端射出，形成射流，并在电场力作用下发生破碎，形成液滴，液滴直径通常比喷嘴直径小得多，能够产生亚微米级的分辨率精度。
当电流体动力学喷印过程采用不同的工艺参数时，射流在喷嘴处可能产生两种情况：(1)当射流产生断裂时，形成垂流、微垂流、纺锤体和微纺锤体等喷射模式；(2)当射流连续喷出时，形成锥射流、震荡射流、旋进射流和多射流模式等喷射模式。除了锥射流以外喷射模式的射流飞行行为比较复杂，难以控制射流的飞行轨迹。因此，电流体动力学喷印过程中所使用的具体工艺参数会对喷印图案质量产生重要影响。
电流体动力学喷印技术很难采用一个成熟的理论方法加以借鉴和指导，实现高质量的图案喷印，在实际的图案喷印过程中需要针对不同的工艺参数对喷印溶液进行大量的实验，得到不同工况的实验结果，实验需要的时间长，同时还要耗费大量的材料，因此存在成本高和效率低的问题。
发明内容
针对上述问题，本发明的目的在于提供一种电流体动力学均匀性喷印图案的工艺参数优化方法，将计算机优化方法和电流体动力学喷印技术机理相结合，根据少量的实验得到不同工艺参数的电流体动力学喷印图案，采用计算机优化方法进行优化计算得到不同喷印图案的工艺参数，提高图案喷印效率，并且节省了大量的人力物力，降低成本。
为实现上述目的，本发明提供一种电流体动力学均匀性喷印图案的工艺参数优化方法，其包括以下步骤：
一、通过若干组实验，获取不同材料的喷印溶液、不同工艺参数的电流体动力学喷印设备进行喷印的图案，获取图案上若干个取样点的厚度，并对若干个取样点的厚度进行均值计算；
二、根据若干组的不同工艺参数值及其得到的相应的图案喷印厚度值，构建电流体动力学喷印图案的原始数据序列：
其中，X′(0)表示电流体动力学喷印图案的原始数据，表示若干组实验得到的图案喷印厚度的数据序列，表示喷头直径的数据序列，表示施加电压的数据序列，表示喷射高度的数据序列，表示注射泵流量的数据序列，表示运动平台移动速度的数据序列，表示喷印溶液粘度数据序列，表示在1,2,…,n次实验所得到的图案喷印厚度序列，表示在1,2,…,n次实验喷头直径序列，表示在1,2,…,n次实验施加电压序列，表示在1,2,…,n次实验喷射高度序列，表示在1,2,…,n次实验注射泵流量序列，表示在1,2,…,n次实验运动平台移动速度序列，表示在1,2,…,n次实验喷印溶液粘度序列；
三、对电流体动力学喷印图案的原始数据序列进行数据处理，生成电流体动力学喷印图案的一阶累加生成数据序列：
其中，X′(1)表示生成电流体动力学喷印图案的一阶累加生成数据序列，表示若干组实验得到的图案喷印厚度的一阶累加生成数据序列，表示喷头直径的一阶累加生成数据序列，表示施加电压的一阶累加生成数据序列，表示喷射高度的一阶累加生成数据序列，表示注射泵流量的一阶累加生成数据序列，表示运动平台移动速度的一阶累加生成数据序列，表示喷印溶液粘度的一阶累加生成数据序列，表示在1,2,…,n次实验所得到的图案喷印厚度的一阶累加生成序列，表示在1,2,…,n次实验喷头直径的一阶累加生成序列，表示在1,2,…,n次实验施加电压的一阶累加生成序列，表示在1,2,…,n次实验喷射高度的一阶累加生成序列，表示在1,2,…,n次实验注射泵流量的一阶累加生成序列，表示在1,2,…,n次实验运动平台移动速度的一阶累加生成序列，表示在1,2,…,n次实验喷印溶液粘度的一阶累加生成序列；
四、根据步骤二的原始数据序列和步骤三的一阶累加生成数据序列，对实验数据进行关联分析，并建立多变量一阶累加生成数据的电流体动力学喷印模型：
其中，表示第i次实验得到的图案喷印厚度，表示第i次实验得到的图案喷印厚度的一阶累加生成序列，表示第i次实验喷头直径的一阶累加生成序列，表示第i次实验施加电压的一阶累加生成序列，表示第i次实验喷射高度的一阶累加生成序列，表示第i次实验注射泵流量的一阶累加生成序列，表示第i次实验运动平台移动速度的一阶累加生成序列，表示第i次实验喷印溶液粘度的一阶累加生成序列，ξ表示电流体动力学喷印模型的背景系数，取值范围为0＜ξ＜1，a表示电流体动力学喷印模型发展系数，bi,i＝2,…,7表示电流体动力学喷印模型的驱动系数，c表示电流体动力学喷印模型的线性修正系数，d表示电流体动力学喷印模型的调整作用量；
五、采用最小二乘法对一阶累加生成数据的电流体动力学喷印模型进行参数计算，得到参数序列表达式：
p＝(ATA)-1ATB
其中，p表示由参数序列组成矢量，其表达式为p＝[b2,b3,b4,b5,b6,b7,a,c,d]T，A由原始实验数据一阶累加生成数据和背景系数组成的9×(n-1)阶矩阵，其表达式为：
B表示在第2次，第3次，…，第n次实验分别得到的图案喷印厚度组成的矢量，其表达式为：
六、获取一阶累加生成数据的电流体动力学喷印模型的时间响应函数：
其中，为一阶累加生成数据实验模型的时间响应函数，k＝2,…,n，μ1表示第一常数项，满足μ2表示第二常数项，满足μ3表示第三常数项，满足μ4表示第四常数项，满足
七、根据步骤六中得到的时间响应函数，确定电流体动力学喷印工艺参数的原始数据序列：
其中，为电流体动力学喷印工艺参数优化模型，且k＝2,…,n，e为大于1的整数，满足e＝2,…,k，r为大于1的整数，满足r＝2,…,k-3；
八、根据步骤七中获取的原始数据序列，还原电流体动力学喷印在不同工艺参数下的n组实验的原始数据序列：
九、根据步骤八中不同工艺参数下的n组实验的原始数据序列，建立电流体动力学喷印图案厚度与工艺参数之间的关系模型：
其中f0表示为截距，f1表示喷头直径的影响系数，f2表示施加电压的影响系数，f3表示喷射高度的影响系数，f4表示注射泵流量的影响系数，f5表示运动平台移动速度的影响系数，f6表示喷印溶液粘度的影响系数；
十、根据步骤八中得到的不同工艺参数的n组实验的原始数据序列，获取电流体动力学喷印工艺参数模型的截距和各项系数：
f＝(ZTZ)-1ZTB，
其中f表示由计算得到的电流体动力学喷印工艺参数模型的截距和各项系数的组成的矢量，其表达式为f＝[f0,f1,f2,f3,f4,f5,f6]T，f0表示为截距，f1表示喷头直径的影响系数，f2表示施加电压的影响系数，f3表示喷射高度的影响系数，f4表示注射泵流量的影响系数，f5表示运动平台移动速度的影响系数，f6表示喷印溶液粘度的影响系数，Z由原始实验数据累加生成数据组成的n×7阶矩阵：
B表示在第2次，第3次，…，第n次实验分别得到的图案喷印厚度组成的矢量
十一、根据电流体动力学喷印时输入的工艺参数，获取对应工艺参数的电流体动力学喷印图案的厚度：
h(k)＝f0+f1y1(k)+f2y2(k)+…+f6y7(k)(k＝1,…,n)，
其中h(k)表示电流体动力学喷印厚度，f0表示为截距，f1表示喷头直径的影响系数，f2表示施加电压的影响系数，f3表示喷射高度的影响系数，f4表示注射泵流量的影响系数，f5表示运动平台移动速度的影响系数，f6表示喷印溶液粘度的影响系数，y1(k)表示喷头直径值，y2(k)表示施加电压值，y3(k)表示喷射高度值，y4(k)表示喷射高度值，表示y5(k)运动平台移动速度运动平台移动速度值，y6(k)表示喷印溶液粘度值。
步骤一中通过至少4组实验获取不同材料的喷印溶液、不同工艺参数的电流体动力学喷印设备进行喷印的图案。
电流体动力学喷印模型的背景系数ξ通过遗传算法、粒子群算法或模拟退火算法进行优化。
所述工艺参数包括喷头直径、施加电压、喷射高度、注射泵流量、运动平台移动速度。
本发明的有益效果：通过优化计算方法得到符合实验结果的电流体动力学喷印图案厚度模型，并且根据图案喷印厚度给出最佳厚度的工艺参数范围，使得电流体动力学喷印图案厚度更加均匀，不但节省了人力物力，还节省了时间成本，同时提高了流体动力学喷印效率和图案喷印质量。
附图说明
图1是本发明的流程图。
具体实施方式
下面针对附图对本发明的实施例作进一步说明：
本发明提供的一种电流体动力学喷印图案的工艺参数优化方法，其包括如下步骤：
(1)针对不同材料的喷印溶液，利用电流体动力学喷印设备进行图案喷印，通过控制喷头直径、施加电压、喷射高度、注射泵流量、运动平台移动速度等工艺参数进行若干组实验，每组实验完成电流体动力学图案喷印后，测量图案上若干个取样点的厚度，并对该若干个取样点的厚度进行均值计算，得到每组实验的图案喷印厚度；完成若干组实验，得到若干组实验的图案喷印厚度；
(2)根据若干组的不同工艺参数值及其得到的相应的图案喷印厚度值，构建电流体动力学喷印图案的原始数据序列：
式中，X′(0)表示电流体动力学喷印图案的原始数据，表示若干组实验得到的图案喷印厚度的数据序列，表示喷头直径的数据序列，表示施加电压的数据序列，表示喷射高度的数据序列，表示注射泵流量的数据序列，表示运动平台移动速度的数据序列，表示喷印溶液粘度数据序列，表示在1,2,…,n次实验所得到的图案喷印厚度序列，表示在1,2,…,n次实验喷头直径序列，表示在1,2,…,n次实验施加电压序列，表示在1,2,…,n次实验喷射高度序列，表示在1,2,…,n次实验注射泵流量序列，表示在1,2,…,n次实验运动平台移动速度序列，表示在1,2,…,n次实验喷印溶液粘度序列；
(3)对电流体动力学喷印图案的原始数据序列进行数据处理，生成电流体动力学喷印图案的一阶累加生成数据序列：
式中，X′(1)表示生成电流体动力学喷印图案的一阶累加生成数据序列，表示若干组实验得到的图案喷印厚度的一阶累加生成数据序列，表示喷头直径的一阶累加生成数据序列，表示施加电压的一阶累加生成数据序列，表示喷射高度的一阶累加生成数据序列，表示注射泵流量的一阶累加生成数据序列，表示运动平台移动速度的一阶累加生成数据序列，表示喷印溶液粘度的一阶累加生成数据序列，表示在1,2,…,n次实验所得到的图案喷印厚度的一阶累加生成序列，表示在1,2,…,n次实验喷头直径的一阶累加生成序列，表示在1,2,…,n次实验施加电压的一阶累加生成序列，表示在1,2,…,n次实验喷射高度的一阶累加生成序列，表示在1,2,…,n次实验注射泵流量的一阶累加生成序列，表示在1,2,…,n次实验运动平台移动速度的一阶累加生成序列，表示在1,2,…,n次实验喷印溶液粘度的一阶累加生成序列；
(4)根据实验的原始数据序列和一阶累加生成数据，对实验数据进行关联分析，建立多变量一阶累加生成数据的电流体动力学喷印模型：
式中，表示第i次实验得到的图案喷印厚度，表示第i次实验得到的图案喷印厚度的一阶累加生成序列，表示第i次实验喷头直径的一阶累加生成序列，表示第i次实验施加电压的一阶累加生成序列，表示第i次实验喷射高度的一阶累加生成序列，表示第i次实验注射泵流量的一阶累加生成序列，表示第i次实验运动平台移动速度的一阶累加生成序列，表示第i次实验喷印溶液粘度的一阶累加生成序列，ξ表示电流体动力学喷印模型的背景系数，取值范围为0＜ξ＜1，a表示电流体动力学喷印模型发展系数，bi,i＝2,…,7表示电流体动力学喷印模型的驱动系数，c表示电流体动力学喷印模型的线性修正系数，d表示电流体动力学喷印模型的调整作用量；
(5)采用最小二乘法，对一阶累加生成数据的电流体动力学喷印模型进行参数计算，得到参数序列表达式：
p＝(ATA)-1ATB (2)
式中，p表示由参数序列组成矢量，其表达式为p＝[b2,b3,b4,b5,b6,b7,a,c,d]T，A由原始实验数据一阶累加生成数据和背景系数组成的9×(n-1)阶矩阵，其表达式为：
B表示在第2次，第3次，…，第n次实验分别得到的图案喷印厚度组成的矢量，其表达式为：
(6)计算一阶累加生成数据实验模型的时间响应函数，其表达式如下：
式中，为一阶累加生成数据实验模型的时间响应函数，k＝2,…,n，μ1表示第一常数项，满足μ2表示第二常数项，满足μ3表示第三常数项，满足μ4表示第四常数项，满足
(7)根据一阶累加生成数据实验模型的时间响应函数，确定电流体动力学喷印工艺参数的原始数据序列的表达式为：
其中，为电流体动力学喷印工艺参数优化模型，且k＝2,…,n，e为大于1的整数，满足e＝2,…,k，r为大于1的整数，满足r＝2,…,k-3；
(8)根据步骤(7)，还原电流体动力学喷印喷头直径、施加电压、喷射高度、注射泵流量、运动平台移动速度等工艺参数的n组实验的原始数据序列，其表示为：
(9)根据电流体动力学喷印工艺参数的n组实验的原始数据序列，建立电流体动力学喷印图案厚度与工艺参数之间的关系模型，其表达式为：
式中，f0表示为截距，f1表示喷头直径的影响系数，f2表示施加电压的影响系数，f3表示喷射高度的影响系数，f4表示注射泵流量的影响系数，f5表示运动平台移动速度的影响系数，f6表示喷印溶液粘度的影响系数；
(10)根据步骤(8)中的得到工艺参数的n组实验的原始数据序列，计算电流体动力学喷印工艺参数模型的截距和各项系数，表示为：
f＝(ZTZ)-1ZTB (7)
其中f表示由计算得到的电流体动力学喷印工艺参数模型的截距和各项系数的组成的矢量，其表达式为f＝[f0,f1,f2,f3,f4,f5,f6]T，f0表示为截距，f1表示喷头直径的影响系数，f2表示施加电压的影响系数，f3表示喷射高度的影响系数，f4表示注射泵流量的影响系数，f5表示运动平台移动速度的影响系数，f6表示喷印溶液粘度的影响系数，Z由原始实验数据累加生成数据组成的n×7阶矩阵，其表达式为：
B表示在第2次，第3次，…，第n次实验分别得到的图案喷印厚度组成的矢量，其表达式为：
(11)根据电流体动力学喷印喷头直径、施加电压、喷射高度、注射泵流量、运动平台移动速度等工艺参数的输入值，得到电流体动力学喷印图案的厚度，其表达式为：
h(k)＝f0+f1y1(k)+f2y2(k)+…+f6y7(k)(k＝1,…,n) (8)
式中，h(k)表示电流体动力学喷印厚度，f0表示为截距，f1表示喷头直径的影响系数，f2表示施加电压的影响系数，f3表示喷射高度的影响系数，f4表示注射泵流量的影响系数，f5表示运动平台移动速度的影响系数，f6表示喷印溶液粘度的影响系数，y1(k)表示喷头直径值，y2(k)表示施加电压值，y3(k)表示喷射高度值，y4(k)表示喷射高度值，表示y5(k)运动平台移动速度运动平台移动速度值，y6(k)表示喷印溶液粘度值；
进一步地，所述步骤(1)中，所述若干组实验为至少4组。
进一步地，所述步骤(4)中，对所述背景系数ξ采用遗传算法、粒子群算法、模拟退火算法等优化方法进行优化，进一步提高工艺参数模型的精度。
实施例不应视为对本发明的限制，但任何基于本发明的精神所作的改进，都应在本发明的保护范围之内。
Claims (4)
1.一种电流体动力学均匀性喷印图案的工艺参数优化方法，其特征在于：其包括以下步骤：
一、通过若干组实验，获取不同材料的喷印溶液、不同工艺参数的电流体动力学喷印设备进行喷印的图案，获取图案上若干个取样点的厚度，并对若干个取样点的厚度进行均值计算；
二、根据若干组的不同工艺参数值及其得到的相应的图案喷印厚度值，构建电流体动力学喷印图案的原始数据序列：
其中，X′(0)表示电流体动力学喷印图案的原始数据，表示若干组实验得到的图案喷印厚度的数据序列，表示喷头直径的数据序列，表示施加电压的数据序列，表示喷射高度的数据序列，表示注射泵流量的数据序列，表示运动平台移动速度的数据序列，表示喷印溶液粘度数据序列，表示在1,2,…,n次实验所得到的图案喷印厚度序列，表示在1,2,…,n次实验喷头直径序列，表示在1,2,…,n次实验施加电压序列，表示在1,2,…,n次实验喷射高度序列，表示在1,2,…,n次实验注射泵流量序列，表示在1,2,…,n次实验运动平台移动速度序列，表示在1,2,…,n次实验喷印溶液粘度序列；
三、对电流体动力学喷印图案的原始数据序列进行数据处理，生成电流体动力学喷印图案的一阶累加生成数据序列：
其中，X′(1)表示生成电流体动力学喷印图案的一阶累加生成数据序列，表示若干组实验得到的图案喷印厚度的一阶累加生成数据序列，表示喷头直径的一阶累加生成数据序列，表示施加电压的一阶累加生成数据序列，表示喷射高度的一阶累加生成数据序列，表示注射泵流量的一阶累加生成数据序列，表示运动平台移动速度的一阶累加生成数据序列，表示喷印溶液粘度的一阶累加生成数据序列，表示在1,2,…,n次实验所得到的图案喷印厚度的一阶累加生成序列，表示在1,2,…,n次实验喷头直径的一阶累加生成序列，表示在1,2,…,n次实验施加电压的一阶累加生成序列，表示在1,2,…,n次实验喷射高度的一阶累加生成序列，表示在1,2,…,n次实验注射泵流量的一阶累加生成序列，表示在1,2,…,n次实验运动平台移动速度的一阶累加生成序列，表示在1,2,…,n次实验喷印溶液粘度的一阶累加生成序列；
四、根据步骤二的原始数据序列和步骤三的一阶累加生成数据序列，对实验数据进行关联分析，并建立多变量一阶累加生成数据的电流体动力学喷印模型：
其中，表示第i次实验得到的图案喷印厚度，表示第i次实验得到的图案喷印厚度的一阶累加生成序列，表示第i次实验喷头直径的一阶累加生成序列，表示第i次实验施加电压的一阶累加生成序列，表示第i次实验喷射高度的一阶累加生成序列，表示第i次实验注射泵流量的一阶累加生成序列，表示第i次实验运动平台移动速度的一阶累加生成序列，表示第i次实验喷印溶液粘度的一阶累加生成序列，ξ表示电流体动力学喷印模型的背景系数，取值范围为0＜ξ＜1，a表示电流体动力学喷印模型发展系数，bi,i＝2,…,7表示电流体动力学喷印模型的驱动系数，c表示电流体动力学喷印模型的线性修正系数，d表示电流体动力学喷印模型的调整作用量；五、采用最小二乘法对一阶累加生成数据的电流体动力学喷印模型进行参数计算，得到参数序列表达式：
p＝(ATA)-1ATB
其中，p表示由参数序列组成矢量，其表达式为p＝[b2,b3,b4,b5,b6,b7,a,c,d]T，A由原始实验数据一阶累加生成数据和背景系数组成的9×(n-1)阶矩阵，其表达式为：
B表示在第2次，第3次，…，第n次实验分别得到的图案喷印厚度组成的矢量，其表达式为：
六、获取一阶累加生成数据的电流体动力学喷印模型的时间响应函数：
其中，为一阶累加生成数据实验模型的时间响应函数，k＝2,…,n，μ1表示第一常数项，满足μ2表示第二常数项，满足μ3表示第三常数项，满足μ4表示第四常数项，满足
七、根据步骤六中得到的时间响应函数，确定电流体动力学喷印工艺参数的原始数据序列：
其中，为电流体动力学喷印工艺参数优化模型，且k＝2,…,n，e为大于1的整数，满足e＝2,…,k，r为大于1的整数，满足r＝2,…,k-3；
八、根据步骤七中获取的原始数据序列，还原电流体动力学喷印在不同工艺参数下的n组实验的原始数据序列：
九、根据步骤八中不同工艺参数下的n组实验的原始数据序列，建立电流体动力学喷印图案厚度与工艺参数之间的关系模型：
其中f0表示为截距，f1表示喷头直径的影响系数，f2表示施加电压的影响系数，f3表示喷射高度的影响系数，f4表示注射泵流量的影响系数，f5表示运动平台移动速度的影响系数，f6表示喷印溶液粘度的影响系数；
十、根据步骤八中得到的不同工艺参数的n组实验的原始数据序列，获取电流体动力学喷印工艺参数模型的截距和各项系数：
f＝(ZTZ)-1ZTB，
其中f表示由计算得到的电流体动力学喷印工艺参数模型的截距和各项系数的组成的矢量，其表达式为f＝[f0,f1,f2,f3,f4,f5,f6]T，f0表示为截距，f1表示喷头直径的影响系数，f2表示施加电压的影响系数，f3表示喷射高度的影响系数，f4表示注射泵流量的影响系数，f5表示运动平台移动速度的影响系数，f6表示喷印溶液粘度的影响系数，Z由原始实验数据累加生成数据组成的n×7阶矩阵B表示在第2次，第3次，…，第n次实验分别得到的电流体动力学图案喷印厚度组成的矢量
十一、根据电流体动力学喷印时输入的工艺参数，获取对应工艺参数的电流体动力学喷印图案的厚度：
h(k)＝f0+f1y1(k)+f2y2(k)+…+f6y7(k)(k＝1,…,n)，
其中h(k)表示电流体动力学喷印厚度，f0表示为截距，f1表示喷头直径的影响系数，f2表示施加电压的影响系数，f3表示喷射高度的影响系数，f4表示注射泵流量的影响系数，f5表示运动平台移动速度的影响系数，f6表示喷印溶液粘度的影响系数，y1(k)表示喷头直径值，y2(k)表示施加电压值，y3(k)表示喷射高度值，y4(k)表示喷射高度值，表示y5(k)运动平台移动速度运动平台移动速度值，y6(k)表示喷印溶液粘度值。
2.根据权利要求1所述的一种电流体动力学均匀性喷印图案的工艺参数优化方法，其特征在于：步骤一中通过至少4组实验获取不同材料的喷印溶液、不同工艺参数的电流体动力学喷印设备进行喷印的图案。
3.根据权利要求1所述的一种电流体动力学均匀性喷印图案的工艺参数优化方法，其特征在于：电流体动力学喷印模型的背景系数ξ通过遗传算法、粒子群算法或模拟退火算法进行优化。
4.根据权利要求1所述的一种电流体动力学均匀性喷印图案的工艺参数优化方法，其特征在于：所述工艺参数包括喷头直径、施加电压、喷射高度、注射泵流量、运动平台移动速度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810952905.7A CN109165439B (zh) | 2018-08-21 | 2018-08-21 | 一种电流体动力学均匀性喷印图案的工艺参数优化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810952905.7A CN109165439B (zh) | 2018-08-21 | 2018-08-21 | 一种电流体动力学均匀性喷印图案的工艺参数优化方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109165439A true CN109165439A (zh) | 2019-01-08 |
| CN109165439B CN109165439B (zh) | 2022-11-25 |
Family
ID=64896234
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810952905.7A Active CN109165439B (zh) | 2018-08-21 | 2018-08-21 | 一种电流体动力学均匀性喷印图案的工艺参数优化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109165439B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110370804A (zh) * | 2019-07-24 | 2019-10-25 | 嘉兴学院 | 图案喷印宽度驱动的电流体动力学喷印控制方法 |
| CN110385912A (zh) * | 2019-08-05 | 2019-10-29 | 嘉兴学院 | 均匀间距的点阵结构图案的电流体动力学喷印控制方法 |
| CN111198497A (zh) * | 2019-12-04 | 2020-05-26 | 嘉兴学院 | 一种柔性电子线路在线智能打印烧结控制系统及控制方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120105528A1 (en) * | 2010-11-01 | 2012-05-03 | Alleyne Andrew | High Resolution Sensing and Control of Electrohydrodynamic Jet Printing |
| CN107284025A (zh) * | 2017-08-09 | 2017-10-24 | 嘉兴学院 | 一种电流体动力学曲面喷印设备及其控制方法 |
| CN107379768A (zh) * | 2017-07-31 | 2017-11-24 | 嘉兴学院 | 一种电流体动力学高分辨率喷墨打印的控制方法 |
-
2018
- 2018-08-21 CN CN201810952905.7A patent/CN109165439B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120105528A1 (en) * | 2010-11-01 | 2012-05-03 | Alleyne Andrew | High Resolution Sensing and Control of Electrohydrodynamic Jet Printing |
| CN107379768A (zh) * | 2017-07-31 | 2017-11-24 | 嘉兴学院 | 一种电流体动力学高分辨率喷墨打印的控制方法 |
| CN107284025A (zh) * | 2017-08-09 | 2017-10-24 | 嘉兴学院 | 一种电流体动力学曲面喷印设备及其控制方法 |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110370804A (zh) * | 2019-07-24 | 2019-10-25 | 嘉兴学院 | 图案喷印宽度驱动的电流体动力学喷印控制方法 |
| CN110385912A (zh) * | 2019-08-05 | 2019-10-29 | 嘉兴学院 | 均匀间距的点阵结构图案的电流体动力学喷印控制方法 |
| CN111198497A (zh) * | 2019-12-04 | 2020-05-26 | 嘉兴学院 | 一种柔性电子线路在线智能打印烧结控制系统及控制方法 |
| CN111198497B (zh) * | 2019-12-04 | 2022-09-30 | 嘉兴南湖学院 | 一种柔性电子线路在线智能打印烧结控制系统及控制方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109165439B (zh) | 2022-11-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109165439B (zh) | 一种电流体动力学均匀性喷印图案的工艺参数优化方法 | |
| US10974442B2 (en) | Set-up and method of electrohydrodynamic jet 3D printing based on resultant effect of electric field and thermal field | |
| Martin et al. | Inkjet printing-the physics of manipulating liquid jets and drops | |
| CN107379768B (zh) | 一种电流体动力学高分辨率喷墨打印的控制方法 | |
| CN112122614B (zh) | 一种基于自激发静电场驱动喷射沉积3d打印装置、工作方法及其应用 | |
| CN104191819A (zh) | 喷嘴喷射独立可控的阵列化电流体喷印头及其实现方法 | |
| Cheng et al. | Study on driving waveform design process for multi-nozzle piezoelectric printhead in material-jetting 3D printing | |
| Wang et al. | Deformation characteristics of a single droplet driven by a piezoelectric nozzle of the drop-on-demand inkjet system | |
| Yang et al. | Drop-on-demand electrohydrodynamic printing of high resolution conductive micro patterns for MEMS repairing | |
| Kim et al. | Drop-on-demand hybrid printing using a piezoelectric MEMS printhead at various waveforms, high voltages and jetting frequencies | |
| WO2005012161A1 (ja) | 立体構造物の製造方法および微細立体構造物 | |
| US8807678B2 (en) | Methods of driving an inkjet printing apparatus | |
| Zhou et al. | Experimental study of the piezoelectric drop-on-demand drop formation in a coaxial airflow | |
| CN101623954A (zh) | 整体转印喷墨喷嘴板及制造其的方法 | |
| CN204123788U (zh) | 喷嘴喷射独立可控的阵列化电流体喷印头 | |
| US20220266513A1 (en) | Drop-on-demand printer having optimized nozzle design | |
| Kim et al. | On-demand electrohydrodynamic jetting with meniscus control by a piezoelectric actuator for ultra-fine patterns | |
| CN107670871A (zh) | 一种变直径喷印装置 | |
| CN113226569B (zh) | 控制喷射器的方法和相关系统 | |
| CN110385912B (zh) | 均匀间距的点阵结构图案的电流体动力学喷印控制方法 | |
| Byun et al. | A hybrid inkjet printer utilizing electrohydrodynamic jetting and piezoelectric actuation | |
| CN105773965A (zh) | 一种倒置式单电极电流体三维喷印装置 | |
| KR101527632B1 (ko) | 표면처리 및 ｅｈｄ 젯을 이용한 기판의 미세 패턴 형성 방법 | |
| CN108819218B (zh) | 一种电流体直写喷嘴和控制方法 | |
| Lee et al. | Microelectromechanical systems print heads for industrial printing |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| 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 |