告别重复劳动!Abaqus Python 自动化建模:手把手教你实现“分层建模”与“节点自动融合”内置函数 mergeNodes

from abaqus import *from abaqusConstants import *import materialimport sectionimport regionToolsetimport displayGroupMdbToolset as dgmimport partimport assemblyimport stepimport interactionimport loadimport meshimport optimizationimport jobimport sketchimport visualizationimport xyPlotimport displayGroupOdbToolset as dgoimport connectorBehaviorimport sysimport osL = 100.0H = 20.0H_split = H / 2.0nx = 20ny1 = 2ny2 = 2dx = L / nxdy1 = H_split / ny1dy2 = H_split / ny2total_force = 1000.0model_name = 'Automated_Plane_Stress_Split'if model_name in mdb.models:   del mdb.models[model_name]my_model = mdb.Model(name=model_name)my_part = my_model.Part(   name='Plate',   dimensionality=TWO_D_PLANAR,   type=DEFORMABLE_BODY)nodes_region1 = []for j in range(ny1 + 1):   for i in range(nx + 1):       x = i * dx       y = j * dy1       n = my_part.Node(coordinates=(x, y, 0.0))       nodes_region1.append(n)nodes_region2 = []for j in range(ny2 + 1):   for i in range(nx + 1):       x = i * dx       y = H_split + j * dy2       n = my_part.Node(coordinates=(x, y, 0.0))       nodes_region2.append(n)all_part_nodes = my_part.nodesnode_dict = {}tolerance = 1e-6for node in all_part_nodes:   coord_key = (       round(node.coordinates[0] / tolerance) * tolerance,       round(node.coordinates[1] / tolerance) * tolerance,       round(node.coordinates[2] / tolerance) * tolerance   )   if coord_key not in node_dict:       node_dict[coord_key] = node.labelnode_merge_map = {}for node in all_part_nodes:   coord_key = (       round(node.coordinates[0] / tolerance) * tolerance,       round(node.coordinates[1] / tolerance) * tolerance,       round(node.coordinates[2] / tolerance) * tolerance   )   node_merge_map[node.label] = node_dict[coord_key]duplicate_count = len(all_part_nodes) - len(node_dict)elements_region1 = []for j in range(ny1):   for i in range(nx):       n1_local = j * (nx + 1) + i       n2_local = j * (nx + 1) + (i + 1)       n3_local = (j + 1) * (nx + 1) + (i + 1)       n4_local = (j + 1) * (nx + 1) + i       node1 = nodes_region1[n1_local]       node2 = nodes_region1[n2_local]       node3 = nodes_region1[n3_local]       node4 = nodes_region1[n4_local]       elem = my_part.Element(           nodes=(node1, node2, node3, node4),           elemShape=QUAD4       )       elements_region1.append(elem)elements_region2 = []for j in range(ny2):   for i in range(nx):       n1_local = j * (nx + 1) + i       n2_local = j * (nx + 1) + (i + 1)       n3_local = (j + 1) * (nx + 1) + (i + 1)       n4_local = (j + 1) * (nx + 1) + i       node1 = nodes_region2[n1_local]       node2 = nodes_region2[n2_local]       node3 = nodes_region2[n3_local]       node4 = nodes_region2[n4_local]       elem = my_part.Element(           nodes=(node1, node2, node3, node4),           elemShape=QUAD4       )       elements_region2.append(elem)all_node_labels = [n.label for n in my_part.nodes]my_part.Set(   nodes=my_part.nodes.sequenceFromLabels(all_node_labels),   name='All-Nodes')my_part.mergeNodes(   nodes=my_part.sets['All-Nodes'],   tolerance=tolerance)mat = my_model.Material(name='Steel')mat.Elastic(table=((210000.0, 0.3), ))my_model.HomogeneousSolidSection(   name='Sec',   material='Steel',   thickness=1.0)p_elems = my_part.elementsmy_part.Set(elements=p_elems, name='All-E')my_part.SectionAssignment(   region=my_part.sets['All-E'],   sectionName='Sec')my_part.setElementType(   regions=my_part.sets['All-E'],   elemTypes=(mesh.ElemType(elemCode=CPS4R),))all_nodes_final = my_part.nodesleft_labels = [n.label for n in all_nodes_final if abs(n.coordinates[0] - 0.0) < tolerance]my_part.Set(   nodes=my_part.nodes.sequenceFromLabels(left_labels),   name='Set-Left')right_labels = [n.label for n in all_nodes_final if abs(n.coordinates[0] - L) < tolerance]my_part.Set(   nodes=my_part.nodes.sequenceFromLabels(right_labels),   name='Set-Right')my_asm = my_model.rootAssemblyinst = my_asm.Instance(name='Inst', part=my_part, dependent=ON)my_model.StaticStep(name='Step-1', previous='Initial')my_model.EncastreBC(   name='Fix',   createStepName='Initial',   region=inst.sets['Set-Left'])f_per_node = total_force / len(right_labels)my_model.ConcentratedForce(   name='Pull',   createStepName='Step-1',   region=inst.sets['Set-Right'],   cf1=f_per_node)job_name = 'Analysis_Task_Split'if job_name in mdb.jobs:   del mdb.jobs[job_name]my_job = mdb.Job(   name=job_name,   model=model_name,   type=ANALYSIS,   description='Automated Plane Stress Analysis - Split Domain')my_job.submit()my_job.waitForCompletion()odb_path = job_name + '.odb'my_odb = session.openOdb(name=odb_path)vp = session.viewports['Viewport: 1']vp.setValues(displayedObject=my_odb)vp.odbDisplay.setPrimaryVariable(   variableLabel='S',   outputPosition=INTEGRATION_POINT,   refinement=(INVARIANT, 'Mises'))vp.odbDisplay.display.setValues(plotState=(CONTOURS_ON_DEF, ))vp.odbDisplay.commonOptions.setValues(deformationScaling=AUTO)