Ehsan 506 Project Report
Ehsan 506 Project Report
} [
] {
}
[
] [
]
Taking the terms with X to the right side and solving for steady-state conditions where
where is the frequency of excitation:
[
]
The matrix
is called the dynamic stiffness matrix. Taking the inverse of this matrix to find X in terms of
:
The matrix
is the dynamic compliance matrix also known as the position FRF matrix. Knowing F and
X we can find this matrix. However, for the system at hand we know X but we do not know F and even
assuming the proportionality we can only find the torque applied on the rotor which cannot be directly
related to F without writing the equations of motion for the whole system. Therefore, we have to consider
the full system to be able to find F.
Figure 3. Free-body diagram for the motor and the first cart
5
2.1. The complete system
The free_body diagram of the rotor and pulley is shown in Figure 3.
The Newtons equation of motion for the rotary coordinate of the pulley is:
As the lower part of the figure shows the net force applied to mass 1 is:
Thus, the relation between the torque and the net force F assumed in the previous section is:
Substituting for
] {
} {
}
[
] {
} {
]
The FRF that we used to find by dividing the output voltage of the encoders by the voltage of excitation is in
fact proportional to the FRF obtained by
.
The output amplitude of oscillations is known and is the same for the simple 2-DoF system and the full
system (including the rotor.) To obtain the correct FRF for the 2-DoF system (without the rotor) however,
we have to divide this amplitude by the force F instead of
We can find
Substituting for
The above equation relates the FRF for the full system and system without the rotor at each frequency.
Figure 4 shows the FRF that is corrected for rotor inertia versus the FRF that is not corrected.
Figure 4. Comparison of calculating the FRF with and without considering the effect of rotor inertia
3.3 System Identification for the Excitation Unit
Next, we assume the moment of inertia of the rotor and its torque constant, as well as the value for the motor
drives command gain given by the factory are not reliable. With this assumption we lose less generality in our
analysis, because we know that each DC motor has an internal impedance and we do not have much
information about the circuit inside the drive. So we look at the drive and motor and pulley system as a black
box for which the voltage is the input and the net force on the first cart applied by the cables under it is the
output and we try to find the transfer function for this black box:
DAQ
AO0
Voltage
Net External
Force
Motor & Drive
System
The Effect of Rotor Inertia not considered
The Effect of Rotor Inertia considered
P
o
s
i
t
i
o
n
F
R
F
(
m
/
N
)
7
We do not have any extensometers or any other means to measure the tension in the cables so in order to
simplify the problem we disconnect the first cart from the springs it is connected to. Considering the friction
forces to be negligible with respect to the inertial forces the free-body diagram for this cart will be:
The equation of motion for this system in steady-state harmonic motion is:
So, for this simplified 1-DoF system without any springs; knowing the X, force can be easily found. After
running a test with this configuration we can have an estimate for the transfer function between the DAC
output voltage and the actual excitation force on the MSD system. However, still this estimation is not
reliable before 3 Hz.
Figure 5. Position frequency response with respect to input voltage
F
N
W
8
Figure 6. Excitation force frequency response with respect to input voltage
We tried to fit rational fraction polynomials to this response (which is also the transfer function for the
excitation unit) but the simplest function which could acceptably follow the trend was 10
th
order in the
numerator and 5
th
order in the denominator. However, we can see that the for frequencies higher than 3 Hz
the response is quite constant and is equal to 0.0141 (N/V) Therefore, to find the frequency response with
respect to the input force to the 2-DoF system we first need to multiply the voltage Fourier transform by this
transfer function and get the force. This method is also implemented in the MATLAB code.
4. Conclusion
Changes to the hardware (wiring, toggle switch) and the software that operate the MSD setup made this setup
safer, more accurate and more reliable. Previously, the FRF was acquired without considering either the effect
that impedance of the motor and drive impose on the input voltage or the effect of the rotor inertia on the
resultant excitation force on the MSD system.
5. References
[1] S200 High Performance Compact Brushless Servo Drives, Reference Manual, Kollmorgen
[2] AKM Series Motors Selection Guide, Daner Motion
[3] Theoritical and Experimental Modal Analysis, Maia and Silva, Research Studies Press