def load_model(self, model_path=None): global nodes, Elements, domains, opt_domains, en_all, plane_strain, plane_stress, axisymmetry, file_name, domain_shells, domain_volumes, elm_states model_path = model_path if (model_path is not None) else file_name [ nodes, Elements, domains, opt_domains, en_all, plane_strain, plane_stress, axisymmetry ] = beso_lib.import_inp(model_path, domains_from_config, domain_optimized, shells_as_composite) domain_shells = {} domain_volumes = {} for dn in domains_from_config: # distinguishing shell elements and volume elements domain_shells[dn] = set(domains[dn]).intersection( list(Elements.tria3.keys()) + list(Elements.tria6.keys()) + list(Elements.quad4.keys()) + list(Elements.quad8.keys())) domain_volumes[dn] = set(domains[dn]).intersection( list(Elements.tetra4.keys()) + list(Elements.tetra10.keys()) + list(Elements.hexa8.keys()) + list(Elements.hexa20.keys()) + list(Elements.penta6.keys()) + list(Elements.penta15.keys())) # initial element state elm_states = {} if continue_from[-4:] == ".frd": elm_states = beso_lib.import_frd_state(continue_from, elm_states, number_of_states, file_name) elif continue_from[-4:] == ".inp": elm_states = beso_lib.import_inp_state(continue_from, elm_states, number_of_states, file_name) else: for dn in domains_from_config: for en in domains[dn]: elm_states[en] = len( domain_density[dn]) - 1 # set to highest state
elm_states = {} if isinstance(continue_from, int): for dn in domains_from_config: if (len(domain_density[dn]) - 1) < continue_from: sn = len(domain_density[dn]) - 1 msg = "\nINFO: elements from the domain " + dn + " were set to the highest state.\n" beso_lib.write_to_log(file_name, msg) print(msg) else: sn = continue_from for en in domains[dn]: elm_states[en] = sn elif continue_from[-4:] == ".frd": elm_states = beso_lib.import_frd_state(continue_from, elm_states, number_of_states, file_name) elif continue_from[-4:] == ".inp": elm_states = beso_lib.import_inp_state(continue_from, elm_states, number_of_states, file_name) elif continue_from[-4:] == ".csv": elm_states = beso_lib.import_csv_state(continue_from, elm_states, file_name) else: for dn in domains_from_config: for en in domains[dn]: elm_states[en] = len(domain_density[dn]) - 1 # set to highest state # computing volume or area, and centre of gravity of each element [cg, cg_min, cg_max, volume_elm, area_elm] = beso_lib.elm_volume_cg(file_name, nodes, Elements) mass = [0.0] mass_full = 0 # sum from initial states TODO make it independent on starting elm_states? for dn in domains_from_config: if domain_optimized[dn] is True: for en in domain_shells[dn]:
elm_states = {} if isinstance(continue_from, int): for dn in domains_from_config: if (len(domain_density[dn]) - 1) < continue_from: sn = len(domain_density[dn]) - 1 msg = "\nINFO: elements from the domain " + dn + " were set to the highest state.\n" beso_lib.write_to_log(file_name, msg) print(msg) else: sn = continue_from for en in domains[dn]: elm_states[en] = sn elif continue_from[-4:] == ".frd": elm_states = beso_lib.import_frd_state(continue_from, elm_states, number_of_states, file_name) elif continue_from[-4:] == ".inp": elm_states = beso_lib.import_inp_state(continue_from, elm_states, number_of_states, file_name) elif continue_from[-4:] == ".csv": elm_states = beso_lib.import_csv_state(continue_from, elm_states, file_name) else: for dn in domains_from_config: for en in domains[dn]: elm_states[en] = len(domain_density[dn]) - 1 # set to highest state # computing volume or area, and centre of gravity of each element [cg, cg_min, cg_max, volume_elm, area_elm] = beso_lib.elm_volume_cg(file_name, nodes, Elements) mass = [0.0] mass_full = 0 # sum from initial states TODO make it independent on starting elm_states? for dn in domains_from_config: if domain_optimized[dn] is True: for en in domain_shells[dn]: