Boundary Control of a Tensioned Elastic Axially Moving String

Size: px
Start display at page:

Download "Boundary Control of a Tensioned Elastic Axially Moving String"

Transcription

1 ICCAS5 June -5 KINTEX Gyeonggi-Do Korea Boundary Conrol of a Tensioned Elasi Aially Moving Sring Chang-Won Kim* Keum-Shik Hong** and Hahn Park* *Deparmen of Mehanial and Inelligen Sysems Engineering Pusan Naional Universiy San 3 Jangjeon-dong Gumjeong-gu Busan Korea. Tel.: wonnyday@pusan.a.kr hpark97@pusan.a.kr **Shool of Mehanial Engineering Pusan Naional Universiy San 3 Jangjeon-dong Gumjeong-gu Busan Korea. Tel.: Fa: kshong@pusan.a.kr Absra: In his paper an aive vibraion onrol of a ensioned elasi aially moving sring is invesigaed. The dynamis of he ranslaing sring are desribed by a non-linear parial differenial equaion oupled wih an ordinary differenial equaion. A ime varying onrol in he form of righ boundary ransverse moions is proposed in sabilizing he ransverse vibraions of he ranslaing oninuum. A onrol law based on yapunov s seond mehod is derived. Eponenial sabiliy of he losed-loop sysem is verified. The effeiveness of he proposed onroller is shown hrough simulaions. Keywords: Aially moving sring boundary onrol eponenial sabiliy hyperboli parial differenial equaion yapunov mehod nonlinear sring. INTRODUCTION The onrol problems of aially moving sysems our in various engineering areas: For eample he srips in hin meal-shee produion lines he ables bels and hains in power ransmission lines he magnei apes in reorders he band saws e. The dynamis of hese sysems an be differenly modeled depending on he lengh fleibiliy and onrol objeives of he sysem onsidered. For insane he dynamis of a moving able of an elevaor an be desribed by a sring equaion bu ha of a rubber bel in he radiional mill an be well represened by a bel equaion. The differene beween a sring and a bel lies in wheher he longiudinal elongaion is onsidered or no. In aially moving sysems he ransverse (laeral vibraion of he moving maerial ofen auses a serious problem in ahieving good qualiy. I is also known ha hese vibraions are ofen aused by he eenriiy of a pulley and/or an irregular speed of he driving moor and/or a non-uniform maerial propery and/or environmenal disurbanes. Sine he qualiy requiremen as well as he produiviy in a produion line is geing srier an aive or a semi-aive vibraion onrol is nowadays seriously onsidered. In his paper he vibraion onrol mehod o redue he vibraion whih ours during a oninuous ho-dip zin galvanizing proess is onsidered. In order o ahieve he uniformiy of he zin deposi on he srip surfaes and o redue he zin onsumpion he srip should pass a equidisane from eah of he air knives. Bu due o he shifing and vibraion of he srip a disrepany beween he averaged deposied masses on he lef and righ srip surfaes and a non-uniformiy of he deposied mass aross he srip our. These variaions in deposied mass will degrade he qualiy of he produ. Depending on he hikness of he srip and he disane beween wo suppor poins Aially moving sysem an be modeled as one ou of hree models: a moving beam a moving sring and a moving bel. In he zin galvanizing line he disane beween wo-suppor poins is quie large ompared o he srip hikness and widh is small. Therefore he modeling as a sring is invesed. In he given sysem in-flu and ou-flu mass eis and a he ime varying boundary whih he onrol fore is given by he moving mass work is ourred. So we use Hamilon s priniple for sysem of hanging mass o derive nonlinear equaion of moion. The previous researhes abou aially moving sring didn onsider he elasiiy of he sring and used more han wo sensors. ee and Moe [] derived an opimal boundary fore onrol law ha dissipaes he vibraion energy of an aially moving sring and in [5] analyzed he wave haraerisis of he beam and derived opimal boundary damping laws as a funion of linear veloiy linear slope and linear fore. The onribuions of his paper are he following. Firs boundary onrol law abou non-linear sring o onsider he onsan ension and elasiiy of he sring is derived. Seond he derived boundary onrol law requires only one sensor o apply iself. Addiionally he damping oeffiien of he auaor is designed.. EQUATIONS OF MOTION Fig. shows a shemai (model of he aially moving sring onsidered whih will be used in deriving equaions of moion and a boundary onrol law. The sring is assumed o ravel a a onsan speed. The lef boundary is fied in he sense ha he boundary iself does no have any verial (ransversal movemen bu i allows he maerial o move longiudinally. However he righ boundary permis a ransversal movemen of he sring under a onrol fore and is ime varying. e be he ime be he spaial oordinae along he longiude of moion v be he aial speed of he sring w ( be he ransversal displaemen of he sring a spaial oordinae and ime and be he lengh of he sring. Then he absolue veloiy of he sring a spaial oordinae is given by v = v i Dw( D j = v i { w ( vw ( j ( where D( / D = ( / v ( / and ( = ( ( = ( denoe he parial derivaives in ime and spaial oordinae respeively. The kinei energy of he aially moving sring inluding he auaor and he poenial energy are wrien as follows: T = ρ A{ v ( w d mw ( = EA U Po d ε ε (3 where T is he kinei energy U is he poenial (srain energy ρ is he mass per uni volume (maerial densiy A is he ross-seional area m is he mass of he auaor

2 ICCAS5 Y June -5 KINTEX Gyeonggi-Do Korea ds v w( d Conroller A a P P r P s w w w d P r Fig. An aially moving srip under a righ boundary onrol fore. (ouh roll E is he elasi modulus of he sring ε is he srain due o he ension P. The poenial energy is proporional o he inrease in sring lengh ds when ompared o he sring a res. For small slopes see Fig. he following relaionship for sring elongaion is valid: ds d ( ε = w w( = d. ( Therefore (3 is rewrien as follows EA U = P w w d. (5 Now o derive he equaions of moion Hamilon s priniple for he sysems wih hanging mass is uilized as follows: δ ( T U W.... = n Wr b d (6 where W n.. is he non-onservaive work W r.b. is he virual momenum ranspor a he righ boundary (no variaions a he lef boundary. The variaions of he non-onservaive work and he virual momenum ranspor a he righ boundary are δwn.. = F ( δ w( dw δ w( (7 δ Wr. b. = { w δ w( ( where F ( is he onrol fore and d is he damping oeffiien of he auaor. Now he variaions of ( and (5 respeively are δ T = ρa ( w ( δw vδw mw δw (9 EA 3 δ U = P wδwd wδwd. ( The subsiuion of (7-( ino (6 yields ( δ T δu δw.. n δwr b d = A( w vw w d d ρ δ ( ρ Avw w δw d d EA 3 P w w δ w d d [ mw δ w δw( ]d d Fig. Shemai of sring elongaion for small slopes. { F ( d w δ w( d { Avw w δw( = ρ. ( And he inegraion of ( by pars yields [( ρ Aw w δw] d ( ρ Aw w δw d d [( ρ Avw w δw] d ( w w δwdd EA 3 P w w w d δ 3EA P w w w wdd δ [ mw δw( ] { mw F ( d w δw( d { Avw w ρ δw( d. ( Noe ha δ w( = beause he lef end is fied (i.e. w ( =. Therefore ( is rewrien as follows: ( ρ Aw w P ρ Av w 3EA δwdd wδwdd { mw F ( d w δw( d EA P w w w( d = δ. (3 Sine δ w is arbirary eep for he requiremen ha he lef end is fied (i.e. δ w( = he following governing equaion and a boundary onsrain a he righ end are derived as follows: 3EA ρaw w P ρ Av w w = ( F ( = mw d w

3 ICCAS5 EA P w w (5 where w is he loal aeleraion in he ransversal direion of he sring w is he Coriolis aeleraion and v w is he enripeal aeleraion. Remark: For a linear sysem whih is he ase ha 3EAw / = he soluion of ( an be obained hrough he mehod of separaion of variables. In his ase he naural nπ frequeny is given by ω n = ( v n = 3 where = P ρa is alled he wave veloiy (see []. The naural frequeny dereases as he raveling speed inreases. If he raveling speed is equal o he wave veloiy he naural frequeny beomes zero and a divergene of he soluion ours. In his sense is alled he riial speed v r. Hene he following is also assumed in his paper. < v < vr = P ρa. (6 If using he parameers in Table we an alulae = P / ρ A = m/se. v r 3. BOUNDARY CONTRO AW In his seion we design a boundary onroller inluding a damper o redue he ransversal vibraion of a longiudinally moving elasi sring. A posiive definie funion in he form of oal mehanial energy of he sring eluding he auaor is firs onsidered as follows: Vs ( = ρ A( w EA P w w d (7 where he subsrip s sands for sring. emma : Consider a funional V ~ ~ V ( = Vs ( V ( ( wih he seond (omplemenary erm is defined by V ( = ρ A w ( w (9 where > is a onsan. Then (7 and ( are equivalen ha is here eis onsans > and < C < saisfying ~ ( C VS ( V ( ( C VS (. ( Proof: The eisene of suh and C will be proved. Firs (9 beomes V ( = ρ A w ( w where ρa w d P ρ A w P ( w d d ( w vw d EA w ρ A EA d C V ( ( S June -5 KINTEX Gyeonggi-Do Korea ρa C = min( P ρa EA >. ( Hene he following holds: CV s ( V ( CV s (. (3 By adding V s ( a boh sides of (3 we obain ( C ~ V S ( V ( ( C V S (. ( In order o C > he range of is resried by min( P ρa EA < <. (5 ρa Verifying he range of using he parameers in Table < < / =.5 is obained. In his paper =.3 is seleed and C =. 6 is alulaed. emma is proved. Now wih emma he following yapunov funion andidae V ( whih is basially equivalen o he oal mehanial energy of he sring and he auaor is proposed. ~ V ( = V ( Va ( (6 where he addiional auaed-relaed erm V a ( is defined as m Va ( = { w ( v w. (7 To derive he ime derivaive of (6 a fied onrol volume is inrodued as in Fig. 3. Volume II represens he par of he sring ha oupies he inner par of he onrol volume a an arbirary ime whereas I and III represen he influ and efflu of he sring a d respeively. Using Reynolds ranspor heorem he ime-derivaive of (6 is given by dv ( / d = V / v V /. ( The firs erm in he righ-hand side of ( means he ime rae of he equivalen energy wihin he onrol volume and he seond erm is he ne energy flu ino he onrol volume. Noing ha V inludes V s in (7 V in (9 and V a in (7 individual erms are evaluaed as follows: Vs ( / = ρ A( w ( w EA P w w w d = = v ( [ ] ( P [ ] P w w w EA 3 [ ] [ ] 3EAv ww w [ w ] onrol volume (9 = I II III v d Fig. 3 Conrol volume of an aially moving sring sysem wih a ime varying righ boundary. d

4 ICCAS5 v Vs ( / = v ρ A( w ( w EAv 3 v P wwd wwd Av vp = ρ [( w ] [ ] EAv w [ w ]. (3 Using (9 and (3 dv s ( / d beomes v ( [ ] ( P [ ] s ( / d = P w w w EA 3 [ w ] Av 3EAv w ρ [ w ] [ ] w Av ρ [( w ] vp vw [ w ] EAv [ w ] = P w w vp { w w ( EA 3 EAv w w { w w ( dv. (3 For V ( he following are derived: V ( / = ρ A w ( w ρ A w ( w (3 v V ( / = vρ A w ( w vρ A w ( w vρ A w ( w. (33 Using (3 and (33 dv ( / d beomes dv ( / d = vρ A ( w w w w w w Av ρ wwd ρ A w w d w ( ρaw w w ρ Av w d. (3 To rewrie (3 he following inegraions by pars are uilized. ( w w w w ww = [ w w ] = w( w (35 wwd = w w d (36 ww d = w w d. (37 Also using he equaion of moion ( he following equaion is derived. w ( ρaw w w 3EA = w P w w d. (3 Therefore he subsiuion of (35-(3 ino (3 yields dv ( / d = w w w w d P P w w d June -5 KINTEX Gyeonggi-Do Korea 3 EA 3EA ρa w w d w ρa w d Av w ρ d. (39 Also he ime-derivaive of (7 is dv a / d = m w v w ( { ( { w ( v w. ( Finally he main par in his paper is saed as follows: Theorem: Consider he following aially moving sysem 3EA ρaw w P ρ Av = w w w ( = mw d w EA P w w( = F ( w = w ( w ( = w& (. ( ( ( If he onrol fore F and he damping oeffiien d are given by F ( = Kw ( ρ A / < d < vρa /( v (3 where K = m( v is he onrol gain hen he losed-loop sysem is eponenially sable. Proof: The subsiuion of ( ino (5 yields mw = d w EA P w w m( v w.( The subsiuion of ( ino ( yields dv ( / d = w v w a { ( EA 3 d w P w w. (5 Therefore ombining (3 (39 and (5 he ime-derivaive of he yapunov funion andidae (6 is given by EA dv ( / d = ( P w w ρa ( P wd w d 3EA E Av w d vp ( w w ( ρa d w ρ Av d ( v w w { ( P w EA w ρa ( P wd w d 3EA E Av w d vp ( w w ( A d w P w ρ ( ρa d w ρ Av d ( v w w. (6 { Sine v is under he riial speed (see (6 P > is saisfied. Noe ha he firs and eighh erms afer he equaliy

5 ICCAS5 sign have been spli ino wo halves. Noe also ha he firs eigh erms followed by he inequaliy sign are all negaive whereas he las hree erms are ombined o make hem negaive as a whole. The following noaions are inrodued. φ = ( P ρ Av / (7 > ( d ρ A { ρ Av d ( v φ = / ( ψ =. (9 From ( and (9 if d saisfies he range d ρ A d and > ( v d d < ρ Av (5 hen he las hree erms in (6 saisfies he following inequaliy. φ w ψ w w φ w min{ φ ψ / φ ( w w. (5 Therefore from (6 and (5 he asympoi sabiliy of he losed-loop sysem is assured. Now he eponenial sabiliy is shown wih furher manipulaion of he erms. (6 an be rewrien by spliing he hird erm ino wo pars as EA dv ( / d w vp w ( E Av w ( ( P wd ρa ( P ( vw w d v 3 EA ( w d P ( w ρa d w min{ φ ψ ϕ{ w ( w EA E Av w vp w ( w ( ρ A ( w d P Av w ρ d ρa ( P min w d ( vw v 3 EA ( w d P ( w ρa d w min{ φ ψ ϕ{ w ( w. (5 By using he following inequaliy w d ( vw ( w (53 and eliminaing he firs four and las (negaive erms (5 an be rewrien as ρa ( P dv ( / d min ( w v 3EA ( P wd w d ( P w ρa d w June -5 KINTEX Gyeonggi-Do Korea ( P ( P min 3 P P EA ρa wd wd ( w ( P ρa min d m( v m m { w ( v w. (5 If using (7 and (6-(7 (5 an be epressed as ( P ( P ~ dv ( / d min 3 V ( P ( P ρa min d V ( A m( v m ( P ( P min 3 P ( P ρa d m( v m ( V ~ ( VA( = λv (. (55 (55 means he following relaionship λ Voal ( V e (56 where V = V ( and λ is given by ( P ( P λ = min 3 P ( P ρa d >. (57 m( v m Therefore all he variables inluded in (6 onverge eponenially o zero.. IMPEMENTATION AND SIMUATION Implemening ( and (3 requires wo hings: measuremen of and saisfaion of he range d < d < d w. Sine he damping range is relaed o a design problem of he auaor i mus be answered ahead. Using he parameers in Table wih =.3 (5 is verified as follows: 7.67 < d < 7.7. (5 Therefore he eisene of suh a damping range is assured. The implemenaion of w an be ahieved by bakwards differening of w measured a eah sep. To demonsrae he performane of he losed loop sysem ompuer simulaions using he finie differene sheme have been performed. The plan parameers used for simulaions are gahered in Table. Wih =.3 given in (5 he onrol gain is alulaed as follows: K = m( v = 5(.3 = 39. (59 For simulaion purpose le d = 5. e he iniial ondiions be w ( = sin(3π [m] w ( = [m/s]. (6

6 ICCAS5 June -5 KINTEX Gyeonggi-Do Korea Fig. The ransverse displaemen wih onrol gain K = 39 and damping oeffiien d = 5 w ( where = m. Fig. 5 The ransverse displaemen wih onrol gain K = 39 and damping oeffiien d = 5 w where = m. Fig. and Fig. 5 show he ransverse displaemen a = / and he onrol fore a = respeively. The designed onroller eliminaes he ransversal vibraion in 3 seonds. Fig. 6 shows he variaion of boundary onrol inpu a =. 5. CONCUSIONS This paper invesigaed a ransverse vibraion suppression sheme of an aially moving non-linear sring and disussed he developmen of an effiien aive onroller and is implemenaion. In previous researhes he sring was mosly modeled as a linear sring. However we derived a non-linear sring PDE equaion wih auaor dynamis. The boundary onrol law was derived by yapunov mehod. Eponenial sabiliy of he losed-loop sysem was proved. The effiieny of he designed onroller was shown by simulaions. ACKNOWEDGEMENT This work was suppored by he Minisry of Siene and Tehnology of Korea under a program of he Naional Researh aboraory gran number NR M J--3-. REFERENCES [] Hong K. S. Kim C. W. and Hong K. T. Boundary Conrol of an Aially Moving Bel Sysem in a Thin-Meal Produion ine Inernaional Journal of Conrol Auomaion and Sysems Vol. No. pp Fig. 6 The onrol fore used in Fig. 5. Table The plan parameers for simulaion. Symbols Definiions Values A Cross-seion area.5. 5 m engh of sring m P Tension of he sring 9 kn m Mass of he auaor 5 kg v Sring moving speed m/s ρ Mass per uni area 7 5 kg/m d Damping oeffiien 5 Ns/m [] Kim C. W. Hong K. S. and Park H. Boundary Conrol of an Aially Moving Sring: Auaor Dynamis Inluded Journal of Mehanial Siene and Tehnology Vol. 9 No. pp [3] Khalil H. Nonlinear Sysems Prenie Hall nd ediion 996. [] ee S. Y. and Moe C. D. Vibraion onrol of an aially moving sring by boundary onrol ASME Journal of Dynami Sysems Measuremen and Conrol Vol. No. pp [5] ee S. Y. and Moe C. D. Wave haraerisis and vibraion onrol of ranslaing beams by opimal boundary damping Journal of Vibraion and Aousis Vol. No. pp [6] ee S. Y. and Park S. G. Free Vibraion Charaerisis of a Sring wih Time-Varying engh Korea Soiey for Noise and Vibraion Engineering Vol. 9 No.5 pp [7] i Y. Aron D. and Rahn C. D. Adapive vibraion isolaion for aially moving srings: heory and eperimen Auomaia Vol. 3 No. 3 pp [] Ryu D. H. and Park Y. P. Transverse vibraion onrol of an aially moving sring by veloiy boundary onrol Transaions of KSME A Vol. 5 No. pp [9] Ying S. and Tan C. A. Aive Vibraion Conrol of he Aially Moving Sring Using Spae Feedfoward and Feedbak Conrollers Journal of Vibraion and Aousis Vol. No. pp

Analysis of Tubular Linear Permanent Magnet Motor for Drilling Application

Analysis of Tubular Linear Permanent Magnet Motor for Drilling Application Analysis of Tubular Linear Permanen Magne Moor for Drilling Appliaion Shujun Zhang, Lars Norum, Rober Nilssen Deparmen of Eleri Power Engineering Norwegian Universiy of Siene and Tehnology, Trondheim 7491

More information

(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)

(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4) Phase Plane Analysis of Linear Sysems Adaped from Applied Nonlinear Conrol by Sloine and Li The general form of a linear second-order sysem is a c b d From and b bc d a Differeniaion of and hen subsiuion

More information

Linear Quadratic Regulator (LQR) - State Feedback Design

Linear Quadratic Regulator (LQR) - State Feedback Design Linear Quadrai Regulaor (LQR) - Sae Feedbak Design A sysem is expressed in sae variable form as x = Ax + Bu n m wih x( ) R, u( ) R and he iniial ondiion x() = x A he sabilizaion problem using sae variable

More information

Amit Mehra. Indian School of Business, Hyderabad, INDIA Vijay Mookerjee

Amit Mehra. Indian School of Business, Hyderabad, INDIA Vijay Mookerjee RESEARCH ARTICLE HUMAN CAPITAL DEVELOPMENT FOR PROGRAMMERS USING OPEN SOURCE SOFTWARE Ami Mehra Indian Shool of Business, Hyderabad, INDIA {Ami_Mehra@isb.edu} Vijay Mookerjee Shool of Managemen, Uniersiy

More information

Hybrid probabilistic interval dynamic analysis of vehicle-bridge interaction system with uncertainties

Hybrid probabilistic interval dynamic analysis of vehicle-bridge interaction system with uncertainties 1 APCOM & SCM 11-14 h Deember, 13, Singapore Hybrid probabilisi inerval dynami analysis of vehile-bridge ineraion sysem wih unerainies Nengguang iu 1, * Wei Gao 1, Chongmin Song 1 and Nong Zhang 1 Shool

More information

mywbut.com Lesson 11 Study of DC transients in R-L-C Circuits

mywbut.com Lesson 11 Study of DC transients in R-L-C Circuits mywbu.om esson Sudy of DC ransiens in R--C Ciruis mywbu.om Objeives Be able o wrie differenial equaion for a d iruis onaining wo sorage elemens in presene of a resisane. To develop a horough undersanding

More information

Generalized The General Relativity Using Generalized Lorentz Transformation

Generalized The General Relativity Using Generalized Lorentz Transformation P P P P IJISET - Inernaional Journal of Innoaie Siene, Engineering & Tehnology, Vol. 3 Issue 4, April 6. www.ijise.om ISSN 348 7968 Generalized The General Relaiiy Using Generalized Lorenz Transformaion

More information

κt π = (5) T surrface k BASELINE CASE

κt π = (5) T surrface k BASELINE CASE II. BASELINE CASE PRACICAL CONSIDERAIONS FOR HERMAL SRESSES INDUCED BY SURFACE HEAING James P. Blanhard Universi of Wisonsin Madison 15 Engineering Dr. Madison, WI 5376-169 68-63-391 blanhard@engr.is.edu

More information

An Inventory Model for Weibull Time-Dependence. Demand Rate with Completely Backlogged. Shortages

An Inventory Model for Weibull Time-Dependence. Demand Rate with Completely Backlogged. Shortages Inernaional Mahemaial Forum, 5, 00, no. 5, 675-687 An Invenory Model for Weibull Time-Dependene Demand Rae wih Compleely Baklogged Shorages C. K. Tripahy and U. Mishra Deparmen of Saisis, Sambalpur Universiy

More information

Inverted Pendulum-type Personal Mobility Considering Human Vibration Sensitivity

Inverted Pendulum-type Personal Mobility Considering Human Vibration Sensitivity (IJACSA) Inernaional Journal of Advaned Compuer Siene and Appliaions, Vol. 5, No. 3, 14 Invered Pendulum-ype Personal Mobiliy Considering Human Vibraion Sensiiviy Misaki Masuda Shool of Siene for Open

More information

Concrete damaged plasticity model

Concrete damaged plasticity model Conree damaged asiiy model Conree damaged asiiy model is a maerial model for he analysis of onree sruures mainly under dynami loads suh as earhquakes(only aes an be analyzed under he dynami loads like

More information

Chapter 8 The Complete Response of RL and RC Circuits

Chapter 8 The Complete Response of RL and RC Circuits Chaper 8 he Complee Response of R and RC Ciruis Exerises Ex 8.3-1 Before he swih loses: Afer he swih loses: 2 = = 8 Ω so = 8 0.05 = 0.4 s. 0.25 herefore R ( ) Finally, 2.5 ( ) = o + ( (0) o ) = 2 + V for

More information

Boyce/DiPrima 9 th ed, Ch 6.1: Definition of. Laplace Transform. In this chapter we use the Laplace transform to convert a

Boyce/DiPrima 9 th ed, Ch 6.1: Definition of. Laplace Transform. In this chapter we use the Laplace transform to convert a Boye/DiPrima 9 h ed, Ch 6.: Definiion of Laplae Transform Elemenary Differenial Equaions and Boundary Value Problems, 9 h ediion, by William E. Boye and Rihard C. DiPrima, 2009 by John Wiley & Sons, In.

More information

EE40 Summer 2005: Lecture 2 Instructor: Octavian Florescu 1. Measuring Voltages and Currents

EE40 Summer 2005: Lecture 2 Instructor: Octavian Florescu 1. Measuring Voltages and Currents Announemens HW # Due oday a 6pm. HW # posed online oday and due nex Tuesday a 6pm. Due o sheduling onflis wih some sudens, lasses will resume normally his week and nex. Miderm enaively 7/. EE4 Summer 5:

More information

Mahgoub Transform Method for Solving Linear Fractional Differential Equations

Mahgoub Transform Method for Solving Linear Fractional Differential Equations Mahgoub Transform Mehod for Solving Linear Fraional Differenial Equaions A. Emimal Kanaga Puhpam 1,* and S. Karin Lydia 2 1* Assoiae Professor&Deparmen of Mahemais, Bishop Heber College Tiruhirappalli,

More information

Energy Momentum Tensor for Photonic System

Energy Momentum Tensor for Photonic System 018 IJSST Volume 4 Issue 10 Prin ISSN : 395-6011 Online ISSN : 395-60X Themed Seion: Siene and Tehnology Energy Momenum Tensor for Phooni Sysem ampada Misra Ex-Gues-Teaher, Deparmens of Eleronis, Vidyasagar

More information

The Contradiction within Equations of Motion with Constant Acceleration

The Contradiction within Equations of Motion with Constant Acceleration The Conradicion wihin Equaions of Moion wih Consan Acceleraion Louai Hassan Elzein Basheir (Daed: July 7, 0 This paper is prepared o demonsrae he violaion of rules of mahemaics in he algebraic derivaion

More information

3. Differential Equations

3. Differential Equations 3. Differenial Equaions 3.. inear Differenial Equaions of Firs rder A firs order differenial equaion is an equaion of he form d() d ( ) = F ( (),) (3.) As noed above, here will in general be a whole la

More information

5.2 Design for Shear (Part I)

5.2 Design for Shear (Part I) 5. Design or Shear (Par I) This seion overs he ollowing opis. General Commens Limi Sae o Collapse or Shear 5..1 General Commens Calulaion o Shear Demand The objeive o design is o provide ulimae resisane

More information

Capacitance and Inductance. The Capacitor

Capacitance and Inductance. The Capacitor apaiane and Induane OUTINE apaiors apaior volage, urren, power, energy Induors eure 9, 9/9/5 Reading Hambley haper 3 (A) EE4 Fall 5 eure 9, Slide The apaior Two onduors (a,b) separaed by an insulaor: differene

More information

A New Formulation of Electrodynamics

A New Formulation of Electrodynamics . Eleromagnei Analysis & Appliaions 1 457-461 doi:1.436/jemaa.1.86 Published Online Augus 1 hp://www.sirp.org/journal/jemaa A New Formulaion of Elerodynamis Arbab I. Arbab 1 Faisal A. Yassein 1 Deparmen

More information

PHYS-3301 Lecture 5. Chapter 2. Announcement. Sep. 12, Special Relativity. What about y and z coordinates? (x - direction of motion)

PHYS-3301 Lecture 5. Chapter 2. Announcement. Sep. 12, Special Relativity. What about y and z coordinates? (x - direction of motion) Announemen Course webpage hp://www.phys.u.edu/~slee/33/ Tebook PHYS-33 Leure 5 HW (due 9/4) Chaper, 6, 36, 4, 45, 5, 5, 55, 58 Sep., 7 Chaper Speial Relaiiy. Basi Ideas. Consequenes of Einsein s Posulaes

More information

8.1. a) For step response, M input is u ( t) Taking inverse Laplace transform. as α 0. Ideal response, K c. = Kc Mτ D + For ramp response, 8-1

8.1. a) For step response, M input is u ( t) Taking inverse Laplace transform. as α 0. Ideal response, K c. = Kc Mτ D + For ramp response, 8-1 8. a For ep repone, inpu i u, U Y a U α α Y a α α Taking invere Laplae ranform a α e e / α / α A α 0 a δ 0 e / α a δ deal repone, α d Y i Gi U i δ Hene a α 0 a i For ramp repone, inpu i u, U Soluion anual

More information

Derivation of longitudinal Doppler shift equation between two moving bodies in reference frame at rest

Derivation of longitudinal Doppler shift equation between two moving bodies in reference frame at rest Deriaion o longiudinal Doppler shi equaion beween wo moing bodies in reerene rame a res Masanori Sao Honda Eleronis Co., d., Oyamazuka, Oiwa-ho, Toyohashi, ihi 44-393, Japan E-mail: msao@honda-el.o.jp

More information

Basic solution to Heat Diffusion In general, one-dimensional heat diffusion in a material is defined by the linear, parabolic PDE

Basic solution to Heat Diffusion In general, one-dimensional heat diffusion in a material is defined by the linear, parabolic PDE New Meio e yd 5 ydrology Program Quaniaie Meods in ydrology Basi soluion o ea iffusion In general one-dimensional ea diffusion in a maerial is defined by e linear paraboli PE or were we assume a is defined

More information

Numerical Dispersion

Numerical Dispersion eview of Linear Numerical Sabiliy Numerical Dispersion n he previous lecure, we considered he linear numerical sabiliy of boh advecion and diffusion erms when approimaed wih several spaial and emporal

More information

Optimal Transform: The Karhunen-Loeve Transform (KLT)

Optimal Transform: The Karhunen-Loeve Transform (KLT) Opimal ransform: he Karhunen-Loeve ransform (KL) Reall: We are ineresed in uniary ransforms beause of heir nie properies: energy onservaion, energy ompaion, deorrelaion oivaion: τ (D ransform; assume separable)

More information

Problem Set 9 Due December, 7

Problem Set 9 Due December, 7 EE226: Random Proesses in Sysems Leurer: Jean C. Walrand Problem Se 9 Due Deember, 7 Fall 6 GSI: Assane Gueye his problem se essenially reviews Convergene and Renewal proesses. No all exerises are o be

More information

Mass Transfer Coefficients (MTC) and Correlations I

Mass Transfer Coefficients (MTC) and Correlations I Mass Transfer Mass Transfer Coeffiiens (MTC) and Correlaions I 7- Mass Transfer Coeffiiens and Correlaions I Diffusion an be desribed in wo ways:. Deailed physial desripion based on Fik s laws and he diffusion

More information

Theory of! Partial Differential Equations-I!

Theory of! Partial Differential Equations-I! hp://users.wpi.edu/~grear/me61.hml! Ouline! Theory o! Parial Dierenial Equaions-I! Gréar Tryggvason! Spring 010! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and

More information

Mathematical Foundations -1- Choice over Time. Choice over time. A. The model 2. B. Analysis of period 1 and period 2 3

Mathematical Foundations -1- Choice over Time. Choice over time. A. The model 2. B. Analysis of period 1 and period 2 3 Mahemaial Foundaions -- Choie over Time Choie over ime A. The model B. Analysis of period and period 3 C. Analysis of period and period + 6 D. The wealh equaion 0 E. The soluion for large T 5 F. Fuure

More information

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities: Mah 4 Eam Review Problems Problem. Calculae he 3rd Taylor polynomial for arcsin a =. Soluion. Le f() = arcsin. For his problem, we use he formula f() + f () + f ()! + f () 3! for he 3rd Taylor polynomial

More information

Theory of! Partial Differential Equations!

Theory of! Partial Differential Equations! hp://www.nd.edu/~gryggva/cfd-course/! Ouline! Theory o! Parial Dierenial Equaions! Gréar Tryggvason! Spring 011! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and

More information

Peakon, pseudo-peakon, and cuspon solutions for two generalized Camassa-Holm equations

Peakon, pseudo-peakon, and cuspon solutions for two generalized Camassa-Holm equations JOURNAL OF MATHEMATICAL PHYSICS 54 13501 (013) Peakon pseudo-peakon and uspon soluions for wo generalized Camassa-Holm equaions Jibin Li 1a) and Zhijun Qiao 3a) 1 Deparmen of Mahemais Zhejiang Normal Universiy

More information

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms

More information

Solutions to Assignment 1

Solutions to Assignment 1 MA 2326 Differenial Equaions Insrucor: Peronela Radu Friday, February 8, 203 Soluions o Assignmen. Find he general soluions of he following ODEs: (a) 2 x = an x Soluion: I is a separable equaion as we

More information

SOME ISSUES ON INERTIA PROPULSION MECHANISMS USING TWO CONTRA-ROTATING MASSES

SOME ISSUES ON INERTIA PROPULSION MECHANISMS USING TWO CONTRA-ROTATING MASSES Преподавание ТММ УДК 61.1 С. G. PROVATIDIS SOE ISSUES ON INERTIA PROPULSION ECHANISS USING TWO CONTRA-ROTATING ASSES 1. INTRODUCTION Among several physial priniples ha are poenially appliable o produe

More information

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon 3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of

More information

Number of modes per unit volume of the cavity per unit frequency interval is given by: Mode Density, N

Number of modes per unit volume of the cavity per unit frequency interval is given by: Mode Density, N SMES404 - LASER PHYSCS (LECTURE 5 on /07/07) Number of modes per uni volume of he aviy per uni frequeny inerval is given by: 8 Mode Densiy, N (.) Therefore, energy densiy (per uni freq. inerval); U 8h

More information

6.2 Transforms of Derivatives and Integrals.

6.2 Transforms of Derivatives and Integrals. SEC. 6.2 Transforms of Derivaives and Inegrals. ODEs 2 3 33 39 23. Change of scale. If l( f ()) F(s) and c is any 33 45 APPLICATION OF s-shifting posiive consan, show ha l( f (c)) F(s>c)>c (Hin: In Probs.

More information

Fractional Method of Characteristics for Fractional Partial Differential Equations

Fractional Method of Characteristics for Fractional Partial Differential Equations Fracional Mehod of Characerisics for Fracional Parial Differenial Equaions Guo-cheng Wu* Modern Teile Insiue, Donghua Universiy, 188 Yan-an ilu Road, Shanghai 51, PR China Absrac The mehod of characerisics

More information

TRANSMISSION LINES AND WAVEGUIDES. Uniformity along the Direction of Propagation

TRANSMISSION LINES AND WAVEGUIDES. Uniformity along the Direction of Propagation TRANSMISSION LINES AND WAVEGUIDES Uniformi along he Direion of Propagaion Definiion: Transmission Line TL is he erm o desribe ransmission ssems wih wo or more mealli onduors eleriall insulaed from eah

More information

Oscillatory Motion Application of Tubular Linear Permanent Magnet Machine

Oscillatory Motion Application of Tubular Linear Permanent Magnet Machine Osillaory Moion Aliaion of Tubular Linear Permanen Magne Mahine Absra-Tubular linear ermanen magne mahine wih gas srings using for drilling aliaions is resened. This mahine has he advanage of direly ransmiing

More information

SIMULATION STUDY OF STOCHASTIC CHANNEL REDISTRIBUTION

SIMULATION STUDY OF STOCHASTIC CHANNEL REDISTRIBUTION Developmens in Business Simulaion and Experienial Learning, Volume 3, 3 SIMULATIO STUDY OF STOCHASTIC CHAEL REDISTRIBUTIO Yao Dong-Qing Towson Universiy dyao@owson.edu ABSTRACT In his paper, we invesigae

More information

KINEMATICS IN ONE DIMENSION

KINEMATICS IN ONE DIMENSION KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec

More information

15. Vector Valued Functions

15. Vector Valued Functions 1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,

More information

CHAPTER 12 DIRECT CURRENT CIRCUITS

CHAPTER 12 DIRECT CURRENT CIRCUITS CHAPTER 12 DIRECT CURRENT CIUITS DIRECT CURRENT CIUITS 257 12.1 RESISTORS IN SERIES AND IN PARALLEL When wo resisors are conneced ogeher as shown in Figure 12.1 we said ha hey are conneced in series. As

More information

Zürich. ETH Master Course: L Autonomous Mobile Robots Localization II

Zürich. ETH Master Course: L Autonomous Mobile Robots Localization II Roland Siegwar Margaria Chli Paul Furgale Marco Huer Marin Rufli Davide Scaramuzza ETH Maser Course: 151-0854-00L Auonomous Mobile Robos Localizaion II ACT and SEE For all do, (predicion updae / ACT),

More information

Section 2.2 Charge and Current 2.6 b) The current direction is designated as the direction of the movement of positive charges.

Section 2.2 Charge and Current 2.6 b) The current direction is designated as the direction of the movement of positive charges. Chaper Soluions Secion. Inroducion. Curren source. Volage source. esisor.4 Capacior.5 Inducor Secion. Charge and Curren.6 b) The curren direcion is designaed as he direcion of he movemen of posiive charges..7

More information

Class Meeting # 10: Introduction to the Wave Equation

Class Meeting # 10: Introduction to the Wave Equation MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

Mathematical Modelling of the Solutions of Partial Differential Equation with Physical Implications

Mathematical Modelling of the Solutions of Partial Differential Equation with Physical Implications IOSR Jornal of Mahemais (IOSR-JM) e-issn: 78-578, p-issn: 39-765X. Volme 4, Isse 3 Ver. II (May - Jne 8), PP 86-97 www.iosrjornals.org Mahemaial Modelling of he Solions of Parial Differenial Eqaion wih

More information

Idealize Bioreactor CSTR vs. PFR... 3 Analysis of a simple continuous stirred tank bioreactor... 4 Residence time distribution... 4 F curve:...

Idealize Bioreactor CSTR vs. PFR... 3 Analysis of a simple continuous stirred tank bioreactor... 4 Residence time distribution... 4 F curve:... Idealize Bioreaor CSTR vs. PFR... 3 Analysis of a simple oninuous sirred ank bioreaor... 4 Residene ime disribuion... 4 F urve:... 4 C urve:... 4 Residene ime disribuion or age disribuion... 4 Residene

More information

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More information

Chapter 8 The Complete Response of RL and RC Circuits

Chapter 8 The Complete Response of RL and RC Circuits Chaper 8 The Complee Response of RL and RC Circuis Seoul Naional Universiy Deparmen of Elecrical and Compuer Engineering Wha is Firs Order Circuis? Circuis ha conain only one inducor or only one capacior

More information

Generalized electromagnetic energy-momentum tensor and scalar curvature of space at the location of charged particle

Generalized electromagnetic energy-momentum tensor and scalar curvature of space at the location of charged particle Generalized eleromagnei energy-momenum ensor and salar urvaure of spae a he loaion of harged parile A.L. Kholmeskii 1, O.V. Missevih and T. Yarman 3 1 Belarus Sae Universiy, Nezavisimosi Avenue, 0030 Minsk,

More information

MATH 31B: MIDTERM 2 REVIEW. x 2 e x2 2x dx = 1. ue u du 2. x 2 e x2 e x2] + C 2. dx = x ln(x) 2 2. ln x dx = x ln x x + C. 2, or dx = 2u du.

MATH 31B: MIDTERM 2 REVIEW. x 2 e x2 2x dx = 1. ue u du 2. x 2 e x2 e x2] + C 2. dx = x ln(x) 2 2. ln x dx = x ln x x + C. 2, or dx = 2u du. MATH 3B: MIDTERM REVIEW JOE HUGHES. Inegraion by Pars. Evaluae 3 e. Soluion: Firs make he subsiuion u =. Then =, hence 3 e = e = ue u Now inegrae by pars o ge ue u = ue u e u + C and subsiue he definiion

More information

Traveling Waves. Chapter Introduction

Traveling Waves. Chapter Introduction Chaper 4 Traveling Waves 4.1 Inroducion To dae, we have considered oscillaions, i.e., periodic, ofen harmonic, variaions of a physical characerisic of a sysem. The sysem a one ime is indisinguishable from

More information

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid

More information

where the coordinate X (t) describes the system motion. X has its origin at the system static equilibrium position (SEP).

where the coordinate X (t) describes the system motion. X has its origin at the system static equilibrium position (SEP). Appendix A: Conservaion of Mechanical Energy = Conservaion of Linear Momenum Consider he moion of a nd order mechanical sysem comprised of he fundamenal mechanical elemens: ineria or mass (M), siffness

More information

Circuit Variables. AP 1.1 Use a product of ratios to convert two-thirds the speed of light from meters per second to miles per second: 1 ft 12 in

Circuit Variables. AP 1.1 Use a product of ratios to convert two-thirds the speed of light from meters per second to miles per second: 1 ft 12 in Circui Variables 1 Assessmen Problems AP 1.1 Use a produc of raios o conver wo-hirds he speed of ligh from meers per second o miles per second: ( ) 2 3 1 8 m 3 1 s 1 cm 1 m 1 in 2.54 cm 1 f 12 in 1 mile

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

The motions of the celt on a horizontal plane with viscous friction

The motions of the celt on a horizontal plane with viscous friction The h Join Inernaional Conference on Mulibody Sysem Dynamics June 8, 18, Lisboa, Porugal The moions of he cel on a horizonal plane wih viscous fricion Maria A. Munisyna 1 1 Moscow Insiue of Physics and

More information

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals

More information

LAPLACE TRANSFORM AND TRANSFER FUNCTION

LAPLACE TRANSFORM AND TRANSFER FUNCTION CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dep. of Chemical and Biological Engineering 5-1 Road Map of he Lecure V Laplace Transform and Transfer funcions

More information

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k Challenge Problems DIS 03 and 0 March 6, 05 Choose one of he following problems, and work on i in your group. Your goal is o convince me ha your answer is correc. Even if your answer isn compleely correc,

More information

MATH 128A, SUMMER 2009, FINAL EXAM SOLUTION

MATH 128A, SUMMER 2009, FINAL EXAM SOLUTION MATH 28A, SUMME 2009, FINAL EXAM SOLUTION BENJAMIN JOHNSON () (8 poins) [Lagrange Inerpolaion] (a) (4 poins) Le f be a funcion defined a some real numbers x 0,..., x n. Give a defining equaion for he Lagrange

More information

Differential Equations

Differential Equations Mah 21 (Fall 29) Differenial Equaions Soluion #3 1. Find he paricular soluion of he following differenial equaion by variaion of parameer (a) y + y = csc (b) 2 y + y y = ln, > Soluion: (a) The corresponding

More information

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem)

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem) Week 1 Lecure Problems, 5 Wha if somehing oscillaes wih no obvious spring? Wha is ω? (problem se problem) Sar wih Try and ge o SHM form E. Full beer can in lake, oscillaing F = m & = ge rearrange: F =

More information

Molecular Motion in Isotropic Turbulence

Molecular Motion in Isotropic Turbulence Moleular Moion in Isoropi Turbulene Jing Fan, Jian-Zheng Jiang, and Fei Fei Laboraory of High Temperaure Gas Dynamis, Insiue of Mehanis Chinese Aademy of Sienes, Being 9, China Absra Moleular moion in

More information

Lecture 20: Riccati Equations and Least Squares Feedback Control

Lecture 20: Riccati Equations and Least Squares Feedback Control 34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he

More information

New Oscillation Criteria For Second Order Nonlinear Differential Equations

New Oscillation Criteria For Second Order Nonlinear Differential Equations Researh Inveny: Inernaional Journal Of Engineering And Siene Issn: 78-47, Vol, Issue 4 (Feruary 03), Pp 36-4 WwwResearhinvenyCom New Osillaion Crieria For Seond Order Nonlinear Differenial Equaions Xhevair

More information

8. Basic RL and RC Circuits

8. Basic RL and RC Circuits 8. Basic L and C Circuis This chaper deals wih he soluions of he responses of L and C circuis The analysis of C and L circuis leads o a linear differenial equaion This chaper covers he following opics

More information

A Numerical Hydraulic Fracture Model Using the Extended Finite Element Method

A Numerical Hydraulic Fracture Model Using the Extended Finite Element Method nernaional Conferene on Mehanial and ndusrial Engineering (CME'2013) Augus 28-29, 2013 Penang (Malaysia) A Numerial Hydrauli Fraure Model Using he Exended Finie Elemen Mehod Ashkan. Mahdavi, and Soheil.

More information

The Arcsine Distribution

The Arcsine Distribution The Arcsine Disribuion Chris H. Rycrof Ocober 6, 006 A common heme of he class has been ha he saisics of single walker are ofen very differen from hose of an ensemble of walkers. On he firs homework, we

More information

THE MYSTERY OF STOCHASTIC MECHANICS. Edward Nelson Department of Mathematics Princeton University

THE MYSTERY OF STOCHASTIC MECHANICS. Edward Nelson Department of Mathematics Princeton University THE MYSTERY OF STOCHASTIC MECHANICS Edward Nelson Deparmen of Mahemaics Princeon Universiy 1 Classical Hamilon-Jacobi heory N paricles of various masses on a Euclidean space. Incorporae he masses in he

More information

AN INVENTORY MODEL FOR DETERIORATING ITEMS WITH EXPONENTIAL DECLINING DEMAND AND PARTIAL BACKLOGGING

AN INVENTORY MODEL FOR DETERIORATING ITEMS WITH EXPONENTIAL DECLINING DEMAND AND PARTIAL BACKLOGGING Yugoslav Journal of Operaions Researh 5 (005) Number 77-88 AN INVENTORY MODEL FOR DETERIORATING ITEMS WITH EXPONENTIAL DECLINING DEMAND AND PARTIAL BACKLOGGING Liang-Yuh OUYANG Deparmen of Managemen Sienes

More information

Teacher Quality Policy When Supply Matters: Online Appendix

Teacher Quality Policy When Supply Matters: Online Appendix Teaher Qualiy Poliy When Supply Maers: Online Appendix Jesse Rohsein July 24, 24 A Searh model Eah eaher draws a single ouside job offer eah year. If she aeps he offer, she exis eahing forever. The ouside

More information

t + t sin t t cos t sin t. t cos t sin t dt t 2 = exp 2 log t log(t cos t sin t) = Multiplying by this factor and then integrating, we conclude that

t + t sin t t cos t sin t. t cos t sin t dt t 2 = exp 2 log t log(t cos t sin t) = Multiplying by this factor and then integrating, we conclude that ODEs, Homework #4 Soluions. Check ha y ( = is a soluion of he second-order ODE ( cos sin y + y sin y sin = 0 and hen use his fac o find all soluions of he ODE. When y =, we have y = and also y = 0, so

More information

Nevertheless, there are well defined (and potentially useful) distributions for which σ 2

Nevertheless, there are well defined (and potentially useful) distributions for which σ 2 M. Meseron-Gibbons: Bioalulus, Leure, Page. The variane. More on improper inegrals In general, knowing only he mean of a isribuion is no as useful as also knowing wheher he isribuion is lumpe near he mean

More information

MOMENTUM CONSERVATION LAW

MOMENTUM CONSERVATION LAW 1 AAST/AEDT AP PHYSICS B: Impulse and Momenum Le us run an experimen: The ball is moving wih a velociy of V o and a force of F is applied on i for he ime inerval of. As he resul he ball s velociy changes

More information

The Asymptotical Behavior of Probability Measures for the Fluctuations of Stochastic Models

The Asymptotical Behavior of Probability Measures for the Fluctuations of Stochastic Models The Asympoial Behavior of Probabiliy Measures for he Fluuaions of Sohasi Models JUN WANG CUINING WEI Deparmen of Mahemais College of Siene Beijing Jiaoong Universiy Beijing Jiaoong Universiy Beijing 44

More information

Economics 202 (Section 05) Macroeconomic Theory Practice Problem Set 7 Suggested Solutions Professor Sanjay Chugh Fall 2013

Economics 202 (Section 05) Macroeconomic Theory Practice Problem Set 7 Suggested Solutions Professor Sanjay Chugh Fall 2013 Deparmen of Eonomis Boson College Eonomis 0 (Seion 05) Maroeonomi Theory Praie Problem Se 7 Suggesed Soluions Professor Sanjay Chugh Fall 03. Lags in Labor Hiring. Raher han supposing ha he represenaive

More information

1.054/1.541 Mechanics and Design of Concrete Structures (3-0-9) Outline 5 Creep and Shrinkage Deformation

1.054/1.541 Mechanics and Design of Concrete Structures (3-0-9) Outline 5 Creep and Shrinkage Deformation 1.54/1.541 Mehanis and Design of Conree ruures pring 24 Prof. Oral Buyukozurk Massahuses Insiue of Tehnology Ouline 5 1.54/1.541 Mehanis and Design of Conree ruures (3--9 Ouline 5 and hrinkage Deformaion

More information

Finite Element Analysis of Structures

Finite Element Analysis of Structures KAIT OE5 Finie Elemen Analysis of rucures Mid-erm Exam, Fall 9 (p) m. As shown in Fig., we model a russ srucure of uniform area (lengh, Area Am ) subjeced o a uniform body force ( f B e x N / m ) using

More information

The Special Theory of Relativity Chapter II

The Special Theory of Relativity Chapter II The Speial Theory of Relaiiy Chaper II 1. Relaiisi Kinemais. Time dilaion and spae rael 3. Lengh onraion 4. Lorenz ransformaions 5. Paradoes? Simulaneiy/Relaiiy If one obserer sees he eens as simulaneous,

More information

1. VELOCITY AND ACCELERATION

1. VELOCITY AND ACCELERATION 1. VELOCITY AND ACCELERATION 1.1 Kinemaics Equaions s = u + 1 a and s = v 1 a s = 1 (u + v) v = u + as 1. Displacemen-Time Graph Gradien = speed 1.3 Velociy-Time Graph Gradien = acceleraion Area under

More information

The fundamental mass balance equation is ( 1 ) where: I = inputs P = production O = outputs L = losses A = accumulation

The fundamental mass balance equation is ( 1 ) where: I = inputs P = production O = outputs L = losses A = accumulation Hea (iffusion) Equaion erivaion of iffusion Equaion The fundamenal mass balance equaion is I P O L A ( 1 ) where: I inpus P producion O oupus L losses A accumulaion Assume ha no chemical is produced or

More information

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3 and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or

More information

Math 334 Fall 2011 Homework 11 Solutions

Math 334 Fall 2011 Homework 11 Solutions Dec. 2, 2 Mah 334 Fall 2 Homework Soluions Basic Problem. Transform he following iniial value problem ino an iniial value problem for a sysem: u + p()u + q() u g(), u() u, u () v. () Soluion. Le v u. Then

More information

Predator - Prey Model Trajectories and the nonlinear conservation law

Predator - Prey Model Trajectories and the nonlinear conservation law Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories

More information

2. Nonlinear Conservation Law Equations

2. Nonlinear Conservation Law Equations . Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear

More information

Lorentz Transformation Properties of Currents for the Particle-Antiparticle Pair Wave Functions

Lorentz Transformation Properties of Currents for the Particle-Antiparticle Pair Wave Functions Open Aess Library Journal 17, Volume 4, e373 ISSN Online: 333-971 ISSN Prin: 333-975 Lorenz Transformaion Properies of Currens for he Parile-Aniparile Pair Wave Funions Raja Roy Deparmen of Eleronis and

More information

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

More information

Problem 1 / 25 Problem 2 / 10 Problem 3 / 15 Problem 4 / 30 Problem 5 / 20 TOTAL / 100

Problem 1 / 25 Problem 2 / 10 Problem 3 / 15 Problem 4 / 30 Problem 5 / 20 TOTAL / 100 Deparmen of Applied Eonomis Johns Hopkins Universiy Eonomis 60 Maroeonomi Theory and Poliy Miderm Exam Suggesed Soluions Professor Sanjay Chugh Summer 0 NAME: The Exam has a oal of five (5) problems and

More information

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan Ground Rules PC11 Fundamenals of Physics I Lecures 3 and 4 Moion in One Dimension A/Prof Tay Seng Chuan 1 Swich off your handphone and pager Swich off your lapop compuer and keep i No alking while lecure

More information

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.

More information

Optimal Path Planning for Flexible Redundant Robot Manipulators

Optimal Path Planning for Flexible Redundant Robot Manipulators 25 WSEAS In. Conf. on DYNAMICAL SYSEMS and CONROL, Venice, Ialy, November 2-4, 25 (pp363-368) Opimal Pah Planning for Flexible Redundan Robo Manipulaors H. HOMAEI, M. KESHMIRI Deparmen of Mechanical Engineering

More information

Sub Module 2.6. Measurement of transient temperature

Sub Module 2.6. Measurement of transient temperature Mechanical Measuremens Prof. S.P.Venkaeshan Sub Module 2.6 Measuremen of ransien emperaure Many processes of engineering relevance involve variaions wih respec o ime. The sysem properies like emperaure,

More information

Dynamic Analysis of Damped Driven Pendulum using Laplace Transform Method

Dynamic Analysis of Damped Driven Pendulum using Laplace Transform Method , ISSN 0974-570X (Online), ISSN 0974-578 (Prin), Vol. 6; Issue No. 3; Year 05, Copyrigh 05 by CESER PUBLICATIONS Dynamic Analysis of Damped Driven Pendulum using Laplace Transform Mehod M.C. Agarana and

More information