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,