Modeling and Simulation Revision III D R. T A R E K A. T U T U N J I P H I L A D E L P H I A U N I V E R S I T Y, J O R D A N 0 1 4
Block Diagrams Block diagram models consist of two fundamental objects: signal blocks and wires. A block is a processing element which operates on input signals and parameters to produce output signals A wire is to transmits a signal from its origination point (usually a block) to its termination point (usually another block). Block diagrams are suitable to represent multidisciplinary models that represent a physical phenomenon. Dr. Tarek A. Tutunji
Block Diagram Example [Ref.] Prof. Shetty
Block Diagrams Manipulation
Block Diagrams Manipulation
Block Diagrams: Direct Method Example Consider the transfer function: We can introduce s state variable, x(t), in order to separate the polynomials
State Equation The differential equation is: Put the needed integrator blocks: Add the required multipliers to obtain the state equation:
Output Equation Repeat the same procedure for the output equation: Connect the two sub-blocks
Block Diagram Modeling: Analogy Approach Physical laws are used to predict the behavior (both static and dynamic) of systems. Electrical engineering relies on Ohm s and Kirchoff s laws Mechanical engineering on Newton s law Electromagnetics on Faradays and Lenz s laws Fluids on continuity and Bernoulli s law Based on electrical analogies, we can derive the fundamental equations of systems in five disciplines of engineering: Electrical, Mechanical, Electromagnetic, Fluid, and Thermal. By using this analogy method to first derive the fundamental relationships in a system, the equations then can be represented in block diagram form, allowing secondary and nonlinear effects to be added. This two-step approach is especially useful when modeling large coupled systems using block diagrams.
Power and Energy Variables: Effort & Flow [Ref.] Raul Longoria
Transitional Mechanical Systems Mechanical movements in a straight line (i.e. linear motion) are called transitional Force displacement Basic Blocks are: Dampers, Masses, and Springs Springs represent the stiffness of the system F kx Dampers (or dashpots) represent the forces opposing to the motion (i.e. friction) F cv velocity Masses represent the inertia F ma mass acceleration Dr. Tarek A. Tutunji
Transitional Mechanical Systems Equations for mechanical systems are based on Newton Laws Free body diagram Spring Mass Force Damper ma F - kx - c dx dt Dr. Tarek A. Tutunji
Example: Mass-Spring-Damper Note: D is Differentiation 1/D is Integration
Example: Two-Mass Mechanical System [Ref.] Prof. Shetty
Example: Two-Mass Mechanical System
Example: Mechanical Model Consider a two carriage train system x1 x Force Mass 1 k(x1-x) k(x1-x) Mass c v1 m m 1 x x 1 f - k(x k(x 1-1 - x x ) - cx ) - cx 1 c v
Example continued Taking the Laplace transform of the equations gives m m 1 s s X X 1 (s) (s) F(s) k(x - 1 k(x (s) - 1 (s) X - X (s)) (s)) - csx - csx Note: Laplace transforms the time domain problem into s-domain (i.e. frequency) L L x(t) X(s) x(t) sx(s) 0 e st x(t)dt (s) 1 (s)
Example continued Manipulating the previous two equations, gives the following transfer function (with F as input and V1 as output) V 1 (s) F(s) m 1 m s 3 c m s cs k m m s km km c s kc 1 1 Note: Transfer function is a frequency domain equation that gives the relationship between a specific input to a specific output
Example continued Simulation using MATLAB m1= 5; m=0.7; k=0.8; c=0.05; num=[m c k]; den=[m1*m c*m1+c*m k*m1+k*m+c*c *k*c]; sys=tf(num,den); % constructs the transfer function impulse(sys); % plots the impulse response step(sys); % plots the step response bode(sys); % plots the Bode plot
Example continued: Impulse response Dr. Tarek A. Tutunji
Example continued: Step response Dr. Tarek A. Tutunji
Dr. Tarek A. Tutunji Example continued: Bode Plot
Example: Motion of Aircraft
Rotational Mechanical Systems Consider a mechanical system that involves rotation Torque, T w dq/dt q T Shaft Side View J dw/dt kq cw The torque, T, replaces the force, F The angle, q, replaces the displacement x The angular velocity, w, replaces velocity v The angular acceleration, a, replaces the acceleration a The moment of inertia J, replaces the mass m Dr. Tarek A. Tutunji Top View
Rotational Mechanical Systems The mechanics equation becomes spring coefficient angle damper coefficient angular velocity moment of inertia angular acceleration Torque T kθ cω Jα dθ T kθ c J dt d dt θ Dr. Tarek A. Tutunji
Example: Rotational-Transitional System Consider a rack-and-pinion system. The rotational motion of the pinion is transformed into transitional motion of the rack w T in pinion r v F rack For simplicity, the spring effects are ignored dω Tin Tout J c1ω dt Dr. Tarek A. Tutunji
Example continued The rotational equation is dω Tin Tout J c1ω dt The transitional equation is F cv m dv dt Using the equations T out rf And manipulating the rotational and transitional equations with the input torque, Tin, as inputs and velocity, v, as output, we get ω v/r c T 1 in c r r v J r mr dv dt Dr. Tarek A. Tutunji
Example continued Let us take a look at the state space equations In general, where x is the states vector, y is the output vector, and u is the input vector x y Ax Bx Cu Du In our example, we will dω use the states: w and v, dv dt the inputs: T in and F dt the output: v v 0 Manipulating the equations in the previous slide, we get c1 J 0 ω v 1 c 0 ω v m 1 J 0 r J T 1 F m in
Conversion: Transitional and Rotational
Gear Trains
where Gear Trains
Electrical Systems: Basic Equations Resistor Ohm s Law Voltage V Ri current Inductor V L di dt Resistance Capacitor Power = Voltage x Current Inductance V i 1 idt C dv C dt Capacitance Dr. Tarek A. Tutunji
Kirchoff Laws Equations for electrical systems are based on Kirchoff s Laws 1. Kirchoff current law: Sum of Input currents at node = Sum of output currents. Kirchoff voltage law: Summation of voltage in closed loop equals zero
Example: RLC circuit R i L C V + - + V R - + V L - + V c - Using Kirchoff voltage law V Ri L di dt 1 C idt Or V Ri L di dt Vc since i C dv dt c Then V RC dv dt c LC d dt V c V A second order differential equation c
RLC MATLAB Code R=1000000; % R = 1MW L=0.001; % L=1 mh C=0.000001; % C= 1mF num=1; den=[l*c R*C 1]; sys=tf(num,den); bode(sys) Impulse(sys) Step(sys)
RLC Simulation: Bode Plot At DC (i.e. frequency = 0), Capacitor is open => Voltage gain is 0dB (i.e. 1 V/V) At high frequency, Capacitor is short => Voltage = 0 Dr. Tarek A. Tutunji
RLC Simulation: Impulse Response Input voltage is pulse => Capacitor stores energy And then releases the energy Dr. Tarek A. Tutunji
RLC Simulation: Step Response At about.3 seconds, the capacitor Voltage becomes 90% of the 1 Volts Dr. Tarek A. Tutunji
Op Amps
PM-DC Motor Modeling R i L V in + - Motor T w The electrical equation is V in Ri L di dt V emf where V emf (Back electromagnetic voltage) = k 1 w dω The mechanical equation is T J bω Tload dt where T = k i
DC Motor Model: Block Diagram T load V in i T + 1/(Ls+R) + - k - 1/(Js+b) w k 1
Simulation Result
Fluid Systems Fluid systems can be divided into two categories: Hydraulic: fluid is a liquid and incompressible Pneumatic: fluid is gas and can be compressed The volumetric rate of flow, q, is equivalent to the current The pressure difference, P 1 -P, is equivalent to voltage The basic building blocks for hydraulic systems are: Hydraulic resistance, capacitance, and inertance
Hydraulic resistance Hydraulic resistance is the resistance to the fluid flow which occurs as a result of valves or pipe diameter changes The relationship between the volume rate of flow, q, and pressure difference, p 1 -p,is given by Ohm s law p p 1 Rq p 1 p p 1 p pipe valve
Hydraulic Capacitance Potential energy stored in a liquid such as height of a liquid in a container q 1 h A p 1 p q Volume Change Volumetric rate of change q V 1 q Ah Cross sectional Area dv dt q 1 q dh A dt height
Hydraulic Capacitance p 1 p p hgρ density pressure height gravity Note that p F / A mg / A p Vg / A p hg dh q q A q1 - q dt d p gρ A dt 1 A gρ dp dt By letting the hydraulic capacitance be We get q q 1 C dp dt C A gρ
Hydraulic Inertance Equivalent to inductance in electrical systems To accelerate a fluid and increase its velocity a force is required F 1 F p1 p Cross sectional Area F mass 1 = p 1 A Length using F = p A Then p p 1 F1 F ma m ALρ q Av dq I dt dv m dt A Where the Inertance is Lρ I A
Hydraulic Example Modeling: an interactive -tank system q in (t) A h 1 (t) A 1 R 1 q h (t) 1 (t) dh1 dt dh dt q (t) q 1 (t) q q h h in 1 1 (t) q (t) q (t) h (t)/r 1 (t) (t) (t) /A /A /R 1 1 R q (t)
Hydraulic Example Modeling: Block Diagram Input: q in 1/A 1 - dh1 h h /dt 1 Sum Integrate Sum 1/R 1 1/A Sum Integrate - - q 1 1/A 1 1/A 1/R Dr. Tarek A. Tutunji Output: q
Hydraulic Example: Simulation Input, q in, is a step Output, q, is taken to a virtual scope Here, we assume all the Cross sectional areas and the resistances equals 1
Hydraulic Example: Simulation
Another Form of Analogies Potential and Flow Variables When systems are in motion, the energy can be Increased by an energy-producing source outside the system Redistributed between components within the system Decreased by energy loss through components out of the system. Therefore, a coupled system becomes synonymous with energy transfer between systems. Potential Variable = PV Flow Variable = FV
Analogies: FV and PV Flow Variable (FV) Potential Variable (PV) Electrical Current Voltage Mechanical Transitional Force Velocity Mechanical Rotational Torque Angular Velocity Hydraulic Volumetric Flow Rate Pressure Pneumatic Mass Flow Rate Pressure Thermal Heat Flow Rate Temperature
Which Analogies to use? Force-Voltage makes more physical sense Graphical Representation: Bond Graphs Force-Current makes mathematical sense Sum of Currents= Zero and Sum of Forces = Zero Graphical Representation: Linear Graphs
Conclusion Mathematical Modeling of physical systems is an essential step in the design process Simulation should follow the modeling in order to investigate the system response Mechatronic systems involve different disciplines and therefore an appropriate modeling technique to use is block diagrams Analogies among disciplines can be used to simplify the understanding of different dynamic behaviors Dr. Tarek A. Tutunji