from sfepy.base.ioutils import InDir from sfepy.homogenization.coefficients import Coefficients from examples.phononic.band_gaps_conf import (BandGapsConf, get_pars, clip, clip_sqrt) clip, clip_sqrt # Make pyflakes happy... incwd = InDir(__file__) filename = data_dir + '/meshes/2d/special/circle_in_square.mesh' output_dir = incwd('output/band_gaps') # aluminium, SI units D_m = get_pars(2, 5.898e10, 2.681e10) density_m = 2799.0 # epoxy, SI units D_c = get_pars(2, 1.798e9, 1.48e9) density_c = 1142.0 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.ioutils import InDir from sfepy.homogenization.coefficients import Coefficients from examples.phononic.band_gaps_conf import (BandGapsConf, get_pars, clip, clip_sqrt) clip, clip_sqrt # Make pyflakes happy... incwd = InDir(__file__) filename = data_dir + '/meshes/2d/special/circle_in_square.mesh' output_dir = incwd('output/band_gaps') # aluminium, SI units D_m = get_pars(2, 5.898e10, 2.681e10) density_m = 2799.0 # epoxy, SI units D_c = get_pars(2, 1.798e9, 1.48e9) density_c = 1142.0 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 = {
from sfepy.base.ioutils import InDir from sfepy.homogenization.coefficients import Coefficients from examples.phononic.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 = {
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',
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, SI units D_m = get_pars(2, 5.898e10, 2.681e10) density_m = 2799.0 # epoxy, SI units D_c = get_pars(2, 1.798e9, 1.48e9) density_c = 1142.0 # lead, SI units, does not matter D_r = get_pars(dim, 4.074e10, 5.556e9) density_r = 11340.0 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',
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.ioutils import InDir from sfepy.homogenization.coefficients import Coefficients from examples.phononic.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'}
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, SI units D_m = get_pars(2, 5.898e10, 2.681e10) density_m = 2799.0 # epoxy, SI units D_c = get_pars(2, 1.798e9, 1.48e9) density_c = 1142.0 # lead, SI units, does not matter D_r = get_pars(dim, 4.074e10, 5.556e9) density_r = 11340.0 mat_pars = Coefficients(D_m=D_m, density_m=density_m, D_c=D_c, density_c=density_c, D_r=D_r,