'estriau1', 'estriau2', 'estrifu1', 'estupper', 'isresau1', 'isresfu1', 'isrupper', 'istupper', 'istriau1', 'istrifu1' ] # element sets of triangular elements on the lower surface # outer_edge_node_nums=[2,1] list_of_tri_lower_elementsets = ['estril_tri', 'istril_tri', 'isresl_tri'] # element sets of triangular elements on the upper surface # outer_edge_node_nums=[3,2] list_of_tri_upper_elementsets = ['estriu_tri', 'istriu_tri', 'isresu_tri'] # import the initial grid object fmt_grid = 'sandia_blade/' + stn_str + '/mesh_' + stn_str + '.abq' g = au.AbaqusGrid(fmt_grid, debug_flag=True, soft_warning=False, auto_parse=True) # manually assign triangular elements into new element sets for num in ES_triax_upper_tri_elem_num: g.list_of_elements[num - 1].element_set = 'estriu_tri' for num in ES_triax_lower_tri_elem_num: g.list_of_elements[num - 1].element_set = 'estril_tri' for num in IS_triax_upper_tri_elem_num: g.list_of_elements[num - 1].element_set = 'istriu_tri' for num in IS_triax_lower_tri_elem_num: g.list_of_elements[num - 1].element_set = 'istril_tri' for num in IS_resin_upper_tri_elem_num: g.list_of_elements[num - 1].element_set = 'isresu_tri' for num in IS_resin_lower_tri_elem_num: g.list_of_elements[num - 1].element_set = 'isresl_tri'
'estrteu1', 'estrteu2', 'estrteu3', 'estrteu4', 'esglteu1', 'esglteu2', 'esglteu3', 'esglteu4', 'esglasw2', 'estrasw2', 'rbtrasw2', 'scupper', 'rbtriscu', 'estriscu', 'is2rsscu', 'esgelscu', 'is2trscu' ] # element sets of triangular elements on the lower surface # outer_edge_node_nums=[2,1] list_of_tri_lower_elementsets = ['is3rtel2_tri', 'tefoaml3_tri'] # element sets of triangular elements on the upper surface # outer_edge_node_nums=[3,2] list_of_tri_upper_elementsets = ['is3rteu2_tri', 'tefoamu3_tri'] # import the initial grid object fmt_grid = 'biplane_blade/' + stn_str + '/mesh_' + stn_str + '.abq' g = au.AbaqusGrid(fmt_grid, debug_flag=True) # manually assign two triangular elements into new element sets g.list_of_elements[IS3_resin_u2_tri_elem_num - 1].element_set = 'is3rteu2_tri' g.list_of_elements[IS3_resin_l2_tri_elem_num - 1].element_set = 'is3rtel2_tri' g.list_of_elements[TE_reinf_foam_u3_tri_elem_num - 1].element_set = 'tefoamu3_tri' g.list_of_elements[TE_reinf_foam_l3_tri_elem_num - 1].element_set = 'tefoaml3_tri' # update the grid object with all the layer plane angles for elem in g.list_of_elements: if elem.element_set in list_of_LE_elementsets: elem.calculate_layer_plane_angle(outer_edge_node_nums=[1, 4], inner_edge_node_nums=[2, 3]) elif elem.element_set in list_of_TE_elementsets: