Residual Modes on Non-linear Resonant Decay Method

Size: px
Start display at page:

Download "Residual Modes on Non-linear Resonant Decay Method"

Transcription

1 Residual Modes on Non-linea Resonant Decay Method D Mehdi Samast, Pofesso Jan R. Wight ABSRAC Non-linea Resonant Decay method (NL-RDM) addesses the identification of multi-degee of feedom non-linea systems. his method offes a pactical appoach to the identification of lumped paamete and continuous systems by poducing a non-linea extension of the classical linea modal model. he method is intoduced and its potential as a pactical identification appoach explained. his pape is concened with the inclusion of esidual modes, above and below any egion of inteest. Any stuctue that is not suppoted to eath will have igid body modes, being modes that have a natual fequency at zeo Hetz; a good example is an aicaft. hese modes ae sometimes ignoed but fo a complete mathematical model obtained fom NL-RDM the effects of igid body modes must be analysed. NL-RDM elies on a method of cuve fitting to geneate modal chaacteistics, in tems of foce, displacement, velocity and acceleation in modal space. he lowe esidual egion was obseved to contain igid body modes, and these wee obseved to affect significantly the esults in the egion of inteest. Monitoing this effect was shown to be too difficult cuently fo the NL-RDM, given technological estictions. A Mass Substitution method was geneated to model the system esponse moe accuately. Its accuacy was demonstated though case studies. Most systems also contain an indefinite numbe of esidual modes, these being modes that occu at fequencies geate than the fequency ange of inteest. he section will attempt to detemine the effects of excluding esidual modes on the accuacy of NL-RDM. he effect of the uppe esidual aea on the modelling of the system within the egion of inteest was obseved though seveal case studies. Simple validation methods have been consideed to peseve the accuacy of the model, taking into account both hamonic and coupled esidual modes and thei effects. Keywods: NL-RDM, Residual Modes, Rigid Body Mode, Out of Range modes. Manuscipt eceived on Mach, 3,. D. Mehdi Samast Autho is with Islamic Azad Univesity, East ehan Banch, ehan, Ian, phone: , Membe of IAENG, msamast@yahoo.com Pofesso Jan R. Wight is with he Univesity of Mancheste, Mancheste, M3 9PL, UK, jan.wight@mancheste.ac.uk I. MEHODOLOGY FOR NL-RDM he method is a development of the Resonant Decay Method, the Restoing Foce Suface method based in modal space and the Foce Appopiation methodology [, ]. It makes use of an appopiated excitation vecto applied in a bust to educe the numbe of modes equied in the cuve fit. In this section, the mathematical fom of the extended modal model fo non-linea systems will be descibed. Conside a lumped paamete system o a continuous one that has been discedited as N degees of feedom (DOF). he equation of motion in physical space fo a dynamic system including stiffness non-lineaities is M x C x K x K x x F(t) () L whee [M] is the mass matix, [C] is the damping matix, [K L ] is the linea stiffness matix, [K NL ] is the non-linea stiffness matix, {F(t)} is the vecto of applied nodal foces, and {x(t)} is the vecto of physical displacements. he tansfomation between physical and modal space is defined by N NL { x( t)} { } p [ ]{ p( t)} whee {p(t)} is a vecto of modal amplitudes, and [] is the modal matix of the N modes {}, =, N of the undelying linea system, which may be obtained by solving the classical eigenvalue poblem fo un-damped fee vibation. he coesponding un-damped natual fequencies ae ( =, N). In pactice, the numbe of degees of feedom equied to model the system with easonable accuacy can be educed by consideing only those modes NR << N having natual fequencies in the fequency ange of inteest, o that ae consideed impotant in the esponse. In this case, the modal matix [ ] is educed to dimension N by NR and the vecto of modal amplitudes {p(t)} to dimension NR by. Substituting the modal expansion into the system equations of motion and pe-multiplying by [] yields: M p( t) C p K K p( t) F L NL Using the othogonality of the modes, this equation of motion in modal space becomes () (3) ISSN: (Pint); ISSN: (Online)

2 m p c p [ k ] p( t) f f (4) L whee the modal mass matix [m] and linea modal stiffness matix [k L ] ae diagonal and the modal damping matix [c] is diagonal fo popotionally damped systems. {f} is the applied modal foce vecto and {f NL } is the vecto of nonlinea estoing foces in modal space. Clealy, in ode to be able to pefom this tansfomation, the modal matix [] is assumed to be known fom low foce level esults. Fo a paticula mode of a popotionally damped non-linea system, equation (4) is in the fom of a single degee of feedom system, namely m p c p k p fnl f,... NR (5) whee p is the th modal displacement, m, c and k ae the th mode modal mass, damping and stiffness, and f is the coesponding applied modal foce. he modal mass and stiffness fo the linea system ae elated by the un-damped natual fequency fom the equation k m [3]. Nonpopotional damping would lead to the pesence of modal damping coupling tems [4]. he tem f NL efes to the th mode non-linea modal estoing foce and in geneal includes othe modal co-odinates to allow fo non-linea coss coupling f ( p, p,.., p NR 3 s i A i ( s, s ) p NR i s ) NL NR NR s ts i 3i j NL A i j s, t p p whee the A tem is the amplitude of the non-linea tem in the th mode equation that is associated with the espective modal displacements p s and p t (fo modes identified by the indices s and t), aised to the powes defined by i and j. Fo example, the diect cubic non-linea stiffness tem associated with the th mode would be A p 3 (s=, t=, i=3, j=) and a typical coss-coupling tem between modes and q might be A p p q (s=, t=q, i=, j=). Note, in this case, that the set of possible basis functions is esticted to tems of 3 d ode involving no moe than two modes. Also, the basis functions in equation (6) ae polynomials wheeas non-polynomial tems may be used to cate fo diffeent types of non-lineaity (e.g. discontinuous) [5]. hee is no eason in pinciple why non-linea modal velocity tems cannot also be included in the extended modal model to cope with damping non-lineaity. A. SIMPLE EXAMPLE Non-lineaity is not easy to define in a geneal manne and so the NL-RDM cannot be applied to a geneal case. hus, it is ecommended that the consideation of the method stat with a simple example such as a two degee of feedom system with some kind of non-lineaity, which could be the common cubic stiffness, to take advantage of a elatively simple case. Some othe type of non-lineaity such i s j t (6) as quadatic damping, fiction and backlash will not be investigated in this pape. II. RIGID BODY AND OU OF RANGE MODES he pevious pape [6] attempted to detemine the effects the use of impefect data has on NL-RDM, assuming use of the coect model. his pape will assume the use of pefect data and test pocess, but will analyse some of the effects of using an incoect model. he analysis will include the effects of igid body modes, of poblems associated with hamonics, and of the esidual pat of the data. he effects will be analysed in tems of deived system paametes and the FRF. Any stuctue that is not suppoted to eath will have igid body modes, being modes that have a natual fequency at zeo Hetz; a good example is an aicaft. hese modes ae sometimes ignoed but fo a complete mathematical model obtained fom NL-RDM the effects of igid body modes must be analysed. NL-RDM elies on a method of cuve fitting to geneate modal chaacteistics, in tems of foce, displacement, velocity and acceleation in modal space. Most systems also contain an indefinite numbe of esidual modes, these being modes that occu at fequencies geate than the fequency ange of inteest. his chapte will attempt to detemine the effects of excluding esidual modes on the accuacy of NL-RDM. Finally, hamonic effects will be analysed. Specifically a system with a cubic stiffness non-lineaity will be obseved, when one mode natual fequency occus at pecisely 3 times the natual fequency of the othe mode. A. IDENIFYING RIGID BODY MODES Rigid body motion occus in many systems and affects the dynamic behaviou, so should be investigated. A igid body mode is a mode with a natual fequency of zeo Hetz and mode shapes which nomally consist of pue tanslation o otation; igid body modes occu in a pefect fee-fee stuctue. Expeimentally, igid body motion is obseved to cause no stain in any elastic element in the system. Rigid body behaviou can include six types of motion, which ae tanslation and otation in 3 dimensions. he way in which the NL-RDM would appoach the analysis of igid body modes needs to be discussed, as a complete model is equied in ode fo the dynamic behaviou of the system to be investigated. he poblems encounteed by this appoach give ise to thee pactical methods of analysing stuctues with igid body modes. Fistly, the mode is identified using the NL-RDM in the same way as fo any othe mode. Secondly, the igid body mode can simply be ignoed. he effects of this on accuate system epesentation ae analysed by compaison between pactical esults fo models with and without igid body modes. hidly, the igid body mode is simulated by a substituted modal mass element without the need fo any specific identification. ISSN: (Pint); ISSN: (Online)

3 MMIF Eigenvalue Phase FRF FRF A.. NON-LINEAR MDOF MODELS WIH RIGID MODY MODES In this section, two systems ae intoduced fo the evaluation of the pactical method, i.e. to detemine the effectiveness of appoximating a fee-fee system by a lightly spung fixed-fixed system. he thee degee of feedom lumped paamete system with stiffness nonlineaity between the second and thid masses (fee-fee stuctue) in Figue, is compaed to anothe stuctue, identical except fo being suppoted at both ends on soft spings (pactical stuctue) in Figue. stuctues intoduced in the pevious section, ae compaed in Figues 3 and FRFs phases in Figue 4. he FRFs show stong ageement, indicating the accuacy of the igid body mode appoximation and justifying the analysis of the suppoted stuctue fo igid body mode identification. Figue 5 shows the MMIF [7] fo the 3 degee of feedom pactical stuctue, whee the igid body mode on soft suppots is included. Clealy, appopiated foce pattens could be obtained fom such an analysis. 3 5 Pactical Pefect 5 5 Figue : Schematic Diagam of 3 Degee of Feedom Feefee System with Stiffness Non-lineaity (pefect fee-fee stuctue) 5 5 Fequency Figue 3: H FRFs fo the Pefect fee-fee Stuctue and the Pactical System 5 Pactical Pefect Figue : Schematic Diagam of 3 Degee of Feedom Fee- Fee System with Stiffness Non-lineaity but suppoted on Soft Spings (pactical stuctue) he FRF esults fo the pefect fee-fee and pactical systems ae shown in Figues 3 and 4 espectively. Compaison between the fist and second modes, i.e. the igid body mode and the fist flexible mode, gives an indication of the suitability of the pactical system fo appoximating the fee-fee system. At a atio of / the fist mode is consideed close enough to zeo Hetz to not have any significant effect on the flexible mode; in this example, the igid body mode natual fequency is.4 Hz and the fist flexible mode is 5.6 Hz. Fo bette compaison between the pefect and pactical system, the mode shapes o modal matix fo the two cases ae as follows.5 ModalMati x (7) Pefect.5 and.9967 ModalMatix P actical Clealy, the suppot has little effect on the modal chaacteistics. he FRFs of the pefect and pactical (8) Fequency Figue 4: H FRFs Phases fo the Pefect fee-fee Stuctue and the Pactical System Fequency (Hz) Figue 5: MMIF fo the 3DOF Pactical Fee-fee System A. IDENIFICAION OF A RIGID BODY MODE USING NL-RDM Having confimed that it is acceptable to conside a fee-fee stuctue to be suppoted on soft spings, one possibility fo identifying the igid body behaviou is to simply apply the NL-RDM [8] appoach just as fo a flexible mode. he NL-RDM is pefectly suited to identify a fee-fee stuctue by simulation, but poblems aise when ISSN: (Pint); ISSN: (Online)

4 H (m/s /N) H (m/s /N) seeking to identify a stuctue in pactice, i.e. expeimentally. Seveal poblems hinde the analysis of igid body modes pactically as the igid body mode may need to be investigated with vey low fequencies. his may ceate a significant poblem with noise, as the level of noise is close to the low acceleation levels expeienced at low fequencies. he vey low fequencies ceate futhe poblems expeimentally when geneating foce and monitoing esponse. Piezo-electic tansduces ae not accuate at low fequency levels as the chaacteistics olloff at low fequencies. Howeve, it would be possible to use acceleometes that opeate down to DC (e.g. piezo-esistive acceleometes) o even displacement tansduces. Commonly, shakes would be used to geneate foce and acceleometes used to monito the motion, so the ability to geneate acceptably low levels of foce and monito the low levels of motion, and hence obtain accuate esults, is limited by the capability of these pieces of equipment. It is well known that it is difficult to use an electodynamic shake to geneate a sinusoidal foce at low fequency. Futhe poblems ae encounteed when using the NL- RDM to identify the igid body mode. Poblems may aise when pepaing fo foce appopiation. A poo MMIF, specifically a pooly identified igid body mode causes an inaccuate foce vecto and poo tuning of the system. hese poblems ae citical with a lowe flexible mode fequency. Fo example, with a fist flexible mode aound 5 Hz, the igid mode fequency would be aound.5 Hz and the pocedue of detecting and identifying the mode becomes much moe difficult than when a stuctue s fist flexible mode fequency appeas aound Hz and the igid body fequency could be moe like Hz. Fo such a stiffe stuctue, thee would be no paticula poblem in using the NL-RDM to identify the igid body mode linea chaacteistics; no non-linea behaviou would be expected. ypically, the instumentation poblems mentioned above become wose below about 4 Hz, coesponding to a flexible mode at 4 Hz if the igid / flexible mode fequency atio is %. heoetically, using impoved foce geneation / ecoding, with bette shakes and acceleometes could poduce bette, moe accuate esults at vey low fequencies. echnologies ae aleady available o being developed that could impove the accuacy of these devices (e.g. lase based measuement systems). Clealy, impoved low-level accuacy of components would have a significant effect on the difficulties associated with igid body mode expeimentation. With accuate enough components, the igid body mode could theoetically be analysed by the NL- RDM in the same way as any othe mode. Howeve, if diect identification of vey low fequency igid body modes is too difficult to pefom accuately, maybe thee ae othe appoaches? Pehaps the igid body mode could simply be ignoed o a theoetically based estimate included? A.3 NEGLECING HE RIGID BODY MODE IN HE Non-IDEALIZED MODEL Anothe method that could be employed to analyse the system esponse in a non-linea system with a igid body mode is to entiely ignoe the igid body mode itself. Fo example, in the flutte analysis of an aicaft, often only the flexible modes ae included. he non-linea system is then analysed like any othe system, in tems of its flexible modes. he ease of this method is clealy its geatest advantage, but what may be a significant mode is missed out. Case studies ae caied out in this section to demonstate the accuacy of this appoach as compaed to identifying the igid body mode fom measued expeimental esults. he same system as intoduced in Figue is chosen fo analysis. Figues 6 and 7 show the calculated FRFs fo the system when the igid body mode is ignoed (Geen colou), and compaes them to esults fom a complete analysis, with the igid body mode on soft suppots included in the pactical system (Blue colou). While the geneal shape of the FRFs match aound the two flexible mode peaks, it is clea that esponse levels ae all somewhat inaccuate when the igid body mode is ignoed. he effect occus not only whee the igid body peak is not modelled, but at othe fequencies too, because the influence of the igid body mode extends acoss the entie fequency ange of the system. hus, while this method can be used as an appoximation to the tue system, its accuacy fo geneating system esponse and thus its application may be seiously limited. Figue 8 implements a flowchat, which indicates the way the compaison between systems is done Pactical neglected 5 5 Fequency (Hz) Figue 6: Compaisons of FRF (H ) Pactical neglected 5 5 Fequency (Hz) Figue 7: Compaisons of FRF (H ) ISSN: (Pint); ISSN: (Online)

5 a modal equation to allow fo the igid body mode contibution. he mass epesentation method involves analysing the igid body mode in tems of an added modal mass, that mass being deived fom a mass summation ove the stuctue, assuming a igid body mode shape. his method has the advantage of a moe complete system analysis, without making any majo changes o additions to the oiginal method. he appoach will be illustated fo the simple case of the 3 DOF system consideed ealie. he extension to moe complex poblems will be outlined late. he igid body modal mass M, which is also in a sense a mass summation, could be calculated as M * M * (9) Physical Figue 8: Compaison between Systems Rigid Body Ignoed o Included A.4 REPRESENING HE RIGID BODY MODE WIHOU SPECIFIC IDENIFICAION Due to the inaccuacy of measuement of low fequency igid body modes, a complete identification using NL-RDM may be consideed too inaccuate to be wokable; cetainly this is the case fo an aicaft whee the flexible modes occu at low fequencies (typically -5 Hz). heefoe a method is equied that will eadily allow fo the igid body mode effects in the identified model; fo example, in cuve fitting an FRF, it is common to include a lowe esidual tem in the model to appoximate the igid body effect [9]. In this section, a method is poposed that will simply add in whee the igid body mode shape fo the 3 DOF system is () and M Physical is the physical mass matix of the system. hus the igid body modal mass is M M * M * M 3 () o M M () M o moe geneally M 3 M M M (3) i otal hus, if the total mass of the stuctue is known by diect measuement o calculation, and if the igid body tanslation mode shape is assumed to have unit values, the igid body modal mass is simply equal to the total mass of the stuctue. In ode to include the igid body mode behaviou as an additional equation to accompany the identified flexible mode equations, then the appopiate modal foce and possibly modal mass and damping values should be estimated. By analogy to the analysis fo mass pefomed above, the igid body modal stiffness K will equal the summation of the suppot stiffnesses (e.g. bungees) and the modal damping C will be the summation of the suppot damping coefficients, both tansfomed into modal space. he igid body modal foce may be detemined though the physical foces to be applied to the system, which means that: F * F (4) Physical F F * F (5) F3 F F (6) F F3 ISSN: (Pint); ISSN: (Online)

6 o moe geneally F Fi Fotal (7) hus the igid body modal foce is simply equal to the summation of the physical foces applied to the stuctue. he final modal equation fo the igid body mode (including suppots if desied) will then be M q Cq Kq F (8) whee q is the igid body modal coodinate. his modal equation is simply added to those identified fo the nonlinea flexible modes. It should be noted that no non-linea modal couplings will exist between igid body and flexible modes, povided the suppot stiffnesses ae linea, because no non-linea foces will be geneated by igid body motion and foces geneated in motion of a flexible mode will not do any wok in the igid body mode. Fo implementation in the simulation pogamme, a system with 3 degee of feedom with cubic stiffness nonlineaity and one igid body mode is consideed. Figue shows the schematic system of the pactical stuctue. Fistly, only the igid body modal mass is included (i.e. effect of suppot stiffness is ignoed) and then the suppot effects ae included though the modal stiffness K. he accuacy is clealy impoved in compaison with the Neglected method, with the FRFs being mostly coincident. he only vaiation is seen in the epesentation of the igid body mode itself, aound which point the methodically calculated FRF levels off, giving a lowe esponse level aound the igid body mode. hus this appoach is a simple and pactical one to including igid body mode effects in the identified nonlinea model without needing to cay out difficult low fequency measuements. he appoach epesents the igid body effects accuately, povided the estimated mass is easonably accuate; indeed, the epesentation of the igid body modes is consideed to be moe accuate than using the lowe esidual tem in a cuve fit. If suppot stiffnesses wee only included to make the test possible, and a model of the tue fee-fee system equied, then only the igid body modal mass and foce need be included in the final model, as the suppot stiffness and damping effects can be ignoed. A.5 PRACICALIIES FOR MORE COMPLEX SRUCURES In eality a stuctue could have as many as six igid body modes. his will affect the mass epesentation method. Fo the 3 tanslation modes, the modal masses will all equal the total mass of the stuctue and the mode shapes will be equal tanslation of all points in the same diection, i.e. simple extension of the above appoach. hus, fo a heave mode, Heave M M M (9) i otal and so on. he othe igid body modes will involve igid body otation about the cente of mass (e.g. oll, pitch and yaw). In this case, the substitution would be an inetia epesentation, as opposed to being mass based fo a tanslation type mode. hus fo a pitch mode, it may be shown that Pitch M Ii I () otal whee I i is the pitch moment of inetia of the i th component of the stuctue about the cente of mass and I otal is the total pitch moment of inetia of the stuctue about the cente of mass. he pitch mode modal foce will be the sum of applied moments of each of the physical foces about the cente of mass. he coesponding otation mode shape must be a adian otation if the modal mass is the moment of inetia. In pactice, if the suppot stiffness effects ae to be included in the igid body modal equations, the values may be found by diect measuement of the sping stiffnesses o by detemining the igid body natual fequencies by exciting o distubing the stuctue into a igid body motion, even though the data may be noisy, and then calculating the modal stiffness fom the modal mass. Flexible modes would be identified as nomal using NL-RDM. B. UPPER RESIDUAL EFFECS It is possible to divide the ange of fequencies into thee egions. Of cental concen is the paticula egion of inteest, which will often include seveal impotant flexible and maybe non-linea modes. If this aea is consideed in isolation, it may be modelled inaccuately, as the enegy fom modes is less than expected in that egion. he peceding section has consideed the low fequency egion below the egion of inteest, containing one o moe igid body modes, and obseved the effect of its pesence on the egion of inteest; appoaches fo allowing fo the igid body effects wee poposed. Additionally, the egion of fequencies highe than the egion of inteest should be consideed and its effects on the egion of inteest obseved. hese highe fequency modes occu at lowe levels of displacement, and so ae moe likely to be linea, and theefoe aguably may be less of a concen. he poblem with such out of ange high fequency modes is that thei effects will appea in the measued data, up to aound the Nyquist fequency [], but the model constucted by the NL-RDM appoach will only cove a limited numbe of modes, likely to be smalle than the numbe of modes in the measued ange. In cuve fitting of an FRF, it is common to include an uppe esidual tem to allow cudely fo the modes above the ange analysed. Howeve, thee appeas to be no obvious way of allowing fo out of ange modes above the ange of inteest; adding a epesentative mode would seem to be athe atificial. It seems the only way is to ensue that the fequency ange ISSN: (Pint); ISSN: (Online)

7 H(m/s /N) analysed is somewhat wide than that fo which the model is eally equied; this is moe impotant fo a non-linea stuctue than fo a linea one as will be shown late. B. EFFECS OF MISSED MODES In a continuous system, thee ae a geat many modes pesent. It would be impactical to model all of them, and so the effect of missed modes should be obseved. hese missed modes may be coupled to included modes and so the effect of those missed modes should not be neglected. o obseve the effect of missed modes on a simple example, a thee degee of feedom non-linea lumped paamete system is modelled, and foce, displacement, velocity and acceleation data deived. hen, in this study case, the thid and highest mode is omitted fom the fitted model (but not its effects fom the data) and the cuve fit would be un fo both systems (3 and DOF). he FRF esult is shown in Figue 9. Results have shown some loss of enegy and inaccuacy on the FRF. Howeve, the estimated paametes fo the two identified modes wee actually the same (i.e. not in eo); this is because thee was no non-linea coupling and 3 shakes wee used so the appopiation was pefect. Eos could be pesent if the appopiation was impefect and if non-linea couplings wee pesent; any eo would then depend upon how lage the missed mode contibution was and how close in fequency Fequency Figue 9: FRF, when hid mode is missed hee is no way to eliminate these eos; the only appoach is to ensue that any modes with a lage esponse ae included and that the fequency ange chosen is adequately lage so that the effect of highe fequency modes is small. B. HARMONIC EFFECS One paticula case of out of ange modes that is of inteest is when an out of ange mode occus at a hamonic of a mode within the ange of inteest. Hamonic effects need to be consideed fo a non-linea system. he hamonic behaviou depends on the ode of non-lineaity, which means that if a system has a quadatic o cubic non-lineaity, hamonic effects may be obseved at two o thee times any natual fequency espectively. he so called hamonic modes outside the egion of inteest may stoe a significant amount of enegy, and so shouldn t be neglected. he effects of these hamonic modes ae consideed though the modelling of a system with two natual fequencies, at 7. and.36 Hz, the second mode occuing at the thid hamonic of the fist mode. he simulation pogamme is un fo cases aound and at exactly thee times the atio between natual fequencies to indicate any effects on the esults. When the system was identified with both modes pesent, the identified paametes wee not affected because all the elevant tems wee included in the analysis; howeve, consideation of this hamonic issue is equied when the 3 times natual fequency is outside the egion of inteest. he esults showed the esponse of the highe fequency (hamonic) mode is clealy much lage in this case, especially fo the acceleation. Anothe compaison based on examining time domain signals showed that the modal acceleation fo mode when seeking to excite mode was doubled when the fequencies diffeed by a facto of 3. What is of moe inteest is that when mode was analysed alone using NL-RDM, with the pesence of mode completely ignoed, the modal paametes identified wee moe in eo when the two modes diffeed by exactly a facto of 3. B.3 PRACICALIIES Pactically, it may be impotant to validate the esponse model geneated within the egion of inteest. his would be best achieved though modelling of the system with one additional mode outside this egion, and obseving the effect of the inclusion of this mode on the oiginally identified paametes and fit to a validation excitation. If thee is no significant change, then the ange chosen is aguably adequate. If thee is a significant diffeence, then the ange needs to be extended. In ecent sections has been concened with the inclusion of esidual modes, above and below any egion of inteest. he lowe esidual egion was obseved to contain igid body modes, and these wee obseved to possibly affect significantly the esults in the egion of inteest. Obtaining igid body paametes expeimentally was shown to be too difficult cuently fo the NL-RDM, given technological estictions. A mass epesentation method was geneated to model the system esponse moe accuately by adding in a modal equation fo each igid body mode. Its accuacy was demonstated though case studies. he effect of the uppe esidual aea on the modelling of the system within the egion of inteest was obseved though seveal case studies. No obvious way of allowing fo highe fequency modes seems available and so it is impotant that an adequate ange of data is consideed. Modes whose fequency is at a hamonic facto of modes of inteest need to be consideed caefully as if they ae omitted, significant eos may occu. Simple validation methods have been consideed to peseve the accuacy of the model, taking into account both hamonic and coupled esidual modes and thei effects. ISSN: (Pint); ISSN: (Online)

8 III. CONCLUSIONS In this pape, poblems with esidual effects have been investigated, including igid body mode effects in the lowe esidual aea, and missed / coupled modes in the uppe esidual aea. Hamonic effects, whee one mode natual fequency was an intege multiple of anothe mode wee also investigated. Measuing any igid body mode is difficult because, even though the stuctue is suppoted on low stiffness mounts, the fequency will usually be vey low and excitation and measuement ae inaccuate. So, some methods wee intoduced in the simulation pogamme fo and 3 degee of feedom systems to detemine o maybe neglect a igid body mode. he fist option could be to neglect the igid body mode. he esults show some inaccuacy not just aound the igid mode that neglected, but also at othe fequencies even when the igid body mode was at a significantly lowe fequency than the flexible mode(s). Anothe option could be consideed when, instead of identifying the igid body mode, the modal mass may be epesented theoetically by a summation of physical mass and inetia paametes and a theoetical modal equation added to the oveall model; a modal stiffness tem may be estimated if the igid body fequencies on the suppots can be estimated. he esults show good ageement when compaed with neglecting the igid body mode. In a eal stuctue, estimating accuate mass and inetia paametes could become a poblem although the FE model could be used. In eal stuctues, because a significant numbe of modes ae pesent in the fequency domain, only egions of inteest ae consideed but then how could uppe esidual effects be allowed fo? Effects of this egion have been investigated. Missing o including a mode in a egion of inteest is compaed fo a 3 degee of feedom system. Unlike linea identification, thee seems no way of including an uppe esidual tem in the non-linea model to allow fo the effect of the out of ange modes. o avoid this poblem, it is ecommended that exta modes be added to the egion of inteest and, the effect of the inceased size model on the validation checked to see if the diffeence is significant; so maybe moe modes need to be included than fo a linea identification. Altenatively, the level of enegy in the out of ange modes could be assessed to see if they need to be included. One moe issue, is about the so-called hamonic effect whee an out of ange mode fequency is at an intege multiple of a mode in the egion of inteest. As explained, the hamonic effect should be consideed caefully in a nonlinea system and depends on the non-linea ode (i.e. is non-lineaity of second o thid ode). Fo example fo cubic stiffness non-lineaity, 3 times fequency should be examined and, depending on the level of enegy in that mode, could be consideed o left out of the egion of inteest. he oveall study showed that the NL-RDM could allow a non-linea modal model to be identified. he levels of inaccuacy pesent fo expected levels of eo was consideed to be acceptable; even though high levels of inaccuacy wee sometimes encounteed, the eos wee highe than might be expected fom moden instumentation and methodologies. his is the natue of non-lineaity fo which thee is no pefect solution to non-linea system identification up to now. he method needs to be modified, but at the end of the day, NL-RDM is still in elatively ealy stages of exploation and development should continue. ACKNOWLEDGMEN he authos ae gateful to D Michael F. Platten fo helpful discussions on the methodology pesented in this pape. REFERENCES [] M. Samast, Identification of Non-Linea Dynamic Systems using the Non-Linea Resonant Decay Method (NL-RDM), PhD hesis, he Univesity of Mancheste, 6. [] M.F. Platten, J.R. Wight, J.E. Coope, M. Samast, Identification of multi-degee-of-feedom non-linea simulated and expeimental systems, in: Poceedings of the Intenational Semina on Modal Analysis (ISMA), Leuven,, pp. 95. [3] J.R. Wight, M.F. Platten, J.E. Coope, M. Samast, Identification of multi-degee-of-feedom weakly nonlinea systems using a model based in modal space, in: Poceedings of the Intenational Confeence on Stuctual System Identification, Kassel,, pp [4] S. Naylo, Identification of Non-popotionally Damped Stuctues using Foce Appopiation, PhD hesis, he Univesity of Mancheste, 998. [5] D. J. Ewins, Modal esting: heoy and Pactice, Reseach Studies Pess 995. [6] Mehdi Samast and Jan R. Wight, Sensitivity of NL- RDM o Eos In Paametes Used Fom he Undelying Linea System, Poceedings of the ASME 9 Intenational Mechanical Engineeing Congess & Exposition, Novembe 3-9, Lake Buena Vista, Floida, [7] R. Williams, J. Cowley and H. Vold, he Multivaiate Mode Indicato Function in Modal Analysis, Poceeding of 4 th Intenational Modal Analysis Confeence, Los Angeles 986. [8] J. R. Wight, M. F. Platten, J. E. Coope and M. Samast, Expeimental Identification of Continuous Non-linea Systems using an Extension of Foce Appopiation, st Intenational Modal Analysis Confeence, Floida, 3. [9] M. I. McEwan, J. R. Wight, J. E. Coope and A.Y.. Leung, A Finite Element echnique fo Non-linea plate and Stiffness Panel Response Pediction, AIAA Confeence Pape 595,. [] K. Woden and G. R. omlinson, Nonlineaity in Stuctual Dynamics, Institute of Physics Publishing,. ISSN: (Pint); ISSN: (Online)

APPLICATION OF MAC IN THE FREQUENCY DOMAIN

APPLICATION OF MAC IN THE FREQUENCY DOMAIN PPLICION OF MC IN HE FREQUENCY DOMIN D. Fotsch and D. J. Ewins Dynamics Section, Mechanical Engineeing Depatment Impeial College of Science, echnology and Medicine London SW7 2B, United Kingdom BSRC he

More information

A NEW VARIABLE STIFFNESS SPRING USING A PRESTRESSED MECHANISM

A NEW VARIABLE STIFFNESS SPRING USING A PRESTRESSED MECHANISM Poceedings of the ASME 2010 Intenational Design Engineeing Technical Confeences & Computes and Infomation in Engineeing Confeence IDETC/CIE 2010 August 15-18, 2010, Monteal, Quebec, Canada DETC2010-28496

More information

LINEAR AND NONLINEAR ANALYSES OF A WIND-TUNNEL BALANCE

LINEAR AND NONLINEAR ANALYSES OF A WIND-TUNNEL BALANCE LINEAR AND NONLINEAR ANALYSES O A WIND-TUNNEL INTRODUCTION BALANCE R. Kakehabadi and R. D. Rhew NASA LaRC, Hampton, VA The NASA Langley Reseach Cente (LaRC) has been designing stain-gauge balances fo utilization

More information

Numerical Integration

Numerical Integration MCEN 473/573 Chapte 0 Numeical Integation Fall, 2006 Textbook, 0.4 and 0.5 Isopaametic Fomula Numeical Integation [] e [ ] T k = h B [ D][ B] e B Jdsdt In pactice, the element stiffness is calculated numeically.

More information

Determining solar characteristics using planetary data

Determining solar characteristics using planetary data Detemining sola chaacteistics using planetay data Intoduction The Sun is a G-type main sequence sta at the cente of the Sola System aound which the planets, including ou Eath, obit. In this investigation

More information

Identification of the degradation of railway ballast under a concrete sleeper

Identification of the degradation of railway ballast under a concrete sleeper Identification of the degadation of ailway ballast unde a concete sleepe Qin Hu 1) and Heung Fai Lam ) 1), ) Depatment of Civil and Achitectual Engineeing, City Univesity of Hong Kong, Hong Kong SAR, China.

More information

PROBLEM SET #1 SOLUTIONS by Robert A. DiStasio Jr.

PROBLEM SET #1 SOLUTIONS by Robert A. DiStasio Jr. POBLM S # SOLUIONS by obet A. DiStasio J. Q. he Bon-Oppenheime appoximation is the standad way of appoximating the gound state of a molecula system. Wite down the conditions that detemine the tonic and

More information

Analysis of high speed machining center spindle dynamic unit structure performance Yuan guowei

Analysis of high speed machining center spindle dynamic unit structure performance Yuan guowei Intenational Confeence on Intelligent Systems Reseach and Mechatonics Engineeing (ISRME 0) Analysis of high speed machining cente spindle dynamic unit stuctue pefomance Yuan guowei Liaoning jidian polytechnic,dan

More information

Physics 221 Lecture 41 Nonlinear Absorption and Refraction

Physics 221 Lecture 41 Nonlinear Absorption and Refraction Physics 221 Lectue 41 Nonlinea Absoption and Refaction Refeences Meye-Aendt, pp. 97-98. Boyd, Nonlinea Optics, 1.4 Yaiv, Optical Waves in Cystals, p. 22 (Table of cystal symmeties) 1. Intoductoy Remaks.

More information

Pearson s Chi-Square Test Modifications for Comparison of Unweighted and Weighted Histograms and Two Weighted Histograms

Pearson s Chi-Square Test Modifications for Comparison of Unweighted and Weighted Histograms and Two Weighted Histograms Peason s Chi-Squae Test Modifications fo Compaison of Unweighted and Weighted Histogams and Two Weighted Histogams Univesity of Akueyi, Bogi, v/noduslód, IS-6 Akueyi, Iceland E-mail: nikolai@unak.is Two

More information

IN SITU SOUND ABSORPTION COEFFICIENT MEASUREMENT OF VARIOUS SURFACES

IN SITU SOUND ABSORPTION COEFFICIENT MEASUREMENT OF VARIOUS SURFACES IN SITU SOUND ABSORPTION COEFFICIENT MEASUREMENT OF VARIOUS SURFACES PACS REFERENCES : 43.20.El, 43.20.Ye, 43.55.Ev, 43.58.Bh Michel Béengie 1 ; Massimo Gaai 2 1 Laboatoie Cental des Ponts et Chaussées

More information

Aalborg Universitet. Load Estimation from Natural input Modal Analysis Aenlle, Manuel López; Brincker, Rune; Canteli, Alfonso Fernández

Aalborg Universitet. Load Estimation from Natural input Modal Analysis Aenlle, Manuel López; Brincker, Rune; Canteli, Alfonso Fernández Aalbog Univesitet Load Estimation fom atual input Modal Analysis Aenlle, Manuel López; Bincke, Rune; Canteli, Alfonso Fenández Published in: Confeence Poceedings Publication date: 005 Document Vesion Publishe's

More information

you of a spring. The potential energy for a spring is given by the parabola U( x)

you of a spring. The potential energy for a spring is given by the parabola U( x) Small oscillations The theoy of small oscillations is an extemely impotant topic in mechanics. Conside a system that has a potential enegy diagam as below: U B C A x Thee ae thee points of stable equilibium,

More information

Light Time Delay and Apparent Position

Light Time Delay and Apparent Position Light Time Delay and ppaent Position nalytical Gaphics, Inc. www.agi.com info@agi.com 610.981.8000 800.220.4785 Contents Intoduction... 3 Computing Light Time Delay... 3 Tansmission fom to... 4 Reception

More information

Hydroelastic Analysis of a 1900 TEU Container Ship Using Finite Element and Boundary Element Methods

Hydroelastic Analysis of a 1900 TEU Container Ship Using Finite Element and Boundary Element Methods TEAM 2007, Sept. 10-13, 2007,Yokohama, Japan Hydoelastic Analysis of a 1900 TEU Containe Ship Using Finite Element and Bounday Element Methods Ahmet Egin 1)*, Levent Kaydıhan 2) and Bahadı Uğulu 3) 1)

More information

Safety variations in steel designed using Eurocode 3

Safety variations in steel designed using Eurocode 3 JCSS Wokshop on eliability Based Code Calibation Safety vaiations in steel designed using Euocode 3 Mike Byfield Canfield Univesity Swindon, SN6 8LA, UK David Nethecot Impeial College London SW7 2BU, UK

More information

OSCILLATIONS AND GRAVITATION

OSCILLATIONS AND GRAVITATION 1. SIMPLE HARMONIC MOTION Simple hamonic motion is any motion that is equivalent to a single component of unifom cicula motion. In this situation the velocity is always geatest in the middle of the motion,

More information

Basic Bridge Circuits

Basic Bridge Circuits AN7 Datafoth Copoation Page of 6 DID YOU KNOW? Samuel Hunte Chistie (784-865) was bon in London the son of James Chistie, who founded Chistie's Fine At Auctionees. Samuel studied mathematics at Tinity

More information

School of Electrical and Computer Engineering, Cornell University. ECE 303: Electromagnetic Fields and Waves. Fall 2007

School of Electrical and Computer Engineering, Cornell University. ECE 303: Electromagnetic Fields and Waves. Fall 2007 School of Electical and Compute Engineeing, Conell Univesity ECE 303: Electomagnetic Fields and Waves Fall 007 Homewok 8 Due on Oct. 19, 007 by 5:00 PM Reading Assignments: i) Review the lectue notes.

More information

Chem 453/544 Fall /08/03. Exam #1 Solutions

Chem 453/544 Fall /08/03. Exam #1 Solutions Chem 453/544 Fall 3 /8/3 Exam # Solutions. ( points) Use the genealized compessibility diagam povided on the last page to estimate ove what ange of pessues A at oom tempeatue confoms to the ideal gas law

More information

2 Governing Equations

2 Governing Equations 2 Govening Equations This chapte develops the govening equations of motion fo a homogeneous isotopic elastic solid, using the linea thee-dimensional theoy of elasticity in cylindical coodinates. At fist,

More information

ASTR415: Problem Set #6

ASTR415: Problem Set #6 ASTR45: Poblem Set #6 Cuan D. Muhlbege Univesity of Mayland (Dated: May 7, 27) Using existing implementations of the leapfog and Runge-Kutta methods fo solving coupled odinay diffeential equations, seveal

More information

The Millikan Experiment: Determining the Elementary Charge

The Millikan Experiment: Determining the Elementary Charge LAB EXERCISE 7.5.1 7.5 The Elementay Chage (p. 374) Can you think of a method that could be used to suggest that an elementay chage exists? Figue 1 Robet Millikan (1868 1953) m + q V b The Millikan Expeiment:

More information

Right-handed screw dislocation in an isotropic solid

Right-handed screw dislocation in an isotropic solid Dislocation Mechanics Elastic Popeties of Isolated Dislocations Ou study of dislocations to this point has focused on thei geomety and thei ole in accommodating plastic defomation though thei motion. We

More information

To Feel a Force Chapter 7 Static equilibrium - torque and friction

To Feel a Force Chapter 7 Static equilibrium - torque and friction To eel a oce Chapte 7 Chapte 7: Static fiction, toque and static equilibium A. Review of foce vectos Between the eath and a small mass, gavitational foces of equal magnitude and opposite diection act on

More information

ME 3600 Control Systems Frequency Domain Analysis

ME 3600 Control Systems Frequency Domain Analysis ME 3600 Contol Systems Fequency Domain Analysis The fequency esponse of a system is defined as the steady-state esponse of the system to a sinusoidal (hamonic) input. Fo linea systems, the esulting steady-state

More information

7.2. Coulomb s Law. The Electric Force

7.2. Coulomb s Law. The Electric Force Coulomb s aw Recall that chaged objects attact some objects and epel othes at a distance, without making any contact with those objects Electic foce,, o the foce acting between two chaged objects, is somewhat

More information

Uniform Circular Motion

Uniform Circular Motion Unifom Cicula Motion Intoduction Ealie we defined acceleation as being the change in velocity with time: a = v t Until now we have only talked about changes in the magnitude of the acceleation: the speeding

More information

J. Electrical Systems 1-3 (2005): Regular paper

J. Electrical Systems 1-3 (2005): Regular paper K. Saii D. Rahem S. Saii A Miaoui Regula pape Coupled Analytical-Finite Element Methods fo Linea Electomagnetic Actuato Analysis JES Jounal of Electical Systems In this pape, a linea electomagnetic actuato

More information

Centripetal Force OBJECTIVE INTRODUCTION APPARATUS THEORY

Centripetal Force OBJECTIVE INTRODUCTION APPARATUS THEORY Centipetal Foce OBJECTIVE To veify that a mass moving in cicula motion expeiences a foce diected towad the cente of its cicula path. To detemine how the mass, velocity, and adius affect a paticle's centipetal

More information

Ultrafast effects in 3D Metamaterials

Ultrafast effects in 3D Metamaterials Ultafast effects in 3D Metamateials Nkoni Katte,2 and Philip G.Evans 3 Electical Engineeing Wilbefoce Univesity, 2 Sensos Laboatoy Wight State Univesity and 3 Oak idge National Laboatoy *Coesponding autho:

More information

Alternative Tests for the Poisson Distribution

Alternative Tests for the Poisson Distribution Chiang Mai J Sci 015; 4() : 774-78 http://epgsciencecmuacth/ejounal/ Contibuted Pape Altenative Tests fo the Poisson Distibution Manad Khamkong*[a] and Pachitjianut Siipanich [b] [a] Depatment of Statistics,

More information

The Substring Search Problem

The Substring Search Problem The Substing Seach Poblem One algoithm which is used in a vaiety of applications is the family of substing seach algoithms. These algoithms allow a use to detemine if, given two chaacte stings, one is

More information

Steady State and Transient Performance Analysis of Three Phase Induction Machine using MATLAB Simulations

Steady State and Transient Performance Analysis of Three Phase Induction Machine using MATLAB Simulations Intenational Jounal of Recent Tends in Engineeing, Vol, No., May 9 Steady State and Tansient Pefomance Analysis of Thee Phase Induction Machine using MATAB Simulations Pof. Himanshu K. Patel Assistant

More information

B. Spherical Wave Propagation

B. Spherical Wave Propagation 11/8/007 Spheical Wave Popagation notes 1/1 B. Spheical Wave Popagation Evey antenna launches a spheical wave, thus its powe density educes as a function of 1, whee is the distance fom the antenna. We

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department. Problem Set 10 Solutions. r s

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department. Problem Set 10 Solutions. r s MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Depatment Physics 8.033 Decembe 5, 003 Poblem Set 10 Solutions Poblem 1 M s y x test paticle The figue above depicts the geomety of the poblem. The position

More information

EM Boundary Value Problems

EM Boundary Value Problems EM Bounday Value Poblems 10/ 9 11/ By Ilekta chistidi & Lee, Seung-Hyun A. Geneal Desciption : Maxwell Equations & Loentz Foce We want to find the equations of motion of chaged paticles. The way to do

More information

ACOUSTIC EXCITATION OF MECHATRONIC SYSTEMS

ACOUSTIC EXCITATION OF MECHATRONIC SYSTEMS Twelfth Intenational Congess on Sound and Vibation ACOUSTIC EXCITATION OF MECHATRONIC SYSTEMS N.B. Roozen, B.T. Vehaa and R.M.G. Rijs Royal Philips Electonics, Philips Applied Technologies, P.O. Box 18,

More information

MATH 415, WEEK 3: Parameter-Dependence and Bifurcations

MATH 415, WEEK 3: Parameter-Dependence and Bifurcations MATH 415, WEEK 3: Paamete-Dependence and Bifucations 1 A Note on Paamete Dependence We should pause to make a bief note about the ole played in the study of dynamical systems by the system s paametes.

More information

THE INFLUENCE OF THE MAGNETIC NON-LINEARITY ON THE MAGNETOSTATIC SHIELDS DESIGN

THE INFLUENCE OF THE MAGNETIC NON-LINEARITY ON THE MAGNETOSTATIC SHIELDS DESIGN THE INFLUENCE OF THE MAGNETIC NON-LINEARITY ON THE MAGNETOSTATIC SHIELDS DESIGN LIVIU NEAMŢ 1, ALINA NEAMŢ, MIRCEA HORGOŞ 1 Key wods: Magnetostatic shields, Magnetic non-lineaity, Finite element method.

More information

SIMULAITION AND ANALYSIS OF PROCESS-INDUCED DISTORTIONS FOR COMPOSITE STRUCTURES

SIMULAITION AND ANALYSIS OF PROCESS-INDUCED DISTORTIONS FOR COMPOSITE STRUCTURES 21 st Intenational Confeence on Composite Mateials Xi an, 20-25 th August 2017 SIMULAITION AND ANALYSIS OF PROCESS-INDUCED DISTORTIONS FOR COMPOSITE STRUCTURES LIU Kai *1, ZHANG Boming 1, YE Jinui 1, Tang

More information

ELASTIC ANALYSIS OF CIRCULAR SANDWICH PLATES WITH FGM FACE-SHEETS

ELASTIC ANALYSIS OF CIRCULAR SANDWICH PLATES WITH FGM FACE-SHEETS THE 9 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS ELASTIC ANALYSIS OF CIRCULAR SANDWICH PLATES WITH FGM FACE-SHEETS R. Sbulati *, S. R. Atashipou Depatment of Civil, Chemical and Envionmental Engineeing,

More information

Physics 181. Assignment 4

Physics 181. Assignment 4 Physics 181 Assignment 4 Solutions 1. A sphee has within it a gavitational field given by g = g, whee g is constant and is the position vecto of the field point elative to the cente of the sphee. This

More information

BIFURCATION ANALYSIS FOR A QUASI-LINEAR CENTRIFUGAL PENDULUM ABSORBER

BIFURCATION ANALYSIS FOR A QUASI-LINEAR CENTRIFUGAL PENDULUM ABSORBER 11 th Intenational Confeence on Vibation Poblems Z. Dimitovová et al. (eds.) Lisbon, Potugal, 9-12 Septembe 213 BIFURCATION ANALYSIS FOR A QUASI-LINEAR CENTRIFUGAL PENDULUM ABSORBER Eugen B. Keme* 1, Mikhail

More information

Fresnel Diffraction. monchromatic light source

Fresnel Diffraction. monchromatic light source Fesnel Diffaction Equipment Helium-Neon lase (632.8 nm) on 2 axis tanslation stage, Concave lens (focal length 3.80 cm) mounted on slide holde, iis mounted on slide holde, m optical bench, micoscope slide

More information

Circular Orbits. and g =

Circular Orbits. and g = using analyse planetay and satellite motion modelled as unifom cicula motion in a univesal gavitation field, a = v = 4π and g = T GM1 GM and F = 1M SATELLITES IN OBIT A satellite is any object that is

More information

Chapter 9 Dynamic stability analysis III Lateral motion (Lectures 33 and 34)

Chapter 9 Dynamic stability analysis III Lateral motion (Lectures 33 and 34) Pof. E.G. Tulapukaa Stability and contol Chapte 9 Dynamic stability analysis Lateal motion (Lectues 33 and 34) Keywods : Lateal dynamic stability - state vaiable fom of equations, chaacteistic equation

More information

Pulse Neutron Neutron (PNN) tool logging for porosity Some theoretical aspects

Pulse Neutron Neutron (PNN) tool logging for porosity Some theoretical aspects Pulse Neuton Neuton (PNN) tool logging fo poosity Some theoetical aspects Intoduction Pehaps the most citicism of Pulse Neuton Neuon (PNN) logging methods has been chage that PNN is to sensitive to the

More information

MODULE 5 ADVANCED MECHANICS GRAVITATIONAL FIELD: MOTION OF PLANETS AND SATELLITES VISUAL PHYSICS ONLINE

MODULE 5 ADVANCED MECHANICS GRAVITATIONAL FIELD: MOTION OF PLANETS AND SATELLITES VISUAL PHYSICS ONLINE VISUAL PHYSICS ONLIN MODUL 5 ADVANCD MCHANICS GRAVITATIONAL FILD: MOTION OF PLANTS AND SATLLITS SATLLITS: Obital motion of object of mass m about a massive object of mass M (m

More information

Experimental Modal Analysis of Reinforced Concrete Square Slabs

Experimental Modal Analysis of Reinforced Concrete Square Slabs Expeimental Modal Analysis of Reinfoced Concete Squae Slabs M. S. Ahmed, F. A. Mohammad Abstact The aim of this pape is to pefom expeimental modal analysis (EMA) of einfoced concete (RC) squae slabs. EMA

More information

Physics 211: Newton s Second Law

Physics 211: Newton s Second Law Physics 211: Newton s Second Law Reading Assignment: Chapte 5, Sections 5-9 Chapte 6, Section 2-3 Si Isaac Newton Bon: Januay 4, 1643 Died: Mach 31, 1727 Intoduction: Kinematics is the study of how objects

More information

r cos, and y r sin with the origin of coordinate system located at

r cos, and y r sin with the origin of coordinate system located at Lectue 3-3 Kinematics of Rotation Duing ou peious lectues we hae consideed diffeent examples of motion in one and seeal dimensions. But in each case the moing object was consideed as a paticle-like object,

More information

LC transfer of energy between the driving source and the circuit will be a maximum.

LC transfer of energy between the driving source and the circuit will be a maximum. The Q of oscillatos efeences: L.. Fotney Pinciples of Electonics: Analog and Digital, Hacout Bace Jovanovich 987, Chapte (AC Cicuits) H. J. Pain The Physics of Vibations and Waves, 5 th edition, Wiley

More information

Lab #4: Newton s Second Law

Lab #4: Newton s Second Law Lab #4: Newton s Second Law Si Isaac Newton Reading Assignment: bon: Januay 4, 1643 Chapte 5 died: Mach 31, 1727 Chapte 9, Section 9-7 Intoduction: Potait of Isaac Newton by Si Godfey Knelle http://www.newton.cam.ac.uk/at/potait.html

More information

Chapter 13 Gravitation

Chapter 13 Gravitation Chapte 13 Gavitation In this chapte we will exploe the following topics: -Newton s law of gavitation, which descibes the attactive foce between two point masses and its application to extended objects

More information

Between any two masses, there exists a mutual attractive force.

Between any two masses, there exists a mutual attractive force. YEAR 12 PHYSICS: GRAVITATION PAST EXAM QUESTIONS Name: QUESTION 1 (1995 EXAM) (a) State Newton s Univesal Law of Gavitation in wods Between any two masses, thee exists a mutual attactive foce. This foce

More information

Absorption Rate into a Small Sphere for a Diffusing Particle Confined in a Large Sphere

Absorption Rate into a Small Sphere for a Diffusing Particle Confined in a Large Sphere Applied Mathematics, 06, 7, 709-70 Published Online Apil 06 in SciRes. http://www.scip.og/jounal/am http://dx.doi.og/0.46/am.06.77065 Absoption Rate into a Small Sphee fo a Diffusing Paticle Confined in

More information

Lab 10: Newton s Second Law in Rotation

Lab 10: Newton s Second Law in Rotation Lab 10: Newton s Second Law in Rotation We can descibe the motion of objects that otate (i.e. spin on an axis, like a popelle o a doo) using the same definitions, adapted fo otational motion, that we have

More information

Duality between Statical and Kinematical Engineering Systems

Duality between Statical and Kinematical Engineering Systems Pape 00, Civil-Comp Ltd., Stiling, Scotland Poceedings of the Sixth Intenational Confeence on Computational Stuctues Technology, B.H.V. Topping and Z. Bittna (Editos), Civil-Comp Pess, Stiling, Scotland.

More information

FUNDAMENTAL PROPERTIES OF LINEAR SHIP STEERING DYNAMIC MODELS

FUNDAMENTAL PROPERTIES OF LINEAR SHIP STEERING DYNAMIC MODELS Jounal of Maine Science and Technology, Vol. 7, No. 2, pp. 79-88 (1999) 79 FUNDAMENTAL PROPERTIES OF LINEAR SHIP STEERING DYNAMIC MODELS Ching-Yaw Tzeng* and Ju-Fen Chen** Keywods: Nomoto model, Contollability

More information

A Three-Dimensional Magnetic Force Solution Between Axially-Polarized Permanent-Magnet Cylinders for Different Magnetic Arrangements

A Three-Dimensional Magnetic Force Solution Between Axially-Polarized Permanent-Magnet Cylinders for Different Magnetic Arrangements Poceedings of the 213 Intenational Confeence on echanics, Fluids, Heat, Elasticity Electomagnetic Fields A Thee-Dimensional agnetic Foce Solution Between Axially-Polaied Pemanent-agnet Cylindes fo Diffeent

More information

Physics 235 Chapter 5. Chapter 5 Gravitation

Physics 235 Chapter 5. Chapter 5 Gravitation Chapte 5 Gavitation In this Chapte we will eview the popeties of the gavitational foce. The gavitational foce has been discussed in geat detail in you intoductoy physics couses, and we will pimaily focus

More information

Graphs of Sine and Cosine Functions

Graphs of Sine and Cosine Functions Gaphs of Sine and Cosine Functions In pevious sections, we defined the tigonometic o cicula functions in tems of the movement of a point aound the cicumfeence of a unit cicle, o the angle fomed by the

More information

Effect of drag on the performance for an efficient wind turbine blade design

Effect of drag on the performance for an efficient wind turbine blade design Available online at www.sciencediect.com Enegy Pocedia 18 (01 ) 404 415 Abstact Effect of dag on the pefomance fo an efficient wind tubine blade design D. Eng. Ali H. Almukhta Univesity of Technology Email-

More information

Electrodynamic Forces in Steel Strip during Induction Heating

Electrodynamic Forces in Steel Strip during Induction Heating Intenational Scientific Colloquium Modelling fo Electomagnetic Pocessing Hannove, Mach 4-6, 3 Electodynamic Foces in Steel Stip duing Induction Heating H. Kasjanow, H. Schülbe,. Nikanoov bstact Tangential

More information

Computational Methods of Solid Mechanics. Project report

Computational Methods of Solid Mechanics. Project report Computational Methods of Solid Mechanics Poject epot Due on Dec. 6, 25 Pof. Allan F. Bowe Weilin Deng Simulation of adhesive contact with molecula potential Poject desciption In the poject, we will investigate

More information

Experiment I Voltage Variation and Control

Experiment I Voltage Variation and Control ELE303 Electicity Netwoks Expeiment I oltage aiation and ontol Objective To demonstate that the voltage diffeence between the sending end of a tansmission line and the load o eceiving end depends mainly

More information

15 B1 1. Figure 1. At what speed would the car have to travel for resonant oscillations to occur? Comment on your answer.

15 B1 1. Figure 1. At what speed would the car have to travel for resonant oscillations to occur? Comment on your answer. Kiangsu-Chekiang College (Shatin) F:EasteHolidaysAssignmentAns.doc Easte Holidays Assignment Answe Fom 6B Subject: Physics. (a) State the conditions fo a body to undego simple hamonic motion. ( mak) (a)

More information

Analysis and Optimization of a Special Type of Dielectric Loaded Resonant Cavity for Mobile Communication Filters

Analysis and Optimization of a Special Type of Dielectric Loaded Resonant Cavity for Mobile Communication Filters 328 Analysis and Optimization of a Special Type of Dielectic Loaded Resonant Cavity fo Mobile Communication Filtes Haold S. Showes, Banmali S. Rawat *, Syam S. Challa Depatment of Electical and Biomedical

More information

Central Coverage Bayes Prediction Intervals for the Generalized Pareto Distribution

Central Coverage Bayes Prediction Intervals for the Generalized Pareto Distribution Statistics Reseach Lettes Vol. Iss., Novembe Cental Coveage Bayes Pediction Intevals fo the Genealized Paeto Distibution Gyan Pakash Depatment of Community Medicine S. N. Medical College, Aga, U. P., India

More information

2. Electrostatics. Dr. Rakhesh Singh Kshetrimayum 8/11/ Electromagnetic Field Theory by R. S. Kshetrimayum

2. Electrostatics. Dr. Rakhesh Singh Kshetrimayum 8/11/ Electromagnetic Field Theory by R. S. Kshetrimayum 2. Electostatics D. Rakhesh Singh Kshetimayum 1 2.1 Intoduction In this chapte, we will study how to find the electostatic fields fo vaious cases? fo symmetic known chage distibution fo un-symmetic known

More information

is the instantaneous position vector of any grid point or fluid

is the instantaneous position vector of any grid point or fluid Absolute inetial, elative inetial and non-inetial coodinates fo a moving but non-defoming contol volume Tao Xing, Pablo Caica, and Fed Sten bjective Deive and coelate the govening equations of motion in

More information

rt () is constant. We know how to find the length of the radius vector by r( t) r( t) r( t)

rt () is constant. We know how to find the length of the radius vector by r( t) r( t) r( t) Cicula Motion Fom ancient times cicula tajectoies hae occupied a special place in ou model of the Uniese. Although these obits hae been eplaced by the moe geneal elliptical geomety, cicula motion is still

More information

The geometric construction of Ewald sphere and Bragg condition:

The geometric construction of Ewald sphere and Bragg condition: The geometic constuction of Ewald sphee and Bagg condition: The constuction of Ewald sphee must be done such that the Bagg condition is satisfied. This can be done as follows: i) Daw a wave vecto k in

More information

Physics 2A Chapter 10 - Moment of Inertia Fall 2018

Physics 2A Chapter 10 - Moment of Inertia Fall 2018 Physics Chapte 0 - oment of netia Fall 08 The moment of inetia of a otating object is a measue of its otational inetia in the same way that the mass of an object is a measue of its inetia fo linea motion.

More information

Rigid Body Dynamics 2. CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2018

Rigid Body Dynamics 2. CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2018 Rigid Body Dynamics 2 CSE169: Compute Animation nstucto: Steve Rotenbeg UCSD, Winte 2018 Coss Poduct & Hat Opeato Deivative of a Rotating Vecto Let s say that vecto is otating aound the oigin, maintaining

More information

Teachers notes. Beyond the Thrills excursions. Worksheets in this book. Completing the worksheets

Teachers notes. Beyond the Thrills excursions. Worksheets in this book. Completing the worksheets Beyond the Thills excusions Teaches notes Physics is the science of how the wold (and Univese) woks. Luna Pak Sydney is a lage hands-on physics laboatoy full of fee falling objects, otating systems and

More information

Conservative Averaging Method and its Application for One Heat Conduction Problem

Conservative Averaging Method and its Application for One Heat Conduction Problem Poceedings of the 4th WSEAS Int. Conf. on HEAT TRANSFER THERMAL ENGINEERING and ENVIRONMENT Elounda Geece August - 6 (pp6-) Consevative Aveaging Method and its Application fo One Heat Conduction Poblem

More information

Rotor Blade Performance Analysis with Blade Element Momentum Theory

Rotor Blade Performance Analysis with Blade Element Momentum Theory Available online at www.sciencediect.com ScienceDiect Enegy Pocedia 5 (7 ) 3 9 The 8 th Intenational Confeence on Applied Enegy ICAE6 Roto Blade Pefomance Analysis with Blade Element Momentum Theoy Faisal

More information

12th WSEAS Int. Conf. on APPLIED MATHEMATICS, Cairo, Egypt, December 29-31,

12th WSEAS Int. Conf. on APPLIED MATHEMATICS, Cairo, Egypt, December 29-31, th WSEAS Int. Conf. on APPLIED MATHEMATICS, Caio, Egypt, Decembe 9-3, 7 5 Magnetostatic Field calculations associated with thick Solenoids in the Pesence of Ion using a Powe Seies expansion and the Complete

More information

Chapter 5 Force and Motion

Chapter 5 Force and Motion Chapte 5 Foce and Motion In Chaptes 2 and 4 we have studied kinematics, i.e., we descibed the motion of objects using paametes such as the position vecto, velocity, and acceleation without any insights

More information

Chapter 5 Linear Equations: Basic Theory and Practice

Chapter 5 Linear Equations: Basic Theory and Practice Chapte 5 inea Equations: Basic Theoy and actice In this chapte and the next, we ae inteested in the linea algebaic equation AX = b, (5-1) whee A is an m n matix, X is an n 1 vecto to be solved fo, and

More information

Chapter 5 Force and Motion

Chapter 5 Force and Motion Chapte 5 Foce and Motion In chaptes 2 and 4 we have studied kinematics i.e. descibed the motion of objects using paametes such as the position vecto, velocity and acceleation without any insights as to

More information

Chapter 3 Optical Systems with Annular Pupils

Chapter 3 Optical Systems with Annular Pupils Chapte 3 Optical Systems with Annula Pupils 3 INTRODUCTION In this chapte, we discuss the imaging popeties of a system with an annula pupil in a manne simila to those fo a system with a cicula pupil The

More information

Application of homotopy perturbation method to the Navier-Stokes equations in cylindrical coordinates

Application of homotopy perturbation method to the Navier-Stokes equations in cylindrical coordinates Computational Ecology and Softwae 5 5(): 9-5 Aticle Application of homotopy petubation method to the Navie-Stokes equations in cylindical coodinates H. A. Wahab Anwa Jamal Saia Bhatti Muhammad Naeem Muhammad

More information

Modeling Fermi Level Effects in Atomistic Simulations

Modeling Fermi Level Effects in Atomistic Simulations Mat. Res. Soc. Symp. Poc. Vol. 717 Mateials Reseach Society Modeling Femi Level Effects in Atomistic Simulations Zudian Qin and Scott T. Dunham Depatment of Electical Engineeing, Univesity of Washington,

More information

DEVIL PHYSICS THE BADDEST CLASS ON CAMPUS IB PHYSICS

DEVIL PHYSICS THE BADDEST CLASS ON CAMPUS IB PHYSICS DEVIL PHYSICS THE BADDEST CLASS ON CAMPUS IB PHYSICS TSOKOS LESSON 10-1 DESCRIBING FIELDS Essential Idea: Electic chages and masses each influence the space aound them and that influence can be epesented

More information

A scaling-up methodology for co-rotating twin-screw extruders

A scaling-up methodology for co-rotating twin-screw extruders A scaling-up methodology fo co-otating twin-scew extudes A. Gaspa-Cunha, J. A. Covas Institute fo Polymes and Composites/I3N, Univesity of Minho, Guimaães 4800-058, Potugal Abstact. Scaling-up of co-otating

More information

( ) [ ] [ ] [ ] δf φ = F φ+δφ F. xdx.

( ) [ ] [ ] [ ] δf φ = F φ+δφ F. xdx. 9. LAGRANGIAN OF THE ELECTROMAGNETIC FIELD In the pevious section the Lagangian and Hamiltonian of an ensemble of point paticles was developed. This appoach is based on a qt. This discete fomulation can

More information

COUPLED MODELS OF ROLLING, SLIDING AND WHIRLING FRICTION

COUPLED MODELS OF ROLLING, SLIDING AND WHIRLING FRICTION ENOC 008 Saint Petesbug Russia June 30-July 4 008 COUPLED MODELS OF ROLLING SLIDING AND WHIRLING FRICTION Alexey Kieenkov Ins ti tu te fo P ob le ms in Me ch an ic s Ru ss ia n Ac ad em y of Sc ie nc es

More information

= 4 3 π( m) 3 (5480 kg m 3 ) = kg.

= 4 3 π( m) 3 (5480 kg m 3 ) = kg. CHAPTER 11 THE GRAVITATIONAL FIELD Newton s Law of Gavitation m 1 m A foce of attaction occus between two masses given by Newton s Law of Gavitation Inetial mass and gavitational mass Gavitational potential

More information

This is a very simple sampling mode, and this article propose an algorithm about how to recover x from y in this condition.

This is a very simple sampling mode, and this article propose an algorithm about how to recover x from y in this condition. 3d Intenational Confeence on Multimedia echnology(icm 03) A Simple Compessive Sampling Mode and the Recovey of Natue Images Based on Pixel Value Substitution Wenping Shao, Lin Ni Abstact: Compessive Sampling

More information

International Journal of Scientific & Engineering Research, Volume 7, Issue 4, April ISSN

International Journal of Scientific & Engineering Research, Volume 7, Issue 4, April ISSN Intenational Jounal of Scientific & Engineeing Reseach, Volume 7, Issue 4, Apil-06 67 ISSN 9-558 Reduction at Payload Using Isolation System in a ypical Launch Vehicle Lekshmi Piya S,D. Sundaaajan, D.K

More information

Encapsulation theory: the transformation equations of absolute information hiding.

Encapsulation theory: the transformation equations of absolute information hiding. 1 Encapsulation theoy: the tansfomation equations of absolute infomation hiding. Edmund Kiwan * www.edmundkiwan.com Abstact This pape descibes how the potential coupling of a set vaies as the set is tansfomed,

More information

A new approach in classical electrodynamics to protect principle of causality

A new approach in classical electrodynamics to protect principle of causality A new appoach in classical electodynamics to potect pinciple of causality Biswaanjan Dikshit * Lase and Plasma Technology Division Bhabha Atomic Reseach Cente, Mumbai-400085 INDIA * Coesponding autho E-mail:

More information

AST 121S: The origin and evolution of the Universe. Introduction to Mathematical Handout 1

AST 121S: The origin and evolution of the Universe. Introduction to Mathematical Handout 1 Please ead this fist... AST S: The oigin and evolution of the Univese Intoduction to Mathematical Handout This is an unusually long hand-out and one which uses in places mathematics that you may not be

More information

Chapter 7-8 Rotational Motion

Chapter 7-8 Rotational Motion Chapte 7-8 Rotational Motion What is a Rigid Body? Rotational Kinematics Angula Velocity ω and Acceleation α Unifom Rotational Motion: Kinematics Unifom Cicula Motion: Kinematics and Dynamics The Toque,

More information

ON THE TWO-BODY PROBLEM IN QUANTUM MECHANICS

ON THE TWO-BODY PROBLEM IN QUANTUM MECHANICS ON THE TWO-BODY PROBLEM IN QUANTUM MECHANICS L. MICU Hoia Hulubei National Institute fo Physics and Nuclea Engineeing, P.O. Box MG-6, RO-0775 Buchaest-Maguele, Romania, E-mail: lmicu@theoy.nipne.o (Received

More information

Swissmetro: design methods for ironless linear transformer

Swissmetro: design methods for ironless linear transformer Swissmeto: design methods fo ionless linea tansfome Nicolas Macabey GESTE Engineeing SA Scientific Pak PSE-C, CH-05 Lausanne, Switzeland Tel (+4) 2 693 83 60, Fax. (+4) 2 693 83 6, nicolas.macabey@geste.ch

More information

Chapter Sixteen: Electric Charge and Electric Fields

Chapter Sixteen: Electric Charge and Electric Fields Chapte Sixteen: Electic Chage and Electic Fields Key Tems Chage Conducto The fundamental electical popety to which the mutual attactions o epulsions between electons and potons ae attibuted. Any mateial

More information