'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:
