NEURAL COMPENSATION AND MODELLING OF A HOT STRIP ROLLING MILL USING RADIAL BASIS FUNCTION

Size: px
Start display at page:

Download "NEURAL COMPENSATION AND MODELLING OF A HOT STRIP ROLLING MILL USING RADIAL BASIS FUNCTION"

Transcription

1 Latn Amercan Apple Research 4:4-48(0) NEURAL COMPENSAION AND MODELLING OF A HO SRIP ROLLING MILL USING RADIAL BASIS FUNCION F.G.ROSSOMANDO, J. DENI F an A.VIGLIOCCO Insttuto e Automátca (INAU), Faculta e Ingenería. Unversa Naconal e San Juan (UNSJ). Av. San Martín 09 Oeste J5400ARL. San Juan Argentna, frosoma@naut.uns.eu.ar Laboratóro e Controle e Instrumentação, Unversae Feeral e Esprto Santo (UFES) C. Av. Fernano Ferrar s/n CEP Vtóra-ES-Brasl,.ent@ele.ufes.br Gerenca e Lamnacón, ERNIUM SIDERAR- San Ncolas-Pca e Bs. As.- Argentna, avglocco@ternum.com.ar Abstract In ths paper a Neural Compensaton Strategy for a hot rollng mll process s propose. he target of ths work s to bult a RBF-NN compensaton approxmaton for the classcal force fee forwar an spee controller. A strategy base on neural networks s propose here, because they are capable of moellng many nonlnear systems an ther neural control va RBF-NN approxmaton. Smulatons emonstrate that the propose soluton eals wth sturbances an moelng errors n a better way than classc solutons o. he analyss of the RBF-NN approxmaton error on the control errors s nclue, an control system performance s verfe through smulatons. Keywors hot strp mll, thckness evaton, neural compensaton, raal bass functons. I. INRODUCION he effcency of a hot strp mll can be ncrease f the amount of reecte materal s reuce. A strp s consere as reecte materal f t oes not meet the requrements set by the customer an, therefore, has to be sol as lower qualty prouct or has to be re-melte. hs last opton mples a tremenous amount of extra materal hanlng an hgh energy costs One way to mprove the effcency of the hot rollng mll s through a better fnshng mll control. In hot steel strp rollng processes, force fee forwar (FFF) control s necessary to reuce the effects of rap strp thckness varatons ue to temperature varatons an other factors such as support an bearngs eformaton as well as long term or slow varatons of thckness n the sx- stan fnshng mll (Huang et al. 004;Pttner an Marwan 006). he am of thckness control s to regulate the ext thckness from the roll gap an the spee output of the strp. Moern an avance control methos can be employe to ecrease the thckness stanar evaton an ecrease the thckness error by keepng the output spee of the strp constant. In ths approach, a complete moel for the system s evelope, whch s entcal to the process use n the facltes of ERNIUM SIDERAR-Argentna HSM, wth real ata an parameters prove by ths company. hs real-system moel s use as to bass from whch to bul a smulator for valatng new control strateges. 4 Some of the controllers esgne so far are base only on lnear moels of the HSM. For nstance, the controllers esgne n Rgler et al. (996) an Peersen et al. (995) compensate the roll gap ynamcs, whereas n Galvez et al. (003) the neural network base prectve control allows for overcomng the exstng tme elays n system ynamcs. In Mahfouf et al. (005) a neural compensaton for roll gap an roll spee s propose usng MLP-NN, but the results shown are base only on off-lne austment of the NN. he work Kug et al. (00) presents a nonlnear thckness controller that s acheve usng a ynamc feeback lnearzaton technque, an mplemente on HSM wth a hyraulc austment system. Hearns et al. (004) emphasze the mportance of the nteracton between loop length an ext thckness, an tre a performance trae off between mass flow control an gauge control by changng the output weghtngs, but the weghts can be change to ensure stable mll operaton at the expense of ext gauge control. Alvarez et al. (00) escrbe the esgn an mplementaton of a supervsory system for real-tme compensaton of uneven thckness on both ses of a rolle strp, snce t s base on vsual nspecton. An n Bouazza an Abbass (007) an off-lne moel-base controller an a process smulator are escrbe. In Rossomano an Dent (006) a lnear parameterzaton of two stans of a HSM s presente along wth the esgn of a thckness an looper arm controller base on ts entre moel. One avantage of ther controller s that ts parameters are rectly relate to the mll stan parameters. However, f the parameters are not correctly entfe or f the change wth tme because of, e.g. loa varaton, the performance of the controller can be severely affecte. In Lee et al. (007) the propose controller s mplemente on the conventonal control system, mprovng the performance of the system wthout any atonal cost. he mprovement was emonstrate by usng the real ata whch are collecte usng the propose control scheme n HSM, but t only austs the roll gap to reuce thckness varaton n the HSM. In Dng et al. (008), an aaptve threang technque was use to prect thckness an materal harness errors so as to mofy the setup of the remanng stans wthout the metal strp rollng n.

2 Latn Amercan Apple Research 4:4-48(0) In ths work, a new strategy that shows a goo performance uner nonlneartes an sturbances s presente, such as the col zones n the metal strp. he am of ths work s to bult a neural compensaton for thckness an spee varatons n the roll system. herefore, the ynamc controller s esgne on the bass of a nomnal moel an a RBF-NN compensaton controller capable of learnng the ynamc fference between the esre an the actual moels. Snce the RBF-NN compensates only for a esre moel error ynamcs, ts computatonal cost s sgnfcantly reuce as compare wth a whole NN ynamcs controller. An analyss s mae to stuy the effects of the RBF-NN approxmaton error on the control error when the entre control system, the aaptve ynamcs controller workng, s apple to a roll gap an spee control task. he paper s organze as follows: Secton II presents the overvew of the propose system, an shows the mathematcal representaton of the complete screw an spee system moel. he neural controller s scusse, n Secton III, as well as the corresponng error analyss (Secton III.A). Secton IV presents some smulaton results on the performance of the aaptve controller. Fnally, conclusons are gven n Secton V. II. PLAN MODELLING A. Process escrpton he rollng mll process conssts of pressng a strp metal between a set of rotatng rolls that reuces the thckness of the passng hot-metal strp. he mll housng, the back-up an work rolls an the roll gap austment system are the man machne sectons that partcpate n the rollng mll process; a fnshng mll secton s a set of rollng mll stans; Fgure shows a ynamc scheme of a smple rollng mll stan wth four rolls each one, two work rolls an two back up rolls, wth austment roll gap systems (screw). he rolls n contact wth the strp are the work rolls, they are the man responsble for the rollng mll physcal phenomena. Conserng a gven stan, the strp plate from the prevous stans (flat prouct or a col), s ntrouce n the gap between work rolls, whch s smaller than the thckness of the strp, ths gap s etermne by the poston of the screw. he work rolls rag an press the strp nto the gap to reuce ts thckness. hs strp has to leave an enter n the next stan untl the wshe thckness s reache an cole at the output of the last Stan. he fnshng mll process n ernum-serar conssts of 6 rollng mll stans (F5 to F0). he complete fnshng mll process s shown n Fg.. Fgure : Fnshng mll process Fgure : Mll stan ynamc moel In ths work, the mathematcal moel of the last Stan F0 of the rollng mll tran to evaluate the propose control technque s presente. he analytcal moelng use s ame at valatng through smulatons the obtane results. he results are compare wth the measure real values n the rollng mll stans controlle by the FFF technque (force fee forwar; Bryant, 973). o obtan an approxmate rollng mll stan moel, Alexaner (97) an Orowan (944) theores were use snce they are more relable to obtan the rollng mll force, an can be use n a we range of rollng mll contons. he mathematcal moel was ft an calbrate wth real values measure n the process, thus attanng a closer response to the real force value. hs moel s represente n a general form by Eq. (). S h t, h t, t, t, S,, t () an the rollng torque s expresse by: S h t, h t, t, t, S,, t, () where s the output strp tenson for Stan F0, s the nput strp tenson for Stan F0, h s the output strp thckness, h s the nput strp thckness, s the rolls gap for the Stan, s the frcton Coeffcent for the Stan, S s the Yel stress for the Stan an s the emperature of the strp. B. Screw system moelng he mplemente mathematcal moel consers last stans process, In ths case, the thckness austment by RBF-NN control technque s mae, usng the mathematcal moel expresse by the next equatons. he output thckness of each stan s a functon of the gap rolls an the stan stretchng. h t t E S (3) where h s the output strp thckness of the stan, s the Rolls gap of the stan, s s the Rollng mll force of the stan an E(.) s the Elastcty functon of the stan. here s a slght fference between the roll gap (wthout mll force) an the thckness of the rolle prouct, cause by the elastcty effect of the mll structure. he elastcty functon can be expresse by: 4

3 F.G.ROSSOMANDO, J. DENI F, A.VIGLIOCCO S S S E C C C. (4) 3 he elastcty coeffcents of Eq. (4) have a etermne value epenng on the locaton of the stan mll. he able I shows the elastcty coeffcents for each stan. he motor prouces the followng torque, whch s proportonal to the motor current: ms t kmsims t, (5) where k ms : s the motor torque constant, I ms s the c motor current an ms s the electrc motor torque. he electrc counter force s proportonal to the motor velocty ms t ems t kes (6) where θ ms s the angular poston of the screw motor, e ms s the electromagnetc counter-force an k es s the electrcal constant of the motor. he electrc equaton of DC motor s expresse by: Ims t Ls RI s ms t ems t Vs t. (7) he mechancal equatons of the DC motor, accorng to Newton s n Law, are ms t ms t ms t Jms JB B... B... K t t (8) B ms t K t t P B ms where J ms s the mechancal nerta of the motor (constant), J B s the nerta of the shaft an B mr s the vscous frcton of the motor. he necessary torque to move the screw ( p ) s gven by the followng formula: f S me me sgn H P F t me sgn H F t able : Elastcty Coeffcents Stan C [mm/ton ] C [mm/ton] C 3 [mm] F5-5.07e F e F e F8 -.67e F9 -.84e F e Fgure 3: Mll force acton on the screw (9) where F(t) s the resultant of the forces apple n the recton of the screw axs, expresse n tons; me s the 43 able : Moel parameters Symbol Parameter Unt a n r Gear rato p p Screw ptch mm/ra. J B Inerta of the screw referre to the motor shaft 60 kgm J ms Inerta of the rotor referre to 6 kgm the motor shaft L s Electrcal nuctance of the Hy motor R s Electrcal Resstance of the ohms. DC Motor K es Electrc-mechancal constant 4.3 of the motor M Weght of the bearngs an 0 ton + 5 ton of the backup rolls M Weght of the bearngs an ton. + 7ton. of the work rolls K el Elastcty nex of the stan 436 ton/mm. F0 (lnear approx. of E) J mr Inerta of the rotor referre to 33 kgm the motor shaft (Roll system) J gb Inerta of the gear box referre to the motor shaft 3093 kgm (Roll system) J r Inerta of the roll referre to 74.4 kgm the motor shaft L r Electrcal nuctance of the 8e-5 Hy motor (Roll system) R r Electrcal Resstance of the ohms. DC Motor (Roll system) K er Electrc-mechancal constant 3 of the motor (Roll system) n gb Gear rato of roll system Stan F5 5,6 Stan F6 3,7 Stan F7,90 Stan F8,00 Stan F9,00 Stan F0,00 Mean ameter of the splacement screw, gven n meters (m); α s the angle of the screw s helx; μ s the frcton coeffcent between the splacement screw an the nut an H s the lea of the screw s helx. Functon f s (θ ) states whether the force F(t) s requre to open or to close the roll gap (see Fg. ). he relaton between torque an mll force s shown n Fgure 3. he translaton movement of the system that nvolves the vertcal movement of the backup rolls s rule by the followng equaton. t t F t M B ca E h (0) where M s the Mass of the backup rolls, work rolls an bearngs (M an M n able I), B ca s the vscous frcton coeffcent between the Stans an bearngs of the backup rolls; E - (.) s the elastcty nverse functon of the stan. Γ s can be erve from Eq. (3) an replacng n Eq. (0) an the resultant of the forces apple n the recton of the screw axs s: t t F t M B ca. () S he relaton between the angular splacement of the shaft an the lnear vertcal splacement of the backup rolls s gven by the followng expresson:

4 Latn Amercan Apple Research 4:4-48(0) pp t t () nr where p p s the ptch of the screw, n r s the ransmsson rato of the gear box. he system responsble n charge of postonng the rolls n the esre gap s known as AGC (Automatc Gage Control) (Dark et al., 989; an Gunzburg, 989). he AGC control loop s close by the feeback system of the gap, an the controller s a PI type (proportonal ntegral), whose parameters are auste to obtan a fast uner ampe response of the gap (γ). (Gunzburg, 989) shows the control system of the gap (AGC). he gap controller s a PI controller (AGC control) an can be efne by the followng equaton: t Vs t satkis ms u (3) KPs. ms t u t where θ u (t) s the nput poston of the screw system, sat s the saturaton functon, V s (t) s the output voltage, K Is s the ntegral gan an K Ps s the proportonal gan. he total process, n whch the propose control shoul be apple, s consttute by the AGC system an by the loa moel of the rollng process. he loa moel epens on the output thckness an at the same tme the AGC system s affecte by the rollng loa (Γ s ). he parameters to bul the moel were prove by ERNIUM SIDERAR an are shown n able I C. Man rve control system In orer to evelop the roll spee response moel, t s necessary to ve the mll nto several mechancal components couple by a flexble shaft lnk, an then evelop the moton equatons for the nvual components usng Newton s law. In ths paper, a moel of three masses (the motor nerta, the reucton gear an the roll mll) structure was ntally aopte for the stan mll, whose smplfe mechancal confguraton s represente as: mr t mr t mr t Jmr B mr Kmr mr t gb t Jgb gb t Bgb gb t Kmr mr t gb t... ngb ngb (4) K gb... gb t r t ngb K gb J r r t B r r t gb t r t M t ngb ngb ngb where J mr, J gb, J r are the mechancal nerta of the motor, gear box, an the flexble shaft, respectvely; B mr, B gb, B r are the vscous frcton of the motor, gear box an the flexble shaft, respectvely an K mr, K gb : are the elastcty of the flexble shaft lnk of the motor an the gear box, respectvely. he DC electrc motor ynamcs s the same as that represente by Eqs. (5),(6) an (7). Equatons parameters are k mr, k er, L r an R r. he motor spee controller s a PI type (proportonal ntegral Spee control) an s etermne by a smlar Eq.(3) (wth K Ir an K Pr, the ntegral an proportonal gans, respectvely); but the nput spee s ω u (t). D. hckness Austment Dynamc he whole system ynamc of the stan s represente n the state space as shown n Eq. (5) an eq. (6) kes BB KB Ims ms t ms t t Jeq Jeq J eq ms t Bca n S t r 0 K B t fs ms t t ms t M ps M M ms t t t Rs kes sat Ims ms t KIss t KPs ms t u t t Ls Ls L s I ms t ms t u t s t mr t ker Bmr Kmr Imr mr t mr t gb t mr mr mr mr t J J J mr t gb t gb mr gb gb gb t n K B K mr t gb t gb t gb t r t Jgb Jgb J gb r t r t gb t I mr t K gb B n r gb gb t r t r t t0 r t Jr Jr Jr r t Rr k I er mr sat mr t KIrr t KPr mr t u t Lr Lr L r mr t u t An the output vector s efne by: ps h ES to n (6) v r r where o r (t) an o s (t) are auxlary varables to calculate the ntegral acton, J eq s the total nerta of the system (J ms +J B ), v s the spee output of strp, r work roll raus an 0 s the elay tme (one step tme)- he torque of rollng τ(.) an the loa force Γ(.) are elaye one samplng tme ( 0 ) to avo algebrac loops. Expressng Eq. (5) an Eq. (6) n compact form. t t, t y hx t x f x u where: x I u an (5) (7) ms ms ms mr mr gb gb r r Imr u y h u v. A basc problem wth the thckness control s the ffculty to measure the thckness when the plate s beng rolle. Instea, the thckness s estmate an the uncertanty ntrouce here s reflecte n the stablty analyss. Some conseratons on the moel are: ) he frcton coeffcent n the moel, the yel stress tenson an the stresses forwar an backwar from the strp are consere constants to ths analyss. 44

5 F.G.ROSSOMANDO, J. DENI F, A.VIGLIOCCO ) he sturbances of the process cross through the mathematcal moel affect the output value of the plant. ) he man conseraton s that the control s apple when a strp exsts nto the stan. hs means that the prevous moment to the entrance of the strp n the last stan s not consere. v) he sample tme for purpose of control s o =0.04 sec. to esgn the elay tme Assumpton: he rollng loa s relate to the strp's temperature, nput an output thckness, the frcton coeffcent, the yel stress an stresses forwar an backwar from the strp. Any varaton n sa varables s consere a sturbance. hese sturbances are boune an affect the process outputs. he whole process consere n ths work s shown on Fgure 4. III. PROCESS CONROL A. Error ynamcs approxmaton A smple control strategy, whch s smlar to the wellknown moel reference aaptve control (MRAC) scheme, s propose n Fg. 4. hs scheme s aopte n ths stuy because t has the avantage of generatng the esre nput sgnal wthout requrng that the network be trane ntally offlne. Moreover, the outputs of the neuro controller u N whch has fault tolerant ablty through on-lne learnng only epen on the past an present value nputs, outputs an esre nput reference of the system. he selecton of the moel shoul not greatly excee the ynamcs of the system uner control whle the sturbances are consere constant, otherwse large control efforts mght be requre. Assumng all the states of the system are accessble, the state error e s efne as e=x-x. he error ynamcs for the overall system s, e x x fx, u fx, u (8) ey yy hxhx Eq. (8) can be expane n the esre reference usng aylor seres. fx, u e xxoxx... x xx uu fx, u u u Ou u u x x uu hx ey xxoxx... x xx (9) where O. represents hgher orer terms. Substtutng xx fxu u xx by B, fxu, x by A, uu, uu Fgure 4: Rollng mll control process scheme hx x x x by C, an neglectng all the hgher orer terms. he Eq. (9) can by expresse by: e AeBuu (0) ey Ce Now applyng the KYP (Kalman-Yacubovch-Popov) Lemma (Vyasagar, 993) n Eq (0). Lemma : Let Z(s) = C(sI A) - B be a x transfer matrx such that Z (s)+ Z (-s) has normal rank, where A s Hurwtz, (A,B) s stablzable, an (C,A) s observable. hen Z(s) s Strctly Postve Real (SPR) f an only f there exst symmetrc postve efnte matrces P an Q such that: PA A P Q () PB P B C () In orer to evelop an algorthm to aust u so that the system becomes asymptotcally stable, n spte of parameter varatons an moelng errors. Now conserng uff u u ff (3) where u ff s the correcton sgnal (roll gap an spee) from the prevous Stan an u s the set pont of the Stan (constant). Conserng the correcton sgnals an makng, uun u ff (4) hus, the total control sgnal to the stan mll s the sum of the FFF austment, spee correcton an the RBF- NN controller sgnals, from Eq. (4) an usng Eq. (0), the error ynamcs can be approxmate as: e Ae BuN uff Bu (5) B. RBF-NN Parameterzaton Settng the RBF-NN s nputs as ζ where y, y, y, un, u, u N N can be approxmate by a RBF-NN through on-lne learnng, m un w exp c n (6) w ξ, c, n where w (m x ) an ξ (m x ) are optmal parameter vectors of weghts w an raal bass functons ξ, respectvely; c an η are optmal parameter vectors of centers c an whs η, respectvely; an ε n s the approxmaton error. 45

6 Latn Amercan Apple Research 4:4-48(0) However, the optmal parameter vectors are unknown, so t s necessary to estmate the values. Defne an estmatve functon. m uˆ ˆ ˆ ˆ N w exp c n (7) wˆ ξˆ, cˆ, ˆ n where ŵ an ˆξ are estmate parameter vectors of w an ξ, respectvely; an ĉ an ˆ are estmate parameter vectors of c an η, respectvely. Wth the RBF-NN compensaton, usng Eq. (7) an replacng u by ˆ N u N, the control nput vector u from (4) s, m uwˆ exp ˆ ˆ c nuff (8) wˆ ξˆ nu ff Defnng w w wˆ an ξ ξ ξˆ, an replacng Eq. (3) an Eq. (8) n Eq. (5) the error ynamcs may be wrtten as, ˆ ˆ ˆ e Ae B w ξ w ξ wˆ ξ w ξ n Bu ff (9) ˆ B w ξ wˆ ξ represents the learnng error ψ l where an conserng ˆ ˆ ff w ξ u an w ξ nto n. herefore, a RBF-NN compensator s use for on-lne learnng of. he Eq. (9) can be expresse by u ff ˆ e Ae B w ξ wˆ ξ n (30) Usng an approxmaton for the functon ξξ, c, ξ ˆ, c ˆ, ˆ. Ïn orer to eal wth ξ, the aylor s expanson of ξ s taken about c c ˆ an ˆ. ˆ ξ, c, ξ, cˆ, ˆ Ξ c Φηδ, c, (3) where δ enotes the hgh-orer arguments n a aylor s seres expanson, an Ξ an Φ are ervatves of ξ, c, wth respect to an η at c ˆ,ˆ. hey are expresse as : ξ, c, Ξ c c cˆ η ηˆ (3) ξ, c, Φ η c cˆ η ηˆ Equaton (3) can be expresse as ξξ c Φη δ, c, (33) From (33) the hgh-orer term δ s boune by δ, c, ξξ c Φη ξ Ξ c Φη (34) c 3 η where κ, κ, an κ 3 are some boune constants ue to the fact that RBF an ts ervatve are always boune by constants (the proof s omtte here to save space). Substtute (33) nto (30), t can obtan that ˆ e Ae B ˆ w ξ w Ξ c Φηδ n (35) ˆ e Ae Bw ξ Bwˆ Ξ cbwˆ Φη ε where the uncertan ε Bwˆ δ n, s reasonable to assume ε s boune by a constant ε H, an sup ε (36) H t t C. Stablty Analyss an neural parameters austment o erve the stable tunng law, the followng Lyapunov functon s chosen, V epe trw w c cη η (37) where P s an x symmetrc postve efnte matrx, an Θ, Λ are m x m non-negatve efnte matrces. he ervatve of the Lyapunov functon s gven by, V tr epe epe w w c c η η (38) From eq. () b, Q A P P A snce A s Hurwtz, the Lyapunov functon can always be foun an has a unque soluton an replacng Eq. (30) nto Eq. (38), ˆ V eqe ξ wb Pe c ΞwB ˆ Pe ηφwb ˆ Pe... ε BPetr w w c c η η (39) Substtutng Eq. () n Eq. (39) ˆ V eqe ξ we ˆ ˆ y c Ξwey ηφwe y... ε BPetr w w c c η η (40) Notng n eq. (39), tr w w w w (4) ˆ ξ we y ξ w ey Eq. (40) becomes, V eqe ε BPec wey c... η Φwe η w e tr w w y y (4) (43) where w s the -th column of matrx w an e y s the - th row of vector e y, If c, η an w are selecte as, w ˆ ξe y,... (44) c ˆ wey (45) η ˆ Φwey (46) hen Eq. (43) becomes, V eqe ε BPe (47) V can be emonstrate negatve accorng to Eq.(47), V mn Q e e P B H (48) Let B H, t can be shown rectly that V s negatve H when e P a H Q (49) mn 46

7 F.G.ROSSOMANDO, J. DENI F, A.VIGLIOCCO Now conserng w w w, c c c η η η an w 0, c 0, η 0 the tunng rules are, w ˆ ξe y,... (50) c ˆ wey (5) η ˆ Φwey (5) he negatveness of the Lyapunov functon s guarantee, resultng n the overall system to be stable. IV. SIMULAION RESULS In orer to make the control technque evaluaton collecte real ata of the rollng mll of the steel col were use, wth N o 98 6 lamnate n the ERNIUM SIDERAR plant n Argentna. he obectve of the smulaton shown n ths secton s to utlze the propose control strategy to reuce varatons of the parameters an moelng errors as well as ecrease any possble sturbance. I n ths smulaton the number of RBF neurons s equal to fve (m=5). he sturbance sgnals are shown n fgure 7 an fgure 8 shows the smulate output strp thckness wth RBF-NN correcton an the strp thckness wth conventonal control (FFF). nput thckness [mm] tme (sec.) 950 Output Spee (m/s.) Fgure 7: Output spee of the strp for the RBF-NN compensaton Output hckness (mm.) me (sec.) Output Spee Espesor RBF-NN FFF Desre RBF-NN FFF Desre me (sec.) emperature [ o C] tme (sec.) Fgure 8: Output thcknesses for the RBF-NN compensaton an FFF Loa Force real estme Fgure 5: Input Dsturbances u NN (mm.) Gap Compensaton me (sec.) Spee Compensaton Loa Force (tons) tme (sec.) u NN (m/s) me (sec.) Fgure 6: Gap an spee compensaton from RBF-NN Fgure 9: Estmaton of the loa force from the rollng moel an the measure loa force. he stues on control smulaton nvolve parameter varatons, moelng error an sturbance reecton cases (temperature an nput thckness varatons), where the output thckness an output spee are controlle usng the RBF-NN roll gap an spee compensaton as the manpulate varable. he RBF-NN control varaton error s lower than 60 m an, n the case of conventonal control (FFF), t reaches up to 00 m. 47

8 Latn Amercan Apple Research 4:4-48(0) Fgure 8 shows the output spee strp for stan F0, beng the output expecte reference spee equal to 5.45 m/s. he same results for fferent ata collecte n the rollng mll process were obtane an the thckness error n the case of the RBF-NN control s always smaller than the error of the control FFF, whch s the actual control technque use n the plant. Fgure 9 shows the loa varaton of the actual rollng loa measure n the process, compare to the estmate loa by the control moel, n whch the estmate loa by the moel has a smlar behavor to the actual process. he error between them coul be cause by an error n the measurements, or because some of the nvolve varables were estmate. V. CONCLUSIONS he HSM moel presente n ths artcle s calbrate wth real values measure n the rollng mll process an constructe wth parameters prove by the frm ERNIUM SIDERAR. In ths way, a ynamc moel for the last rollng mll stan s obtane wth a behavor whch s also smlar to the real process wth non-lneartes that was consere. Lyapunov's stablty theory s use to erve a stable tunng rule to upate all parameters n the RBF-NN, thus ensurng the local stablty of the overall system. he results show the strp thckness varaton obtane wth the two control technques. he RBF-NN scheme of the controlle process evaluate by computatonal smulaton shows a smaller thckness varaton than the real system force fee forwar FFF that was taken for comparson presents greater thckness sperson than the RBF-NN compensaton. Smulaton stues base on HSM moel emonstrate that the propose scheme can tolerate parameters varatons, an ths scheme has the avantage that t can be use n parallel, for system-performance verfcaton uner real operaton contons, allowng for a future probable complementng of the system's force fee forwar (FFF) n HSM stans. REFERENCES Alexaner, J.M., On the heory of Rollng, Proc. R. Soc. Lonon, A. 36, (97). Alvarez, J.C., A.B. Díez, D. Alvarez, J.A. González an F. Obeso, hck unevenness compensaton n a hot rollng mll havng automatc gage control, IEEE transactons on nustry applcatons, 38, (00). Bouazza, S.-E. an H.A. Abbass, Moel Base Control System for Hot Steel Strp Rollng Mll Stans, Asan Journal of Informaton echnology, 6, (007). Bryant, G.F., Automaton of anem Mlls, he Iron an Steel Insttute Carlton House errace, SWY5BD, Lonon (973). Dark, O., H. Mabush, I. Degawa an H. Nakamura, Progress of AGC utlzaton technques at plate 48 mll of ota works, Nppon Steel echncal Report, 4, (989). Dng, J., X.L. Hu, J.M. Jao, G.F. She an X.H. Lu, Applcaton of Aaptve hreang echnque to Hot Strp Mll, Int. Journal of Iron an Steel Research, 5, 9-3 (008). Gálvez, J.M., L.E. Záratell an H. Helman, A moelbase prectve control scheme for steal rollng mlls usng neural networks, J. Braz. Soc. Mech. Sc. & Eng., 5, ISSN (003) Gnzburg, V.B., Steel Rollng echnology: heory an Practce, Marcel Dekker. New York an Basel (989). Hearns, G., P. Reeve,.S. Blkhu an P. Smth, Multvarable gauge an mass flow control for hot strp mlls, IFAC Automaton n mnng, mneral an metal proc. (MMM'04), France, 3-34 (004). Huang, M., X. Fang, J. Wang an S. Gu, Roll eccentrcty compensaton control for strp rollng mlls base on wavelet packet e-nosng theory, Intellgent Control an Automaton, 4, (004). Kug, A., K. Schlacher an R. Novak, Nonlnear Control n Rollng Mlls: A New Perspectve, IEEE transactons on nustry applcatons, 37, (00). Lee, Y.K., Y.J. Jang an S.W. Km, Aaptve Feeforwar Automatc Gauge Control n Hot Strp Fnshng Mll, ISIJ Int., 47, (007). Mahfouf, M., Y.Y. Yang, M.A. Gama an D.A. Lnkens, Roll Spee an Roll Gap Control wth Neural Networks compensaton, ISIJ Internatonal, 45, (005). Orowan, E.. he Calculaton of Roll Pressure n Hot an Col Rollng, Proc. Inst. of Mechancal Engneers, 50, (944). Peersen, L.M., Moelng an entfcaton of a hot rollng mll, Amercan Control Conference, Seatle, Washngton, (995). Pttner, J. an A.S. Marwan, Control of a contnuous tanem col metal rollng process, Control Engneerng Practce, 6, (007). Rgler, G.W., H.R. Aberl, W.A. Stauer, K. Astlener, an K.H. Wenberg, Improve Rollng Mll automaton by means of avance control technques an Dynamc Smulaton, IEEE transactons on nustry applcatons, 3, (996). Rossomano, F.G. an J. Dent Flho, Moellng an control of a hot rollng mll, Latn Amercan Apple Research, 36, (006). Vyasagar, M., Nonlnear Systems Analyss, Prentce- Hall, N.J. (993). Receve: December 6, 009 Accepte: September 7, 00 Recommene by Subect Etor: José Guvant

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng

More information

Yukawa Potential and the Propagator Term

Yukawa Potential and the Propagator Term PHY304 Partcle Physcs 4 Dr C N Booth Yukawa Potental an the Propagator Term Conser the electrostatc potental about a charge pont partcle Ths s gven by φ = 0, e whch has the soluton φ = Ths escrbes the

More information

High-Order Hamilton s Principle and the Hamilton s Principle of High-Order Lagrangian Function

High-Order Hamilton s Principle and the Hamilton s Principle of High-Order Lagrangian Function Commun. Theor. Phys. Bejng, Chna 49 008 pp. 97 30 c Chnese Physcal Socety Vol. 49, No., February 15, 008 Hgh-Orer Hamlton s Prncple an the Hamlton s Prncple of Hgh-Orer Lagrangan Functon ZHAO Hong-Xa an

More information

New Liu Estimators for the Poisson Regression Model: Method and Application

New Liu Estimators for the Poisson Regression Model: Method and Application New Lu Estmators for the Posson Regresson Moel: Metho an Applcaton By Krstofer Månsson B. M. Golam Kbra, Pär Sölaner an Ghaz Shukur,3 Department of Economcs, Fnance an Statstcs, Jönköpng Unversty Jönköpng,

More information

Neuro-Adaptive Design - I:

Neuro-Adaptive Design - I: Lecture 36 Neuro-Adaptve Desgn - I: A Robustfyng ool for Dynamc Inverson Desgn Dr. Radhakant Padh Asst. Professor Dept. of Aerospace Engneerng Indan Insttute of Scence - Bangalore Motvaton Perfect system

More information

ENGI9496 Lecture Notes Multiport Models in Mechanics

ENGI9496 Lecture Notes Multiport Models in Mechanics ENGI9496 Moellng an Smulaton of Dynamc Systems Mechancs an Mechansms ENGI9496 Lecture Notes Multport Moels n Mechancs (New text Secton 4..3; Secton 9.1 generalzes to 3D moton) Defntons Generalze coornates

More information

SIMPLIFIED MODEL-BASED OPTIMAL CONTROL OF VAV AIR- CONDITIONING SYSTEM

SIMPLIFIED MODEL-BASED OPTIMAL CONTROL OF VAV AIR- CONDITIONING SYSTEM Nnth Internatonal IBPSA Conference Montréal, Canaa August 5-8, 2005 SIMPLIFIED MODEL-BASED OPTIMAL CONTROL OF VAV AIR- CONDITIONING SYSTEM Nabl Nassf, Stanslaw Kajl, an Robert Sabourn École e technologe

More information

TORQUE-SPEED ADAPTIVE OBSERVER AND INERTIA IDENTIFICATION WITHOUT CURRENT TRANSDUCERS FOR CONTROL OF ELECTRIC DRIVES

TORQUE-SPEED ADAPTIVE OBSERVER AND INERTIA IDENTIFICATION WITHOUT CURRENT TRANSDUCERS FOR CONTROL OF ELECTRIC DRIVES PUBLISHING HOUSE PROCEEDINGS OF THE ROMANIAN ACADEMY, Seres A, OF THE ROMANIAN ACADEMY Volume 4, Number 3/2003, pp.000-000 TORQUE-SPEED ADAPTIVE OBSERVER AND INERTIA IDENTIFICATION WITHOUT CURRENT TRANSDUCERS

More information

( ) = : a torque vector composed of shoulder torque and elbow torque, corresponding to

( ) = : a torque vector composed of shoulder torque and elbow torque, corresponding to Supplementary Materal for Hwan EJ, Donchn O, Smth MA, Shamehr R (3 A Gan-Fel Encon of Lmb Poston an Velocty n the Internal Moel of Arm Dynamcs. PLOS Boloy, :9-. Learnn of ynamcs usn bass elements he nternal

More information

Neuro-Adaptive Design II:

Neuro-Adaptive Design II: Lecture 37 Neuro-Adaptve Desgn II: A Robustfyng Tool for Any Desgn Dr. Radhakant Padh Asst. Professor Dept. of Aerospace Engneerng Indan Insttute of Scence - Bangalore Motvaton Perfect system modelng s

More information

Design of Optimum Controllers for Gas Turbine Engines

Design of Optimum Controllers for Gas Turbine Engines Desgn of Optmum Controllers for Gas Turbne Engnes Junxa Mu 1, Dav Rees 1, Cer Evans 1 an Neophytos Chras 1 School of Electroncs, Unversty of Glamorgan Pontypr, CF37 1DL, Wales, UK Phone: +44() 1443 4859

More information

Responsiveness Improvement of Idling Speed Control for Automotive Using SMC

Responsiveness Improvement of Idling Speed Control for Automotive Using SMC J. Software Engneerng & Applcatons, 009, : 309-315 o:10.436/jsea.009.5040 Publshe Onlne December 009 (http://www.scrp.org/journal/jsea) 309 Responsveness Improvement of Ilng Spee Control for Automotve

More information

Week 9 Chapter 10 Section 1-5

Week 9 Chapter 10 Section 1-5 Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,

More information

Analytical classical dynamics

Analytical classical dynamics Analytcal classcal ynamcs by Youun Hu Insttute of plasma physcs, Chnese Acaemy of Scences Emal: yhu@pp.cas.cn Abstract These notes were ntally wrtten when I rea tzpatrck s book[] an were later revse to

More information

Variable Structure Control ~ Motor Control

Variable Structure Control ~ Motor Control Varable Structure Control ~ Motor Control Harry G. Kwatny Department of Mechancal Engneerng & Mechancs Drexel Unversty Outlne Moels of ac Motors Synchronous motors Inucton motors Brushless c motors VS

More information

Time dependent weight functions for the Trajectory Piecewise-Linear approach?

Time dependent weight functions for the Trajectory Piecewise-Linear approach? Tme epenent weght functons for the Trajectory Pecewse-Lnear approach? Juan Pablo Amorocho an Heke Faßbener Abstract Moel orer reucton (MOR) has become an ubqutous technque n the smulaton of large-scale

More information

An Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors

An Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors An Algorthm to Solve the Inverse Knematcs Problem of a Robotc Manpulator Based on Rotaton Vectors Mohamad Z. Al-az*, Mazn Z. Othman**, and Baker B. Al-Bahr* *AL-Nahran Unversty, Computer Eng. Dep., Baghdad,

More information

Precision Tracking Control of a Piezoelectric-Actuated System

Precision Tracking Control of a Piezoelectric-Actuated System Proceengs of the 5th Meterranean Conference on Control & Automaton, July 7-9, 7, Athens - Greece T3- Precson Trackng Control of a Pezoelectrc-Actuate System Jng-Chung Shen, Wen-Yuh Jywe, *Huan-Keng Chang

More information

On the Coordinated Control of Multiple HVDC Links: Modal Analysis Approach

On the Coordinated Control of Multiple HVDC Links: Modal Analysis Approach R. Erksson an V. Knazkns / GMSARN nternatonal Journal 2 (2008) 15-20 On the Coornate Control of Multple HVDC Lnks: Moal Analyss Approach Robert Erksson an Valerjs Knazkns Abstract There are several possbltes

More information

A Hybrid Variational Iteration Method for Blasius Equation

A Hybrid Variational Iteration Method for Blasius Equation Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method

More information

Finite Element Modelling of truss/cable structures

Finite Element Modelling of truss/cable structures Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures

More information

Competitive Experimentation and Private Information

Competitive Experimentation and Private Information Compettve Expermentaton an Prvate Informaton Guseppe Moscarn an Francesco Squntan Omtte Analyss not Submtte for Publcaton Dervatons for te Gamma-Exponental Moel Dervaton of expecte azar rates. By Bayes

More information

Study on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI

Study on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI 2017 2nd Internatonal Conference on Electrcal and Electroncs: echnques and Applcatons (EEA 2017) ISBN: 978-1-60595-416-5 Study on Actve Mcro-vbraton Isolaton System wth Lnear Motor Actuator Gong-yu PAN,

More information

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

More information

Dynamic Modeling of a Synchronous Generator Using T-S Fuzzy Approach

Dynamic Modeling of a Synchronous Generator Using T-S Fuzzy Approach e-issn : 0975-0 Hee-Jn Lee / Internatonal Journal of Engneerng an echnology (IJE) Dynamc oelng of a Synchronous Generator Usng -S Fuzzy Approach Hee-Jn Lee Department of Electronc Engneerng Kumoh Natonal

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

WHY NOT USE THE ENTROPY METHOD FOR WEIGHT ESTIMATION?

WHY NOT USE THE ENTROPY METHOD FOR WEIGHT ESTIMATION? ISAHP 001, Berne, Swtzerlan, August -4, 001 WHY NOT USE THE ENTROPY METHOD FOR WEIGHT ESTIMATION? Masaak SHINOHARA, Chkako MIYAKE an Kekch Ohsawa Department of Mathematcal Informaton Engneerng College

More information

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1 C/CS/Phy9 Problem Set 3 Solutons Out: Oct, 8 Suppose you have two qubts n some arbtrary entangled state ψ You apply the teleportaton protocol to each of the qubts separately What s the resultng state obtaned

More information

Field and Wave Electromagnetic. Chapter.4

Field and Wave Electromagnetic. Chapter.4 Fel an Wave Electromagnetc Chapter.4 Soluton of electrostatc Problems Posson s s an Laplace s Equatons D = ρ E = E = V D = ε E : Two funamental equatons for electrostatc problem Where, V s scalar electrc

More information

Supporting Information

Supporting Information Supportng Informaton The neural network f n Eq. 1 s gven by: f x l = ReLU W atom x l + b atom, 2 where ReLU s the element-wse rectfed lnear unt, 21.e., ReLUx = max0, x, W atom R d d s the weght matrx to

More information

Adaptive sliding mode reliable excitation control design for power systems

Adaptive sliding mode reliable excitation control design for power systems Acta Technca 6, No. 3B/17, 593 6 c 17 Insttute of Thermomechancs CAS, v.v.. Adaptve sldng mode relable exctaton control desgn for power systems Xuetng Lu 1, 3, Yanchao Yan Abstract. In ths paper, the problem

More information

Parameter Estimation for Dynamic System using Unscented Kalman filter

Parameter Estimation for Dynamic System using Unscented Kalman filter Parameter Estmaton for Dynamc System usng Unscented Kalman flter Jhoon Seung 1,a, Amr Atya F. 2,b, Alexander G.Parlos 3,c, and Klto Chong 1,4,d* 1 Dvson of Electroncs Engneerng, Chonbuk Natonal Unversty,

More information

Introduction. - The Second Lyapunov Method. - The First Lyapunov Method

Introduction. - The Second Lyapunov Method. - The First Lyapunov Method Stablty Analyss A. Khak Sedgh Control Systems Group Faculty of Electrcal and Computer Engneerng K. N. Toos Unversty of Technology February 2009 1 Introducton Stablty s the most promnent characterstc of

More information

Modeling of Dynamic Systems

Modeling of Dynamic Systems Modelng of Dynamc Systems Ref: Control System Engneerng Norman Nse : Chapters & 3 Chapter objectves : Revew the Laplace transform Learn how to fnd a mathematcal model, called a transfer functon Learn how

More information

1.050 Content overview Engineering Mechanics I Content overview. Outline and goals. Lecture 28

1.050 Content overview Engineering Mechanics I Content overview. Outline and goals. Lecture 28 .5 Content overvew.5 Engneerng Mechancs I Lecture 8 Introucton: Energy bouns n lnear elastcty (cont I. Dmensonal analyss. On monsters, mce an mushrooms Lectures -. Smlarty relatons: Important engneerng

More information

Calculation of Coherent Synchrotron Radiation in General Particle Tracer

Calculation of Coherent Synchrotron Radiation in General Particle Tracer Calculaton of Coherent Synchrotron Raaton n General Partcle Tracer T. Myajma, Ivan V. Bazarov KEK-PF, Cornell Unversty 9 July, 008 CSR n GPT D CSR wake calculaton n GPT usng D. Sagan s formula. General

More information

Visualization of 2D Data By Rational Quadratic Functions

Visualization of 2D Data By Rational Quadratic Functions 7659 Englan UK Journal of Informaton an Computng cence Vol. No. 007 pp. 7-6 Vsualzaton of D Data By Ratonal Quaratc Functons Malk Zawwar Hussan + Nausheen Ayub Msbah Irsha Department of Mathematcs Unversty

More information

Large-Scale Data-Dependent Kernel Approximation Appendix

Large-Scale Data-Dependent Kernel Approximation Appendix Large-Scale Data-Depenent Kernel Approxmaton Appenx Ths appenx presents the atonal etal an proofs assocate wth the man paper [1]. 1 Introucton Let k : R p R p R be a postve efnte translaton nvarant functon

More information

On a one-parameter family of Riordan arrays and the weight distribution of MDS codes

On a one-parameter family of Riordan arrays and the weight distribution of MDS codes On a one-parameter famly of Roran arrays an the weght strbuton of MDS coes Paul Barry School of Scence Waterfor Insttute of Technology Irelan pbarry@wte Patrck Ftzpatrck Department of Mathematcs Unversty

More information

EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES

EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES Manuel J. C. Mnhoto Polytechnc Insttute of Bragança, Bragança, Portugal E-mal: mnhoto@pb.pt Paulo A. A. Perera and Jorge

More information

Linear Feature Engineering 11

Linear Feature Engineering 11 Lnear Feature Engneerng 11 2 Least-Squares 2.1 Smple least-squares Consder the followng dataset. We have a bunch of nputs x and correspondng outputs y. The partcular values n ths dataset are x y 0.23 0.19

More information

Chapter 7: Conservation of Energy

Chapter 7: Conservation of Energy Lecture 7: Conservaton o nergy Chapter 7: Conservaton o nergy Introucton I the quantty o a subject oes not change wth tme, t means that the quantty s conserve. The quantty o that subject remans constant

More information

CHAPTER 14 GENERAL PERTURBATION THEORY

CHAPTER 14 GENERAL PERTURBATION THEORY CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves

More information

PHYS 705: Classical Mechanics. Newtonian Mechanics

PHYS 705: Classical Mechanics. Newtonian Mechanics 1 PHYS 705: Classcal Mechancs Newtonan Mechancs Quck Revew of Newtonan Mechancs Basc Descrpton: -An dealzed pont partcle or a system of pont partcles n an nertal reference frame [Rgd bodes (ch. 5 later)]

More information

In this section is given an overview of the common elasticity models.

In this section is given an overview of the common elasticity models. Secton 4.1 4.1 Elastc Solds In ths secton s gven an overvew of the common elastcty models. 4.1.1 The Lnear Elastc Sold The classcal Lnear Elastc model, or Hooean model, has the followng lnear relatonshp

More information

AGC Introduction

AGC Introduction . Introducton AGC 3 The prmary controller response to a load/generaton mbalance results n generaton adjustment so as to mantan load/generaton balance. However, due to droop, t also results n a non-zero

More information

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

Topological Sensitivity Analysis for Three-dimensional Linear Elasticity Problem

Topological Sensitivity Analysis for Three-dimensional Linear Elasticity Problem 6 th Worl Congress on Structural an Multscplnary Optmzaton Ro e Janero, 30 May - 03 June 2005, Brazl Topologcal Senstvty Analyss for Three-mensonal Lnear Elastcty Problem A.A. Novotny 1, R.A. Fejóo 1,

More information

Lecture 21: Numerical methods for pricing American type derivatives

Lecture 21: Numerical methods for pricing American type derivatives Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)

More information

CHAPTER 4 HYDROTHERMAL COORDINATION OF UNITS CONSIDERING PROHIBITED OPERATING ZONES A HYBRID PSO(C)-SA-EP-TPSO APPROACH

CHAPTER 4 HYDROTHERMAL COORDINATION OF UNITS CONSIDERING PROHIBITED OPERATING ZONES A HYBRID PSO(C)-SA-EP-TPSO APPROACH 77 CHAPTER 4 HYDROTHERMAL COORDINATION OF UNITS CONSIDERING PROHIBITED OPERATING ZONES A HYBRID PSO(C)-SA-EP-TPSO APPROACH 4.1 INTRODUCTION HTC consttutes the complete formulaton of the hyrothermal electrc

More information

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,

More information

ENGN 40 Dynamics and Vibrations Homework # 7 Due: Friday, April 15

ENGN 40 Dynamics and Vibrations Homework # 7 Due: Friday, April 15 NGN 40 ynamcs and Vbratons Homework # 7 ue: Frday, Aprl 15 1. Consder a concal hostng drum used n the mnng ndustry to host a mass up/down. A cable of dameter d has the mass connected at one end and s wound/unwound

More information

Inductance Calculation for Conductors of Arbitrary Shape

Inductance Calculation for Conductors of Arbitrary Shape CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors

More information

Non-negative Matrices and Distributed Control

Non-negative Matrices and Distributed Control Non-negatve Matrces an Dstrbute Control Yln Mo July 2, 2015 We moel a network compose of m agents as a graph G = {V, E}. V = {1, 2,..., m} s the set of vertces representng the agents. E V V s the set of

More information

Chapter 9: Statistical Inference and the Relationship between Two Variables

Chapter 9: Statistical Inference and the Relationship between Two Variables Chapter 9: Statstcal Inference and the Relatonshp between Two Varables Key Words The Regresson Model The Sample Regresson Equaton The Pearson Correlaton Coeffcent Learnng Outcomes After studyng ths chapter,

More information

ENTROPIC QUESTIONING

ENTROPIC QUESTIONING ENTROPIC QUESTIONING NACHUM. Introucton Goal. Pck the queston that contrbutes most to fnng a sutable prouct. Iea. Use an nformaton-theoretc measure. Bascs. Entropy (a non-negatve real number) measures

More information

One-sided finite-difference approximations suitable for use with Richardson extrapolation

One-sided finite-difference approximations suitable for use with Richardson extrapolation Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,

More information

Appendix B. The Finite Difference Scheme

Appendix B. The Finite Difference Scheme 140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton

More information

The Study of Teaching-learning-based Optimization Algorithm

The Study of Teaching-learning-based Optimization Algorithm Advanced Scence and Technology Letters Vol. (AST 06), pp.05- http://dx.do.org/0.57/astl.06. The Study of Teachng-learnng-based Optmzaton Algorthm u Sun, Yan fu, Lele Kong, Haolang Q,, Helongang Insttute

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

EEE 241: Linear Systems

EEE 241: Linear Systems EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they

More information

On Liu Estimators for the Logit Regression Model

On Liu Estimators for the Logit Regression Model CESIS Electronc Workng Paper Seres Paper No. 59 On Lu Estmators for the Logt Regresson Moel Krstofer Månsson B. M. Golam Kbra October 011 The Royal Insttute of technology Centre of Excellence for Scence

More information

COMPLEX NUMBERS AND QUADRATIC EQUATIONS

COMPLEX NUMBERS AND QUADRATIC EQUATIONS COMPLEX NUMBERS AND QUADRATIC EQUATIONS INTRODUCTION We know that x 0 for all x R e the square of a real number (whether postve, negatve or ero) s non-negatve Hence the equatons x, x, x + 7 0 etc are not

More information

1 Convex Optimization

1 Convex Optimization Convex Optmzaton We wll consder convex optmzaton problems. Namely, mnmzaton problems where the objectve s convex (we assume no constrants for now). Such problems often arse n machne learnng. For example,

More information

Summary. Introduction

Summary. Introduction Sesmc reflecton stuy n flu-saturate reservor usng asymptotc Bot s theory Yangun (Kevn) Lu* an Gennay Goloshubn Unversty of Houston Dmtry Sln Lawrence Bereley Natonal Laboratory Summary It s well nown that

More information

Spring 2002 Lecture #13

Spring 2002 Lecture #13 44-50 Sprng 00 ecture # Dr. Jaehoon Yu. Rotatonal Energy. Computaton of oments of nerta. Parallel-as Theorem 4. Torque & Angular Acceleraton 5. Work, Power, & Energy of Rotatonal otons Remember the md-term

More information

The Concept of Beamforming

The Concept of Beamforming ELG513 Smart Antennas S.Loyka he Concept of Beamformng Generc representaton of the array output sgnal, 1 where w y N 1 * = 1 = w x = w x (4.1) complex weghts, control the array pattern; y and x - narrowband

More information

Kinematics of Fluid Motion

Kinematics of Fluid Motion Knematcs of Flu Moton R. Shankar Subramanan Department of Chemcal an Bomolecular Engneerng Clarkson Unversty Knematcs s the stuy of moton wthout ealng wth the forces that affect moton. The scusson here

More information

Fall 2012 Analysis of Experimental Measurements B. Eisenstein/rev. S. Errede

Fall 2012 Analysis of Experimental Measurements B. Eisenstein/rev. S. Errede Fall 0 Analyss of Expermental easurements B. Esensten/rev. S. Errede We now reformulate the lnear Least Squares ethod n more general terms, sutable for (eventually extendng to the non-lnear case, and also

More information

Section 8.3 Polar Form of Complex Numbers

Section 8.3 Polar Form of Complex Numbers 80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the

More information

p(z) = 1 a e z/a 1(z 0) yi a i x (1/a) exp y i a i x a i=1 n i=1 (y i a i x) inf 1 (y Ax) inf Ax y (1 ν) y if A (1 ν) = 0 otherwise

p(z) = 1 a e z/a 1(z 0) yi a i x (1/a) exp y i a i x a i=1 n i=1 (y i a i x) inf 1 (y Ax) inf Ax y (1 ν) y if A (1 ν) = 0 otherwise Dustn Lennon Math 582 Convex Optmzaton Problems from Boy, Chapter 7 Problem 7.1 Solve the MLE problem when the nose s exponentally strbute wth ensty p(z = 1 a e z/a 1(z 0 The MLE s gven by the followng:

More information

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD Journal of Appled Mathematcs and Computatonal Mechancs 7, 6(3), 7- www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.3. e-issn 353-588 THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS

More information

Electrical Circuits 2.1 INTRODUCTION CHAPTER

Electrical Circuits 2.1 INTRODUCTION CHAPTER CHAPTE Electrcal Crcuts. INTODUCTION In ths chapter, we brefly revew the three types of basc passve electrcal elements: resstor, nductor and capactor. esstance Elements: Ohm s Law: The voltage drop across

More information

Principles of Food and Bioprocess Engineering (FS 231) Solutions to Example Problems on Heat Transfer

Principles of Food and Bioprocess Engineering (FS 231) Solutions to Example Problems on Heat Transfer Prncples of Food and Boprocess Engneerng (FS 31) Solutons to Example Problems on Heat Transfer 1. We start wth Fourer s law of heat conducton: Q = k A ( T/ x) Rearrangng, we get: Q/A = k ( T/ x) Here,

More information

Physics 111: Mechanics Lecture 11

Physics 111: Mechanics Lecture 11 Physcs 111: Mechancs Lecture 11 Bn Chen NJIT Physcs Department Textbook Chapter 10: Dynamcs of Rotatonal Moton q 10.1 Torque q 10. Torque and Angular Acceleraton for a Rgd Body q 10.3 Rgd-Body Rotaton

More information

Digital Signal Processing

Digital Signal Processing Dgtal Sgnal Processng Dscrete-tme System Analyss Manar Mohasen Offce: F8 Emal: manar.subh@ut.ac.r School of IT Engneerng Revew of Precedent Class Contnuous Sgnal The value of the sgnal s avalable over

More information

ALTERNATIVE METHODS FOR RELIABILITY-BASED ROBUST DESIGN OPTIMIZATION INCLUDING DIMENSION REDUCTION METHOD

ALTERNATIVE METHODS FOR RELIABILITY-BASED ROBUST DESIGN OPTIMIZATION INCLUDING DIMENSION REDUCTION METHOD Proceengs of IDETC/CIE 00 ASME 00 Internatonal Desgn Engneerng Techncal Conferences & Computers an Informaton n Engneerng Conference September 0-, 00, Phlaelpha, Pennsylvana, USA DETC00/DAC-997 ALTERATIVE

More information

Parametric fractional imputation for missing data analysis. Jae Kwang Kim Survey Working Group Seminar March 29, 2010

Parametric fractional imputation for missing data analysis. Jae Kwang Kim Survey Working Group Seminar March 29, 2010 Parametrc fractonal mputaton for mssng data analyss Jae Kwang Km Survey Workng Group Semnar March 29, 2010 1 Outlne Introducton Proposed method Fractonal mputaton Approxmaton Varance estmaton Multple mputaton

More information

Output Feedback Adaptive Robust Control of Uncertain Linear Systems With Disturbances

Output Feedback Adaptive Robust Control of Uncertain Linear Systems With Disturbances Bn Yao 1 School of Mechancal Engneerng, Purdue Unversty, West Lafayette, IN 47907 e-mal: byao@purdue.edu L Xu Research and Advanced Engneerng, Ford Motor Company, Dearborn, MI 48124 Output Feedback Adaptve

More information

Speeding up Computation of Scalar Multiplication in Elliptic Curve Cryptosystem

Speeding up Computation of Scalar Multiplication in Elliptic Curve Cryptosystem H.K. Pathak et. al. / (IJCSE) Internatonal Journal on Computer Scence and Engneerng Speedng up Computaton of Scalar Multplcaton n Ellptc Curve Cryptosystem H. K. Pathak Manju Sangh S.o.S n Computer scence

More information

IMPROVED TRAJECTORY CONTROL FOR AN INTERIOR PERMANENT MAGNET SYNCHRONOUS MOTOR DRIVE WITH EXTENDED OPERATING LIMIT

IMPROVED TRAJECTORY CONTROL FOR AN INTERIOR PERMANENT MAGNET SYNCHRONOUS MOTOR DRIVE WITH EXTENDED OPERATING LIMIT IMPROVED TRAJECTORY CONTRO FOR AN INTERIOR PERMANENT MAGNET SYNCHRONOUS MOTOR DRIVE WITH EXTENDED OPERATING IMIT Abstract M. E. Haue,. Zhong an M. F. Rahman School of Electrcal Engneerng an Telecommuncaton

More information

Part C Dynamics and Statics of Rigid Body. Chapter 5 Rotation of a Rigid Body About a Fixed Axis

Part C Dynamics and Statics of Rigid Body. Chapter 5 Rotation of a Rigid Body About a Fixed Axis Part C Dynamcs and Statcs of Rgd Body Chapter 5 Rotaton of a Rgd Body About a Fxed Axs 5.. Rotatonal Varables 5.. Rotaton wth Constant Angular Acceleraton 5.3. Knetc Energy of Rotaton, Rotatonal Inerta

More information

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

Operating conditions of a mine fan under conditions of variable resistance

Operating conditions of a mine fan under conditions of variable resistance Paper No. 11 ISMS 216 Operatng condtons of a mne fan under condtons of varable resstance Zhang Ynghua a, Chen L a, b, Huang Zhan a, *, Gao Yukun a a State Key Laboratory of Hgh-Effcent Mnng and Safety

More information

RNN Based Auto-tuning of PID Control Gains in Hot Strip Looper Controller

RNN Based Auto-tuning of PID Control Gains in Hot Strip Looper Controller Memors of the Faculty of Engneerng, Okayama Unversty, Vol.40, pp.16-22, January, 2006 RNN Based Auto-tunng of PID Control Gans n Hot Strp Looper Controller Yoshhro Abe Dept. of Electrcal and Electronc

More information

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method Appled Mathematcal Scences, Vol. 7, 0, no. 47, 07-0 HIARI Ltd, www.m-hkar.com Comparson of the Populaton Varance Estmators of -Parameter Exponental Dstrbuton Based on Multple Crtera Decson Makng Method

More information

A Robust Method for Calculating the Correlation Coefficient

A Robust Method for Calculating the Correlation Coefficient A Robust Method for Calculatng the Correlaton Coeffcent E.B. Nven and C. V. Deutsch Relatonshps between prmary and secondary data are frequently quantfed usng the correlaton coeffcent; however, the tradtonal

More information

A particle in a state of uniform motion remain in that state of motion unless acted upon by external force.

A particle in a state of uniform motion remain in that state of motion unless acted upon by external force. The fundamental prncples of classcal mechancs were lad down by Galleo and Newton n the 16th and 17th centures. In 1686, Newton wrote the Prncpa where he gave us three laws of moton, one law of gravty,

More information

Matrix Approximation via Sampling, Subspace Embedding. 1 Solving Linear Systems Using SVD

Matrix Approximation via Sampling, Subspace Embedding. 1 Solving Linear Systems Using SVD Matrx Approxmaton va Samplng, Subspace Embeddng Lecturer: Anup Rao Scrbe: Rashth Sharma, Peng Zhang 0/01/016 1 Solvng Lnear Systems Usng SVD Two applcatons of SVD have been covered so far. Today we loo

More information

829. An adaptive method for inertia force identification in cantilever under moving mass

829. An adaptive method for inertia force identification in cantilever under moving mass 89. An adaptve method for nerta force dentfcaton n cantlever under movng mass Qang Chen 1, Mnzhuo Wang, Hao Yan 3, Haonan Ye 4, Guola Yang 5 1,, 3, 4 Department of Control and System Engneerng, Nanng Unversty,

More information

Indeterminate pin-jointed frames (trusses)

Indeterminate pin-jointed frames (trusses) Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all

More information

Numerical modeling of a non-linear viscous flow in order to determine how parameters in constitutive relations influence the entropy production

Numerical modeling of a non-linear viscous flow in order to determine how parameters in constitutive relations influence the entropy production Technsche Unverstät Berln Fakultät für Verkehrs- un Maschnensysteme, Insttut für Mechank Lehrstuhl für Kontnuumsmechank un Materaltheore, Prof. W.H. Müller Numercal moelng of a non-lnear vscous flow n

More information

Pop-Click Noise Detection Using Inter-Frame Correlation for Improved Portable Auditory Sensing

Pop-Click Noise Detection Using Inter-Frame Correlation for Improved Portable Auditory Sensing Advanced Scence and Technology Letters, pp.164-168 http://dx.do.org/10.14257/astl.2013 Pop-Clc Nose Detecton Usng Inter-Frame Correlaton for Improved Portable Audtory Sensng Dong Yun Lee, Kwang Myung Jeon,

More information

Lab 2e Thermal System Response and Effective Heat Transfer Coefficient

Lab 2e Thermal System Response and Effective Heat Transfer Coefficient 58:080 Expermental Engneerng 1 OBJECTIVE Lab 2e Thermal System Response and Effectve Heat Transfer Coeffcent Warnng: though the experment has educatonal objectves (to learn about bolng heat transfer, etc.),

More information

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

More information

Wavelet chaotic neural networks and their application to continuous function optimization

Wavelet chaotic neural networks and their application to continuous function optimization Vol., No.3, 04-09 (009) do:0.436/ns.009.307 Natural Scence Wavelet chaotc neural networks and ther applcaton to contnuous functon optmzaton Ja-Ha Zhang, Yao-Qun Xu College of Electrcal and Automatc Engneerng,

More information

NMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582

NMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582 NMT EE 589 & UNM ME 48/58 ROBOT ENGINEERING Dr. Stephen Bruder NMT EE 589 & UNM ME 48/58 7. Robot Dynamcs 7.5 The Equatons of Moton Gven that we wsh to fnd the path q(t (n jont space) whch mnmzes the energy

More information

Plasticity of Metals Subjected to Cyclic and Asymmetric Loads: Modeling of Uniaxial and Multiaxial Behavior

Plasticity of Metals Subjected to Cyclic and Asymmetric Loads: Modeling of Uniaxial and Multiaxial Behavior Plastcty of Metals Subjecte to Cyclc an Asymmetrc Loas: Moelng of Unaxal an Multaxal Behavor Dr Kyrakos I. Kourouss Captan, Hellenc Ar Force 1/16 Abstract Strength analyss of materals submtte to cyclc

More information

A Particle Swarm approach for the Design of Variable Structure Stabilizer for a Nonlinear Model of SMIB System

A Particle Swarm approach for the Design of Variable Structure Stabilizer for a Nonlinear Model of SMIB System A Partcle Swarm approach for the Desgn of Varable Structure Stablzer for a Nonlnear Moel of SMIB System NAI A AL-MUSABI**, ZAKARIYA M AL-HAMOUZ*, HUSSAIN N AL-DUWAISH* ** The Petroleum Insttute, Electrcal

More information