<?xml-model href='http://www.tei-c.org/release/xml/tei/custom/schema/relaxng/tei_all.rng' schematypens='http://relaxng.org/ns/structure/1.0'?><TEI xmlns="http://www.tei-c.org/ns/1.0">
	<teiHeader>
		<fileDesc>
			<titleStmt><title level='a'>Sequential Quadratic Programming Iterative Learning Control for a Roll-to-Roll Manufacturing Process</title></titleStmt>
			<publicationStmt>
				<publisher>ASME</publisher>
				<date>10/01/2025</date>
			</publicationStmt>
			<sourceDesc>
				<bibl> 
					<idno type="par_id">10679081</idno>
					<idno type="doi">10.1115/1.4069342</idno>
					<title level='j'>ASME Letters in Dynamic Systems and Control</title>
<idno>2689-6117</idno>
<biblScope unit="volume">5</biblScope>
<biblScope unit="issue">4</biblScope>					

					<author>Christopher Martin</author><author>Wei Li</author><author>Dongmei Chen</author>
				</bibl>
			</sourceDesc>
		</fileDesc>
		<profileDesc>
			<abstract><ab><![CDATA[Roll-to-roll (R2R) mechanical dry transfer is an enabling technology for high-throughput, environmentally friendly fabrication of advanced thin-film devices. However, precise control is required to ensure high-quality transfer, presenting a significant challenge due to nonlinear peeling dynamics, abrupt material property changes, and input constraints. This study proposes a sequential quadratic programming iterative learning control (SQP ILC) approach to regulate the R2R mechanical dry transfer process. The method leverages the system's iterative structure to improve the performance across successive transfer tasks while rigorously accounting for nonlinear dynamics and input constraints. Experimental validation on a lab-scale testbed and a case study on the chemical vapor deposition-grown graphene transfer show that the SQP ILC significantly improves transfer quality with minimal online computation, making it a scalable solution for industrial R2R applications.</p>]]></ab></abstract>
		</profileDesc>
	</teiHeader>
	<text><body xmlns="http://www.tei-c.org/ns/1.0" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:xlink="http://www.w3.org/1999/xlink">
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="1.">INTRODUCTION</head><p>Roll-to-roll (R2R) manufacturing is a class of continuous industrial processes that fabricate thin-film devices by moving flexible substrates through a series of rollers. R2R techniques enable efficient, high-throughput production of advanced technologies such as lithium-ion batteries <ref type="bibr">[1]</ref>, flexible solar cells <ref type="bibr">[2]</ref>, and fuel cell membranes <ref type="bibr">[3]</ref>. A key step in these processes is transferring functional materials or devices from a donor to a target substrate for device integration. Traditional batch-style transfer methods are inefficient <ref type="bibr">[4]</ref>, and while R2R-compatible methods exist, many have limited throughput <ref type="bibr">[5]</ref> or use hazardous chemical etchants that destroy the donor substrate <ref type="bibr">[6]</ref>.</p><p>Recent advancements in R2R mechanical dry transfer provide an efficient, eco-friendly alternative to liquid-based or discrete transfer methods <ref type="bibr">[7]</ref>, <ref type="bibr">[8]</ref>. As shown in Figure <ref type="figure">1</ref>, the process involves laminating a 2D material or device between donor and target substrates, peeling them apart so the material adheres to the target, and rewinding the functional material or device for further processing and the donor substrate for recycling. This process can be repeated for a continuous series of 2D devices or material samples.</p><p>Successful R2R dry transfer requires precise control of peeling conditions, particularly web tensions and speeds <ref type="bibr">[7]</ref>, <ref type="bibr">[8]</ref>, <ref type="bibr">[9]</ref>. Maintaining optimal web tension is critical but challenging due to nonlinear peeling dynamics and abrupt material property changes from device patterning <ref type="bibr">[10]</ref>, <ref type="bibr">[11]</ref>-challenges absent in typical R2R processes without peeling. Furthermore, control design must account for practical constraints such as motor torque limits <ref type="bibr">[11]</ref> and online computational capacity. Alongside these challenges, R2R dry transfer has an iterative structure that can be leveraged to improve control. Each 2D sample or device transfer represents one "transfer task" iteration, meaning data from previous iterations can be used to improve performance in future transfers within the continuous manufacturing line.</p><p>Researchers have tried to develop controllers for the R2R mechanical dry transfer process. Tuned PI controllers ignore the nonlinear process model <ref type="bibr">[10]</ref>, <ref type="bibr">[12]</ref>, while linear differential inclusion (LDI)-based optimal controllers are overly conservative <ref type="bibr">[13]</ref>, <ref type="bibr">[14]</ref>, <ref type="bibr">[15]</ref>. Both methods neglect input constraints. Nonlinear MPC offers high performance but incurs significant online computational costs <ref type="bibr">[11]</ref>. In addition, these existing controllers for R2R mechanical dry transfer do not take advantage of its iterative structure. Thus, a control method is needed that incorporates the nonlinear process model, considers input constraints, has minimal online computation, and uses previous task data to improve performance.</p><p>This last consideration motivates an iterative learning control (ILC) approach. ILC, illustrated in Figure <ref type="figure">2</ref>, improves system performance by adjusting control inputs based on measurements from previous task iterations. The most common form, classical ILC, uses a constant learning gain to adjust the control based on the tracking error from the previous task but does not account for nonlinear models or input constraints <ref type="bibr">[16]</ref>. Norm-optimal ILC uses a constant linear-time-varying process *Corresponding author. Email: cbmartin129@uexas.edu 2 MECC25_39, Martin</p><p>model to improve the controller after each task iteration with respect to a quadratic cost function. If there are no input constraints, this optimal update takes the form of a constant learning gain <ref type="bibr">[17]</ref>. Otherwise, a quadratic program (QP) must be solved after each iteration <ref type="bibr">[18]</ref>. When a constant linear approximation of the nonlinear model is insufficient, Sequential Quadratic Programming ILC (SQP ILC) becomes necessary. SQP ILC formulates the ILC problem within the SQP framework, which is a robust method for nonlinear constrained optimization offering reliable performance and superlinear convergence when initialized near a feasible solution <ref type="bibr">[19]</ref>. Each iteration of SQP ILC involves performing the task, linearizing the nonlinear process model, and then using these measurements and linear dynamics to solve a QP that improves the control <ref type="bibr">[19]</ref>.</p><p>Given the nonlinear dynamics and input constraints in R2R mechanical dry transfer, SQP ILC is a suitable approach. We thus propose an SQP ILC approach to regulate R2R mechanical dry transfer, incorporating nonlinear system dynamics and input constraints while maintaining low online computational cost. The key contributions of this study are as follows: 1) SQP ILC is used to control an R2R system for the first time. 2) R2R mechanical dry transfer is formulated within the ILC paradigm for the first time. 3) The practicality of this approach is validated experimentally. 4) A case study suggests that the proposed method could significantly improve the R2R transfer quality of advanced 2D materials.</p><p>The paper is structured as follows: Section 2 introduces the dynamic model of R2R mechanical dry transfer, followed by the control problem formulation in Section 3. Section 4 presents the SQP ILC method, while Section 5 validates this approach through an experiment. Finally, Section 6 demonstrates the utility of the approach on an advanced manufacturing example.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="2.">R2R MECHANICAL TRANSFER DYNAMIC MODEL</head><p>This section provides an overview of the R2R mechanical dry transfer dynamic model presented in <ref type="bibr">[20]</ref>, <ref type="bibr">[21]</ref>. The variables and parameters of the system are listed in Table <ref type="table">1</ref>, and the process is illustrated in Figure <ref type="figure">3</ref>.</p><p>The difference equations for the web velocities and unstretched lengths are given as follows:</p><p>&#119905; + 1) = &#119907; &#119894; (&#119905;) + -&#119877; &#119894; 2 &#119879; &#119894; (&#119905;) -&#119891; &#119894; &#119907; &#119894; (&#119905;) + &#119877; &#119894; &#119906; &#119894; (&#119905;) &#119869; &#119894; &#8226; &#119879; &#119878; (1) &#119897; 1 (&#119905; + 1) = &#119897; 1 (&#119905;) + &#119907; 1 -&#119907; &#119901; (&#119905;) 1 + &#120576; 1 (&#119905;) &#8226; &#119879; &#119878; (2. &#119860;) &#119897; &#119894; (&#119905; + 1) = &#119897; &#119894; (&#119905;) + ( &#119907; &#119901; (&#119905;) 1 + &#120576; 1 (&#119905;) -&#119907; &#119894; (&#119905;) 1 + &#120576; &#119894; (&#119905;) ) &#8226; &#119879; &#119878; (2. &#119861;) where &#119905; is the discrete timestep, &#119879; &#119904; is the sampling time, and &#119894; = 2, 3. Equation (2) can be interpreted as a mass conservation, where &#120576; &#119894; = &#119879; &#119894; (&#119905;) &#119860; &#119894; &#119864; &#119894;</p><p>and &#119907; &#119901; is the peeling front velocity. &#119907; &#119901; is a nonlinear function of the web tensions, web speeds, and energy balance at the peeling front <ref type="bibr">[20]</ref>. This energy balance can be summarized as follows <ref type="bibr">[21]</ref>:</p><p>where the time dependence in Eq. ( <ref type="formula">3</ref>) has been omitted for brevity. Thus, there is a nonlinear relationship between the three web tensions and the bending and adhesion energy at the peeling front, parameterized by &#119880;(&#119905;).</p><p>Using the geometry in Figure <ref type="figure">3</ref>, the three web tensions, &#119879; &#119894; (&#119894; = 1, 2, 3), can be uniquely determined from the corresponding unstretched web lengths, &#119897; &#119894; , and vice-versa. This relationship is computed through an iterative root-finding algorithm, represented as <ref type="bibr">[20]</ref>: Initial Parameters Perform Task Control Improvement Measurements Updated Control Parameters = + 1 = FIGURE 3: SIMPLIFIED CARTOON OF R2R MECHANICAL DRY TRANSFER (LEFT); ZOOMED PEELING FRONT (RIGHT). </p><p>Moment of inertia roller &#119894; (kg m 2 ) &#119891; &#119894; , &#119894; = 2,3</p><p>Rotational friction for roller &#119894; (m/s) &#119887; Width of the contact surface (m)</p><p>*Corresponding author. Email: cbmartin129@uexas.edu 3 MECC25_39, Martin</p><p>Thus, using Eqs. ( <ref type="formula">1</ref>)-( <ref type="formula">5</ref>), it is possible to form a state space model of the R2R mechanical dry transfer system: &#119961;(&#119905; + 1) = &#119943;(&#119961;(&#119905;), &#119958;(&#119905;), &#119960;(&#119905;)) <ref type="bibr">(6)</ref> where &#119961;(&#119905;) &#8712; &#8477; &#119899; &#119909; = [&#119907; 2 (&#119905;), &#119907; 3 (&#119905;), &#119879; 1 (&#119905;), &#119879; 2 (&#119905;), &#119879; 3 (&#119905;)] &#119879; , &#119958;(&#119905;) &#8712; &#8477; &#119899; &#119906; = [&#119906; 2 (&#119905;), &#119906; 3 (&#119905;)] &#119879; , and &#119960;(&#119905;) &#8712; &#8477; &#119899; &#119908; = &#119880;(&#119905;). This model has been experimentally validated <ref type="bibr">[20]</ref>, <ref type="bibr">[21]</ref>, and, with small adjustments, is fast enough for online control <ref type="bibr">[11]</ref>.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="3.">CONTROL PROBLEM FORMULATION</head><p>R2R mechanical dry transfer can be modeled as an iterative discrete-time system where each iteration, corresponding to the transfer of one sample from the donor to the receiver substrate, consists of &#119873; timesteps with a step size &#119879; &#119904; . The values of &#119873; and &#119879; &#119904; are determined by the pre-set unwinding velocity and control system sampling rates. The system dynamics follow:</p><p>&#119961; &#119896; (&#119905; + 1) = &#119943;(&#119961; &#119896; (&#119905;), &#119958; &#119896; (&#119905;), &#119905;), &#119905; = 1, 2, &#8230; , &#119873; -1 <ref type="bibr">(7)</ref> where &#119961; &#119896; (&#119905;) &#8712; &#8477; &#119899; &#119909; and &#119958; &#119896; (&#119905;) &#8712; &#8477; &#119899; &#119906; are the state and control at timestep &#119905; in iteration . The initial condition &#119961; &#119896; (1) = &#119961; 0 is assumed constant, as the R2R line typically reaches steady-state operation between transfers. Note that Eq. ( <ref type="formula">7</ref>) is a reformulation of Eq. ( <ref type="formula">6</ref>) within an iterative framework. The time-dependence in Eq. ( <ref type="formula">7</ref>) accounts for &#119960;(&#119905;) &#8712; &#8477; &#119899; &#119908; , which is assumed to be known a priori <ref type="bibr">[22]</ref>.</p><p>Next, suppose the control law has the following form:</p><p>where &#119958; &#773; &#119896; (&#119905;) is the feedforward control signal, &#119922; is a stabilizing feedback gain, &#119961; &#770;&#119896;(&#119905;) is the measured state, and &#119961; &#119903; (&#119905;) is the desired or reference state. The control input is constrained by motor torque limits: &#119940; &#119958; &#8804; &#119958; &#119896; (&#119905;) &#8804; &#119940; &#119958; , &#119940; &#119941;&#119958; &#8804; &#119889;&#119958; &#119896; (&#119905;) &#8804; &#119940; &#119941;&#119958; <ref type="bibr">(9)</ref> where &#119889;&#119963;(&#119905;) = &#119963;(&#119905;) -&#119963;(&#119905; -1).</p><p>Let the process of obtaining measurements from one transfer iteration using the control law (8) be represented as:</p><p>where &#119917;: &#8477; &#119899; &#119906; &#215;(&#119873;-1) &#8594; &#8477; &#119899; &#119909; &#215;&#119873; maps the feedforward control sequence &#119932; &#773; &#119896; = [&#119958; &#773; &#119896; <ref type="bibr">(1)</ref>, &#119958; &#773; &#119896; (2), &#8230; , &#119958; &#773; &#119896; (&#119873; -1) ] to the measured system states &#119935; &#770;&#119896; = [&#119961; &#770;&#119896; <ref type="bibr">(1)</ref>, &#119961; &#770;&#119896;(2), &#8230; , &#119961; &#770;&#119896;(&#119873;) ]. The goal of each iteration is to refine the feedforward control &#119932; &#773; &#119896; using past measurements (&#119935; &#770;1, &#119935; &#770;2, &#8230; , &#119935; &#770;&#119896;-1 ) such that the input constraints (9) are met, and the following cost function is minimized:</p><p>+&#8214;&#119889;&#119958; &#773; &#119896; (&#119905;)&#8214; &#119930; + &#8214;&#119961; &#119896; (&#119873;) -&#119961; &#119903; (&#119873;)&#8214; &#119928; &#119891; <ref type="bibr">(11)</ref> where &#119928;, &#119929;, and &#119930; are weighting matrices that penalize tracking error, control effort, and control variation, respectively; and &#8214;&#119963;&#8214; &#119924; denotes the quadratic form &#119963; &#119879; &#119924;&#119963;. This control problem formulation is summarized in Figure <ref type="figure">4</ref>.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="4.">SQP ILC FOR R2R MECHANICAL DRY TRANSFER</head><p>The key next question is: how to use the system measurements to improve the feedforward trajectory in each iteration according to the objective function <ref type="bibr">(11)</ref>, while ensuring compliance with the dynamic constraints <ref type="bibr">(7)</ref> and the input constraints ( <ref type="formula">9</ref>)? This endeavor presents a challenging nonconvex optimization problem where measurements from new task iterations need to be integrated into the control updates.</p><p>These considerations motivate the use of the SQP ILC framework, first introduced in <ref type="bibr">[19]</ref>. Assuming a quadratic objective, SQP iteratively repeats four steps: 1) execute and measure the controlled task, 2) linearize the dynamic and input constraints around a "guess" of the optimal trajectory, 3) solve a QP to find the optimal update direction, and 4) update the control policy using an appropriate step size. The remainder of this section presents a formulation of this approach crafted to fit R2R mechanical dry transfer.</p><p>After obtaining measurements from the th transfer task (Step 1), the linearized dynamics are updated as follows <ref type="bibr">(</ref>Step 2): &#119912; &#119905;,&#119896; &#8776; &#120597;&#119943; &#120597;&#119961; (&#119961; &#119896; &#119892; (&#119905;), &#119958; &#119896; &#119892; (&#119905;), &#119905;) (12. &#119860;) &#119913; &#119905;,&#119896; &#8776; &#120597;&#119943; &#120597;&#119958; (&#119961; &#119896; &#119892; (&#119905;), &#119958; &#119896; &#119892; (&#119905;), &#119905;) (12. &#119861;) &#119941; &#119905;,&#119896; = &#119961; &#770;&#119896;(&#119905; + 1) -&#119961; &#119896; &#119892; (&#119905; + 1) (12. &#119862;) where &#119961; &#119896; &#119892; (&#119905;) and &#119958; &#119896; &#119892; (&#119905;) represent the "guess" of the optimal state and control. The following QP is then solved to find the update direction for &#119932; &#773; &#119896; (Step 3): min &#916;&#119935; &#119896; ,&#916;&#119932; &#773; &#119896; &#8721;&#8214;&#119961; &#119896; &#119892; (&#119905;) + &#916;&#119961; &#119896; (&#119905;) -&#119961; &#119903; (&#119905;)&#8214; &#119928; &#119873;-1 &#119905;=1 +&#8214;&#119958; &#773; &#119896; (&#119905;) + &#916;&#119958; &#773; &#119896; (&#119905;)&#8214; &#119929; + &#8214;&#119889;(&#119958; &#773; &#119896; (&#119905;) + &#916;&#119958; &#773; &#119896; (&#119905;))&#8214; &#119930; +&#8214;&#119961; &#119896; &#119892; (&#119873;) + &#916;&#119961; &#119896; (&#119873;) -&#119961; &#119903; (&#119873;)&#8214; &#119928; &#119891; (13. &#119860;) subject to the dynamic constraints: &#916;&#119961; &#119896; (&#119905; + 1) = &#119912; &#119905;,&#119896; &#916;&#119961; &#119896; (&#119905;) + &#119913; &#119905;,&#119896; &#916;&#119958; &#119896; (&#119905;) + &#119941; &#119905;,&#119896; (13. &#119861;) and the input constraints: &#119940; &#119958; &#8804; &#119958; &#119896; &#119892; (&#119905;) + &#916;&#119958; &#773; &#119896; (&#119905;) + &#119922;&#916;&#119961; &#119896; (&#119905;) &#8804; &#119940; &#119958; &#119940; &#119941;&#119958; &#8804; &#119889; (&#119958; &#119896; &#119892; (&#119905;) + &#916;&#119958; &#773; &#119896; (&#119905;) + &#119922;&#916;&#119961; &#119896; (&#119905;)) &#8804; &#119940; &#119941;&#119958; (13. &#119862;) FIGURE 4: SUMMARY OF THE ILC PROBLEM FORMULATION FOR R2R MECHANICAL DRY TRANSFER. min &#119869; &#119935; &#119896; , &#119932; &#773; &#119896; such that: &#119961; &#119905; + 1 = &#119943; &#119961; &#119905; , &#119958; &#119905; , &#119905; , input constraints in Eq. (9) hold. &#119932; &#773; &#119896; &#119935; &#770;&#119896; &#119932; &#773; &#119896; 1 = + 1 Transfer one sample/device: &#119935; &#770;&#119896; = &#119917; &#119932; &#773; &#119896; . For &#119905; = 1, 2, &#8230; , &#119873;: R2R Transfer Control law Eq. (8) Motors &#119960; &#119905; &#119961; &#770;&#119896; &#119905; &#119879; &#119878; &#119958; &#119896; &#119905; &#119958; &#773; &#119896; &#119905; Controller with saturation Discrete-time system Improve feedforward traj.: &#119932; &#773; 0 , = *Corresponding author. Email: cbmartin129@uexas.edu 4 MECC25_39, Martin where &#916;&#119961; &#119896; (&#119905;) and &#916;&#119958; &#773; &#119896; (&#119905;) are the update directions, &#916;&#119935; &#119896; = [&#916;&#119961; &#119896; (1), &#916;&#119961; &#119896; (2), &#8230; , &#916;&#119961; &#119896; (&#119873;)], &#916;&#119961; &#119896; (1) = , &#916;&#119932; &#773; &#119896; = [&#916;&#119958; &#773; &#119896; (1), &#916;&#119958; &#773; &#119896; (2), &#8230; , &#916;&#119958; &#773; &#119896; (&#119873; -1)], &#916;&#119958; &#773; &#119896; ( ) = , and &#119958; &#119896; &#119892; ( ) = &#119958; 0 &#119892; ( ) is a user-defined constant. Finally, the feedforward trajectory and guess optimal trajectory are updated as follows (Step 4): &#119935; &#119896; 1 &#119892; &#8592; &#119935; &#119896; &#119892; + &#120572; &#119896; &#916;&#119935; &#119896; &#119892; (14. &#119860;) &#119932; &#773; &#119896; 1 &#8592; &#119932; &#773; &#119896; + &#120572; &#119896; &#916;&#119932; &#773; &#119896; (14. &#119861;) &#119932; &#119896; 1 &#119892; &#8592; &#119932; &#773; &#119896; 1 + &#119922;(&#119935; &#119896; 1 &#119892; -&#119935; &#119903; ) (14. &#119862;) where &#119935; &#119903; = [&#119961; &#119903; (1), &#119961; &#119903; (2), &#8230; , &#119961; &#119903; (&#119873;)], &#119935; &#119896; &#119892; = [&#119961; &#119896; &#119892; (1), &#119961; &#119896; &#119892; (2), &#8230; , &#119961; &#119896; &#119892; (&#119873;)], &#119932; &#119896; &#119892; = [&#119958; &#119896; &#119892; (1), &#119958; &#119896; &#119892; (2), &#8230; , &#119958; &#119896; &#119892; (&#119873; -1)], and &#120572; &#119896; is the step size, given by: &#120572; &#119896; = &#120578; -&#119888; (15) where &#120578; and &#119888; are user-defined. Typically, &#119888; &#8712; [ ,1] and &#120578; &#8712; ( ,1]. This SQP ILC framework is summarized in Algorithm 1. Algorithm 1: SQP ILC Input: &#119935; 0 &#119892; , &#119932; 0 &#119892; , &#119932; &#773; 0 while termination criteria not met do 1. Measurement: &#119935; &#770;&#119896; &#8592; &#119917;(&#119932; &#773; &#119896; ) 2. Linearize: &#119912; &#119905;,&#119896; , &#119913; &#119905;,&#119896; , &#119941; &#119905;,&#119896; &#8592; Expand around &#119961; &#119896; &#119892; (&#119905;), &#119958; &#119896; &#119892; (&#119905;) 3. Optimal control: &#916;&#119935; &#119896; , &#916;&#119932; &#773; &#119896; &#8592; Solution of QP (13) 4. Update: &#119935; &#119896; 1 &#119892; &#8592; &#119935; &#119896; &#119892; + &#120572; &#119896; &#916;&#119935; &#119896; ; &#119932; &#773; &#119896; 1 &#8592; &#119932; &#773; &#119896; + &#120572; &#119896; &#916;&#119932; &#773; &#119896; ; &#119932; &#119896; 1 &#119892; &#8592; &#119932; &#773; &#119896; 1 + &#119922;(&#119935; &#119896; 1 &#119892; -&#119935; &#119903; ) &#8592; + 1</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="5.">EXPERIMENTAL VALIDATION</head><p>The SQP ILC framework for R2R mechanical dry transfer is validated on an experimental testbed, shown in Figure <ref type="figure">5</ref>. The experiment involves peeling Trimaco masking tape, chosen for its suitable adhesion properties, from a PET transfer web. R2R peeling of an adhesive strip shares key characteristics with R2R mechanical dry transfer, including nonlinear peeling dynamics and motor torque constraints <ref type="bibr">[14]</ref>. Therefore, this experiment serves as a practical test of SQP ILC's feasibility and effectiveness for R2R mechanical dry transfer. The physical parameters of this R2R peeling testbed are given in Table <ref type="table">2</ref>. All system states were measured online using loadcells and encoders, and the motor torques were determined using state feedback, as described in Eq. ( <ref type="formula">8</ref>).</p><p>The testing process followed the sequence outlined in Algorithm 1. Each th iteration began with measuring the system's response to a step increase in reference tension, a standard approach for evaluating control regulation performance. State measurements from the th peeling task were then used to update the control parameters according to Eqs. ( <ref type="formula">12</ref>)- <ref type="bibr">(15)</ref>. The new feedforward control was then applied in the ( + 1) th peeling task. Thus, the experiment aligned with the proposed ILC flow for an industrial R2R line: each peeling task generated new measurements, which were used to improve the control policy before the next iteration.</p><p>The feedback gain &#119922; from Eq. ( <ref type="formula">8</ref>) was chosen to be a linear quadratic regulator (LQR) <ref type="bibr">[14]</ref>, while the initial feedforward trajectory was derived from a linearized system model. The first peeling task iteration, which used this LQR feedback and na&#239;ve feedforward trajectory, served as a benchmark for evaluating subsequent iterations. The &#119928;, &#119929;, and &#119930; matrices in the cost function (13.A) were determined using Bryson's rule <ref type="bibr">[23]</ref>, where &#119928; &#119894;,&#119894; ~&#119909;&#119894; -2 , &#119929; &#119894;,&#119894; ~&#119906;&#119894; -2 , and &#119930; &#119894;,&#119894; ~&#119906;&#119894; -2 . Additional weight was given to the elements of &#119928; associated with &#119879; As can be seen in Figure <ref type="figure">7</ref>, the tension RMSE in the first iteration, using LQR feedback and a na&#239;ve feedforward trajectory, is 2.36 N. After 20 iterations of the peeling task within the SQP ILC framework, this RMSE converges to 0.67 N, a 1.69 N or 71% improvement. Notably, 82% of this improvement comes in the first 5 iterations. These results demonstrate that the    *Corresponding author. Email: cbmartin129@uexas.edu 5 MECC25_39, Martin SQP ILC method is both implementable and useful for R2R mechanical peeling dry transfer.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="6.">CASE STUDY: R2R MECHANICAL DRY TRANSFER OF LARGE-SCALE GRAPHENE</head><p>With the practicality and effectiveness of the proposed method demonstrated experimentally, its applicability to advanced manufacturing is further evaluated. Specifically, SQP ILC is utilized to control the R2R mechanical dry transfer of chemical vapor deposition (CVD)-grown graphene, where the donor is the copper growth substrate <ref type="bibr">[7]</ref>, <ref type="bibr">[8]</ref>. This R2R transfer setup, illustrated in Figure <ref type="figure">8</ref>, consists of regularly spaced graphene-copper samples sandwiched between two PET carrier films. The adhesion and bending properties, parameterized by &#119880;, of the PET-PET and graphene-copper interfaces are also detailed in Figure <ref type="figure">8</ref>. Additional physical parameters, sourced from previous experiments <ref type="bibr">[20]</ref>, <ref type="bibr">[21]</ref>, are listed in Table <ref type="table">3</ref>. To ensure minimal damage to the transferred graphene, the objective is to maintain an optimal web tension setpoint during peeling despite abrupt changes in &#119880; when the peeling interface changes.</p><p>As a comparison controller, the nonlinear MPC scheme from <ref type="bibr">[11]</ref> was also used on this task. While this MPC achieves high performance without requiring iterative learning, it must solve a nonlinear optimization over a prediction horizon at each timestep, resulting in drastically higher online computational cost than SQP ILC, which only computes state feedback. N-m] T for both controllers. For the proposed SQP ILC scheme, &#120578; = 1, &#119888; = 0.5, and &#119873; = 60. Additionally, &#119922; will once again be synthesized using an LQR approach. For the nonlinear MPC, the prediction horizon, &#119873; &#119867; , was 20 timesteps.</p><p>A total of 100 graphene transfer tasks were simulated, corresponding to 100 iterations of the SQP ILC framework. The R2R mechanical dry transfer model, used for these simulations, was developed in <ref type="bibr">[20]</ref>. Figure <ref type="figure">9</ref> gives a representative selection of the web tension and control input trajectories of the SQP ILC method as well as the trajectory of the nonlinear MPC for comparison. Note that the nonlinear MPC is only run for one graphene transfer task, since it does not update from iteration to iteration. Figure <ref type="figure">10</ref> gives the average tension RMSE signals for each of the 100 iterations controlled using SQP ILC.</p><p>As can be seen in Figure <ref type="figure">9</ref>, the tension RMSE in the first iteration, using LQR feedback with a na&#239;ve feedforward trajectory, is 2.24 N. After 100 iterations of the transfer task within the SQP ILC framework, this RMSE converges to 0.411 N, a 1.83 N or 82% improvement. Notably, 78% of this performance improvement comes in the first 10 iterations. This reduction in tension RMSE corresponds to an estimated 1.5 k&#937; decrease in electrical sheet resistance, a key metric of graphene transfer quality <ref type="bibr">[7]</ref>. These results suggest that the SQP ILC method can significantly improve R2R mechanical dry transfer quality for advanced applications like CVD graphene.</p><p>Next, the performance of the SQP ILC scheme on this task is compared to that of the nonlinear MPC. The average tension RMSE of the MPC was 0.415 N. The SQP ILC surpassed this performance after 93 iterations, got within 20% of this performance after 25 iterations, and got within 50% of this performance after 10 iterations. Importantly however, a key advantage of SQP ILC is its minimal computational burden, requiring only a single "offline" optimization step between graphene transfer tasks, as defined in Eqs. ( <ref type="formula">12</ref>)-( <ref type="formula">15</ref>). This computation converged within a few seconds on a Dell G7 7500 (Intel i7-10750 CPU, 16GB RAM), making it feasible for highspeed industrial R2R processes with throughputs of 10-100 devices per minute. In contrast, nonlinear MPC must solve a nonconvex optimization problem at every timestep of each peeling task, making it computationally impractical for most real-world industrial applications. A quantitative summary of this performance comparison is presented in Table <ref type="table">4</ref>. Tape PVA Graphene Copper Tape PET PET-PET Graphene-Copper &#119880; = 6.5 J m -1 Graphene-Copper &#119880; = 10 J m -1 PET-PET &#119880; for each peeling section: Sample Sample +1 FIGURE 8: SELECTED TRAJECTORIES OF THE R2R MECHANICAL DRY TRANSFER OF CVD GRAPHENE SIMULATIONS. The grey shaded region corresponds to the time segment when the graphene sample is being peeled. *Corresponding author. Email: cbmartin129@uexas.edu 6 MECC25_39, Martin</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="7.">CONCLUSION</head><p>This study presents an SQP ILC method for controlling R2R mechanical dry transfer. This method achieves excellent performance after about 5-10 iterations by accounting for the nonlinear dynamics, input constraints, and measurements from previous transfer tasks. Moreover, since the control policy updates are calculated between iterations of the transfer task, the method has minimal online computation. Experimental results show that the proposed method is both practical and effective for a lab-scale R2R transfer task. Over 20 iterations, the tension RMSE was reduced by 71%, with most of the improvement occurring in the first several iterations. In addition, a case study indicates that the method can significantly improve the transfer quality of CVD graphene relative to a baseline optimal controller, and that it can achieve equivalent regulation performance compared to an advanced nonlinear MPC with essentially no online computation cost. </p></div></body>
		</text>
</TEI>
