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

Transcription

1 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

2 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

3 Block Diagram Example [Ref.] Prof. Shetty

4 Block Diagrams Manipulation

5 Block Diagrams Manipulation

6 Block Diagrams: Direct Method Example Consider the transfer function: We can introduce s state variable, x(t), in order to separate the polynomials

7 State Equation The differential equation is: Put the needed integrator blocks: Add the required multipliers to obtain the state equation:

8 Output Equation Repeat the same procedure for the output equation: Connect the two sub-blocks

9 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.

10 Power and Energy Variables: Effort & Flow [Ref.] Raul Longoria

11 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

12 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

13 Example: Mass-Spring-Damper Note: D is Differentiation 1/D is Integration

14 Example: Two-Mass Mechanical System [Ref.] Prof. Shetty

15 Example: Two-Mass Mechanical System

16 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 x x ) - cx ) - cx 1 c v

17 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)

18 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

19 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

20 Example continued: Impulse response Dr. Tarek A. Tutunji

21 Example continued: Step response Dr. Tarek A. Tutunji

22 Dr. Tarek A. Tutunji Example continued: Bode Plot

23 Example: Motion of Aircraft

24 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

25 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

26 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

27 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

28 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

29 Conversion: Transitional and Rotational

30 Gear Trains

31 where Gear Trains

32 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

33 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

34 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

35 RLC MATLAB Code R= ; % R = 1MW L=0.001; % L=1 mh C= ; % C= 1mF num=1; den=[l*c R*C 1]; sys=tf(num,den); bode(sys) Impulse(sys) Step(sys)

36 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

37 RLC Simulation: Impulse Response Input voltage is pulse => Capacitor stores energy And then releases the energy Dr. Tarek A. Tutunji

38 RLC Simulation: Step Response At about.3 seconds, the capacitor Voltage becomes 90% of the 1 Volts Dr. Tarek A. Tutunji

39 Op Amps

40 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

41 DC Motor Model: Block Diagram T load V in i T + 1/(Ls+R) + - k - 1/(Js+b) w k 1

42 Simulation Result

43 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

44 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

45 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

46 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ρ

47 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

48 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)

49 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

50 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

51 Hydraulic Example: Simulation

52 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

53 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

54 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

55 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