<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">ICA</journal-id><journal-title-group><journal-title>Intelligent Control and Automation</journal-title></journal-title-group><issn pub-type="epub">2153-0653</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ica.2015.62014</article-id><article-id pub-id-type="publisher-id">ICA-56542</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Computer Science&amp;Communications</subject></subj-group></article-categories><title-group><article-title>
 
 
  An Electrothermal Model Based Adaptive Control of Resistance Spot Welding Process
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>iyad</surname><given-names>Kas</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Manohar</surname><given-names>Das</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Department of Electrical and Computer Engineering, Oakland University, Rochester, USA</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>zrkas@oakland.edu(IK)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>30</day><month>03</month><year>2015</year></pub-date><volume>06</volume><issue>02</issue><fpage>134</fpage><lpage>146</lpage><history><date date-type="received"><day>26</day>	<month>February</month>	<year>2015</year></date><date date-type="rev-recd"><day>accepted</day>	<month>18</month>	<year>May</year>	</date><date date-type="accepted"><day>22</day>	<month>May</month>	<year>2015</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  Resistance Spot Welding (RSW) is a process commonly used for joining a stack of two or three metal sheets at desired spots. The weld is accomplished by holding the metallic workpieces together by applying pressure through the tips of a pair of electrodes and then passing a strong electric current for a short duration. Inconsistent weld and insufficient nugget size are some of the common problems associated with RSW. To overcome these problems, a new adaptive control scheme is proposed in this paper. It is based on an electrothermal dynamical model of the RSW process, and utilizes the principle of adaptive one-step-ahead control. It is basically a tracking controller that adjusts the weld current continuously to make sure that the temperature of the workpieces or the weld nugget tracks a desired reference temperature profile. The proposed control scheme is expected to reduce energy consumption by 5% or more per weld, which can result in significant energy savings for any application requiring a high volume of spot welds. The design steps are discussed in details. Also, results of some simulation studies are presented.
 
</p></abstract><kwd-group><kwd>Resistance Spot Welding</kwd><kwd> Adaptive Control</kwd><kwd> Nugget Formation</kwd><kwd> Energy Saving</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>In resistance spot welding, the welding process begins by applying pressure on a stack of metal sheets, held together between a pair of electrodes. A weld current is then passed through the electrodes, causing resistive heating of the metal workpieces and the formation of a welded joint or nugget, as shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>. The formation of a weld nugget strongly depends on the electrical and thermal properties of the sheet and coating materials [<xref ref-type="bibr" rid="scirp.56542-ref1">1</xref>] . Since the contact resistance near the faying surface is much higher than the resistance of the sheets and electrodes, most of the heating is concentrated near the faying surface, causing melting and formation of a nugget</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Resistance spot welding system</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-7900405x5.png"/></fig><p>there. Depending on the thickness and type of material, welding current ranges from 1,000 to 20,000 amperes or more, while the voltage typically is between 1 and 30 volts [<xref ref-type="bibr" rid="scirp.56542-ref2">2</xref>] .</p><p>A Resistance Spot Welding cycle consists of three main stages as follows:</p><p>Stage 1: Squeeze time, which is the time when electrodes press the welded workpieces together.</p><p>Stage 2: Weld time, which is the time when welding current is applied producing heat at the faying surface of the workpieces and thus creating a weld nugget.</p><p>Stage 3: Hold time, which is the time when electrode force still presses the workpieces together and cools the weld down after the welding current is switched off.</p><p>One of the most common applications of resistance spot welding is in the automobile manufacturing industry, where it is used almost universally to weld the sheet metals to form the car body and parts. A typical automotive vehicle today requires about 4000 - 6000 spot welds per vehicle. Considering a worldwide annual production volume of 80 million automotive vehicles, an energy saving RSW controller can result in significant energy savings and reduce carbon footprint accordingly.</p><p>During the past two decades, a number of studies have been carried out to improve the RSW process, which focuses on monitoring and control of weld parameters to improve weld quality. The RSW control techniques proposed to date include Proportional-Integral (PI) [<xref ref-type="bibr" rid="scirp.56542-ref3">3</xref>] , Proportional-Derivative (PD) [<xref ref-type="bibr" rid="scirp.56542-ref4">4</xref>] , Proportional-Integral- Derivative (PID) [<xref ref-type="bibr" rid="scirp.56542-ref5">5</xref>] , Fuzzy [<xref ref-type="bibr" rid="scirp.56542-ref6">6</xref>] - [<xref ref-type="bibr" rid="scirp.56542-ref8">8</xref>] , Neural Networks (NN) [<xref ref-type="bibr" rid="scirp.56542-ref9">9</xref>] [<xref ref-type="bibr" rid="scirp.56542-ref10">10</xref>] , or a combination of Fuzzy and NN [<xref ref-type="bibr" rid="scirp.56542-ref11">11</xref>] . The main drawback of these techniques is that they do not take into account the thermal dynamics of the RSW process, i.e. they do not utilize dynamical models that govern the heat transfer and nugget formation in the RSW process. Also, these systems don’t take into account any welding process variations, such as variations in coating materials, electrode degradation, and weld force variations.</p><p>In this paper, a novel approach to RSW control is presented. This approach has not been explored by other researchers. We start with a simplified heat balance model of a RSW process proposed in [<xref ref-type="bibr" rid="scirp.56542-ref12">12</xref>] and [<xref ref-type="bibr" rid="scirp.56542-ref13">13</xref>] , and then use it to design a controller. This thermal model of the heat balance is a function of nugget growth and it determines the temperature variation during welding time. This model is used later to design an adaptive-one-step- ahead (AOSA) controller and an adaptive-weighted one-step-ahead (AWOSA) controller that compensate for unknown process variations and track a desired reference temperature profile. Finally, some simulation results that show the performance of the proposed controllers are presented and compared to the performance of a PID controller. Simulation results show that AOSA and AWOSA controllers are capable of tracking a reference temperature profile when the weld parameters are unknown, as well as reduce the energy needed to make a weld by 6%.</p><p>The organization of this paper is as follows. Section 2 presents a simplified electrothermal dynamical model of a RSW nugget formation process. The design of adaptive OSA and WOSA controllers is discussed in Section 3. Section 4 presents the results of some simulation studies, and finally some concluding results are provided in Section 5.</p></sec><sec id="s2"><title>2. Electrothermal Dynamical Model of a RSW Nugget Formation Process</title><p>To start with, we consider a simplified heat balance model of a RSW process, presented in [<xref ref-type="bibr" rid="scirp.56542-ref13">13</xref>] . The simplified dynamical model of a RSW process determines the heat balance in the system as a function of nugget temperature. For a simplified nugget model, shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>, the heat balance can be described by the following equations:</p><p>The total heat generation rate, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x6.png" xlink:type="simple"/></inline-formula>is given by</p><disp-formula id="scirp.56542-formula1584"><label>(1a)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x7.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.56542-formula1585"><label>(1b)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x8.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x9.png" xlink:type="simple"/></inline-formula> denotes the welding current, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x10.png" xlink:type="simple"/></inline-formula> denotes the total resistance consisting of the resistance of work pieces, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x11.png" xlink:type="simple"/></inline-formula>, contact resistance, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x12.png" xlink:type="simple"/></inline-formula>, and electrode resistance,<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x13.png" xlink:type="simple"/></inline-formula>. Since <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x14.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x15.png" xlink:type="simple"/></inline-formula> are very small compared to the total contact resistance<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x16.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x17.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x18.png" xlink:type="simple"/></inline-formula> can be neglected in (1b).</p><p>The total contact resistance can then be described as,</p><disp-formula id="scirp.56542-formula1586"><label>(1c)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x19.png"  xlink:type="simple"/></disp-formula><p>A linear relationship between the resistance and temperature is assumed to model the heat generated as a function of temperature. Thus,</p><disp-formula id="scirp.56542-formula1587"><label>(1d)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x20.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.56542-formula1588"><label>(1e)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x21.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.56542-formula1589"><label>(1f)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x22.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x23.png" xlink:type="simple"/></inline-formula> denotes the resistivity of the material, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x24.png" xlink:type="simple"/></inline-formula>denotes the distance from the melting interface to electrode contact surface, p denotes the penetration, A is the cross sectional area, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x25.png" xlink:type="simple"/></inline-formula>denotes the resistivity at reference temperature<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x25.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x26.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x25.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x27.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x25.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x28.png" xlink:type="simple"/></inline-formula> are the temperature to be controlled and the temperature coefficient respectively.</p><p>Substituting (1f) in (1d) and (1e) we get</p><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> A simplified model of a weld nugget</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-7900405x29.png"/></fig><disp-formula id="scirp.56542-formula1590"><label>(1g)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x30.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.56542-formula1591"><label>(1h)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x31.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.56542-formula1592"><label>(1i)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x32.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.56542-formula1593"><label>(1j)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x33.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.56542-formula1594"><label>(1k)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x34.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.56542-formula1595"><label>(1l)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x35.png"  xlink:type="simple"/></disp-formula><p>Substituting (1g) and (1h) in (1a) we get</p><disp-formula id="scirp.56542-formula1596"><label>(1m)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x36.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.56542-formula1597"><label>(1n)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x37.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.56542-formula1598"><label>(1o)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x38.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.56542-formula1599"><label>(1p)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x39.png"  xlink:type="simple"/></disp-formula><p>The heat of fusion required for nugget formation is given by:</p><disp-formula id="scirp.56542-formula1600"><label>(2a)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x40.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.56542-formula1601"><label>(2b)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x41.png"  xlink:type="simple"/></disp-formula><p>where H denotes the heat of fusion per unit volume, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x42.png" xlink:type="simple"/></inline-formula>denotes the nugget volume, and p, a denote the penetration and nugget radius respectively. Substituting (2b) in (2a) and normalizing over the weld duration, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x43.png" xlink:type="simple"/></inline-formula>, we get the heat of fusion per unit time:</p><disp-formula id="scirp.56542-formula1602"><label>(2c)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x44.png"  xlink:type="simple"/></disp-formula><p>Neglecting the heat loss in the surroundings and the electrodes, the heat required to raise temperature by <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x45.png" xlink:type="simple"/></inline-formula> is given by</p><disp-formula id="scirp.56542-formula1603"><label>(3a)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x46.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x47.png" xlink:type="simple"/></inline-formula> denotes the density, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x48.png" xlink:type="simple"/></inline-formula>denotes the specific heat, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x49.png" xlink:type="simple"/></inline-formula>is the volume, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x50.png" xlink:type="simple"/></inline-formula> is the temperature rise. We rewrite (3a) as:</p><disp-formula id="scirp.56542-formula1604"><label>(3b)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x51.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.56542-formula1605"><label>(3c)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x52.png"  xlink:type="simple"/></disp-formula><p>The total heat loss rate is given by</p><disp-formula id="scirp.56542-formula1606"><label>(4a)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x53.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.56542-formula1607"><graphic  xlink:href="http://html.scirp.org/file/4-7900405x54.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.56542-formula1608"><graphic  xlink:href="http://html.scirp.org/file/4-7900405x55.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.56542-formula1609"><label>(4b)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x56.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.56542-formula1610"><label>(4c)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x57.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.56542-formula1611"><label>(4d)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x58.png"  xlink:type="simple"/></disp-formula><p>In the above equations, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x59.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x59.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x60.png" xlink:type="simple"/></inline-formula> denote the axial and radial loss rates, respectively; <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x59.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x60.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x61.png" xlink:type="simple"/></inline-formula>repre-</p><p>sents thermal conductivity, a is the nugget radius;<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x62.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x63.png" xlink:type="simple"/></inline-formula>, represent the melting temperature and the interface temperature at the work piece respectively; <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x64.png" xlink:type="simple"/></inline-formula>is the distance from the melting interface to the electrodes contact area; <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x65.png" xlink:type="simple"/></inline-formula>represents the final penetration to work piece thickness ratio; L is the sheet thickness; <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x66.png" xlink:type="simple"/></inline-formula>represent the electrode radius and thermal diffusivity of work piece respectively.</p><p>The heat balance equation over time <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x67.png" xlink:type="simple"/></inline-formula> is given by</p><disp-formula id="scirp.56542-formula1612"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x68.png"  xlink:type="simple"/></disp-formula><p>Substituting (1n), (2c), (3b), and (4b) in (5) and rearranging it, we get</p><disp-formula id="scirp.56542-formula1613"><label>(6a)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x69.png"  xlink:type="simple"/></disp-formula><p>or, equivalently,</p><disp-formula id="scirp.56542-formula1614"><label>(6b)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x70.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.56542-formula1615"><label>(6c)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x71.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.56542-formula1616"><label>(6b)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x72.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.56542-formula1617"><label>(6c)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x73.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.56542-formula1618"><label>(6d)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x74.png"  xlink:type="simple"/></disp-formula><p>For the sake of notational convenience, let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x75.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x75.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x76.png" xlink:type="simple"/></inline-formula>. Then (6b) can rewritten as</p><disp-formula id="scirp.56542-formula1619"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x77.png"  xlink:type="simple"/></disp-formula><p>Equation (7) represents a bilinear electrothermal dynamical model of a RSW process. Note that this simplified model neglects the heat required to raise the temperature of the electrodes and the nugget surroundings. Also, it assumes that most of the heating occurs near the faying surface due to its high contact resistance. The size of the workpieces is assumed to be infinite in the radial direction and the nugget shape is assumed to be a disk growing radially and axially in the same proportions. The nominal nugget diameter is assumed to be<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x78.png" xlink:type="simple"/></inline-formula>, where L is the sheet thickness.</p><p>Using a first order Euler approximation for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x79.png" xlink:type="simple"/></inline-formula> with a sampling period<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x79.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x80.png" xlink:type="simple"/></inline-formula>, the following discrete time equation is derived from the system Equation (7):</p><disp-formula id="scirp.56542-formula1620"><label>(8a)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x81.png"  xlink:type="simple"/></disp-formula><p>or</p><disp-formula id="scirp.56542-formula1621"><label>(8b)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x82.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.56542-formula1622"><label>(8c)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x83.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.56542-formula1623"><label>(8d)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x84.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.56542-formula1624"><label>(8e)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x85.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.56542-formula1625"><label>(8f)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x86.png"  xlink:type="simple"/></disp-formula><p>Also, k denotes the discrete time index <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x87.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x87.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x88.png" xlink:type="simple"/></inline-formula> denote the sampling instances. The above electrothermal model is characterized by four unknown parameters, namely, A, B, C, and D.</p></sec><sec id="s3"><title>3. Design of a RSW Controller</title><p>To develop a control scheme for controlling the nugget temperature of the RSW model presented by Equation (8a), we realize that it presents a bilinear system characterized by some unknown parameters. These parameters can vary from weld to weld, and in most cases we have no prior knowledge of the parameter values. In view of this, we propose to use an adaptive OSA and WOSA controllers.</p><p>The proposed adaptive control scheme involves measurement of the inputs and outputs of the system, estimation of unknown system parameters using a recursive least squares (RLS) parameter estimation algorithm, and computation of a control signal based on the estimated parameter values. Also, the temperature of the weld nugget is monitored indirectly by assuming it to be proportional to the contact resistance.</p><sec id="s3_1"><title>3.1. Adaptive OSA and WOSA Controllers</title><p>In an adaptive controller, the sampled measurements, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x89.png" xlink:type="simple"/></inline-formula>and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x89.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x90.png" xlink:type="simple"/></inline-formula>, are used to estimate the model parameters, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x89.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x90.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x91.png" xlink:type="simple"/></inline-formula>and D in Equation (8b), using a recursive parameter estimation method, such as recursive least square (RLS). The estimated values of these parameters are then used to compute the OSA/WOSA control signals.</p></sec><sec id="s3_2"><title>3.2. Parameter Estimation</title><p>First we write model Equation (7) in the following form:</p><disp-formula id="scirp.56542-formula1626"><label>(9a)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x92.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.56542-formula1627"><label>(9b)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x93.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.56542-formula1628"><label>(9c)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x94.png"  xlink:type="simple"/></disp-formula><p>Next, the estimated value of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x95.png" xlink:type="simple"/></inline-formula> is computed recursively using the following RLS algorithm:</p><disp-formula id="scirp.56542-formula1629"><label>(10a)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x96.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.56542-formula1630"><label>(10b)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x97.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.56542-formula1631"><label>(10c)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x98.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.56542-formula1632"><label>(10d)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x99.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x100.png" xlink:type="simple"/></inline-formula> is a small number and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x100.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x101.png" xlink:type="simple"/></inline-formula> is chosen to be large. Also, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x100.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x101.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x102.png" xlink:type="simple"/></inline-formula>is always constrained to be non-negative, i.e.,</p><disp-formula id="scirp.56542-formula1633"><label>(10e)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x103.png"  xlink:type="simple"/></disp-formula><p>Given an estimate <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x104.png" xlink:type="simple"/></inline-formula> of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x104.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x105.png" xlink:type="simple"/></inline-formula>, we define the predicted output at time <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x104.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x105.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x106.png" xlink:type="simple"/></inline-formula> as:</p><disp-formula id="scirp.56542-formula1634"><label>(11)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x107.png"  xlink:type="simple"/></disp-formula></sec><sec id="s3_3"><title>3.3. Adaptive-One-Step-Ahead Tracking Controller</title><p>One-step-ahead (OSA) control scheme for linear systems has been well investigated in [<xref ref-type="bibr" rid="scirp.56542-ref14">14</xref>] . An OSA controller attempts to bring the predicted output, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x108.png" xlink:type="simple"/></inline-formula>at time<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x108.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x109.png" xlink:type="simple"/></inline-formula>, to the desired value, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x108.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x109.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x110.png" xlink:type="simple"/></inline-formula>in one step. Thus, it minimizes the following cost function:</p><disp-formula id="scirp.56542-formula1635"><label>(12)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x111.png"  xlink:type="simple"/></disp-formula><p>The corresponding OSA control law is given by [<xref ref-type="bibr" rid="scirp.56542-ref14">14</xref>] :</p><disp-formula id="scirp.56542-formula1636"><label>(13)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x112.png"  xlink:type="simple"/></disp-formula><p>The above control signal needs to be constrained by the maximum current delivery capacity of the controller, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x113.png" xlink:type="simple"/></inline-formula>, as follows:</p><disp-formula id="scirp.56542-formula1637"><label>(14)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x114.png"  xlink:type="simple"/></disp-formula><p>The adaptive OSA controller uses the estimate, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x115.png" xlink:type="simple"/></inline-formula>in Equation (11) to compute the control signal, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x115.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x116.png" xlink:type="simple"/></inline-formula>, from the following adaptive version of Equation (13) above:</p><disp-formula id="scirp.56542-formula1638"><label>(15)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x117.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x118.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x118.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x119.png" xlink:type="simple"/></inline-formula> denote the estimated values of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x118.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x119.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x120.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x118.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x119.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x120.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x121.png" xlink:type="simple"/></inline-formula> respectively, at time <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x118.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x119.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x120.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x121.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x122.png" xlink:type="simple"/></inline-formula></p><p>One of the potential drawbacks of OSA controllers is excessive control efforts that often result from attempting to bring <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x123.png" xlink:type="simple"/></inline-formula> to <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x123.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x124.png" xlink:type="simple"/></inline-formula> in one step. To address this potential problem, an AWOSA controller is discussed below.</p></sec><sec id="s3_4"><title>3.4. Adaptive Weighted One-Step-Ahead Controller</title><p>The excessive effort to bring the output <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x125.png" xlink:type="simple"/></inline-formula> to the desired value <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x125.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x126.png" xlink:type="simple"/></inline-formula> in one step using AOSA may</p><p>result in an unfavorable saturation of the input. The adaptive weighted one-step-ahead controller attempts to seek a tradeoff between tracking accuracy and control effort by considering a slight generalization of the cost function (12) to the form (16) given below. Thus, it minimizes the following cost function:</p><disp-formula id="scirp.56542-formula1639"><label>(16)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x127.png"  xlink:type="simple"/></disp-formula><p>where, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x128.png" xlink:type="simple"/></inline-formula>is chosen to provide a desired tradeoff.</p><p>The minimization of the cost function in (16) leads to the weighted one-step-ahead control law [<xref ref-type="bibr" rid="scirp.56542-ref14">14</xref>] :</p><disp-formula id="scirp.56542-formula1640"><label>(17)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x129.png"  xlink:type="simple"/></disp-formula><p>The above control law is also constrained by the maximum current delivery capacity, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x130.png" xlink:type="simple"/></inline-formula>, as shown in Equation (14) above. The choice of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x130.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x131.png" xlink:type="simple"/></inline-formula> provides a desired tradeoff between tracking accuracy and control effort. A small <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x130.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x131.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x132.png" xlink:type="simple"/></inline-formula> results in good tracking but requires high level of control effort. A large<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x130.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x131.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x132.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x133.png" xlink:type="simple"/></inline-formula>, on the other hand, reduces control efforts at the cost of tracking accuracy.</p><p>The adaptive WOSA controller uses the estimate, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x134.png" xlink:type="simple"/></inline-formula>in Equation (11) to compute the control signal, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x134.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x135.png" xlink:type="simple"/></inline-formula>from the following adaptive version of Equation (17) above:</p><disp-formula id="scirp.56542-formula1641"><label>(18)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x136.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x137.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x137.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x138.png" xlink:type="simple"/></inline-formula> denote the estimated values of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x137.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x138.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x139.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x137.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x138.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x139.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x140.png" xlink:type="simple"/></inline-formula> respectively, at time <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x137.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x138.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x139.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x140.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x141.png" xlink:type="simple"/></inline-formula></p></sec></sec><sec id="s4"><title>4. Simulation Results and Discussion</title><p>This section presents the results of a simulation study showing the performance of the system with the proposed AOSA and AWOSA controllers and also compare them with a PID controller. Each controller is designed for tracking a reference temperature profile.</p><p>The reference temperature profile is a good indicator of the weld quality. Therefore, it is desirable to keep the temperature variation close to a desired variation curve, which may be experimentally predetermined for the good welds. A typical reference temperature profile for good weld is shown in <xref ref-type="fig" rid="fig3">Figure 3</xref> below [<xref ref-type="bibr" rid="scirp.56542-ref1">1</xref>] . Basically, such a curve is characterized by a fast rise of temperature to melting point, melting of the workpieces at the faying surface area which causes a slight drop in temperature, followed by a cooling zone that results from removal of weld current. The actual nugget temperature is measured during the weld cycle using the relationship described by Equation (1f). Depending on the tracking error signal, the welding current is adjusted so as to reduce the temperature error.</p><p>For these simulations, we have selected two sheets of mild steel with the same thickness as the materials to be welded. The force variation and electrode wear are considered as unknown process variables that impact the nugget size (diameter and penetration). The Figures below show the performance of the AOSA, AWOSA, and PID controllers due to 20% increase in nugget diameter and 50% increase in indentation from their desired values.</p><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> Desired reference temperature profile</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-7900405x142.png"/></fig><p><xref ref-type="fig" rid="fig4">Figure 4</xref> shows the performance of the AOSA controller using<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x143.png" xlink:type="simple"/></inline-formula>, where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x143.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x144.png" xlink:type="simple"/></inline-formula> denotes the maximum current delivery capacity of the weld controller. We can see that the AOSA controller adapts to the parameter change and force the output temperature profile to follow the desired temperature profile. Also, we can see that the energy required for the weld is lower than that of the PID controller.</p><p><xref ref-type="fig" rid="fig5">Figure 5</xref> and <xref ref-type="fig" rid="fig6">Figure 6</xref> show the performance of AWOSA controller using <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x145.png" xlink:type="simple"/></inline-formula> with <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x145.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x146.png" xlink:type="simple"/></inline-formula> and</p><p>1, respectively. Here we notice that when <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x147.png" xlink:type="simple"/></inline-formula> is high, the output temperature profile does not follow the desired output temperature profile well. However, increasing <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x147.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x148.png" xlink:type="simple"/></inline-formula> results in decreasing the total energy required for the weld.</p><p><xref ref-type="fig" rid="fig7">Figure 7</xref> shows the performance of the PID controller prior to any parameter change using<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x149.png" xlink:type="simple"/></inline-formula>. After multiple trial and error attempts to get satisfactory results, the parameters of the PID controllers are: Proportional (P) = 0.5, Integral (I) = 26.56, Derivative (D) = 0.</p><p>In <xref ref-type="fig" rid="fig8">Figure 8</xref> we see that the PID controller looses track of the reference temperature profile due to weld parameters change. Also, we can see that PID controller requires more energy for the weld comparing to AOSA and AWOSA.</p><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> Performance of AOSA Controller with 20% increase in nugget diameter and 50% increase in indentation;<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x151.png" xlink:type="simple"/></inline-formula></title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-7900405x150.png"/></fig><fig id="fig5"  position="float"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Performance of AWOSA Controller with 20% increase in nugget diameter and 50% increase in indentation;<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x153.png" xlink:type="simple"/></inline-formula></title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-7900405x152.png"/></fig><fig id="fig6"  position="float"><label><xref ref-type="fig" rid="fig6">Figure 6</xref></label><caption><title> Performance of AWOSA Controller with 20% increase in nugget diameter and 50% increase in indentation;<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x155.png" xlink:type="simple"/></inline-formula></title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-7900405x154.png"/></fig><fig id="fig7"  position="float"><label><xref ref-type="fig" rid="fig7">Figure 7</xref></label><caption><title> Performance of PID Controller prior to unknown parameter variations; I<sub>max</sub> =<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x158.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x158.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x157.png" xlink:type="simple"/></inline-formula></title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-7900405x156.png"/></fig><fig id="fig8"  position="float"><label><xref ref-type="fig" rid="fig8">Figure 8</xref></label><caption><title> Performance of PID Controller with 20% increase in nugget diameter and 50% increase in indentation;<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x161.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x161.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x160.png" xlink:type="simple"/></inline-formula></title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-7900405x159.png"/></fig><p>Comparing the simulation results for the three controllers, we can see that AOSA and AWOSA controllers compensate for the parameter variations and track the reference temperature profile quite well. Simulation results in <xref ref-type="fig" rid="fig5">Figure 5</xref> for the AWOSA controller show satisfactory performance and a good tradeoff between tracking error and total energy required for the weld regardless of change in weld parameters. The output temperature profile follows the desired temperature profile reasonably well during the heating stage prior to the melting point. Also, we can see that the total energy required to make a weld using AWOSA is reduced by 6% comparing to the PID controller when<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x162.png" xlink:type="simple"/></inline-formula>. This can result in significant energy savings for applications requiring a high volume of spot welds, such as manufacturing of automotive vehicles.</p></sec><sec id="s5"><title>5. Conclusion</title><p>This paper presents a new approach for designing adaptive OSA and WOSA controllers for resistance spot welding processes by utilizing a simplified electrothermal dynamical model of the process. Simulation results of AOSA and AWOSA performance are compared with those of a PID controller. These results indicate that using the proposed AOSA and AWOSA controllers, the nugget temperature profile is forced to track a desired reference temperature profile in presence of unknown parameter variations. Also, these controllers reduce the energy consumed to perform a spot weld, which can result in significant energy savings for applications requiring a high volume of spot welds, such as manufacturing of automotive vehicles.</p></sec><sec id="s6"><title>Appendix</title>Boundedness of Nugget Temperature<p>Since a RSW is a time limited process (<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x163.png" xlink:type="simple"/></inline-formula>usually), establishing a proof of asymptotic tracking would be meaningless. However, it is important to make sure that the nugget temperature remains bounded during time<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x163.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x164.png" xlink:type="simple"/></inline-formula>. A theoretical upper bound of the nugget temperature rise, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x163.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x164.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x165.png" xlink:type="simple"/></inline-formula>during time, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x163.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x164.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x165.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x166.png" xlink:type="simple"/></inline-formula>can be established as follows.</p><p>Notice the amount of heat absorbed = the amount of heat supplied ? the amount of heat loss</p><p>Suppose</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x167.png" xlink:type="simple"/></inline-formula>= rise in temperature during time, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x167.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x168.png" xlink:type="simple"/></inline-formula></p><p>Thus,</p><disp-formula id="scirp.56542-formula1642"><label>(19)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x169.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x170.png" xlink:type="simple"/></inline-formula> is a constant.</p><disp-formula id="scirp.56542-formula1643"><label>(20a)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x171.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.56542-formula1644"><label>(20b)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x172.png"  xlink:type="simple"/></disp-formula><p>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x173.png" xlink:type="simple"/></inline-formula> denotes the maximum weld current.</p><disp-formula id="scirp.56542-formula1645"><label>(21)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x174.png"  xlink:type="simple"/></disp-formula><p>where, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x175.png" xlink:type="simple"/></inline-formula>is a constant.</p><p>Thus,</p><disp-formula id="scirp.56542-formula1646"><label>(22a)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x176.png"  xlink:type="simple"/></disp-formula><p>or,</p><disp-formula id="scirp.56542-formula1647"><label>(22b)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x177.png"  xlink:type="simple"/></disp-formula><p>or,</p><disp-formula id="scirp.56542-formula1648"><label>(22c)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-7900405x178.png"  xlink:type="simple"/></disp-formula><p>which proves the boundedness of the nugget temperature rise during weld time,<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-7900405x179.png" xlink:type="simple"/></inline-formula></p></sec></body><back><ref-list><title>References</title><ref id="scirp.56542-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Zhang, H. and Senkara, J. (2012) Resistance Welding Fundamentals and Applications. Taylor &amp; Francis Group, Boca Raton.</mixed-citation></ref><ref id="scirp.56542-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Govik, A. (2009) Modeling of the Resistance Spot Welding Process. M.S. Thesis, Institute of Technology, Linkopings University, Linkoping.</mixed-citation></ref><ref id="scirp.56542-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Won, Y.J., Cho, H.S. and Lee, C.W. (1983) A Microprocessor-Based Control System for Resistance Spot Welding Process. Proceedings of ACC, San Francisco, 22-24 June 1983, 734-738.</mixed-citation></ref><ref id="scirp.56542-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Zhou, K. and Cai, L. (2014) A Nonlinear Current Control Method for Resistance Spot Welding. Proceedings of ASME Transactions on Mechatronics, 19, 559-569. http://dx.doi.org/10.1109/TMECH.2013.2251351</mixed-citation></ref><ref id="scirp.56542-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Salem, M. and Brown, L.J. (2011) Improved Consistency of Resistance Spot Welding with Tip Voltage Control. Proceedings of CCECE, Niagara Falls, 8-11 May 2011, 548-551. http://dx.doi.org/10.1109/ccece.2011.6030511</mixed-citation></ref><ref id="scirp.56542-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Chen, X., Araki, K. and Mizuno, T. (1997) Modeling and Fuzzy Control of the Resistance Spot Welding Process. Proceedings of SICE, Tokushima, 29-31 July 1997, 898-994.</mixed-citation></ref><ref id="scirp.56542-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">El-Banna, M., Filev, D. and Chinnam, R.B. (2006) Intelligent Constant Current Control for Resistance Spot Welding. Proceedings of IEEE Conference on Fuzzy Systems, Vancouver, 16-21 July 2006, 1570-1577.http://dx.doi.org/10.1109/fuzzy.2006.1681917</mixed-citation></ref><ref id="scirp.56542-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Chen, X. and Araki, K. (1997) Fuzzy Adaptive Process Control of Resistance Spot Welding with a Current Reference Model. Proceedings of IEEE Conference on Intelligent Processing Systems, Beijing, 28-31 October 1997, 190-194.</mixed-citation></ref><ref id="scirp.56542-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Shriver, J., Peng, H. and Hu, S.J. (1999) Control of Resistance Spot Welding. Proceedings of ACC, San Diego, 2-4 June 1999, 187-191. http://dx.doi.org/10.1109/acc.1999.782766</mixed-citation></ref><ref id="scirp.56542-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Ivezic, N., Allen Jr, J.D. and Zacharia, T. (1999) Neural Network-Based Resistance Spot Welding Control and Quality Prediction. Proceedings of IPMM, Honolulu, 10-15 July 1999, 989-994. http://dx.doi.org/10.1109/ipmm.1999.791516</mixed-citation></ref><ref id="scirp.56542-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Messler Jr, R.W., Jou, M. and Li, C.J. (1995) An Intelligent Control System for Resistance Spot Welding Using a Neural Network and Fuzzy Logic. Proceeding of IAC, Orlando, October 1995, 1757-1763.http://dx.doi.org/10.1109/ias.1995.530518</mixed-citation></ref><ref id="scirp.56542-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Kim, E.W. and Eagar, T.W. (1988) Parametric Analysis of Resistance Spot Welding Lobe Curve. SAE Technical Paper Series, Warrendale.</mixed-citation></ref><ref id="scirp.56542-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Kas, Z. and Das, M. (2014) A Thermal Dynamical Model Based Control of Resistance Spot Welding. Proceedings of IEEE EIT 2014, Milwaukee, 5-7 June 2014, 264-269. http://dx.doi.org/10.1109/eit.2014.6871774</mixed-citation></ref><ref id="scirp.56542-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Goodwin, G.C. and Sin, K.S. (1983) Adaptive Filtering Prediction and Control. Prentice-Hall, Englewood Cliffs.</mixed-citation></ref></ref-list></back></article>