[b]Hi everyone,
I hope to do some simulations on a 2MW wind turbine by using FAST.
As the first step, I changed some parameters of the NREL WP-1.5MW (test.fst) to that of my 2MW wind turbine, as follows,
The primary input file:[/b]
[code]-------------------------------------------------------------------------------------- FAST INPUT FILE --------------------------------------------------------
FAST certification Test #13: WindPACT 1.5 MW Baseline with many DOFs with VS and VP and FF turbulence.
Model properties from “InputData1.5A08V07adm.xls” (from C. Hansen) with bugs removed. Compatible with FAST v7.02.00.
---------------------- SIMULATION CONTROL --------------------------------------
False Echo - Echo input data to “echo.out” (flag)
1 ADAMSPrep - ADAMS preprocessor mode {1: Run FAST, 2: use FAST as a preprocessor to create an ADAMS model, 3: do both} (switch)
1 AnalMode - Analysis mode {1: Run a time-marching simulation, 2: create a periodic linearized model} (switch)
3 NumBl - Number of blades (-)
700.0 TMax - Total run time (s)
0.005 DT - Integration time step (s)
---------------------- ENVIRONMENTAL CONDITIONS --------------------------------
9.80665 Gravity - Gravitational acceleration (m/s^2)
---------------------- FEATURE FLAGS -------------------------------------------
False FlapDOF1 - First flapwise blade mode DOF (flag)
False FlapDOF2 - Second flapwise blade mode DOF (flag)
False EdgeDOF - First edgewise blade mode DOF (flag)
False TeetDOF - Rotor-teeter DOF (flag) [unused for 3 blades]
True DrTrDOF - Drivetrain rotational-flexibility DOF (flag)
True GenDOF - Generator DOF (flag)
False YawDOF - Yaw DOF (flag)
False TwFADOF1 - First fore-aft tower bending-mode DOF (flag)
False TwFADOF2 - Second fore-aft tower bending-mode DOF (flag)
False TwSSDOF1 - First side-to-side tower bending-mode DOF (flag)
False TwSSDOF2 - Second side-to-side tower bending-mode DOF (flag)
True CompAero - Compute aerodynamic forces (flag)
False CompNoise - Compute aerodynamic noise (flag)
---------------------- TURBINE CONFIGURATION -----------------------------------
46.451 TipRad - The distance from the rotor apex to the blade tip (meters)!!!
1.17 HubRad - The distance from the rotor apex to the blade root (meters)!!!
1 PSpnElN - Number of the innermost blade element which is still part of the pitchable portion of the blade for partial-span pitch control [1 to BldNodes] [CURRENTLY IGNORED] (-)
0.0 UndSling - Undersling length [distance from teeter pin to the rotor apex] (meters) [unused for 3 blades]
0.0 HubCM - Distance from rotor apex to hub mass [positive downwind] (meters)
-3.492 OverHang - Distance from yaw axis to rotor apex [3 blades] (meters)
-1 NacCMxn - Downwind distance from the tower-top to the nacelle CM (meters)
0.0 NacCMyn - Lateral distance from the tower-top to the nacelle CM (meters)
1.430 NacCMzn - Vertical distance from the tower-top to the nacelle CM (meters)
77.99 TowerHt - Height of tower above ground level [onshore] or MSL [offshore] (meters)
2.01 Twr2Shft - Vertical distance from the tower-top to the rotor shaft (meters)
0.0 TwrRBHt - Tower rigid base height (meters)
-5.0 ShftTilt - Rotor shaft tilt angle (degrees)
0.0 Delta3 - Delta-3 angle for teetering rotors (degrees) [unused for 3 blades]
3.5 PreCone(1) - Blade 1 cone angle (degrees)
3.5 PreCone(2) - Blade 2 cone angle (degrees)
3.5 PreCone(3) - Blade 3 cone angle (degrees) [unused for 2 blades]
0.0 AzimB1Up - Azimuth value to use for I/O when blade 1 points up (degrees)
---------------------- MASS AND INERTIA ----------------------------------------
0.0 YawBrMass - Yaw bearing mass (kg)
87E3 NacMass - Nacelle mass (kg)
20E3 HubMass - Hub mass (kg)
0.0 TipMass(1) - Tip-brake mass, blade 1 (kg)
0.0 TipMass(2) - Tip-brake mass, blade 2 (kg)
0.0 TipMass(3) - Tip-brake mass, blade 3 (kg)
100E3 NacYIner - Nacelle inertia about yaw axis (kg m^2)
89 GenIner - Generator inertia about HSS (kg m^2)
10E3 HubIner - Hub inertia about rotor axis [3 blades](kg m^2)
[/code]
The Aerodyn input file:
1.5 MW baseline aerodynamic parameters for FAST certification test #13.
SI SysUnits - System of units for used for input and output [must be SI for FAST] (unquoted string)
BEDDOES StallMod - Dynamic stall included [BEDDOES or STEADY] (unquoted string)
NO_CM UseCm - Use aerodynamic pitching moment model? [USE_CM or NO_CM] (unquoted string)
DYNIN InfModel - Inflow model [DYNIN or EQUIL] (unquoted string)
SWIRL IndModel - Induction-factor model [NONE or WAKE or SWIRL] (unquoted string)
0.005 AToler - Induction-factor tolerance (convergence criteria) (-)
PRANDtl TLModel - Tip-loss model (EQUIL only) [PRANDtl, GTECH, or NONE] (unquoted string)
PRANDtl HLModel - Hub-loss model (EQUIL only) [PRANdtl or NONE] (unquoted string)
"Wind/WP_Baseline/Shr12_02.wnd" WindFile - Name of file containing wind data (quoted string)
80 HH - Wind reference (hub) height [TowerHt+Twr2Shft+OverHang*SIN(ShftTilt)] (m)
0.0 TwrShad - Tower-shadow velocity deficit (-)
9999.9 ShadHWid - Tower-shadow half width (m)
9999.9 T_Shad_Refpt - Tower-shadow reference point (m)
1.225 AirDens - Air density (kg/m^3)
1.4639E-5 KinVisc - Kinematic air viscosity (m^2/sec)
0.005 DTAero - Time interval for aerodynamic calculations (sec)
4 NumFoil - Number of airfoil files (-)
"AeroData/WP_Baseline/cylinder.dat" FoilNm - Names of the airfoil files [NumFoil lines] (quoted strings)
"AeroData/WP_Baseline/s818_2703.dat"
"AeroData/WP_Baseline/s825_2103.dat"
"AeroData/WP_Baseline/s826_1603.dat"
26 BldNodes - Number of blade nodes used for analysis (-)!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
RNodes AeroTwst DRNodes Chord NFoil PrnElm
2.0140 10.8200 1.6880 2.2520 1
3.7015 13.0300 1.6870 2.5820 1
5.3890 14.8200 1.6880 2.9590 2
7.0765 15.1300 1.6870 3.3360 2
8.7640 14.4200 1.6880 3.5700 2
10.4515 11.7300 1.6870 3.5120 2
12.1390 9.6900 1.6880 3.2960 2
13.8265 8.0800 1.6870 3.0590 3
15.5140 6.7800 1.6880 2.8140 3
17.2015 5.7200 1.6870 2.5800 3
18.8890 4.8500 1.6880 2.3720 3
20.5765 4.1100 1.6870 2.1920 3
22.2640 3.4800 1.6880 2.0350 3
23.9515 2.9300 1.6870 1.8990 3
25.6390 2.4500 1.6880 1.7800 3
27.3265 2.0300 1.6870 1.6750 3
29.0140 1.6000 1.6880 1.5810 3
30.7015 1.2200 1.6870 1.4950 3
32.3890 0.8600 1.6880 1.4130 3
34.0765 0.5400 1.6870 1.3320 3
35.7640 0.2400 1.6880 1.2500 4
37.4515 -0.0400 1.6870 1.1660 4
39.1390 -0.3000 1.6880 1.0810 4
40.8265 -0.5400 1.6870 0.9960 4
42.5140 -0.7400 1.6880 0.9110 4
44.2015 -0.8600 1.6870 0.8260 4
45.7480 5.1400 1.4060 0.0220 4
I disabled all degrees-of-freedoms expect GenDOF and DrTrDOF. Then I run FAST. The following error appears:
Running AeroDyn (v13.00.02a-bjj, 20-Feb-2013).
Heading of the AeroDyn input file: 1.5 MW baseline aerodynamic parameters for FAST certification
test #13.
TipRadius must be equal to HubRadius + SUM( DR(:) ).
Could not initialize AeroDyn.
Aborting FAST.
[b]I checked the parameters. But found that my parameters absolutely meet the requirement of TipRadius=HubRadius + SUM( DR(
):
HubRadius + SUM( DR( )=1.17+sum(1.688 1.687 1.688 1.687 1.688 1.687 1.688 1.687 1.688 1.687 1.688 1.687 1.688 1.687 1.688 1.687 1.688 1.687 1.688 1.687 1.688 1.687 1.688 1.687 1.688 1.687 1.406)=46.451=TipRadius
It is appreciated if someone could give me some advices about how to check where the problem is? Thanks in advance[/b]
Best Regards,
Chao