- Critical parts are tested (PS) - The remaining parts are simulated (NS).

Size: px
Start display at page:

Download "- Critical parts are tested (PS) - The remaining parts are simulated (NS)."

Transcription

1

2 Hybrid simulatio of a structural system Loads (Disp.) simulatig deck motio due to earthquake K V T B CBdV Deck modeled aalytically by F.E.M. Piers tested physically - Critical parts are tested (PS) - The remaiig parts are simulated (NS). 1) Physical Substrucutre should move cotiuously, with a costat timescale betwee the seismic evet ad the testig duratio; 2) I the curret applicatio, a refied FE model is itroduced for the Numerical Substructure.

3 Servo-Hydraulic Actuators The implemetatio scheme developed at JRC, Ispra Force trasducers Displacemet Trasducers Referece frame A digital implemetatio of PsD cotrol scheme: Dt s = 1 or 2 msec Substructurig capacity Measuremet of reactio force R... dt Imposed displacemet d.. M d Cd. R d MIa g t Accelerogram x g Numerical itegratio Numerical model

4 The classical moolithic time itegratio:.. M d Cd. Whe the NS is refied, a moolithic itegratio step caot be complete withi the cotrol samplig time. R d MIa g è 1) Use of Partitioed itegratio ad differet time steps betwee the Ns ad the PS. è 2) Explicit scheme is preferable for PS (without iteratio), whereas implicit scheme is suitable to the NS (stability cosideratio). t K. Park, C. Felippa, 1983; T. Hughes, W. Liu, 1978; W. Liu, T. Belytschko, 1982

5 A Multi-time-step Explicit-implicit Method (Gravouil ad Combescure 2001) Steps are sequetial. Steps A 1 ad B 1 are cocurret.

6 Time steppig schemes applied to liear problems ca be recast i a recursive form as where X is the state vector, A is the amplificatio matrix ad L is the load vector. For PM method, X reads with X AX L 1 X ( X, X, X ) A A B T 1 X ( u, u, u, u ), X ( u, u, u ) A A A A A T B B B B T, free X has 11 variables for a partitioed SDoF system. Boelli et.al, 2011

7 Photograph courtesy of JRC Isolator: the physical substructure. The emulated structure. Bridge deck: the umerical substructure.

8 Problems related to umerical dissipatio I FE aalysis, the higher modes are ofte spurious due to the artifacts of the umerical discretizatio process. These spurious modes could be excited due to experimetal errors durig hybrid simulatio. T. Hughes, 1987; Bursi, Shig, 1997; Hulbert, Chag, 1993/1996; C. Farhat, F.-X.Roux,

9 A ovel oe-step four-stage form of the G-alpha method Predictor: Equilibrium equatio (Im.): Corrector: u u tu t u t a 2 f 2 1 D D D ( ) 1 m 2 1 m f u u Dt u Dt(1 ) a 1 1 m m (Ex.): Mu Cu Ku F ( t ) ext Mu Cu Ku F ( t ) ext u u t u u u t u 2 f f 1 D 1; 1 D 1 1 m 1 m The equilibrium equatio is idetical to those of Newmark schemes. The itegratio parameters are chose as those of the origial implicit/explicit G-alpha method, that is, r as the spectral radius at the high frequecy limit for the implicit method ad r b as the spectral radius at the bifurcatio limit for the explicit method. He, 2008; Bursi, He 2012

10 Represetative Numerical Simulatio (I) 0.04 PM m2 u2 (rad) 0.00 m1 la ma u1 Aalytical Part mb time (s) lb Experimetal Part 0.04 PM-alpha (rad) 0.00 Partitioed three-dof system f 1 =0.24 Hz, f 2 =1.6 Hz, f 3 =35 Hz time (s) Rotatioal respose i free vibratio for the Three-DoF system (Dt A =5ms, ss = 100): (a) the PM method; (b) the PM-alpha method, r = r b = 0.5

11 A three-actuator testig setup of the previous 3DOF model Simulatio cosiderig the cotrol errors of the actuators i udershootig type (Mosqueda et.al 2007) ad measured errors of load cells as white oise (0.06%) Respose of partitioed three-dof system (Texp./Ts.e.=50; Dt A =2ms ; Dt B =0.02ms ;ss=100; )

12 This work presets partitioed time itegratio methods for solvig dyamic problems of hybrid structural systems, typical of pseudo-dyamic testig with substructurig or hybrid simulatio, as decomposed by computatioal/physical cosideratios. The methods implemet the iterfield parallel aalysis procedure suitable to maitai a costat time scale betwee the umerical ad physical substructures, by usig the implicit ad explicit schemes of differet time steps. The characteristics of the proposed time schemes are aalyzed with respect to their covergece (stability ad accuracy). High-frequecy umerical dissipatio based o a ew form of the G-alpha methods is itroduced, which is show to be helpful to hybrid simulatio with the presece of both the cotrol errors ad the measuremet errors.

13 Thak you for your attetio! Ackowledgemets: Prof. O.S. UNITN Dr. P. JRC, Ispra

2C09 Design for seismic and climate changes

2C09 Design for seismic and climate changes 2C09 Desig for seismic ad climate chages Lecture 02: Dyamic respose of sigle-degree-of-freedom systems I Daiel Grecea, Politehica Uiversity of Timisoara 10/03/2014 Europea Erasmus Mudus Master Course Sustaiable

More information

Mechatronics. Time Response & Frequency Response 2 nd -Order Dynamic System 2-Pole, Low-Pass, Active Filter

Mechatronics. Time Response & Frequency Response 2 nd -Order Dynamic System 2-Pole, Low-Pass, Active Filter Time Respose & Frequecy Respose d -Order Dyamic System -Pole, Low-Pass, Active Filter R 4 R 7 C 5 e i R 1 C R 3 - + R 6 - + e out Assigmet: Perform a Complete Dyamic System Ivestigatio of the Two-Pole,

More information

Spectral Analysis. This week in lab. Next classes: 3/26 and 3/28. Your next experiment Homework is to prepare

Spectral Analysis. This week in lab. Next classes: 3/26 and 3/28. Your next experiment Homework is to prepare Spectral Aalysis This week i lab Your ext experimet Homework is to prepare Next classes: 3/26 ad 3/28 Aero Testig, Fracture Toughess Testig Read the Experimets 5 ad 7 sectios of the course maual Spectral

More information

MONOLITHIC AND PARTITIONED L-STABLE ROSENBROCK METHODS FOR DYNAMIC SUBSTRUCTURE TESTS. Oreste S. Bursi, Chuanguo Jia and Zhen Wang

MONOLITHIC AND PARTITIONED L-STABLE ROSENBROCK METHODS FOR DYNAMIC SUBSTRUCTURE TESTS. Oreste S. Bursi, Chuanguo Jia and Zhen Wang MONOLITHIC AND PARTITIONED L-STABLE ROSENBROCK METHODS FOR DYNAMIC SUBSTRUCTURE TESTS Oreste S. Bursi, Chuanguo Jia and Zhen Wang University of Trento/ Department of Mechanical and Structural Engineering,

More information

Basics of Dynamics. Amit Prashant. Indian Institute of Technology Gandhinagar. Short Course on. Geotechnical Aspects of Earthquake Engineering

Basics of Dynamics. Amit Prashant. Indian Institute of Technology Gandhinagar. Short Course on. Geotechnical Aspects of Earthquake Engineering Basics of yamics Amit Prashat Idia Istitute of Techology Gadhiagar Short Course o Geotechical Aspects of Earthquake Egieerig 4 8 March, 213 Our ear Pedulum Revisited g.si g l s Force Equilibrium: Cord

More information

FIXED-FREE AGAINST FREE-FREE BEAMS FOR DYNAMIC YOUNG S MODULUS OF WOOD By: Mehran Roohnia

FIXED-FREE AGAINST FREE-FREE BEAMS FOR DYNAMIC YOUNG S MODULUS OF WOOD By: Mehran Roohnia FIXED-FREE AGAINST FREE-FREE BEAMS FOR DYNAMIC YOUNG S MODULUS OF WOOD By: Mehra Roohia Itroductio: Itroductio: Modulus of Elasticity Hook's Law: s : Stress (e.g. Normal Stress) E: Modulus of Elasticity,

More information

2C09 Design for seismic and climate changes

2C09 Design for seismic and climate changes C9 Desig for seismic ad climate chages Lecture 3: Dyamic respose of sigle-degree-of-freedom systems II Daiel Grecea, Politehica Uiversity of Timisoara 11/3/14 Europea Erasmus Mudus Master Course Sustaiable

More information

Elastic Plastic Behavior of Geomaterials: Modeling and Simulation Issues

Elastic Plastic Behavior of Geomaterials: Modeling and Simulation Issues Elastic Plastic Behavior of Geomaterials: Modelig ad Simulatio Issues Boris Zhaohui Yag (UA), Zhao Cheg (EarthMechaics Ic.), Mahdi Taiebat (UBC) Departmet of Civil ad Evirometal Egieerig Uiversity of Califoria,

More information

Vibratory Motion. Prof. Zheng-yi Feng NCHU SWC. National CHung Hsing University, Department of Soil and Water Conservation

Vibratory Motion. Prof. Zheng-yi Feng NCHU SWC. National CHung Hsing University, Department of Soil and Water Conservation Vibratory Motio Prof. Zheg-yi Feg NCHU SWC 1 Types of vibratory motio Periodic motio Noperiodic motio See Fig. A1, p.58 Harmoic motio Periodic motio Trasiet motio impact Trasiet motio earthquake A powerful

More information

577. Estimation of surface roughness using high frequency vibrations

577. Estimation of surface roughness using high frequency vibrations 577. Estimatio of surface roughess usig high frequecy vibratios V. Augutis, M. Sauoris, Kauas Uiversity of Techology Electroics ad Measuremets Systems Departmet Studetu str. 5-443, LT-5368 Kauas, Lithuaia

More information

Natural frequencies, mode shapes and modal damping values are the most important parameters to

Natural frequencies, mode shapes and modal damping values are the most important parameters to Abstract Natural frequecies, mode shapes ad modal dampig values are the most importat parameters to describe the oise ad vibratio behaviour of a mechaical system. For rotatig machiery, however, the directivity

More information

Research Article Health Monitoring for a Structure Using Its Nonstationary Vibration

Research Article Health Monitoring for a Structure Using Its Nonstationary Vibration Advaces i Acoustics ad Vibratio Volume 2, Article ID 69652, 5 pages doi:.55/2/69652 Research Article Health Moitorig for a Structure Usig Its Nostatioary Vibratio Yoshimutsu Hirata, Mikio Tohyama, Mitsuo

More information

A PROCEDURE TO MODIFY THE FREQUENCY AND ENVELOPE CHARACTERISTICS OF EMPIRICAL GREEN'S FUNCTION. Lin LU 1 SUMMARY

A PROCEDURE TO MODIFY THE FREQUENCY AND ENVELOPE CHARACTERISTICS OF EMPIRICAL GREEN'S FUNCTION. Lin LU 1 SUMMARY A POCEDUE TO MODIFY THE FEQUENCY AND ENVELOPE CHAACTEISTICS OF EMPIICAL GEEN'S FUNCTION Li LU SUMMAY Semi-empirical method, which divides the fault plae of large earthquake ito mets ad uses small groud

More information

EXPERIMENT OF SIMPLE VIBRATION

EXPERIMENT OF SIMPLE VIBRATION EXPERIMENT OF SIMPLE VIBRATION. PURPOSE The purpose of the experimet is to show free vibratio ad damped vibratio o a system havig oe degree of freedom ad to ivestigate the relatioship betwee the basic

More information

Provläsningsexemplar / Preview TECHNICAL REPORT INTERNATIONAL SPECIAL COMMITTEE ON RADIO INTERFERENCE

Provläsningsexemplar / Preview TECHNICAL REPORT INTERNATIONAL SPECIAL COMMITTEE ON RADIO INTERFERENCE TECHNICAL REPORT CISPR 16-4-3 2004 AMENDMENT 1 2006-10 INTERNATIONAL SPECIAL COMMITTEE ON RADIO INTERFERENCE Amedmet 1 Specificatio for radio disturbace ad immuity measurig apparatus ad methods Part 4-3:

More information

Research Article Health Monitoring for a Structure Using Its Nonstationary Vibration

Research Article Health Monitoring for a Structure Using Its Nonstationary Vibration Hidawi Publishig Corporatio Advaces i Acoustics ad Vibratio Volume 2, Article ID 69652, 5 pages doi:.55/2/69652 Research Article Health Moitorig for a Structure Usig Its Nostatioary Vibratio Yoshimutsu

More information

1. Linearization of a nonlinear system given in the form of a system of ordinary differential equations

1. Linearization of a nonlinear system given in the form of a system of ordinary differential equations . Liearizatio of a oliear system give i the form of a system of ordiary differetial equatios We ow show how to determie a liear model which approximates the behavior of a time-ivariat oliear system i a

More information

School of Mechanical Engineering Purdue University. ME375 Frequency Response - 1

School of Mechanical Engineering Purdue University. ME375 Frequency Response - 1 Case Study ME375 Frequecy Respose - Case Study SUPPORT POWER WIRE DROPPERS Electric trai derives power through a patograph, which cotacts the power wire, which is suspeded from a cateary. Durig high-speed

More information

THE APPEARANCE OF FIBONACCI AND LUCAS NUMBERS IN THE SIMULATION OF ELECTRICAL POWER LINES SUPPLIED BY TWO SIDES

THE APPEARANCE OF FIBONACCI AND LUCAS NUMBERS IN THE SIMULATION OF ELECTRICAL POWER LINES SUPPLIED BY TWO SIDES THE APPEARANCE OF FIBONACCI AND LUCAS NUMBERS IN THE SIMULATION OF ELECTRICAL POWER LINES SUPPLIED BY TWO SIDES Giuseppe Ferri Dipartimeto di Ihgegeria Elettrica-FacoM di Igegeria, Uiversita di L'Aquila

More information

μ are complex parameters. Other

μ are complex parameters. Other A New Numerical Itegrator for the Solutio of Iitial Value Problems i Ordiary Differetial Equatios. J. Suday * ad M.R. Odekule Departmet of Mathematical Scieces, Adamawa State Uiversity, Mubi, Nigeria.

More information

Introduction to Signals and Systems, Part V: Lecture Summary

Introduction to Signals and Systems, Part V: Lecture Summary EEL33: Discrete-Time Sigals ad Systems Itroductio to Sigals ad Systems, Part V: Lecture Summary Itroductio to Sigals ad Systems, Part V: Lecture Summary So far we have oly looked at examples of o-recursive

More information

AN OPEN-PLUS-CLOSED-LOOP APPROACH TO SYNCHRONIZATION OF CHAOTIC AND HYPERCHAOTIC MAPS

AN OPEN-PLUS-CLOSED-LOOP APPROACH TO SYNCHRONIZATION OF CHAOTIC AND HYPERCHAOTIC MAPS http://www.paper.edu.c Iteratioal Joural of Bifurcatio ad Chaos, Vol. 1, No. 5 () 119 15 c World Scietific Publishig Compay AN OPEN-PLUS-CLOSED-LOOP APPROACH TO SYNCHRONIZATION OF CHAOTIC AND HYPERCHAOTIC

More information

CUMULATIVE DAMAGE ESTIMATION USING WAVELET TRANSFORM OF STRUCTURAL RESPONSE

CUMULATIVE DAMAGE ESTIMATION USING WAVELET TRANSFORM OF STRUCTURAL RESPONSE CUMULATIVE DAMAGE ESTIMATION USING WAVELET TRANSFORM OF STRUCTURAL RESPONSE Ryutaro SEGAWA 1, Shizuo YAMAMOTO, Akira SONE 3 Ad Arata MASUDA 4 SUMMARY Durig a strog earthquake, the respose of a structure

More information

Sample Size Estimation in the Proportional Hazards Model for K-sample or Regression Settings Scott S. Emerson, M.D., Ph.D.

Sample Size Estimation in the Proportional Hazards Model for K-sample or Regression Settings Scott S. Emerson, M.D., Ph.D. ample ie Estimatio i the Proportioal Haards Model for K-sample or Regressio ettigs cott. Emerso, M.D., Ph.D. ample ie Formula for a Normally Distributed tatistic uppose a statistic is kow to be ormally

More information

Added mass estimation of open flat membranes vibrating in still air

Added mass estimation of open flat membranes vibrating in still air Added mass estimatio of ope flat membraes vibratig i still air *Yua-qi Li 1),Yi Zhou ), Akihito YOHIDA 3) ad Yukio TAMURA 4) 1), ) Departmet of Buildig Egieerig, Togji Uiversity, haghai, 9, Chia 1) liyq@togji.edu.c

More information

Dynamic Response of Linear Systems

Dynamic Response of Linear Systems Dyamic Respose of Liear Systems Liear System Respose Superpositio Priciple Resposes to Specific Iputs Dyamic Respose of st Order Systems Characteristic Equatio - Free Respose Stable st Order System Respose

More information

Dynamic Response of Second Order Mechanical Systems with Viscous Dissipation forces

Dynamic Response of Second Order Mechanical Systems with Viscous Dissipation forces Hadout #b (pp. 4-55) Dyamic Respose o Secod Order Mechaical Systems with Viscous Dissipatio orces M X + DX + K X = F t () Periodic Forced Respose to F (t) = F o si( t) ad F (t) = M u si(t) Frequecy Respose

More information

CS284A: Representations and Algorithms in Molecular Biology

CS284A: Representations and Algorithms in Molecular Biology CS284A: Represetatios ad Algorithms i Molecular Biology Scribe Notes o Lectures 3 & 4: Motif Discovery via Eumeratio & Motif Represetatio Usig Positio Weight Matrix Joshua Gervi Based o presetatios by

More information

567. Research of Dynamics of a Vibration Isolation Platform

567. Research of Dynamics of a Vibration Isolation Platform 567. Research of Dyamics of a Vibratio Isolatio Platform A. Kilikevičius, M. Jurevičius 2, M. Berba 3 Vilius Gedimias Techical Uiversity, Departmet of Machie buildig, J. Basaavičiaus str. 28, LT-03224

More information

TR Use of the Implicit HHT-I3 and the Explicit ADAMS Methods with the Absolute Nodal Coordinate Formulation

TR Use of the Implicit HHT-I3 and the Explicit ADAMS Methods with the Absolute Nodal Coordinate Formulation R-007-05 Use of the Implicit HH-I3 ad the Explicit ADAMS Methods with the Absolute Nodal Coordiate Formulatio Bassam Hussei, Da Negrut, Ahmed A. Shabaa August 007 Abstract his ivestigatio is cocered with

More information

Lecture 2: Monte Carlo Simulation

Lecture 2: Monte Carlo Simulation STAT/Q SCI 43: Itroductio to Resamplig ethods Sprig 27 Istructor: Ye-Chi Che Lecture 2: ote Carlo Simulatio 2 ote Carlo Itegratio Assume we wat to evaluate the followig itegratio: e x3 dx What ca we do?

More information

PROPOSING INPUT-DEPENDENT MODE CONTRIBUTION FACTORS FOR SIMPLIFIED SEISMIC RESPONSE ANALYSIS OF BUILDING SYSTEMS

PROPOSING INPUT-DEPENDENT MODE CONTRIBUTION FACTORS FOR SIMPLIFIED SEISMIC RESPONSE ANALYSIS OF BUILDING SYSTEMS he 4 th World Coferece o Earthquake Egieerig October -7, 008, Beiig, Chia PROPOSING INPU-DEPENDEN ODE CONRIBUION FACORS FOR SIPLIFIED SEISIC RESPONSE ANALYSIS OF BUILDING SYSES ahmood Hosseii ad Laya Abbasi

More information

DYNAMIC ANALYSIS OF BEAM-LIKE STRUCTURES SUBJECT TO MOVING LOADS

DYNAMIC ANALYSIS OF BEAM-LIKE STRUCTURES SUBJECT TO MOVING LOADS DYNAMIC ANALYSIS OF BEAM-LIKE STRUCTURES SUBJECT TO MOVING LOADS Ivaa Štimac 1, Ivica Kožar 1 M.Sc,Assistat, Ph.D. Professor 1, Faculty of Civil Egieerig, Uiverity of Rieka, Croatia INTRODUCTION The vehicle-iduced

More information

Correlation. Two variables: Which test? Relationship Between Two Numerical Variables. Two variables: Which test? Contingency table Grouped bar graph

Correlation. Two variables: Which test? Relationship Between Two Numerical Variables. Two variables: Which test? Contingency table Grouped bar graph Correlatio Y Two variables: Which test? X Explaatory variable Respose variable Categorical Numerical Categorical Cotigecy table Cotigecy Logistic Grouped bar graph aalysis regressio Mosaic plot Numerical

More information

DETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS.

DETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS. ICSV4 Cairs Australia 9- July 7 DTRMINATION OF MCHANICAL PROPRTIS OF A NON- UNIFORM BAM USING TH MASURMNT OF TH XCITD LONGITUDINAL LASTIC VIBRATIONS Pavel Aokhi ad Vladimir Gordo Departmet of the mathematics

More information

EECE 301 Signals & Systems

EECE 301 Signals & Systems EECE 301 Sigals & Systems Prof. Mark Fowler Note Set #8 D-T Covolutio: The Tool for Fidig the Zero-State Respose Readig Assigmet: Sectio 2.1-2.2 of Kame ad Heck 1/14 Course Flow Diagram The arrows here

More information

(all terms are scalars).the minimization is clearer in sum notation:

(all terms are scalars).the minimization is clearer in sum notation: 7 Multiple liear regressio: with predictors) Depedet data set: y i i = 1, oe predictad, predictors x i,k i = 1,, k = 1, ' The forecast equatio is ŷ i = b + Use matrix otatio: k =1 b k x ik Y = y 1 y 1

More information

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense, 3. Z Trasform Referece: Etire Chapter 3 of text. Recall that the Fourier trasform (FT) of a DT sigal x [ ] is ω ( ) [ ] X e = j jω k = xe I order for the FT to exist i the fiite magitude sese, S = x [

More information

Application of an Innovative Precise Integration Method in Solving Equilibrium Equations of Motion for Structural Dynamic Problems

Application of an Innovative Precise Integration Method in Solving Equilibrium Equations of Motion for Structural Dynamic Problems Applicatio of a Iovative Precise Itegratio Method i Solvig Equilibrium Equatios of Motio for Structural Dyamic Problems Chiu-li Wu Natioal Ceter for Research o Earthquake Egieerig,Taipei City, Taiwa Chig-Chiag

More information

Chapter 6 Sampling Distributions

Chapter 6 Sampling Distributions Chapter 6 Samplig Distributios 1 I most experimets, we have more tha oe measuremet for ay give variable, each measuremet beig associated with oe radomly selected a member of a populatio. Hece we eed to

More information

Nonlinear Vibrations of Aerospace Structures

Nonlinear Vibrations of Aerospace Structures Noliear Vibratios of Aerospace Structures Uiversity of Liège, Belgium L Itroductio Course objectives Review of liear theory Istructors: G. Kersche, J.P. Noel, T. Detroux Cotact details: Space Structures

More information

Proceedings of the 2nd International Conference on Computer Science and Electronics Engineering (ICCSEE 2013)

Proceedings of the 2nd International Conference on Computer Science and Electronics Engineering (ICCSEE 2013) Proceedigs of the 2d Iteratioal Coferece o Computer Sciece ad Electroics Egieerig (ICCE 20) GPU Implemetatio of the Electric Power System Model for Real-Time Simulatio of Electromagetic Trasiets Slawomir

More information

Microscopic traffic flow modeling

Microscopic traffic flow modeling Chapter 34 Microscopic traffic flow modelig 34.1 Overview Macroscopic modelig looks at traffic flow from a global perspective, whereas microscopic modelig, as the term suggests, gives attetio to the details

More information

Exponential Moving Average Pieter P

Exponential Moving Average Pieter P Expoetial Movig Average Pieter P Differece equatio The Differece equatio of a expoetial movig average lter is very simple: y[] x[] + (1 )y[ 1] I this equatio, y[] is the curret output, y[ 1] is the previous

More information

Computational Fluid Dynamics. Lecture 5

Computational Fluid Dynamics. Lecture 5 Time differecig cotiued. Three level schemes. B. Modified L-F schemes. C. Higher order methods. Three level schemes Computatioal Fluid Dyamics Lecture 5 ψ = α α β tf β cosistet schemes if α α = ad β β

More information

Difference Equation Construction (1) ENGG 1203 Tutorial. Difference Equation Construction (2) Grow, baby, grow (1)

Difference Equation Construction (1) ENGG 1203 Tutorial. Difference Equation Construction (2) Grow, baby, grow (1) ENGG 03 Tutorial Differece Equatio Costructio () Systems ad Cotrol April Learig Objectives Differece Equatios Z-trasform Poles Ack.: MIT OCW 6.0, 6.003 Newto s law of coolig states that: The chage i a

More information

Solution of Differential Equation from the Transform Technique

Solution of Differential Equation from the Transform Technique Iteratioal Joural of Computatioal Sciece ad Mathematics ISSN 0974-3189 Volume 3, Number 1 (2011), pp 121-125 Iteratioal Research Publicatio House http://wwwirphousecom Solutio of Differetial Equatio from

More information

System Diagnostics using Kalman Filter Estimation Error

System Diagnostics using Kalman Filter Estimation Error System Diagostics usig Kalma Filter Estimatio Error Prof. Seugchul Lee, Seugtae Park, Heechag Kim, Hyusuk Huh ICMR2015 (11/25/2015) Machie Health Diagostics Desig Desig DNA + + Blue Prit phagocytes lymphocytes

More information

REAL-TIME MONITORING OF POWER AND NEUTRON CAPTURE CROSS SECTION IN NUCLEAR RESEARCH REACTOR

REAL-TIME MONITORING OF POWER AND NEUTRON CAPTURE CROSS SECTION IN NUCLEAR RESEARCH REACTOR D5 REAL-TIME MONITORING OF POWER AND NEUTRON CAPTURE CROSS SECTION IN NUCLEAR RESEARCH REACTOR D.A.P. PALMA Brazilia Nuclear Eergy Commissio (CNEN), Rio de Jaeiro, Brazil A.Z. MESQUITA, R.M.G.P. SOUZA

More information

Dynamic System Response

Dynamic System Response Solutio of Liear, Costat-Coefficiet, Ordiary Differetial Equatios Classical Operator Method Laplace Trasform Method Laplace Trasform Properties 1 st -Order Dyamic System Time ad Frequecy Respose d -Order

More information

FREE VIBRATION RESPONSE OF A SYSTEM WITH COULOMB DAMPING

FREE VIBRATION RESPONSE OF A SYSTEM WITH COULOMB DAMPING Mechaical Vibratios FREE VIBRATION RESPONSE OF A SYSTEM WITH COULOMB DAMPING A commo dampig mechaism occurrig i machies is caused by slidig frictio or dry frictio ad is called Coulomb dampig. Coulomb dampig

More information

Finite Element Analysis of Cable-Stayed Bridges with Appropriate Initial Shapes Under Seismic Excitations Focusing on Deck-Stay Interaction

Finite Element Analysis of Cable-Stayed Bridges with Appropriate Initial Shapes Under Seismic Excitations Focusing on Deck-Stay Interaction Chapter 9 Fiite Elemet Aalysis of Cable-Stayed Bridges with Appropriate Iitial Shapes Uder Seismic Excitatios Focusig o Deck-Stay Iteractio Mig-Yi Liu ad Pao-Hsii Wag Additioal iformatio is available at

More information

Vibration Absorbers/Neutralisers. Professor Mike Brennan

Vibration Absorbers/Neutralisers. Professor Mike Brennan Vibratio Absorbers/Neutralisers Professor Mike Brea Vibratio absorbers/eutralisers m k c Vibratio absorber/eutraliser What does it do? Absorber or eutraliser Mobility approach effectiveess of the device

More information

APPLICATION OF NETSLAB FOR REMOTE TESTING OF BRIDGE SYSTEMS ABSTRACT

APPLICATION OF NETSLAB FOR REMOTE TESTING OF BRIDGE SYSTEMS ABSTRACT APPLICATION OF NETSLAB FOR REMOTE TESTING OF BRIDGE SYSTEMS Yulei Fa, Ya Xiao, ad Yurog Guo ABSTRACT A Iteret based remote hybrid testig platform was developed for bridge systems. The testig platform,

More information

CALCULATION OF STIFFNESS AND MASS ORTHOGONAL VECTORS

CALCULATION OF STIFFNESS AND MASS ORTHOGONAL VECTORS 14. CALCULATION OF STIFFNESS AND MASS ORTHOGONAL VECTORS LDR Vectors are Always More Accurate tha Usig the Exact Eigevectors i a Mode Superpositio Aalysis 14.1 INTRODUCTION The major reaso to calculate

More information

Microscopic Theory of Transport (Fall 2003) Lecture 6 (9/19/03) Static and Short Time Properties of Time Correlation Functions

Microscopic Theory of Transport (Fall 2003) Lecture 6 (9/19/03) Static and Short Time Properties of Time Correlation Functions .03 Microscopic Theory of Trasport (Fall 003) Lecture 6 (9/9/03) Static ad Short Time Properties of Time Correlatio Fuctios Refereces -- Boo ad Yip, Chap There are a umber of properties of time correlatio

More information

MOST PEOPLE WOULD RATHER LIVE WITH A PROBLEM THEY CAN'T SOLVE, THAN ACCEPT A SOLUTION THEY CAN'T UNDERSTAND.

MOST PEOPLE WOULD RATHER LIVE WITH A PROBLEM THEY CAN'T SOLVE, THAN ACCEPT A SOLUTION THEY CAN'T UNDERSTAND. XI-1 (1074) MOST PEOPLE WOULD RATHER LIVE WITH A PROBLEM THEY CAN'T SOLVE, THAN ACCEPT A SOLUTION THEY CAN'T UNDERSTAND. R. E. D. WOOLSEY AND H. S. SWANSON XI-2 (1075) STATISTICAL DECISION MAKING Advaced

More information

Complex Algorithms for Lattice Adaptive IIR Notch Filter

Complex Algorithms for Lattice Adaptive IIR Notch Filter 4th Iteratioal Coferece o Sigal Processig Systems (ICSPS ) IPCSIT vol. 58 () () IACSIT Press, Sigapore DOI:.7763/IPCSIT..V58. Complex Algorithms for Lattice Adaptive IIR Notch Filter Hog Liag +, Nig Jia

More information

Load Dependent Ritz Vector Algorithm and Error Analysis

Load Dependent Ritz Vector Algorithm and Error Analysis Load Depedet Ritz Vector Algorithm ad Error Aalysis Writte by Ed Wilso i 006. he Complete eigealue subspace I the aalysis of structures subected to three base acceleratios there is a requiremet that oe

More information

The Phi Power Series

The Phi Power Series The Phi Power Series I did this work i about 0 years while poderig the relatioship betwee the golde mea ad the Madelbrot set. I have fially decided to make it available from my blog at http://semresearch.wordpress.com/.

More information

Finite Element Modeling of Seismic Response of Field Fabricated Liquefied Natural Gas (LNG) Spherical Storage Vessels

Finite Element Modeling of Seismic Response of Field Fabricated Liquefied Natural Gas (LNG) Spherical Storage Vessels Egieerig, 013, 5, 543-550 doi:10.436/eg.013.56065 Published Olie Jue 013 (http://www.scirp.org/joural/eg) Fiite Elemet Modelig of Seismic Respose of Field Fabricated Liquefied Natural Gas (LNG) Spherical

More information

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would

More information

Application of Dynamic Relaxation in Thermo-Elastic Structural Analysis of Highway Pavement Structures

Application of Dynamic Relaxation in Thermo-Elastic Structural Analysis of Highway Pavement Structures 9 th Iteratioal LS-DYNA Users Coferece Simulatio Techology () Applicatio of Dyamic Relaxatio i Thermo-Elastic Structural Aalysis of Highway Pavemet Structures Samir N. Shoukry, Gergis W. William, Mourad

More information

Principle Of Superposition

Principle Of Superposition ecture 5: PREIMINRY CONCEP O RUCUR NYI Priciple Of uperpositio Mathematically, the priciple of superpositio is stated as ( a ) G( a ) G( ) G a a or for a liear structural system, the respose at a give

More information

Dynamic Programming. Sequence Of Decisions

Dynamic Programming. Sequence Of Decisions Dyamic Programmig Sequece of decisios. Problem state. Priciple of optimality. Dyamic Programmig Recurrece Equatios. Solutio of recurrece equatios. Sequece Of Decisios As i the greedy method, the solutio

More information

Dynamic Programming. Sequence Of Decisions. 0/1 Knapsack Problem. Sequence Of Decisions

Dynamic Programming. Sequence Of Decisions. 0/1 Knapsack Problem. Sequence Of Decisions Dyamic Programmig Sequece Of Decisios Sequece of decisios. Problem state. Priciple of optimality. Dyamic Programmig Recurrece Equatios. Solutio of recurrece equatios. As i the greedy method, the solutio

More information

ADVANCED DIGITAL SIGNAL PROCESSING

ADVANCED DIGITAL SIGNAL PROCESSING ADVANCED DIGITAL SIGNAL PROCESSING PROF. S. C. CHAN (email : sccha@eee.hku.hk, Rm. CYC-702) DISCRETE-TIME SIGNALS AND SYSTEMS MULTI-DIMENSIONAL SIGNALS AND SYSTEMS RANDOM PROCESSES AND APPLICATIONS ADAPTIVE

More information

Modification of Arrhenius Model for Numerical Modelling of Turbulent Flames

Modification of Arrhenius Model for Numerical Modelling of Turbulent Flames J. Basic. Appl. Sci. Res., (5)480-486, 01 01, TextRoad Publicatio ISSN 090-4304 Joural of Basic ad Applied Scietific Research www.textroad.com Modificatio of Arrheius Model for Numerical Modellig of Turbulet

More information

A STUDY OF VIBRATION MEASURING AND FATIGUE ANALYSIS FOR CANTILEVER BEAMS

A STUDY OF VIBRATION MEASURING AND FATIGUE ANALYSIS FOR CANTILEVER BEAMS Joural of Techology, Vol. 3, No., pp. 47-56 (07) 47, * LabView ANSYS A STUDY OF VIBRATION MEASURING AND FATIGUE ANALYSIS FOR CANTILEVER BEAMS Yuo-Ter Tsai, * Hsie-Yag Li Departmet of Mechaical Egieerig

More information

THE ENERGY BALANCE ERROR FOR CIRCUIT TRANSIENT ANALYSIS

THE ENERGY BALANCE ERROR FOR CIRCUIT TRANSIENT ANALYSIS THE ENERGY BALANCE ERROR FOR CIRCUIT TRANSIENT ANALYSIS FLORIN CONSTANTINESCU, ALEXANDRU GABRIEL GHEORGHE, MIRUNA NIŢESCU Key words: Trasiet aalysis, Eergy balace error, Time step coice. Two algoritms

More information

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would

More information

Analysis of Algorithms. Introduction. Contents

Analysis of Algorithms. Introduction. Contents Itroductio The focus of this module is mathematical aspects of algorithms. Our mai focus is aalysis of algorithms, which meas evaluatig efficiecy of algorithms by aalytical ad mathematical methods. We

More information

Chapter 5 Vibrational Motion

Chapter 5 Vibrational Motion Fall 4 Chapter 5 Vibratioal Motio... 65 Potetial Eergy Surfaces, Rotatios ad Vibratios... 65 Harmoic Oscillator... 67 Geeral Solutio for H.O.: Operator Techique... 68 Vibratioal Selectio Rules... 7 Polyatomic

More information

8. СОВЕТУВАЊЕ. Охрид, септември ANALYSIS OF NO LOAD APPARENT POWER AND FREQUENCY SPECTRUM OF MAGNETIZING CURRENT FOR DIFFERENT CORE TYPES

8. СОВЕТУВАЊЕ. Охрид, септември ANALYSIS OF NO LOAD APPARENT POWER AND FREQUENCY SPECTRUM OF MAGNETIZING CURRENT FOR DIFFERENT CORE TYPES 8. СОВЕТУВАЊЕ Охрид, 22 24 септември Leoardo Štrac Frajo Keleme Kočar Power Trasformers Ltd. ANALYSS OF NO LOAD APPARENT POWER AND FREQENCY SPECTRM OF MAGNETZNG CRRENT FOR DFFERENT CORE TYPES ABSTRACT

More information

STA6938-Logistic Regression Model

STA6938-Logistic Regression Model Dr. Yig Zhag STA6938-Logistic Regressio Model Topic -Simple (Uivariate) Logistic Regressio Model Outlies:. Itroductio. A Example-Does the liear regressio model always work? 3. Maximum Likelihood Curve

More information

Orthogonal Gaussian Filters for Signal Processing

Orthogonal Gaussian Filters for Signal Processing Orthogoal Gaussia Filters for Sigal Processig Mark Mackezie ad Kiet Tieu Mechaical Egieerig Uiversity of Wollogog.S.W. Australia Abstract A Gaussia filter usig the Hermite orthoormal series of fuctios

More information

New Exponential Strengthening Buffer Operators and Numerical Simulation

New Exponential Strengthening Buffer Operators and Numerical Simulation Sesors & Trasducers, Vol. 59, Issue, November 0, pp. 7-76 Sesors & Trasducers 0 by IFSA http://www.sesorsportal.com New Expoetial Stregtheig Buffer Operators ad Numerical Simulatio Cuifeg Li, Huajie Ye,

More information

EECS 442 Computer vision. Multiple view geometry Affine structure from Motion

EECS 442 Computer vision. Multiple view geometry Affine structure from Motion EECS 442 Computer visio Multiple view geometry Affie structure from Motio - Affie structure from motio problem - Algebraic methods - Factorizatio methods Readig: [HZ] Chapters: 6,4,8 [FP] Chapter: 2 Some

More information

Analysis of Foundation-Layered Soil Interaction Using Propagating Wave-modes Analysis

Analysis of Foundation-Layered Soil Interaction Using Propagating Wave-modes Analysis Aalysis of Foudatio-Layered Soil Iteractio Usig Propagatig Wave-modes Aalysis IGOR ŠPACAPAN & MIROSLAV PREMROV Faculty of Civil Egieerig Uiversity of Maribor Smetaova 7, Maribor SLOVENIA Abstract: - I

More information

Visco-acoustic Reverse-time Migration Using the Pseudo-analytical Method

Visco-acoustic Reverse-time Migration Using the Pseudo-analytical Method Visco-acoustic Reverse-time Migratio Usig the Pseudo-aalytical Method J. Ramos-Martiez* (Petroleum Geo-Services), A.A. Valeciao (Petroleum Geo-Services) & M. Orlovich (Petroleum Geo-Services) SUMMARY We

More information

A Simplified Derivation of Scalar Kalman Filter using Bayesian Probability Theory

A Simplified Derivation of Scalar Kalman Filter using Bayesian Probability Theory 6th Iteratioal Workshop o Aalysis of Dyamic Measuremets Jue -3 0 Göteorg Swede A Simplified Derivatio of Scalar Kalma Filter usig Bayesia Proaility Theory Gregory Kyriazis Imetro - Electrical Metrology

More information

DECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS. Park Road, Islamabad, Pakistan

DECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS. Park Road, Islamabad, Pakistan Mathematical ad Computatioal Applicatios, Vol. 9, No. 3, pp. 30-40, 04 DECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS Muhammad Aslam Noor, Khalida Iayat Noor ad Asif Waheed

More information

Chapter 12 EM algorithms The Expectation-Maximization (EM) algorithm is a maximum likelihood method for models that have hidden variables eg. Gaussian

Chapter 12 EM algorithms The Expectation-Maximization (EM) algorithm is a maximum likelihood method for models that have hidden variables eg. Gaussian Chapter 2 EM algorithms The Expectatio-Maximizatio (EM) algorithm is a maximum likelihood method for models that have hidde variables eg. Gaussia Mixture Models (GMMs), Liear Dyamic Systems (LDSs) ad Hidde

More information

BIOSTATISTICS. Lecture 5 Interval Estimations for Mean and Proportion. dr. Petr Nazarov

BIOSTATISTICS. Lecture 5 Interval Estimations for Mean and Proportion. dr. Petr Nazarov Microarray Ceter BIOSTATISTICS Lecture 5 Iterval Estimatios for Mea ad Proportio dr. Petr Nazarov 15-03-013 petr.azarov@crp-sate.lu Lecture 5. Iterval estimatio for mea ad proportio OUTLINE Iterval estimatios

More information

TESTING OF THE FORCES IN CABLE OF SUSPENSION STRUCTURE AND BRIDGES

TESTING OF THE FORCES IN CABLE OF SUSPENSION STRUCTURE AND BRIDGES TSTING OF TH FORCS IN CABL OF SUSPNSION STRUCTUR AND BRIDGS Zhou, M. 1, Liu, Z. ad Liu, J. 1 College of the Muicipal Techology, Guagzhou Uiversity, Guagzhou. Guagzhou Muicipal ad Ladscape gieerig Quality

More information

Numerical Solution of the Two Point Boundary Value Problems By Using Wavelet Bases of Hermite Cubic Spline Wavelets

Numerical Solution of the Two Point Boundary Value Problems By Using Wavelet Bases of Hermite Cubic Spline Wavelets Australia Joural of Basic ad Applied Scieces, 5(): 98-5, ISSN 99-878 Numerical Solutio of the Two Poit Boudary Value Problems By Usig Wavelet Bases of Hermite Cubic Splie Wavelets Mehdi Yousefi, Hesam-Aldie

More information

Using SAS to Evaluate Integrals and Reverse Functions in Power and Sample Size Calculations

Using SAS to Evaluate Integrals and Reverse Functions in Power and Sample Size Calculations Usig SAS to Evaluate Itegrals ad everse Fuctios i Power ad Sample Size Calculatios Xigshu Zhu, Merck &Co., Ic. Blue Bell, PA 94 Shupig Zhag, Merck &Co., Ic. Blue Bell, PA 94 ABSTACT We have recetly created

More information

Bayesian Control Charts for the Two-parameter Exponential Distribution

Bayesian Control Charts for the Two-parameter Exponential Distribution Bayesia Cotrol Charts for the Two-parameter Expoetial Distributio R. va Zyl, A.J. va der Merwe 2 Quitiles Iteratioal, ruaavz@gmail.com 2 Uiversity of the Free State Abstract By usig data that are the mileages

More information

Harmonic balance approaches to the nonlinear oscillators in which the restoring force is inversely proportional to the dependent variable

Harmonic balance approaches to the nonlinear oscillators in which the restoring force is inversely proportional to the dependent variable Joural of Soud ad Vibratio, Vol. 314, Nº 3-5, 775-78 (008) Harmoic balace approaches to the oliear oscillators i which the restorig force is iversely proportioal to the depedet variable A. Belédez, D.

More information

Stability Analysis and Bifurcation Control of Hysteresis Current Controlled Ćuk Converter Using Filippov s Method

Stability Analysis and Bifurcation Control of Hysteresis Current Controlled Ćuk Converter Using Filippov s Method Stability Aalysis ad Bifurcatio otrol of Hysteresis urret otrolled Ću overter Usig Filippov s Method I. Daho*, D. Giaouris * (Member IE, B. Zahawi*(Member IE, V. Picer*(Member IE ad S. Baerjee** *School

More information

Signal Processing in Mechatronics

Signal Processing in Mechatronics Sigal Processig i Mechatroics Zhu K.P. AIS, UM. Lecture, Brief itroductio to Sigals ad Systems, Review of Liear Algebra ad Sigal Processig Related Mathematics . Brief Itroductio to Sigals What is sigal

More information

Some New Iterative Methods for Solving Nonlinear Equations

Some New Iterative Methods for Solving Nonlinear Equations World Applied Scieces Joural 0 (6): 870-874, 01 ISSN 1818-495 IDOSI Publicatios, 01 DOI: 10.589/idosi.wasj.01.0.06.830 Some New Iterative Methods for Solvig Noliear Equatios Muhammad Aslam Noor, Khalida

More information

LOSS-MINIMIZATION CONTROL OF SCALAR- CONTROLLED INDUCTION MOTOR DRIVES

LOSS-MINIMIZATION CONTROL OF SCALAR- CONTROLLED INDUCTION MOTOR DRIVES LOSS-MINIMIZATION CONTROL OF SCALAR- CONTROLLED INDUCTION MOTOR DRIVES Hussei Sarha, Rateb Al-Issa, ad Qazem Jaber Departmet of Mechatroics Egieerig, Faculty of Egieerig Techology Al-Balqa Applied Uiversity,

More information

Addition: Property Name Property Description Examples. a+b = b+a. a+(b+c) = (a+b)+c

Addition: Property Name Property Description Examples. a+b = b+a. a+(b+c) = (a+b)+c Notes for March 31 Fields: A field is a set of umbers with two (biary) operatios (usually called additio [+] ad multiplicatio [ ]) such that the followig properties hold: Additio: Name Descriptio Commutativity

More information

POWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS

POWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS Joural of Applied Mathematics ad Computatioal Mechaics 4 3(3) 3-8 POWER SERIES SOLUION OF FIRS ORDER MARIX DIFFERENIAL EQUAIONS Staisław Kukla Izabela Zamorska Istitute of Mathematics Czestochowa Uiversity

More information

ASYMPTOTIC CLOSED-LOOP DESIGN OF ERROR RESILIENT PREDICTIVE COMPRESSION SYSTEMS. Sina Zamani, Tejaswi Nanjundaswamy, Kenneth Rose

ASYMPTOTIC CLOSED-LOOP DESIGN OF ERROR RESILIENT PREDICTIVE COMPRESSION SYSTEMS. Sina Zamani, Tejaswi Nanjundaswamy, Kenneth Rose ASYMPTOTIC OSED-LOOP DESIGN OF ERROR RESILIENT PREDICTIVE COMPRESSION SYSTEMS Sia Zamai, Tejaswi Najudaswamy, Keeth Rose Departmet of Electrical ad Computer Egieerig, Uiversity of Califoria Sata Barbara,

More information

Damped Vibration of a Non-prismatic Beam with a Rotational Spring

Damped Vibration of a Non-prismatic Beam with a Rotational Spring Vibratios i Physical Systems Vol.6 (0) Damped Vibratio of a No-prismatic Beam with a Rotatioal Sprig Wojciech SOCHACK stitute of Mechaics ad Fudametals of Machiery Desig Uiversity of Techology, Czestochowa,

More information

Notes on iteration and Newton s method. Iteration

Notes on iteration and Newton s method. Iteration Notes o iteratio ad Newto s method Iteratio Iteratio meas doig somethig over ad over. I our cotet, a iteratio is a sequece of umbers, vectors, fuctios, etc. geerated by a iteratio rule of the type 1 f

More information

Modeling of Finite Gradient Elasticity Using Optimal Transportation Method

Modeling of Finite Gradient Elasticity Using Optimal Transportation Method 9 th Iteratioal Cogress of Croatia Society of Mechaics 18-22 September 2018 Split, Croatia Modelig of Fiite Gradiet Elasticity Usig Optimal Trasportatio Method Boris JALUŠIĆ *, Christia WEIßENFELS +, Jurica

More information