Results extraction and analysis from LS-DYNA sources#

This example provides an overview of the LS-DYNA results providers.

Note

This example requires DPF 6.1 (ansys-dpf-server-2023-2-pre0) or above. For more information, see Compatibility.

import matplotlib.pyplot as plt
from ansys.dpf import core as dpf
from ansys.dpf.core import examples

d3plot file results extraction#

Create the model and print its contents. This LS-DYNA d3plot file contains several individual results, each at different times. The d3plot file does not contain information related to Units. In this case, as the simulation was run through Mechanical, a file.actunits file is produced. If this file is supplemented in the data_sources, the units will be correctly fetched for all results in the file as well as for the mesh.

d3plot = examples.download_d3plot_beam()
ds = dpf.DataSources()
ds.set_result_file_path(d3plot[0], "d3plot")
ds.add_file_path(d3plot[3], "actunits")
model = dpf.Model(ds)
print(model)
DPF Model
------------------------------
Unknown analysis
Unit system: NMM: mm, ton, N, s, mV, mA, degC
Physics Type: Unknown
Available results:
     -  global_kinetic_energy: TimeFreq_steps Global Kinetic Energy
     -  global_internal_energy: TimeFreq_steps Global Internal Energy
     -  global_total_energy: TimeFreq_steps Global Total Energy
     -  global_velocity: TimeFreq_steps Global Velocity
     -  initial_coordinates: Nodal Initial Coordinates
     -  coordinates: Nodal Coordinates
     -  velocity: Nodal Velocity
     -  acceleration: Nodal Acceleration
     -  stress: Elemental Stress
     -  stress_von_mises: Elemental Stress Von Mises
     -  plastic_strain_eqv: Elemental Plastic Strain Eqv
     -  total_strain: Elemental Total Strain
     -  beam_axial_force: Elemental Beam Axial Force
     -  beam_s_shear_force: Elemental Beam S Shear Force
     -  beam_t_shear_force: Elemental Beam T Shear Force
     -  beam_s_bending_moment: Elemental Beam S Bending Moment
     -  beam_t_bending_moment: Elemental Beam T Bending Moment
     -  beam_torsional_moment: Elemental Beam Torsional Moment
     -  beam_axial_stress: Elemental Beam Axial Stress
     -  beam_rs_shear_stress: Elemental Beam Rs Shear Stress
     -  beam_tr_shear_stress: Elemental Beam Tr Shear Stress
     -  beam_axial_plastic_strain: Elemental Beam Axial Plastic Strain
     -  beam_axial_total_strain: Elemental Beam Axial Total Strain
     -  history_variablesihv__[1__10]: Elemental History Variables(ihv: [1, 10])
     -  displacement: Nodal Displacement
------------------------------
DPF  Meshed Region:
  1940 nodes
  2056 elements
  Unit: mm
  With solid (3D) elements, beam (1D) elements
------------------------------
DPF  Time/Freq Support:
  Number of sets: 12
Cumulative     Time (s)       LoadStep       Substep
1              0.000000       1              1
2              0.997460       2              1
3              1.997818       3              1
4              2.998237       4              1
5              3.998657       5              1
6              4.999076       6              1
7              5.999496       7              1
8              6.999915       8              1
9              7.997509       9              1
10             8.997929       10             1
11             9.998348       11             1
12             10.001174      12             1

The model has solid (3D) elements and beam (1D) elements. Some of the results only apply to one type of elements (such as the stress tensor for solids, or the axial force for beams, for example).

Let’s extract beam axial force for the last time step and plot it, deforming the mesh according to the displacement field at the same time step.

N = model.results.beam_axial_force(time_scoping=[12]).eval()
u = model.results.displacement(time_scoping=[12]).eval()

sargs = dict(title="N", fmt="%.2e", title_font_size=30, label_font_size=20)
N[0].plot(deform_by=u[0], scalar_bar_args=sargs)
00 lsdyna operators

The axial force has only been computed for the beam elements (the bottom frame), whereas the top sphere, which is comprised by solid elements, has only been deformed by the displacement field.

PyDPF also allows you to animate the results in a FieldsContainer. Thus, if all time steps are extracted, an animation can be produced.

N_all = model.results.beam_axial_force.on_all_time_freqs.eval()
u_all = model.results.displacement.on_all_time_freqs.eval()
N_all.animate(deform_by=u_all, save_as="falling_ball.gif", scalar_bar_args=sargs)
00 lsdyna operators

Some of the results are marked as global. They are not scoped over any mesh entity, but are global variables of the model.

K = model.results.global_kinetic_energy().eval()
print(K)
DPF  Fields Container
  with 1 field(s)
  defined on labels:

  with:
  - field 0 {} with TimeFreq_steps location, 1 components and 12 entities.

Energy plots over time (the sphere was released with some initial velocity).

U = model.results.global_internal_energy().eval()
H = model.results.global_total_energy().eval()

plt.plot(K.time_freq_support.time_frequencies.data, K[0].data, label="Kinetic")
plt.plot(U.time_freq_support.time_frequencies.data, U[0].data, label="Internal")
plt.plot(H.time_freq_support.time_frequencies.data, H[0].data, label="Total")
plt.xlabel("Time ({:s})".format(K.time_freq_support.time_frequencies.unit))
plt.ylabel("Energies ({:s})".format(K[0].unit))
plt.legend()
plt.show()
00 lsdyna operators

binout file results extraction#

Create the model and dprint its contents. This LS-DYNA binout file contains several branches (glstat, matsum and rcforc).

binout = examples.download_binout_matsum()
ds = dpf.DataSources()
ds.set_result_file_path(binout, "binout")
model = dpf.Model(ds)
print(model)
DPF Model
------------------------------
Unknown analysis
Unit system: Undefined
Physics Type: Unknown
Available results:
     -  global_time_step: TimeFreq_steps Global Time Step
     -  global_kinetic_energy: TimeFreq_steps Global Kinetic Energy
     -  global_internal_energy: TimeFreq_steps Global Internal Energy
     -  global_spring_damper_energy: TimeFreq_steps Global Spring Damper Energy
     -  global_joint_internal_energy: TimeFreq_steps Global Joint Internal Energy
     -  global_system_damping_energy: TimeFreq_steps Global System Damping Energy
     -  global_sliding_interface_energy: TimeFreq_steps Global Sliding Interface Energy
     -  global_external_work: TimeFreq_steps Global External Work
     -  global_eroded_kinetic_energy: TimeFreq_steps Global Eroded Kinetic Energy
     -  global_eroded_internal_energy: TimeFreq_steps Global Eroded Internal Energy
     -  global_total_energy: TimeFreq_steps Global Total Energy
     -  global_energy_ratio: TimeFreq_steps Global Energy Ratio
     -  global_energy_ratio_wo_eroded: TimeFreq_steps Global Energy Ratio Wo Eroded
     -  global_velocity: TimeFreq_steps Global Velocity
     -  global_added_mass: TimeFreq_steps Global Added Mass
     -  global_added_mass_percentage: TimeFreq_steps Global Added Mass Percentage
     -  part_internal_energy: TimeFreq_steps Part Internal Energy
     -  part_kinetic_energy: TimeFreq_steps Part Kinetic Energy
     -  part_eroded_internal_energy: TimeFreq_steps Part Eroded Internal Energy
     -  part_eroded_kinetic_energy: TimeFreq_steps Part Eroded Kinetic Energy
     -  part_added_mass: TimeFreq_steps Part Added Mass
     -  part_momentum: TimeFreq_steps Part Momentum
     -  part_rigid_body_velocity: TimeFreq_steps Part Rigid Body Velocity
     -  interface_contact_mass: TimeFreq_steps Interface Contact Mass
     -  interface_contact_force: TimeFreq_steps Interface Contact Force
     -  secforc_force_x: TimeFreq_steps Secforc Force X
     -  secforc_force_y: TimeFreq_steps Secforc Force Y
     -  secforc_force_z: TimeFreq_steps Secforc Force Z
     -  secforc_force_total: TimeFreq_steps Secforc Force Total
     -  secforc_moment_x: TimeFreq_steps Secforc Moment X
     -  secforc_moment_y: TimeFreq_steps Secforc Moment Y
     -  secforc_moment_z: TimeFreq_steps Secforc Moment Z
     -  secforc_moment_total: TimeFreq_steps Secforc Moment Total
     -  secforc_centroid_x: TimeFreq_steps Secforc Centroid X
     -  secforc_centroid_y: TimeFreq_steps Secforc Centroid Y
     -  secforc_centroid_z: TimeFreq_steps Secforc Centroid Z
     -  secforc_area: TimeFreq_steps Secforc Area
     -  nodout_acceleration_rx: TimeFreq_steps Nodout Acceleration Rx
     -  nodout_acceleration_ry: TimeFreq_steps Nodout Acceleration Ry
     -  nodout_acceleration_rz: TimeFreq_steps Nodout Acceleration Rz
     -  nodout_acceleration_x: TimeFreq_steps Nodout Acceleration X
     -  nodout_acceleration_y: TimeFreq_steps Nodout Acceleration Y
     -  nodout_acceleration_z: TimeFreq_steps Nodout Acceleration Z
     -  nodout_displacement_rx: TimeFreq_steps Nodout Displacement Rx
     -  nodout_displacement_ry: TimeFreq_steps Nodout Displacement Ry
     -  nodout_displacement_rz: TimeFreq_steps Nodout Displacement Rz
     -  nodout_displacement_x: TimeFreq_steps Nodout Displacement X
     -  nodout_displacement_y: TimeFreq_steps Nodout Displacement Y
     -  nodout_displacement_z: TimeFreq_steps Nodout Displacement Z
     -  nodout_velocity_rx: TimeFreq_steps Nodout Velocity Rx
     -  nodout_velocity_ry: TimeFreq_steps Nodout Velocity Ry
     -  nodout_velocity_rz: TimeFreq_steps Nodout Velocity Rz
     -  nodout_velocity_x: TimeFreq_steps Nodout Velocity X
     -  nodout_velocity_y: TimeFreq_steps Nodout Velocity Y
     -  nodout_velocity_z: TimeFreq_steps Nodout Velocity Z
     -  nodout_coordinate_x: TimeFreq_steps Nodout Coordinate X
     -  nodout_coordinate_y: TimeFreq_steps Nodout Coordinate Y
     -  nodout_coordinate_z: TimeFreq_steps Nodout Coordinate Z
     -  abstat_density: TimeFreq_steps Abstat Density
     -  abstat_input_mass_flow_rate_from_inflator: TimeFreq_steps Abstat Input Mass Flow Rate From Inflator
     -  abstat_leak_out_mass_flow_rate_from_vent_holes_and_fabric_porous_leakage: TimeFreq_steps Abstat Leak Out Mass Flow Rate From Vent Holes And Fabric Porous Leakage
     -  abstat_leak_out_mass_flow_rate_from_vent_holes: TimeFreq_steps Abstat Leak Out Mass Flow Rate From Vent Holes
     -  abstat_leak_out_mass_flow_rate_from_fabric_porous_leakage: TimeFreq_steps Abstat Leak Out Mass Flow Rate From Fabric Porous Leakage
     -  abstat_gas_temperature: TimeFreq_steps Abstat Gas Temperature
     -  abstat_internal_energy: TimeFreq_steps Abstat Internal Energy
     -  abstat_pressure: TimeFreq_steps Abstat Pressure
     -  abstat_reaction_force: TimeFreq_steps Abstat Reaction Force
     -  abstat_area_surface: TimeFreq_steps Abstat Area Surface
     -  abstat_area_total_mass: TimeFreq_steps Abstat Area Total Mass
     -  abstat_area_volume: TimeFreq_steps Abstat Area Volume
     -  sbtout_belt_force: TimeFreq_steps Sbtout Belt Force
     -  sbtout_belt_length: TimeFreq_steps Sbtout Belt Length
     -  sbtout_retractor_force: TimeFreq_steps Sbtout Retractor Force
     -  sbtout_retractor_pullout: TimeFreq_steps Sbtout Retractor Pullout
     -  sbtout_slipring: TimeFreq_steps Sbtout Slipring
     -  jntforc_force_x: TimeFreq_steps Jntforc Force X
     -  jntforc_force_y: TimeFreq_steps Jntforc Force Y
     -  jntforc_force_z: TimeFreq_steps Jntforc Force Z
     -  jntforc_moment_x: TimeFreq_steps Jntforc Moment X
     -  jntforc_moment_y: TimeFreq_steps Jntforc Moment Y
     -  jntforc_moment_z: TimeFreq_steps Jntforc Moment Z
     -  jntforc_resultant_force: TimeFreq_steps Jntforc Resultant Force
     -  jntforc_resultant_moment: TimeFreq_steps Jntforc Resultant Moment
     -  jntforc_joint_energy: TimeFreq_steps Jntforc Joint Energy
     -  bndout_rbdout_dircos_11: TimeFreq_steps Bndout Rbdout Dircos 11
     -  bndout_rbdout_dircos_12: TimeFreq_steps Bndout Rbdout Dircos 12
     -  bndout_rbdout_dircos_13: TimeFreq_steps Bndout Rbdout Dircos 13
     -  bndout_rbdout_dircos_21: TimeFreq_steps Bndout Rbdout Dircos 21
     -  bndout_rbdout_dircos_22: TimeFreq_steps Bndout Rbdout Dircos 22
     -  bndout_rbdout_dircos_23: TimeFreq_steps Bndout Rbdout Dircos 23
     -  bndout_rbdout_dircos_31: TimeFreq_steps Bndout Rbdout Dircos 31
     -  bndout_rbdout_dircos_32: TimeFreq_steps Bndout Rbdout Dircos 32
     -  bndout_rbdout_dircos_33: TimeFreq_steps Bndout Rbdout Dircos 33
     -  bndout_rbdout_global_acceleration: TimeFreq_steps Bndout Rbdout Global Acceleration
     -  bndout_rbdout_rotation_acceleration: TimeFreq_steps Bndout Rbdout Rotation Acceleration
     -  bndout_rbdout_dircos_rotation_displacement: TimeFreq_steps Bndout Rbdout Dircos Rotation Displacement
     -  bndout_rbdout_rotation_velocity: TimeFreq_steps Bndout Rbdout Rotation Velocity
     -  bndout_rbdout_global_velocity: TimeFreq_steps Bndout Rbdout Global Velocity
     -  bndout_rbdout_global_coordinate: TimeFreq_steps Bndout Rbdout Global Coordinate
     -  bndout_rbdout_location_acceleration: TimeFreq_steps Bndout Rbdout Location Acceleration
     -  bndout_rbdout_local_displacement: TimeFreq_steps Bndout Rbdout Local Displacement
     -  bndout_rbdout_local_rotation_acceleration: TimeFreq_steps Bndout Rbdout Local Rotation Acceleration
     -  bndout_rbdout_local_rotation_displacement: TimeFreq_steps Bndout Rbdout Local Rotation Displacement
     -  bndout_rbdout_local_velocity: TimeFreq_steps Bndout Rbdout Local Velocity
     -  bndout_rbdout_local_acceleration: TimeFreq_steps Bndout Rbdout Local Acceleration
     -  bndout_sleout_friction_energy: TimeFreq_steps Bndout Sleout Friction Energy
     -  bndout_sleout_master: TimeFreq_steps Bndout Sleout Master
     -  bndout_sleout_slave: TimeFreq_steps Bndout Sleout Slave
     -  bndout_sleout_slavemaster: TimeFreq_steps Bndout Sleout Slavemaster
     -  bndout_sleout_total_energy: TimeFreq_steps Bndout Sleout Total Energy
     -  bndout_sleout_total_friction: TimeFreq_steps Bndout Sleout Total Friction
     -  bndout_sleout_total_master: TimeFreq_steps Bndout Sleout Total Master
     -  bndout_sleout_total_slave: TimeFreq_steps Bndout Sleout Total Slave
------------------------------
DPF  Meshed Region:
  1 nodes
  Unit:
------------------------------
DPF  Time/Freq Support:
  Number of sets: 2002
Cumulative     Frequency ()   LoadStep       Substep
1              0.000000       1              1
2              0.008897       2              1
3              0.097868       3              1
4              0.106765       4              1
5              0.195737       5              1
6              0.204634       6              1
7              0.293605       7              1
8              0.302502       8              1
9              0.391474       9              1
10             0.400371       10             1
11             0.498239       11             1
12             0.507136       12             1
13             0.596107       13             1
14             0.605004       14             1
15             0.693976       15             1
16             0.702873       16             1
17             0.791844       17             1
18             0.800741       18             1
19             0.898610       19             1
20             0.907507       20             1
21             0.996478       21             1
22             1.005375       22             1
23             1.094346       23             1
24             1.103243       24             1
25             1.192215       25             1
26             1.201112       26             1
27             1.298980       27             1
28             1.307877       28             1
29             1.396849       29             1
30             1.405746       30             1
31             1.494717       31             1
32             1.503614       32             1
33             1.592585       33             1
34             1.601483       34             1
35             1.699351       35             1
36             1.708248       36             1
37             1.797219       37             1
38             1.806116       38             1
39             1.895088       39             1
40             1.903985       40             1
41             1.992956       41             1
42             2.001853       42             1
43             2.099721       43             1
44             2.108618       44             1
45             2.197590       45             1
46             2.206487       46             1
47             2.295462       47             1
48             2.304362       48             1
49             2.393364       49             1
50             2.402261       50             1
51             2.491304       51             1
52             2.500201       52             1
53             2.598073       53             1
54             2.606983       54             1
55             2.696023       55             1
56             2.704934       56             1
57             2.794000       57             1
58             2.802911       58             1
59             2.891951       59             1
60             2.900861       60             1
61             2.998879       61             1
62             3.007790       62             1
63             3.096866       63             1
64             3.105777       64             1
65             3.194883       65             1
66             3.203794       66             1
67             3.292901       67             1
68             3.301811       68             1
69             3.399788       69             1
70             3.408699       70             1
71             3.497805       71             1
72             3.506716       72             1
73             3.595823       73             1
74             3.604733       74             1
75             3.693840       75             1
76             3.702751       76             1
77             3.791857       77             1
78             3.800768       78             1
79             3.898755       79             1
80             3.907665       80             1
81             3.996773       81             1
82             4.005683       82             1
83             4.094749       83             1
84             4.103660       84             1
85             4.192767       85             1
86             4.201678       86             1
87             4.299624       87             1
88             4.308535       88             1
89             4.397628       89             1
90             4.406539       90             1
91             4.495646       91             1
92             4.504557       92             1
93             4.593663       93             1
94             4.602574       94             1
95             4.691681       95             1
96             4.700592       96             1
97             4.798609       97             1
98             4.807520       98             1
99             4.896627       99             1
100            4.905538       100            1
101            4.994594       101            1
102            5.003504       102            1
103            5.092611       103            1
104            5.101522       104            1
105            5.199547       105            1
106            5.208459       106            1
107            5.297572       107            1
108            5.306483       108            1
109            5.395590       109            1
110            5.404501       110            1
111            5.493581       111            1
112            5.502492       112            1
113            5.591599       113            1
114            5.600510       114            1
115            5.698528       115            1
116            5.707439       116            1
117            5.796547       117            1
118            5.805458       118            1
119            5.894565       119            1
120            5.903475       120            1
121            5.992583       121            1
122            6.001493       122            1
123            6.099511       123            1
124            6.108422       124            1
125            6.197529       125            1
126            6.206440       126            1
127            6.295547       127            1
128            6.304458       128            1
129            6.393566       129            1
130            6.402476       130            1
131            6.491584       131            1
132            6.500494       132            1
133            6.598444       133            1
134            6.607341       134            1
135            6.696421       135            1
136            6.705332       136            1
137            6.794439       137            1
138            6.803350       138            1
139            6.892458       139            1
140            6.901369       140            1
141            6.999387       141            1
142            7.008299       142            1
143            7.097407       143            1
144            7.106318       144            1
145            7.195426       145            1
146            7.204337       146            1
147            7.293445       147            1
148            7.302355       148            1
149            7.391466       149            1
150            7.400377       150            1
151            7.498406       151            1
152            7.507318       152            1
153            7.596435       153            1
154            7.605347       154            1
155            7.694464       155            1
156            7.703375       156            1
157            7.792492       157            1
158            7.801404       158            1
159            7.899434       159            1
160            7.908345       160            1
161            7.997462       161            1
162            8.006374       162            1
163            8.095469       163            1
164            8.104369       164            1
165            8.193441       165            1
166            8.202342       166            1
167            8.291460       167            1
168            8.300357       168            1
169            8.398328       169            1
170            8.407240       170            1
171            8.496346       171            1
172            8.505246       172            1
173            8.594317       173            1
174            8.603230       174            1
175            8.692324       175            1
176            8.701236       176            1
177            8.799255       177            1
178            8.808166       178            1
179            8.897273       179            1
180            8.906184       180            1
181            8.995292       181            1
182            9.004202       182            1
183            9.093309       183            1
184            9.102221       184            1
185            9.191328       185            1
186            9.200238       186            1
187            9.298257       187            1
188            9.307167       188            1
189            9.396275       189            1
190            9.405186       190            1
191            9.494293       191            1
192            9.503203       192            1
193            9.592311       193            1
194            9.601222       194            1
195            9.699240       195            1
196            9.708151       196            1
197            9.797258       197            1
198            9.806169       198            1
199            9.895276       199            1
200            9.904187       200            1
201            9.993294       201            1
202            10.002205      202            1
203            10.091312      203            1
204            10.100224      204            1
205            10.198241      205            1
206            10.207151      206            1
207            10.296259      207            1
208            10.305170      208            1
209            10.394279      209            1
210            10.403189      210            1
211            10.492304      211            1
212            10.501216      212            1
213            10.599245      213            1
214            10.608157      214            1
215            10.697274      215            1
216            10.706186      216            1
217            10.795304      217            1
218            10.804215      218            1
219            10.893333      219            1
220            10.902245      220            1
221            10.991363      221            1
222            11.000275      222            1
223            11.098304      223            1
224            11.107215      224            1
225            11.196333      225            1
226            11.205245      226            1
227            11.294362      227            1
228            11.303274      228            1
229            11.392391      229            1
230            11.401303      230            1
231            11.499333      231            1
232            11.508245      232            1
233            11.597363      233            1
234            11.606274      234            1
235            11.695392      235            1
236            11.704304      236            1
237            11.793422      237            1
238            11.802333      238            1
239            11.891449      239            1
240            11.900359      240            1
241            11.998384      241            1
242            12.007294      242            1
243            12.096347      243            1
244            12.105259      244            1
245            12.194366      245            1
246            12.203276      246            1
247            12.292382      247            1
248            12.301292      248            1
249            12.399310      249            1
250            12.408221      250            1
251            12.497328      251            1
252            12.506238      252            1
253            12.595345      253            1
254            12.604255      254            1
255            12.693352      255            1
256            12.702262      256            1
257            12.791368      257            1
258            12.800280      258            1
259            12.898229      259            1
260            12.907125      260            1
261            12.996218      261            1
262            13.005129      262            1
263            13.094235      263            1
264            13.103146      264            1
265            13.192252      265            1
266            13.201163      266            1
267            13.299180      267            1
268            13.308091      268            1
269            13.397198      269            1
270            13.406108      270            1
271            13.495215      271            1
272            13.504128      272            1
273            13.593246      273            1
274            13.602159      274            1
275            13.691278      275            1
276            13.700190      276            1
277            13.798222      277            1
278            13.807134      278            1
279            13.896254      279            1
280            13.905166      280            1
281            13.994215      281            1
282            14.003116      282            1
283            14.092141      283            1
284            14.101042      284            1
285            14.198971      285            1
286            14.207870      286            1
287            14.296902      287            1
288            14.305802      288            1
289            14.394816      289            1
290            14.403717      290            1
291            14.492726      291            1
292            14.501623      292            1
293            14.599604      293            1
294            14.608504      294            1
295            14.697553      295            1
296            14.706450      296            1
297            14.795456      297            1
298            14.804366      298            1
299            14.893403      299            1
300            14.902300      300            1
301            14.991365      301            1
302            15.000277      302            1
303            15.098253      303            1
304            15.107150      304            1
305            15.196204      305            1
306            15.205114      306            1
307            15.294127      307            1
308            15.303024      308            1
309            15.392003      309            1
310            15.400900      310            1
311            15.498770      311            1
312            15.507667      312            1
313            15.596642      313            1
314            15.605539      314            1
315            15.694518      315            1
316            15.703419      316            1
317            15.792427      317            1
318            15.801324      318            1
319            15.899206      319            1
320            15.908117      320            1
321            15.997162      321            1
322            16.006063      322            1
323            16.095060      323            1
324            16.103958      324            1
325            16.192940      325            1
326            16.201836      326            1
327            16.299715      327            1
328            16.308613      328            1
329            16.397587      329            1
330            16.406485      330            1
331            16.495459      331            1
332            16.504356      332            1
333            16.593330      333            1
334            16.602228      334            1
335            16.691200      335            1
336            16.700098      336            1
337            16.797968      337            1
338            16.806864      338            1
339            16.895941      339            1
340            16.904840      340            1
341            16.993937      341            1
342            17.002836      342            1
343            17.091888      343            1
344            17.100786      344            1
345            17.198776      345            1
346            17.207687      346            1
347            17.296793      347            1
348            17.305695      348            1
349            17.394678      349            1
350            17.403576      350            1
351            17.492554      351            1
352            17.501465      352            1
353            17.599354      353            1
354            17.608252      354            1
355            17.697252      355            1
356            17.706154      356            1
357            17.795135      357            1
358            17.804033      358            1
359            17.893005      359            1
360            17.901901      360            1
361            17.999775      361            1
362            18.008673      362            1
363            18.097645      363            1
364            18.106543      364            1
365            18.195517      365            1
366            18.204412      366            1
367            18.293386      367            1
368            18.302284      368            1
369            18.391256      369            1
370            18.400154      370            1
371            18.498024      371            1
372            18.506922      372            1
373            18.595894      373            1
374            18.604790      374            1
375            18.693762      375            1
376            18.702660      376            1
377            18.791632      377            1
378            18.800529      378            1
379            18.898399      379            1
380            18.907295      380            1
381            18.996269      381            1
382            19.005167      382            1
383            19.094139      383            1
384            19.103037      384            1
385            19.192009      385            1
386            19.200907      386            1
387            19.298775      387            1
388            19.307673      388            1
389            19.396645      389            1
390            19.405542      390            1
391            19.494514      391            1
392            19.503412      392            1
393            19.592386      393            1
394            19.601282      394            1
395            19.699152      395            1
396            19.708050      396            1
397            19.797024      397            1
398            19.805920      398            1
399            19.894894      399            1
400            19.903791      400            1
401            19.992764      401            1
402            20.001661      402            1
403            20.099531      403            1
404            20.108429      404            1
405            20.197401      405            1
406            20.206299      406            1
407            20.295271      407            1
408            20.304169      408            1
409            20.393141      409            1
410            20.402037      410            1
411            20.499907      411            1
412            20.508804      412            1
413            20.597776      413            1
414            20.606674      414            1
415            20.695646      415            1
416            20.704544      416            1
417            20.793518      417            1
418            20.802414      418            1
419            20.891388      419            1
420            20.900286      420            1
421            20.998158      421            1
422            21.007055      422            1
423            21.096027      423            1
424            21.104925      424            1
425            21.193899      425            1
426            21.202797      426            1
427            21.291771      427            1
428            21.300667      428            1
429            21.398764      429            1
430            21.400885      430            1
431            21.498735      431            1
432            21.500870      432            1
433            21.599386      433            1
434            21.601534      434            1
435            21.698536      435            1
436            21.700699      436            1
437            21.798298      437            1
438            21.800472      438            1
439            21.898647      439            1
440            21.900835      440            1
441            21.999575      441            1
442            22.001776      442            1
443            22.098860      443            1
444            22.101072      444            1
445            22.198702      445            1
446            22.200928      446            1
447            22.299099      447            1
448            22.301336      448            1
449            22.400043      449            1
450            22.402294      450            1
451            22.499275      451            1
452            22.501537      452            1
453            22.599031      453            1
454            22.601305      454            1
455            22.699329      455            1
456            22.701616      456            1
457            22.797865      457            1
458            22.800161      458            1
459            22.899208      459            1
460            22.901518      460            1
461            22.998747      461            1
462            23.001068      462            1
463            23.098791      463            1
464            23.101124      464            1
465            23.199358      465            1
466            23.201704      466            1
467            23.298098      467            1
468            23.300455      468            1
469            23.399715      469            1
470            23.402084      470            1
471            23.499483      471            1
472            23.501863      472            1
473            23.599775      473            1
474            23.602169      474            1
475            23.698185      475            1
476            23.700592      476            1
477            23.799513      477            1
478            23.801933      478            1
479            23.898933      479            1
480            23.901365      480            1
481            23.998865      481            1
482            24.001308      482            1
483            24.099319      483            1
484            24.101776      484            1
485            24.197836      485            1
486            24.200304      486            1
487            24.299360      487            1
488            24.301844      488            1
489            24.398956      489            1
490            24.401453      490            1
491            24.499113      491            1
492            24.501623      492            1
493            24.599850      493            1
494            24.602375      494            1
495            24.698679      495            1
496            24.701221      496            1
497            24.798164      497            1
498            24.800724      498            1
499            24.898348      499            1
500            24.900925      500            1
501            24.999241      501            1
502            25.001837      502            1
503            25.098253      503            1
504            25.100870      504            1
505            25.198011      505            1
506            25.200647      506            1
507            25.298561      507            1
508            25.301218      508            1
509            25.399933      509            1
510            25.402613      510            1
511            25.499472      511            1
512            25.502174      512            1
513            25.599909      513            1
514            25.602638      514            1
515            25.698580      515            1
516            25.701336      516            1
517            25.798267      517            1
518            25.801052      518            1
519            25.899038      519            1
520            25.901854      520            1
521            25.998114      521            1
522            26.000961      522            1
523            26.098392      523            1
524            26.101276      524            1
525            26.199966      525            1
526            26.202888      526            1
527            26.299913      527            1
528            26.302872      528            1
529            26.398115      529            1
530            26.401110      530            1
531            26.497496      531            1
532            26.500526      532            1
533            26.598053      533            1
534            26.601118      534            1
535            26.699776      535            1
536            26.702877      536            1
537            26.799532      537            1
538            26.802668      538            1
539            26.897236      539            1
540            26.900406      540            1
541            26.999231      541            1
542            27.002438      542            1
543            27.099184      543            1
544            27.102427      544            1
545            27.197046      545            1
546            27.200327      546            1
547            27.299398      547            1
548            27.302721      548            1
549            27.399677      549            1
550            27.403040      550            1
551            27.497784      551            1
552            27.501188      552            1
553            27.597071      553            1
554            27.600515      554            1
555            27.697582      555            1
556            27.701069      556            1
557            27.799368      557            1
558            27.802900      558            1
559            27.898903      559            1
560            27.902481      560            1
561            27.999710      561            1
562            28.003336      562            1
563            28.098152      563            1
564            28.101822      564            1
565            28.197824      565            1
566            28.201540      566            1
567            28.298733      567            1
568            28.302494      568            1
569            28.397083      569            1
570            28.400888      570            1
571            28.496624      571            1
572            28.500477      572            1
573            28.597397      573            1
574            28.601299      574            1
575            28.699442      575            1
576            28.703392      576            1
577            28.798777      577            1
578            28.802776      578            1
579            28.899324      579            1
580            28.903372      580            1
581            28.997013      581            1
582            29.001108      582            1
583            29.100006      583            1
584            29.104153      584            1
585            29.195889      585            1
586            29.200081      586            1
587            29.297098      587            1
588            29.301342      588            1
589            29.399496      589            1
590            29.403788      590            1
591            29.498739      591            1
592            29.503080      592            1
593            29.599081      593            1
594            29.603468      594            1
595            29.696087      595            1
596            29.700520      596            1
597            29.798573      597            1
598            29.803053      598            1
599            29.897617      599            1
600            29.902142      600            1
601            29.997641      601            1
602            30.002209      602            1
603            30.098639      603            1
604            30.103252      604            1
605            30.195957      605            1
606            30.200613      606            1
607            30.298866      607            1
608            30.303568      608            1
609            30.398010      609            1
610            30.402754      610            1
611            30.498035      611            1
612            30.502819      612            1
613            30.598740      613            1
614            30.603500      614            1
615            30.697105      615            1
616            30.701708      616            1
617            30.796635      617            1
618            30.801075      618            1
619            30.896976      619            1
620            30.901255      620            1
621            30.996698      621            1
622            31.001333      622            1
623            31.095596      623            1
624            31.100388      624            1
625            31.197756      625            1
626            31.202702      626            1
627            31.298031      627            1
628            31.303120      628            1
629            31.395048      629            1
630            31.400173      630            1
631            31.497917      631            1
632            31.503080      632            1
633            31.596348      633            1
634            31.601547      634            1
635            31.695446      635            1
636            31.700680      636            1
637            31.795191      637            1
638            31.800457      638            1
639            31.895548      639            1
640            31.900846      640            1
641            31.996508      641            1
642            32.001839      642            1
643            32.098083      643            1
644            32.103447      644            1
645            32.194904      645            1
646            32.200298      646            1
647            32.297726      647            1
648            32.303154      648            1
649            32.395702      649            1
650            32.401161      650            1
651            32.499683      651            1
652            32.505169      652            1
653            32.598660      653            1
654            32.604172      654            1
655            32.698044      655            1
656            32.703575      656            1
657            32.797794      657            1
658            32.803349      658            1
659            32.897907      659            1
660            32.903481      660            1
661            32.998405      661            1
662            33.003998      662            1
663            33.099312      663            1
664            33.104931      664            1
665            33.195011      665            1
666            33.200653      666            1
667            33.296768      667            1
668            33.302433      668            1
669            33.398956      669            1
670            33.404648      670            1
671            33.495861      671            1
672            33.501572      672            1
673            33.598877      673            1
674            33.604610      674            1
675            33.696568      675            1
676            33.702324      676            1
677            33.794651      677            1
678            33.800434      678            1
679            33.898952      679            1
680            33.904762      680            1
681            33.997894      681            1
682            34.003731      682            1
683            34.097282      683            1
684            34.103146      684            1
685            34.197136      685            1
686            34.203026      686            1
687            34.297470      687            1
688            34.303387      688            1
689            34.398304      689            1
690            34.404251      690            1
691            34.499653      691            1
692            34.505630      692            1
693            34.595516      693            1
694            34.601524      694            1
695            34.697891      695            1
696            34.703930      696            1
697            34.794743      697            1
698            34.800812      698            1
699            34.898186      699            1
700            34.904289      700            1
701            34.996063      701            1
702            35.002201      702            1
703            35.094463      703            1
704            35.100632      704            1
705            35.199596      705            1
706            35.205799      706            1
707            35.299114      707            1
708            35.305351      708            1
709            35.399197      709            1
710            35.405468      710            1
711            35.499851      711            1
712            35.506161      712            1
713            35.594730      713            1
714            35.601074      714            1
715            35.695652      715            1
716            35.701862      716            1
717            35.799515      717            1
718            35.805504      718            1
719            35.899506      719            1
720            35.905266      720            1
721            35.995525      721            1
722            36.001049      722            1
723            36.098106      723            1
724            36.103367      724            1
725            36.195683      725            1
726            36.200676      726            1
727            36.297588      727            1
728            36.302284      728            1
729            36.397541      729            1
730            36.401917      730            1
731            36.496498      731            1
732            36.501205      732            1
733            36.598366      733            1
734            36.603374      734            1
735            36.698254      735            1
736            36.700443      736            1
737            36.799324      737            1
738            36.801632      738            1
739            36.898586      739            1
740            36.901009      740            1
741            36.997723      741            1
742            37.000256      742            1
743            37.098717      743            1
744            37.101364      744            1
745            37.198433      745            1
746            37.201183      746            1
747            37.299068      747            1
748            37.301918      748            1
749            37.397396      749            1
750            37.400333      750            1
751            37.498909      751            1
752            37.501942      752            1
753            37.597363      753            1
754            37.600483      754            1
755            37.698647      755            1
756            37.701859      756            1
757            37.799484      757            1
758            37.802784      758            1
759            37.899643      759            1
760            37.903023      760            1
761            37.998882      761            1
762            38.002346      762            1
763            38.096966      763            1
764            38.100510      764            1
765            38.197273      765            1
766            38.200897      766            1
767            38.299793      767            1
768            38.303497      768            1
769            38.396961      769            1
770            38.400734      770            1
771            38.499851      771            1
772            38.503700      772            1
773            38.596916      773            1
774            38.600834      774            1
775            38.699688      775            1
776            38.703678      776            1
777            38.796211      777            1
778            38.800270      778            1
779            38.898472      779            1
780            38.902599      780            1
781            38.998302      781            1
782            39.002499      782            1
783            39.099781      783            1
784            39.104046      784            1
785            39.198582      785            1
786            39.202911      786            1
787            39.298496      787            1
788            39.302895      788            1
789            39.395855      789            1
790            39.400311      790            1
791            39.499107      791            1
792            39.503632      792            1
793            39.599281      793            1
794            39.603867      794            1
795            39.696175      795            1
796            39.700817      796            1
797            39.798985      797            1
798            39.803692      798            1
799            39.898350      799            1
800            39.903111      800            1
801            39.998882      801            1
802            40.003696      802            1
803            40.095692      803            1
804            40.100555      804            1
805            40.198425      805            1
806            40.203342      806            1
807            40.297260      807            1
808            40.302223      808            1
809            40.396984      809            1
810            40.401989      810            1
811            40.497517      811            1
812            40.502563      812            1
813            40.598789      813            1
814            40.603870      814            1
815            40.695621      815            1
816            40.700733      816            1
817            40.798138      817            1
818            40.803280      818            1
819            40.896027      819            1
820            40.901188      820            1
821            40.999531      821            1
822            41.004719      822            1
823            41.098259      823            1
824            41.103466      824            1
825            41.197323      825            1
826            41.202545      826            1
827            41.296677      827            1
828            41.301914      828            1
829            41.396275      829            1
830            41.401524      830            1
831            41.496067      831            1
832            41.501324      832            1
833            41.596027      833            1
834            41.601295      834            1
835            41.696152      835            1
836            41.701424      836            1
837            41.796444      837            1
838            41.801731      838            1
839            41.896938      839            1
840            41.902233      840            1
841            41.997654      841            1
842            42.002960      842            1
843            42.097473      843            1
844            42.102619      844            1
845            42.198185      845            1
846            42.203098      846            1
847            42.298985      847            1
848            42.303658      848            1
849            42.399212      849            1
850            42.403641      850            1
851            42.496552      851            1
852            42.501148      852            1
853            42.595505      853            1
854            42.600342      854            1
855            42.699448      855            1
856            42.704525      856            1
857            42.797859      857            1
858            42.803154      858            1
859            42.894871      859            1
860            42.900318      860            1
861            42.998627      861            1
862            43.004097      862            1
863            43.097332      863            1
864            43.102829      864            1
865            43.196476      865            1
866            43.201996      866            1
867            43.296085      867            1
868            43.301632      868            1
869            43.396194      869            1
870            43.401772      870            1
871            43.496834      871            1
872            43.502441      872            1
873            43.598022      873            1
874            43.603661      874            1
875            43.699772      875            1
876            43.705441      876            1
877            43.796383      877            1
878            43.802082      878            1
879            43.899239      879            1
880            43.904968      880            1
881            43.996918      881            1
882            44.002682      882            1
883            44.095146      883            1
884            44.100941      884            1
885            44.199768      885            1
886            44.205601      886            1
887            44.299168      887            1
888            44.305035      888            1
889            44.399162      889            1
890            44.405064      890            1
891            44.499771      891            1
892            44.505707      892            1
893            44.595058      893            1
894            44.601032      894            1
895            44.696953      895            1
896            44.702965      896            1
897            44.799515      897            1
898            44.805569      898            1
899            44.896626      899            1
900            44.902714      900            1
901            44.994236      901            1
902            45.000305      902            1
903            45.095139      903            1
904            45.100925      904            1
905            45.196796      905            1
906            45.202286      906            1
907            45.298306      907            1
908            45.303482      908            1
909            45.398643      909            1
910            45.403484      910            1
911            45.496727      911            1
912            45.501213      912            1
913            45.599728      913            1
914            45.604412      914            1
915            45.696220      915            1
916            45.701199      916            1
917            45.798355      917            1
918            45.803604      918            1
919            45.894894      919            1
920            45.900387      920            1
921            45.995911      921            1
922            46.001659      922            1
923            46.095600      923            1
924            46.101597      924            1
925            46.199570      925            1
926            46.205822      926            1
927            46.294743      927            1
928            46.301140      928            1
929            46.397125      929            1
930            46.403522      930            1
931            46.499519      931            1
932            46.505920      932            1
933            46.595528      933            1
934            46.601925      934            1
935            46.697941      935            1
936            46.704342      936            1
937            46.793964      937            1
938            46.800365      938            1
939            46.898586      939            1
940            46.900734      940            1
941            46.999065      941            1
942            47.001286      942            1
943            47.098663      943            1
944            47.100971      944            1
945            47.199722      945            1
946            47.202118      946            1
947            47.299709      947            1
948            47.302193      948            1
949            47.398182      949            1
950            47.400749      950            1
951            47.497311      951            1
952            47.499962      952            1
953            47.599579      953            1
954            47.602310      954            1
955            47.699291      955            1
956            47.702103      956            1
957            47.798820      957            1
958            47.801701      958            1
959            47.897900      959            1
960            47.900848      960            1
961            47.999321      961            1
962            48.002342      962            1
963            48.096905      963            1
964            48.099987      964            1
965            48.199684      965            1
966            48.202835      966            1
967            48.298237      967            1
968            48.301449      968            1
969            48.398777      969            1
970            48.402054      970            1
971            48.498035      971            1
972            48.501381      972            1
973            48.599335      973            1
974            48.602749      974            1
975            48.699295      975            1
976            48.702782      976            1
977            48.797794      977            1
978            48.801350      978            1
979            48.898296      979            1
980            48.901924      980            1
981            48.997124      981            1
982            49.000820      982            1
983            49.097820      983            1
984            49.101585      984            1
985            49.196495      985            1
986            49.200325      986            1
987            49.296780      987            1
988            49.300667      988            1
989            49.398594      989            1
990            49.402538      990            1
991            49.497856      991            1
992            49.501854      992            1
993            49.598370      993            1
994            49.602413      994            1
995            49.695961      995            1
996            49.700050      996            1
997            49.798687      997            1
998            49.802818      998            1
999            49.898251      999            1
1000           49.902416      1000           1
1001           49.998661      1001           1
1002           50.002861      1002           1
1003           50.095600      1003           1
1004           50.099827      1004           1
1005           50.197403      1005           1
1006           50.201656      1006           1
1007           50.299770      1007           1
1008           50.304047      1008           1
1009           50.398323      1009           1
1010           50.402618      1010           1
1011           50.497257      1011           1
1012           50.501564      1012           1
1013           50.596539      1013           1
1014           50.600864      1014           1
1015           50.696148      1015           1
1016           50.700485      1016           1
1017           50.796036      1017           1
1018           50.800385      1018           1
1019           50.896156      1019           1
1020           50.900513      1020           1
1021           50.996483      1021           1
1022           51.000851      1022           1
1023           51.097012      1023           1
1024           51.101387      1024           1
1025           51.197742      1025           1
1026           51.202129      1026           1
1027           51.298710      1027           1
1028           51.303104      1028           1
1029           51.395569      1029           1
1030           51.399979      1030           1
1031           51.497196      1031           1
1032           51.501625      1032           1
1033           51.599281      1033           1
1034           51.603733      1034           1
1035           51.697464      1035           1
1036           51.701942      1036           1
1037           51.796288      1037           1
1038           51.800797      1038           1
1039           51.895859      1039           1
1040           51.900402      1040           1
1041           51.996254      1041           1
1042           52.000835      1042           1
1043           52.097549      1043           1
1044           52.102177      1044           1
1045           52.195145      1045           1
1046           52.199814      1046           1
1047           52.298397      1047           1
1048           52.303116      1048           1
1049           52.397980      1049           1
1050           52.402748      1050           1
1051           52.498501      1051           1
1052           52.503304      1052           1
1053           52.599464      1053           1
1054           52.604271      1054           1
1055           52.695415      1055           1
1056           52.700203      1056           1
1057           52.795685      1057           1
1058           52.800449      1058           1
1059           52.899094      1059           1
1060           52.903641      1060           1
1061           52.997910      1061           1
1062           53.002590      1062           1
1063           53.097610      1063           1
1064           53.102375      1064           1
1065           53.197834      1065           1
1066           53.202618      1066           1
1067           53.298611      1067           1
1068           53.303429      1068           1
1069           53.395298      1069           1
1070           53.400158      1070           1
1071           53.497799      1071           1
1072           53.502708      1072           1
1073           53.596516      1073           1
1074           53.601482      1074           1
1075           53.696472      1075           1
1076           53.701508      1076           1
1077           53.797825      1077           1
1078           53.802933      1078           1
1079           53.895523      1079           1
1080           53.900707      1080           1
1081           53.994690      1081           1
1082           53.999950      1082           1
1083           54.095367      1083           1
1084           54.100708      1084           1
1085           54.197590      1085           1
1086           54.203014      1086           1
1087           54.295849      1087           1
1088           54.301350      1088           1
1089           54.395443      1089           1
1090           54.401016      1090           1
1091           54.496315      1091           1
1092           54.501953      1092           1
1093           54.598404      1093           1
1094           54.604111      1094           1
1095           54.695896      1095           1
1096           54.701664      1096           1
1097           54.794392      1097           1
1098           54.800217      1098           1
1099           54.899723      1099           1
1100           54.905605      1100           1
1101           54.994213      1101           1
1102           55.000149      1102           1
1103           55.099247      1103           1
1104           55.104866      1104           1
1105           55.196472      1105           1
1106           55.201630      1106           1
1107           55.299496      1107           1
1108           55.304134      1108           1
1109           55.399052      1109           1
1110           55.403740      1110           1
1111           55.495045      1111           1
1112           55.499928      1112           1
1113           55.598309      1113           1
1114           55.603260      1114           1
1115           55.697903      1115           1
1116           55.702923      1116           1
1117           55.799137      1117           1
1118           55.804249      1118           1
1119           55.897148      1119           1
1120           55.902359      1120           1
1121           55.997131      1121           1
1122           56.002453      1122           1
1123           56.099277      1123           1
1124           56.104717      1124           1
1125           56.198238      1125           1
1126           56.203804      1126           1
1127           56.299503      1127           1
1128           56.305199      1128           1
1129           56.397381      1129           1
1130           56.403206      1130           1
1131           56.497566      1131           1
1132           56.503532      1132           1
1133           56.594013      1133           1
1134           56.600109      1134           1
1135           56.698742      1135           1
1136           56.704941      1136           1
1137           56.797749      1137           1
1138           56.803925      1138           1
1139           56.896503      1139           1
1140           56.902672      1140           1
1141           56.995148      1141           1
1142           57.001312      1142           1
1143           57.093777      1143           1
1144           57.099941      1144           1
1145           57.198627      1145           1
1146           57.204800      1146           1
1147           57.297436      1147           1
1148           57.303619      1148           1
1149           57.396423      1149           1
1150           57.402618      1150           1
1151           57.495628      1151           1
1152           57.501835      1152           1
1153           57.595066      1153           1
1154           57.601288      1154           1
1155           57.694752      1155           1
1156           57.700989      1156           1
1157           57.794678      1157           1
1158           57.800930      1158           1
1159           57.894829      1159           1
1160           57.901096      1160           1
1161           57.995201      1161           1
1162           58.001484      1162           1
1163           58.095810      1163           1
1164           58.102104      1164           1
1165           58.196674      1165           1
1166           58.202988      1166           1
1167           58.297829      1167           1
1168           58.304161      1168           1
1169           58.399292      1169           1
1170           58.405643      1170           1
1171           58.497662      1171           1
1172           58.505173      1172           1
1173           58.596508      1173           1
1174           58.602882      1174           1
1175           58.693451      1175           1
1176           58.701885      1176           1
1177           58.792042      1177           1
1178           58.800476      1178           1
1179           58.893253      1179           1
1180           58.901688      1180           1
1181           58.993610      1181           1
1182           59.001785      1182           1
1183           59.097088      1183           1
1184           59.104786      1184           1
1185           59.195107      1185           1
1186           59.200550      1186           1
1187           59.294846      1187           1
1188           59.299881      1188           1
1189           59.398270      1189           1
1190           59.402866      1190           1
1191           59.496750      1191           1
1192           59.501377      1192           1
1193           59.595139      1193           1
1194           59.601757      1194           1
1195           59.698166      1195           1
1196           59.701157      1196           1
1197           59.798031      1197           1
1198           59.801189      1198           1
1199           59.898396      1199           1
1200           59.901714      1200           1
1201           59.996727      1201           1
1202           60.000195      1202           1
1203           60.099335      1203           1
1204           60.102951      1204           1
1205           60.198711      1205           1
1206           60.202461      1206           1
1207           60.297974      1207           1
1208           60.301865      1208           1
1209           60.396893      1209           1
1210           60.400921      1210           1
1211           60.499100      1211           1
1212           60.503250      1212           1
1213           60.595608      1213           1
1214           60.599850      1214           1
1215           60.698551      1215           1
1216           60.702885      1216           1
1217           60.799133      1217           1
1218           60.803547      1218           1
1219           60.899158      1219           1
1220           60.903641      1220           1
1221           60.998318      1221           1
1222           61.002853      1222           1
1223           61.098686      1223           1
1224           61.103275      1224           1
1225           61.195549      1225           1
1226           61.200184      1226           1
1227           61.295708      1227           1
1228           61.300388      1228           1
1229           61.399109      1229           1
1230           61.403828      1230           1
1231           61.498520      1231           1
1232           61.503273      1232           1
1233           61.598572      1233           1
1234           61.603352      1234           1
1235           61.699268      1235           1
1236           61.702904      1236           1
1237           61.798626      1237           1
1238           61.802292      1238           1
1239           61.898388      1239           1
1240           61.903275      1240           1
1241           61.996387      1241           1
1242           62.001305      1242           1
1243           62.097511      1243           1
1244           62.102463      1244           1
1245           62.196812      1245           1
1246           62.201797      1246           1
1247           62.296928      1247           1
1248           62.301945      1248           1
1249           62.397636      1249           1
1250           62.402691      1250           1
1251           62.498623      1251           1
1252           62.502464      1252           1
1253           62.598114      1253           1
1254           62.601910      1254           1
1255           62.698845      1255           1
1256           62.702484      1256           1
1257           62.797821      1257           1
1258           62.801239      1258           1
1259           62.897034      1259           1
1260           62.900238      1260           1
1261           62.997009      1261           1
1262           63.000065      1262           1
1263           63.099350      1263           1
1264           63.102325      1264           1
1265           63.197098      1265           1
1266           63.200054      1266           1
1267           63.297955      1267           1
1268           63.300938      1268           1
1269           63.397018      1269           1
1270           63.400043      1270           1
1271           63.497719      1271           1
1272           63.500797      1272           1
1273           63.597069      1273           1
1274           63.600201      1274           1
1275           63.698181      1275           1
1276           63.701366      1276           1
1277           63.797832      1277           1
1278           63.801075      1278           1
1279           63.899082      1279           1
1280           63.902374      1280           1
1281           63.998466      1281           1
1282           64.001808      1282           1
1283           64.099220      1283           1
1284           64.102600      1284           1
1285           64.198006      1285           1
1286           64.201439      1286           1
1287           64.298363      1287           1
1288           64.301857      1288           1
1289           64.396996      1289           1
1290           64.400558      1290           1
1291           64.497734      1291           1
1292           64.501373      1292           1
1293           64.597267      1293           1
1294           64.601006      1294           1
1295           64.699341      1295           1
1296           64.703163      1296           1
1297           64.797836      1297           1
1298           64.802208      1298           1
1299           64.895775      1299           1
1300           64.900726      1300           1
1301           64.999107      1301           1
1302           65.004463      1302           1
1303           65.096153      1303           1
1304           65.101578      1304           1
1305           65.194511      1305           1
1306           65.200012      1306           1
1307           65.294067      1307           1
1308           65.299629      1308           1
1309           65.394691      1309           1
1310           65.400307      1310           1
1311           65.496292      1311           1
1312           65.501961      1312           1
1313           65.598442      1313           1
1314           65.604019      1314           1
1315           65.695824      1315           1
1316           65.701057      1316           1
1317           65.797157      1317           1
1318           65.802040      1318           1
1319           65.896500      1319           1
1320           65.901062      1320           1
1321           65.998512      1321           1
1322           66.003067      1322           1
1323           66.095955      1323           1
1324           66.100761      1324           1
1325           66.196014      1325           1
1326           66.200500      1326           1
1327           66.297638      1327           1
1328           66.301964      1328           1
1329           66.398659      1329           1
1330           66.402740      1330           1
1331           66.497978      1331           1
1332           66.501846      1332           1
1333           66.596672      1333           1
1334           66.600403      1334           1
1335           66.696022      1335           1
1336           66.699661      1336           1
1337           66.797050      1337           1
1338           66.800629      1338           1
1339           66.897072      1339           1
1340           66.900635      1340           1
1341           66.996864      1341           1
1342           67.000427      1342           1
1343           67.096977      1343           1
1344           67.100563      1344           1
1345           67.197739      1345           1
1346           67.201355      1346           1
1347           67.299347      1347           1
1348           67.302994      1348           1
1349           67.398209      1349           1
1350           67.401886      1350           1
1351           67.498039      1351           1
1352           67.501755      1352           1
1353           67.599007      1353           1
1354           67.602768      1354           1
1355           67.697502      1355           1
1356           67.701317      1356           1
1357           67.797569      1357           1
1358           67.801453      1358           1
1359           67.895668      1359           1
1360           67.899635      1360           1
1361           67.996025      1361           1
1362           68.000092      1362           1
1363           68.098991      1363           1
1364           68.103165      1364           1
1365           68.195976      1365           1
1366           68.200233      1366           1
1367           68.299057      1367           1
1368           68.303391      1368           1
1369           68.399063      1369           1
1370           68.403427      1370           1
1371           68.495331      1371           1
1372           68.499718      1372           1
1373           68.596405      1373           1
1374           68.600807      1374           1
1375           68.697884      1375           1
1376           68.702301      1376           1
1377           68.795311      1377           1
1378           68.799751      1378           1
1379           68.897575      1379           1
1380           68.902031      1380           1
1381           68.995827      1381           1
1382           69.000305      1382           1
1383           69.099022      1383           1
1384           69.103523      1384           1
1385           69.194160      1385           1
1386           69.199623      1386           1
1387           69.295479      1387           1
1388           69.301437      1388           1
1389           69.396385      1389           1
1390           69.401413      1390           1
1391           69.495796      1391           1
1392           69.500328      1392           1
1393           69.598282      1393           1
1394           69.602631      1394           1
1395           69.695694      1395           1
1396           69.699799      1396           1
1397           69.799469      1397           1
1398           69.803352      1398           1
1399           69.898781      1399           1
1400           69.902550      1400           1
1401           69.996033      1401           1
1402           69.999756      1402           1
1403           70.096512      1403           1
1404           70.100235      1404           1
1405           70.197105      1405           1
1406           70.200829      1406           1
1407           70.297966      1407           1
1408           70.301704      1408           1
1409           70.399284      1409           1
1410           70.403053      1410           1
1411           70.497482      1411           1
1412           70.501274      1412           1
1413           70.596344      1413           1
1414           70.600159      1414           1
1415           70.695747      1415           1
1416           70.699577      1416           1
1417           70.799263      1417           1
1418           70.803101      1418           1
1419           70.899025      1419           1
1420           70.902863      1420           1
1421           70.999016      1421           1
1422           71.002869      1422           1
1423           71.095772      1423           1
1424           71.099663      1424           1
1425           71.197426      1425           1
1426           71.201363      1426           1
1427           71.296288      1427           1
1428           71.300262      1428           1
1429           71.396271      1429           1
1430           71.400299      1430           1
1431           71.497467      1431           1
1432           71.501541      1432           1
1433           71.595886      1433           1
1434           71.600014      1434           1
1435           71.695580      1435           1
1436           71.699753      1436           1
1437           71.796181      1437           1
1438           71.800385      1438           1
1439           71.897385      1439           1
1440           71.901619      1440           1
1441           71.999176      1441           1
1442           72.003433      1442           1
1443           72.097397      1443           1
1444           72.101685      1444           1
1445           72.196312      1445           1
1446           72.200630      1446           1
1447           72.295929      1447           1
1448           72.300278      1448           1
1449           72.396133      1449           1
1450           72.400505      1450           1
1451           72.496887      1451           1
1452           72.501282      1452           1
1453           72.598228      1453           1
1454           72.602646      1454           1
1455           72.695869      1455           1
1456           72.700325      1456           1
1457           72.798912      1457           1
1458           72.803421      1458           1
1459           72.898674      1459           1
1460           72.903236      1460           1
1461           72.997437      1461           1
1462           73.001984      1462           1
1463           73.095413      1463           1
1464           73.100021      1464           1
1465           73.196739      1465           1
1466           73.201279      1466           1
1467           73.298820      1467           1
1468           73.303131      1468           1
1469           73.398743      1469           1
1470           73.402756      1470           1
1471           73.496010      1471           1
1472           73.499779      1472           1
1473           73.598969      1473           1
1474           73.602562      1474           1
1475           73.698006      1475           1
1476           73.701500      1476           1
1477           73.798454      1477           1
1478           73.801895      1478           1
1479           73.898056      1479           1
1480           73.901489      1480           1
1481           73.997803      1481           1
1482           74.001251      1482           1
1483           74.098282      1483           1
1484           74.101761      1484           1
1485           74.196243      1485           1
1486           74.199760      1486           1
1487           74.298874      1487           1
1488           74.302437      1488           1
1489           74.399155      1489           1
1490           74.402763      1490           1
1491           74.497017      1491           1
1492           74.500664      1492           1
1493           74.596191      1493           1
1494           74.599892      1494           1
1495           74.696945      1495           1
1496           74.700706      1496           1
1497           74.795631      1497           1
1498           74.799461      1498           1
1499           74.895981      1499           1
1500           74.899872      1500           1
1501           74.997833      1501           1
1502           75.001778      1502           1
1503           75.097069      1503           1
1504           75.101067      1504           1
1505           75.197525      1505           1
1506           75.201569      1506           1
1507           75.299156      1507           1
1508           75.303246      1508           1
1509           75.397606      1509           1
1510           75.401726      1510           1
1511           75.496498      1511           1
1512           75.500618      1512           1
1513           75.595306      1513           1
1514           75.599419      1514           1
1515           75.698036      1515           1
1516           75.702148      1516           1
1517           75.796631      1517           1
1518           75.800735      1518           1
1519           75.899315      1519           1
1520           75.903419      1520           1
1521           75.998055      1521           1
1522           76.002182      1522           1
1523           76.097130      1523           1
1524           76.101273      1524           1
1525           76.196556      1525           1
1526           76.200706      1526           1
1527           76.296249      1527           1
1528           76.300407      1528           1
1529           76.396255      1529           1
1530           76.400436      1530           1
1531           76.496902      1531           1
1532           76.501122      1532           1
1533           76.598694      1533           1
1534           76.602966      1534           1
1535           76.697784      1535           1
1536           76.702141      1536           1
1537           76.798668      1537           1
1538           76.803085      1538           1
1539           76.896393      1539           1
1540           76.900848      1540           1
1541           76.998802      1541           1
1542           77.003242      1542           1
1543           77.095863      1543           1
1544           77.100235      1544           1
1545           77.195580      1545           1
1546           77.199867      1546           1
1547           77.297424      1547           1
1548           77.301613      1548           1
1549           77.396805      1549           1
1550           77.400887      1550           1
1551           77.497551      1551           1
1552           77.501526      1552           1
1553           77.595688      1553           1
1554           77.599564      1554           1
1555           77.699280      1555           1
1556           77.703079      1556           1
1557           77.797119      1557           1
1558           77.800842      1558           1
1559           77.897072      1559           1
1560           77.900742      1560           1
1561           77.995789      1561           1
1562           77.999428      1562           1
1563           78.097328      1563           1
1564           78.100937      1564           1
1565           78.198364      1565           1
1566           78.201965      1566           1
1567           78.299225      1567           1
1568           78.302826      1568           1
1569           78.396614      1569           1
1570           78.400230      1570           1
1571           78.498070      1571           1
1572           78.501701      1572           1
1573           78.596519      1573           1
1574           78.600182      1574           1
1575           78.699333      1575           1
1576           78.703018      1576           1
1577           78.799156      1577           1
1578           78.802864      1578           1
1579           78.895920      1579           1
1580           78.899651      1580           1
1581           78.997055      1581           1
1582           79.000816      1582           1
1583           79.098854      1583           1
1584           79.102638      1584           1
1585           79.197380      1585           1
1586           79.201180      1586           1
1587           79.296219      1587           1
1588           79.300026      1588           1
1589           79.399139      1589           1
1590           79.402962      1590           1
1591           79.498566      1591           1
1592           79.502396      1592           1
1593           79.598297      1593           1
1594           79.602135      1594           1
1595           79.698227      1595           1
1596           79.702072      1596           1
1597           79.798187      1597           1
1598           79.802032      1598           1
1599           79.897949      1599           1
1600           79.901779      1600           1
1601           79.997330      1601           1
1602           80.001144      1602           1
1603           80.096260      1603           1
1604           80.100060      1604           1
1605           80.198669      1605           1
1606           80.202461      1606           1
1607           80.297554      1607           1
1608           80.301369      1608           1
1609           80.397369      1609           1
1610           80.401230      1610           1
1611           80.498566      1611           1
1612           80.502487      1612           1
1613           80.597290      1613           1
1614           80.601265      1614           1
1615           80.697189      1615           1
1616           80.701202      1616           1
1617           80.797874      1617           1
1618           80.801910      1618           1
1619           80.898842      1619           1
1620           80.902878      1620           1
1621           80.995499      1621           1
1622           80.999519      1622           1
1623           81.095490      1623           1
1624           81.099472      1624           1
1625           81.198395      1625           1
1626           81.202332      1626           1
1627           81.295998      1627           1
1628           81.299866      1628           1
1629           81.395844      1629           1
1630           81.399643      1630           1
1631           81.497581      1631           1
1632           81.501312      1632           1
1633           81.597313      1633           1
1634           81.600975      1634           1
1635           81.698814      1635           1
1636           81.702408      1636           1
1637           81.798645      1637           1
1638           81.802177      1638           1
1639           81.897079      1639           1
1640           81.900574      1640           1
1641           81.997871      1641           1
1642           82.001335      1642           1
1643           82.097832      1643           1
1644           82.101265      1644           1
1645           82.197243      1645           1
1646           82.200668      1646           1
1647           82.296303      1647           1
1648           82.299713      1648           1
1649           82.398483      1649           1
1650           82.401886      1650           1
1651           82.497002      1651           1
1652           82.500397      1652           1
1653           82.598633      1653           1
1654           82.602013      1654           1
1655           82.696518      1655           1
1656           82.699890      1656           1
1657           82.797386      1657           1
1658           82.800743      1658           1
1659           82.897896      1659           1
1660           82.901245      1660           1
1661           82.998161      1661           1
1662           83.001503      1662           1
1663           83.098282      1663           1
1664           83.101624      1664           1
1665           83.198311      1665           1
1666           83.201645      1666           1
1667           83.298271      1667           1
1668           83.301598      1668           1
1669           83.398163      1669           1
1670           83.401489      1670           1
1671           83.497871      1671           1
1672           83.501190      1672           1
1673           83.597122      1673           1
1674           83.600418      1674           1
1675           83.698853      1675           1
1676           83.702118      1676           1
1677           83.796303      1677           1
1678           83.799530      1678           1
1679           83.899185      1679           1
1680           83.902390      1680           1
1681           83.998177      1681           1
1682           84.001358      1682           1
1683           84.097023      1683           1
1684           84.100212      1684           1
1685           84.196312      1685           1
1686           84.199532      1686           1
1687           84.296509      1687           1
1688           84.299759      1688           1
1689           84.397743      1689           1
1690           84.401024      1690           1
1691           84.496628      1691           1
1692           84.499939      1692           1
1693           84.596344      1693           1
1694           84.599678      1694           1
1695           84.696762      1695           1
1696           84.700119      1696           1
1697           84.797684      1697           1
1698           84.801056      1698           1
1699           84.898834      1699           1
1700           84.902206      1700           1
1701           84.996559      1701           1
1702           84.999931      1702           1
1703           85.097366      1703           1
1704           85.100723      1704           1
1705           85.197601      1705           1
1706           85.200928      1706           1
1707           85.296936      1707           1
1708           85.300224      1708           1
1709           85.398407      1709           1
1710           85.401657      1710           1
1711           85.498558      1711           1
1712           85.501770      1712           1
1713           85.597412      1713           1
1714           85.600586      1714           1
1715           85.698296      1715           1
1716           85.701431      1716           1
1717           85.798248      1717           1
1718           85.801361      1718           1
1719           85.897491      1719           1
1720           85.900581      1720           1
1721           85.999207      1721           1
1722           86.002274      1722           1
1723           86.097252      1723           1
1724           86.100304      1724           1
1725           86.197800      1725           1
1726           86.200836      1726           1
1727           86.297852      1727           1
1728           86.300873      1728           1
1729           86.397476      1729           1
1730           86.400490      1730           1
1731           86.496696      1731           1
1732           86.499695      1732           1
1733           86.598503      1733           1
1734           86.601494      1734           1
1735           86.696907      1735           1
1736           86.699883      1736           1
1737           86.797882      1737           1
1738           86.800850      1738           1
1739           86.898445      1739           1
1740           86.901398      1740           1
1741           86.998558      1741           1
1742           87.001495      1742           1
1743           87.098228      1743           1
1744           87.101151      1744           1
1745           87.197449      1745           1
1746           87.200356      1746           1
1747           87.299103      1747           1
1748           87.301994      1748           1
1749           87.397354      1749           1
1750           87.400238      1750           1
1751           87.498009      1751           1
1752           87.500877      1752           1
1753           87.598175      1753           1
1754           87.601028      1754           1
1755           87.697823      1755           1
1756           87.700661      1756           1
1757           87.796829      1757           1
1758           87.799644      1758           1
1759           87.897827      1759           1
1760           87.900620      1760           1
1761           87.997841      1761           1
1762           88.000610      1762           1
1763           88.096748      1763           1
1764           88.099480      1764           1
1765           88.197258      1765           1
1766           88.199959      1766           1
1767           88.296837      1767           1
1768           88.299522      1768           1
1769           88.398468      1769           1
1770           88.401131      1770           1
1771           88.497032      1771           1
1772           88.499695      1772           1
1773           88.598061      1773           1
1774           88.600716      1774           1
1775           88.698975      1775           1
1776           88.701630      1776           1
1777           88.797234      1777           1
1778           88.799889      1778           1
1779           88.898125      1779           1
1780           88.900780      1780           1
1781           88.998909      1781           1
1782           89.001556      1782           1
1783           89.096809      1783           1
1784           89.099449      1784           1
1785           89.196999      1785           1
1786           89.199631      1786           1
1787           89.296768      1787           1
1788           89.299385      1788           1
1789           89.398796      1789           1
1790           89.401405      1790           1
1791           89.497932      1791           1
1792           89.500534      1792           1
1793           89.596710      1793           1
1794           89.599304      1794           1
1795           89.697548      1795           1
1796           89.700119      1796           1
1797           89.797577      1797           1
1798           89.800133      1798           1
1799           89.899040      1799           1
1800           89.901566      1800           1
1801           89.999130      1801           1
1802           90.001610      1802           1
1803           90.097702      1803           1
1804           90.100151      1804           1
1805           90.197128      1805           1
1806           90.199532      1806           1
1807           90.297394      1807           1
1808           90.299759      1808           1
1809           90.398682      1809           1
1810           90.401016      1810           1
1811           90.498909      1811           1
1812           90.501228      1812           1
1813           90.598305      1813           1
1814           90.600609      1814           1
1815           90.696991      1815           1
1816           90.699280      1816           1
1817           90.797256      1817           1
1818           90.799522      1818           1
1819           90.898933      1819           1
1820           90.901184      1820           1
1821           90.997368      1821           1
1822           90.999596      1822           1
1823           91.099159      1823           1
1824           91.101357      1824           1
1825           91.197655      1825           1
1826           91.199829      1826           1
1827           91.297287      1827           1
1828           91.299438      1828           1
1829           91.398048      1829           1
1830           91.400177      1830           1
1831           91.497818      1831           1
1832           91.499931      1832           1
1833           91.598694      1833           1
1834           91.600784      1834           1
1835           91.698486      1835           1
1836           91.700554      1836           1
1837           91.799164      1837           1
1838           91.801201      1838           1
1839           91.898529      1839           1
1840           91.900543      1840           1
1841           91.994476      1841           1
1842           91.999649      1842           1
1843           92.093620      1843           1
1844           92.099930      1844           1
1845           92.194557      1845           1
1846           92.200867      1846           1
1847           92.295479      1847           1
1848           92.301788      1848           1
1849           92.396393      1849           1
1850           92.402695      1850           1
1851           92.497284      1851           1
1852           92.503593      1852           1
1853           92.598167      1853           1
1854           92.604469      1854           1
1855           92.699028      1855           1
1856           92.705330      1856           1
1857           92.793571      1857           1
1858           92.799873      1858           1
1859           92.894409      1859           1
1860           92.900711      1860           1
1861           92.995239      1861           1
1862           93.001541      1862           1
1863           93.096054      1863           1
1864           93.102356      1864           1
1865           93.196861      1865           1
1866           93.203163      1866           1
1867           93.297661      1867           1
1868           93.303963      1868           1
1869           93.398453      1869           1
1870           93.404747      1870           1
1871           93.492920      1871           1
1872           93.499222      1872           1
1873           93.593674      1873           1
1874           93.599968      1874           1
1875           93.694130      1875           1
1876           93.700340      1876           1
1877           93.798470      1877           1
1878           93.804520      1878           1
1879           93.894218      1879           1
1880           93.900131      1880           1
1881           93.993454      1881           1
1882           93.999214      1882           1
1883           94.095642      1883           1
1884           94.101227      1884           1
1885           94.194656      1885           1
1886           94.200058      1886           1
1887           94.295464      1887           1
1888           94.300659      1888           1
1889           94.397217      1889           1
1890           94.402184      1890           1
1891           94.499069      1891           1
1892           94.503792      1892           1
1893           94.595764      1893           1
1894           94.600235      1894           1
1895           94.695816      1895           1
1896           94.700272      1896           1
1897           94.796463      1897           1
1898           94.801163      1898           1
1899           94.897522      1899           1
1900           94.902451      1900           1
1901           94.998192      1901           1
1902           95.003334      1902           1
1903           95.097672      1903           1
1904           95.103004      1904           1
1905           95.195129      1905           1
1906           95.200623      1906           1
1907           95.295425      1907           1
1908           95.301079      1908           1
1909           95.398483      1909           1
1910           95.404289      1910           1
1911           95.498398      1911           1
1912           95.504356      1912           1
1913           95.594734      1913           1
1914           95.600830      1914           1
1915           95.693306      1915           1
1916           95.699539      1916           1
1917           95.794113      1917           1
1918           95.800484      1918           1
1919           95.896873      1919           1
1920           95.903358      1920           1
1921           95.994904      1921           1
1922           96.001495      1922           1
1923           96.094551      1923           1
1924           96.101250      1924           1
1925           96.195740      1925           1
1926           96.202538      1926           1
1927           96.298393      1927           1
1928           96.305290      1928           1
1929           96.395454      1929           1
1930           96.402435      1930           1
1931           96.493706      1931           1
1932           96.500771      1932           1
1933           96.593170      1933           1
1934           96.600319      1934           1
1935           96.693840      1935           1
1936           96.701073      1936           1
1937           96.795700      1937           1
1938           96.803024      1938           1
1939           96.898735      1939           1
1940           96.906143      1940           1
1941           96.995483      1941           1
1942           97.002975      1942           1
1943           97.098320      1943           1
1944           97.104607      1944           1
1945           97.197815      1945           1
1946           97.204727      1946           1
1947           97.295593      1947           1
1948           97.301880      1948           1
1949           97.391281      1949           1
1950           97.399094      1950           1
1951           97.493439      1951           1
1952           97.501350      1952           1
1953           97.596687      1953           1
1954           97.604630      1954           1
1955           97.691605      1955           1
1956           97.699478      1956           1
1957           97.797691      1957           1
1958           97.804039      1958           1
1959           97.898926      1959           1
1960           97.905212      1960           1
1961           97.993690      1961           1
1962           98.000038      1962           1
1963           98.094635      1963           1
1964           98.100922      1964           1
1965           98.195473      1965           1
1966           98.201759      1966           1
1967           98.296059      1967           1
1968           98.302345      1968           1
1969           98.396629      1969           1
1970           98.402916      1970           1
1971           98.497223      1971           1
1972           98.503510      1972           1
1973           98.597824      1973           1
1974           98.604111      1974           1
1975           98.698425      1975           1
1976           98.704712      1976           1
1977           98.799034      1977           1
1978           98.805321      1978           1
1979           98.893364      1979           1
1980           98.899651      1980           1
1981           98.993980      1981           1
1982           99.000267      1982           1
1983           99.094604      1983           1
1984           99.100891      1984           1
1985           99.195229      1985           1
1986           99.201515      1986           1
1987           99.295860      1987           1
1988           99.302155      1988           1
1989           99.396500      1989           1
1990           99.402786      1990           1
1991           99.497131      1991           1
1992           99.503426      1992           1
1993           99.597763      1993           1
1994           99.604057      1994           1
1995           99.698402      1995           1
1996           99.704689      1996           1
1997           99.799042      1997           1
1998           99.805336      1998           1
1999           99.893394      1999           1
2000           99.899689      2000           1
2001           99.993752      2001           1
2002           99.999969      2002           1

In this case, the Unit System is not attached to the data_source, but it can be directly assigned to the results. As we are employing the dpf.Model API, we only need to assign the Unit System once (the Model will assign it for the rest of the results).

Results from the matsum branch of the binout file are a FieldsContainer on a LabelSpace comprised by part IDs. Extract part kinetic energy for all parts:

PKE_op = model.results.part_kinetic_energy()
PKE_op.inputs.unit_system.connect(dpf.unit_systems.solver_mks)
PKE = PKE_op.eval()
print(PKE)
DPF  Fields Container
  with 129 field(s)
  defined on labels: part

  with:
  - field 0 {part:  1} with TimeFreq_steps location, 1 components and 101 entities.
  - field 1 {part:  2} with TimeFreq_steps location, 1 components and 101 entities.
  - field 2 {part:  3} with TimeFreq_steps location, 1 components and 101 entities.
  - field 3 {part:  4} with TimeFreq_steps location, 1 components and 101 entities.
  - field 4 {part:  5} with TimeFreq_steps location, 1 components and 101 entities.
  - field 5 {part:  6} with TimeFreq_steps location, 1 components and 101 entities.
  - field 6 {part:  7} with TimeFreq_steps location, 1 components and 101 entities.
  - field 7 {part:  8} with TimeFreq_steps location, 1 components and 101 entities.
  - field 8 {part:  9} with TimeFreq_steps location, 1 components and 101 entities.
  - field 9 {part:  10} with TimeFreq_steps location, 1 components and 101 entities.
  - field 10 {part:  11} with TimeFreq_steps location, 1 components and 101 entities.
  - field 11 {part:  12} with TimeFreq_steps location, 1 components and 101 entities.
  - field 12 {part:  13} with TimeFreq_steps location, 1 components and 101 entities.
  - field 13 {part:  14} with TimeFreq_steps location, 1 components and 101 entities.
  - field 14 {part:  15} with TimeFreq_steps location, 1 components and 101 entities.
  - field 15 {part:  16} with TimeFreq_steps location, 1 components and 101 entities.
  - field 16 {part:  17} with TimeFreq_steps location, 1 components and 101 entities.
  - field 17 {part:  18} with TimeFreq_steps location, 1 components and 101 entities.
  - field 18 {part:  19} with TimeFreq_steps location, 1 components and 101 entities.
  - field 19 {part:  20} with TimeFreq_steps location, 1 components and 101 entities.
  - field 20 {part:  21} with TimeFreq_steps location, 1 components and 101 entities.
  - field 21 {part:  22} with TimeFreq_steps location, 1 components and 101 entities.
  - field 22 {part:  23} with TimeFreq_steps location, 1 components and 101 entities.
  - field 23 {part:  24} with TimeFreq_steps location, 1 components and 101 entities.
  - field 24 {part:  25} with TimeFreq_steps location, 1 components and 101 entities.
  - field 25 {part:  26} with TimeFreq_steps location, 1 components and 101 entities.
  - field 26 {part:  27} with TimeFreq_steps location, 1 components and 101 entities.
  - field 27 {part:  28} with TimeFreq_steps location, 1 components and 101 entities.
  - field 28 {part:  29} with TimeFreq_steps location, 1 components and 101 entities.
  - field 29 {part:  30} with TimeFreq_steps location, 1 components and 101 entities.
  - field 30 {part:  31} with TimeFreq_steps location, 1 components and 101 entities.
  - field 31 {part:  32} with TimeFreq_steps location, 1 components and 101 entities.
  - field 32 {part:  33} with TimeFreq_steps location, 1 components and 101 entities.
  - field 33 {part:  34} with TimeFreq_steps location, 1 components and 101 entities.
  - field 34 {part:  35} with TimeFreq_steps location, 1 components and 101 entities.
  - field 35 {part:  36} with TimeFreq_steps location, 1 components and 101 entities.
  - field 36 {part:  37} with TimeFreq_steps location, 1 components and 101 entities.
  - field 37 {part:  38} with TimeFreq_steps location, 1 components and 101 entities.
  - field 38 {part:  39} with TimeFreq_steps location, 1 components and 101 entities.
  - field 39 {part:  40} with TimeFreq_steps location, 1 components and 101 entities.
  - field 40 {part:  41} with TimeFreq_steps location, 1 components and 101 entities.
  - field 41 {part:  42} with TimeFreq_steps location, 1 components and 101 entities.
  - field 42 {part:  43} with TimeFreq_steps location, 1 components and 101 entities.
  - field 43 {part:  44} with TimeFreq_steps location, 1 components and 101 entities.
  - field 44 {part:  45} with TimeFreq_steps location, 1 components and 101 entities.
  - field 45 {part:  46} with TimeFreq_steps location, 1 components and 101 entities.
  - field 46 {part:  47} with TimeFreq_steps location, 1 components and 101 entities.
  - field 47 {part:  48} with TimeFreq_steps location, 1 components and 101 entities.
  - field 48 {part:  49} with TimeFreq_steps location, 1 components and 101 entities.
  - field 49 {part:  50} with TimeFreq_steps location, 1 components and 101 entities.
  - field 50 {part:  51} with TimeFreq_steps location, 1 components and 101 entities.
  - field 51 {part:  52} with TimeFreq_steps location, 1 components and 101 entities.
  - field 52 {part:  53} with TimeFreq_steps location, 1 components and 101 entities.
  - field 53 {part:  54} with TimeFreq_steps location, 1 components and 101 entities.
  - field 54 {part:  55} with TimeFreq_steps location, 1 components and 101 entities.
  - field 55 {part:  56} with TimeFreq_steps location, 1 components and 101 entities.
  - field 56 {part:  57} with TimeFreq_steps location, 1 components and 101 entities.
  - field 57 {part:  58} with TimeFreq_steps location, 1 components and 101 entities.
  - field 58 {part:  59} with TimeFreq_steps location, 1 components and 101 entities.
  - field 59 {part:  60} with TimeFreq_steps location, 1 components and 101 entities.
  - field 60 {part:  61} with TimeFreq_steps location, 1 components and 101 entities.
  - field 61 {part:  62} with TimeFreq_steps location, 1 components and 101 entities.
  - field 62 {part:  63} with TimeFreq_steps location, 1 components and 101 entities.
  - field 63 {part:  64} with TimeFreq_steps location, 1 components and 101 entities.
  - field 64 {part:  65} with TimeFreq_steps location, 1 components and 101 entities.
  - field 65 {part:  66} with TimeFreq_steps location, 1 components and 101 entities.
  - field 66 {part:  67} with TimeFreq_steps location, 1 components and 101 entities.
  - field 67 {part:  68} with TimeFreq_steps location, 1 components and 101 entities.
  - field 68 {part:  69} with TimeFreq_steps location, 1 components and 101 entities.
  - field 69 {part:  70} with TimeFreq_steps location, 1 components and 101 entities.
  - field 70 {part:  71} with TimeFreq_steps location, 1 components and 101 entities.
  - field 71 {part:  72} with TimeFreq_steps location, 1 components and 101 entities.
  - field 72 {part:  73} with TimeFreq_steps location, 1 components and 101 entities.
  - field 73 {part:  74} with TimeFreq_steps location, 1 components and 101 entities.
  - field 74 {part:  75} with TimeFreq_steps location, 1 components and 101 entities.
  - field 75 {part:  76} with TimeFreq_steps location, 1 components and 101 entities.
  - field 76 {part:  77} with TimeFreq_steps location, 1 components and 101 entities.
  - field 77 {part:  78} with TimeFreq_steps location, 1 components and 101 entities.
  - field 78 {part:  79} with TimeFreq_steps location, 1 components and 101 entities.
  - field 79 {part:  80} with TimeFreq_steps location, 1 components and 101 entities.
  - field 80 {part:  81} with TimeFreq_steps location, 1 components and 101 entities.
  - field 81 {part:  82} with TimeFreq_steps location, 1 components and 101 entities.
  - field 82 {part:  83} with TimeFreq_steps location, 1 components and 101 entities.
  - field 83 {part:  84} with TimeFreq_steps location, 1 components and 101 entities.
  - field 84 {part:  85} with TimeFreq_steps location, 1 components and 101 entities.
  - field 85 {part:  86} with TimeFreq_steps location, 1 components and 101 entities.
  - field 86 {part:  87} with TimeFreq_steps location, 1 components and 101 entities.
  - field 87 {part:  88} with TimeFreq_steps location, 1 components and 101 entities.
  - field 88 {part:  89} with TimeFreq_steps location, 1 components and 101 entities.
  - field 89 {part:  90} with TimeFreq_steps location, 1 components and 101 entities.
  - field 90 {part:  91} with TimeFreq_steps location, 1 components and 101 entities.
  - field 91 {part:  92} with TimeFreq_steps location, 1 components and 101 entities.
  - field 92 {part:  93} with TimeFreq_steps location, 1 components and 101 entities.
  - field 93 {part:  301} with TimeFreq_steps location, 1 components and 101 entities.
  - field 94 {part:  303} with TimeFreq_steps location, 1 components and 101 entities.
  - field 95 {part:  304} with TimeFreq_steps location, 1 components and 101 entities.
  - field 96 {part:  305} with TimeFreq_steps location, 1 components and 101 entities.
  - field 97 {part:  306} with TimeFreq_steps location, 1 components and 101 entities.
  - field 98 {part:  307} with TimeFreq_steps location, 1 components and 101 entities.
  - field 99 {part:  308} with TimeFreq_steps location, 1 components and 101 entities.
  - field 100 {part:  309} with TimeFreq_steps location, 1 components and 101 entities.
  - field 101 {part:  310} with TimeFreq_steps location, 1 components and 101 entities.
  - field 102 {part:  311} with TimeFreq_steps location, 1 components and 101 entities.
  - field 103 {part:  312} with TimeFreq_steps location, 1 components and 101 entities.
  - field 104 {part:  718} with TimeFreq_steps location, 1 components and 101 entities.
  - field 105 {part:  1500} with TimeFreq_steps location, 1 components and 101 entities.
  - field 106 {part:  1501} with TimeFreq_steps location, 1 components and 101 entities.
  - field 107 {part:  1502} with TimeFreq_steps location, 1 components and 101 entities.
  - field 108 {part:  1503} with TimeFreq_steps location, 1 components and 101 entities.
  - field 109 {part:  1504} with TimeFreq_steps location, 1 components and 101 entities.
  - field 110 {part:  1505} with TimeFreq_steps location, 1 components and 101 entities.
  - field 111 {part:  1506} with TimeFreq_steps location, 1 components and 101 entities.
  - field 112 {part:  1507} with TimeFreq_steps location, 1 components and 101 entities.
  - field 113 {part:  1508} with TimeFreq_steps location, 1 components and 101 entities.
  - field 114 {part:  1509} with TimeFreq_steps location, 1 components and 101 entities.
  - field 115 {part:  1510} with TimeFreq_steps location, 1 components and 101 entities.
  - field 116 {part:  1511} with TimeFreq_steps location, 1 components and 101 entities.
  - field 117 {part:  1512} with TimeFreq_steps location, 1 components and 101 entities.
  - field 118 {part:  1513} with TimeFreq_steps location, 1 components and 101 entities.
  - field 119 {part:  1514} with TimeFreq_steps location, 1 components and 101 entities.
  - field 120 {part:  1515} with TimeFreq_steps location, 1 components and 101 entities.
  - field 121 {part:  1516} with TimeFreq_steps location, 1 components and 101 entities.
  - field 122 {part:  1517} with TimeFreq_steps location, 1 components and 101 entities.
  - field 123 {part:  1518} with TimeFreq_steps location, 1 components and 101 entities.
  - field 124 {part:  1519} with TimeFreq_steps location, 1 components and 101 entities.
  - field 125 {part:  1520} with TimeFreq_steps location, 1 components and 101 entities.
  - field 126 {part:  1521} with TimeFreq_steps location, 1 components and 101 entities.
  - field 127 {part:  1522} with TimeFreq_steps location, 1 components and 101 entities.
  - field 128 {part:  1523} with TimeFreq_steps location, 1 components and 101 entities.

Extract part internal energy for only a selected number of parts.

part_sco = dpf.Scoping(ids=[50, 1522], location="part")
PIE_op = model.results.part_internal_energy()
PIE_op.inputs.entity_scoping.connect(part_sco)
PIE = PIE_op.eval()

Plot part kinetic and internal energy for a selection of parts. In this case, the TimeFreqSupport of the matsum branch does not have all the time steps in the binout file. Thus, a rescoping operation is needed:

rescope_op = dpf.operators.scoping.rescope()
rescope_op.inputs.fields.connect(PKE.time_freq_support.time_frequencies)
rescope_op.inputs.mesh_scoping.connect(PKE[0].scoping)
t_field = rescope_op.outputs.fields_as_field()
t_vals = t_field.data

plt.plot(t_vals, PKE.get_field({"part": 50}).data, label="Kinetic, Part 50")
plt.plot(t_vals, PKE.get_field({"part": 1522}).data, label="Kinetic, Part 1522")
plt.plot(t_vals, PIE.get_field({"part": 50}).data, label="Internal, Part 50")
plt.plot(t_vals, PIE.get_field({"part": 1522}).data, label="Internal, Part 1522")
plt.xlabel("Time ({:s})".format(t_field.unit))
plt.ylabel("Energy ({:s})".format(PIE.get_field({"part": 50}).unit))
plt.legend()
plt.show()
00 lsdyna operators

Similarly, results from the rcforc branch of the binout file are a FieldsContainer on a LabelSpace comprised by interface IDs. Extract interface contact force for only one interface.

interface_sco = dpf.Scoping(ids=[19], location="interface")
FC_op = model.results.interface_contact_force()
FC_op.inputs.entity_scoping.connect(interface_sco)
FC = FC_op.eval()
print(FC)
DPF  Fields Container
  with 2 field(s)
  defined on labels: idtype interface

  with:
  - field 0 {interface:  19, idtype:  0} with TimeFreq_steps location, 3 components and 1001 entities.
  - field 1 {interface:  19, idtype:  1} with TimeFreq_steps location, 3 components and 1001 entities.

In addition to interface, the FieldsContainer is scoped on idtype (0 for the master side of the interface, 1 for the slave). Contact force is a vector, and the three components are available.

rescope_op = dpf.operators.scoping.rescope()
rescope_op.inputs.fields.connect(FC.time_freq_support.time_frequencies)
rescope_op.inputs.mesh_scoping.connect(FC[0].scoping)
t_field = rescope_op.outputs.fields_as_field()
t_vals = t_field.data

FX = FC.select_component(0)
FY = FC.select_component(1)
FZ = FC.select_component(2)

plt.plot(t_vals, FX.get_field({"interface": 19, "idtype": 0}).data, label="FX, slave")
plt.plot(t_vals, FX.get_field({"interface": 19, "idtype": 1}).data, label="FX, master")
plt.plot(t_vals, FY.get_field({"interface": 19, "idtype": 0}).data, label="FY, slave")
plt.plot(t_vals, FY.get_field({"interface": 19, "idtype": 1}).data, label="FY, master")
plt.plot(t_vals, FZ.get_field({"interface": 19, "idtype": 0}).data, label="FZ, slave")
plt.plot(t_vals, FZ.get_field({"interface": 19, "idtype": 1}).data, label="FZ, master")
plt.xlabel("Time ({:s})".format(t_field.unit))
plt.xlim([0, 10])
plt.ylabel("Contact Force ({:s})".format(FX.get_field({"interface": 19, "idtype": 0}).unit))
plt.legend()
plt.show()
00 lsdyna operators

Total running time of the script: (0 minutes 12.710 seconds)

Gallery generated by Sphinx-Gallery