Python二次开发——findAt方法的正确食用方法

#!/user/bin/python# -* - coding:UTF-8 -*-# 导入脚本中使用的各个模块from abaqus import *from abaqusConstants import *from caeModules import *#创建新模型 Model-1myModel = mdb.models['Model-1']# 创建新的视口myViewport=session.Viewport(name='Region syntax',                origin=(20, 20), width=200, height=100)# 创建一个草图，并绘制两个矩形 mySketch = myModel.ConstrainedSketch(name='Sketch A',    sheetSize=200.0)mySketch.rectangle(point1=(-40.0, 30.0),    point2=(-10.0, 0.0))mySketch.rectangle(point1=(10.0, 30.0),    point2=(40.0, 0.0))# 对上面创建的两个矩形进行拉伸操作，生成三维部件door = myModel.Part(name='Door',    dimensionality=THREE_D, type=DEFORMABLE_BODY)door.BaseSolidExtrude(sketch=mySketch, depth=20.0)# 创建部件实例myAssembly = myModel.rootAssemblydoorInstance = myAssembly.Instance(name='Door-1',    part=door, dependent=OFF)# 选择两个顶点pillarVertices = doorInstance.vertices.findAt(    ((-40,30,0),), ((40,0,0),) )# 创建静力分析步（static）myModel.StaticStep(name='impact',    previous='Initial', initialInc=1, timePeriod=1)# 在选择的顶点上施加集中力myPillarLoad = myModel.ConcentratedForce(    name='pillarForce', createStepName='impact',    region=(pillarVertices,), cf1=12.50E4)# 选择两个面topFace = doorInstance.faces.findAt(((-25,30,10),))bottomFace = doorInstance.faces.findAt(((-25,0,10),))# 在选择的面上施加压力（pressure）# 对于同一部件实例的相同类型的实体，可以使用 + 号myFenderLoad = myModel.Pressure(    name='pillarPressure', createStepName='impact',    region=((topFace+bottomFace, SIDE1),),    magnitude=10E4)# 在同一部件实例上选择两条边myEdge1 = doorInstance.edges.findAt(((10,15,20),))myEdge2 = doorInstance.edges.findAt(((10,15,0),))# 对一个面、两条边和两个顶点施加边界条件myDisplacementBc= myModel.DisplacementBC(    name='xBC', createStepName='impact',    region=(pillarVertices, myEdge1+myEdge2,    topFace), u1=5.0)# 使用面上的任意点选择两个面faceRegion = doorInstance.faces.findAt(    ((-30,15,20), ), ((30,15,20),))# 创建包含两个面（face）的表面（surface）mySurface = myModel.rootAssembly.Surface(    name='exterior', side1Faces=faceRegion)# 使用这个表面来创建弹性地基（elastic foundation）myFoundation = myModel.ElasticFoundation(     name='elasticFloor', createStepName='Initial',     surface=mySurface, stiffness=1500)# 显示施加荷载和边界条件后的装配件myViewport.setValues(displayedObject=myAssembly)myViewport.assemblyDisplay.setValues(step='impact',     loads=ON, bcs=ON, predefinedFields=ON)

a1 = mdb.models['Model-1'].rootAssemblyv1 = a1.instances['Part-1-1'].verticesverts1 = v1.findAt(((5.0,5.0,200.0),))region = a1.Set(vertices=verts1, name='Set-2')

import mathcood_x=5.0*math.sin(45.0/180.0*math.pi)cood_y=5.0*math.cos(45.0/180.0*math.pi)side1Edges3 =s1.findAt(((cood_x,cood_y,0.0),))region3 = a.Surface(side1Edges=side1Edges3, name='Surf-1')mdb.models['Model-1'].Pressure(name='Load-1', createStepName='Step-1',   region=region3, distributionType=UNIFORM, field='', magnitude=-pp,   amplitude=UNSET)

a1 = mdb.models['Model-1'].rootAssemblyv1 = a1.instances['Part-1-1'].verticesverts1 = v1.findAt(((5.0,5.0,200.0),),((5.0,-5.0,200.0),),((-5.0,-5.0,200.0),),((-5.0,5.0,200.0),))region = a1.Set(vertices=verts1, name='Set-2')mdb.models['Model-1'].ConcentratedForce(name='Load-1', createStepName='Step-1',     region=region, cf3=1000.0, distributionType=UNIFORM, field='',     localCsys=None)