I'm analyzing a wind turbine blade dynamic analysis to get the mode shapes and natural frequencies,
so after I inserted all required airfoil cross section data (0m to 5m, having 10 cross-section) into PreVABS and run it at VABS,
now I need to put VABS output files all together as a (.dat) file to be run in GEBT.
Any example or tutorial of such case would be very much appreciated, noting that each cross-section is different from the other, regarding cord length pitch axis, etc...
PS: as we all know VABS has the following data regarding the airfoil cross section data:The 6X6 Mass Matrix, The Mass Center of the Cross Section, The 6X6 Mass Matrix at the Mass Center, The Mass Properties with respect to Principal Inertial Axes, Mass Per Unit Span, Mass Moments of Intertia about x1, x2 and x3 axises, The Principal Inertial Axes Rotated from User Coordinate System, Degree about the positive direction of x1 axis, The mass-weighted radius of gyration, The Geometric Center of the Cross Section Xg2 and Xg3, Classical Stiffness Matrix (1-extension; 2-twist; 3,4-bending), Classical Flexibility Matrix (1-extension; 2-twist; 3,4-bending), The Neutral Axes (or Tension Center) of the Cross Section Xt2 and Xt3, Timoshenko Stiffness Matrix (1-extension; 2,3-shear, 4-twist; 5,6-bending), Timoshenko Flexibility Matrix (1-extension; 2,3-shear, 4-twist; 5,6-bending), and The Generalized Shear Center of the Cross Section in the User Coordinate System Xs2 and Xs3.
(Aside,this work is a part of my master of science in aerospace engineering thesis, so it has a high importance for me)
Dear sirs,
I have managed such a file to run a modal analysis for my wind turbine blade having 10 different crossections but still Im having an error I dont know the reason,
can anybody tell please, very urgently
However, I think key point 2 is not a boundary point; change this key point
number to 10 (which I think this is the point at the tip) and try it again.
Chadi, Please use attachment instead on the board. Our team usually will not be responsive to email questions as we can Q&A to be shared with all other users.
Attached please find the corrected input file for this analysis. There are two mistakes in the original file:
The boundary condition should be applied to boundary points; for this case, the boundary points are 1 and 10.
The nframe tag in your case should be "0". Everything is defined in the global coordinates, thus, no local frame is needed.
Moreover, you really don't need so many elements for this case. In your original input file, you meshed each segment by 100 elements ( and there are 9 segments in total).
The file u fixed run now perfectly at GEBT now, but I have two questions:
1- why the eigenvalues shows two values one is positive while the other is negative? and why the first 6 modes are not around zero value Hz since I run the dynamic analysis flag 3 for freefree mode? is it correct that their values are in Hertz although we input the Rpm in rd/s ?
2- when I run the PostGEBT file for that it shows an error in line 33 which I couldnt fix to see the graphs of the mode shapes? (please see attachment) (while no such error for other examples = no installation mistake)
1 The eigensolver in GEBT returns all the "complex eigen pairs". PostGEBT will generate a file named "NaturalFrequency.pgt" in which only the natural frequencies are documented.
2 I run the PostGEBT successfully: the natural frequencies are documented in the ".pgt" file and all the corresponding eigenvectors are plotted in ".png" files.
1-why each member has its own mode for Disp and Force sepeerately but not all the beam as a whole ?because my professor asked me to show all the beam mode shapes related to each natural frequency.
2- why the limits of the RPM in rd/s are, 0.000001 and 83 for my case?
The PostGEBT will plot the eigenmodes for each beam at this moment. But you can output the quantities of your interested only by changing the "IO.f90" subroutine in your GEBT download package.
I don't quite understand what does "the limits" mean in your description.
The PostGEBT is plotting now the mode shapes (eigenmodes) for each member (of the same beam) separately as we saw,
- What quantities did you mean here sir? (if number of frequencies and related modeshapes of the output, I know where to edit in the main GEBT input)
- for (IO.f90) on which line of its subroutines we are able to fix such issue?
(why when Im doing the free-free modal analysis I do not see around 0 Hz values for the first 6 natural frequencies? I understood that only flexible mode shapes shows up in such case.)
my request is mode shapes for the whole beam all together.
chadi amer @ on — Edited @ on
Dear all,
I'm analyzing a wind turbine blade dynamic analysis to get the mode shapes and natural frequencies,
so after I inserted all required airfoil cross section data (0m to 5m, having 10 cross-section) into PreVABS and run it at VABS,
now I need to put VABS output files all together as a (.dat) file to be run in GEBT.
Any example or tutorial of such case would be very much appreciated, noting that each cross-section is different from the other, regarding cord length pitch axis, etc...
PS: as we all know VABS has the following data regarding the airfoil cross section data:The 6X6 Mass Matrix, The Mass Center of the Cross Section, The 6X6 Mass Matrix at the Mass Center, The Mass Properties with respect to Principal Inertial Axes, Mass Per Unit Span, Mass Moments of Intertia about x1, x2 and x3 axises, The Principal Inertial Axes Rotated from User Coordinate System, Degree about the positive direction of x1 axis, The mass-weighted radius of gyration, The Geometric Center of the Cross Section Xg2 and Xg3, Classical Stiffness Matrix (1-extension; 2-twist; 3,4-bending), Classical Flexibility Matrix (1-extension; 2-twist; 3,4-bending), The Neutral Axes (or Tension Center) of the Cross Section Xt2 and Xt3, Timoshenko Stiffness Matrix (1-extension; 2,3-shear, 4-twist; 5,6-bending), Timoshenko Flexibility Matrix (1-extension; 2,3-shear, 4-twist; 5,6-bending), and The Generalized Shear Center of the Cross Section in the User Coordinate System Xs2 and Xs3.
(Aside,this work is a part of my master of science in aerospace engineering thesis, so it has a high importance for me)
Kindest Regards,
B.S Aeronautical Engineer
M.Sc Aerospace Engineering Student
Su Tian @ on
Hi,
There are some supporting documents under the GEBT Tool webpage, including a GEBT manual and several examples with sample input files.
I am working with VABS and GEBT right now as well. So if you have further questions, we can discuss later.
Best regards,
Su Tian
chadi amer @ on
ok Mr. Tian
I have read the manual actually but I still find it hard to insert my data corresponding
would u please send me an email so we can discuss about it
mine is chadi_amer@hotmail.com
there I can show u my files I am working on or u may send me some examples simulating a similar case
thanks a lot
chadi amer @ on
Dear sirs, I have managed such a file to run a modal analysis for my wind turbine blade having 10 different crossections but still Im having an error I dont know the reason, can anybody tell please, very urgently
0.00E+005.54E-037.45E-043.56E-010.00E+000.00E+00 -5.54E-030.00E+000.00E+000.00E+001.76E-01-8.23E-04 -7.45E-040.00E+000.00E+000.00E+00-8.23E-041.80E-01
3.99E-096.07E-13-2.09E-131.24E-11-1.58E-081.41E-09 6.07E-132.59E-08-1.41E-092.62E-09-3.96E-11-1.31E-11 -2.09E-13-1.41E-093.43E-08-8.01E-09-1.85E-11-7.53E-12 -1.58E-08-3.96E-11-1.85E-11-2.05E-118.58E-072.40E-07
-1.0405803491E-13 -1.9771112026E-08 9.3979911387E-08 -1.2698921987E-08 -2.8328015753E-11 -1.0123873627E-11 -4.5694757922E-08 -4.4892344385E-11 -2.8328015753E-11 -1.1351405996E-10 2.0322494782E-06 6.2063368541E-07 -8.2146533760E-10 -1.2262430452E-11 -1.0123873627E-11 5.7687841030E-11 6.2063368541E-07 3.4401597992E-07
-1.9045753866E+00 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00 -2.8330114837E-01 7.8456325172E-01
-9.62E-08-1.76E-10-3.30E-11-3.03E-105.23E-061.12E-06 -1.12E-08-3.91E-11-2.74E-11-1.73E-101.12E-066.86E-07
-5.02E-010.00E+000.00E+000.00E+00-7.77E-023.22E-01
-1.14E-12-2.65E-082.77E-07-5.89E-081.87E-122.74E-12 -1.71E-07-3.48E-101.87E-12-1.08E-091.07E-051.26E-06
-6.15E-010.00E+000.00E+000.00E+00-3.73E-022.32E-01
-7.60E-12-3.48E-084.51E-07-4.34E-07-4.66E-12-1.20E-11 -3.67E-07-1.46E-09-4.66E-12-4.19E-094.12E-052.78E-06
-1.78E-12-6.67E-087.10E-07-1.44E-069.74E-10-7.04E-11 -1.07E-06-3.00E-099.74E-105.98E-091.16E-043.89E-06
-1.23E-010.00E+000.00E+000.00E+00-1.96E-032.71E-02
-3.92E-11-9.11E-088.90E-067.42E-062.42E-09-2.75E-09 -1.91E-06-7.78E-092.42E-09-2.82E-082.54E-044.33E-06
-7.87E-020.00E+000.00E+000.00E+00-8.02E-041.36E-02
-8.52E-07-4.83E-061.52E-04-9.27E-049.61E-05-2.46E-05 -1.95E-04-1.05E-059.61E-05-1.44E-031.00E-01-1.32E-03-4.73E-030.00E+000.00E+000.00E+00-7.29E-061.42E-04
The Error Message was as such at the end of the .txt.ech file
============================================
Failure of MA28AD with iflag=-1
chadi amer @ on
and some times if gives me
frame number unmatch
what does the frames resembles here?
any body knows ?
Thanks!
Qi Wang @ on — Edited @ on
Hi Chadi,
My first guess on the this bug is that your boundary conditions are not correct (ncond_pt). In your input file, it reads:
1 #point conditions 1 2 3 4 5 6 0 0 0 0 0 0 # corresponding value 0 0 0 0 0 0 # corresponding time functions 0 0 0 0 0 0 # indicate whether it is a follower condition
2 7 8 9 10 11 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
However, I think key point 2 is not a boundary point; change this key point number to 10 (which I think this is the point at the tip) and try it again.
Qi
Emailed forum response from aaron.qi.wang@gmail.com
chadi amer @ on
thanks for your opinion sir
but the error is as such
I/O error: read number of member frames
I dont really understand what does frames stands for
Qi Wang @ on — Edited @ on
Can you upload your input file as attachment please?
The data in the post is in poor format that I can’t help debug.
Qi
Emailed forum response from aaron.qi.wang@gmail.com
chadi amer @ on
can i have your email to send u my file directly and u may check it more carefully ?
mine is
chadi_amer@hotmail.com
Wenbin Yu @ on
Chadi, Please use attachment instead on the board. Our team usually will not be responsive to email questions as we can Q&A to be shared with all other users.
chadi amer @ on
ok sir, I didnt realize the icon of attachment at first. thanks
Qi Wang @ on
Hi Chadi,
Attached please find the corrected input file for this analysis. There are two mistakes in the original file:
Moreover, you really don't need so many elements for this case. In your original input file, you meshed each segment by 100 elements ( and there are 9 segments in total).
Qi
temp.dat
9 KBClick to download
chadi amer @ on
Thanks a lot Dr. Qi Wang,
The file u fixed run now perfectly at GEBT now, but I have two questions:
1- why the eigenvalues shows two values one is positive while the other is negative? and why the first 6 modes are not around zero value Hz since I run the dynamic analysis flag 3 for freefree mode? is it correct that their values are in Hertz although we input the Rpm in rd/s ?
2- when I run the PostGEBT file for that it shows an error in line 33 which I couldnt fix to see the graphs of the mode shapes? (please see attachment) (while no such error for other examples = no installation mistake)
Regards
Qi Wang @ on
Hi Chadi,
1 The eigensolver in GEBT returns all the "complex eigen pairs". PostGEBT will generate a file named "NaturalFrequency.pgt" in which only the natural frequencies are documented.
2 I run the PostGEBT successfully: the natural frequencies are documented in the ".pgt" file and all the corresponding eigenvectors are plotted in ".png" files.
Qi
chadi amer @ on
:) everything works perfectly now, much appreciation for all your time and patience Dr. Qi
and I hope that my thesis namely ''HIGH FIDELITY FEM FOR WIND TURBINE BLADE WITH MODAL TESTING VERIFICATION''
presents something new to the natural and applied science when its published this February 2015
Best Regards,
Chadi Amer
chadi amer @ on
two small question to add here,
1-why each member has its own mode for Disp and Force sepeerately but not all the beam as a whole ?because my professor asked me to show all the beam mode shapes related to each natural frequency.
2- why the limits of the RPM in rd/s are, 0.000001 and 83 for my case?
Qi Wang @ on
Hi Chadi,
The PostGEBT will plot the eigenmodes for each beam at this moment. But you can output the quantities of your interested only by changing the "IO.f90" subroutine in your GEBT download package.
I don't quite understand what does "the limits" mean in your description.
Qi
chadi amer @ on
Dear Sir,
The PostGEBT is plotting now the mode shapes (eigenmodes) for each member (of the same beam) separately as we saw,
- What quantities did you mean here sir? (if number of frequencies and related modeshapes of the output, I know where to edit in the main GEBT input)
- for (IO.f90) on which line of its subroutines we are able to fix such issue?
(why when Im doing the free-free modal analysis I do not see around 0 Hz values for the first 6 natural frequencies? I understood that only flexible mode shapes shows up in such case.)
my request is mode shapes for the whole beam all together.
IO.f90
23 KBClick to download