Coverage: through entire course material (updated July 18, 2017) Terminology/miscellaneous. for a linear instrument:
|
|
- Gwendolyn Chapman
- 5 years ago
- Views:
Transcription
1 Boston University ME 310: Summer 017 Course Notes Prof. C. Farny DISCLAIMER: This document is meant to serve as a reference and an overview, and not an exclusive study guide for the course. Coverage: through entire course material (updated July 18, 017) Terminology/miscellaneous Precision Accuracy Resolution Full Scale Input, Output Significant figures Static sensitivity for a linear instrument: K = f(x) x K = output input where f(x) is the output (the functional response of the instrument) and x is the input parameter. Probability density function p(x) p(x) = 1 σ π exp [ ] (x x) σ The absolute probability percentage P is defined over a defined interval of deviation ±δx about the average value x: P ( x δx x x + δx) = x+δx x δx (1) () (3) p(x)dx (4) Convert to a simplified scale by defining β, the normal variate for all samples x, and z 1, the normal variate for a particular sample value x 1 : β = x x σ ; z 1 = x 1 x ; σ dx = σdβ (5) and substituting into Eqn 4, P ( z 1 β z 1 ) = 1 π z1 z 1 exp( β /) dβ (6) or the more common version to reflect the symmetrical distribution assumption (note change in the integration bounds and constant term): P ( z 1 β z 1 ) = π z1 1 0 exp( β /) dβ (7)
2 Analog-Digital Conversion (ADC) The ADC process discretizes the analog signal at discrete time points (sampling frequency) and discrete voltage amplitudes. Amplitude discretization: The resulting amplitude discretization corresponds to n = V F SR V res (8) where n is the number of bits (an integer value, commonly even). The full scale range must be greater than the amplitude of the analog signal in order to properly resolve the signal. Temporal discretization: The ADC can only sample incoming data at a certain rate. This rate is generally adjustable, and should be set according to the Nyquist frequency. The Nyquist sampling criterium states that at least two data points are required to in order to resolve the highest frequency present. This results in: f s f max (9) (10) f Nyq = f s = 1 (11) t s Aliasing: Effect occurs when f Nyquist is lower than one or more frequencies in the analog signal. The energy from these higher analog frequencies appears as energy at a frequency lower than f Nyquist. ADC Uncertainty: Due to its discrete nature, there are two types of uncertainty related to ADC, both of which act like resolution-based bias sources: Temporal: The digitizer is sampling the signal at discrete time points, so the amplitude at any given time has an inherent uncertainty: U B,t = ± t s Magnitude resolution: The vertical resolution due to the number of bits and full scale range introduces a resolution error: (1) U B,V = ± V res (13)
3 Uncertainty analysis Precision Uncertainty: Inherent scatter in a measurable x about a mean x due to uncontrolled fluctuations. 1. Infinite Statistics: Large sample sets (sample size N > 63) σ x U P,x = ±z 1 ; (14) N z 1 = x 1 x ; (15) σ (xi x) σ x =. Finite Statistics: Small sample sets. Degrees of freedom ν = N-1 N (16) S x U P,x = ±t ν,95% ; (17) N (xi x) S x = (18) ν Note: The subscript P stands for precision (not pressure!). If the measurement were for a pressure, it would be denoted U P,p and if it were for a voltage, it would be denoted U P,V. Bias Uncertainty: Systematic uncertainties from instrumentation. Affects every data point in a data set of measurable x. Common bias sources: Resolution: Arises from rounding error. U B,resolution,x = ±resolution/ (19) Hysteresis: Arises from discrepancy in output due to direction-dependent loading effects. Commonly found in elastic systems. U B,hysteresis,x = ±%FS (0) Linearity (also referred to as Nonlinearity): Arises from instrument response, when output deviates from a linear response over the full scale range. U B,linearity,x = ±%FS (1) Sensitivity (also referred to as a Calibration Error): Arises from drift in sensitivity over time. Total bias uncertainty: U B,sensitivity,x = ±%FS () U B,x,total = (U B,x,i ) (3) 3
4 Total Uncertainty: Combine bias and precision sources for each measurable x using root-sum-square (RSS) method U total,x = (U P,x,i ) + (U B,x,i ) = U P,x,total + U B,x,total (4) Class 5 Error propagation: Combining uncertainty from multiple measurables x, y, z, where f(...) represents the generic analytical expression that relates x, y, z to the resultant parameter R: U total,r = U total,r = f(u x, U y, U z ) (5) ( f ) ( ) ( ) f f x U x + y U y + z U z (6) where each individual measurable uncertainty (U x, etc) is the total uncertainty for that parameter. Error propagation can be used in a single-instrument setting for the purpose of converting its input and output units. Measurement analysis Input/output impedance: Recall the voltage divider scenario that the function generator and oscilloscope describe: [ ] R L V L = V S (7) R S + R L where V is the generic instrument voltage, R is the generic instrument impedance, L represents the load (the downstream measurement instrument, or the oscilloscope), and S represents the source (the signal source instrument, or function generator). Linear system modeling Given 3 linear instruments 1,, and 3 in line with each other, an instrument-specific static sensitivity K i and static offset O i, and an input parameter x, the instrument I output is y, instrument II output is z and instrument III output is W. In generic functional response terms: y = f 1 (x); z = f (y); W = f 3 (z) = g(x); (8) Plugging in for the instrument-specific sensitivities: W = f 3 (f (f 1 (x))) = g(x); (9) y = K 1 x + O 1 ; z = K y + O ; W = K 3 z + O 3 ; (30) W = K sys x + O sys ; (31) where K sys = K 1 K K 3, O sys = O 3 + K 3 O + K K 3 O 1 (3) 4
5 Resistance measurement methods RTDs Strain gauge/wheatstone bridge Bridge constant Compensation techniques Linear System Calibration Linear Regression Analysis: The method of least squares error gives the following relationship for a linear system (with additional terms for an appropriately-nonlinear relationship): y fit = f(x) = a 0 + a 1 x (33) x a 0 = i ym,i x i xi y m,i ; a 1 = N x i y m,i x i ym,i (34) B B B = N x i ( xi ) (35) where x i is each control point (input) value, y m,i is each measured response (output) point value, y fit is the calibration value for a given input x, and N is the number of unique control points involved in the calibration. Note that it is helpful at times to invert the functional relationship to instead solve for x in terms of a measurement of y: x = f 1 (y) = y a 0 a 1 (36) For an exercise involving multiple output measurements per control point, a weighted scheme is used: wi x i wi y m,i w i x i wi x i y m,i a 0 = (37) B w wi wi x i y m,i w i x i wi y m,i a 1 = (38) B w B w = w i wi x i ( wi x i ) (39) w i = 1 U i where U i is the precision uncertainty value per control point. Calibration uncertainty: Potentially the most confusing of the uncertainty calculations. The goal is to understand the relationship between an output y (likely a voltage) and a physical input parameter x for an uncalibrated instrument A. Remember that A is being used to measure x, so x is unknown, and y is simply a voltage (and thus not helpful on its own). The functional response of A may only be known via use of a calibrated 5 (40)
6 instrument C that is capable of knowing what the input x actually is. C tells us what x is, and a regression analysis can then reveal how to get the constants in Eq 34 for instrument A. It gets confusing because both instruments can have bias uncertainties, and both may have their own precision uncertainties. We re dealing with multiple units (for x and for y), so a final, total uncertainty value can only be obtained via propagation of the multiple uncertainties to a final uncertainty in units of x or of y. Try to clearly differentiate which uncertainties belong with each instrument; keeping careful track of units helps a good deal. Here U meas refers to the uncertainties from the uncalibrated instrument A and U cal refers to those from the calibrated instrument C. Obtaining a total uncertainty from C is familiar and straightforward: U cal,x = ± ( U B,x,cal,i) + U p,x,cal,max (41) Remember that these will all be in units of x. Here the precision uncertainty would come from the scenario where you ve taken multiple measurements of x i at each control point value, for the purpose of accounting for scatter in the measurement. This would give rise to N values of U p,x,cal ; since only one such value would be considered for the uncertainty, you should take the maximum such value over the N that you ve calculated. For the uncalibrated instrument A, the various uncertainty parameters will likely be in units of y, or voltage. The new component of the analysis here involves the precision uncertainty of instrument A. The precision uncertainty is now a scatter of the data about a linear trend (the calibration curve), known as the precision uncertainty of the fit U p,fit that relies on a standard error of the fit S yx : U p,y,fit = ±t ν,95% Syx (4) N (yfit,i y m,i ) S yx = (43) ν The number of statistical degrees of freedom in this case (specifically, calibration uncertainty analysis) ν is N (I sometimes notate ν here as ν fit as a reminder of the calibration aspect), and N is the number of control points on the x axis. The parameter y fit,i represents the known/calibrated output value at each control point. The uncertainty from instrument A may be then found per the usual method: U meas,y = ± ( U B,y,i) + U p,y,fit (44) Since the values of the calibration for instrument A depend on the inherent uncertainties of instrument C, the two uncertainties U meas and U cal must be combined, to obtain a total uncertainty. Error propagation must be used to properly combine the different units of x and y: ( f ) U total,y = ± x U cal,x + Umeas,y (45) 6
7 ( f 1 or: U total,x = ± y U meas,y ) + U cal,x (46) where f represents the calibration equation y = a 1 x + a 0 and f 1 is the inverse of that equation, or x = y a 0 a 1. Signal conditioning Signal-to-Noise Ratio: SNR = 0 log 10 ( Vsig V noise ) Noise sources: Electrical (broadband), frequency-specific (60 Hz power/ground loops, environmental sources) Filtering techniques Hardware: Passive (RC); Active (op amp, powered) Software: Sampling issues Amplification Measurement Systems: The response of a measurement device to a dynamically-changing physical parameter (pressure, temperature, etc) that it is seeking to measure. All measurement systems must interact with the physical parameter in order to measure it, so the inherent response of the measurement device will be involved with the measurement signal output. The output signal will be affected by how the measurement is performed (ie, the physical characteristics of the measurement instrument). The input parameter may be thought of as a driving, or stimulus energy, since it is causing the instrument to respond to it. We are now incorporating the effect of time in the measurement response dynamic signals. We have studied types of measurement systems (1st and nd-order) and their response to different types of driving parameters: Step input: where H(t) is the Heaviside function. Sinusoidal input: where x(t) = AH(t) (47) x(t) = A sin(ωt + Φ 0 ) + A 0 (48) T = 1 ; ω = πf (49) f 7
8 First Order Systems 1st-Order systems only have a storage component. One example is a thermocouple or thermometer. They have no inertia, so there is no oscillatory reaction to a driving system. The storage aspect means that the system will have some inherent response time to a change in input parameter. The system (regardless of the input parameter x(t)) may be modeled as: τẏ(t) + y(t) = Kx(t) (50) where K is the static sensitivity and τ is the time constant. Response to input parameter x(t): Step input: x(t) = AH(t) The solution to the model is: y(t) = KA + (y 0 KA)e t/τ + O (51) Sinusoidal input: x(t) = A sin(ωt + Φ 0 ) + A 0 The steady-state solution to the model is: y(t) = KAM(ω) sin(ωt + Φ(ω) + Φ 0 ) + KA 0 + O (5) where M(ω) is the Magnitude Ratio, Φ(ω) is the phase lag, and O is the instrument-specific DC offset. The frequency-dependent terms may be modeled as: M(ω) = (ωτ) (53) (54) Φ(ω) = tan 1 (ωτ) (55) Quantifying the response Dynamic error Γ(t): Frequency-based dynamic error δ(ω): Γ(t) = y(t) y y 0 y = e t/τ (56) δ(ω) = M(ω) 1 (57) Rise time, t rise : Time that it takes to reach 90% of the steady state value: y(t rise ) = 0.9(y y 0 ). This occurs when t =.3τ. For 1st-order, step input functions, t rise = t response. The system reaches 99.3% of the steady-state value when t = 5τ At t = τ, the system has reached 63% of the steady-state value 8
9 Second Order Systems nd-order systems inherently have a mass, spring, and damping, which are all involved in the response of the system to the driving signal. The system (regardless of the input parameter x(t)) may be modeled as: mÿ(t) + cẏ(t) + ky(t) = Kx(t) (58) where m is the mass, c is the damping coefficient, k is the spring constant, and K is the static sensitivity. Since the system has a second-order dependence on time, its response will depend on highly on the values of these three parameters (m, c, k). We can simplify the equation slightly, by dividing through by m, and restating as: ÿ(t) + c mẏ(t) + ω ny(t) = Kx(t) (59) where ω n = k/m. The response is based on the value of the generic solution to the homogeneous equation (x(t) = 0), Be λt : λ 1, = c m ± c 4km m (60) The solution is imaginary for a critical value of c: c cr = km. This defines the damping ratio: ζ = c c cr = c. km Simplifying the solution: λ 1, = ζω n ± ω n ζ 1 (61) It s now easy to see the boundary for the importance of ζ, since it dictates the solution to the model. 3 solutions arise: ζ < 1: Underdamped solution y(t) = Ce ζωnt sin(ω d t + Φ) (6) where ω d = ω n 1 ζ : the damped natural frequency. ζ = 1: Critically damped solution y(t) = C 1 e λ 1t + C te λ t (63) ζ > 1: Overdamped solution y(t) = C 1 e λ 1t + C e λ t (64) Not much happens in the critically and overdamped cases: most of our attention has been focused on the underdamped scenarios. 9
10 Response to input parameter x(t): We ve mostly focused on a sinusoidal input, for underdamped systems. Sinusoidal input: x(t) = A sin(ωt + Φ 0 ) + A 0. The steady-state solution is the same as 1st-order systems: y(t) = KAM(ω) sin(ωt + Φ(ω) + Φ 0 ) + KA 0 + O (65) The difference between systems is the magnitude ratio and phase lag form: M(ω) = 1 ( ( ) ) ω 1 ω n + ( ζ ω ) Φ(ω) = tan 1 ω n 1 (ω/ω n ) ( ζ ω ω n ) (66) Resonance: When output amplitude exceeds input amplitude. Resonance occurs at the maximum amplitude, M(ω = ω r ) = M max, where ω r is the resonance frequency: (67) (68) ω r = ω n 1 ζ (69) Q = ω nτ relax = ω r ω ω 1 = The Quality Factor Q defines the sharpness of resonance. (70) 1 ζ 1 ζ = M max (71) Step input: x(t) = AH(t) + A 0 The solution to the model for an underdamped system, where A 0 is the previous amplitude of the input and A is the new change in input amplitude: y(t) = KA + O + (KA 0 KA) exp( ζω n t) [ ζ 1 ζ sin(ω dt) + cos(ω d t) Rise time, t rise : Time that it takes to first reach 90% of the steady state value Response time, t response, or t settle : Time that it takes to stay within 10% of the steady-state value Think of as the free decay: No driving signal is applied. The system s natural response takes over, and its damping governs the decay. The system responds at its damped natural frequency ω d : ] (7) ω d = ω n 1 ζ (73) Special case: IF the scenario involves an initial input amplitude A 0 and a final amplitude of A = 0, its envelope has an exponential decay: y(t) = KA 0 exp( ζω n t) (74) 10
11 τ relax is the time that it takes the undriven signal to decay from its previous excitation state to y(t = τ relax ) = 0.37y 0. Quantifying the response Time-based dynamic error Γ(t): Miscellaneous Frequency-based dynamic error δ(ω): Optimal value: ζ = Γ(t) = y(t) y y 0 y (75) δ(ω) = M(ω) 1 (76) First order systems asymptote to Φ = π/ for large frequencies Second order systems asymptote to Φ = π for large frequencies First order systems pass through Φ = π/4 when ω = 1/τ Second order systems pass through Φ = π/ when ω = ω n The measurement always lags (ie, is behind) the driving signal, so Φ(ω) 0. Because the system has a physical inertial response, it cannot mechanically respond instantaneously! Quality Factor, Q A unitless parameter that describes the sharpness of the instrument s resonance response to an oscillatory input: where and Q = M max = 1 ζ 1 ζ = ω r ω BW (77) M max = M(ω r ), (78) ω BW = ω ω 1, (79) and ω 1, ω correspond to the frequencies at which M(ω) = 0.707M max, or when the response is at its half-power relative to its maximum amplitude (at the resonant frequency). The Quality Factor may also be visualized during a step input response, or free decay scenario. In this scenario, the instrument s natural response takes over, and its damping governs the decay. The system responds at its damped natural frequency ω d : ω d = ω n 1 ζ (80) 11
12 and its envelope has an exponential decay: y(t) = A 0 exp( ζω n t) (81) The time constant τ relax is the time that it takes the instrument output to decay from its previous excitation state to y(t = τ relax ) = 0.37A 0. The Quality Factor is related to this time constant: Q = ω nτ relax Frequency analysis Note: See the course text (Chapter: The Analog Measurand: Time-Dependent Characteristics) for the full derivation of the equations. Premise: If a measured signal y(t) contains energy at multiple frequencies, it will likely be difficult to visually determine the amplitude of each frequency. To better understand the processes that are going on in the signal, we can use Fourier Analysis. Fourier analysis allows us to decompose the measured time-domain signal into all of its frequency components. This is ultimately accomplished by multiplying the time-domain signal by its sin and cos components and integrating the signal over some time window, in an operation called the Fourier Transform: Y (ω) = F(y(t)) = (8) y(t) exp iωt dt (83) Y (ω) = A(ω) ib(ω) = Y (ω) exp iφ(ω) (84) 1 Im(Y (ω)) Φ(ω) = tan Re(Y (ω)) (85) where Y (ω) is the complex signal at each frequency ω in the signal spectrum, A represents the real (amplitude) component and B represents the imaginary (phase) component. Note that these equations can all be expressed in terms of f instead of ω. Fourier analysis and instrumentation: To apply the Fourier Transform to digitized data, you will want to use a Fast Fourier Transform (FFT) in Matlab. There are a few important things to remember for this process: The FFT returns a complex vector that has the same number of points as the data y(t). The absolute value represents the amplitude and the imaginary values represent the phase angle. In the context of ME310, we only care about the amplitude. The maximum frequency represented by Y (f) is equal to the Nyquist frequency, or f s /, and the frequency resolution f res is equal to f s /N, where N is the number of data points. The FFT returns the negative and positive frequency components, of which we only need to use the positive components. Practically speaking, this means we only need half of the points contained in Y (f). In Matlab, the first index value in Y (f) is the DC offset, or A 0. The second index point corresponds to the lowest frequency f min, the third point corresponds to f min + f res, and so on. This is important to remember when plotting the frequency amplitude or accessing individual index values. 1
13 The amplitude Y (f) does not have the same units as y(t)! The amplitude scaling in Matlab involves a complicated relationship between the y(t) units, frequency, and the number of data points, so it is best to simply think of the amplitude as a relative scale. The time-domain data can be recovered from Y (f) by applying the inverse FFT, or IFFT ( ifft in Matlab). This involves both the complex value Y (f), not just the amplitude. 13
14 Boston University ME 310: nd exam review questions Prof. C. Farny Review problems; solutions on next page: 1. A pressure transducer behaves as a second-order system. If the natural frequency is 4 khz and the damping ratio is 0.75, determine the resonance frequency, as well as the frequency range(s) over which the measurement error is not greater than 5%. Repeat this exercise with a damping ratio of A mercury glass thermometer that initially reads 5 is suddenly immersed into a liquid that is maintained at 100. After a time interval of.0 s, the thermometer reads 76. Assuming the thermometer is a first-order system, estimate its time constant. 3. The output from a temperature sensor is expected to vary from.5 mv to 3.5 mv. Its output is sampled by a 1-bit AD converter with a ±5 V range. a) Estimate the resolution of the AD system. How does this compare to the temperature magnitude from the sensor? b) You implement a preamplifier to improve the measurement accuracy/reduce uncertainty. What standard gain setting (state in both linear and in db scale) is appropriate for such a preamplifier? Most preamplifiers offer steps of either 6 or 0 db. 14
15 Solutions: 1. For ζ = 0.75, the system has no resonance, so there is no concept of a resonance frequency. Solve for the frequency range by plugging into Eqn 76 and solving for M. M will be flat up until some frequency, and this is the frequency you re interested in finding. Plug into Eqn 66 and solve for f. For ζ = 0.75, f = 1900 Hz. For ζ = 0.5, the system will be resonant, so there will be a resonant frequency (see Eqn 69), where ω r = krad/s, or f r =.88 khz. There will exist two frequency regions that satisfy the criteria for the error: 0 < f < 1300 Hz; 3785 < f < 4195 Hz.. Solve by plugging into Eqn 51 and solving for τ: 76 = (5 100) exp ( /τ) (86) τ = 1.76 s (87) 3. a) Use Eqn 8, which returns a resolution of.44 mv on the same magnitude as the temperature output you re looking to measure. b) Comparing the V F S of the AD converter with the maximum output from the sensor, the ratio is 149. Choosing a standard amplification value, you should choose a gain of 1000, or 60 db. 15
Coverage: through class #7 material (updated June 14, 2017) Terminology/miscellaneous. for a linear instrument:
Boston University ME 310: Summer 017 Course otes to-date Prof. C. Farny DISCLAIMER: This document is meant to serve as a reference and an overview, and not an exclusive study guide for the course. Coverage:
More informationChapter 3. Second Order Measurement System Behaviour
Chapter 3 Second Order Measurement System Behaviour 1 GENERAL MEASUREMENT SYSTEM MODEL INITIAL CONDITIONS y(0) SIGNAL INPUT F(t) MEASUREMENT SYSTEM SYSTEM OUTPUT y(t) Understanding the theoretical response
More informationDepartment of Mechanical and Aerospace Engineering. MAE334 - Introduction to Instrumentation and Computers. Final Examination.
Name: Number: Department of Mechanical and Aerospace Engineering MAE334 - Introduction to Instrumentation and Computers Final Examination December 12, 2003 Closed Book and Notes 1. Be sure to fill in your
More information(A) (B) (D) (C) 1.5. Amplitude (volts) 1.5. Amplitude (volts) Time (seconds) Time (seconds)
Reminder: Lab #1 : Limitations of A/D conversion Lab #2 : Thermocouple, static and dynamic calibration Lab #3 : Conversion of work into heat Lab #4 : Pressure transducer, static and dynamic calibration
More informationDynamic circuits: Frequency domain analysis
Electronic Circuits 1 Dynamic circuits: Contents Free oscillation and natural frequency Transfer functions Frequency response Bode plots 1 System behaviour: overview 2 System behaviour : review solution
More informationChapter 3. 1 st Order Sine Function Input. General Solution. Ce t. Measurement System Behavior Part 2
Chapter 3 Measurement System Behavior Part 2 1 st Order Sine Function Input Examples of Periodic: vibrating structure, vehicle suspension, reciprocating pumps, environmental conditions The frequency of
More informationVibrations: Second Order Systems with One Degree of Freedom, Free Response
Single Degree of Freedom System 1.003J/1.053J Dynamics and Control I, Spring 007 Professor Thomas Peacock 5//007 Lecture 0 Vibrations: Second Order Systems with One Degree of Freedom, Free Response Single
More informationDepartment of Mechanical and Aerospace Engineering. MAE334 - Introduction to Instrumentation and Computers. Final Examination.
Name: Number: Department of Mechanical and Aerospace Engineering MAE334 - Introduction to Instrumentation and Computers Final Examination December 17, 2001 Closed Book and Notes 1. Be sure to fill in your
More informationControl Engineering BDA30703
Control Engineering BDA30703 Lecture 3: Performance characteristics of an instrument Prepared by: Ramhuzaini bin Abd. Rahman Expected Outcomes At the end of this lecture, students should be able to; 1)
More informationDamped Oscillation Solution
Lecture 19 (Chapter 7): Energy Damping, s 1 OverDamped Oscillation Solution Damped Oscillation Solution The last case has β 2 ω 2 0 > 0. In this case we define another real frequency ω 2 = β 2 ω 2 0. In
More informationLaboratory III: Operational Amplifiers
Physics 33, Fall 2008 Lab III - Handout Laboratory III: Operational Amplifiers Introduction Operational amplifiers are one of the most useful building blocks of analog electronics. Ideally, an op amp would
More informationLab 11 - Free, Damped, and Forced Oscillations
Lab 11 Free, Damped, and Forced Oscillations L11-1 Name Date Partners Lab 11 - Free, Damped, and Forced Oscillations OBJECTIVES To understand the free oscillations of a mass and spring. To understand how
More informationNotes on the Periodically Forced Harmonic Oscillator
Notes on the Periodically orced Harmonic Oscillator Warren Weckesser Math 38 - Differential Equations 1 The Periodically orced Harmonic Oscillator. By periodically forced harmonic oscillator, we mean the
More informationProblem 1: Lagrangians and Conserved Quantities. Consider the following action for a particle of mass m moving in one dimension
105A Practice Final Solutions March 13, 01 William Kelly Problem 1: Lagrangians and Conserved Quantities Consider the following action for a particle of mass m moving in one dimension S = dtl = mc dt 1
More informationThe Harmonic Oscillator
The Harmonic Oscillator Math 4: Ordinary Differential Equations Chris Meyer May 3, 008 Introduction The harmonic oscillator is a common model used in physics because of the wide range of problems it can
More informationIntroduction to Vibration. Mike Brennan UNESP, Ilha Solteira São Paulo Brazil
Introduction to Vibration Mike Brennan UNESP, Ilha Solteira São Paulo Brazil Vibration Most vibrations are undesirable, but there are many instances where vibrations are useful Ultrasonic (very high
More informationSource-Free RC Circuit
First Order Circuits Source-Free RC Circuit Initial charge on capacitor q = Cv(0) so that voltage at time 0 is v(0). What is v(t)? Prof Carruthers (ECE @ BU) EK307 Notes Summer 2018 150 / 264 First Order
More informationBiomedical Instrumentation System
BME/EECS 458 - Biomedical Instrumentation & Design Matt O Donnell I.0 Introduction What is a biomedical instrument? To many it s an EKG machine, to others it s a chemical biosensor, and to some it s a
More informationNUMERICAL EVALUATION OF DISPLACEMENT AND ACCELERATION FOR A MASS, SPRING, DASHPOT SYSTEM
American Society for Engineering Education Salt Lake City, Utah June 24 NUMERICAL EVALUATION OF DISPLACEMENT AND ACCELERATION FOR A MASS, SPRING, DASHPOT SYSTEM ACCELEROMETER MASS SPRING FIXED SUPPORT
More informationDeveloping a Multisemester Interwoven Dynamic Systems Project to Foster Learning and Retention of STEM Material
100 10 1 ζ=0.1% ζ=1% ζ=2% ζ=5% ζ=10% ζ=20% 0-90 -180 ζ=20% ζ=10% ζ=5% ζ=2% ζ=1% ζ=0.1% Proceedings of 2004 IMECE: November 13-19, 2004, Anaheim, CA 2004 ASME International Mechanical Engineering Congress
More informationDamped harmonic motion
Damped harmonic motion March 3, 016 Harmonic motion is studied in the presence of a damping force proportional to the velocity. The complex method is introduced, and the different cases of under-damping,
More informationClassical Mechanics Phys105A, Winter 2007
Classical Mechanics Phys5A, Winter 7 Wim van Dam Room 59, Harold Frank Hall vandam@cs.ucsb.edu http://www.cs.ucsb.edu/~vandam/ Phys5A, Winter 7, Wim van Dam, UCSB Midterm New homework has been announced
More informationOscillations. Simple Harmonic Motion of a Mass on a Spring The equation of motion for a mass m is attached to a spring of constant k is
Dr. Alain Brizard College Physics I (PY 10) Oscillations Textbook Reference: Chapter 14 sections 1-8. Simple Harmonic Motion of a Mass on a Spring The equation of motion for a mass m is attached to a spring
More informationDepartment of Mechanical and Aerospace Engineering MAE334 - Introduction to Instrumentation and Computers. Midterm Examination.
Department of Mechanical and Aerospace Engineering MAE334 - Introduction to Instrumentation and Computers Midterm Examination October 19, 2005 o Closed Book and Notes o Fill in your name on your scoring
More information4. Complex Oscillations
4. Complex Oscillations The most common use of complex numbers in physics is for analyzing oscillations and waves. We will illustrate this with a simple but crucially important model, the damped harmonic
More informationCoupled Electrical Oscillators Physics Advanced Physics Lab - Summer 2018 Don Heiman, Northeastern University, 5/24/2018
Coupled Electrical Oscillators Physics 3600 - Advanced Physics Lab - Summer 08 Don Heiman, Northeastern University, 5/4/08 I. INTRODUCTION The objectives of this experiment are: () explore the properties
More informationMODULE I. Transient Response:
Transient Response: MODULE I The Transient Response (also known as the Natural Response) is the way the circuit responds to energies stored in storage elements, such as capacitors and inductors. If a capacitor
More informationModeling and Experimentation: Mass-Spring-Damper System Dynamics
Modeling and Experimentation: Mass-Spring-Damper System Dynamics Prof. R.G. Longoria Department of Mechanical Engineering The University of Texas at Austin July 20, 2014 Overview 1 This lab is meant to
More informationTHE FOURIER TRANSFORM (Fourier series for a function whose period is very, very long) Reading: Main 11.3
THE FOURIER TRANSFORM (Fourier series for a function whose period is very, very long) Reading: Main 11.3 Any periodic function f(t) can be written as a Fourier Series a 0 2 + a n cos( nωt) + b n sin n
More informationDynamic Response. Assoc. Prof. Enver Tatlicioglu. Department of Electrical & Electronics Engineering Izmir Institute of Technology.
Dynamic Response Assoc. Prof. Enver Tatlicioglu Department of Electrical & Electronics Engineering Izmir Institute of Technology Chapter 3 Assoc. Prof. Enver Tatlicioglu (EEE@IYTE) EE362 Feedback Control
More informationME EXAM #1 Tuesday, October 7, :30 7:30 pm PHYS 112 and 114
ME 36500 EXAM # Tuesday, October 7, 04 6:30 7:30 pm PHYS and 4 Division: Chiu(0:30) / Shelton(:30) / Bae(:30) (circle one) HW ID: Name: Solution Instructions () This is a closed book examination, but you
More informationProcess Control & Instrumentation (CH 3040)
First-order systems Process Control & Instrumentation (CH 3040) Arun K. Tangirala Department of Chemical Engineering, IIT Madras January - April 010 Lectures: Mon, Tue, Wed, Fri Extra class: Thu A first-order
More informationAspects of Continuous- and Discrete-Time Signals and Systems
Aspects of Continuous- and Discrete-Time Signals and Systems C.S. Ramalingam Department of Electrical Engineering IIT Madras C.S. Ramalingam (EE Dept., IIT Madras) Networks and Systems 1 / 45 Scaling the
More informationDeveloping a Virtual Model of a Second Order System to Simulation Real Laboratory Measurement Problems
Proceedings of 24 IMECE: November 13-19, 24, Anaheim, CA 24 ASME International Mechanical Engineering Congress and RD&D Expo Developing a Virtual Model of a Second Order System to Simulation Real Laboratory
More informationEE 230 Lecture 40. Data Converters. Amplitude Quantization. Quantization Noise
EE 230 Lecture 40 Data Converters Amplitude Quantization Quantization Noise Review from Last Time: Time Quantization Typical ADC Environment Review from Last Time: Time Quantization Analog Signal Reconstruction
More information1 Pushing your Friend on a Swing
Massachusetts Institute of Technology MITES 017 Physics III Lecture 05: Driven Oscillations In these notes, we derive the properties of both an undamped and damped harmonic oscillator under the influence
More informationIntroduction to Vibration. Professor Mike Brennan
Introduction to Vibration Professor Mie Brennan Introduction to Vibration Nature of vibration of mechanical systems Free and forced vibrations Frequency response functions Fundamentals For free vibration
More informationPatrick F. Dunn 107 Hessert Laboratory Department of Aerospace and Mechanical Engineering University of Notre Dame Notre Dame, IN 46556
Learning Objectives to Accompany MEASUREMENT AND DATA ANALYSIS FOR ENGINEERING AND SCIENCE Second Edition, Taylor and Francis/CRC Press, c 2010 ISBN: 9781439825686 Patrick F. Dunn pdunn@nd.edu 107 Hessert
More informationThe Torsion Pendulum (One or two weights)
The Torsion Pendulum (One or two weights) Exercises I through V form the one-weight experiment. Exercises VI and VII, completed after Exercises I -V, add one weight more. Preparatory Questions: 1. The
More informationR-L-C Circuits and Resonant Circuits
P517/617 Lec4, P1 R-L-C Circuits and Resonant Circuits Consider the following RLC series circuit What's R? Simplest way to solve for is to use voltage divider equation in complex notation. X L X C in 0
More informationProblem Weight Total 100
EE 350 Problem Set 3 Cover Sheet Fall 2016 Last Name (Print): First Name (Print): ID number (Last 4 digits): Section: Submission deadlines: Turn in the written solutions by 4:00 pm on Tuesday September
More informationChaos and R-L diode Circuit
Chaos and R-L diode Circuit Rabia Aslam Chaudary Roll no: 2012-10-0011 LUMS School of Science and Engineering Thursday, December 20, 2010 1 Abstract In this experiment, we will use an R-L diode circuit
More informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Illinois Institute of Technology Lecture 8: Response Characteristics Overview In this Lecture, you will learn: Characteristics of the Response Stability Real
More informationOutline. Introduction
Outline Periodic and Non-periodic Dipartimento di Ingegneria Civile e Ambientale, Politecnico di Milano March 5, 4 A periodic loading is characterized by the identity p(t) = p(t + T ) where T is the period
More informationTransient Response of a Second-Order System
Transient Response of a Second-Order System ECEN 830 Spring 01 1. Introduction In connection with this experiment, you are selecting the gains in your feedback loop to obtain a well-behaved closed-loop
More informationTrial version. Temperature Sensing. How does the temperature sensor work and how can it be used to control the temperature of a refrigerator?
Temperature Sensing How does the temperature sensor work and how can it be used to control the temperature of a refrigerator? Temperature Sensing page: 1 of 13 Contents Initial Problem Statement 2 Narrative
More informationFigure 1.1.1: Magnitude ratio as a function of time using normalized temperature data
9Student : Solutions to Technical Questions Assignment : Lab 3 Technical Questions Date : 1 November 212 Course Title : Introduction to Measurements and Data Analysis Course Number : AME 2213 Instructor
More informationOscillatory Motion SHM
Chapter 15 Oscillatory Motion SHM Dr. Armen Kocharian Periodic Motion Periodic motion is motion of an object that regularly repeats The object returns to a given position after a fixed time interval A
More informationEE 40: Introduction to Microelectronic Circuits Spring 2008: Midterm 2
EE 4: Introduction to Microelectronic Circuits Spring 8: Midterm Venkat Anantharam 3/9/8 Total Time Allotted : min Total Points:. This is a closed book exam. However, you are allowed to bring two pages
More informationLecture 4: R-L-C Circuits and Resonant Circuits
Lecture 4: R-L-C Circuits and Resonant Circuits RLC series circuit: What's V R? Simplest way to solve for V is to use voltage divider equation in complex notation: V X L X C V R = in R R + X C + X L L
More informationLaboratory I: Impedance
Physics 331, Fall 2008 Lab I - Handout 1 Laboratory I: Impedance Reading: Simpson Chapter 1 (if necessary) & Chapter 2 (particularly 2.9-2.13) 1 Introduction In this first lab we review the properties
More informationElectronic Circuits. Prof. Dr. Qiuting Huang Integrated Systems Laboratory
Electronic Circuits Prof. Dr. Qiuting Huang 6. Transimpedance Amplifiers, Voltage Regulators, Logarithmic Amplifiers, Anti-Logarithmic Amplifiers Transimpedance Amplifiers Sensing an input current ii in
More informationLinear second-order differential equations with constant coefficients and nonzero right-hand side
Linear second-order differential equations with constant coefficients and nonzero right-hand side We return to the damped, driven simple harmonic oscillator d 2 y dy + 2b dt2 dt + ω2 0y = F sin ωt We note
More informationResponse to Periodic and Non-periodic Loadings. Giacomo Boffi. March 25, 2014
Periodic and Non-periodic Dipartimento di Ingegneria Civile e Ambientale, Politecnico di Milano March 25, 2014 Outline Introduction Fourier Series Representation Fourier Series of the Response Introduction
More informationLecture V : Oscillatory motion and spectral analysis
Lecture V : Oscillatory motion and spectral analysis I. IDEAL PENDULUM AND STABILITY ANALYSIS Let us remind ourselves of the equation of motion for the pendulum. Remembering that the external torque applied
More informationSolved Problems. Electric Circuits & Components. 1-1 Write the KVL equation for the circuit shown.
Solved Problems Electric Circuits & Components 1-1 Write the KVL equation for the circuit shown. 1-2 Write the KCL equation for the principal node shown. 1-2A In the DC circuit given in Fig. 1, find (i)
More informationSinusoidal Steady-State Analysis
Chapter 4 Sinusoidal Steady-State Analysis In this unit, we consider circuits in which the sources are sinusoidal in nature. The review section of this unit covers most of section 9.1 9.9 of the text.
More information1. (10 points) Find the general solution to the following second-order differential equation:
Math 307A, Winter 014 Midterm Solutions Page 1 of 8 1. (10 points) Find the general solution to the following second-order differential equation: 4y 1y + 9y = 9t. To find the general solution to this nonhomogeneous
More informationENGR 2405 Chapter 8. Second Order Circuits
ENGR 2405 Chapter 8 Second Order Circuits Overview The previous chapter introduced the concept of first order circuits. This chapter will expand on that with second order circuits: those that need a second
More informationDynamics of structures
Dynamics of structures 1.2 Viscous damping Luc St-Pierre October 30, 2017 1 / 22 Summary so far We analysed the spring-mass system and found that its motion is governed by: mẍ(t) + kx(t) = 0 k y m x x
More informationChapter 3: Capacitors, Inductors, and Complex Impedance
hapter 3: apacitors, Inductors, and omplex Impedance In this chapter we introduce the concept of complex resistance, or impedance, by studying two reactive circuit elements, the capacitor and the inductor.
More informationCircuits Advanced Topics by Dr. Colton (Fall 2016)
ircuits Advanced Topics by Dr. olton (Fall 06). Time dependence of general and L problems General and L problems can always be cast into first order ODEs. You can solve these via the particular solution
More informationTime Varying Circuit Analysis
MAS.836 Sensor Systems for Interactive Environments th Distributed: Tuesday February 16, 2010 Due: Tuesday February 23, 2010 Problem Set # 2 Time Varying Circuit Analysis The purpose of this problem set
More informationCMPT 889: Lecture 3 Fundamentals of Digital Audio, Discrete-Time Signals
CMPT 889: Lecture 3 Fundamentals of Digital Audio, Discrete-Time Signals Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University October 6, 2005 1 Sound Sound waves are longitudinal
More informationDifferential Equations
Differential Equations A differential equation (DE) is an equation which involves an unknown function f (x) as well as some of its derivatives. To solve a differential equation means to find the unknown
More informationVibrations of Single Degree of Freedom Systems
Vibrations of Single Degree of Freedom Systems CEE 541. Structural Dynamics Department of Civil and Environmental Engineering Duke University Henri P. Gavin Fall, 16 This document describes free and forced
More information2.3 Damping, phases and all that
2.3. DAMPING, PHASES AND ALL THAT 107 2.3 Damping, phases and all that If we imagine taking our idealized mass on a spring and dunking it in water or, more dramatically, in molasses), then there will be
More informationSensor Characteristics
Lecture (3) Sensor Characteristics (Part Two) Prof. Kasim M. Al-Aubidy Philadelphia University-Jordan AMSS-MSc Prof. Kasim Al-Aubidy 1 3. Computation of Stimulus: The main objective of sensing is to determine
More informationBasic Procedures for Common Problems
Basic Procedures for Common Problems ECHE 550, Fall 2002 Steady State Multivariable Modeling and Control 1 Determine what variables are available to manipulate (inputs, u) and what variables are available
More informationPhysics 8 Monday, December 4, 2017
Physics 8 Monday, December 4, 2017 HW12 due Friday. Grace will do a review session Dec 12 or 13. When? I will do a review session: afternoon Dec 17? Evening Dec 18? Wednesday, I will hand out the practice
More informationagree w/input bond => + sign disagree w/input bond => - sign
1 ME 344 REVIEW FOR FINAL EXAM LOCATION: CPE 2.204 M. D. BRYANT DATE: Wednesday, May 7, 2008 9-noon Finals week office hours: May 6, 4-7 pm Permitted at final exam: 1 sheet of formulas & calculator I.
More informationResponse to Periodic and Non-periodic Loadings. Giacomo Boffi.
Periodic and Non-periodic http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano March 7, 218 Outline Undamped SDOF systems
More informationForced Oscillation and Resonance
Chapter Forced Oscillation and Resonance The forced oscillation problem will be crucial to our understanding of wave phenomena Complex exponentials are even more useful for the discussion of damping and
More informationSTABILITY. Have looked at modeling dynamic systems using differential equations. and used the Laplace transform to help find step and impulse
SIGNALS AND SYSTEMS: PAPER 3C1 HANDOUT 4. Dr David Corrigan 1. Electronic and Electrical Engineering Dept. corrigad@tcd.ie www.sigmedia.tv STABILITY Have looked at modeling dynamic systems using differential
More informationik () uk () Today s menu Last lecture Some definitions Repeatability of sensing elements
Last lecture Overview of the elements of measurement systems. Sensing elements. Signal conditioning elements. Signal processing elements. Data presentation elements. Static characteristics of measurement
More informationMSE 383, Unit 3-3. Joshua U. Otaigbe Iowa State University Materials Science & Engineering Dept.
Dynamic Mechanical Behavior MSE 383, Unit 3-3 Joshua U. Otaigbe Iowa State University Materials Science & Engineering Dept. Scope Why DMA & TTS? DMA Dynamic Mechanical Behavior (DMA) Superposition Principles
More informationIndex. Index. More information. in this web service Cambridge University Press
A-type elements, 4 7, 18, 31, 168, 198, 202, 219, 220, 222, 225 A-type variables. See Across variable ac current, 172, 251 ac induction motor, 251 Acceleration rotational, 30 translational, 16 Accumulator,
More informationTable of contents. d 2 y dx 2, As the equation is linear, these quantities can only be involved in the following manner:
M ath 0 1 E S 1 W inter 0 1 0 Last Updated: January, 01 0 Solving Second Order Linear ODEs Disclaimer: This lecture note tries to provide an alternative approach to the material in Sections 4. 4. 7 and
More informationChapter 15. Oscillatory Motion
Chapter 15 Oscillatory Motion Part 2 Oscillations and Mechanical Waves Periodic motion is the repeating motion of an object in which it continues to return to a given position after a fixed time interval.
More informationRaktim Bhattacharya. . AERO 422: Active Controls for Aerospace Vehicles. Frequency Response-Design Method
.. AERO 422: Active Controls for Aerospace Vehicles Frequency Response- Method Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. ... Response to
More informationTopic 5 Notes Jeremy Orloff. 5 Homogeneous, linear, constant coefficient differential equations
Topic 5 Notes Jeremy Orloff 5 Homogeneous, linear, constant coefficient differential equations 5.1 Goals 1. Be able to solve homogeneous constant coefficient linear differential equations using the method
More information5.74 Introductory Quantum Mechanics II
MIT OpenCourseWare http://ocw.mit.edu 5.74 Introductory Quantum Mechanics II Spring 9 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Andrei Tokmakoff,
More informationLab 1: Damped, Driven Harmonic Oscillator
1 Introduction Lab 1: Damped, Driven Harmonic Oscillator The purpose of this experiment is to study the resonant properties of a driven, damped harmonic oscillator. This type of motion is characteristic
More informationOutline of parts 1 and 2
to Harmonic Loading http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano March, 6 Outline of parts and of an Oscillator
More informationMathematical Physics
Mathematical Physics MP205 Vibrations and Waves Lecturer: Office: Lecture 9-10 Dr. Jiří Vala Room 1.9, Mathema
More informationIn-class exercises Day 1
Physics 4488/6562: Statistical Mechanics http://www.physics.cornell.edu/sethna/teaching/562/ Material for Week 11 Exercises due Mon Apr 16 Last correction at April 16, 2018, 11:19 am c 2018, James Sethna,
More informationChapter 2: Linear Constant Coefficient Higher Order Equations
Chapter 2: Linear Constant Coefficient Higher Order Equations The wave equation is a linear partial differential equation, which means that sums of solutions are still solutions, just as for linear ordinary
More informationEE 242 EXPERIMENT 8: CHARACTERISTIC OF PARALLEL RLC CIRCUIT BY USING PULSE EXCITATION 1
EE 242 EXPERIMENT 8: CHARACTERISTIC OF PARALLEL RLC CIRCUIT BY USING PULSE EXCITATION 1 PURPOSE: To experimentally study the behavior of a parallel RLC circuit by using pulse excitation and to verify that
More informationBIOEN 302, Section 3: AC electronics
BIOEN 3, Section 3: AC electronics For this section, you will need to have very present the basics of complex number calculus (see Section for a brief overview) and EE5 s section on phasors.. Representation
More informationChapter 4 Transients. Chapter 4 Transients
Chapter 4 Transients Chapter 4 Transients 1. Solve first-order RC or RL circuits. 2. Understand the concepts of transient response and steady-state response. 1 3. Relate the transient response of first-order
More informationBasic Electronics. Introductory Lecture Course for. Technology and Instrumentation in Particle Physics Chicago, Illinois June 9-14, 2011
Basic Electronics Introductory Lecture Course for Technology and Instrumentation in Particle Physics 2011 Chicago, Illinois June 9-14, 2011 Presented By Gary Drake Argonne National Laboratory Session 2
More informationPHY217: Vibrations and Waves
Assessed Problem set 1 Issued: 5 November 01 PHY17: Vibrations and Waves Deadline for submission: 5 pm Thursday 15th November, to the V&W pigeon hole in the Physics reception on the 1st floor of the GO
More informationSection 3.4. Second Order Nonhomogeneous. The corresponding homogeneous equation
Section 3.4. Second Order Nonhomogeneous Equations y + p(x)y + q(x)y = f(x) (N) The corresponding homogeneous equation y + p(x)y + q(x)y = 0 (H) is called the reduced equation of (N). 1 General Results
More informationTheory and Design for Mechanical Measurements
Theory and Design for Mechanical Measurements Third Edition Richard S. Figliola Clemson University Donald E. Beasley Clemson University John Wiley & Sons, Inc. New York / Chichester / Weinheim / Brisbane
More informationLaboratory handouts, ME 340
Laboratory handouts, ME 340 This document contains summary theory, solved exercises, prelab assignments, lab instructions, and report assignments for Lab 4. 2014-2016 Harry Dankowicz, unless otherwise
More informationLaboratory notes. Torsional Vibration Absorber
Titurus, Marsico & Wagg Torsional Vibration Absorber UoB/1-11, v1. Laboratory notes Torsional Vibration Absorber Contents 1 Objectives... Apparatus... 3 Theory... 3 3.1 Background information... 3 3. Undamped
More informationC(s) R(s) 1 C(s) C(s) C(s) = s - T. Ts + 1 = 1 s - 1. s + (1 T) Taking the inverse Laplace transform of Equation (5 2), we obtain
analyses of the step response, ramp response, and impulse response of the second-order systems are presented. Section 5 4 discusses the transient-response analysis of higherorder systems. Section 5 5 gives
More information4.2 Homogeneous Linear Equations
4.2 Homogeneous Linear Equations Homogeneous Linear Equations with Constant Coefficients Consider the first-order linear differential equation with constant coefficients a 0 and b. If f(t) = 0 then this
More informationPhysics 106a/196a Problem Set 7 Due Dec 2, 2005
Physics 06a/96a Problem Set 7 Due Dec, 005 Version 3, Nov 7, 005 In this set we finish up the SHO and study coupled oscillations/normal modes and waves. Problems,, and 3 are for 06a students only, 4, 5,
More informationWaves 2006 Physics 23. Armen Kocharian Lecture 3: Sep
Waves 2006 Physics 23 Armen Kocharian Lecture 3: Sep 12. 2006 Last Time What is a wave? A "disturbance" that moves through space. Mechanical waves through a medium. Transverse vs. Longitudinal e.g., string
More information