Next Article in Journal
Editorial for the Special Issue on Micromachines for Dielectrophoresis, Volume II
Next Article in Special Issue
Research of Vertical via Based on Silicon, Ceramic and Glass
Previous Article in Journal
Magnetostriction Enhancement in Midrange Modulus Magnetorheological Elastomers for Sensor Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Waveform Design Method for Piezoelectric Print-Head Based on Iterative Learning and Equivalent Circuit Model

School of Mechano-Electronic Engineering, Xidian University, Xi’an 710071, China
*
Author to whom correspondence should be addressed.
Micromachines 2023, 14(4), 768; https://doi.org/10.3390/mi14040768
Submission received: 2 January 2023 / Revised: 25 March 2023 / Accepted: 28 March 2023 / Published: 30 March 2023
(This article belongs to the Special Issue Advanced Packaging for Microsystem Applications, 2nd Edition)

Abstract

:
Piezoelectric print-heads (PPHs) are used with a variety of fluid materials with specific functions. Thus, the volume flow rate of the fluid at the nozzle determines the formation process of droplets, which is used to design the drive waveform of the PPH, control the volume flow rate at the nozzle, and effectively improve droplet deposition quality. In this study, based on the iterative learning and the equivalent circuit model of the PPHs, we proposed a waveform design method to control the volume flow rate at the nozzle. Experimental results show that the proposed method can accurately control the volume flow of the fluid at the nozzle. To verify the practical application value of the proposed method, we designed two drive waveforms to suppress residual vibration and produce smaller droplets. The results are exceptional, indicating that the proposed method has good practical application value.

1. Introduction

Piezoelectric droplet deposition technology has been widely used in many fields, including additive manufacturing [1] and electronic devices [2], because of its low operating cost and high material deposition accuracy. However, new application fields such as printing electrochemical sensors and optical microelectromechanical systems devices require better droplet deposition quality. A typical single-nozzle PPH structure is shown in Figure 1, and it mainly consists of a piezoelectric tube and glass tubes. When a voltage is applied to the electrode, the piezoelectric tube is polarized and causes the glass tube to shift radially. The deformation in the glass tube structure can cause pressure vibrations within the fluid. The pressure drives the liquid to form droplets at the nozzle. A detailed description of the droplet formation process can be found in the literature [3].
The droplet-forming quality is determined by the physical properties of the liquid, the structure of the print-head, and the drive waveform [4]. For a specific fluid and PPH, it is difficult to adjust the structure of the PPH and the physical properties of the liquid. Thus, the appropriate design of the drive waveform is the most convenient way to improve the droplet-forming quality. Currently, several drive waveforms exist for PPH according to different application requirements. The most classic drive waveform is the STDW [5,6,7], which drives the droplet produced with the PPH to have the same radius as that of the nozzle. The jetting ability of the PPH could be improved by adjusting the pulse width of the STDW. However, after the droplet is pushed from the PPH, the residual pressure vibrations in the fluid cannot be quickly dissipated. If the next droplet is jetted before the residual pressure vibrations settle, the resulting droplet properties (radius and velocity) will be different from those of the previous droplet. This can degrade the droplet-forming quality. To suppress residual pressure vibrations, a negative trapezoid pulse, called the DTDW, is added after the STDW [8,9,10]. The effect of the first trapezoid pulse is to produce droplets, and that of the second trapezoid pulse is to suppress residual pressure vibrations. Although the DTDW can suppress the residual pressure vibrations to some extent, some residual pressure vibrations remain. Moreover, some additional complex waveforms can be used to produce small droplets without reducing the nozzle radius [11,12]. In general, these drive waveforms ultimately change the volume flow rate of the liquid at the nozzle of the PPH. Therefore, controlling the volume flow rate of the fluid at the nozzle can improve the droplet-forming quality.
It is optimal to use an iterative learning method to improve the droplet-forming quality of PPHs because the droplet-forming process is executed repeatedly during material deposition. Wassink first applied the control method to PPHs to improve their jet frequencies [13]. Given the lack of a theoretical model for PPHs at the time, he identified a frequency response function as the theoretical model of the PPH based on the self-sensing signal [14]. However, since the self-sensing signal is usually composed of structural vibrations and pressure waves, the accuracy of the identified theoretical model was relatively low. In this study, based on the equivalent circuit model proposed in [15], we devise a new iterative, more practical learning waveform design method. Moreover, the existing iterative learning law for PPHs cannot be applied directly. This is because PPHs have two requirements for the drive waveform. The first requirement is that the drive waveform must have a limited driving cycle, and the second requirement mandates that the starting and ending voltages of the drive waveform must be zero. To satisfy these two requirements, some auxiliary constraints are added to the iterative learning process without changing the convergence. To verify the effect of this method, for the MJ-AL-80 PPH, we designed two types of injection volume flow rates of the fluid at the nozzle to suppress residual vibrations and generate smaller droplets. The optimized waveforms corresponding to the two flow rates are obtained using iterative learning. Finally, residual vibration suppression and smaller droplet-generating experiments were conducted based on the DWS [16]. The experimental results indicate that the proposed method can accurately control the volume flow rate at the nozzle.
The remainder of this paper is organized as follows. In Section 2, the equivalent circuit model and parameter estimation method are presented. In Section 3, we describe the iterative learning process. In Section 4, a schematic of the experimental system is introduced. Moreover, to suppress residual pressure vibrations and produce smaller droplets, two drive waveforms were designed based on the method proposed in this study, and the effectiveness of the method was verified using DWS. Finally, Section 5 summarizes this paper.

2. Model and Desired Volume Flow Rate

2.1. Equivalent Circuit Model

The structure of a circle squeeze-mode print-head MJ-AL-80 (a type of PPH fabricated by MicroFab; “MJ-AL” is the name of a series of PPHs, wherein “80” indicates the diameter of the nozzle in µm) is illustrated in Figure 1. The working principle of this PPH is that the drive waveform polarizes the piezoelectric tube to change the diameter of the pressure chamber. The volume change in the pressure chamber leads to a change in the inner fluid density, which causes pressure vibrations in the glass tube. This pressure vibration propagates to the nozzle and drives the fluid to produce droplets. In this process, the diameter change in the glass tube is very small (nanometer scale). To simplify the modeling process, the change in liquid density caused by the diameter change in the pressure chamber can be equivalent to the change in density caused by the increase or decrease in the inner fluid without a change in diameter.
The fluid injection process inside the PPH contains complex flow characteristics. Although the CFD theory can accurately simulate this process, the calculation process is very time-consuming, and the theoretical model in the differential form is difficult to combine with the control theory. The equivalent circuit model corresponds the fluid pressure and volume flow rate to the voltage and current in the circuit. Using a more mature circuit theory, the fluid injection dynamics can be expressed as a state equation. Not only is the solution efficiency very high, but the state equation can be fleetly combined with the control theory; therefore, the equivalent circuit model of the PPH is prepared in this study. Figure 2 illustrates the equivalent circuit of the MJ-AL-80 print-head. Capacitance C1 represents the compression in the fluid volume at the pressure chamber. The compressibility of the liquid in the right channel is equivalent to capacitance C2. The compressibility of the liquid in the left channel is ignored because the flow resistance at the exit is very small. The flow resistance and inertial force of the left channel and left half of the pressure chamber are equivalent to resistance R1 and inductance L1, respectively. Analogously, the flow resistance and inertial force of the right channel and right half of the pressure chamber are equivalent to resistance R2 and inductance L2, respectively. The flow resistance and inertial force of the nozzle are equivalent to resistance R3 and inductance L3, respectively. The ambient pressure is equivalent to voltage source Us (typically Us = 1 [atm]) [17]. The voltage source Ud is used to replace a pressure constant derived from the equation of state. Capacitance C3 represents the fluid stored by the meniscus surface tension [18]. The current source is is used to simulate the increase or decrease in the fluid inside the pressure chamber.
According to KVL and KCL, the equivalent circuit illustrated in Figure 2 can be described using the following state equation:
{ x ˙ = A x + B u y = C x
A = [ 0 0 0 1 C 1 1 C 1 0 0 0 0 0 1 C 2 1 C 2 0 0 0 0 0 1 C 3 1 L 1 0 0 R 1 L 1 0 0 1 L 2 1 L 2 0 0 R 2 L 2 0 0 1 L 3 1 L 3 0 0 R 3 L 3 ] , B = [ 1 C 1 0 0 0 0 0 0 1 L 1 0 0 0 1 L 3 ] , C = [ 0 0 0 0 0 1 ] T , u = [ U s + U d i s ] T , x = [ u 1 u 2 u 3 i 1 i 2 i 3 ] T
In Equation (1), A is the state matrix, x is the state variable, B is the input matrix, u is the system input, C is the output matrix, and y is the system output.

2.2. Parameter Estimation

According to electroacoustic theory, the resistance, inductance, capacitance, and voltage source in the equivalent circuit shown in Figure 2 represent the viscous force, inertial force, compressibility, and pressure in the fluid, respectively [17]. The parameter estimations of resistance and inductance can be undertaken using Equations (2) and (3), respectively,
R f = 2 c ρ μ l 2 l T π r 6
L f = l ρ π r 2
where c is the speed of sound in the fluid, ρ is the density of the fluid, µ is the viscosity of the fluid, l is the length of the estimated portion of the channel, lT is the total length of the channel, r is the inner radius of the channel, and Rf and Lf are the equivalent resistance and inductance of the estimated portion of the channel, respectively.
As the droplet deposition process is very transient, it can be assumed to be an adiabatic process [19], and the pressure change in the pressure chamber can be described by the state equation (dp = c2) under adiabatic conditions. The state equation of the pressure chamber can be described as
d p = c 2 ρ ( t 0 ) c 2 ρ ( t 0 ) [ V ( t 0 ) + t 0 t q ( t ) d t ] V ( t )
where dp is the increase in pressure, ρ(t0) is the initial density of the liquid in the pressure chamber, q(t) is the volume flow rate of the fluid in the pressure chamber, and V(t) is the volume of the pressure chamber. As the volume change in the pressure chamber is very small, the change in liquid density caused by the diameter change in the pressure chamber can be simulated using the fluid source qs(t). Then, Equation (4) can be written as follows:
d p = c 2 ρ ( t 0 ) c 2 ρ ( t 0 ) ( V ( t 0 ) + t 0 t q ( t ) d t + t 0 t q s ( t ) d t ) V ( t 0 )
where V(t0) is the initial volume of the pressure chamber. To establish the relationship between V(t) and qs(t), let Equations (4) and (5) have the same pressure increment. Therefore, the relationship between V(t) and qs(t) can be described as follows:
q s ( t ) = V ( t ) t
According to Equation (5), the pressure in the pressure chamber can be expressed as the sum of the initial pressure and the pressure increment:
p ( t ) = p ( t 0 ) + c 2 ρ ( t 0 ) V ( t 0 ) + t 0 t q ( t ) d t + t 0 t q s ( t ) d t V ( t 0 ) c 2 ρ ( t 0 )
where p(t0) is the initial pressure in the pressure chamber. In the equivalent circuit illustrated in Figure 2, voltage u1 represents the pressure change in the pressure chamber. According to the Kirchhoff Voltage Law (KVL), u1 can be written as
u 1 ( t ) = U s U d + Q 1 ( t 0 ) + i c 1 d t + i s d t C 1
where ic1 is the current through capacitance C1, is is the current from the current source, and Q1(t0) is the initial charge on capacitor C1. This current source releases the charge to capacitance C1, causing voltage uc1 to increase. The current source absorbs the charge from capacitance C1, causing voltage uc1 to decrease. This process is very similar to that causing the pressure changes in the pressure chamber. Therefore, the parameter estimation method for capacitance and the voltage source Ud in the equivalent circuit can be respectively written as
C f = V ( t 0 ) c 2 ρ ( t 0 )
U d = c 2 ρ ( t 0 )
where Cf is the equivalent capacitance of the liquid in the pressure chamber. A capacitor Cn is equivalent to the Laplace pressure caused by surface tension. To estimate the expression of the capacitor parameter, the average values of the flow-out/-in volumes at the nozzle are roof-shaped [17]. The equivalent capacitance of the surface tension can be described as
C n = V m e n p m e n = π r n 4 3 σ
where Vmen is the roof-shaped volume of the fluid with nozzle radius rn, pmen is the Laplace pressure, σ is the coefficient of surface tension, and Cn is the equivalent capacitance of the meniscus. Based on the above equation, the estimation method of component parameters in the equivalent circuit illustrated in Figure 2 is listed in Table 1. A more detailed method of parameter estimation can be found in Reference [15].

3. Iterative Learning Process

3.1. Iterative Learning Method

As the PPH system is a causal system, the input of the drive waveform at any given time will always affect the output thereafter. It is difficult to obtain a good drive waveform using the classical proportional integral derivative (PID) iterative learning law. However, as the time scale of the drive waveform is limited to tens of microseconds, we can consider the system input vector as the object of iteration. This framework is called the lifted-system framework [20]. The discrete form of state Equation (1) can be written as
y ( k ) = G u ( k ) + d 0
where d 0 = [ C A x ( 0 ) , C A 2 x ( 0 ) , , C A N x ( 0 ) ] T is the initial conditional response of the system, y ( k ) = [ y ( 1 , k ) , y ( 2 , k ) , , y ( N , k ) ] T is the system output vector corresponding to iteration k, u ( k ) = [ u ( 1 , k ) , u ( 2 , k ) , , u ( N , k ) ] T is the system input vector corresponding to iteration k, and G = [ C B 0 0 C A B C B 0 C A N 1 B C A N 2 B C B ] is a Toeplitz matrix. A, B, and C are the discrete system state-transition matrix, input matrix, and output matrix, respectively. For a print-head system, the goal of the iterative learning waveform design method is to find system input u within a given time interval [0, N] so that the system output is consistent with the reference output (desired volume flow rate). According to Equation (12), the error vector e(k) corresponding to iteration k between the system output and reference output can be written as
e ( k ) = y d G u ( k ) d 0
where y d = [ y d ( 1 ) , y d ( 2 ) , , y d ( N ) ] T is the system reference output vector. The quadratic form of the error vector corresponding to iteration k can be written as
E ( k ) = 1 2 e ( k ) T e ( k )
where E(k) is the total error corresponding to iteration k. The purpose of the iterative learning waveform design is to find a system input u such that lim k E ( k ) = 0 . To rapidly reduce the total error E, the steepest descent method is used to design the iterative learning laws [21,22]. The normalized negative gradient vector of the total error can be written as
{ v = [ E ( k ) u ( 0 , k ) , E ( k ) u ( 1 , k ) , , E ( k ) u ( N 1 , k ) ] T Δ u ( k ) = v v 2
where Δu is the normalized negative gradient vector of the total error and v is the negative gradient vector without normalization. Therefore, the learning law of the iterative learning system can be written as
u ( k + 1 ) = u ( k ) + η Δ u ( k )
where η is the iterative step length.

3.2. Optimal Iteration Step Length

In Equation (16), when η is too large, the learning process will not converge, and when η is too small, the convergence rate will be reduced. Combining Equations (13) and (16), the total error at time k + 1 can be written as
E ( k + 1 ) = E ( k ) + 1 2 [ Δ u ( k ) T G T G Δ u ( k ) η 2 2 e ( k ) T G Δ u ( k ) η ]
The second term on the right-hand side of Equation (17) is denoted as the total error increment ΔE(k) at iteration k. The expression for the total error increment ΔE(k) can be written as
Δ E ( k ) = 1 2 Δ u ( k ) T G T G Δ u ( k ) η 2 e ( k ) T G Δ u ( k ) η
According to Equation (18), the total error increment is a quadratic function with the iterative step length η as the independent variable. An iterative step length η must exist in each iteration cycle to minimize the total error increment ΔE. As E(k + 1) = E(k) + ΔE(k), the total error will decrease as fast as possible when the total error increment ΔE is equal to the minimum value in all the iteration steps. According to the extremum formula of the quadratic function, the optimal iterative step length η can be expressed as [21]
η = e ( k ) T G Δ u ( k ) Δ u ( k ) T G T G Δ u ( k )

3.3. Adding to Iterative Learning Constraint

The aforementioned iterative learning process can be used to obtain a good drive waveform. However, the drive waveform of the print-head must meet some requirements in practical applications. The first requirement is that the drive waveform should be loaded within a finite amount of time. The second requirement is that the waveform starting and ending voltages must be zero. To satisfy the first requirement, a weight function, as shown in Figure 3, is used to ensure that the system input is valid for a limited period of time. The iterative learning law described in Equation (16) can be rewritten as
u ( k + 1 ) = u ( k ) + η Δ u ( k ) w T
where W is also a vector defined in n-dimensional real space. This weight function makes iterative learning effective only for a limited amount of time.
The system input in Equation (1) consists of two parts: the first part is a pressure constant, and the second part is the volume flow rate of the fluid in the pressure chamber. The task of iterative learning is to find a change rule of the volume flow rate of the fluid in the pressure chamber to control the volume flow rate of the fluid at the nozzle. The volume flow rate of the fluid in the pressure chamber is equal to the derivative of the drive waveform multiplied by a coefficient. Therefore, the drive waveform can be obtained by integrating the volume flow rate obtained using iterative learning. To satisfy the second requirement, the second term on the right-hand-side of Equation (20) is modified using the following equation:
{ ( m ) P = ( m ) P ( ( m ) N d t ( m ) P d t ) ( m ) P d t > ( m ) N d t ( m ) N = ( m ) N ( ( m ) P d t ( m ) N d t ) ( m ) P d t < ( m ) N d t
where m = ηΔu(k)wT is the weighted iterative learning system input correction, and (m)P and (m)N represent the elements of the m vector that are greater than zero and less than zero, respectively. The purpose of Equation (21) is to set the integral of the system input correction ηΔu(k)wT calculated in each iteration step to zero, which implies that the starting and ending voltages of the drive waveform obtained by integrating the volume flow rate are zero.

3.4. Calculated Waveform with the Iterative Learning Method

Based on the algorithm principle mentioned above, the control input corresponding to the desired injection volume flow rate can be iteratively computed according to the calculation process shown in Figure 4. Firstly, the expected injection volume flow rate at the nozzle is designed and set as the iterative optimization target yd. Secondly, the model output y(k) corresponding to any initial system input u(0) is calculated, and the correction amount of the control input is calculated according to Equations (14)–(21). This process is repeated until the total error E(k) converges to the set threshold. Finally, the corresponding control input is converted into a voltage waveform, which generates the expected injection behavior when driving the PPH and thus optimizes the fluid injection performance.

4. Results and Discussion

4.1. Experimental Setup

We verify the actual effect of droplet formation by optimizing waveform generation. A schematic of the DWS is illustrated in Figure 5. This system mainly consists of a USB camera (MV-VD040SM, Microvision, Redmond, WA, USA), light-emitting diode (LED) strobe, PPH (MJ-AL-80, MicroFab, Plano, TX, USA), droplet watch system controller, pressure controller, and data processing software that runs on a personal computer. A USB camera with an external trigger interface is used to record the droplet images at different times. The LED strobe is used to control the exposure time of the USB camera because the default exposure time of the camera is so long that the edge of the captured droplet image is blurred. To reduce the effective exposure time of the USB camera, the lighting time of the LED strobe is controlled with the droplet watch system controller. When the LED strobe is dark, the USB camera cannot be effectively exposed even during the exposure time, because the incident light is insufficient. The exposure time of the USB camera can be adjusted by setting the lighting time of the LED strobe. To ensure that the edge of the droplet image captured with the USB camera is as clear as possible without reducing the contrast, the lighting time of the LED strobe should be as short as possible. The pressure controller generates negative pressure inside the PPH to balance the weight of the fluid. The data processing software is used to analyze the photos captured with the USB camera and extract motion information of the meniscus. The detailed experimental theory can be found in Reference [16].

4.2. Suppressed Residual Vibration

Residual vibration during droplet deposition is the primary factor limiting the deposition efficiency and quality. For ethanol, the volume flow rate corresponding to STDW is shown in the magenta triangle line (the line is obtained by solving the state equation) in Figure 6a, and the results show that there are multiple peaks.
Except for the first peak, other peaks are redundant. To ensure that the new droplet ejection behavior is not affected by residual oscillation, it is necessary to wait for the residual peak energy to decay to zero before the new droplet ejection, which limits the ejection frequency of microdroplets. Therefore, removing residual peaks is crucial to maximizing injection frequency. In this study, we use a weighted window function (the blue line in Figure 6a) to remove residual peaks and obtain the desired injection volume flow rate (the red circle line in Figure 6a).
The desired injection volume flow rate in Figure 6a is used as the optimization target for iterative learning, and the iterative algorithms are started. Figure 6b shows the injection volume flow rate (solid blue lines) of each iteration step in the iterative learning process, the results show that as the iterations increase, the injection volume flow rate keeps approaching the optimization target, and at the 50th iteration step, the residual peak has been well suppressed (as shown by the magenta circle line in Figure 6b). Figure 6c shows the change in total error, E(k), and iteration step length, η(k), with the iterations. The results show that the total error in the iteration process is strictly convergent (it can be seen from the figure that the total error, E(k), basically converges after 10 iterations, which shows that the iterative algorithm has fast convergence characteristics and can improve the real-time performance of the controller), and the optimal iteration step length is automatically calculated according to the error. Figure 6d shows the optimized drive waveform. Compared with the STDW, the optimized drive waveform has a new negative voltage pulse, which is mainly used to suppress the residual pressure vibrations.
The experimental system described in Section 4.1 is used to verify the actual effects of suppressing residual vibrations. Figure 7a shows the meniscus vibration in the ethanol when the STDW is used. The results show that the vibration period of the meniscus is about 46 μs, and the volume is the largest at around 60, 106, 153, and 200 μs, respectively. This is highly consistent with the prediction results in Figure 6a, which proves that the model established in this study has high accuracy. Figure 7b shows the vibration in the meniscus when the optimized drive waveform is used. The results indicate that the optimized drive waveform presented can effectively suppress the residual vibrations. Moreover, when the STDW is used for MJ-AL-80 PPH, and the frequency is higher than 1 kHz, an unstable ejection phenomenon will occur (as shown in Figure 7c). When the optimized drive waveform is used, the droplet ejection process is still stable at the frequency of 7 kHz (as shown in Figure 7d). This shows that the optimized waveform can significantly improve the droplet ejection frequency.

4.3. Production of Smaller Droplets

Typically, a droplet produced with the PPH has the same diameter as that of the nozzle. However, some special drive waveforms can be used to reduce the droplet diameter without reducing the nozzle diameter. These drive waveforms produce smaller drops by forming a tongue on the meniscus [11]. However, it is difficult to produce smaller droplets when the viscosity of the fluid is too low [12]. This is because it is difficult for low-viscosity fluids to form a stable tongue on the meniscus. However, as shown in Figure 8a, increasing the drive force of a low-viscosity fluid can produce multiple smaller droplets, and the second smaller droplet is clearly formed by trailing tails. One way to eliminate the second/following droplet is to create a new meniscus that pulls the tail back inside the nozzle. To create a new meniscus at the nozzle after the main droplet is formed, the volume flow rate with two volume flow rate peaks shown in Figure 8b is used. The first peak is used to produce smaller droplets, and the second peak is used to form a new meniscus at the nozzle to pull the tail back inside the PPH. Similar to the reference volume flow rate to suppress residual pressure vibration, to ensure the accuracy of calculation, the discretization period of this process is 1 μs. The volume flow rate in Figure 8b is used for iterative learning. Similar to Figure 6, the optimized waveform shown in Figure 8c can be obtained using iterative learning. Figure 8d shows the formation of smaller droplets when the optimized waveform in Figure 8c is used. When the main droplet breaks, a new meniscus is formed to pull the tail back inside the nozzle. The diameter of the PPH nozzle is 80 µm and that of the smaller droplet is approximately 26 µm. The droplet diameter decreases by 67.5%.
Taken together, the experimental results indicate that the proposed method can accurately control the volume flow rate of fluid at the nozzle. It can be used to suppress residual vibrations and produce smaller droplets. Moreover, this method can be applied to various other applications such as controlling the speed and size of the droplets.

5. Conclusions

This paper proposes a more practical waveform design method based on the equivalent circuit model and iterative learning, which can quickly calculate the drive waveform corresponding to the desired volume flow rate. To verify the practical effect of the proposed method, we designed the drive waveforms for suppressing residual pressure vibration and generating smaller droplets, respectively. For the suppression of residual pressure vibration, the experimental results show that the drive waveform designed with the proposed method can well remove the redundant pressure peak, and the injection frequency can be increased by 7 times. For the generation of smaller droplets, the experimental results show that the drive waveform designed with the proposed method can reduce the droplet volume of low-viscosity fluid by 67%. The experimental results indicate that this method can accurately control the volume flow rate of the fluid at the nozzle. However, the surface behavior of the fluid at the nozzle is complex, and the surface behaviors will change with different volume flow rates. For more complex applications, it is difficult to find a corresponding reference volume flow rate. Future studies will focus on this aspect.

Author Contributions

Conceptualization, J.W. and J.H.; methodology, J.W. and C.X.; software, J.W.; validation, C.X., J.P., J.Z. and P.Z.; formal analysis, J.W.; writing—original draft preparation, J.W. and C.X.; project administration, J.H.; funding acquisition, J.W. and J.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by [Natural Science Basic Research Program of Shaanxi] grant number [2023-JC-QN-0407 and 2023-JC-YB-320], [Natural Science Foundation of China] grant number [52035010 and 52275372], [Fundamental Research Funds for the Central Universities] grant number [ZYTS23012]. And the APC was funded by [ZYTS23012].

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

AbbreviationMeaning
PPHPiezoelectric print-heads
STDWSingle trapezoidal drive waveform
DTDWDouble trapezoidal drive waveform
DWSDroplet watch system
KVLKirchhoff’s voltage law
KCLKirchhoff’s current law
PIDproportional integral derivative

References

  1. Kyobula, M.; Adedeji, A.; Alexander, M.R.; Saleh, E.; Wildman, R.; Ashcroft, I.; Gellert, P.R.; Roberts, C.J. 3D inkjet printing of tablets exploiting bespoke complex geometries for controlled and tuneable drug release. J. Control Release 2017, 261, 207–215. [Google Scholar] [CrossRef] [PubMed]
  2. Roh, J.; Kim, H.; Park, M.; Kwak, J.; Lee, C. Improved electron injection in all-solution-processed n-type organic field-effect transistors with an inkjet-printed ZnO electron injection layer. Appl. Surf. Sci. 2017, 420, 100–104. [Google Scholar] [CrossRef]
  3. Bogy, D.B.; Talke, F.E. Experimental and Theoretical Study of Wave Propagation Phenomena in Drop-on-Demand Ink Jet Devices. IBM J. Res. Dev. 1984, 28, 314–321. [Google Scholar] [CrossRef]
  4. Liou, T.-M.; Chan, C.-Y.; Shih, K.-C. Effects of actuating waveform, ink property, and nozzle size on piezoelectrically driven inkjet droplets. Microfluid. Nanofluid. 2009, 8, 575–586. [Google Scholar] [CrossRef]
  5. Kwon, K.-S. Waveform Design Methods for Piezo Inkjet Dispensers Based on Measured Meniscus Motion. J. Microelectromech. Syst. 2009, 18, 1118–1125. [Google Scholar] [CrossRef]
  6. Kwon, K.-S. Experimental analysis of waveform effects on satellite and ligament behavior via in situ measurement of the drop-on-demand drop formation curve and the instantaneous jetting speed curve. J. Micromech. Microeng. 2010, 20, 115005. [Google Scholar] [CrossRef] [Green Version]
  7. Wu, C.-H.; Hwang, W.-S. The effect of the echo-time of a bipolar pulse waveform on molten metallic droplet formation by squeeze mode piezoelectric inkjet printing. Microelectron. Reliab. 2015, 55, 630–636. [Google Scholar] [CrossRef]
  8. Khalate, A.A.; Bombois, X.; Babuška, R.; Wijshoff, H.; Waarsing, R. Performance improvement of a drop-on-demand inkjet printhead using an optimization-based feedforward control method. Control Eng. Pract. 2011, 19, 771–781. [Google Scholar] [CrossRef]
  9. Khalate, A.A.; Bombois, X.; Scorletti, G.; Babuska, R.; Koekebakker, S.; de Zeeuw, W. A Waveform Design Method for a Piezo Inkjet Printhead Based on Robust Feedforward Control. J. Microelectromech. Syst. 2012, 21, 1365–1374. [Google Scholar] [CrossRef]
  10. He, M.-W.; Sun, L.-L.; Hu, K.-Y.; Zhu, Y.-L.; Chen, H.-N. Analysis of DoD inkjet printhead performance for printable electronics fabrication using dynamic lumped element modeling and swarm intelligence based optimal prediction. J. Central South Univ. 2015, 22, 3925–3934. [Google Scholar] [CrossRef]
  11. Rensch, C. Creation of Small Microdrops; MicroFab Technote; MicroFab Technologies: Plano, TX, USA, 2006. [Google Scholar]
  12. Chen, A.U.; Basaran, O.A. A new method for significantly reducing drop radius without reducing nozzle radius in drop-on-demand drop production. Phys. Fluids 2002, 14, L1–L4. [Google Scholar] [CrossRef]
  13. Wassink, M.G.; Bosch, N.J.M.; Bosgra, O.H.; Koekebakker, S. Enabling higher jet frequencies for an inkjet printhead using Iterative Learning Control. In Proceedings of the 2005 IEEE Conference on Control Applications, Toronto, ON, Canada, 28–31 August 2005. [Google Scholar]
  14. Kwon, K.-S.; Kim, W. A waveform design method for high-speed inkjet printing based on self-sensing measurement. Sens. Actuators A Phys. 2007, 140, 75–83. [Google Scholar] [CrossRef]
  15. Wang, J.; Huang, J.; Peng, J. Hydrodynamic response model of a piezoelectric inkjet print-head. Sens. Actuators A Phys. 2019, 285, 50–58. [Google Scholar] [CrossRef]
  16. Kwon, K.-S.; Jang, M.-H.; Park, H.Y.; Ko, H.-S. An inkjet vision measurement technique for high-frequency jetting. Rev. Sci. Instrum. 2014, 85, 065101. [Google Scholar] [CrossRef] [PubMed]
  17. Zhu, Z.M.; Du, G.H.; Gong, X.F. Theory of Acoustics, 3rd ed.; Nanjing University Press: Nanjing, China, 2012. [Google Scholar]
  18. Kim, B.H.; Lee, H.S.; Kim, S.W.; Kang, P.; Park, Y.S. Hydrodynamic responses of a piezoelectric driven MEMS inkjet print-head. Sens. Actuators A Phys. 2014, 210, 131–140. [Google Scholar] [CrossRef]
  19. Poozesh, S.; Akafuah, N.; Saito, K. New criteria for filament breakup in droplet-on-demand inkjet printing using volume of fluid (VOF) method. Korean J. Chem. Eng. 2016, 33, 775–781. [Google Scholar] [CrossRef]
  20. Dijkstra, B.G. Iterative Learning Control with Applications to a Wafer Stage. Ph.D. Thesis, T.U. Delft, Delft, The Netherlands, September 2004. [Google Scholar]
  21. Owens, D.; Rogers, E.; Moore, K.L. Analysis of Linear Iterative Learning Control Schemes Using Repetitive Process Theory. Asian J. Control 2008, 4, 68–89. [Google Scholar] [CrossRef]
  22. Rogers, E.; Owens, D.H.; Werner, H.; Freeman, C.T.; Lewin, P.L.; Kichhoff, S.; Lichtenberg, G. Norm-Optimal It-erative Learning Control with Application to Problems in Accelerator-Based Free Electron Lasers and Rehabilitation Ro-botics. Eur. J. Control 2010, 16, 497–522. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Structure of the MJ-AL-80 PPH.
Figure 1. Structure of the MJ-AL-80 PPH.
Micromachines 14 00768 g001
Figure 2. Equivalent circuit model of the inkjet print-head.
Figure 2. Equivalent circuit model of the inkjet print-head.
Micromachines 14 00768 g002
Figure 3. Weight function controlling the effective interval of iterative learning.
Figure 3. Weight function controlling the effective interval of iterative learning.
Micromachines 14 00768 g003
Figure 4. Principles of the iterative learning method.
Figure 4. Principles of the iterative learning method.
Micromachines 14 00768 g004
Figure 5. Schematic of the experimental system.
Figure 5. Schematic of the experimental system.
Micromachines 14 00768 g005
Figure 6. (a). Desired injection volume flow rate designed to suppress residual vibration. (b). The iterative learning process of the injection volume flow rate that inhibits the residual vibrations (c). The change in total error and optimal step length with iterations. (d). Optimal waveform for suppressing residual pressure vibrations.
Figure 6. (a). Desired injection volume flow rate designed to suppress residual vibration. (b). The iterative learning process of the injection volume flow rate that inhibits the residual vibrations (c). The change in total error and optimal step length with iterations. (d). Optimal waveform for suppressing residual pressure vibrations.
Micromachines 14 00768 g006
Figure 7. (a). Meniscus vibration process corresponding to STDW. (b). Meniscus vibration process corresponding to the optimized waveform. (c). At 2 kHz frequency, the droplet ejection process corresponding to STDW is unstable. (d). At 7 kHz frequency, the droplet ejection process corresponding to the optimized waveform is stable.
Figure 7. (a). Meniscus vibration process corresponding to STDW. (b). Meniscus vibration process corresponding to the optimized waveform. (c). At 2 kHz frequency, the droplet ejection process corresponding to STDW is unstable. (d). At 7 kHz frequency, the droplet ejection process corresponding to the optimized waveform is stable.
Micromachines 14 00768 g007
Figure 8. (a). Process of forming two smaller droplets (ethanol is used as the fluid). (b). Desired volume flow rate to produce smaller droplet. (c). Optimized waveform obtained using iterative learning. (d). Process of generating a smaller droplet using the optimized waveform.
Figure 8. (a). Process of forming two smaller droplets (ethanol is used as the fluid). (b). Desired volume flow rate to produce smaller droplet. (c). Optimized waveform obtained using iterative learning. (d). Process of generating a smaller droplet using the optimized waveform.
Micromachines 14 00768 g008
Table 1. Parameters mapping relation.
Table 1. Parameters mapping relation.
Equivalent Circuit
Components
Mapping ExpressionsParameters for
MJ-AL-80
Capacitance C1, C2, C3 π r 2 L 2 c 2 ρ ( t 0 ) π r 2 L 3 c 2 ρ ( t 0 ) ,   π r o 4 3 σ C1 = 1.4374 × 10−18
C2 = 8.2565 × 10−19
C3 = 1.1888 × 10−16
Voltage source Ud c 2 ρ ( t 0 ) Ud = 9.897216 × 108
Resistance R1, R2, R3 2 c ρ μ ( L 1 + L 2 / 2 ) 2 l T π r 6 , 2 c ρ μ ( L 3 + L 2 / 2 ) 2 l T π r 6 , 2 c ρ μ L o 2 l T π r e 6 R1 = 2.3018 × 1011
R2 = 1.5635 × 1011
R3 = 3.7675 × 1011
Inductance L1, L2, L3 ρ ( L 1 + L 2 / 2 ) π r 2 , ρ ( L 3 + L 2 / 2 ) π r 2 , ρ L o π r e 2 ,R1 = 5.8984 × 107
R2 = 4.0065 × 107
R3 = 3.4867 × 107
Current source is i s = V ( t ) t
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, J.; Xiong, C.; Huang, J.; Peng, J.; Zhang, J.; Zhao, P. Waveform Design Method for Piezoelectric Print-Head Based on Iterative Learning and Equivalent Circuit Model. Micromachines 2023, 14, 768. https://doi.org/10.3390/mi14040768

AMA Style

Wang J, Xiong C, Huang J, Peng J, Zhang J, Zhao P. Waveform Design Method for Piezoelectric Print-Head Based on Iterative Learning and Equivalent Circuit Model. Micromachines. 2023; 14(4):768. https://doi.org/10.3390/mi14040768

Chicago/Turabian Style

Wang, Jianjun, Chuqing Xiong, Jin Huang, Ju Peng, Jie Zhang, and Pengbing Zhao. 2023. "Waveform Design Method for Piezoelectric Print-Head Based on Iterative Learning and Equivalent Circuit Model" Micromachines 14, no. 4: 768. https://doi.org/10.3390/mi14040768

APA Style

Wang, J., Xiong, C., Huang, J., Peng, J., Zhang, J., & Zhao, P. (2023). Waveform Design Method for Piezoelectric Print-Head Based on Iterative Learning and Equivalent Circuit Model. Micromachines, 14(4), 768. https://doi.org/10.3390/mi14040768

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop