The objective of this lab is to review of the behavior of second-order systems. Students will use gain a better understanding of the importance system identification. Students will also develop a hands-on understanding of the concept of hardware gain and why it will play a crucial role in controller design. System Identification will be implemented in LabVIEW using the System Identification Toolkit.
System identification for the rectilinear plant
Objectives
Review of the behavior of second-order systems
Understand the importance of system identification
Understand the concept of hardware gain and why it will play a crucial role in controller design
Pre-lab
Derive the equations of motion for the 1DOF spring-mass system with friction damping shown below and find the transfer function. Simulate the unforced, natural response in LabVIEW to an initial displacement of
$3.0cm$ (look for the Initial Response function in the Control Design Toolkit). Determine the natural frequency from the plot. Use the following parameters:
$M=3.0kg$$k=350N/m$$c=4.0\mathrm{Nsec}/m$
Lab procedure
Power on the PXI, and while it's booting up configure the plant as explained in the next step.
Clamp the second and third mass carriages, and attach the medium stiffness spring between the first and second carriage as shown in Fig. 2. Secure four
$0.5kg$ brass weights to the first and second carriages.
Start LabVIEW on your host PC and target the RT system.
Open the
Lab 2 - Rectilinear System Identification.vi file, enter a loop time of 0.005 seconds, and then run the VI. Also, turn on the amplifier now.
Reset the encoders if necessary and enter the parameters to execute a zero amplitude step with a dwell time of
$4000ms$ . Prepare to manually displace the first mass carriage. When you are ready, execute the step command, and manually displace the first carriage approximately
$3cm$ ,
then release it. Observe the natural response of the carriage, and wait for the system to finish executing the trajectory and acquiring the data.
Recall that for underdamped systems, the damped natural frequency is related to the undamped natural frequency through the following relationship:
For the purposes of plant identification, you may assume that the damping due to friction in the system is sufficiently small to satisfy the above condition.
Plot the encoder 1 position data and calculate the natural frequency of this system in both
$Hz$ and
$rad/\mathrm{sec}$ . Call this value
${\omega}_{{n}_{\mathrm{m11}}}$ (the subscript
$\mathrm{m11}$ denotes mass #1 trial #1). Use only oscillations with amplitude greater than 1000 counts in your calculation. This is because smaller amplitudes start to become dominated by nonlinear friction effects.
Remove the four brass weights and repeat steps 5 and 6 to obtain
${\omega}_{{n}_{\mathrm{m12}}}$ for the unloaded carriage. Decrease the execution time if necessary.
Measure the initial cycle amplitude
${X}_{0}$ and the last cycle amplitude
${X}_{n}$ for the n cycles measured in Step 6. Using relationships associated with the logarithmic decrement find the damping ratio (call it
${\zeta}_{\mathrm{m12}}$ ) as
Configure the plant as shown in Fig. 2b where the first mass carriage is now clamped and the second is free. Repeat steps 5 through 8 for the second mass carriage to obtain
${\omega}_{{\mathrm{n}}_{\mathrm{m21}}}$ ,
${\omega}_{{\mathrm{n}}_{\mathrm{m22}}}$ and
${\zeta}_{\mathrm{m22}}$ . How does this damping ratio compare with that for the first mass? Be sure to save this plotted data as it will be used in the next experiment.
Connect the mass carriage extension bracket and dashpot to the second mass as shown in Fig. 2c, and load four brass weights onto the second mass carriage. Open the damping (air flow) adjustment knob 2.0 turns from the fully closed position. Perform the necessary steps to obtain
${\zeta}_{d}$ where the "d" subscript denotes the dashpot. For this calculation use amplitudes
$\ge $ 500 counts.
Each brass weight has a mass of
$0.5\pm 0.01kg$ . (You may weigh the pieces if a more precise value is desired). Call the mass of the four weights combined
${m}_{w}$ . Recall that the dynamics for a second-order system in terms of the Laplace variable
$s$ is given by
$${s}^{2}+2\zeta {\omega}_{n}s+{\omega}_{n}^{2}$$ Compare the above with the equations of motion you derived in terms of mass, damping coefficient, and spring constant (
$m$ ,
$c$ , and
$k$ ). Use your knowledge of second-order systems along with your recorded values of
${\omega}_{{n}_{\mathrm{m11}}}$ and
${\omega}_{{n}_{\mathrm{m12}}}$ to solve for the unloaded carriage mass
${m}_{\mathrm{c1}}$ , and spring constant
$k$ .
$$$$ Repeat these calculations for the trials involving the second mass carriage and spring, mc2 and k.
$$$$ Since the same spring was used in the experiments, the values for
$k$ that result from the above calculations should be very close. You may use the average of the two for all future calculations and experiments. Call this value
${k}_{med}$ to denote the medium stiffness spring.
$$$$ Now with your recorded values for
$\zeta $ calculate
${c}_{1}$ , the damping coefficient for the first mass carriage.
$$$$ Repeat this calculation to find the damping coefficient for the second mass carriage,
${c}_{2}$ .
$$$$ Finally, calculate the damping coefficient of the dashpot,
${c}_{d}$ .
Remove the carriage extension bracket and dashpot from the second mass carriage, replace the medium stiffness spring with the high stiffness spring, and perform the necessary steps to obtain the resulting natural frequency
${\omega}_{\mathrm{m23}}$ . Repeat this frequency measurement using the least stiff spring to obtain
${\omega}_{\mathrm{m24}}$ .Now you can compute
${k}_{\mathrm{high}}$ and
${k}_{low}$ .
$$$$ All dynamic parameters have been identified! Values for
${m}_{\mathrm{c1}}$ and
${m}_{\mathrm{c2}}$ for any configuration of masses may be found by adding the calculated mass contribution of the weights to that of the unloaded carriages. Be sure that you have filled in all of the system parameters in the table provided before proceeding.
The following is necessary to establish the hardware gain for control modeling purposes.
Remove the spring connecting the first and second carriages and secure four
$0.5kg$ brass weights onto the first carriage. Using the stop blocks, secure the second mass carriage clear from the first, and remove any stop blocks that are constraining the first carriage. Verify that the masses are secure and that the carriage slides freely. Position the first mass approximately
$3cm$ to the left (negative
${x}_{1}$ position) of its center of travel.
Perform a
$2.0V$ bidirectional step with a dwell time of
$75ms$ (recall that on the PXI,
$\pm {2}^{15}$ counts corresponds to a voltage of
$\pm 10V$ ).
Plot the velocity data and observe the four segments of the profile with nominal shapes of: linear increase (constant acceleration), nearly constant velocity (zero acceleration), linear decrease (constant deceleration), and constant. Obtain the acceleration from the linear, positive-sloped segment of the velocity profile. Repeat this for the negative-sloped segment. You may call these two values
${\stackrel{..}{x}}_{1{e}_{acc}}$ and
${\stackrel{..}{x}}_{1{e}_{dec}}$ , respectively. Be sure that you account for the scaling of velocity data that is shown on the front panel graph indicator, and remember to include units in your calculation.Calculate the average of the magnitudes of the positive and negative accelerations (call this value
${\stackrel{..}{x}}_{1e}$ ) to be used in obtaining
${k}_{hw}$ below.
Recall that the hardware gain for this system is given by
$${k}_{hw}={k}_{c}{k}_{a}{k}_{t}{k}_{mp}{k}_{ep}{k}_{e}$$ where
$$$$${k}_{c}$ , the DAC gain, =
$\mathrm{10}V/32,768$ DAC counts
$$$$${k}_{a}$ , the servo amp gain, = to be determined
$\mathrm{Amps}/V$$$$$${k}_{t}$ , the servo motor torque constant, = to be determined
$Nm/Amp$$$$$${k}_{mp}$ , the motor pinion pitch radius inverse, =
$26.25{m}^{-1}$$$$$${k}_{ep}$ , the encoder pinion pitch radius inverse, =
$89{m}^{-1}$$$$$${k}_{e}$ , the encoder gain, =
$16,000$ pulses/2
$\pi $ radians
$$$$ Using what you have learned about hardware gain, calculate a
value for the product
${k}_{a}{k}_{t}$ . Use this value to compute the hardware gain
${k}_{hw}$ .
Save any files or plots of interest. Stop and exit the VI,
power off the amplifier and PXI, and return all plant materials tothe instructor.
Post-lab
Complete the table below making sure to include units.
How do your values of
${m}_{\mathrm{c1}}$ and
${c}_{1}$ compare with
${m}_{\mathrm{c2}}$ and
${c}_{2}$ . Can
you explain why they may be different?
What were your values for
${\zeta}_{\mathrm{m12}}$ and
${\zeta}_{\mathrm{m22}}$ Are these values
sufficiently small for the approximations made in steps 6 and 8 tobe valid?
Explain the importance of accurately identifying the plant
parameters and hardware gain in terms of how they will affectcontroller design.
Questions & Answers
anyone know any internet site where one can find nanotechnology papers?
fullerene is a bucky ball aka Carbon 60 molecule. It was name by the architect Fuller. He design the geodesic dome. it resembles a soccer ball.
Tarell
what is the actual application of fullerenes nowadays?
Damian
That is a great question Damian. best way to answer that question is to Google it. there are hundreds of applications for buck minister fullerenes, from medical to aerospace. you can also find plenty of research papers that will give you great detail on the potential applications of fullerenes.
Tarell
what is the Synthesis, properties,and applications of carbon nano chemistry
Yeah, it is a pain to say the least. You basically have to heat the substarte up to around 1000 degrees celcius then pass phosphene gas over top of it, which is explosive and toxic by the way, under very low pressure.
Harper
Do you know which machine is used to that process?