def create_rbe3s_between_close_nodes(bdf_filename, bdf_filename_out, tol, renumber_nodes=False, neq_max=4, xref=True, node_set=None, size=8, is_double=False, debug=True): """ Creates semi-rigid RBE3 elements between two nodes within tolerance. Parameters ---------- bdf_filename : str / BDF str : bdf file path BDF : a BDF model that is fully valid (see xref) bdf_filename_out : str a bdf_filename to write tol : float the spherical tolerance renumber_nodes : bool should the nodes be renumbered (default=False) neq_max : int the number of "close" points (default=4) xref bool: bool does the model need to be cross_referenced (default=True; only applies to model option) node_set : List[int] / (n, ) ndarray the list/array of nodes to consider (not supported with renumber_nodes=True) size : int; {8, 16}; default=8 the bdf write precision is_double : bool; default=False the field precision to write Returns ------- model : BDF() The BDF model corresponding to bdf_filename_out .. warning:: I doubt SPOINTs/EPOINTs work correctly .. warning:: xref not fully implemented (assumes cid=0) .. todo:: node_set stil does work on the all the nodes in the big kdtree loop, which is very inefficient """ nodes_xyz, model, nids, inew = _eq_nodes_setup( bdf_filename, tol, renumber_nodes=renumber_nodes, xref=xref, node_set=node_set, debug=debug) eid = 1 if len(model.rigid_elements): eid = max(model.rigid_elements) + 1 ieq, slots = _eq_nodes_build_tree(nodes_xyz, nids, tol, inew=inew, node_set=node_set, neq_max=neq_max)[1:] nid_pairs = _eq_nodes_find_pairs(nids, slots, ieq, node_set=node_set) for (nid1, nid2) in nid_pairs: node1 = model.nodes[nid1] node2 = model.nodes[nid2] # TODO: doesn't use get position... distance = norm(node1.xyz - node2.xyz) if distance > tol: continue refgrid = nid1 refc = '123456' weights = [1.0] comps = ['123456'] gijs = [[nid2]] rbe3 = RBE3(eid, refgrid, refc, weights, comps, gijs, comment='') model.rigid_elements[eid] = rbe3 eid += 1 #_eq_nodes_final(nid_pairs, model, tol, node_set=node_set) if bdf_filename_out is not None: model.write_bdf(bdf_filename_out, size=size, is_double=is_double) return model
def create_rbe3s_between_close_nodes(bdf_filename, bdf_filename_out, tol, renumber_nodes=False, neq_max=4, xref=True, node_set=None, size=8, is_double=False, debug=True): """ Creates semi-rigid RBE3 elements between two nodes within tolerance. Parameters ---------- bdf_filename : str / BDF str : bdf file path BDF : a BDF model that is fully valid (see xref) bdf_filename_out : str a bdf_filename to write tol : float the spherical tolerance renumber_nodes : bool should the nodes be renumbered (default=False) neq_max : int the number of "close" points (default=4) xref bool: bool does the model need to be cross_referenced (default=True; only applies to model option) node_set : List[int] / (n, ) ndarray the list/array of nodes to consider (not supported with renumber_nodes=True) size : int; {8, 16}; default=8 the bdf write precision is_double : bool; default=False the field precision to write Returns ------- model : BDF() The BDF model corresponding to bdf_filename_out .. warning:: I doubt SPOINTs/EPOINTs work correctly .. warning:: xref not fully implemented (assumes cid=0) .. todo:: node_set stil does work on the all the nodes in the big kdtree loop, which is very inefficient """ nodes_xyz, model, nids, inew = _eq_nodes_setup( bdf_filename, tol, renumber_nodes=renumber_nodes, xref=xref, node_set=node_set, debug=debug) eid = 1 if len(model.rigid_elements): eid = max(model.rigid_elements) + 1 ieq, slots = _eq_nodes_build_tree(nodes_xyz, nids, tol, inew=inew, node_set=node_set, neq_max=neq_max)[1:] nid_pairs = _eq_nodes_find_pairs(nids, slots, ieq, node_set=node_set) for (nid1, nid2) in nid_pairs: node1 = model.nodes[nid1] node2 = model.nodes[nid2] # TODO: doesn't use get position... distance = norm(node1.xyz - node2.xyz) if distance > tol: continue refgrid = nid1 refc = '123456' weights = [1.0] comps = ['123456'] gijs = [[nid2]] rbe3 = RBE3(eid, refgrid, refc, weights, comps, gijs, comment='') model.rigid_elements[eid] = rbe3 eid += 1 #_eq_nodes_final(nid_pairs, model, tol, node_set=node_set) if bdf_filename_out is not None: model.write_bdf(bdf_filename_out, size=size, is_double=is_double) return model