'innerdonut', 'outerdonut', 'random', 'trans'] H = 9 el = 21 n_pc_tot = 20 loading = "1" prop = "yield" wrt_file = 'log_%s.txt' % (time.strftime("%Y-%m-%d_h%Hm%M")) f = h5py.File("spatial_stats.hdf5", 'w') f.close() """Gather data from vtk files""" for ii in xrange(len(set_id_cal)): vtk.read_euler(el, strt_cal[ii], ns_cal[ii], set_id_cal[ii], dir_cal[ii], wrt_file, 0) for ii in xrange(len(set_id_val)): vtk.read_euler(el, strt_val[ii], ns_val[ii], set_id_val[ii], dir_val[ii], wrt_file, 0) """get the data for the linkage""" f = h5py.File("responses.hdf5", 'w') f.close() f = h5py.File("responses_for_linkage.hdf5", 'w') f.close() for ii in xrange(len(set_id_cal)): # gr.Eeff(el, ns_cal[ii], set_id_cal[ii], step, dir_cal[ii], wrt_file) # gr.Eeff(el, ns_val[ii], set_id_val[ii], step, dir_val[ii], wrt_file) gr.FIP(el, ns_cal[ii], set_id_cal[ii], step, dir_cal[ii], wrt_file)
for step in xrange(1, 2): wrt_file = 'log_step%s_%s.txt' % (step, time.strftime("%Y-%m-%d_h%Hm%M")) ## The tensorID determines the type of tensor data read from the .vtk file ## if tensorID == 0, we read the stress tensor ## if tensorID == 1, we read the strain tensor ## if tensorID == 2, we read the plastic strain tensor tensor_ID = 1 ## Gather data from calibration vtk files dir_cal = 'cal_equi' vtk.read_euler(el, ns_cal, set_id_cal, step, dir_cal, wrt_file, 0) for comp in compl: vtk.read_meas(el, ns_cal, set_id_cal, step, comp, tensor_ID, dir_cal, wrt_file) ## Gather data from validation vtk files dir_val = 'val_basal' vtk.read_euler(el, ns_val, set_id_val, step, dir_val, wrt_file, 0) for comp in compl: vtk.read_meas(el, ns_val, set_id_val, step, comp, tensor_ID, dir_val, wrt_file)
set_id_cal = C['set_id_cal'] strt_cal = C['strt_cal'] ns_cal = C['ns_cal'] dir_cal = C['dir_cal'] set_id_val = C['set_id_val'] strt_val = C['strt_val'] ns_val = C['ns_val'] dir_val = C['dir_val'] f = h5py.File("spatial.hdf5", 'w') f.close() """Gather data from vtk files""" for ii in xrange(len(set_id_cal)): vtk.read_euler(strt_cal[ii], ns_cal[ii], set_id_cal[ii], dir_cal[ii], 1) for ii in xrange(len(set_id_val)): vtk.read_euler(strt_val[ii], ns_val[ii], set_id_val[ii], dir_val[ii], 1) """Compute GSH coefficients to create microstructure function in real and fourier space""" for ii in xrange(len(set_id_cal)): get_M.get_M(ns_cal[ii], set_id_cal[ii]) for ii in xrange(len(set_id_val)): get_M.get_M(ns_val[ii], set_id_val[ii]) """Compute the periodic statistics for the microstructures""" for ii in xrange(len(set_id_cal)): corr.correlate(ns_cal[ii], set_id_cal[ii]) for ii in xrange(len(set_id_val)): corr.correlate(ns_val[ii], set_id_val[ii]) """Perform PCA on correlations"""
# if tensorID == 1, we read the strain tensor # if tensorID == 2, we read the plastic strain tensor tensor_ID = 1 # # Gather data from calibration vtk files # vtk_r.read_euler(el, ns_cal, set_id_cal, step, dir_cal, wrt_file, 1) # for comp in compl: # vtk_r.read_meas(el, ns_cal, set_id_cal, step, comp, tensor_ID, dir_cal, # wrt_file) # Gather data from validation vtk files vtk_r.read_euler(el, ns_val, set_id_val, step, dir_val, wrt_file, 1) # for comp in compl: # vtk_r.read_meas(el, ns_val, set_id_val, step, comp, tensor_ID, dir_val, # wrt_file) # read grain IDs from the VTK files vtk_r.read_scalar(el, ns_val, set_id_val, step, dir_val, wrt_file) # # Convert the orientations from the calibration datasets from bunge euler # # angles to GSH coefficients # gsh.euler_to_gsh(el, H, ns_cal, set_id_cal, step, wrt_file) # Convert the orientations from the validation datasets from bunge euler # angles to GSH coefficients gsh.euler_to_gsh(el, H, ns_val, set_id_val, step, wrt_file)
tensor_ID = 0 # ## Gather data from calibration vtk files # # dir_cal = 'cal' # # vtk.read_euler(ns_cal,set_id_cal,step,vtk_filename,dir_cal, wrt_file, 0) # # for comp in compl: # vtk.read_meas(ns_cal,set_id_cal,step,comp,vtk_filename,tensor_ID,dir_cal,wrt_file) ## Gather data from validation vtk files dir_val = 'val' vtk.read_euler(ns_val, set_id_val, step, vtk_filename, dir_val, wrt_file) for comp in compl: vtk.read_meas(ns_val, set_id_val, step, comp, vtk_filename, tensor_ID, dir_val, wrt_file, 1) # ## Convert the orientations from the calibration datasets from bunge euler angles # ## to GSH coefficients # gsh.euler_to_gsh(ns_cal,set_id_cal,step,wrt_file) ## Convert the orientations from the validation datasets from bunge euler angles ## to GSH coefficients gsh.euler_to_gsh(ns_val, set_id_val, step, wrt_file) # ## Generate the fftn of the calibration microstructure function # msf.micr_func(ns_cal,set_id_cal,step,wrt_file)
step = 1 n_pc_tot = 20 wrt_file = 'log_%s.txt' % (time.strftime("%Y-%m-%d_h%Hm%M")) f = h5py.File("spatial_stats.hdf5", 'w') f.close() """ The tensorID determines the type of tensor data read from the .vtk file if tensorID == 0, we read the stress tensor if tensorID == 1, we read the strain tensor if tensorID == 2, we read the plastic strain tensor """ """Gather data from vtk files""" for ii in xrange(len(set_id_cal)): vtk.read_euler(el, ns_cal[ii], set_id_cal[ii], step, dir_cal[ii], wrt_file, 1) vtk.read_euler(el, ns_val[ii], set_id_val[ii], step, dir_val[ii], wrt_file, 1) """get the data for the linkage""" f = h5py.File("responses.hdf5", 'w') f.close() f = h5py.File("responses_for_linkage.hdf5", 'w') f.close() for ii in xrange(len(set_id_cal)): # gr.Eeff(el, ns_cal[ii], set_id_cal[ii], step, dir_cal[ii], wrt_file) # gr.Eeff(el, ns_val[ii], set_id_val[ii], step, dir_val[ii], wrt_file) gr.FIP(el, ns_cal[ii], set_id_cal[ii], step, dir_cal[ii], wrt_file) gr.FIP(el, ns_val[ii], set_id_val[ii], step, dir_val[ii], wrt_file) """Compute GSH coefficients to create microstructure function in real and
C = const() ns_cal = C['ns_cal'] set_id_cal = C['set_id_cal'] dir_cal = C['dir_cal'] ns_val = C['ns_val'] set_id_val = C['set_id_val'] dir_val = C['dir_val'] f = h5py.File("spatial.hdf5", 'w') f.close() """Gather data from vtk files""" vr.read_euler(ns_cal, set_id_cal, dir_cal, 1) vr.read_euler(ns_val, set_id_val, dir_val, 1) """get the data for the linkage""" f = h5py.File("responses.hdf5", 'w') f.close() vr.read_fip(ns_cal, set_id_cal, dir_cal) vr.read_fip(ns_cal, set_id_val, dir_val) """Compute GSH coefficients to create microstructure function in real and fourier space""" gm.get_mf(ns_cal, set_id_cal) gm.get_mf(ns_val, set_id_val)
compl = ['11', '22', '33', '12', '13', '23'] # compl = ['11'] for step in xrange(1, 2): wrt_file = 'log_step%s_%s.txt' % (step, time.strftime("%Y-%m-%d_h%Hm%M")) """ The tensorID determines the type of tensor data read from the .vtk file if tensorID == 0, we read the stress tensor if tensorID == 1, we read the strain tensor if tensorID == 2, we read the plastic strain tensor """ """Gather data from validation vtk files""" tensorID = 1 vtk.read_euler(el, ns_val, set_id_val, step, dir_val, wrt_file, tensorID) for comp in compl: vtk.read_meas(el, ns_val, set_id_val, step, comp, tensorID, dir_val, wrt_file) """get the plastic strain tensor from CPFEM""" for comp in compl: vtk.read_meas(el, ns_val, set_id_val, step, comp, 2, dir_val, wrt_file) # """Convert the orientations from the calibration datasets from bunge euler # angles to GSH coefficients""" # gsh.euler_to_gsh(el, H, ns_val, set_id_val, step, wrt_file) # """Perform the validation""" # for comp in compl: # validation.validation(el, H, ns_cal, ns_val, set_id_cal, set_id_val,
el = 21 step = 6 wrt_file = 'log_%s.txt' % (time.strftime("%Y-%m-%d_h%Hm%M")) f = h5py.File("spatial_stats.hdf5", 'w') f.close() """ The tensorID determines the type of tensor data read from the .vtk file if tensorID == 0, we read the stress tensor if tensorID == 1, we read the strain tensor if tensorID == 2, we read the plastic strain tensor """ """Gather data from calibration vtk files""" for ii in xrange(len(set_id_D3D)): vtk.read_euler(el, ns_D3D[ii], set_id_D3D[ii], step, dir_D3D[ii], wrt_file, 1) # gen.delta(el, ns_cal[1], set_id_cal[1], step, wrt_file) # vfrac = np.array([.1, .05]) # gen.inclusion(el, ns_cal[2], set_id_cal[2], step, wrt_file, vfrac) # gen.bicrystal(el, ns_cal[0], set_id_cal[0], step, wrt_file) """Compute GSH coefficients to create microstructure function in real and fourier space""" for ii in xrange(len(set_id_cal)): get_M.get_M(el, H, ns_cal[ii], set_id_cal[ii], step, wrt_file) """Compute the periodic statistics for the microstructures""" for ii in xrange(len(set_id_cal)): corr.correlate(el, ns_cal[ii], H, set_id_cal[ii], step, wrt_file) """Perform PCA on correlations"""