|
| string | sp_statistics.file_name = 'sp_data.hdf5' |
| |
| | sp_statistics.f = h5py.File(file_name, 'r') |
| |
| | sp_statistics.test_id = f.attrs['test_id'] |
| |
| | sp_statistics.internal_radius = f.attrs['internal_radius'] |
| |
| | sp_statistics.external_radius = f.attrs['external_radius'] |
| |
| | sp_statistics.interface_radius = f.attrs['interface_radius'] |
| |
| | sp_statistics.thickness = f.attrs['thickness'] |
| |
| | sp_statistics.volume = f.attrs['volume'] |
| |
| | sp_statistics.real_probe_height = f.attrs['real_probe_height'] |
| |
| | sp_statistics.target_porosity = f.attrs['target_porosity'] |
| |
| | sp_statistics.porosity = f.attrs['porosity'] |
| |
| | sp_statistics.density = f.attrs['density'] |
| |
| | sp_statistics.height_factor = real_probe_height / thickness |
| |
| tuple | sp_statistics.porosity_factor = (1.0-target_porosity)/(1.0-porosity) |
| |
| float | sp_statistics.gram_factor = 1000.0 |
| |
| | sp_statistics.initial_radii = np.array(f['0'].get('radius')) |
| |
| | sp_statistics.initial_radii_3 = np.power(initial_radii,3) |
| |
| float | sp_statistics.initial_masses = 4.0/3.0 * np.pi * density * gram_factor * height_factor * porosity_factor * initial_radii_3 |
| |
| | sp_statistics.initial_mass = np.sum(initial_masses) |
| |
| list | sp_statistics.all_times = [] |
| |
| list | sp_statistics.all_pressures = [] |
| |
| list | sp_statistics.all_sps = [] |
| |
| int | sp_statistics.failure_step = 175 |
| |
| int | sp_statistics.max_step = 2*failure_step |
| |
| | sp_statistics.times = np.zeros(max_step+1) |
| |
| int | sp_statistics.p_rate = 5e10 |
| |
| float | sp_statistics.psi_factor = 0.000145038 |
| |
| list | sp_statistics.pressures = [p_rate*t*psi_factor for t in times] |
| |
| | sp_statistics.t_f = times[failure_step] |
| |
| float | sp_statistics.max_num_bonds = 6.0 |
| |
| | sp_statistics.sps = np.zeros(max_step+1) |
| |
| | sp_statistics.all_radii = np.array(f[str(i)].get('radius')) |
| |
| | sp_statistics.continuum_bonds = np.array(f[str(i)].get('current_continuum_bonds')) |
| |
| | sp_statistics.xs = np.array(f[str(i)].get('x')) |
| |
| | sp_statistics.ys = np.array(f[str(i)].get('y')) |
| |
| | sp_statistics.xs_2 = np.power(xs,2) |
| |
| | sp_statistics.ys_2 = np.power(ys,2) |
| |
| | sp_statistics.distance_2 = xs_2 + ys_2 |
| |
| | sp_statistics.weak_radius = interface_radius |
| |
| | sp_statistics.internal_radii = np.where(distance_2<weak_radius**2,all_radii,0.0) |
| |
| | sp_statistics.external_radii = np.where(distance_2>=weak_radius**2,all_radii,0.0) |
| |
| | sp_statistics.cont_internal_radii = np.where(continuum_bonds>numbonds,internal_radii,0.0) |
| |
| | sp_statistics.cont_external_radii = np.where(continuum_bonds>0.0,external_radii,0.0) |
| |
| | sp_statistics.total_radii_3 = np.power(cont_internal_radii+cont_external_radii,3) |
| |
| float | sp_statistics.masses = 4.0/3.0 * np.pi * density * gram_factor * height_factor * porosity_factor * total_radii_3 |
| |
| | sp_statistics.current_mass = np.sum(masses) |
| |
| | sp_statistics.current_sp = initial_mass - current_mass |
| |
| list | sp_statistics.exp_times = [] |
| |
| list | sp_statistics.exp_pressures = [] |
| |
| list | sp_statistics.exp_sps = [] |
| |
| | sp_statistics.fields = line.strip().split() |
| |
| string | sp_statistics.graph_name = 'sp_bonds_t.pdf' |
| |
| list | sp_statistics.graph_labels |
| |
| | sp_statistics.productions |
| |
| | sp_statistics.label |
| |
| | sp_statistics.loc |
| |
| | sp_statistics.prop |
| |
| | sp_statistics.bbox_inches |
| |
1.9.1