# 0.4 System identification for the rectilinear plant

 Page 1 / 1
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.

## 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:
${\omega }_{d}={\omega }_{n}\sqrt{1-{\zeta }^{2}}$
and it follows that
${\omega }_{n}=\frac{{\omega }_{d}}{\sqrt{1-{\zeta }^{2}}}$
Therefore, for light damping (small values of ζ)
${\omega }_{n}\approx {\omega }_{d}$
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 }_{m12}$ ) as
$\frac{\zeta }{\sqrt{1-{\zeta }^{2}}}=\frac{1}{2\mathrm{pi}n}ln\frac{{X}_{0}}{{X}_{n}}-\zeta \approx \frac{1}{2\mathrm{pi}n}ln\frac{{X}_{0}}{{X}_{n}}$
• 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 }_{{n}_{m21}}$ , ${\omega }_{{n}_{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±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, $±{2}^{15}$ counts corresponds to a voltage of $±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

Is there any normative that regulates the use of silver nanoparticles?
what king of growth are you checking .?
Renato
What fields keep nano created devices from performing or assimulating ? Magnetic fields ? Are do they assimilate ?
why we need to study biomolecules, molecular biology in nanotechnology?
?
Kyle
yes I'm doing my masters in nanotechnology, we are being studying all these domains as well..
why?
what school?
Kyle
biomolecules are e building blocks of every organics and inorganic materials.
Joe
anyone know any internet site where one can find nanotechnology papers?
research.net
kanaga
sciencedirect big data base
Ernesto
Introduction about quantum dots in nanotechnology
what does nano mean?
nano basically means 10^(-9). nanometer is a unit to measure length.
Bharti
do you think it's worthwhile in the long term to study the effects and possibilities of nanotechnology on viral treatment?
absolutely yes
Daniel
how to know photocatalytic properties of tio2 nanoparticles...what to do now
it is a goid question and i want to know the answer as well
Maciej
characteristics of micro business
Abigail
for teaching engĺish at school how nano technology help us
Anassong
Do somebody tell me a best nano engineering book for beginners?
there is no specific books for beginners but there is book called principle of nanotechnology
NANO
what is fullerene does it is used to make bukky balls
are you nano engineer ?
s.
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
Mostly, they use nano carbon for electronics and for materials to be strengthened.
Virgil
is Bucky paper clear?
CYNTHIA
carbon nanotubes has various application in fuel cells membrane, current research on cancer drug,and in electronics MEMS and NEMS etc
NANO
so some one know about replacing silicon atom with phosphorous in semiconductors device?
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?
s.
how to fabricate graphene ink ?
for screen printed electrodes ?
SUYASH
What is lattice structure?
of graphene you mean?
Ebrahim
or in general
Ebrahim
in general
s.
Graphene has a hexagonal structure
tahir
On having this app for quite a bit time, Haven't realised there's a chat room in it.
Cied
what is biological synthesis of nanoparticles
how did you get the value of 2000N.What calculations are needed to arrive at it
Privacy Information Security Software Version 1.1a
Good
Got questions? Join the online conversation and get instant answers!  By By   By   By 