WT_perf使用手册

WT_Perf主要用来预测机组功率

 以下是输入文件

-----  WT_Perf Input File  -----------------------------------------------------
WT_Perf Test01 input file.  UAE Phase 3 turbine (Non-dimen, English, Space, PROP-PC).
Compatible with WT_Perf v3.00f
-----  Input Configuration  ----------------------------------------------------
False                Echo:                      Echo input parameters to "<rootname>.ech"?
False                DimenInp:                  Turbine parameters are dimensional?
False                Metric:                    Turbine parameters are Metric (MKS vs FPS)?
-----  Model Configuration  ----------------------------------------------------
16                   NumSect:                   Number of circumferential sectors.                 ##叶素个数,至少4个,没有上限##
5000                 MaxIter:                   Max number of iterations for induction factor. ##诱导系数最大迭代计算次数,超出迭代次数,输出功率设置为-999.99##
1.0e-6               ATol:                      Error tolerance for induction iteration.              ##诱导系数计算最大容差##
1.0e-6               SWTol:                     Error tolerance for skewed-wake iteration.    ##旋转尾流计算容差##
-----  Algorithm Configuration  ------------------------------------------------
True                 TipLoss:                   Use the Prandtl tip-loss model?
False                HubLoss:                   Use the Prandtl hub-loss model?
True                 Swirl:                     Include Swirl effects?
True                 SkewWake:                  Apply skewed-wake correction?
True                 AdvBrake:                  Use the advanced brake-state model?
True                 IndProp:                   Use PROP-PC instead of PROPX induction algorithm?
False                AIDrag:                    Use the drag term in the axial induction calculation?
False                TIDrag:                    Use the drag term in the tangential induction calculation?
-----  Turbine Data  -----------------------------------------------------------
3                    NumBlade:                  Number of blades.
16.5                 RotorRad:                  Rotor radius [length].
0.2                  HubRad:                    Hub radius [length or div by radius].
3.5                  PreCone:                   Precone angle, positive downwind [deg].
0.0                  Tilt:                      Shaft tilt [deg].
10.0                 Yaw:                       Yaw error [deg].
3.3333               HubHt:                     Hub height [length or div by radius].
16                   NumSeg:                    Number of blade segments (entire rotor radius).
   RElm      Twist      Chord  AFfile  PrntElem     ###弦长扭角分布###
  0.225      0.000     0.0911     1      False
  0.275      0.000     0.0911     1      False
  0.325      0.000     0.0911     1      False
  0.375      0.000     0.0911     1      False
  0.425      0.000     0.0911     1      False
  0.475      0.000     0.0911     1      False
  0.525      0.000     0.0911     1      False
  0.575      0.000     0.0911     1      False
  0.625      0.000     0.0911     1      False
  0.675      0.000     0.0911     1      False
  0.725      0.000     0.0911     1      False
  0.775      0.000     0.0911     1      False
  0.825      0.000     0.0911     1      False
  0.875      0.000     0.0911     1      False
  0.925      0.000     0.0911     1      False
  0.975      0.000     0.0911     1      t False
-----  Aerodynamic Data  -------------------------------------------------------
0.0019749                                Rho:                 Air density [mass/volume].
0.0001625                                KinVisc:             Kinematic air viscosity
0.143                                    ShearExp:            Wind shear exponent (1/7 law = 0.143).
False                                    UseCm:               Are Cm data included in the airfoil tables?
1                                        NumAF:               Number of airfoil files.
"airfoils/unsteadyaeroexp/s809_cln.dat"  AF_File:             List of NumAF airfoil files.
-----  I/O Settings  -----------------------------------------------------------
True                 TabDel:                    Make output tab-delimited (fixed-width otherwise).
True                 KFact:                     Output dimensional parameters in K (e.g., kN instead on N)
True                 WriteBED:                  Write out blade element data to "<rootname>.bed"?
False                InputTSR:                  Input speeds as TSRs?
"fps"                SpdUnits:                  Wind-speed units (mps, fps, mph).
-----  Combined-Case Analysis  -------------------------------------------------
0                    NumCases:                  Number of cases to run.  Enter zero for parametric analysis.
WS or TSR   RotSpd   Pitch                      Remove following block of lines if NumCases is zero.
-----  Parametric Analysis (Ignored if NumCases > 0 )  -------------------------
3                    ParRow:                    Row parameter    (1-rpm, 2-pitch, 3-tsr/speed).
2                    ParCol:                    Column parameter (1-rpm, 2-pitch, 3-tsr/speed).
1                    ParTab:                    Table parameter  (1-rpm, 2-pitch, 3-tsr/speed).
True                 OutPwr:                    Request output of rotor power?
True                 OutCp:                     Request output of Cp?
True                 OutTrq:                    Request output of shaft torque?
True                 OutFlp:                    Request output of flap bending moment?
True                 OutThr:                    Request output of rotor thrust?
3, 4, 1              PitSt, PitEnd, PitDel:     First, last, delta blade pitch (deg).
72, 73, 1            OmgSt, OmgEnd, OmgDel:     First, last, delta rotor speed (rpm).
15, 75, 1            SpdSt, SpdEnd, SpdDel:     First, last, delta speeds.

接下来就要定义翼型的升阻力系数分布

AeroDyn airfoil file.  Compatible with AeroDyn v13.0.
S809 Airfoil, OSU data at Re=.75 Million, Clean roughness
NREL/TP-442-7817 Appendix B, Viterna used aspect ratio=11
  1            Number of airfoil tables in this file
  9.0e9        Table ID parameter (Reynolds number in milllions).  For efficiency, make very large if only one table.
 15.3          Stall angle (deg)
 -0.38         Zero lift angle of attack (deg)
  7.12499      Cn slope for zero lift (dimensionless)
  1.9408       Cn at stall value for positive angle of attack
 -0.8000       Cn at stall value for negative angle of attack
  2.0000       Angle of attack for minimum CD (deg)
  0.0116       Minimum CD value
  -180.00      .000    .1748    .0000
  -170.00      .230    .2116    .4000
  -160.00      .460    .3172    .1018
  -150.00      .494    .4784    .1333
  -140.00      .510    .6743    .1727
  -130.00      .486    .8799    .2132
  -120.00      .415   1.0684    .2498
  -110.00      .302   1.2148    .2779
  -100.00      .159   1.2989    .2933
   -90.00      .000   1.3080    .2936
   -80.00     -.159   1.2989    .2933
   -70.00     -.302   1.2148    .2779
   -60.00     -.415   1.0684    .2498
   -50.00     -.486    .8799    .2132
   -40.00     -.510    .6743    .1727
   -30.00     -.494    .4784    .1333
   -20.10     -.560    .3027    .0612
   -18.10     -.670    .3069    .0904
   -16.10     -.790    .1928    .0293
   -14.20     -.840    .0898   -.0090
   -12.20     -.700    .0553   -.0045
   -10.10     -.630    .0390   -.0044
    -8.20     -.560    .0233   -.0051
    -6.10     -.640    .0131    .0018
    -4.10     -.420    .0134   -.0216
    -2.10     -.210    .0119   -.0282
      .10      .050    .0122   -.0346
     2.00      .300    .0116   -.0405
     4.10      .540    .0144   -.0455
     6.20      .790    .0146   -.0507
     8.10      .900    .0162   -.0404
    10.20      .930    .0274   -.0321
    11.30      .920    .0303   -.0281
    12.10      .950    .0369   -.0284
    13.20      .990    .0509   -.0322
    14.20     1.010    .0648   -.0361
    15.30     1.020    .0776   -.0363
    16.30     1.000    .0917   -.0393
    17.10      .940    .0994   -.0398
    18.10      .850    .2306   -.0983
    19.10      .700    .3142   -.1242
    20.10      .660    .3186   -.1155
    30.00      .705    .4784   -.2459
    40.00      .729    .6743   -.2813
    50.00      .694    .8799   -.3134
    60.00      .593   1.0684   -.3388
    70.00      .432   1.2148   -.3557
    80.00      .227   1.2989   -.3630
    90.00      .000   1.3080   -.3604
   100.00     -.159   1.2989   -.3600
   110.00     -.302   1.2148   -.3446
   120.00     -.415   1.0684   -.3166
   130.00     -.486    .8799   -.2800
   140.00     -.510    .6743   -.2394
   150.00     -.494    .4784   -.2001
   160.00     -.460    .3172   -.1685
   170.00     -.230    .2116   -.5000
   180.00      .000    .1748    .0000
EOT

纠正代码:trainsets = pd.read_csv('/Users/zhangxinyu/Desktop/trainsets82.csv') testsets = pd.read_csv('/Users/zhangxinyu/Desktop/testsets82.csv') y_train_forced_turnover_nolimited = trainsets['m3_forced_turnover_nolimited'] X_train = trainsets.drop(['m3_P_perf_ind_all_1','m3_P_perf_ind_all_2','m3_P_perf_ind_all_3','m3_P_perf_ind_allind_1',\ 'm3_P_perf_ind_allind_2','m3_P_perf_ind_allind_3','m3_P_perf_ind_year_1','m3_P_perf_ind_year_2',\ 'm3_P_perf_ind_year_3','m3_forced_turnover_nolimited','m3_forced_turnover_3mon',\ 'm3_forced_turnover_6mon','m3_forced_turnover_1year','m3_forced_turnover_3year',\ 'm3_forced_turnover_5year','m3_forced_turnover_10year',\ 'CEOid','CEO_turnover_N','year','Firmid','appo_year'],axis=1) y_test_forced_turnover_nolimited = testsets['m3_forced_turnover_nolimited'] X_test = testsets.drop(['m3_P_perf_ind_all_1','m3_P_perf_ind_all_2','m3_P_perf_ind_all_3','m3_P_perf_ind_allind_1',\ 'm3_P_perf_ind_allind_2','m3_P_perf_ind_allind_3','m3_P_perf_ind_year_1','m3_P_perf_ind_year_2',\ 'm3_P_perf_ind_year_3','m3_forced_turnover_nolimited','m3_forced_turnover_3mon',\ 'm3_forced_turnover_6mon','m3_forced_turnover_1year','m3_forced_turnover_3year',\ 'm3_forced_turnover_5year','m3_forced_turnover_10year',\ 'CEOid','CEO_turnover_N','year','Firmid','appo_year'],axis=1) from sklearn.ensemble import RandomForestClassifier rfc = RandomForestClassifier(n_estimators=100, max_depth=10, random_state=42) rfc.fit(X_train, y_train_forced_turnover_nolimited) y_pred = rfc.predict_proba(X_test) # 计算AUC值 auc = roc_auc_score(y_test_forced_turnover_nolimited, y_pred) # 输出AUC值 print('测试集AUC值为:', auc)
06-03
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值