On the Development of a Walking Rehabilitation Device with a Large Workspace

Size: px
Start display at page:

Download "On the Development of a Walking Rehabilitation Device with a Large Workspace"

Transcription

1 2011 IEEE Internationa Conference on Rehabiitation Robotics Rehab Week Zurich, ETH Zurich Science City, Switzerand, June 29 - Juy 1, 2011 On the Deveopment of a Waking Rehabiitation Device with a Large Workspace Cément Gossein and Thierry Laiberté Abstract Rehabiitation therapy aiming at heping a person to regain or improve the abiity to wak is a abour-intensive activity. In this context, the patient is often imited to very restricted waking spaces. This paper presents a device that can support the weight of a person waking freey over a arge workspace. The device is based on an overhead gantry mechanism combined with a cabe routing that decoupes the vertica motion from the horizonta dispacements of the person. The mechanica architecture is first presented and it is shown that the principe can be appied to the design of a competey passive device in which the portion of the weight to be supported can be adjusted. A simpified dynamic mode is aso derived in order to highight the characteristics of the device. A powered version of the device is then discussed. Finay, a prototype of a passive device buit at fu scae is presented and discussed. A video accompanying the paper iustrates the experimenta tests underway with the prototype. I. INTRODUCTION Bipeda waking is a mechanicay compex activity that invoves the coordination of many musces based on a variety of sensing inputs. Therefore, the rehabiitation of patients with waking disabiities requires earning, which can in turn transate into ong hours of exercises. Moreover, assisting a person that requires support to wak is a difficut task. If the person is assisted ony by therapists, then severa of them are needed which becomes very costy. On the other hand, if the person is supported by a harness, the waking space is usuay very imited (e.g. a treadmi is used) which is a significant drawback, especiay in advanced stages of the earning process. The deveopment of devices to assist in the rehabiitation of patients with waking disabiities has been addressed by severa research groups. Two main issues have been considered, namey the support of the patient and the guidance of his/her imbs. Most of the proposed concepts are based on a the use of a treadmi (see [1] and [2] for instance). Athough treadmis offer some advantages in terms of footprint and contro, they cannot reproduce the rea dynamics of waking because the inertia of the upper body of the patient is not incuded in the waking scenario. Other simiar devices have aso been deveoped for appications in the support of edery peope (see for instance [3]). In [4], a robotic assistive device is used to support the weight of the patient. The portion of the weight to be supported can be adjusted. However, the workspace avaiabe This work was supported by The Natura Sciences and Engineering Research Counci of Canada (NSERC) as we as the Canada Research Chair program. The authors are with the Department of Mechanica Engineering, Université Lava, Québec, Québec, Canada, G1V 0A6, gossein@gmc.uava.ca, thierry@gmc.uava.ca /11/$ IEEE 263 to the patient is imited and the foor footprint of the device is rather arge. Mobie robotic devices that foow the patient have aso been deveoped (see for instance [5] and [6]) in order to increase the workspace and aow the patient to wak naturay on a foor. However, this approach invoves a arge mobie robot that requires significant foor space and raises issues regarding the safe active contro of the robot. In other approaches, robotic patforms are attached to the feet of the user and controed directy (see for instance [7], [8] and [9]). Such systems are simiar to the use of treadmis regarding the waking simuation, athough they provide more fexibiity and aow the simuation of obstaces, stairs and others. In [10], a waking rehabiitation device based on a cabe-driven parae mechanism is proposed. This concept eads to an effective ight-weight apparatus and aows the appication of controed forces to the person. However, it aso requires the use of a treadmi with the patient waking on the spot. In [11], a robotic system suspended from a ceiing-mounted rai is proposed. This device aows the patient to wak more naturay on a foor. However, deviations from a inear path and changes of direction are not possibe. Other commercia systems based on some of the concepts mentioned above are aso avaiabe (see for instance [12] and [13]). This paper presents a new device that aims at providing a arge waking space whie supporting the weight of the patient. The proposed concept aows the user to wak freey on a foor in a directions, incuding changes of rotation, over a arge workspace. The device naturay eads itsef to a passive mechanism but can aso be impemented as an active robotic system. After presenting the basic principe of the device, its impementation as a passive mechanism or as an active robotic system is discussed. A dynamic mode is then studied to highight the design chaenges. Finay, the experimenta vaidation on a fu-scae prototype is presented and some concusions are drawn. II. MECHANICAL ARCHITECTURE A. Four-degree-of-freedom gantry mechanism The device proposed in this work is based on an overhead gantry mechanism simiar to those commony found in industria settings. Two parae rais attached to the ceiing provide the motion trave in one horizonta direction. A bridge is suspended to these rais. A troey is suspended to the bridge on a rai orthogona to the first set of rais. The troey provides the motion trave in the second horizonta direction. A vertica sider mechanism is attached to the troey and provides the third degree of freedom, i.e., vertica

2 Y Z X Bridge Troey Y Z X Bridge x θ y Vertica Actuation and Baancing Troey End effector z Payoad Fig. 1. Schematic representation of the 4-dof gantry structure. Fig. 2. Cabe routing for the weight compensation mechanism. transation. Finay, a revoute joint with a vertica axis can be added in order to aow the end body to rotate freey around a vertica axis. Aternativey, the vertica sider and the rotation can be repaced by a cabe hoist. The four degrees of freedom of the mechanism correspond to the so-caed Schönfies motions or SCARA motions. The kinematic architecture of the gantry mechanism is iustrated schematicay in Fig. 1. One advantage of an overhead gantry mechanism is that it does not occupy any of the foor space. B. Weight compensation The overhead gantry mechanism described in the preceding subsection provides the degrees of freedom required to aow a person to wak freey in a directions, incuding a rotation around a vertica axis. However, the vertica degree of freedom must provide a force to baance the weight (or part of the weight) of the person. In other words, the main objective of this mechanism is to be abe to foow the person with an overhead support and compensate part of its weight. One effective means of supporting the weight of the person is to use an adjustabe counterweight mounted on the vertica sider mechanism. However, mounting the counterweight directy on the vertica sider woud significanty increase the inertia of the moving bodies of the mechanism, thereby eading to a arge inertia in the horizonta directions. In a passive system, this inertia woud be highy detrimenta because it woud require the person to drag a significant additiona mass. Therefore, it is desired to move the counterweight away from the person and mount it on the fixed frame. A cabe transmission using the routing iustrated in Fig. 2 can be used to move the counterweight to a fixed ocation. This arrangement basicay connects the counterweight to the person via a cabe transmission that cances out the horizonta motions. Additionay, it is noted that the tota ength of cabe in the routing is constant and therefore no spoo is required. The motion of the counterweight is competey decouped from the horizonta motions which means that horizonta motions of the mechanism eave the counterweight at rest. A motion of the counterweight is induced ony when a 20 kg Fig kg 80 kg Adjustabe weight compensation mechanism. 100 kg vertica motion of the end-effector takes pace. Therefore, the counterweight does not contribute to the horizonta inertia. The weight compensation mechanism must be adjustabe in order to accommodate variations on the weight of the user. Moreover, for a given user, the portion of the weight to be compensated may vary depending on the person s condition. It may aso be desirabe to change this parameter during a working session in order to test different waking conditions or vary the stimui. Therefore, the counterweight is mounted on a rotating ever. Siding the counterweight aong the ever and ocking it in different positions provides a simpe and effective means of adjusting the weight compensation. The adjustabe compensation mechanism is iustrated in Fig. 3. It is cear from the above that the proposed device is wesuited for a passive assistive waking mechanism. Indeed, horizonta motions are provided by the bridge and troey gantry mechanism which can be designed to induce imited friction and does not require arge forces to be activated. Aso, because of the cabe transmission, the counterweight is mounted on a fixed pivoting ever and does not contribute to increasing the inertia in the horizonta directions. Finay, the adjustment of the mass to be compensated is easiy performed by moving the counterweight aong the ever arm. III. DYNAMIC MODELLING A dynamic mode is derived in order to determine the forces required to acceerate the mechanism when performing waking exercises. One can write f = f r + A(p) p + b(p, ṗ) + c(p) (1) where f is the externa force vector defined as f = [f x f y f z τ] T, (2) 264

3 f r is the friction force vector defined as f r = [f rx f ry f rz τ r ] T, (3) p is the Cartesian position vector defined as p = [x y z θ] T, (4) in which x, y, z and θ are the position components in a directions, where θ is the rotation around the vertica axis, ṗ is the Cartesian veocity vector and p is the Cartesian acceeration vector defined as [ ṗ = ẋ ẏ ż θ ] T [ θ] T, p = ẍ ÿ z. (5) Matrix A describes the inertia of the system which can be written as M A = 0 M M 33 0 (6) I z with M 11 = m 1 + m 2 + m 3 + m h (7) M 22 = m 2 + m 3 + m h (8) and where m 1, m 2, m 3 and m h are respectivey the mass of the bridge, the mass of the troey, the mass of the endeffector and the mass of the human subject whie I z is the inertia of the end-effector pus the human subject around the vertica axis. Vector b represents the so-caed centrifuga and Coriois terms b = [0 0 B 3 0] T, (9) whose components are a zero except for the vertica motion for which there is a noninear reationship between the motion of the end-effector and the counterweight. Vector c represents the terms arising from the potentia energy of the system and can be written as c = [0 0 C 3 0] T, (10) whose components are a zero except for the vertica motion, which invoves gravitationa potentia energy. The expressions for M 33, B 3 and C 3 wi be given after the kinematics of the baancing system is estabished. The geometry of the baancing system is described in Fig. 4. It can be seen that with From the reationships z = ρ o ρ (11) ρ o = h 2 + ( o ) 2. (12) ρ sin φ = o cos ψ (13) h ρ cos φ = sin ψ (14) and eiminating φ, one can obtain ρ as ρ = h o + 2 2(h sin ψ + o cos ψ). (15) 265 Fig. 4. m + m 3 h z φ h ρ o ρ m Geometry of the baancing system. The time derivative of eq. (11) gives c o ż = (h cos ψ o sin ψ) ψ. (16) ρ In order to simpify the dynamic modeing and provide more insight, the foowing approximations are introduced. The engths o and are assumed to be approximatey equa, h is assumed to be significanty arger than o and, and ψ is kept reativey sma, say ess than 30 degrees. As a resut, ange φ is sma and cos φ 1. Then, eq.(11) is approximated by z = sin ψ (17) and eq.(16) is approximated by ψ ż = ψ cos ψ (18) The dynamic mode is obtained using the Lagrangian approach ( ) d (T V ) (T V ) = f z f rz (19) dt ż z The kinetic energy T is T = 1 2 (m 3 + m h )ż I c ψ 2 (20) where I c is the inertia of the counterweight around its pivot. The potentia energy V is V = m c gc sin ψ (m 3 + m h )gz = ( m cc c m 3 m h )gz (21) where m c is the mass of the counterweight. The resut can then be incuded in eq. (1) which gives M 33 = m 3 + m h + I c 2 cos 2 ψ (22) B 3 = ż2 I c tan ψ 3 cos 3 (23) ψ ( mc c ) C 3 = m 3 m h g. (24) In the mechanica architecture described in the preceding section, the troey is mounted on the bridge and the vertica

4 sider is mounted on the troey. Therefore one has, by inspection of eqs.(7) and (8), M 11 > M 22 (25) Hence, considering eqs.(6) and (25), it is cear that the inertia is not equa in the two horizonta directions. This is a drawback of the bridge and troey arrangement. In a passive system, it means that the spurious inertia fet by the user wi not be the same depending on the direction of waking. Therefore, it is important that the design of the mechanism aims at minimizing m 1. Ceary, a masses shoud be minimized but the atter has the most impact on the performance of a passive device. Simuations can be performed based on the simpified dynamic mode presented above in order to estimate the performances in terms of the horizonta forces required for waking. However, at ow veocities, friction may be dominant and shoud therefore be estimated in a practica design. Approximating the counterweight as a point mass, its inertia around the fixed pivot, noted I c in the above equations, can be written as I c = m c c 2. (26) The baancing ratio, noted γ, is defined as the ratio of the static moment of the counterweight with respect to the fixed pivot over the static moment generated by the payoad on the pivot. In other words, when γ is equa to 1, the payoad is exacty baanced. From eq.(24), one obtains m c c γ = (m 3 + m h ). (27) The inertia ratio, noted α, is defined as the ratio of the effective inertia of the counterweight over the inertia of the payoad. A sma vaue of α is desirabe in order to minimize the spurious inertia. From eqs.(22) and (26), one obtains α = m c c 2 2 cos 2 ψ(m 3 + m h ) = cγ cos 2 ψ. (28) It can be observed that by using a arger counterweight mounted coser to the pivot, the tota inertia associated with the vertica motion is reduced. This is so because the mass of the counterweight required to obtain static baancing increases ineary with the ever arm whie the inertia increases quadraticay. Therefore, in practice a arge counterweight mounted cose to the pivot is preferabe and wi ead to smaer spurious inertia (smaer vaues of α). It can aso be observed that it is preferabe to keep the ever around the horizonta configuration in order to avoid excessive inertia. This is obtained if the ever is sufficienty ong compared to the tota vertica dispacement. Aso, from eq.(23), it can be seen that the effect of I c is imited if ψ and ż are kept sma. IV. POWERED ASSISTIVE DEVICE Athough the device described above eads naturay to a passive mechanism, it is aso possibe to introduce actuation in order to improve the performances and reduce the spurious forces to be provided by the user due to the inertia and friction in the mechanism. To this end, actuators coud be 266 incuded directy on the joints of the gantry mechanism. However, this approach woud further increase the inertia of the moving components. In order to circumvent this probem, cosed-oop bet-driven transmissions are used. Severa cosed-oop transmissions were proposed in [14] for appication in four-dof gantry mechanisms. The transmissions used here are based on the arrangements proposed in the atter reference. First, it is observed that because of the mechanica architecture of the grantry device, the inertia induced by the mechanism around the vertica axis of rotation is rather negigeabe. Therefore, this degree of freedom is eft passive and no actuation is required. Then, it is pointed out that the cabe routing used to attach the payoad to the counterweight can be used to ocate the actuator providing the vertica motion on the base. Indeed, the actuator for the vertica motion can be mounted directy at the counterweight. Finay, the horizonta motions (bridge and troey) are considered. Actuation is most reevant for these degrees of freedom since they are more ikey to produce undesirabe effects on the user. The actuation of the bridge is based on the cosed-oop transmission iustrated in Fig. 5-eft. In this figure, the bridge is represented by a rectange and the pueys by circes. The s mark the paces where the bet is attached to the bridge and the square around one of the circes denotes that the corresponding puey is actuated. As it can be observed on the figure, the ength of the cosed oop is constant and therefore no spoo is needed. Aso, the forces appied on the bridge are we baanced and provide stabiity. The actuation of the troey is based on the cosedoop transmission iustrated in Fig. 5-right. Simiary to the preceding routing, the cosed oop for the actuation of the troey aso has a constant ength and does not require a spoo. Moreover, it is pointed out that the two transmissions are competey decouped. In other words, a motion of one of the actuators does not have any impact on the other degree of freedom and vice-versa. One important advantage of decouping is that no antagonistic forces are generated, which makes the contro easier and which reduces the power required at the actuators. Using the above actuation schemes, it is possibe to actuate the horizonta motions of the gantry mechanism and compensate for the inertia and friction in these directions by tracking the motion of the torso of the person. V. EXPERIMENTAL VALIDATION A passive 3-dof version of the concept described above was impemented experimentay for demonstration purposes. The mechanism buit in the aboratory shown in Fig. 6 is based on a device that was deveoped for other appications (payoad manipuation). The main modification is the addition of a chair as an interface between the human subject and the end-effector. Therefore, athough it provides a proof of concept, the experimenta set-up is not taiored to the appication described here. First, the mechanism has ony three degrees of freedom and does not aow rotations around the vertica axis, which is definitey a imitation

5 θ1 θ 2 Fig. 5. Cosed-oop bet-driven transmissions. Left: for the actuation of the bridge motion. Right: for the actuation of the troey motion. for waking devices. Additionay, the inertia of the bridge, troey and end-effector coud be significanty reduced since the prototype mechanism was deveoped for arger payoads. Despite the aforementioned imitations, the experimenta impementation performs rather we, as shown in the accompanying video. In particuar, it can be observed that the gravity compensation is effective and we decouped from horizonta motions. It is observed that for arge baancing ratios (γ > 0.8), the static forces between the human subject and the foor are so sma that it is difficut to perform horizonta acceerations by waking. The user tends to jump when trying to appy more force on the foor. It is aso observed that for sma baancing ratios (γ < 0.2) the prototype coud not aways foow the user, i.e., the chair separates from the subject when arge upwards acceerations are performed. In order to avoid this phenomenon, the unoaded device (without the person) shoud be abe to perform arge upward acceerations. In other words, with m h = 0, C 3 shoud be as arge as possibe compared to M 33. The use of a chair aowed to adapt a system intended for another appication. However, the waking motion was not fuy natura and comfort was imited. It is cear that a better interface, such as a harness, shoud be used in the future. Horizonta motions are aso we demonstrated, athough it is cear that friction and inertia shoud be reduced in a rea impementation of the concept. As mentioned above, such a reduction shoud be possibe in a design based on reaistic requirements. In order to verify the dynamic mode aong the vertica motion, position mesasurements were performed. Numerica differentiation of the position was then performed in order to obtain the veocity and acceeration. The resuts obtained can then be compared with the acceeration computed from eq.(1). The parameters of the prototype are as foows: = 1.0m, o = 0.96m, h = 2.2m, m 3 = 12kg and m h = 70kg. For reasons reated to other appications of the prototype, the range of motion of the ever is 15 o < ψ < 30 o. As a resut, a person in a standing position with their feet on the foor corresponds to z = 0.31m or ψ = 18 o. Aso, the vaue of φ is comprised between 1 o and 1 o for 80% of the range of motion, with a maximum of 3 o. Therefore, the geometric 267 Fig. 6. Passive 3-dof assistive waking device buit in the aboratory. approximation used for the reation between z and ψ is appicabe. From measurements taken with a dynamometer, the friction is estimated to be approximatey f rz = 20N. Six baancing cases are studied. They are presented in Tabe I. Typica experimenta curves are potted in Fig. 7. Since the forces appied by the person on the foor are not known, the mode can be verified ony during the free-faing trajectory of the person jumping in the air. For the ast two cases given in the tabe, since the baancing ratio is sma, it is possibe to manuay throw the chair downwards and measure the acceeration of the unoaded device. However, it is difficut for the person to stay in the air ong enough to obtain effective measurements with the setup used. TABLE I THE SIX BALANCING CASES. THE ACCELERATIONS z ARE CALCULATED AROUND THE POINT WHERE THE VELOCITY IS ZERO. THE CORRESPONDING MEASURED VALUES ARE GIVEN IN PARENTHESES. Case m c c γ α z z (ψ = 0) with person chair ony (kg) (m) (m/s 2 ) (m/s 2 ) C (1.8) -5.6 C (2.0) -9.4 C (4.0) -4.2 C (4.4) -6.6 C (-2.2) C (-3.0)

6 Position (m) and Veocity (m/s) Veocity Position Time (s) 3 2 Acceeration (m\s2) Time (s) appication is ceary iustrated in the accompanying video, where the effect of gravity compensation is demonstrated using eaping motions that are characteristic of astronauts evoving in a ow-gravity environment. VII. CONCLUSION This paper presented a waking rehabiitation device that aows a person to wak freey on a foor over a arge workspace. The basic principe of the device consists of a cabe-routing system that foows the user in order to provide gravity compensation without hindering waking motions. The device naturay eads itsef to a passive mechanism but can aso be impemented as an active robotic system. Using a dynamic mode, it was shown that a rea impementation can be deveoped based on reaistic requirements. Aso, it was shown that the baancing principe based on the use of a counterweight can be expoited to minimize the spurious inertia. Finay, the experimenta vaidation on a fu-scae prototype provided design hints. Athough there were some imitations on the experimenta impementation, the prototype and the accompanying video ceary demonstrate the potentia of the concept in a context of rehabiitation. Fig. 7. Exampe of measured position, veocity and acceeration aong z. Acceeration curves obtained from eq.1 are potted during free-faing. The steps in the curves correspond to the change of direction of the friction force. The sight changes at the ends of the curves are due to the variation of ψ. From the six cases studied, the effect of the baancing ratio (γ) on the acceeration is very cear. Aso, the effect of the inertia ratio (α) is significant, especiay on the unoaded device. Therefore, it is of interest to minimize the inertia in order to obtain better performances during the upwards acceeration (user pushing up). The effect of friction, which pushes against the direction of motion, is sma but visibe. It is noted that, with the setup used, the free fa of the person occurs around ψ = 0 o and the upwards acceeration of the chair occurs around ψ = 25 o. A arge vaue of ψ produces an increase of the effective inertia (see eq. (22)). Therefore, it is advantageous to keep ψ sma to minimize this effect. For instance, in a future impementation, the baancing ever wi be designed such that it is horizonta when the person is standing. Simuations with and without taking into account B 3 were performed, and the resuts were aways very cose for the present setup. Therefore, the effect of the veocity which rarey exceeded 1m/s in the experiments on the effective inertia seems to be very sma. VI. OTHER APPLICATIONS The concept presented in this paper can aso be used in other appications such as the deveopment of training faciities for astronauts (see [15] for instance). By propery adjusting the baancing ratio, the device can emuate the gravity fied of other panets or the moon and the trainees can wak in a arge workspace without being hindered by foormounted mechanisms. Aso, the device does not require the use of springs, which raise safety and reiabiity issues. This 268 REFERENCES [1] M. Bouri, Y. Stauffer, C. Schmit, Y. Aemand, S. Gnemmi, R. Cave, P. Metraier and R. Brodard, The WakTrainer: A Robotic System for Waking Rehabiitation in IEEE Int. Conference on Robotics and Biomimetics, pp , [2] S. Jezernik, G.Coombo, T. Keer, H. Frueh and M. Morani, Robotic orthosis Lokomat: a rehabiitation and research too, Neuromoduation: Technoogy at the Neura Interface, Vo. 6, No. 2, pp , [3] J. Jang, S. Yu, J. Han and C. Han, Deveopment of a waking assistive service robot for rehabiitation of edery peope, in Service Robot Appications (Takahashi, ed.), In-Tech, August [4] N. Siddiqi, F. Gazzani, J. Des Jardins and E.Y.S. Chao, The use of a robotic device for gait training and rehabiitation Medicine Meets Virtua Reaity, K.S. Morgan et a. (Eds.), IOS Press, [5] C.Y. Lee and J.J. Lee, Waking-support robot system for waking rehabiitation: design and contro, Artificia Life and Robotics, Vo. 4, No. 4, pp , [6] J. Patton, D.A. Brown, M. Peshkin, J.J. Santos-Munne, A. Makhin, E. Lewis, E.J. Cogate and D. Schwandt, Design and Deveopment of a Robotic Overground Gait and Baance Therapy Device, Topics in Stroke Rehabiitation, Vo. 15, No. 2, pp , [7] H. Schmidt, D. Sorowka, S. Hesse and R. Bernhardt, Deveopment of a robotic waking simuator for gait rehabiitation, Biomed tech, Vo. 48, No. 10, pp , Berin, [8] S. Perreaut and C. Gossein, Cabe-driven parae mechanisms: appication to a ocomotion interface, ASME Journa of Mechanica Design, Vo. 130, Artice , [9] J. Yoon and J. Ryu, A nove ocomotion interface with two 6-DOF parae manipuators that aows human waking on various virtua terrains, Internationa Journa of Robotics Research, Vo. 25, No. 7, pp , [10] D. Surdiovic and R. Bernhardt, STRING-MAN-Wire robot for gait rehabiitation: further deveopment and testing, European Symposium on Technica Aids for Rehabiitation, January 25 26, Technica University of Berin, [11] (2011). [12] (2011). [13] (2011). [14] T. Laiberté, C. Gossein and D. Gao, Cosed-oop actuation routings for Cartesian SCARA-type manipuators, Proceedings of the ASME IDETC/CIE Conferences, August 15 18, Montrea, Canada, [15] O. Ma, Q. Lu, J. McAvoy and K. Rube, Concept study of a passive reduced-gravity simuator for training astronauts, Proceedings of the ASME IDETC/CIE Conferences, August 15 18, Montrea, Canada, 2010.

Spring Gravity Compensation Using the Noncircular Pulley and Cable For the Less-Spring Design

Spring Gravity Compensation Using the Noncircular Pulley and Cable For the Less-Spring Design The 14th IFToMM Word Congress, Taipei, Taiwan, October 25-30, 2015 DOI Number: 10.6567/IFToMM.14TH.WC.PS3.010 Spring Gravity Compensation Using the Noncircuar Puey and Cabe For the Less-Spring Design M.C.

More information

Measurement of acceleration due to gravity (g) by a compound pendulum

Measurement of acceleration due to gravity (g) by a compound pendulum Measurement of acceeration due to gravity (g) by a compound penduum Aim: (i) To determine the acceeration due to gravity (g) by means of a compound penduum. (ii) To determine radius of gyration about an

More information

12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes

12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes Maima and Minima 1. Introduction In this Section we anayse curves in the oca neighbourhood of a stationary point and, from this anaysis, deduce necessary conditions satisfied by oca maima and oca minima.

More information

Technical Data for Profiles. Groove position, external dimensions and modular dimensions

Technical Data for Profiles. Groove position, external dimensions and modular dimensions Technica Data for Profies Extruded Profie Symbo A Mg Si 0.5 F 25 Materia number.206.72 Status: artificiay aged Mechanica vaues (appy ony in pressing direction) Tensie strength Rm min. 245 N/mm 2 Yied point

More information

PHYSICS LOCUS / / d dt. ( vi) mass, m moment of inertia, I. ( ix) linear momentum, p Angular momentum, l p mv l I

PHYSICS LOCUS / / d dt. ( vi) mass, m moment of inertia, I. ( ix) linear momentum, p Angular momentum, l p mv l I 6 n terms of moment of inertia, equation (7.8) can be written as The vector form of the above equation is...(7.9 a)...(7.9 b) The anguar acceeration produced is aong the direction of appied externa torque.

More information

Traffic data collection

Traffic data collection Chapter 32 Traffic data coection 32.1 Overview Unike many other discipines of the engineering, the situations that are interesting to a traffic engineer cannot be reproduced in a aboratory. Even if road

More information

1) For a block of mass m to slide without friction up a rise of height h, the minimum initial speed of the block must be

1) For a block of mass m to slide without friction up a rise of height h, the minimum initial speed of the block must be v m 1) For a bock of mass m to side without friction up a rise of height h, the minimum initia speed of the bock must be a ) gh b ) gh d ) gh e ) gh c ) gh P h b 3 15 ft 3) A man pus a pound crate up a

More information

An Algorithm for Pruning Redundant Modules in Min-Max Modular Network

An Algorithm for Pruning Redundant Modules in Min-Max Modular Network An Agorithm for Pruning Redundant Modues in Min-Max Moduar Network Hui-Cheng Lian and Bao-Liang Lu Department of Computer Science and Engineering, Shanghai Jiao Tong University 1954 Hua Shan Rd., Shanghai

More information

Adaptive Fuzzy Sliding Control for a Three-Link Passive Robotic Manipulator

Adaptive Fuzzy Sliding Control for a Three-Link Passive Robotic Manipulator Adaptive Fuzzy Siding Contro for a hree-link Passive Robotic anipuator Abstract An adaptive fuzzy siding contro (AFSC scheme is proposed to contro a passive robotic manipuator. he motivation for the design

More information

SEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l

SEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l Probem 7. Simpe Penduum SEMINAR. PENDULUMS A simpe penduum means a mass m suspended by a string weightess rigid rod of ength so that it can swing in a pane. The y-axis is directed down, x-axis is directed

More information

Numerical simulation of javelin best throwing angle based on biomechanical model

Numerical simulation of javelin best throwing angle based on biomechanical model ISSN : 0974-7435 Voume 8 Issue 8 Numerica simuation of javein best throwing ange based on biomechanica mode Xia Zeng*, Xiongwei Zuo Department of Physica Education, Changsha Medica University, Changsha

More information

Lecture 6: Moderately Large Deflection Theory of Beams

Lecture 6: Moderately Large Deflection Theory of Beams Structura Mechanics 2.8 Lecture 6 Semester Yr Lecture 6: Moderatey Large Defection Theory of Beams 6.1 Genera Formuation Compare to the cassica theory of beams with infinitesima deformation, the moderatey

More information

arxiv: v1 [physics.flu-dyn] 2 Nov 2007

arxiv: v1 [physics.flu-dyn] 2 Nov 2007 A theoretica anaysis of the resoution due to diffusion and size-dispersion of partices in deterministic atera dispacement devices arxiv:7.347v [physics.fu-dyn] 2 Nov 27 Martin Heer and Henrik Bruus MIC

More information

PhysicsAndMathsTutor.com

PhysicsAndMathsTutor.com . Two points A and B ie on a smooth horizonta tabe with AB = a. One end of a ight eastic spring, of natura ength a and moduus of easticity mg, is attached to A. The other end of the spring is attached

More information

Tracking Control of Multiple Mobile Robots

Tracking Control of Multiple Mobile Robots Proceedings of the 2001 IEEE Internationa Conference on Robotics & Automation Seou, Korea May 21-26, 2001 Tracking Contro of Mutipe Mobie Robots A Case Study of Inter-Robot Coision-Free Probem Jurachart

More information

STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION

STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION Journa of Sound and Vibration (996) 98(5), 643 65 STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM G. ERDOS AND T. SINGH Department of Mechanica and Aerospace Engineering, SUNY at Buffao,

More information

A Statistical Framework for Real-time Event Detection in Power Systems

A Statistical Framework for Real-time Event Detection in Power Systems 1 A Statistica Framework for Rea-time Event Detection in Power Systems Noan Uhrich, Tim Christman, Phiip Swisher, and Xichen Jiang Abstract A quickest change detection (QCD) agorithm is appied to the probem

More information

Module 22: Simple Harmonic Oscillation and Torque

Module 22: Simple Harmonic Oscillation and Torque Modue : Simpe Harmonic Osciation and Torque.1 Introduction We have aready used Newton s Second Law or Conservation of Energy to anayze systems ike the boc-spring system that osciate. We sha now use torque

More information

c 2007 Society for Industrial and Applied Mathematics

c 2007 Society for Industrial and Applied Mathematics SIAM REVIEW Vo. 49,No. 1,pp. 111 1 c 7 Society for Industria and Appied Mathematics Domino Waves C. J. Efthimiou M. D. Johnson Abstract. Motivated by a proposa of Daykin [Probem 71-19*, SIAM Rev., 13 (1971),

More information

PHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased

PHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased PHYS 110B - HW #1 Fa 2005, Soutions by David Pace Equations referenced as Eq. # are from Griffiths Probem statements are paraphrased [1.] Probem 6.8 from Griffiths A ong cyinder has radius R and a magnetization

More information

LECTURE 10. The world of pendula

LECTURE 10. The world of pendula LECTURE 0 The word of pendua For the next few ectures we are going to ook at the word of the pane penduum (Figure 0.). In a previous probem set we showed that we coud use the Euer- Lagrange method to derive

More information

THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS

THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS ECCM6-6 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Sevie, Spain, -6 June 04 THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS M. Wysocki a,b*, M. Szpieg a, P. Heström a and F. Ohsson c a Swerea SICOMP

More information

Solution Set Seven. 1 Goldstein Components of Torque Along Principal Axes Components of Torque Along Cartesian Axes...

Solution Set Seven. 1 Goldstein Components of Torque Along Principal Axes Components of Torque Along Cartesian Axes... : Soution Set Seven Northwestern University, Cassica Mechanics Cassica Mechanics, Third Ed.- Godstein November 8, 25 Contents Godstein 5.8. 2. Components of Torque Aong Principa Axes.......................

More information

Forces of Friction. through a viscous medium, there will be a resistance to the motion. and its environment

Forces of Friction. through a viscous medium, there will be a resistance to the motion. and its environment Forces of Friction When an object is in motion on a surface or through a viscous medium, there wi be a resistance to the motion This is due to the interactions between the object and its environment This

More information

APPENDIX C FLEXING OF LENGTH BARS

APPENDIX C FLEXING OF LENGTH BARS Fexing of ength bars 83 APPENDIX C FLEXING OF LENGTH BARS C.1 FLEXING OF A LENGTH BAR DUE TO ITS OWN WEIGHT Any object ying in a horizonta pane wi sag under its own weight uness it is infinitey stiff or

More information

Candidate Number. General Certificate of Education Advanced Level Examination June 2010

Candidate Number. General Certificate of Education Advanced Level Examination June 2010 Centre Number Surname Candidate Number For Examiner s Use Other Names Candidate Signature Examiner s Initias Genera Certificate of Education Advanced Leve Examination June 2010 Question 1 2 Mark Physics

More information

Separation of Variables and a Spherical Shell with Surface Charge

Separation of Variables and a Spherical Shell with Surface Charge Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation

More information

Turbo Codes. Coding and Communication Laboratory. Dept. of Electrical Engineering, National Chung Hsing University

Turbo Codes. Coding and Communication Laboratory. Dept. of Electrical Engineering, National Chung Hsing University Turbo Codes Coding and Communication Laboratory Dept. of Eectrica Engineering, Nationa Chung Hsing University Turbo codes 1 Chapter 12: Turbo Codes 1. Introduction 2. Turbo code encoder 3. Design of intereaver

More information

Problem Set 6: Solutions

Problem Set 6: Solutions University of Aabama Department of Physics and Astronomy PH 102 / LeCair Summer II 2010 Probem Set 6: Soutions 1. A conducting rectanguar oop of mass M, resistance R, and dimensions w by fas from rest

More information

Lecture Notes for Math 251: ODE and PDE. Lecture 34: 10.7 Wave Equation and Vibrations of an Elastic String

Lecture Notes for Math 251: ODE and PDE. Lecture 34: 10.7 Wave Equation and Vibrations of an Elastic String ecture Notes for Math 251: ODE and PDE. ecture 3: 1.7 Wave Equation and Vibrations of an Eastic String Shawn D. Ryan Spring 212 ast Time: We studied other Heat Equation probems with various other boundary

More information

Non-linear robust control for inverted-pendulum 2D walking

Non-linear robust control for inverted-pendulum 2D walking Non-inear robust contro for inverted-penduum 2D waking Matthew Key and Andy Ruina 2 Abstract We present an approach to high-eve contro for bipeda waking exempified with a 2D point-mass inextensibeegs inverted-penduum

More information

Term Test AER301F. Dynamics. 5 November The value of each question is indicated in the table opposite.

Term Test AER301F. Dynamics. 5 November The value of each question is indicated in the table opposite. U N I V E R S I T Y O F T O R O N T O Facuty of Appied Science and Engineering Term Test AER31F Dynamics 5 November 212 Student Name: Last Name First Names Student Number: Instructions: 1. Attempt a questions.

More information

XI PHYSICS. M. Affan Khan LECTURER PHYSICS, AKHSS, K. https://promotephysics.wordpress.com

XI PHYSICS. M. Affan Khan LECTURER PHYSICS, AKHSS, K. https://promotephysics.wordpress.com XI PHYSICS M. Affan Khan LECTURER PHYSICS, AKHSS, K affan_414@ive.com https://promotephysics.wordpress.com [TORQUE, ANGULAR MOMENTUM & EQUILIBRIUM] CHAPTER NO. 5 Okay here we are going to discuss Rotationa

More information

PREPARED BY: ER. VINEET LOOMBA (B.TECH. IIT ROORKEE)

PREPARED BY: ER. VINEET LOOMBA (B.TECH. IIT ROORKEE) Cass XI TARGET : JEE Main/Adv PREPARED BY: ER. VINEET LOOMBA (B.TECH. IIT ROORKEE) ALP ADVANCED LEVEL LPROBLEMS ROTATION- Topics Covered: Rigid body, moment of inertia, parae and perpendicuar axes theorems,

More information

In-plane shear stiffness of bare steel deck through shell finite element models. G. Bian, B.W. Schafer. June 2017

In-plane shear stiffness of bare steel deck through shell finite element models. G. Bian, B.W. Schafer. June 2017 In-pane shear stiffness of bare stee deck through she finite eement modes G. Bian, B.W. Schafer June 7 COLD-FORMED STEEL RESEARCH CONSORTIUM REPORT SERIES CFSRC R-7- SDII Stee Diaphragm Innovation Initiative

More information

UI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE

UI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE UI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE Juan Huang, Ronghui Wang and Tao Tang Coege of Traffic and Communications, South China University of Technoogy, Guangzhou, Guangdong 51641,

More information

Acoustic Nondestructive Testing and Measurement of Tension for Steel Reinforcing Members

Acoustic Nondestructive Testing and Measurement of Tension for Steel Reinforcing Members Acoustic Nondestructive Testing and Measurement of Tension for Stee Reinforcing Members Part 1 Theory by Michae K. McInerney PURPOSE: This Coasta and Hydrauics Engineering Technica Note (CHETN) describes

More information

Combining reaction kinetics to the multi-phase Gibbs energy calculation

Combining reaction kinetics to the multi-phase Gibbs energy calculation 7 th European Symposium on Computer Aided Process Engineering ESCAPE7 V. Pesu and P.S. Agachi (Editors) 2007 Esevier B.V. A rights reserved. Combining reaction inetics to the muti-phase Gibbs energy cacuation

More information

THE THREE POINT STEINER PROBLEM ON THE FLAT TORUS: THE MINIMAL LUNE CASE

THE THREE POINT STEINER PROBLEM ON THE FLAT TORUS: THE MINIMAL LUNE CASE THE THREE POINT STEINER PROBLEM ON THE FLAT TORUS: THE MINIMAL LUNE CASE KATIE L. MAY AND MELISSA A. MITCHELL Abstract. We show how to identify the minima path network connecting three fixed points on

More information

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is

More information

The influence of temperature of photovoltaic modules on performance of solar power plant

The influence of temperature of photovoltaic modules on performance of solar power plant IOSR Journa of Engineering (IOSRJEN) ISSN (e): 2250-3021, ISSN (p): 2278-8719 Vo. 05, Issue 04 (Apri. 2015), V1 PP 09-15 www.iosrjen.org The infuence of temperature of photovotaic modues on performance

More information

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische

More information

$, (2.1) n="# #. (2.2)

$, (2.1) n=# #. (2.2) Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier

More information

arxiv: v1 [math.ca] 6 Mar 2017

arxiv: v1 [math.ca] 6 Mar 2017 Indefinite Integras of Spherica Besse Functions MIT-CTP/487 arxiv:703.0648v [math.ca] 6 Mar 07 Joyon K. Boomfied,, Stephen H. P. Face,, and Zander Moss, Center for Theoretica Physics, Laboratory for Nucear

More information

1. Measurements and error calculus

1. Measurements and error calculus EV 1 Measurements and error cacuus 11 Introduction The goa of this aboratory course is to introduce the notions of carrying out an experiment, acquiring and writing up the data, and finay anayzing the

More information

A Brief Introduction to Markov Chains and Hidden Markov Models

A Brief Introduction to Markov Chains and Hidden Markov Models A Brief Introduction to Markov Chains and Hidden Markov Modes Aen B MacKenzie Notes for December 1, 3, &8, 2015 Discrete-Time Markov Chains You may reca that when we first introduced random processes,

More information

Research on liquid sloshing performance in vane type tank under microgravity

Research on liquid sloshing performance in vane type tank under microgravity IOP Conference Series: Materias Science and Engineering PAPER OPEN ACCESS Research on iquid soshing performance in vane type tan under microgravity Reated content - Numerica simuation of fuid fow in the

More information

8 Digifl'.11 Cth:uits and devices

8 Digifl'.11 Cth:uits and devices 8 Digif'. Cth:uits and devices 8. Introduction In anaog eectronics, votage is a continuous variabe. This is usefu because most physica quantities we encounter are continuous: sound eves, ight intensity,

More information

Statistical Learning Theory: A Primer

Statistical Learning Theory: A Primer Internationa Journa of Computer Vision 38(), 9 3, 2000 c 2000 uwer Academic Pubishers. Manufactured in The Netherands. Statistica Learning Theory: A Primer THEODOROS EVGENIOU, MASSIMILIANO PONTIL AND TOMASO

More information

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s

More information

Some Measures for Asymmetry of Distributions

Some Measures for Asymmetry of Distributions Some Measures for Asymmetry of Distributions Georgi N. Boshnakov First version: 31 January 2006 Research Report No. 5, 2006, Probabiity and Statistics Group Schoo of Mathematics, The University of Manchester

More information

Lobontiu: System Dynamics for Engineering Students Website Chapter 3 1. z b z

Lobontiu: System Dynamics for Engineering Students Website Chapter 3 1. z b z Chapter W3 Mechanica Systems II Introduction This companion website chapter anayzes the foowing topics in connection to the printed-book Chapter 3: Lumped-parameter inertia fractions of basic compiant

More information

International Journal of Mass Spectrometry

International Journal of Mass Spectrometry Internationa Journa of Mass Spectrometry 280 (2009) 179 183 Contents ists avaiabe at ScienceDirect Internationa Journa of Mass Spectrometry journa homepage: www.esevier.com/ocate/ijms Stark mixing by ion-rydberg

More information

T.C. Banwell, S. Galli. {bct, Telcordia Technologies, Inc., 445 South Street, Morristown, NJ 07960, USA

T.C. Banwell, S. Galli. {bct, Telcordia Technologies, Inc., 445 South Street, Morristown, NJ 07960, USA ON THE SYMMETRY OF THE POWER INE CHANNE T.C. Banwe, S. Gai {bct, sgai}@research.tecordia.com Tecordia Technoogies, Inc., 445 South Street, Morristown, NJ 07960, USA Abstract The indoor power ine network

More information

CONCHOID OF NICOMEDES AND LIMACON OF PASCAL AS ELECTRODE OF STATIC FIELD AND AS WAVEGUIDE OF HIGH FREQUENCY WAVE

CONCHOID OF NICOMEDES AND LIMACON OF PASCAL AS ELECTRODE OF STATIC FIELD AND AS WAVEGUIDE OF HIGH FREQUENCY WAVE Progress In Eectromagnetics Research, PIER 30, 73 84, 001 CONCHOID OF NICOMEDES AND LIMACON OF PASCAL AS ELECTRODE OF STATIC FIELD AND AS WAVEGUIDE OF HIGH FREQUENCY WAVE W. Lin and Z. Yu University of

More information

Candidate Number. General Certificate of Education Advanced Level Examination January 2012

Candidate Number. General Certificate of Education Advanced Level Examination January 2012 entre Number andidate Number Surname Other Names andidate Signature Genera ertificate of Education dvanced Leve Examination January 212 Physics PHY4/1 Unit 4 Fieds and Further Mechanics Section Tuesday

More information

Parallel-Axis Theorem

Parallel-Axis Theorem Parae-Axis Theorem In the previous exampes, the axis of rotation coincided with the axis of symmetry of the object For an arbitrary axis, the paraeaxis theorem often simpifies cacuations The theorem states

More information

Interconnect effects on performance of Field Programmable Analog Array

Interconnect effects on performance of Field Programmable Analog Array nterconnect effects on performance of Fied Programmabe Anaog Array D. Anderson,. Bir, O. A. Pausinsi 3, M. Spitz, K. Reiss Motoroa, SPS, Phoenix, Arizona, USA, University of Karsruhe, Karsruhe, Germany,

More information

Bayesian Learning. You hear a which which could equally be Thanks or Tanks, which would you go with?

Bayesian Learning. You hear a which which could equally be Thanks or Tanks, which would you go with? Bayesian Learning A powerfu and growing approach in machine earning We use it in our own decision making a the time You hear a which which coud equay be Thanks or Tanks, which woud you go with? Combine

More information

Lecture 6 Povh Krane Enge Williams Properties of 2-nucleon potential

Lecture 6 Povh Krane Enge Williams Properties of 2-nucleon potential Lecture 6 Povh Krane Enge Wiiams Properties of -nuceon potentia 16.1 4.4 3.6 9.9 Meson Theory of Nucear potentia 4.5 3.11 9.10 I recommend Eisberg and Resnik notes as distributed Probems, Lecture 6 1 Consider

More information

More Scattering: the Partial Wave Expansion

More Scattering: the Partial Wave Expansion More Scattering: the Partia Wave Expansion Michae Fower /7/8 Pane Waves and Partia Waves We are considering the soution to Schrödinger s equation for scattering of an incoming pane wave in the z-direction

More information

Approximated MLC shape matrix decomposition with interleaf collision constraint

Approximated MLC shape matrix decomposition with interleaf collision constraint Approximated MLC shape matrix decomposition with intereaf coision constraint Thomas Kainowski Antje Kiese Abstract Shape matrix decomposition is a subprobem in radiation therapy panning. A given fuence

More information

Physics 127c: Statistical Mechanics. Fermi Liquid Theory: Collective Modes. Boltzmann Equation. The quasiparticle energy including interactions

Physics 127c: Statistical Mechanics. Fermi Liquid Theory: Collective Modes. Boltzmann Equation. The quasiparticle energy including interactions Physics 27c: Statistica Mechanics Fermi Liquid Theory: Coective Modes Botzmann Equation The quasipartice energy incuding interactions ε p,σ = ε p + f(p, p ; σ, σ )δn p,σ, () p,σ with ε p ε F + v F (p p

More information

Quantum Electrodynamical Basis for Wave. Propagation through Photonic Crystal

Quantum Electrodynamical Basis for Wave. Propagation through Photonic Crystal Adv. Studies Theor. Phys., Vo. 6, 01, no. 3, 19-133 Quantum Eectrodynamica Basis for Wave Propagation through Photonic Crysta 1 N. Chandrasekar and Har Narayan Upadhyay Schoo of Eectrica and Eectronics

More information

Nonlinear Analysis of Spatial Trusses

Nonlinear Analysis of Spatial Trusses Noninear Anaysis of Spatia Trusses João Barrigó October 14 Abstract The present work addresses the noninear behavior of space trusses A formuation for geometrica noninear anaysis is presented, which incudes

More information

Approximation and Fast Calculation of Non-local Boundary Conditions for the Time-dependent Schrödinger Equation

Approximation and Fast Calculation of Non-local Boundary Conditions for the Time-dependent Schrödinger Equation Approximation and Fast Cacuation of Non-oca Boundary Conditions for the Time-dependent Schrödinger Equation Anton Arnod, Matthias Ehrhardt 2, and Ivan Sofronov 3 Universität Münster, Institut für Numerische

More information

Position Control of Rolling Skateboard

Position Control of Rolling Skateboard osition Contro of Roing kateboard Baazs Varszegi enes Takacs Gabor tepan epartment of Appied Mechanics, Budapest University of Technoogy and Economics, Budapest, Hungary (e-mai: varszegi@mm.bme.hu) MTA-BME

More information

Related Topics Maxwell s equations, electrical eddy field, magnetic field of coils, coil, magnetic flux, induced voltage

Related Topics Maxwell s equations, electrical eddy field, magnetic field of coils, coil, magnetic flux, induced voltage Magnetic induction TEP Reated Topics Maxwe s equations, eectrica eddy fied, magnetic fied of cois, coi, magnetic fux, induced votage Principe A magnetic fied of variabe frequency and varying strength is

More information

Approximated MLC shape matrix decomposition with interleaf collision constraint

Approximated MLC shape matrix decomposition with interleaf collision constraint Agorithmic Operations Research Vo.4 (29) 49 57 Approximated MLC shape matrix decomposition with intereaf coision constraint Antje Kiese and Thomas Kainowski Institut für Mathematik, Universität Rostock,

More information

Paper presented at the Workshop on Space Charge Physics in High Intensity Hadron Rings, sponsored by Brookhaven National Laboratory, May 4-7,1998

Paper presented at the Workshop on Space Charge Physics in High Intensity Hadron Rings, sponsored by Brookhaven National Laboratory, May 4-7,1998 Paper presented at the Workshop on Space Charge Physics in High ntensity Hadron Rings, sponsored by Brookhaven Nationa Laboratory, May 4-7,998 Noninear Sef Consistent High Resoution Beam Hao Agorithm in

More information

A simple reliability block diagram method for safety integrity verification

A simple reliability block diagram method for safety integrity verification Reiabiity Engineering and System Safety 92 (2007) 1267 1273 www.esevier.com/ocate/ress A simpe reiabiity bock diagram method for safety integrity verification Haitao Guo, Xianhui Yang epartment of Automation,

More information

(This is a sample cover image for this issue. The actual cover is not yet available at this time.)

(This is a sample cover image for this issue. The actual cover is not yet available at this time.) (This is a sampe cover image for this issue The actua cover is not yet avaiabe at this time) This artice appeared in a journa pubished by Esevier The attached copy is furnished to the author for interna

More information

Simplified analysis of EXAFS data and determination of bond lengths

Simplified analysis of EXAFS data and determination of bond lengths Indian Journa of Pure & Appied Physics Vo. 49, January 0, pp. 5-9 Simpified anaysis of EXAFS data and determination of bond engths A Mishra, N Parsai & B D Shrivastava * Schoo of Physics, Devi Ahiya University,

More information

17 Lecture 17: Recombination and Dark Matter Production

17 Lecture 17: Recombination and Dark Matter Production PYS 652: Astrophysics 88 17 Lecture 17: Recombination and Dark Matter Production New ideas pass through three periods: It can t be done. It probaby can be done, but it s not worth doing. I knew it was

More information

Steepest Descent Adaptation of Min-Max Fuzzy If-Then Rules 1

Steepest Descent Adaptation of Min-Max Fuzzy If-Then Rules 1 Steepest Descent Adaptation of Min-Max Fuzzy If-Then Rues 1 R.J. Marks II, S. Oh, P. Arabshahi Λ, T.P. Caude, J.J. Choi, B.G. Song Λ Λ Dept. of Eectrica Engineering Boeing Computer Services University

More information

High Efficiency Development of a Reciprocating Compressor by Clarification of Loss Generation in Bearings

High Efficiency Development of a Reciprocating Compressor by Clarification of Loss Generation in Bearings Purdue University Purdue e-pubs Internationa Compressor Engineering Conference Schoo of Mechanica Engineering 2010 High Efficiency Deveopment of a Reciprocating Compressor by Carification of Loss Generation

More information

Lecture Note 3: Stationary Iterative Methods

Lecture Note 3: Stationary Iterative Methods MATH 5330: Computationa Methods of Linear Agebra Lecture Note 3: Stationary Iterative Methods Xianyi Zeng Department of Mathematica Sciences, UTEP Stationary Iterative Methods The Gaussian eimination (or

More information

ASummaryofGaussianProcesses Coryn A.L. Bailer-Jones

ASummaryofGaussianProcesses Coryn A.L. Bailer-Jones ASummaryofGaussianProcesses Coryn A.L. Baier-Jones Cavendish Laboratory University of Cambridge caj@mrao.cam.ac.uk Introduction A genera prediction probem can be posed as foows. We consider that the variabe

More information

FOURIER SERIES ON ANY INTERVAL

FOURIER SERIES ON ANY INTERVAL FOURIER SERIES ON ANY INTERVAL Overview We have spent considerabe time earning how to compute Fourier series for functions that have a period of 2p on the interva (-p,p). We have aso seen how Fourier series

More information

Induction and Inductance

Induction and Inductance Induction and Inductance How we generate E by B, and the passive component inductor in a circuit. 1. A review of emf and the magnetic fux. 2. Faraday s Law of Induction 3. Lentz Law 4. Inductance and inductor

More information

Unit 48: Structural Behaviour and Detailing for Construction. Deflection of Beams

Unit 48: Structural Behaviour and Detailing for Construction. Deflection of Beams Unit 48: Structura Behaviour and Detaiing for Construction 4.1 Introduction Defection of Beams This topic investigates the deformation of beams as the direct effect of that bending tendency, which affects

More information

A Simple and Efficient Algorithm of 3-D Single-Source Localization with Uniform Cross Array Bing Xue 1 2 a) * Guangyou Fang 1 2 b and Yicai Ji 1 2 c)

A Simple and Efficient Algorithm of 3-D Single-Source Localization with Uniform Cross Array Bing Xue 1 2 a) * Guangyou Fang 1 2 b and Yicai Ji 1 2 c) A Simpe Efficient Agorithm of 3-D Singe-Source Locaization with Uniform Cross Array Bing Xue a * Guangyou Fang b Yicai Ji c Key Laboratory of Eectromagnetic Radiation Sensing Technoogy, Institute of Eectronics,

More information

Copyright information to be inserted by the Publishers. Unsplitting BGK-type Schemes for the Shallow. Water Equations KUN XU

Copyright information to be inserted by the Publishers. Unsplitting BGK-type Schemes for the Shallow. Water Equations KUN XU Copyright information to be inserted by the Pubishers Unspitting BGK-type Schemes for the Shaow Water Equations KUN XU Mathematics Department, Hong Kong University of Science and Technoogy, Cear Water

More information

Efficient Visual-Inertial Navigation using a Rolling-Shutter Camera with Inaccurate Timestamps

Efficient Visual-Inertial Navigation using a Rolling-Shutter Camera with Inaccurate Timestamps Efficient Visua-Inertia Navigation using a Roing-Shutter Camera with Inaccurate Timestamps Chao X. Guo, Dimitrios G. Kottas, Ryan C. DuToit Ahmed Ahmed, Ruipeng Li and Stergios I. Roumeiotis Mutipe Autonomous

More information

Session : Electrodynamic Tethers

Session : Electrodynamic Tethers Session : Eectrodynaic Tethers Eectrodynaic tethers are ong, thin conductive wires depoyed in space that can be used to generate power by reoving kinetic energy fro their orbita otion, or to produce thrust

More information

Section 6: Magnetostatics

Section 6: Magnetostatics agnetic fieds in matter Section 6: agnetostatics In the previous sections we assumed that the current density J is a known function of coordinates. In the presence of matter this is not aways true. The

More information

Compression Ratio Expansion for Complementary Code Set Compressing a Signal to a Width of Several Sub-pulses

Compression Ratio Expansion for Complementary Code Set Compressing a Signal to a Width of Several Sub-pulses Compression Ratio Expansion for Compementary Code Set Compressing a Signa to a Width of Severa Sub-puses Reiji Sato 2nd Research Center, Technica R and D Institute Japan Defense Agency 1-2-24 Ikejiri,

More information

SE-514 (OPTIMAL CONTROL) OPTIMAL CONTROL FOR SINGLE AND DOUBLE INVERTED PENDULUM. DONE BY: Fatai Olalekan ( Ayman Abdallah (973610)

SE-514 (OPTIMAL CONTROL) OPTIMAL CONTROL FOR SINGLE AND DOUBLE INVERTED PENDULUM. DONE BY: Fatai Olalekan ( Ayman Abdallah (973610) SE-54 (OPTIAL CONTROL OPTIAL CONTROL FOR SINGLE AND DOUBLE INVERTED PENDULU DONE BY: Fatai Oaekan (363 Ayman Abdaah (9736 PREPARED FOR: Dr. Sami E-Ferik Tabe of contents Abstract... 3 Introduction... 3

More information

Extended Simulations for the Observer-based Control of Position and Tension for an Aerial Robot Tethered to a Moving Platform

Extended Simulations for the Observer-based Control of Position and Tension for an Aerial Robot Tethered to a Moving Platform Extended Simuations for the Observer-based Contro of Position and Tension for an eria obot Tethered to a Moving Patform Marco Tognon, ntonio Franchi To cite this version: Marco Tognon, ntonio Franchi.

More information

Introduction. Figure 1 W8LC Line Array, box and horn element. Highlighted section modelled.

Introduction. Figure 1 W8LC Line Array, box and horn element. Highlighted section modelled. imuation of the acoustic fied produced by cavities using the Boundary Eement Rayeigh Integra Method () and its appication to a horn oudspeaer. tephen Kirup East Lancashire Institute, Due treet, Bacburn,

More information

Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Delft University of Technology. Marijn Drillenburg

Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Delft University of Technology. Marijn Drillenburg Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Deft University of Technoogy Marijn Drienburg October 2017 Contents 1 Introduction 2 1.1 Hand Cacuation....................................

More information

1 Equations of Motion 3: Equivalent System Method

1 Equations of Motion 3: Equivalent System Method 8 Mechanica Vibrations Equations of Motion : Equivaent System Method In systems in which masses are joined by rigid ins, evers, or gears and in some distributed systems, various springs, dampers, and masses

More information

IMPLEMENT OF VARIOUS MODELS The coordinate system of the wind-bridge interaction system for anaysis in a these modes is shown in Fig.. A unit-ength se

IMPLEMENT OF VARIOUS MODELS The coordinate system of the wind-bridge interaction system for anaysis in a these modes is shown in Fig.. A unit-ength se The Seventh Internationa Cooquium on Buff Body Aerodynamics and Appications (BBAA7) Shanghai, China; September -6, Comparison of various modeing schemes for bridge aerodynamics and aeroeasticity Teng Wu

More information

Physics 235 Chapter 8. Chapter 8 Central-Force Motion

Physics 235 Chapter 8. Chapter 8 Central-Force Motion Physics 35 Chapter 8 Chapter 8 Centra-Force Motion In this Chapter we wi use the theory we have discussed in Chapter 6 and 7 and appy it to very important probems in physics, in which we study the motion

More information

Transverse Anisotropy in Softwoods

Transverse Anisotropy in Softwoods Transverse Anisotropy in Softwoods Modeing and Experiments CARL MODÉN Licenciate Thesis Stockhom, Sweden 2006 TRITA-AVE 2006:30 ISSN 1651-7660 ISBN 91-7178-385-7 KTH Engineering Sciences Department of

More information

On a geometrical approach in contact mechanics

On a geometrical approach in contact mechanics Institut für Mechanik On a geometrica approach in contact mechanics Aexander Konyukhov, Kar Schweizerhof Universität Karsruhe, Institut für Mechanik Institut für Mechanik Kaiserstr. 12, Geb. 20.30 76128

More information

The basic equation for the production of turbulent kinetic energy in clouds is. dz + g w

The basic equation for the production of turbulent kinetic energy in clouds is. dz + g w Turbuence in couds The basic equation for the production of turbuent kinetic energy in couds is de TKE dt = u 0 w 0 du v 0 w 0 dv + g w q 0 q 0 e The first two terms on the RHS are associated with shear

More information

Experimental Investigation and Numerical Analysis of New Multi-Ribbed Slab Structure

Experimental Investigation and Numerical Analysis of New Multi-Ribbed Slab Structure Experimenta Investigation and Numerica Anaysis of New Muti-Ribbed Sab Structure Jie TIAN Xi an University of Technoogy, China Wei HUANG Xi an University of Architecture & Technoogy, China Junong LU Xi

More information

Numerical Simulation for Optimizing Temperature Gradients during Single Crystal Casting Process

Numerical Simulation for Optimizing Temperature Gradients during Single Crystal Casting Process ISIJ Internationa Vo 54 (2014) No 2 pp 254 258 Numerica Simuation for Optimizing Temperature Gradients during Singe Crysta Casting Process Aeksandr Aeksandrovich INOZEMTSEV 1) Aeksandra Sergeevna DUBROVSKAYA

More information

Chapter 4. Moving Observer Method. 4.1 Overview. 4.2 Theory

Chapter 4. Moving Observer Method. 4.1 Overview. 4.2 Theory Chapter 4 Moving Observer Method 4.1 Overview For a compete description of traffic stream modeing, one woud reuire fow, speed, and density. Obtaining these parameters simutaneousy is a difficut task if

More information