approx_p = '3_4_P0' order = 2 filename_mesh = data_dir + '/meshes/3d/cylinder.mesh' #! Regions #! ------- #! Whole domain 'Omega', left and right ends. regions = { 'Omega' : 'all', 'Left' : ('vertices in (x < 0.001)', 'facet'), 'Right' : ('vertices in (x > 0.099)', 'facet'), } #! Materials #! --------- #! The linear elastic material model is used. materials = { 'solid' : ({'D' : stiffness_from_youngpoisson_mixed(dim, 0.7e9, 0.4), 'gamma' : 1.0/bulk_from_youngpoisson(0.7e9, 0.4)},), } #! Fields #! ------ #! A field is used to define the approximation on a (sub)domain fields = { 'displacement': ('real', 'vector', 'Omega', 1), 'pressure' : ('real', 'scalar', 'Omega', 0), } #! Integrals #! --------- #! Define the integral type Volume/Surface and quadrature rule. integrals = { 'i' : order, }
filename_mesh = data_dir + '/meshes/3d/matrix_fiber.mesh' region_lbn = (0, 0, 0) region_rtf = (1, 1, 1) #! Regions #! ------- #! Regions, edges, ... regions = { 'Y' : 'all', 'Ym' : 'cells of group 1', 'Yc' : 'cells of group 2', } regions.update( define_box_regions( dim, region_lbn, region_rtf ) ) #! Materials #! --------- materials = { 'mat' : ({'D' : {'Ym': stiffness_from_youngpoisson_mixed(dim, 7.0e9, 0.4), 'Yc': stiffness_from_youngpoisson_mixed(dim, 70.0e9, 0.2)}, 'gamma': {'Ym': 1.0/bulk_from_youngpoisson(7.0e9, 0.4), 'Yc': 1.0/bulk_from_youngpoisson(70.0e9, 0.2)}},), } #! Fields #! ------ #! Scalar field for corrector basis functions. fields = { 'corrector_u' : ('real', dim, 'Y', 1), 'corrector_p' : ('real', 1, 'Y', 0), } #! Variables #! --------- #! Unknown and corresponding test variables. Parameter fields #! used for evaluation of homogenized coefficients.
dim = 3 filename_mesh = data_dir + '/meshes/3d/matrix_fiber.mesh' region_lbn = (0, 0, 0) region_rtf = (1, 1, 1) regions = { 'Y': 'all', 'Ym': 'cells of group 1', 'Yc': 'cells of group 2', } regions.update(define_box_regions(dim, region_lbn, region_rtf)) materials = { 'mat': ({'D': {'Ym': stiffness_from_youngpoisson_mixed(dim, 7.0e9, 0.4), 'Yc': stiffness_from_youngpoisson_mixed(dim, 70.0e9, 0.2)}, 'gamma': {'Ym': 1.0/bulk_from_youngpoisson(7.0e9, 0.4), 'Yc': 1.0/bulk_from_youngpoisson(70.0e9, 0.2)}},), } fields = { 'corrector_u': ('real', dim, 'Y', 1), 'corrector_p': ('real', 1, 'Y', 0), } variables = { 'u': ('unknown field', 'corrector_u'), 'v': ('test field', 'corrector_u', 'u'), 'p': ('unknown field', 'corrector_p'), 'q': ('test field', 'corrector_p', 'p'),
filename_mesh = data_dir + '/meshes/3d/matrix_fiber.mesh' region_lbn = (0, 0, 0) region_rtf = (1, 1, 1) #! Regions #! ------- #! Regions, edges, ... regions = { 'Y' : ('all', {}), 'Ym' : ('elements of group 1', {}), 'Yc' : ('elements of group 2', {}), } regions.update( define_box_regions( dim, region_lbn, region_rtf ) ) #! Materials #! --------- materials = { 'matrix' : ({'D' : stiffness_from_youngpoisson_mixed(dim, 0.7e9, 0.4), 'gamma' : 1.0/bulk_from_youngpoisson(0.7e9, 0.4)},), 'reinf' : ({'D' : stiffness_from_youngpoisson_mixed(dim, 70.0e9, 0.2), 'gamma' : 1.0/bulk_from_youngpoisson(70.0e9, 0.2)},), } #! Fields #! ------ #! Scalar field for corrector basis functions. fields = { 'corrector_u' : ('real', dim, 'Y', 1), 'corrector_p' : ('real', 1, 'Y', 0), } #! Variables #! --------- #! Unknown and corresponding test variables. Parameter fields #! used for evaluation of homogenized coefficients.
dim = 3 approx_u = "3_4_P1" approx_p = "3_4_P0" order = 2 filename_mesh = data_dir + "/meshes/3d/cylinder.mesh" #! Regions #! ------- #! Whole domain 'Omega', left and right ends. regions = {"Omega": ("all", {}), "Left": ("nodes in (x < 0.001)", {}), "Right": ("nodes in (x > 0.099)", {})} #! Materials #! --------- #! The linear elastic material model is used. materials = { "solid": ( {"D": stiffness_from_youngpoisson_mixed(dim, 0.7e9, 0.4), "gamma": 1.0 / bulk_from_youngpoisson(0.7e9, 0.4)}, ) } #! Fields #! ------ #! A field is used to define the approximation on a (sub)domain fields = {"displacement": ("real", "vector", "Omega", 1), "pressure": ("real", "scalar", "Omega", 0)} #! Integrals #! --------- #! Define the integral type Volume/Surface and quadrature rule. integrals = {"i1": ("v", order)} #! Variables #! --------- #! Define displacement and pressure fields and corresponding fields #! for test variables. variables = {