img
Comments...

img PD Control Example

A  common control scheme is called PID or Proportional-Derivative-Integral control. PID controller are used in everything from process controls to motor controls. The following example illustrates the use of a PD, Proportional-Derivative, control or a rotating mass.

6.1 Physical Model

A simple DC motor model is shown below. It consists of a rotating mass, an input torque and viscous friction.

image/svg+xmlTorqueθJFriction
where J=10
and b=20
The input is torque (T) and the output is motor position (θ)

The system transfer can be found as follows:

J
d2θ_________dt2
= ∑ Torque
J
d2θ_________dt2
= T - b
d θ_________dt

d2θ_________dt2
+
b_________J
_________dt
=
1_________J
T
substituting in value for b and J gives
d2θ_________dt2
+2
_________dt
= 0.1 T
and the transfer function is
T_________θ
=
0.10_________s2 + 2 s

We wish to control the motor position (θ) and have a damping ratio (ζ) of 0.707

6.2 Open-Loop Analysis

The open characteristics of the system, shown below, don't meet the design specification. First the system does not have a ζ of 0.707. Second, the output is an integral of the input (pole at the origin). We wanted the output to follow the input.

image/svg+xmlReImXX-2
Fig 6.1 Open-Loop Analysis  imgimgimgimgimg0:00

6.2 Proportional Control

Use proportional feedback to control the motor position. Proportional control has the form
T = kpdesired - θ)

Substituting into the differential equation for the system gives

d2θ_________dt2
+2
_________dt
= 0.1 kpdesired - θ)
d2θ_________dt2
+2
_________dt
+ 0.1 kp θ = 0.1 kp θdesired

This can also be done using block diagrams and transfer functions

image/svg+xmlkp_+θdesiredθ0.1s2+2sT
where θ is the measured position, θd is the desired position and kp is the feedback gain.

Now simplify the transfer function

G(s)_________1+G(s)
=
0.1 kp_________s2 + 2 s
─────────
1+
0.1 kp_________s2 + 2 s
=
0.1 kp_________s2 + 2 s + 0.1 kp

Note that when analyzing the closed loop system using root locus the characteristic equation is

s2 + 2 s + 0.1 kp
which d(s) = s2 + 2 s and n(s) = 0.1

We can find kp using either root locus or, if we know the desired closed-loop poles, we can find kp directly through algebra.

Finding kp Using Root Locus

Given the closed loop characteristic equation

s2 + 2 s + 0.1 kp = 0
where d(s) = s2 + 2 s and n(s) = 0.1
plot the poles (roots of d(s)) and zeros (roots of n(s)) as 0 ≤ kp ≤ ∞.
image/svg+xmlReImXX-2-1Ts = 4/1 = 4sζ=0.707-1±idesired pole locationXX

From all of the possible solutions choose kp that results in ζ=0.707. Do this using brute force by selecting a kp and finding the roots of the resulting quadratic equation. As you change kp the root locations are indicated with an x on the root locus plot.

s2 + 2 s + 0.1 kp = 0    where   kp =
0.000

roots are: s = 0, -2

Solving for kp directly

If you know the form of the characteristic equation for the closed-loop system, you can often solve directly for kp. In our problem closed-loop characteristic equation is:

s2 + 2 s + 0.1 kp=0

Compare this to the standard 2nd order form

s2 + 2 ζ ωn + ωn2 = 0

For the two equations to be equal each terms in the two equations must match. Set ζ=0.707, match terms and solve:

2 = 2 ζ ωn=2 (0.707)ωn
and
0.1 kp = ωn2
Solve for kp
ωn = 2/(2 (0.707)) = 1.414
and
0.1 k = ωn2 k = 1.4142/.1 = 20.0

Alternately if we know the desired close-loop pole locations we can solve for kp directly. The key is that you must know the exact value for the pole locations and those must lie on the root locus. The desired pole locations from the root-locus graph are: s= -0.025±-0.025 i. Therefore the desired closed-loop equation must be

(s + 1 + 1 i)(s + 1 - 1 i)=0
s2 + 2 s + 2
The actual closed-loop characteristic equation is
s2 + 2 s + 0.1 kp
matching terms with the closed-loop equation yields
s2 = s2
2 s = 2 s
2 = 0.1 kp     kp = 20

Verifying the Solution

Use the following simulation to verify the solution. Select a value for k and change the value for θdesired. The settling time and damping ratio for the system should match the calculate settling time (remind me). Experiment with different value for k and verify that you calculated closed-loop poles (Tsettling and ζ) match the simulated response.

When kp = 20 the system has closed-loop poles at -1 ± 1 i. This corresponds to a settling time of
Tsettling = -4/(real part of pole location)=4 sec

The damping ratio (ζ) of 0.707 is a response that has a slight overshoot but does not oscillate.

In addition to tracking θdesired the system also rejects disturbances to the motor. Click and drag on the wheel. Note how the torque increases to resist the disturbance.

0
°
kp =
20.0

The final check in the design is to verify that the motor torque requirements don't exceed the physical limits of the system. Starting at 0° set θdesired to 90° (or some value that represents the maximum change in position) and observe the maximum torque required by the control system. If this value exceeds the physical limits of the system you may need to reduce kp, limit the maximum change in the system, or limit how fast the system is allowed to change (this technique is discussed later in this text).