From 783211c142ce4412c7663cc4b0a62b2bbf7cf284 Mon Sep 17 00:00:00 2001 From: pibo Date: Mon, 24 May 2021 21:32:23 -0600 Subject: [PATCH 1/5] improve a few things to run olaf from weis --- .../IEA-15-240-RWT/IEA-15-240-RWT_BeamDyn.dat | 2 +- .../IEA-15-240-RWT_BeamDyn_blade_VABS.dat | 325 ------------------ examples/01_aeroelasticse/run_OLAF.py | 2 +- weis/aeroelasticse/FAST_reader.py | 21 +- weis/aeroelasticse/FAST_writer.py | 2 + weis/aeroelasticse/openmdao_openfast.py | 49 ++- 6 files changed, 43 insertions(+), 358 deletions(-) delete mode 100644 examples/01_aeroelasticse/OpenFAST_models/IEA-15-240-RWT/IEA-15-240-RWT/IEA-15-240-RWT_BeamDyn_blade_VABS.dat diff --git a/examples/01_aeroelasticse/OpenFAST_models/IEA-15-240-RWT/IEA-15-240-RWT/IEA-15-240-RWT_BeamDyn.dat b/examples/01_aeroelasticse/OpenFAST_models/IEA-15-240-RWT/IEA-15-240-RWT/IEA-15-240-RWT_BeamDyn.dat index 9efc0c850..801760d4d 100644 --- a/examples/01_aeroelasticse/OpenFAST_models/IEA-15-240-RWT/IEA-15-240-RWT/IEA-15-240-RWT_BeamDyn.dat +++ b/examples/01_aeroelasticse/OpenFAST_models/IEA-15-240-RWT/IEA-15-240-RWT/IEA-15-240-RWT_BeamDyn.dat @@ -75,7 +75,7 @@ True RotStates - Orient states in the rotating frame during linea ---------------------- MESH PARAMETER ------------------------------------------ 5 order_elem - Order of interpolation (basis) function (-) ---------------------- MATERIAL PARAMETER -------------------------------------- -"IEA 15MW Offshore Reference Turbine, with taped chord tip design_BeamDyn_Blade.dat" BldFile - Name of file containing properties for blade (quoted string) +"../IEA-15-240-RWT/IEA-15-240-RWT_BeamDyn_blade.dat" BldFile - Name of file containing properties for blade (quoted string) ---------------------- PITCH ACTUATOR PARAMETERS ------------------------------- False UsePitchAct - Whether a pitch actuator should be used (flag) 200 PitchJ - Pitch actuator inertia (kg-m^2) [used only when UsePitchAct is true] diff --git a/examples/01_aeroelasticse/OpenFAST_models/IEA-15-240-RWT/IEA-15-240-RWT/IEA-15-240-RWT_BeamDyn_blade_VABS.dat b/examples/01_aeroelasticse/OpenFAST_models/IEA-15-240-RWT/IEA-15-240-RWT/IEA-15-240-RWT_BeamDyn_blade_VABS.dat deleted file mode 100644 index 094d64635..000000000 --- a/examples/01_aeroelasticse/OpenFAST_models/IEA-15-240-RWT/IEA-15-240-RWT/IEA-15-240-RWT_BeamDyn_blade_VABS.dat +++ /dev/null @@ -1,325 +0,0 @@ - ------- BEAMDYN V1.00.* INDIVIDUAL BLADE INPUT FILE -------------------------- - IEA 15MW Offshore Reference Turbine, with taped chord tip design blade. Generated with Sonata/VABS (includes initial twist and curvature) by Roland Feil et al (A Cross-Sectional Aeroelastic Analysis and Structural Optimization Tool for Slender Composite Structures) - ---------------------- BLADE PARAMETERS -------------------------------------- -21 station_total - Number of blade input stations (-) - 1 damp_type - Damping type: 0: no damping; 1: damped - ---------------------- DAMPING COEFFICIENT------------------------------------ - mu1 mu2 mu3 mu4 mu5 mu6 - (-) (-) (-) (-) (-) (-) -1.0E-02 1.0E-02 1.0E-02 1.0E-02 1.0E-02 1.0E-02 - ---------------------- DISTRIBUTED PROPERTIES--------------------------------- - 0.000000 - 6.6908231106999998e+09 4.9875429303000001e+06 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 1.5897781766000000e+08 - 4.9875429303000001e+06 6.6810923405000000e+09 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 -1.3543091622000000e+07 - 0.0000000000000000e+00 0.0000000000000000e+00 4.6040677402000000e+10 -1.0884485653000000e+09 3.5931854910000002e+08 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 -1.0884485653000000e+09 1.4902591801000000e+11 -3.8536788920999996e+07 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 3.5931854910000002e+08 -3.8536788920999996e+07 1.4915276481000000e+11 0.0000000000000000e+00 - 1.5897781766000000e+08 -1.3543091622000000e+07 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 8.7446200981000000e+10 - - 3.1271965618999998e+03 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 7.3931936436000001e+01 - 0.0000000000000000e+00 3.1271965618999998e+03 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 1.2083518807000001e+01 - 0.0000000000000000e+00 0.0000000000000000e+00 3.1271965618999998e+03 -7.3931936436000001e+01 -1.2083518807000001e+01 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 -7.3931936436000001e+01 1.0163145202000000e+04 6.0929376303000005e-01 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 -1.2083518807000001e+01 6.0929376303000005e-01 1.0164748801000000e+04 0.0000000000000000e+00 - 7.3931936436000001e+01 1.2083518807000001e+01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 2.0327894003000001e+04 - - 0.050000 - 4.2273972884000001e+09 -1.9412711552999999e+07 -1.3749566443000000e+07 -1.6386787851000000e+07 -8.9952999974000006e+04 -4.3159278444000000e+08 - -1.9412711552999999e+07 5.1900302736000004e+09 5.8887223583000002e+05 -1.7151776432000000e+06 8.9412399243999999e+06 9.3299144861000001e+07 - -1.3749566443000000e+07 5.8887223583000002e+05 3.3748426993000000e+10 3.9718481230000000e+09 -1.8751484559999999e+08 9.1321652960000001e+06 - -1.6386787851000000e+07 -1.7151776432000000e+06 3.9718481230000000e+09 1.0261402304000000e+11 -2.2858830309000001e+09 2.5070575329000000e+07 - -8.9952999974000006e+04 8.9412399243999999e+06 -1.8751484559999999e+08 -2.2858830309000001e+09 9.1923974113000000e+10 6.0070288247000001e+07 - -4.3159278444000000e+08 9.3299144861000001e+07 9.1321652960000001e+06 2.5070575329000000e+07 6.0070288247000001e+07 5.4176806420000000e+10 - - 2.2041158952999999e+03 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 -2.7142512751999999e+02 - 0.0000000000000000e+00 2.2041158952999999e+03 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 2.4170221126000001e+01 - 0.0000000000000000e+00 0.0000000000000000e+00 2.2041158952999999e+03 2.7142512751999999e+02 -2.4170221126000001e+01 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 2.7142512751999999e+02 6.9571871693000003e+03 -3.3373259730999997e+01 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 -2.4170221126000001e+01 -3.3373259730999997e+01 5.8238848565999997e+03 0.0000000000000000e+00 - -2.7142512751999999e+02 2.4170221126000001e+01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 1.2781072026000000e+04 - - 0.100000 - 2.4031873314000001e+09 7.4245162653999999e+07 -5.8923454737999998e+07 -1.3127022075000000e+08 1.2832049021000000e+07 -1.4313381863000000e+08 - 7.4245162653999999e+07 3.9657475948000002e+09 1.8557857279999999e+06 -9.2277206769000005e+05 -5.0133546035000002e+06 8.3989388888999999e+07 - -5.8923454737999998e+07 1.8557857279999999e+06 2.6858344563000000e+10 6.3640009755000000e+09 -4.8262776811000001e+08 1.2976757009999999e+08 - -1.3127022075000000e+08 -9.2277206769000005e+05 6.3640009755000000e+09 6.6747454769000000e+10 -5.1295607331999998e+09 2.3615395219999999e+08 - 1.2832049021000000e+07 -5.0133546035000002e+06 -4.8262776811000001e+08 -5.1295607331999998e+09 5.2568641475000000e+10 -5.6793124452999998e+06 - -1.4313381863000000e+08 8.3989388888999999e+07 1.2976757009999999e+08 2.3615395219999999e+08 -5.6793124452999998e+06 2.4960016204000000e+10 - - 1.6712070776000000e+03 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 -4.4390931882000001e+02 - 0.0000000000000000e+00 1.6712070776000000e+03 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 2.9926378696000000e+01 - 0.0000000000000000e+00 0.0000000000000000e+00 1.6712070776000000e+03 4.4390931882000001e+02 -2.9926378696000000e+01 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 4.4390931882000001e+02 4.6542513771000004e+03 -1.1944138900999999e+02 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 -2.9926378696000000e+01 -1.1944138900999999e+02 2.7817167724000001e+03 0.0000000000000000e+00 - -4.4390931882000001e+02 2.9926378696000000e+01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 7.4359681495000004e+03 - - 0.150000 - 9.0679270074000001e+08 7.1783373003000006e+07 -7.4044817680999994e+07 -1.9515534075000000e+08 2.8552463930000000e+07 1.7060326012000000e+08 - 7.1783373003000006e+07 2.3314893804000001e+09 -8.1419598595000003e+06 1.1336559359999999e+07 -6.2489422181999996e+07 3.4018349965000004e+07 - -7.4044817680999994e+07 -8.1419598595000003e+06 2.1210763465000000e+10 6.0287195731999998e+09 -1.6872879959000000e+08 2.3595622538999999e+08 - -1.9515534075000000e+08 1.1336559359999999e+07 6.0287195731999998e+09 3.9997060590000000e+10 -4.8997378675000000e+09 4.1415048487000000e+08 - 2.8552463930000000e+07 -6.2489422181999996e+07 -1.6872879959000000e+08 -4.8997378675000000e+09 2.7814376678000000e+10 -1.9179116116000000e+07 - 1.7060326012000000e+08 3.4018349965000004e+07 2.3595622538999999e+08 4.1415048487000000e+08 -1.9179116116000000e+07 7.6704118448000002e+09 - - 9.8068419680000000e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 -4.2712470944000000e+02 - 0.0000000000000000e+00 9.8068419680000000e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 1.8552555125000001e+01 - 0.0000000000000000e+00 0.0000000000000000e+00 9.8068419680000000e+02 4.2712470944000000e+02 -1.8552555125000001e+01 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 4.2712470944000000e+02 2.7369259900000002e+03 -1.0140747589000000e+02 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 -1.8552555125000001e+01 -1.0140747589000000e+02 9.4045127984999999e+02 0.0000000000000000e+00 - -4.2712470944000000e+02 1.8552555125000001e+01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 3.6773772697999998e+03 - - 0.200000 - 4.4875261644000000e+08 4.3616584689000003e+07 -6.1190618108000003e+07 -1.2386407687000000e+08 1.6631414532000000e+07 7.0570049469999999e+07 - 4.3616584689000003e+07 1.0149256158000000e+09 -7.3133623804000001e+06 1.4910947710999999e+07 -9.0627753312999994e+07 3.4323628478000000e+07 - -6.1190618108000003e+07 -7.3133623804000001e+06 1.8883258348000000e+10 4.5758212826999998e+09 -2.8225843602999997e+08 1.8791743661000001e+08 - -1.2386407687000000e+08 1.4910947710999999e+07 4.5758212826999998e+09 2.4308699290000000e+10 -2.8853517304000001e+09 3.0756077491000003e+08 - 1.6631414532000000e+07 -9.0627753312999994e+07 -2.8225843602999997e+08 -2.8853517304000001e+09 2.1638926709000000e+10 -9.6189813441000003e+06 - 7.0570049469999999e+07 3.4323628478000000e+07 1.8791743661000001e+08 3.0756077491000003e+08 -9.6189813441000003e+06 2.6113972210000000e+09 - - 6.0666066332000003e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 -3.3340866762000002e+02 - 0.0000000000000000e+00 6.0666066332000003e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 1.5788648778000001e+01 - 0.0000000000000000e+00 0.0000000000000000e+00 6.0666066332000003e+02 3.3340866762000002e+02 -1.5788648778000001e+01 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 3.3340866762000002e+02 1.6309254443000000e+03 -6.3602442058999998e+01 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 -1.5788648778000001e+01 -6.3602442058999998e+01 4.6981683225000000e+02 0.0000000000000000e+00 - -3.3340866762000002e+02 1.5788648778000001e+01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 2.1007422765000001e+03 - - 0.250000 - 2.8072081675999999e+08 2.7429716309999999e+07 -5.3791011803999998e+07 -9.0577889878999993e+07 8.5773351194000002e+06 4.5293945458999999e+07 - 2.7429716309999999e+07 5.5356632377999997e+08 -6.9236061853999998e+06 3.2666786148000001e+06 -6.1733702818999998e+07 4.0401546254000001e+07 - -5.3791011803999998e+07 -6.9236061853999998e+06 1.8928721569000000e+10 4.1569322419000001e+09 -2.3595885088999999e+08 1.2823289395999999e+08 - -9.0577889878999993e+07 3.2666786148000001e+06 4.1569322419000001e+09 1.9772889899000000e+10 -1.7693786825000000e+09 2.4644634903999999e+08 - 8.5773351194000002e+06 -6.1733702818999998e+07 -2.3595885088999999e+08 -1.7693786825000000e+09 1.5353044678000000e+10 -5.3454390472999997e+06 - 4.5293945458999999e+07 4.0401546254000001e+07 1.2823289395999999e+08 2.4644634903999999e+08 -5.3454390472999997e+06 1.0260318231000000e+09 - - 4.7690739313000000e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 -2.9310667175999998e+02 - 0.0000000000000000e+00 4.7690739313000000e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 1.4287126102000000e+01 - 0.0000000000000000e+00 0.0000000000000000e+00 4.7690739313000000e+02 2.9310667175999998e+02 -1.4287126102000000e+01 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 2.9310667175999998e+02 1.2474296161000000e+03 -4.5066773380999997e+01 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 -1.4287126102000000e+01 -4.5066773380999997e+01 2.6450926969000000e+02 0.0000000000000000e+00 - -2.9310667175999998e+02 1.4287126102000000e+01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 1.5119388858000000e+03 - - 0.300000 - 2.2334704344999999e+08 1.3447678317000000e+07 -5.5554370652999997e+07 -8.0694192193000004e+07 4.4905992836999996e+06 2.7094596614000000e+07 - 1.3447678317000000e+07 3.8419832579000002e+08 -3.3688856924000001e+06 1.6643441961000001e+06 -4.9407763880999997e+07 3.0155587206999999e+07 - -5.5554370652999997e+07 -3.3688856924000001e+06 1.9825241678000000e+10 4.4667581271999998e+09 -2.1648425953000000e+08 1.0940317166000000e+08 - -8.0694192193000004e+07 1.6643441961000001e+06 4.4667581271999998e+09 1.8680234954000000e+10 -9.8732231813999999e+08 2.1787883534000000e+08 - 4.4905992836999996e+06 -4.9407763880999997e+07 -2.1648425953000000e+08 -9.8732231813999999e+08 1.3039262973000000e+10 -4.4524942631999999e+06 - 2.7094596614000000e+07 3.0155587206999999e+07 1.0940317166000000e+08 2.1787883534000000e+08 -4.4524942631999999e+06 6.1932076416999996e+08 - - 4.3416105263999998e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 -2.8425098801000001e+02 - 0.0000000000000000e+00 4.3416105263999998e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 1.1547358941000001e+01 - 0.0000000000000000e+00 0.0000000000000000e+00 4.3416105263999998e+02 2.8425098801000001e+02 -1.1547358941000001e+01 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 2.8425098801000001e+02 1.0850901613000001e+03 -3.1091777955000001e+01 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 -1.1547358941000001e+01 -3.1091777955000001e+01 1.9975678597999999e+02 0.0000000000000000e+00 - -2.8425098801000001e+02 1.1547358941000001e+01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 1.2848469473000000e+03 - - 0.350000 - 1.8750843983000001e+08 4.7434107164000003e+06 -4.8343553248000003e+07 -7.0828055276999995e+07 1.7982204338000000e+06 2.1832437488000002e+07 - 4.7434107164000003e+06 3.4830745445999998e+08 6.9280961918000004e+05 4.5234330219999999e+06 -3.4693151843000002e+07 2.6664997728000000e+07 - -4.8343553248000003e+07 6.9280961918000004e+05 2.0703459547000000e+10 5.1155476706999998e+09 -1.8631430261000001e+08 8.8788258746999994e+07 - -7.0828055276999995e+07 4.5234330219999999e+06 5.1155476706999998e+09 1.8899173593000000e+10 -4.3880708856999999e+08 1.8538557906000000e+08 - 1.7982204338000000e+06 -3.4693151843000002e+07 -1.8631430261000001e+08 -4.3880708856999999e+08 1.0556796249000000e+10 -3.0353536727999998e+06 - 2.1832437488000002e+07 2.6664997728000000e+07 8.8788258746999994e+07 1.8538557906000000e+08 -3.0353536727999998e+06 4.5670607506999999e+08 - - 4.2421580561000002e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 -2.9520298902000002e+02 - 0.0000000000000000e+00 4.2421580561000002e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 9.3446667771000005e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 4.2421580561000002e+02 2.9520298902000002e+02 -9.3446667771000005e+00 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 2.9520298902000002e+02 1.0240564349000001e+03 -2.1507019011000001e+01 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 -9.3446667771000005e+00 -2.1507019011000001e+01 1.5310149835000001e+02 0.0000000000000000e+00 - -2.9520298902000002e+02 9.3446667771000005e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 1.1771579333000000e+03 - - 0.400000 - 1.6047782522999999e+08 2.5388412543000001e+06 -3.4062172715000004e+07 -4.8177568817000002e+07 8.8888365310000000e+05 1.8596924901000001e+07 - 2.5388412543000001e+06 3.2888406347000003e+08 4.3539888265000004e+06 1.9604119210000001e+07 -2.2460588557999998e+07 1.3088871514000000e+07 - -3.4062172715000004e+07 4.3539888265000004e+06 1.9104411739000000e+10 4.4971248598000002e+09 -1.1877964955000000e+08 5.6048441723999999e+07 - -4.8177568817000002e+07 1.9604119210000001e+07 4.4971248598000002e+09 1.6095765482000000e+10 -2.7837423794000000e+08 1.1917741537000000e+08 - 8.8888365310000000e+05 -2.2460588557999998e+07 -1.1877964955000000e+08 -2.7837423794000000e+08 8.3926194080000000e+09 -1.7869603999000001e+06 - 1.8596924901000001e+07 1.3088871514000000e+07 5.6048441723999999e+07 1.1917741537000000e+08 -1.7869603999000001e+06 3.5480450149000001e+08 - - 4.0064426630999998e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 -2.6586696054999999e+02 - 0.0000000000000000e+00 4.0064426630999998e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 6.4097879477999999e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 4.0064426630999998e+02 2.6586696054999999e+02 -6.4097879477999999e+00 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 2.6586696054999999e+02 8.5829076979000001e+02 -1.6722454831000000e+01 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 -6.4097879477999999e+00 -1.6722454831000000e+01 1.1827514899000001e+02 0.0000000000000000e+00 - -2.6586696054999999e+02 6.4097879477999999e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 9.7656591877000005e+02 - - 0.450000 - 1.3713355400000000e+08 3.4135088078000001e+06 -2.0558824640999999e+07 -2.7409318067000002e+07 4.1782630551999999e+05 1.6474908615000000e+07 - 3.4135088078000001e+06 3.3123565902999997e+08 6.8684566156000004e+05 4.7790759845000003e+06 -1.3311475300000001e+07 1.2700634419000000e+07 - -2.0558824640999999e+07 6.8684566156000004e+05 1.9278068946000000e+10 3.5505433059000001e+09 -3.2730426942000002e+07 3.3926953747000001e+07 - -2.7409318067000002e+07 4.7790759845000003e+06 3.5505433059000001e+09 1.1926544805000000e+10 -1.8407909019999999e+08 6.3091197277000003e+07 - 4.1782630551999999e+05 -1.3311475300000001e+07 -3.2730426942000002e+07 -1.8407909019999999e+08 6.4821239165000000e+09 -9.0437065335000004e+05 - 1.6474908615000000e+07 1.2700634419000000e+07 3.3926953747000001e+07 6.3091197277000003e+07 -9.0437065335000004e+05 2.7535130629000002e+08 - - 3.6547262165000001e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 -2.0996524237000000e+02 - 0.0000000000000000e+00 3.6547262165000001e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 4.6589532624999999e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 3.6547262165000001e+02 2.0996524237000000e+02 -4.6589532624999999e+00 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 2.0996524237000000e+02 6.3408896159000005e+02 -1.2548039064999999e+01 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 -4.6589532624999999e+00 -1.2548039064999999e+01 8.9470594765000001e+01 0.0000000000000000e+00 - -2.0996524237000000e+02 4.6589532624999999e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 7.2355955634999998e+02 - - 0.500000 - 1.1520796803000000e+08 3.5740442618999998e+06 -1.3609537523000000e+07 -1.7047753719999999e+07 2.5120666803999999e+05 1.4679461338000000e+07 - 3.5740442618999998e+06 3.3753098569999999e+08 2.2512669880000001e+05 -5.1162871535000001e+05 -8.7679376871000007e+06 1.5587995848999999e+07 - -1.3609537523000000e+07 2.2512669880000001e+05 1.8548933390000000e+10 2.5866000739000001e+09 -7.3669102905000001e+07 2.2220613824000001e+07 - -1.7047753719999999e+07 -5.1162871535000001e+05 2.5866000739000001e+09 8.3210520526000004e+09 -1.0817572931000000e+08 3.5538366060000002e+07 - 2.5120666803999999e+05 -8.7679376871000007e+06 -7.3669102905000001e+07 -1.0817572931000000e+08 4.8271284071000004e+09 -6.0465612887999997e+05 - 1.4679461338000000e+07 1.5587995848999999e+07 2.2220613824000001e+07 3.5538366060000002e+07 -6.0465612887999997e+05 2.0906418147000000e+08 - - 3.2928828972999997e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 -1.5777473846999999e+02 - 0.0000000000000000e+00 3.2928828972999997e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 4.8197585656999999e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 3.2928828972999997e+02 1.5777473846999999e+02 -4.8197585656999999e+00 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 1.5777473846999999e+02 4.4591133460999998e+02 -8.9904436487999995e+00 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 -4.8197585656999999e+00 -8.9904436487999995e+00 6.5504149738999999e+01 0.0000000000000000e+00 - -1.5777473846999999e+02 4.8197585656999999e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 5.1141548434999999e+02 - - 0.550000 - 9.6467011657000005e+07 2.2457118941000002e+06 -8.4274022239999995e+06 -1.0250371884000000e+07 1.2415670527999999e+05 1.2696741321000000e+07 - 2.2457118941000002e+06 3.4469413475000000e+08 1.7646903403000000e+06 3.2750044122000001e+06 -5.5308002078000000e+06 4.0429699517999999e+06 - -8.4274022239999995e+06 1.7646903403000000e+06 1.6622628598000000e+10 1.8631698103000000e+09 -1.5986116706000000e+08 1.2325566801000001e+07 - -1.0250371884000000e+07 3.2750044122000001e+06 1.8631698103000000e+09 5.9971336136999998e+09 -4.8048857623999998e+07 1.9193325936000001e+07 - 1.2415670527999999e+05 -5.5308002078000000e+06 -1.5986116706000000e+08 -4.8048857623999998e+07 3.5701540331999998e+09 -4.0869668411999999e+05 - 1.2696741321000000e+07 4.0429699517999999e+06 1.2325566801000001e+07 1.9193325936000001e+07 -4.0869668411999999e+05 1.5796873102000001e+08 - - 2.9981748379999999e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 -1.2112823803000001e+02 - 0.0000000000000000e+00 2.9981748379999999e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 4.9683646858000001e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 2.9981748379999999e+02 1.2112823803000001e+02 -4.9683646858000001e+00 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 1.2112823803000001e+02 3.2226227340000003e+02 -6.4207541478000003e+00 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 -4.9683646858000001e+00 -6.4207541478000003e+00 4.7775347357000001e+01 0.0000000000000000e+00 - -1.2112823803000001e+02 4.9683646858000001e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 3.7003762075999998e+02 - - 0.600000 - 8.0313388391000003e+07 3.0576877716000001e+06 -4.9340406523000002e+06 -6.0312837122999998e+06 1.0184250015000001e+05 1.1425447971000001e+07 - 3.0576877716000001e+06 3.8030394695999998e+08 2.3476064942999999e+06 1.0519831812000000e+06 -3.4178701115999999e+06 8.8125248970999997e+06 - -4.9340406523000002e+06 2.3476064942999999e+06 1.6372608371000000e+10 1.3027459067000000e+09 -2.6148757911000001e+08 7.9322361441000002e+06 - -6.0312837122999998e+06 1.0519831812000000e+06 1.3027459067000000e+09 4.0953627031999998e+09 -3.6862523983999997e+07 9.6808672774999999e+06 - 1.0184250015000001e+05 -3.4178701115999999e+06 -2.6148757911000001e+08 -3.6862523983999997e+07 2.6344077321999998e+09 -3.0745473697000003e+05 - 1.1425447971000001e+07 8.8125248970999997e+06 7.9322361441000002e+06 9.6808672774999999e+06 -3.0745473697000003e+05 1.1998899775000000e+08 - - 2.7097991674999997e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 -8.6675510721999999e+01 - 0.0000000000000000e+00 2.7097991674999997e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 5.8836984569000004e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 2.7097991674999997e+02 8.6675510721999999e+01 -5.8836984569000004e+00 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 8.6675510721999999e+01 2.1930711817000000e+02 -4.6977436150000003e+00 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 -5.8836984569000004e+00 -4.6977436150000003e+00 3.4836089972000003e+01 0.0000000000000000e+00 - -8.6675510721999999e+01 5.8836984569000004e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 2.5414320814999999e+02 - - 0.650000 - 6.6406987479999997e+07 4.0806485926999999e+06 -2.9783046140999999e+06 -3.9330139145999998e+06 8.6190954702999996e+04 1.0442748527000001e+07 - 4.0806485926999999e+06 4.4522080201999998e+08 3.8356519218000001e+06 3.0171741957999999e+05 -2.2312910797999999e+06 1.4838901524000000e+07 - -2.9783046140999999e+06 3.8356519218000001e+06 1.5390339910000000e+10 8.0494743775000000e+08 -3.4225879641000003e+08 5.6423046731000002e+06 - -3.9330139145999998e+06 3.0171741957999999e+05 8.0494743775000000e+08 2.7256609914000001e+09 -2.8280843895000000e+07 4.9944657255999995e+06 - 8.6190954702999996e+04 -2.2312910797999999e+06 -3.4225879641000003e+08 -2.8280843895000000e+07 1.9070773187000000e+09 -2.4832681609000001e+05 - 1.0442748527000001e+07 1.4838901524000000e+07 5.6423046731000002e+06 4.9944657255999995e+06 -2.4832681609000001e+05 9.0687763310000002e+07 - - 2.4291095505000001e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 -5.9130394832000000e+01 - 0.0000000000000000e+00 2.4291095505000001e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 6.4887480899999996e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 2.4291095505000001e+02 5.9130394832000000e+01 -6.4887480899999996e+00 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 5.9130394832000000e+01 1.4493425407999999e+02 -3.3303070787000002e+00 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 -6.4887480899999996e+00 -3.3303070787000002e+00 2.5025257034999999e+01 0.0000000000000000e+00 - -5.9130394832000000e+01 6.4887480899999996e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 1.6995951110999999e+02 - - 0.700000 - 5.4994082805000000e+07 4.8184355719999997e+06 -1.9060491206000000e+06 -3.0422898568000002e+06 6.7566693723999997e+04 9.4616666366000008e+06 - 4.8184355719999997e+06 5.3578712462000000e+08 7.7137833485000003e+06 7.3779884759000002e+05 -1.6113454435000001e+06 9.2224275599000007e+06 - -1.9060491206000000e+06 7.7137833485000003e+06 1.3857898560000000e+10 4.3414821520999998e+08 -4.0817478385000002e+08 4.3601461103999997e+06 - -3.0422898568000002e+06 7.3779884759000002e+05 4.3414821520999998e+08 1.8016640499000001e+09 -1.5971826897000000e+07 2.6683363961999998e+06 - 6.7566693723999997e+04 -1.6113454435000001e+06 -4.0817478385000002e+08 -1.5971826897000000e+07 1.3472383192000000e+09 -2.1109678135000000e+05 - 9.4616666366000008e+06 9.2224275599000007e+06 4.3601461103999997e+06 2.6683363961999998e+06 -2.1109678135000000e+05 6.8021703585999995e+07 - - 2.1312784884000001e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 -3.8501453607000002e+01 - 0.0000000000000000e+00 2.1312784884000001e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 6.7693416248000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 2.1312784884000001e+02 3.8501453607000002e+01 -6.7693416248000000e+00 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 3.8501453607000002e+01 9.4046538154999993e+01 -2.2090212465999999e+00 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 -6.7693416248000000e+00 -2.2090212465999999e+00 1.7646972125000001e+01 0.0000000000000000e+00 - -3.8501453607000002e+01 6.7693416248000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 1.1169351028000000e+02 - - 0.750000 - 4.4532689961000003e+07 3.7855724704999998e+06 -1.5672704434000000e+06 -3.0235377555999998e+06 5.1909411982999998e+04 8.2055638525999999e+06 - 3.7855724704999998e+06 4.8748771911000001e+08 1.0253502158000000e+07 6.9563826246000000e+05 -1.2525839791000001e+06 1.3147165897000000e+07 - -1.5672704434000000e+06 1.0253502158000000e+07 1.0962725249000000e+10 3.1619657647000003e+08 -3.8508825074000001e+08 4.2235395332000004e+06 - -3.0235377555999998e+06 6.9563826246000000e+05 3.1619657647000003e+08 1.4079102655999999e+09 -6.0152521610000003e+06 2.3596887472999999e+06 - 5.1909411982999998e+04 -1.2525839791000001e+06 -3.8508825074000001e+08 -6.0152521610000003e+06 8.3311139392999995e+08 -2.0981475915000000e+05 - 8.2055638525999999e+06 1.3147165897000000e+07 4.2235395332000004e+06 2.3596887472999999e+06 -2.0981475915000000e+05 5.0462376886000000e+07 - - 1.7380336410999999e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 -2.9814273096000001e+01 - 0.0000000000000000e+00 1.7380336410999999e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 6.2885469630999999e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 1.7380336410999999e+02 2.9814273096000001e+01 -6.2885469630999999e+00 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 2.9814273096000001e+01 7.1021270795000007e+01 -1.5419452683999999e+00 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 -6.2885469630999999e+00 -1.5419452683999999e+00 1.1183459569000000e+01 0.0000000000000000e+00 - -2.9814273096000001e+01 6.2885469630999999e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 8.2204730364000000e+01 - - 0.800000 - 3.8356631333999999e+07 3.3395198176000002e+06 -5.8233507851000002e+05 -1.3995981921999999e+06 2.0926886187000000e+04 8.3315559623999996e+06 - 3.3395198176000002e+06 4.0804832088000000e+08 4.2730116397000002e+06 1.4832464395000000e+05 -3.7879300182000000e+05 1.5420843998000000e+07 - -5.8233507851000002e+05 4.2730116397000002e+06 7.2907120376000004e+09 1.7317709434000000e+08 -2.5857627581999999e+08 1.8291024441000000e+06 - -1.3995981921999999e+06 1.4832464395000000e+05 1.7317709434000000e+08 1.0152418234000000e+09 -2.3232957521000002e+06 9.3622078683000000e+05 - 2.0926886187000000e+04 -3.7879300182000000e+05 -2.5857627581999999e+08 -2.3232957521000002e+06 4.5975706705000001e+08 -8.0702550625999997e+04 - 8.3315559623999996e+06 1.5420843998000000e+07 1.8291024441000000e+06 9.3622078683000000e+05 -8.0702550625999997e+04 4.1023688162000000e+07 - - 1.3018754444000001e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 -1.9640799484999999e+01 - 0.0000000000000000e+00 1.3018754444000001e+02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 4.6392499095000002e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 1.3018754444000001e+02 1.9640799484999999e+01 -4.6392499095000002e+00 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 1.9640799484999999e+01 5.1908965705999996e+01 -1.0818028347999999e+00 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 -4.6392499095000002e+00 -1.0818028347999999e+00 6.6180715242000003e+00 0.0000000000000000e+00 - -1.9640799484999999e+01 4.6392499095000002e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 5.8527037229999998e+01 - - 0.850000 - 2.7037664366000000e+07 1.8052328012000001e+06 -2.5623480766000001e+04 -5.2892151980000002e+04 7.1417807922999998e+02 3.9289505153000001e+06 - 1.8052328012000001e+06 2.6151739238999999e+08 1.5130030108999999e+05 1.1177750855000000e+04 -1.1109355777999999e+04 6.6767695229000002e+06 - -2.5623480766000001e+04 1.5130030108999999e+05 4.4728871703999996e+09 1.5710167424000001e+08 -1.4782098644999999e+08 6.4668055868000003e+04 - -5.2892151980000002e+04 1.1177750855000000e+04 1.5710167424000001e+08 5.7073460941999996e+08 -1.8406614483000000e+06 3.4020500175000001e+04 - 7.1417807922999998e+02 -1.1109355777999999e+04 -1.4782098644999999e+08 -1.8406614483000000e+06 2.4050708547000000e+08 -2.5035405931999999e+03 - 3.9289505153000001e+06 6.6767695229000002e+06 6.4668055868000003e+04 3.4020500175000001e+04 -2.5035405931999999e+03 2.2859208517999999e+07 - - 8.2157678352000005e+01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 -1.4401316204000000e+01 - 0.0000000000000000e+00 8.2157678352000005e+01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 2.7594394443999999e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 8.2157678352000005e+01 1.4401316204000000e+01 -2.7594394443999999e+00 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 1.4401316204000000e+01 2.8824418436999999e+01 -5.9589664765999995e-01 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 -2.7594394443999999e+00 -5.9589664765999995e-01 3.5760270598999999e+00 0.0000000000000000e+00 - -1.4401316204000000e+01 2.7594394443999999e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 3.2400445497000000e+01 - - 0.900000 - 1.9381963243999999e+07 1.2022489213000000e+06 8.5162952946000005e+04 1.6884186919999999e+05 -2.6218548114000000e+03 2.2564549682000000e+06 - 1.2022489213000000e+06 1.6659039418000001e+08 -4.2393765831999999e+05 -2.6881476429999999e+04 1.9931926370000001e+04 4.5307964637000002e+06 - 8.5162952946000005e+04 -4.2393765831999999e+05 2.6083808581999998e+09 8.6127992585999995e+07 -7.9225450568000004e+07 -1.9291306426000001e+05 - 1.6884186919999999e+05 -2.6881476429999999e+04 8.6127992585999995e+07 3.0656377656999999e+08 -1.5495510183000001e+06 -8.7398993895000007e+04 - -2.6218548114000000e+03 1.9931926370000001e+04 -7.9225450568000004e+07 -1.5495510183000001e+06 1.1218305862000000e+08 6.3732981298000004e+03 - 2.2564549682000000e+06 4.5307964637000002e+06 -1.9291306426000001e+05 -8.7398993895000007e+04 6.3732981298000004e+03 1.3040168952000000e+07 - - 5.1585262628000002e+01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 -8.1387186000000007e+00 - 0.0000000000000000e+00 5.1585262628000002e+01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 1.5775199094000001e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 5.1585262628000002e+01 8.1387186000000007e+00 -1.5775199094000001e+00 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 8.1387186000000007e+00 1.5461936073000000e+01 -3.1382283168000003e-01 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 -1.5775199094000001e+00 -3.1382283168000003e-01 1.7610994122999999e+00 0.0000000000000000e+00 - -8.1387186000000007e+00 1.5775199094000001e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 1.7223035486000001e+01 - - 0.950000 - 1.3794705691000000e+07 6.0897303269999998e+05 1.2636423659000000e+05 1.8950537085000001e+05 -4.4924517849000003e+03 1.4322905356999999e+06 - 6.0897303269999998e+05 9.8206222994000003e+07 -3.1433398976999999e+05 -1.9452934992999999e+04 7.8952360753000003e+03 2.9592355263999999e+06 - 1.2636423659000000e+05 -3.1433398976999999e+05 1.1385780670000000e+09 5.8589137399999999e+07 -3.0646733372000001e+07 -1.9244792384999999e+05 - 1.8950537085000001e+05 -1.9452934992999999e+04 5.8589137399999999e+07 1.5056671188999999e+08 -1.9984993668000000e+06 -1.0513157257999999e+05 - -4.4924517849000003e+03 7.8952360753000003e+03 -3.0646733372000001e+07 -1.9984993668000000e+06 3.6006313774999999e+07 5.4597793146000004e+03 - 1.4322905356999999e+06 2.9592355263999999e+06 -1.9244792384999999e+05 -1.0513157257999999e+05 5.4597793146000004e+03 7.3791966996999998e+06 - - 3.2226395472000000e+01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 -5.4999391553999999e+00 - 0.0000000000000000e+00 3.2226395472000000e+01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 8.7531323553999996e-01 - 0.0000000000000000e+00 0.0000000000000000e+00 3.2226395472000000e+01 5.4999391553999999e+00 -8.7531323553999996e-01 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 5.4999391553999999e+00 9.3422586507999998e+00 -2.0413469801000000e-01 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 -8.7531323553999996e-01 -2.0413469801000000e-01 7.5684283515999995e-01 0.0000000000000000e+00 - -5.4999391553999999e+00 8.7531323553999996e-01 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 1.0099101486000000e+01 - - 1.000000 - 1.0698098555000001e+06 2.5725377725000001e+05 4.0396003678000002e+03 2.2014856988000001e+03 -6.0977746271000001e+01 8.1254414661000003e+04 - 2.5725377725000001e+05 1.6520197075999999e+07 -1.7886818336000000e+04 -2.9370253184000001e+02 -1.3330425389000000e+01 1.0806092268000000e+05 - 4.0396003678000002e+03 -1.7886818336000000e+04 1.2014872833000000e+08 1.9141905636000000e+06 -7.0470167760000005e+05 -2.4532477761000000e+03 - 2.2014856988000001e+03 -2.9370253184000001e+02 1.9141905636000000e+06 1.5788371902000001e+06 -3.3245422585000000e+04 -2.8659551955000001e+02 - -6.0977746271000001e+01 -1.3330425389000000e+01 -7.0470167760000005e+05 -3.3245422585000000e+04 1.8154183775000001e+05 2.7223926702000000e+01 - 8.1254414661000003e+04 1.0806092268000000e+05 -2.4532477761000000e+03 -2.8659551955000001e+02 2.7223926702000000e+01 7.4377927549999993e+04 - - 5.6225720411999998e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 -2.2219707415000001e-01 - 0.0000000000000000e+00 5.6225720411999998e+00 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 3.4429766177999997e-02 - 0.0000000000000000e+00 0.0000000000000000e+00 5.6225720411999998e+00 2.2219707415000001e-01 -3.4429766177999997e-02 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 2.2219707415000001e-01 1.0876819022000001e-01 -2.8261744220999998e-03 0.0000000000000000e+00 - 0.0000000000000000e+00 0.0000000000000000e+00 -3.4429766177999997e-02 -2.8261744220999998e-03 7.0468259357999996e-03 0.0000000000000000e+00 - -2.2219707415000001e-01 3.4429766177999997e-02 0.0000000000000000e+00 0.0000000000000000e+00 0.0000000000000000e+00 1.1581501616000001e-01 - diff --git a/examples/01_aeroelasticse/run_OLAF.py b/examples/01_aeroelasticse/run_OLAF.py index 1bccf5fea..8f21e07fd 100644 --- a/examples/01_aeroelasticse/run_OLAF.py +++ b/examples/01_aeroelasticse/run_OLAF.py @@ -58,7 +58,7 @@ case_inputs[("ElastoDyn","BlPitch1")] = {'vals': pitch_init, 'group': 1} case_inputs[("ElastoDyn","BlPitch2")] = case_inputs[("ElastoDyn","BlPitch1")] case_inputs[("ElastoDyn","BlPitch3")] = case_inputs[("ElastoDyn","BlPitch1")] -dt_wanted, nNWPanel, nFWPanel, nFWPanelFree = OLAFParams(omega_init) +dt_wanted, tMax, nNWPanel, nFWPanel, nFWPanelFree = OLAFParams(omega_init) dt_olaf = np.zeros_like(dt_wanted) dt = case_inputs[("Fst","DT")]["vals"] n_dt = dt_wanted / dt diff --git a/weis/aeroelasticse/FAST_reader.py b/weis/aeroelasticse/FAST_reader.py index 011238ed7..538e91748 100644 --- a/weis/aeroelasticse/FAST_reader.py +++ b/weis/aeroelasticse/FAST_reader.py @@ -428,8 +428,9 @@ def execute(self): if self.fst_vt['Fst']['CompMooring'] == 3: # MoorDyn implimented self.read_MoorDyn() - if self.fst_vt['Fst']['CompElast'] == 2: # BeamDyn read assumes all 3 blades are the same - self.read_BeamDyn() + bd_file = os.path.join(self.FAST_directory, self.fst_vt['Fst']['BDBldFile(1)']) + if os.path.exists(bd_file): + self.read_BeamDyn(bd_file) def read_MainInput(self): # Main FAST v8.16-v8.17 Input File @@ -705,12 +706,9 @@ def read_ElastoDyn(self): f.close() - def read_BeamDyn(self): + def read_BeamDyn(self, bd_file): # BeamDyn Input File - - bd_file = os.path.join(self.FAST_directory, self.fst_vt['Fst']['BDBldFile(1)']) f = open(bd_file) - f.readline() f.readline() f.readline() @@ -1126,10 +1124,9 @@ def read_AeroDyn15(self): self.read_AeroDyn15Blade() self.read_AeroDyn15Polar() self.read_AeroDyn15Coord() - if self.fst_vt['AeroDyn15']['WakeMod'] == 3: - if self.fst_vt['AeroDyn15']['AFAeroMod'] == 2: - raise Exception('OLAF is called with unsteady airfoil aerodynamics, but OLAF currently only supports AFAeroMod == 1') - self.read_AeroDyn15OLAF() + olaf_filename = os.path.join(self.FAST_directory, self.fst_vt['AeroDyn15']['OLAFInputFileName']) + if os.path.exists(olaf_filename): + self.read_AeroDyn15OLAF(olaf_filename) def read_AeroDyn15Blade(self): # AeroDyn v5.00 Blade Definition File @@ -1274,10 +1271,9 @@ def read_AeroDyn15Coord(self): f.close() - def read_AeroDyn15OLAF(self): + def read_AeroDyn15OLAF(self, olaf_filename): self.fst_vt['AeroDyn15']['OLAF'] = {} - olaf_filename = os.path.join(self.FAST_directory, self.fst_vt['AeroDyn15']['OLAFInputFileName']) f = open(olaf_filename) f.readline() f.readline() @@ -1323,7 +1319,6 @@ def read_AeroDyn15OLAF(self): f.readline() f.close() - def read_ServoDyn(self): # ServoDyn v1.05 Input File # Currently no differences between FASTv8.16 and OpenFAST. diff --git a/weis/aeroelasticse/FAST_writer.py b/weis/aeroelasticse/FAST_writer.py index c48fcadf7..49214b52c 100644 --- a/weis/aeroelasticse/FAST_writer.py +++ b/weis/aeroelasticse/FAST_writer.py @@ -878,6 +878,8 @@ def write_AeroDyn15(self): self.write_AeroDyn15Coord() if self.fst_vt['AeroDyn15']['WakeMod'] == 3: + if self.fst_vt['AeroDyn15']['AFAeroMod'] == 2: + raise Exception('OLAF is called with unsteady airfoil aerodynamics, but OLAF currently only supports AFAeroMod == 1') self.write_OLAF() # Generate AeroDyn v15.03 input file diff --git a/weis/aeroelasticse/openmdao_openfast.py b/weis/aeroelasticse/openmdao_openfast.py index 21466ecb9..8b2161166 100644 --- a/weis/aeroelasticse/openmdao_openfast.py +++ b/weis/aeroelasticse/openmdao_openfast.py @@ -1994,36 +1994,49 @@ def writeCpsurfaces(self, inputs): def RayleighCDF(x, xbar=10.): return 1.0 - np.exp(-np.pi/4.0*(x/xbar)**2) -def OLAFParams(omega_rpm, deltaPsiDeg=6, nNWrot=2, nFWrot=10, nFWrotFree=3, nPerAzimuth=None): +def OLAFParams(omega_rpm, deltaPsiDeg=6, nNWrot=2, nFWrot=10, nFWrotFree=3, nPerRot=None, totalRot=None, show=False): """ Computes recommended time step and wake length based on the rotational speed in RPM - + INPUTS: - omega_rpm: rotational speed in RPM - deltaPsiDeg : azimuthal discretization in deg - nNWrot : number of near wake rotations - nFWrot : total number of far wake rotations - nFWrotFree : number of far wake rotations that are free - - deltaPsiDeg - nPerAzimuth - 5 72 - 6 60 - 7 51.5 - 8 45 + + deltaPsiDeg - nPerRot + 5 72 + 6 60 + 7 51.5 + 8 45 """ + omega_rpm = np.asarray(omega_rpm) omega = omega_rpm*2*np.pi/60 T = 2*np.pi/omega - if nPerAzimuth is not None: - dt_wanted = np.around(T/nPerAzimuth,2) + if nPerRot is not None: + dt_wanted = np.around(T/nPerRot,4) else: - dt_wanted = np.around(deltaPsiDeg/(6*omega_rpm), 4) - nPerAzimuth = int(2*np.pi /(deltaPsiDeg*np.pi/180)) - - nNWPanel = nNWrot*nPerAzimuth - nFWPanel = nFWrot*nPerAzimuth - nFWPanelFree = nFWrotFree*nPerAzimuth - - return dt_wanted, nNWPanel, nFWPanel, nFWPanelFree + dt_wanted = np.around(deltaPsiDeg/(6*omega_rpm),4) + nPerRot = int(2*np.pi /(deltaPsiDeg*np.pi/180)) + + nNWPanel = nNWrot*nPerRot + nFWPanel = nFWrot*nPerRot + nFWPanelFree = nFWrotFree*nPerRot + + if totalRot is None: + totalRot = (nNWrot + nFWrot)*3 # going three-times through the entire wake + + tMax = dt_wanted*nPerRot*totalRot + + if show: + print(dt_wanted , ' dt') + print(int (nNWPanel ), ' nNWPanel ({} rotations)'.format(nNWrot)) + print(int (nFWPanel ), ' FarWakeLength ({} rotations)'.format(nFWrot)) + print(int (nFWPanelFree), ' FreeFarWakeLength ({} rotations)'.format(nFWrotFree)) + print(tMax , ' Tmax ({} rotations)'.format(totalRot)) + + return dt_wanted, tMax, nNWPanel, nFWPanel, nFWPanelFree class ModesElastoDyn(ExplicitComponent): """ From 9e9ea0b6f3a7d6eeaa96557ab7aca890b9c22320 Mon Sep 17 00:00:00 2001 From: pibo Date: Mon, 24 May 2021 22:02:32 -0600 Subject: [PATCH 2/5] isfile instead of exists --- weis/aeroelasticse/FAST_reader.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/weis/aeroelasticse/FAST_reader.py b/weis/aeroelasticse/FAST_reader.py index 538e91748..d28ef46ea 100644 --- a/weis/aeroelasticse/FAST_reader.py +++ b/weis/aeroelasticse/FAST_reader.py @@ -429,7 +429,7 @@ def execute(self): self.read_MoorDyn() bd_file = os.path.join(self.FAST_directory, self.fst_vt['Fst']['BDBldFile(1)']) - if os.path.exists(bd_file): + if os.path.isfile(bd_file): self.read_BeamDyn(bd_file) def read_MainInput(self): @@ -1125,7 +1125,7 @@ def read_AeroDyn15(self): self.read_AeroDyn15Polar() self.read_AeroDyn15Coord() olaf_filename = os.path.join(self.FAST_directory, self.fst_vt['AeroDyn15']['OLAFInputFileName']) - if os.path.exists(olaf_filename): + if os.path.isfile(olaf_filename): self.read_AeroDyn15OLAF(olaf_filename) def read_AeroDyn15Blade(self): From 9de575105017b1553336cd8d374c265e00da368f Mon Sep 17 00:00:00 2001 From: pibo Date: Thu, 27 May 2021 16:43:43 -0600 Subject: [PATCH 3/5] one more output from olaf params --- weis/aeroelasticse/openmdao_openfast.py | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/weis/aeroelasticse/openmdao_openfast.py b/weis/aeroelasticse/openmdao_openfast.py index 8b2161166..8b7afd0b8 100644 --- a/weis/aeroelasticse/openmdao_openfast.py +++ b/weis/aeroelasticse/openmdao_openfast.py @@ -628,7 +628,7 @@ def update_FAST_model(self, fst_vt, inputs, discrete_inputs): # Update OLAF if fst_vt['AeroDyn15']['WakeMod'] == 3: - _, nNWPanel, nFWPanel, nFWPanelFree = OLAFParams(fst_vt['ElastoDyn']['RotSpeed']) + _, _, nNWPanel, nFWPanel, nFWPanelFree = OLAFParams(fst_vt['ElastoDyn']['RotSpeed']) fst_vt['AeroDyn15']['OLAF']['nNWPanel'] = nNWPanel fst_vt['AeroDyn15']['OLAF']['nFWPanel'] = nFWPanel fst_vt['AeroDyn15']['OLAF']['nFWPanelFree'] = nFWPanelFree @@ -1154,7 +1154,7 @@ def DLC_creation_IEC(self, inputs, discrete_inputs, fst_vt, powercurve=False): # Set DT according to OLAF guidelines if fst_vt['AeroDyn15']['WakeMod'] == 3: - dt_wanted, _, _, _ = OLAFParams(inputs['Omega_init']) + dt_wanted, _, _, _, _ = OLAFParams(inputs['Omega_init']) iec.init_cond[("Fst","DT")] = {'U':inputs['U_init']} iec.init_cond[("Fst","DT")]['val'] = dt_wanted @@ -1273,7 +1273,7 @@ def DLC_creation_powercurve(self, inputs, discrete_inputs, fst_vt): # Set DT according to OLAF guidelines if fst_vt['AeroDyn15']['WakeMod'] == 3: - dt_wanted, _, _, _ = OLAFParams(inputs['Omega_init']) + dt_wanted, _, _, _, _ = OLAFParams(inputs['Omega_init']) case_inputs[("Fst","DT")] = {'vals':dt_wanted, 'group':1} case_list, case_name = CaseGen_General(case_inputs, self.FAST_runDirectory, self.FAST_namingOut + '_powercurve') From e112d52e396fe5bcc47cf65a12cb9bb40140f125 Mon Sep 17 00:00:00 2001 From: pibo Date: Tue, 1 Jun 2021 16:18:32 -0600 Subject: [PATCH 4/5] raise grid center point when running olaf --- weis/aeroelasticse/CaseGen_IEC.py | 5 +++++ weis/aeroelasticse/pyIECWind.py | 7 ++++--- 2 files changed, 9 insertions(+), 3 deletions(-) diff --git a/weis/aeroelasticse/CaseGen_IEC.py b/weis/aeroelasticse/CaseGen_IEC.py index d0ea4233b..e21371c16 100644 --- a/weis/aeroelasticse/CaseGen_IEC.py +++ b/weis/aeroelasticse/CaseGen_IEC.py @@ -134,6 +134,11 @@ def execute(self, case_inputs={}): iecwind.dir_change = self.transient_dir_change iecwind.shear_orient = self.transient_shear_orientation iecwind.z_hub = self.z_hub + # if running olaf, raise the center of the grid above hub height because the wake will expand + if case_inputs_i['AeroDyn15', 'WakeMod']['vals'][0] == 3: + iecwind.z_grid_center = iecwind.z_hub * 1.5 + else: + iecwind.z_grid_center = iecwind.z_hub iecwind.D = self.D iecwind.PLExp = alpha diff --git a/weis/aeroelasticse/pyIECWind.py b/weis/aeroelasticse/pyIECWind.py index 32d0b42c6..8e9cee878 100644 --- a/weis/aeroelasticse/pyIECWind.py +++ b/weis/aeroelasticse/pyIECWind.py @@ -432,6 +432,7 @@ def __init__(self): self.seed = np.random.uniform(1, 1e8) self.Turbulence_Class = 'B' # IEC Turbulance Class self.z_hub = 90. # wind turbine hub height (m) + self.z_grid_center = 90 # height of grid center from the ground (m) self.D = 126. # rotor diameter (m) self.PLExp = 0.2 self.AnalysisTime = 720. @@ -443,9 +444,9 @@ def setup(self): turbsim_vt.runtime_options.RandSeed1 = self.seed turbsim_vt.runtime_options.WrADTWR = False turbsim_vt.tmspecs.AnalysisTime = self.AnalysisTime - turbsim_vt.tmspecs.HubHt = self.z_hub # wind grid centered at hub height - turbsim_vt.tmspecs.GridHeight = (self.z_hub - 1.) * 2.0 # wind grid stops 1 meter above the ground - turbsim_vt.tmspecs.GridWidth = (self.z_hub - 1.) * 2.0 # squared wind grid + turbsim_vt.tmspecs.HubHt = self.z_grid_center # wind grid may be centered at hub height, but it may also go higher. RefHt stays at hub height + turbsim_vt.tmspecs.GridHeight = (self.z_grid_center - 1.) * 2.0 # wind grid stops 1 meter above the ground + turbsim_vt.tmspecs.GridWidth = (self.z_grid_center - 1.) * 2.0 # squared wind grid turbsim_vt.tmspecs.NumGrid_Z = 25 turbsim_vt.tmspecs.NumGrid_Y = 25 turbsim_vt.tmspecs.HFlowAng = 0.0 From 77ca4f7f19dc2ccdd01dce0fb068d0a0b1dff350 Mon Sep 17 00:00:00 2001 From: pibo Date: Tue, 1 Jun 2021 16:39:34 -0600 Subject: [PATCH 5/5] better flags --- weis/aeroelasticse/CaseGen_IEC.py | 7 ++----- weis/aeroelasticse/openmdao_openfast.py | 5 ++++- 2 files changed, 6 insertions(+), 6 deletions(-) diff --git a/weis/aeroelasticse/CaseGen_IEC.py b/weis/aeroelasticse/CaseGen_IEC.py index e21371c16..9f69ea456 100644 --- a/weis/aeroelasticse/CaseGen_IEC.py +++ b/weis/aeroelasticse/CaseGen_IEC.py @@ -53,6 +53,7 @@ def __init__(self): self.TStart = 30. self.uniqueSeeds = False self.uniqueWaveSeeds = False + self.grid_center_over_hh = 1. # Ratio between turbsim grid center and hub height self.debug_level = 2 self.parallel_windfile_gen = False @@ -134,11 +135,7 @@ def execute(self, case_inputs={}): iecwind.dir_change = self.transient_dir_change iecwind.shear_orient = self.transient_shear_orientation iecwind.z_hub = self.z_hub - # if running olaf, raise the center of the grid above hub height because the wake will expand - if case_inputs_i['AeroDyn15', 'WakeMod']['vals'][0] == 3: - iecwind.z_grid_center = iecwind.z_hub * 1.5 - else: - iecwind.z_grid_center = iecwind.z_hub + iecwind.z_grid_center = self.grid_center_over_hh * iecwind.z_hub iecwind.D = self.D iecwind.PLExp = alpha diff --git a/weis/aeroelasticse/openmdao_openfast.py b/weis/aeroelasticse/openmdao_openfast.py index 8b7afd0b8..40714294b 100644 --- a/weis/aeroelasticse/openmdao_openfast.py +++ b/weis/aeroelasticse/openmdao_openfast.py @@ -1152,11 +1152,14 @@ def DLC_creation_IEC(self, inputs, discrete_inputs, fst_vt, powercurve=False): iec.init_cond[("ElastoDyn","BlPitch2")] = iec.init_cond[("ElastoDyn","BlPitch1")] iec.init_cond[("ElastoDyn","BlPitch3")] = iec.init_cond[("ElastoDyn","BlPitch1")] - # Set DT according to OLAF guidelines + # If running OLAF... if fst_vt['AeroDyn15']['WakeMod'] == 3: + # Set DT according to OLAF guidelines dt_wanted, _, _, _, _ = OLAFParams(inputs['Omega_init']) iec.init_cond[("Fst","DT")] = {'U':inputs['U_init']} iec.init_cond[("Fst","DT")]['val'] = dt_wanted + # Raise the center of the grid 50% above hub height because the wake will expand + iec.grid_center_over_hh = 1.5 # Todo: need a way to handle Metocean conditions for Offshore # if offshore: