Dear Jason
when i write
RotTripletIndicesStates = 3
RotTripletIndicesStates =
3
Trial>> ndof =15
ndof =
15
Trial>> mbc3
Running mbc3 (v1.00.00a-gbir, 15-Oct-2008)
**ERROR: the number of colmn vectors in RotTripletIndicesStates must equal 3, the num of blades
Dear Rana,
You should not have to write “RotTripletIndicesStates = 3”; instead, RotTripletIndicesStates should be a matrix written by GetMats_f8.m.
Best regards,
Dear Jason
I very very thankful for your patience, I apologize if i bothering you.
Finally it works.
my regrads
Hi all,
How can I get the averaged state space matrixes for the linearized periodic model? I got the periodic state matrices from the FAST linearization output file. But I don’t know how to use it based on it is periodic.
Thank you!
Yanhua
Dear Yanhua,
For a 3-bladed rotor with states in the rotating frame, I would normally expect that you’d linearize a sizeable number of times (e.g. 36 times, linearizing once every 10 degrees in azimuth angle, to get the full 360-degree periodic model). Then, you should apply the MBC tool (nwtc.nrel.gov/MBC) to transform the rotating states into the fixed frame and azimuth-average the transformed linear model to get linear time-invariant (LTI) matrices. If you don’t have states in the rotating frame (i.e. by disabling the blade degrees of freedom), then you wouldn’t need to apply MBC and you could azimuth-average the period matrices directly. See the documentation supplied with MBC for more information.
Best regards,
Hi all,
How can I get the averaged state space matrixes for the linearized periodic model? I got the periodic state matrices from the FAST v7 linearization output file. I guess that Eigenanalysis.m can get the azimuth-averaged matrices. But I can’t find the file Eigenanalysis.m in the CertTest file. Can the linearized periodic model from FAST V7 use the MBC procedure for FAST V8?
Thank you!
Yanhua
Dear Yanhua,
Use the MBC tool (nwtc.nrel.gov/MBC), which works with linearization files from both FAST v7 and FAST v8.
Best regards,
Dear Jason,
Thank you so much for your nice answer. I tried to use MBC to get the LTI model from the periodic model. I am using the first-order system matrix. However, I found that only the A matrix to be azimuth-averaged. The matrix B, Bd, C, D, Dd are still time-varying. How can I get these corresponding azimuth-averaged matrix for B, Bd, C, D, Dd?
Thank you.
All the best.
Yanhua
Dear Yanhua,
I’m not sure why the MBC3 script does not automatically azimuth-average the other MBC-transformed matrices (B, Bd, C, D, Dd), but the calculation is easy to implement – simply copy the logic used to calculate MBC_AvgA.
Best regards,
Dear Jason,
Thank you so much. I find the reason why I only got the average-azimuthed A. It is because I run ccm.m which is only for A matrix. Actually, I got 36 linearized models because it’s based on the azimuth angle (10 degrees increase each time) and I attached the .lin file (changed into .txt). After the MBC3.mat, 36 linearized models are still periodic to azimuth angle. So what is the meaning of the MBC transformation? After the MBC3.mat, how to use these 36 linearized models for controller design?
All the best.
Yanhua
init5MW.lin.txt (214 KB)
Dear Yanhua,
I have not looked at your attachment, but I would not use the cce.m script. It is far better to apply MBC3 and then azimuth-average the MBC-transformed matrices. Please see the MBC3 documentation for information on what it is doing. As discussed in my post dated May 07, 2010 in the following forum topic: Fast Linearized Models - #4 by Jason.Jonkman, MBC does not eliminate the periodicity, it simply transforms the rotating states of a model into the nonrotating frame and further azimuth-averaging of the lightly-periodic MBC-transformed matrices (i.e. averaging the 36 matrices into 1 single matrix) has been shown to be quite adequate for controls design.
Best regards,
Dear Jason,
Thanks for your reply. I linearized the model at 36 azimuth angles (10 degrees increase) then obtained the first-order average-azimuthed linear model using MBC. I activated 4 DOFs (Generator DOF and First flapwise blade mode DOF). Now I am confusing about the state vector.
As I know, q is the vector of DOF displacements. So in my case q will be [xr, xbfl1, xbfl2, xbfl3], which means rotor displacement, the flapwise tip displacement of the i-th blade, respectively. Then the state vector for the first-order linearized model will be [xr, xbfl1, xbfl2, xbfl3,xr’,xbfl1’, xbfl2’, xbfl3’]. And I think xr’ is the rotor speed. After MBC transformation, the state vector changes into [xr, xbflo,xbflc, xbfls, xr’,xbflo’,xbflc’, xbfls’], where xbflo,xbflc, xbfls are the zero, cosine, sine components respectively.
I am not sure what the DOF displacement of Generator DOF is? Is it the rotor displacement which can be obtained from the integration of rotor speed. Is it right? If it is not, what it will be? Hope you can help me. Thank you so much.
All the best.
Yanhua
Dear Yanhua,
Your understanding is correct, but I would call xr the perturbation of the rotor-azimuth angle, not the integral of the rotor speed. Keep in mind that the states of the linear model are perturbations of the degrees of freedom (DOFs) about the operating point, not the DOFs themselves.
Best regards,
Dear Jason,
Thanks for your timely reply. I am trying to compare the simulation results between the nonlinear model and the obtained linearized model. But I am not sure how to compare. My linearization is done under wind speed 18m/s with 0.2 vertical shear, 14.75-degree pitch angle. According to you, after MBC transformation, the state vector changes into [azimuth-r, xbflo,xbflc, xbfls, azimuth-r’,xbflo’,xbflc’, xbfls’], where azimuth-r means the rotor-azimuth angle. How can I get the perturbations of the degrees of freedom (DOFs) about the operating point? For example, xbflo1=xbflo-xbflop, where xbflo1 means the perturbation of the flapwise tip displacement zero component, xbflo is the real time value from the sensor, xbflop is the value in this operating point. Can you give me some hints? I am really confusing how to evaluate the obtained linearized model. Thank you.
All the best.
Yanhua
Dear Yanhua,
It sounds like you understand how the perturbations are defined. I’m not what else I can explain. Perhaps searching the forum for other FAST linearization-related discussion would be of use to you?
Best regards,
Dear Jason,
I intend to linearise the NREL 5MW turbine model in region 3 of operation. About the linearization of the file I put the annex.I have used MAC software to calculate the matrix of the state space model and build the model in simulink.But I don’t know how to determine the value of input and wind disturbance. Can you help me?
Dear Zhe.Zhang,
The inputs (u) and disturbances (ud) are the perturbations of the controller input (collective blade-pitch angle) and wind disturbance (wind speed) about their operating point values. I would expect that you’d set the wind disturbance to be what you want (e.g. 1 m/s steps) and that you’d design/implement a controller to calculate the controller input derived from the system outputs.
Best regards,
Dear Jason,
Thanks for your timely reply.Are the initial values for these states all zero?
Dear Zhe.Zhang,
You can set the initial states in your Simulink model. See the MATLAB/Simulink documentation for more information.
Best regards,
Dear Jason,
I know that the dimensions of these state vectors are set according to the DOF in FAST, but how are the initial values of these states determined, they represent the deviation from the set point, so are they all 0?
Best regards,