dim = 2 if dim == 3: filename = data_dir + '/meshes/3d/special/cube_sphere.mesh' else: filename = data_dir + '/meshes/2d/special/circle_in_square.mesh' output_dir = incwd('output/band_gaps_rigid') # Rigid inclusion diameter. yr_diameter = 0.125 # aluminium, in 1e+10 Pa D_m = get_pars(dim, 5.898, 2.681) density_m = 0.2799 # in 1e4 kg/m3 # epoxy, in 1e+10 Pa D_c = get_pars(dim, 0.1798, 0.148) density_c = 0.1142 # in 1e4 kg/m3 # lead, in 1e+10 Pa, does not matter D_r = get_pars(dim, 4.074, 0.5556) density_r = 1.1340 # in 1e4 kg/m3 mat_pars = Coefficients(D_m=D_m, density_m=density_m, D_c=D_c, density_c=density_c, D_r=D_r,
dim = 2 if dim == 3: filename = data_dir + '/meshes/3d/special/cube_sphere.mesh' else: filename = data_dir + '/meshes/2d/special/circle_in_square.mesh' output_dir = incwd('output/band_gaps_rigid') # Rigid inclusion diameter. yr_diameter = 0.125 # aluminium, in 1e+10 Pa D_m = get_pars(dim, 5.898, 2.681) density_m = 0.2799 # in 1e4 kg/m3 # epoxy, in 1e+10 Pa D_c = get_pars(dim, 0.1798, 0.148) density_c = 0.1142 # in 1e4 kg/m3 # lead, in 1e+10 Pa, does not matter D_r = get_pars(dim, 4.074, 0.5556) density_r = 1.1340 # in 1e4 kg/m3 mat_pars = Coefficients(D_m=D_m, density_m=density_m, D_c=D_c, density_c=density_c, D_r=D_r, density_r=density_r) region_selects = Struct(matrix='cells of group 1',
from sfepy.base.base import Struct from sfepy.base.ioutils import InDir from sfepy.homogenization.coefficients import Coefficients from band_gaps_conf import BandGapsConf, get_pars, clip_sqrt, normalize clip_sqrt, normalize # Make pyflakes happy... incwd = InDir(__file__) filename = data_dir + '/meshes/2d/special/circle_in_square.mesh' output_dir = incwd('output/band_gaps') # aluminium, in 1e+10 Pa D_m = get_pars(2, 5.898, 2.681) density_m = 0.2799 # in 1e4 kg/m3 # epoxy, in 1e+10 Pa D_c = get_pars(2, 0.1798, 0.148) density_c = 0.1142 # in 1e4 kg/m3 mat_pars = Coefficients(D_m=D_m, density_m=density_m, D_c=D_c, density_c=density_c) region_selects = Struct(matrix='cells of group 1', inclusion='cells of group 2') corrs_save_names = {'evp': 'evp', 'corrs_rs': 'corrs_rs'}
from sfepy.base.base import Struct from sfepy.base.ioutils import InDir from sfepy.homogenization.coefficients import Coefficients from band_gaps_conf import BandGapsConf, get_pars, clip_sqrt, normalize clip_sqrt, normalize # Make pyflakes happy... incwd = InDir(__file__) filename = data_dir + '/meshes/2d/special/circle_in_square.mesh' output_dir = incwd('output/band_gaps') # aluminium, in 1e+10 Pa D_m = get_pars(2, 5.898, 2.681) density_m = 0.2799 # in 1e4 kg/m3 # epoxy, in 1e+10 Pa D_c = get_pars(2, 0.1798, 0.148) density_c = 0.1142 # in 1e4 kg/m3 mat_pars = Coefficients(D_m=D_m, density_m=density_m, D_c=D_c, density_c=density_c) region_selects = Struct(matrix='cells of group 1', inclusion='cells of group 2') corrs_save_names = {'evp' : 'evp', 'corrs_rs' : 'corrs_rs'} options = {