Problems

Problem 4.1 (TRIANGLE)

Respond to the following questions and imperatives with one or two sentences and, if needed, equations and/or sketch.

  1. Why do we include a resistor in lumped-parameter motor models?

  2. How are brushes used in brushed DC motors?

  3. With regard to standard motor curves, why do we say the “braking power” is equivalent to the power that could be successfully transferred by the motor to the mechanical system?

  4. In terms of electrical and mechanical processes, why does an efficiency versus torque motor curve have a peak?

  5. As a DC motor’s bearings wear down, how will its efficiency curve be affected?

  1. There is no physical resistor in a motor, but its rotor and (if applicable) stator windings have a small amount of resistance. This resistance is important to model both because it is the dominant factor in power losses in a motor and because it can thereby generate such heat that damage to the windings is a concern.

  2. Brushes reverse the direction of current flow through the rotor windings and enable electrical contact between stationary and a rotating conductors.

  3. In this context, “braking power” refers to the power dissipated by the brake in the experimental setup used to measure motor curves. Although this power is dissipated in the experiment, it is potentially delivered to a load (in steady state) in application.

  4. The efficiency motor curve has a peak because there are two processes dissipating power: the electrical resistance and the mechanical rotational damping, which dominate separately at different ends of the operating spectrum. Specifically, at high-speed and low-torque, the mechanical damping dominates; conversely, at low-speed and high-torque, the electrical resistance dominates. In between these two ends of the spectrum, the trade-off peaks.

  5. As a DC motor’s bearings wear down, the mechanical damping energy dissipation will increase, lowering the overall efficiency curve and shifting the peak toward higher-torque and lower-speed, due to the trade-off described in the previous part of this solution.

Problem 4.2 (SQUARE)

schematic of an electromechanical system for .

Consider the system presented in the schematic of . Let the DC motor have motor constant Ka (units N-m/A) and let the motor be driven by an ideal current source IS. Assume the motor inertia has been lumped into J1 and motor damping lumped into B1.

  1. Draw a linear graph model.

  2. Draw a normal tree.

  3. Identify any dependent energy storage elements. If the motor was driven by an ideal voltage source instead, how would this change?

Draw a linear graph model and normal tree.

The linear graph (a) and normal tree (b) are shown in .

a normal tree on a linear graph model of the electromechanical system of .

We had no choice but to include the inductor L in the normal tree because we cannot include IS, needed the tree to reach the node between R and IS (therefore R must be in the tree), needed the tree to be connected (therefore L must be in the tree), and also needed the tree to reach the electrical ground (therefore 1 must be in the tree).

Since L is in the normal tree and it is a T-type energy storage element, it is dependent.

This would change if the source was instead a voltage source VS, which can and should be included in the normal tree, which now can and must exclude L (so as not to form a loop). Therefore, L would be an independent energy storage element in this case.

Problem 4.3 (RECTANGLE)

Consider the system presented in the schematic of . From the linear graph model and normal tree derived in , derive a state-space model in standard form. Let the outputs be θJ1 and θJ2, the angular positions of the flywheels.

We can use the results of the solution to , including its linear graph and normal tree of .

We begin by identifying the state, input, and output vectors: x=[ΩJ1Tk1ΩJ2Tk2],u=[IS],y=[θJ1θJ2]. Note that iL is not a state variable because L is not an independent energy storage element.

Now we write the elemental equation and continuity or compatibility equation for each element in the following table.

Now we can eliminate all secondary variables via substitution of the continuity or compatibility equation corresponding to each element into its elemental equation, which results in the following table.

element
J1 dtΩJ1=1J1(Tk1+TB1+T2)
k1 dtTk1=k1(ΩJ1ΩJ2)
J2 dtΩJ2=1J2(Tk1TB2Tk2)
k2 dtTk2=k2ΩJ2
B1 TB1=B1ΩJ1
B2 TB2=B2ΩJ2
1 T2=KaIS
2 v1=KaΩJ1
L vL=LdtIS
R vR=ISR

The first four equations yield the state equations when we substitute some of the last six.

element
J1 dtΩJ1=1J1(Tk1+B1ΩJ1+KaIS)
k1 dtTk1=k1(ΩJ1ΩJ2)
J2 dtΩJ2=1J2(Tk1B2ΩJ2Tk2)
k2 dtTk2=k2ΩJ2

The standard form of the state equation is dxdt=[B1/J11/J100k10k1001/J1B2/J21/J200k20]x+[Ka/J1000]u. The output angular positions can be related to the spring torques (which are state variables) as follows: θJ1=θk1+θk2=Tk1/k1+Tk2/k2 and θJ2=θk2=Tk2/k2. Therefore the output equation in standard form is y=[01/k101/k20001/k2]x+[00]u.

Problem 4.4 (QUADRILATERAL)

a linear graph model of the electromechanical system.

Consider the linear graph model of a motor coupled to a rotational mechanical system shown in . This is similar to the model from the , but includes the flexibility of the shaft coupler. An ideal voltage source drives the motor—modeled as an ideal transducer with armature resistance R and inductance L, given by the manufacturer in . The ideal transducer’s rotational mechanical side (2) is connected to a moment of inertia Jm modeling the rotor inertia and damping Bb modeling the internal motor damping, both values given in the motor specifications. Take Bb=Bm and Jf=0.324103 kg-m2. Assume the shaft coupling has a torsional stiffness of k=100 N-m/rad.

  1. Derive a state-space model for the system with outputs i1 and ΩJf.

  2. Create a Matlab model for the system and simulate its response from rest to an input voltage VS=10 V.

  3. Plot the outputs through time until they reach steady state.

a linear graph model with normal tree of the electromechanical system.

The normal tree is shown in . The state variables are ΩJm, ΩJf, iL, and Tk, so the system order is n=4. The state, input, and output vectors are (different orderings are valid) x=[ΩJmΩJfiLTk]u=[VS]y=[i1ΩJf]

The following are the elementary, continuity, and compatibility equations.

We now use StateMint: statemint.stmartin.edu The .rnd file is at the following url: ricopic.one/assets/quadrilateral.rnd

StateMint gives a state-space model that has standard form dxdt=Ax+Buy=Cx+Du with A=[BmJm0KaJm1JmBbJf001JfKaL0RL0kk00]B=[001L0]C=[00100100]D=[00].

Matlab

We now turn to Matlab.

First, we define each of the parameters. The following, symbolic approach that stores paramter values in a params structure array is a little more work, but when we start designing, it will pay off to hold off on defining parameters as long as possible, so we get in the practice, now.

syms R L Ka k Bm Jm Bb Jf

Kv_spec = 10.2; % V/krpm ... voltage constant spec 
rads_krpm = 1e3*2*pi/60; % (rad/s)/krpm
Kv_si = Kv_spec/rads_krpm; % V/(rad/s)
params.Ka = Kv_si; % N-m/A ... TF/motor constant
m_in = 0.0254; % m/in
d = 2.5*m_in; % m ... flywheel diameter
thick = 1*m_in; % m ... flywheel thickness
vol = pi*(d/2)^2*thick; % flywheel volume
rho = 8000; % kg/m^3 ... flywheel density (304) 
m = rho*vol; % kg ... flywheel mass
params.Jf = 1/2*m*(d/2)^2; % kg-m^2 ... flywheel in.
params.Jm = 56.5e-6; % kg-m^2 ... inertia of rotor
params.Bm = 16.9e-6; % N-m/s^2 ... motor damping coef
params.Bb = params.Bm; % N-m/s^2 ... brng damp coef
params.k = 100; % N-m/rad ... torsional stiffness
params.R = 1.6; % Ohm ... armature resistance
params.L = 4.1e-3; % H ... armature inductance

Now we can define the state-space model standard matrices.

A = [...
  -Bm./Jm 0 Ka./Jm -1./Jm; ...
  -Bb./Jf 0 0 1./Jf; ...
  -Ka./L 0 -R./L 0; ...
  k -k 0 0 ...
];
B = [0; 0; 1./L; 0];
C = [0 0 1 0; 0 1 0 0];
D = [0; 0];

Now we can create an LTI system model using ss. Note that Matlab doesn’t support symbolic LTI models, so we have to substitute our parameters.

A_s = double(subs(A,params)); % sym2num
B_s = double(subs(B,params)); % sym2num
C_s = double(subs(C,params)); % sym2num
D_s = double(subs(D,params)); % sym2num
sys = ss(A_s,B_s,C_s,D_s);

Now we can simulate with lsim (or even step in this case).

t = linspace(0,.3,400);
u = 10*ones(size(t)); % V ... step input
y = lsim(sys,u,t);

Now let’s plot the results.

figure
[ax,l1,l2] = plotyy(t,y(:,1),t,y(:,2));
set(l1,'linewidth',1.5); set(l2,'linewidth',1.5);
ylabel(ax(1),'i_1 (A)')
ylabel(ax(2),'\Omega_{J_f} (rad/s)')
xlabel('time (s)')

The “wiggle” of is a contribution of the flexible shaft coupler. It does damp out over time – try plotting (with more points) up to 1 sec – but relatively slowly. The contribution is small, but potentially not negligible.

Note the two vertical axes in the plot. This is good form when we plot two outputs on the same plot because they generally have different physical units! (In our case, they are current with units A and angular velocity with units rad/s.) It is challenging to plot more than two quantities with different units on the same plot, but such an example is the motor curves of . Note that Matlab recommends using yyaxis now, but plotyy will continue to work, although it is more limited.

Problem 4.5 (MRPOTATOHEAD)

Consider the linear graph model (with normal tree) of . This is a model of a motor with constant Ka connected to a pair of meshing gears with transformer ratio N, the output over input gear ratio. An ideal voltage source drives the motor—modeled as an ideal transducer with armature resistance R and inductance L. The motor’s rotational mechanical side (2) is connected to a moment of inertia J2 modeling the rotor and drive gear combined inertia. The damping element B2 models the internal motor damping and the drive gear bearing damping. The output side of the gear transducer (4) is connected to a moment of inertia J4 modeling the output gear and load combined inertia. The damping element B4 models the internal motor damping and the drive gear bearing damping. Use the parameter values given in .

  1. Derive a state-space model for the system with outputs ΩJ2 and ΩJ4.

  2. Create a Matlab model for the system and simulate its response from rest to an input voltage VS=20 V.

  3. Plot the outputs through time until they reach steady state.

R 2 Ω
L 8 mH
Ka 0.2 N-m/A
J2 0.1103 kg-m2
B2 50 μN-m/(rad/s)
N 5
J4 1103 kg-m2
B4 70 μN-m/(rad/s)
A linear graph model with normal tree in green of an electromechanical system with a gear reduction.
Figure 4.1: A linear graph model with normal tree in green of an electromechanical system with a gear reduction.

The state variables are across-variables on A-type branches and through-variables on T-type links, so ΩJ4 and iL, so the system order is n=2.1 The state, input, and output vectors are (different orderings are valid) x=[ΩJ4ΩiL]u=[VS]y=[ΩJ2ΩJ4]

The following are the elementary, continuity, and compatibility equations.

We now use StateMint: statemint.stmartin.edu The .rnd file is at the following url: ricopic.one/dynamic_systems/source/mrpotatohead.rnd

StateMint gives a state-space model that has standard form dxdt=Ax+Buy=Cx+Du with A=[B2N2+B4J2N2+J4KaNJ2N2+J4KaNLRL]B=[01L]C=[N010]D=[00].

Matlab

We now turn to Matlab.

First, we define each of the parameters. The following, symbolic approach that stores paramter values in a params structure array is a little more work, but when we start designing, it will pay off to hold off on defining parameters as long as possible, so we get in the practice, now.

syms R L Ka B2 J2 N B4 J4

params.R = 2; % Ohm ... armature resistance
params.L = 8e-3; % H ... armature inductance
params.Ka = 0.2; % N-m/A ... motor constant
params.J2 = .1e-3; % kg-m^2 ... inertia
params.B2 = 50e-6; % N-m/(rad/s) ... damping
params.N = 5; % gear ratio
params.J4 = 1e-3; % kg-m^2 ... inertia
params.B4 = 70e-6; % N-m/(rad/s) ... damping

Now we can define the state-space model standard matrices.

A = [...
  -(B2.*N.^2 + B4)./(J2.*N.^2 + J4), ...
  Ka.*N./(J2.*N.^2 + J4); ...
  -Ka.*N./L, ...
  -R./L ...
];
B = [0; 1./L];
C = [N 0; 1 0];
D = [0; 0];

Now we can create an LTI system model using ss. Note that Matlab doesn’t support symbolic LTI models, so we have to substitute our parameters.

A_s = double(subs(A,params)); % sym2num
B_s = double(subs(B,params)); % sym2num
C_s = double(subs(C,params)); % sym2num
D_s = double(subs(D,params)); % sym2num
sys = ss(A_s,B_s,C_s,D_s);

Now we can simulate with lsim (or even step in this case).

t = linspace(0,.23,300);
u = 20*ones(size(t)); % V ... step input
y = lsim(sys,u,t);

Now let’s plot the results.

figure
[ax,l1,l2] = plotyy(t,y(:,1),t,y(:,2));
set(l1,'linewidth',1.5); 
set(l2,'linewidth',1.5,'linestyle','--');
set(ax(1),'YLim',[0,120])
set(ax(2),'YLim',[0,24])
ylabel(ax(1),'\Omega_{J_2} (rad/s)')
ylabel(ax(2),'\Omega_{J_4} (rad/s)')
xlabel('time (s)')

Problem 4.6 (CLUNKER)

Draw a linear graph, a normal tree, identify state variables, identify system order, and denote any dependent energy storage elements for each of the following schematics.

  1. The electronic system of figure 4.21, voltage and current sources, and transformer with transformer ratio N.

  2. The electromechanical system of figure 4.22 with motor model parameters shown, coordinate arrow in green. Model the propeller as a moment of inertia J2 and damping B2.

  3. The translational mechanical system of figure 4.23, force source, coordinate arrow in green.

a circuit diagram.
Figure 4.21: a circuit diagram.
Sketch of a motor coupled to a fan.
Figure 4.22: Sketch of a motor coupled to a fan.
Schematic of a mechanical system.
Figure 4.23: Schematic of a mechanical system.

problem 4.6 part (a) linear graphs and normal trees in green – four possible choices!{#clunker4 .figure h=“clunker4” caption_plain=” part (a) linear graphs and normal trees in green – four possible choices!“}

problem 4.6 part (b) linear graphs and normal tree in green
Figure : problem 4.6 part (b) linear graphs and normal tree in green
problem 4.6 part (c) linear graphs and normal tree in green
Figure : problem 4.6 part (c) linear graphs and normal tree in green
  1. The electronic system of figure 4.21 has linear graph and normal trees of . The state variables are vC2 and either iL1 or iL2. The system order is n=2. Capacitor C1 and either L1 or L2 are dependent energy storage elements.

  2. The electromechanical system of figure 4.22 has linear graph and normal trees of . The state variables are ΩJm, ΩJ2, iL, and Tk1. The system order is n=4. There are no dependent energy storage elements.

  3. The translational mechanical system of figure 4.23 has linear graph and normal trees of . The state variables are vm1, vm2, vm3, fk1, and fk2. The system order is n=5. There are no dependent energy storage elements.

Problem 4.7 (CURVY)

Consider the DC motor curves of figure 4.13, reproduced in figure 4.24.

  1. At peak efficiency, what is the steady-state motor speed?

  2. At peak efficiency, what is the steady-state motor torque?

  3. You are to use this motor to drive a load at a constant angular speed of 100 rad/s with at least 1 N-m of torque. You wisely choose to use a gear reduction between the motor and load. What should the gear ratio be to meet the above requirements and optimize efficiency? Justify your answer in terms of the motor curves of figure 4.24.

Figure 4.24: The motor curve figure 4.13.
  1. From the motor curves, peak efficiency corresponds to about 580 rad/s.

  2. From the motor curves, 0.2 N-m.

  3. The gear reduction allows us to trade speed and torque. Without the gear reduction, you can easily meet the 100 rad/s requirement with ample torque, around 3 N-m. However, in those operating conditions, efficiency is around 20 percent (due to resistive electronic dissipation, mostly). In order to optimize efficiency, we’d like to operate as closely as possible to 0.2 N-m and 580 rad/s. If we gear down the full 580 to 100 rad/s (N=5.8), we get a torque of 0.25.8=1.16 N-m. This is more than enough torque to meet our 1 N-m torque requirement, therefore N5.8 is the optimal gear ratio.

Problem 4.8 (CHAIR)

Consider the opamp circuit of figure 4.25, which will be used to drive a PMDC motor. The input can supply a variable VS[0,10] V, the motor has constant Ka=0.05 V/(rad/s) and coil resistance Rm=1 Ω, and the opamp has differential supplies ±24 V. Assume the maximum torque magnitude required from the motor at top speed is |T2|=0.1 N-m and ignore any voltage drop in the motor due to the coil inductance.2 Select R1 and R2 to demonstrably meet the following design requirements:

  1. drivable motor speeds of at least [0,400] rad/s,

  2. no saturation of the opamp (i.e. |vo|<24 V), and

  3. a maximum combined power dissipation by R1 and R2 less than 300 mW.

An opamp circuit.
Figure 4.25: An opamp circuit.

Hint: start with the elemental equations of the DC motor to determine the necessary amplifier output vo, then constrain R1 and R2 to meet the gain requirements, and finally further constrain R1 and R2 to meet the power dissipation requirement.

The motor has elemental equations Ω2=v1/Kaand T2=Kai1. The first elemental equation gives at top speed v1=KaΩ2=(0.05)(400)=20 V. And here we need a current that meets the torque requirement i1=T2/Ka=0.1/.05=2 A, which corresponds to a resistive coil voltage drop of vRm=(2)(1)=2 V. So we need 22 V at top 400 rad/s. This is within the saturation of the opamp (24 V), so if vo=22 V, we meet all requirements so far.

We recognize the non-inverting opamp circuit, which delivers an output voltage vo=R1+R2R2VS as long as it doesn’t saturate (|vo|<24 V). To produce a top vo=22 V with max VS=10 V, the resistance ratio must be R1+R2R2=22/10. The last requirement limits power dissipation by R1 and R2 to 300 mW. Since these share a current, they are effectively in series and the total voltage drop across them is \tag not allowed in aligned environment From the circuit diagram, vR12=vo. The power dissipated by the resistors then is \tag not allowed in aligned environment which gives the constraint vo2R1+R2<0.3 WR1+R2>2220.31.613 kΩ. Combining this with the yields 2.2R2>1613R2>733.3 Ω or R1>880.0 Ω and R1=1.2R2. So, for instance, a valid solution is R1,R2=12,10 kΩ. So it turned out to be relatively easy to meet the power dissipation requirement by selecting large enough resistors.

Problem 4.9 (ONOMATOPOEIA)

Consider the DC motor curves of figure 4.13, reproduced in figure 4.24. If this motor is running at 400 rads,

  1. How much torque is produced?

  2. What is the output power?

  3. What is the input power?

  4. Why are the input and output power the same or different?

  1. At 400 rads the motor produces approximately 1.3 Nm of torque.

  2. At 1.3 Nm the output power is approximately 520 W

  3. At 1.3 Nm the efficiency is approximately 65%. This means the output power is 65% lower than the input power. Therefore, 520 W=0.65Pin, so Pin=5200.65=800 W.

  4. The output power is lower than the input power because the resistance of the motor armature and friction of the bearings cause energy to be lost. This ends up causing the motor to heat up as it is run.

Problem 4.10 (DEGLAZIFICATION)

Explain in your own words what lumped parameter elements should be used when modeling an electric motor and why.

Problem 4.11 (CONFUZZLED)

In the linear graph below a system is depicted consisting of a motor with its related damping and inertia driven by a voltage source and connected to a set of gears driving a second inertia. A rotary spring is attached between the two inertias.
Given this linear graph:

  1. draw a normal tree,

  2. determine the state variables and system order, and

  3. list any dependent energy storage elements and explain what this implies.

Problem 4.12 (LEVITATION)

In the linear graph and normal tree below a system is depicted consisting of a motor driven by a voltage source VS with inertia J driving a rotary damper and spring connected in series. Let the motor constant be Ka, and outputs of the system be the rotational velocity of the inertia, ΩJ, and the change in rotational velocity across the rotational damper, ΩB.
Given this linear graph and normal tree:

  1. determine the state variables,

  2. define the state, input, and output vectors,

  3. write the elemental, continuity, and compatibility equations, and

  4. solve for the state and output equations.


  1. Do not ignore the voltage drop across Rm, though. Note that this amounts to an assumption of steady-state operation at top speed. By requiring a specific T2, we are also implicitly ignoring torque losses due to motor bearing damping.↩︎

  2. Do not ignore the voltage drop across Rm, though. Note that this amounts to an assumption of steady-state operation at top speed. By requiring a specific T2, we are also implicitly ignoring torque losses due to motor bearing damping.↩︎

Online Resources for Section 4.10

No online resources.