|
| string | sp_statistics_2.file_name = 'sp_data.hdf5' |
| |
| | sp_statistics_2.f = h5py.File(file_name, 'r') |
| |
| | sp_statistics_2.test_id = f.attrs['test_id'] |
| |
| | sp_statistics_2.internal_radius = f.attrs['internal_radius'] |
| |
| | sp_statistics_2.external_radius = f.attrs['external_radius'] |
| |
| | sp_statistics_2.interface_radius = f.attrs['interface_radius'] |
| |
| | sp_statistics_2.thickness = f.attrs['thickness'] |
| |
| | sp_statistics_2.volume = f.attrs['volume'] |
| |
| | sp_statistics_2.real_probe_height = f.attrs['real_probe_height'] |
| |
| | sp_statistics_2.target_porosity = f.attrs['target_porosity'] |
| |
| | sp_statistics_2.porosity = f.attrs['porosity'] |
| |
| | sp_statistics_2.density = f.attrs['density'] |
| |
| int | sp_statistics_2.failure_step = 300 |
| |
| int | sp_statistics_2.max_step = failure_step |
| |
| | sp_statistics_2.times = np.zeros(max_step) |
| |
| tuple | sp_statistics_2.delta_probe_radius = (interface_radius-internal_radius)/5.0 |
| |
| tuple | sp_statistics_2.probe_radius_1 = (internal_radius+delta_probe_radius*1)**2 |
| |
| tuple | sp_statistics_2.probe_radius_2 = (internal_radius+delta_probe_radius*2)**2 |
| |
| tuple | sp_statistics_2.probe_radius_3 = (internal_radius+delta_probe_radius*3)**2 |
| |
| tuple | sp_statistics_2.probe_radius_4 = (internal_radius+delta_probe_radius*4)**2 |
| |
| int | sp_statistics_2.interface_radius_2 = interface_radius**2 |
| |
| | sp_statistics_2.avg_num_intact_bonds_0 = np.zeros(max_step) |
| |
| | sp_statistics_2.avg_num_intact_bonds_1 = np.zeros(max_step) |
| |
| | sp_statistics_2.avg_num_intact_bonds_2 = np.zeros(max_step) |
| |
| | sp_statistics_2.avg_num_intact_bonds_3 = np.zeros(max_step) |
| |
| | sp_statistics_2.avg_num_intact_bonds_4 = np.zeros(max_step) |
| |
| | sp_statistics_2.avg_num_broken_bonds_0 = np.zeros(max_step) |
| |
| | sp_statistics_2.avg_num_broken_bonds_1 = np.zeros(max_step) |
| |
| | sp_statistics_2.avg_num_broken_bonds_2 = np.zeros(max_step) |
| |
| | sp_statistics_2.avg_num_broken_bonds_3 = np.zeros(max_step) |
| |
| | sp_statistics_2.avg_num_broken_bonds_4 = np.zeros(max_step) |
| |
| | sp_statistics_2.continuum_bonds = np.array(f[str(i)].get('current_continuum_bonds')) |
| |
| | sp_statistics_2.initial_continuum_bonds = np.array(f[str(i)].get('initial_continuum_bonds')) |
| |
| | sp_statistics_2.xs = np.array(f[str(i)].get('x')) |
| |
| | sp_statistics_2.ys = np.array(f[str(i)].get('y')) |
| |
| | sp_statistics_2.xs_2 = np.power(xs,2) |
| |
| | sp_statistics_2.ys_2 = np.power(ys,2) |
| |
| | sp_statistics_2.distance_2 = xs_2 + ys_2 |
| |
| | sp_statistics_2.num_intact_bonds_0 = np.where(distance_2<probe_radius_1,continuum_bonds,-100.0) |
| |
| | sp_statistics_2.index_to_delete = np.where(num_intact_bonds_0<0.0)[0] |
| |
| | sp_statistics_2.num_init_bonds_0 = np.where(distance_2<probe_radius_1,initial_continuum_bonds,-100.0) |
| |
| | sp_statistics_2.num_broken_bonds_0 = num_init_bonds_0 - num_intact_bonds_0 |
| |
| | sp_statistics_2.num_intact_bonds_1 = np.where(distance_2>=probe_radius_1,continuum_bonds,-100) |
| |
| | sp_statistics_2.num_init_bonds_1 = np.where(distance_2>=probe_radius_1,initial_continuum_bonds,-100) |
| |
| | sp_statistics_2.num_broken_bonds_1 = num_init_bonds_1 - num_intact_bonds_1 |
| |
| | sp_statistics_2.num_intact_bonds_2 = np.where(distance_2>=probe_radius_2,continuum_bonds,-100) |
| |
| | sp_statistics_2.num_init_bonds_2 = np.where(distance_2>=probe_radius_2,initial_continuum_bonds,-100) |
| |
| | sp_statistics_2.num_broken_bonds_2 = num_init_bonds_2 - num_intact_bonds_2 |
| |
| | sp_statistics_2.num_intact_bonds_3 = np.where(distance_2>=probe_radius_3,continuum_bonds,-100) |
| |
| | sp_statistics_2.num_init_bonds_3 = np.where(distance_2>=probe_radius_3,initial_continuum_bonds,-100) |
| |
| | sp_statistics_2.num_broken_bonds_3 = num_init_bonds_3 - num_intact_bonds_3 |
| |
| | sp_statistics_2.num_intact_bonds_4 = np.where(distance_2>probe_radius_4,continuum_bonds,-100.0) |
| |
| | sp_statistics_2.num_init_bonds_4 = np.where(distance_2>probe_radius_4,initial_continuum_bonds,-100.0) |
| |
| | sp_statistics_2.num_broken_bonds_4 = num_init_bonds_4 - num_intact_bonds_4 |
| |
| | sp_statistics_2.fig1 |
| |
| | sp_statistics_2.axs1 |
| |
| | sp_statistics_2.figsize |
| |
| string | sp_statistics_2.graph_name = 'average_number_intact_bonds_probe_radius.pdf' |
| |
| | sp_statistics_2.fig2 |
| |
| | sp_statistics_2.axs2 |
| |
1.9.1