CC_3D_Y_1.rar_Abaqus 周期边界_abaqus 脚本_abaqus 边界条件_abaqus周期性_python

# Subroutine for Periodic Boundary conditions from abaqus import * from abaqusConstants import * #from Numeric import * from math import * #from numpy import * import interaction a = mdb.models['Model-1'].rootAssembly p1 = mdb.models['Model-1'].rootAssembly.sets['Y_face01'] n1 = p1.nodes p2 = mdb.models['Model-1'].rootAssembly.sets['Y_face02'] n2 = p2.nodes #tol=0.003 tol=0.018 ii=0 for i1 in n1 : for i2 in n2 : b1 = i1.coordinates[0] b2 = i1.coordinates[2] c1 = i2.coordinates[0] c2 = i2.coordinates[2] d1 = (c1 - b1) d2 = (c2 - b2) d3 = sqrt(pow(d1,2)+pow(d2,2)) #d2 = sqrt(pow(d1,2)) if (d3<tol) : label1=i1.label label2=i2.label nodeset=mdb.models['Model-1'].rootAssembly.SetFromNodeLabels(name='Y_PF1_setN1_%03d' %(ii),nodeLabels=(('toron_x-2', (label1,),),),) nodeset=mdb.models['Model-1'].rootAssembly.SetFromNodeLabels(name='Y_PF1_setN2_%03d' %(ii),nodeLabels=(('toron_x-1', (label2,),),),) mdb.models['Model-1'].Equation(name='Y_PF1_EQN1_%03d' %(ii), terms=((1.0, 'Y_PF1_setN1_%03d' %(ii) , 1), (-1.0, 'Y_PF1_setN2_%03d' %(ii) , 1))) mdb.models['Model-1'].Equation(name='Y_PF1_EQN2_%03d' %(ii), terms=((1.0, 'Y_PF1_setN1_%03d' %(ii) , 2), (-1.0, 'Y_PF1_setN2_%03d' %(ii) , 2))) mdb.models['Model-1'].Equation(name='Y_PF1_EQN3_%03d' %(ii), terms=((1.0, 'Y_PF1_setN1_%03d' %(ii) , 3), (-1.0, 'Y_PF1_setN2_%03d' %(ii) , 3))) #mdb.models['WovenCell2D'].Equation(name='EQN3_%03d' %(ii), terms=((1.0, 'setN1_%03d' %(ii) , 3), (-1.0, 'setN2_%03d' %(ii) , 3), (-1.0, 'RP1', 3))) ii=ii+1 break print ii