simlab python二次开发2-一键生成轴瓦并设定节点号

1、节点坐标计算并建立
- 1.1、建坐标原点节点，并得到Model-1.gda
- 1.2、轴瓦节点计算并建立
2、由节点建面
- 2.1、由4个节点建面得到3个面单元Body
- 2.2、得到Bodies名称
- 2.3、根据==Bodies名称选面特征==（放入Group_1组）
- 2.4、镜像选择的面
3、旋转面得到体SLBodies
4、==simlab.getBodiesWithSubString==（）得到SLBodies名称
5、MergeBodies
6、平移拷贝轴瓦
7、根据节点坐标重定义节点号

1、节点坐标计算并建立

1.1、建坐标原点节点，并得到Model-1.gda

#***************************************************************
#SimLab Version 2022 
#Created at Mon Apr 22 19:05:48 2024
#***************************************************************
#For debugging this python script,Please comment(#) out the line "from hwx import simlab" and uncomment the line "import simlab"
#import simlab
from hwx import simlabUnitSystem=''' <UnitSystem UUID="3aca8564-4d38-4b0b-887c-6a542d4001c6"><SetCurrentDisplaySystem Name="MPA (mm t N s)"/></UnitSystem>''';
simlab.execute(UnitSystem);NodeByXYZ=''' <NodeByXYZ UUID="F200B5A2-D615-4d01-8DE2-25596B3B1EB8"><tag Value="-1"/><Name Value=""/><LocalCoordinateSystem Value="0"/><Position Value="0 mm,0 mm,0 mm"/><Node Value="1"/><UniqueNodeId Value="0"/><ModelName Value=""/></NodeByXYZ>''';
simlab.execute(NodeByXYZ);

1.2、轴瓦节点计算并建立

#***************************************************************
#SimLab Version 2022 
#Created at Sun Apr 21 20:29:55 2024
#***************************************************************
#For debugging this python script,Please comment(#) out the line "from hwx import simlab" and uncomment the line "import simlab"
#import simlab
from hwx import simlab
import numpy as np
import csvUnitSystem=''' <UnitSystem UUID="3aca8564-4d38-4b0b-887c-6a542d4001c6"><SetCurrentDisplaySystem Name="MPA (mm t N s)"/></UnitSystem>''';
simlab.execute(UnitSystem);Ni=7 #轴瓦数
Nj=7 #单轴瓦轴向x节点数
Nk=60 #周向节点数
R=[30.,32.] #轴瓦内外半径
L=163 #轴瓦间距
x=[0.,2.,6.,12.] #第一个半轴瓦节点轴向坐标
xx=[-x[3],-x[2],-x[1],x[0],x[1],x[2],x[3]] #第一个轴瓦节点轴向坐标
ang=np.linspace(0.,np.pi,Nk,False)
y=np.cos(ang)
z=np.sin(ang)
with open('C:/Users/Admin/Documents/SimLab/nodes.csv', 'w', newline='') as file:writer = csv.writer(file)for i in range(Ni):for j in range(Nj):for k in range(Nk):writer.writerow((i*L+xx[j],y[k],z[k],str(i+1)+str(j+1)+str(k+1).zfill(3))xij=[[[x[0],R[0],0.],[x[1],R[0],0.],[x[2],R[0],0.],[x[3],R[0],0.]],[[x[0],R[1],0.],[x[1],R[1],0.],[x[2],R[1],0.],[x[3],R[1],0.]]]
for i in range(2):for j in range(4):        NodeByXYZ=''' <NodeByXYZ UUID="F200B5A2-D615-4d01-8DE2-25596B3B1EB8"><tag Value="-1"/><Name Value=""/><LocalCoordinateSystem Value="0"/><Position Value="'''+str(xij[i][j][0])+''' mm,'''+str(xij[i][j][1])+''' mm,'''+str(xij[i][j][2])+''' mm"/><Node Value="'''+str(i+1)+str(j+1)+'''"/><UniqueNodeId Value="0"/><ModelName Value="Model_1.gda"/></NodeByXYZ>''';simlab.execute(NodeByXYZ);

在这里插入图片描述

2、由节点建面

2.1、由4个节点建面得到3个面单元Body

for i in range(3):CreateFaceFromFourNodes=''' <FaceUsingFourNodes gda="" UUID="6C69767C-8731-433c-8474-BE95AE45FCDE"><SupportEntities PythonDoc="The four nodes defining the vertices of the face"><Entities><Model>Model_1.gda</Model><Node>1'''+str(i+1)+''',1'''+str(i+2)+''',2'''+str(i+1)+''',2'''+str(i+2)+'''</Node></Entities></SupportEntities><x_no.ofElements Value="1"/><y_no.ofElements Value="1"/><Element_type Value="Quad"/><tag Value="-1"/></FaceUsingFourNodes>''';simlab.execute(CreateFaceFromFourNodes);

在这里插入图片描述

2.2、得到Bodies名称

在这里插入图片描述

2.3、根据Bodies名称选面特征（放入Group_1组）

SelectFeatures=''' <SelectFeatures CheckBox="ON" UUID="CF82E8FB-9B3E-4c02-BA93-9466C1342C6E"><SupportEntities><Entities><Model>Model_1.gda</Model><Body>'''+str(outputBodies).replace("'",'"').strip('()')+'''</Body></Entities></SupportEntities><Arcs Value="0" MinValue="0 mm" MaxValue="0 mm"/><ArcsAll Value="216"/><Circles Value="0" MinValue="0 mm" MaxValue="0 mm"/><CirclesAll Value="0"/><Cones Value="0" MinValue="0 mm" MaxValue="0 mm"/><ConeAll Value="0"/><FullCone Value="0"/><ClosedPartialCone Value="0"/><OpenPartialCone Value="0"/><TaperAngle Angle="0 deg" Value="0"/><Dics Value="0" MinValue="0 mm" MaxValue="0 mm"/><DicsAll Value="0"/><HollowDics Value="0" MinValue="0 mm" MaxValue="0 mm"/><HollowDicsAll Value="0"/><Cylinders Value="0" MinValue="0 mm" MaxValue="0 mm"/><CylindersAll Value="0"/><FullCylinder Value="0"/><ClosedPartialCylinder Value="0"/><OpenPartialCylinder Value="0"/><Fillets Value="0" MinValue="0 mm" MaxValue="0 mm"/><FilletsOption Value="1"/><PlanarFaces Value="1"/><FourEdgedFaces Value="0"/><ConnectedCoaxialFaces Value="0"/><ThroughBoltHole Value="0" MinValue="0 mm" MaxValue="0 mm"/><BlindBoltHole Value="0" MinValue="0 mm" MaxValue="0 mm"/><BlindBoltHoleDepth Value="0" MinValue="0 mm" MaxValue="0 mm"/><SlotEdges Value="0" MinValue="0 mm" MaxValue="0 mm"/><SlotEdgesAll Value="0"/><CreateGrp Value="1" Name="Group_1"/><ArcLengthBased Value=""/><AngleBased Value=""/><SharpEdges Option="" Angle="" Value=""/></SelectFeatures>''';
simlab.execute(SelectFeatures);

2.4、镜像选择的面

Mirroring=''' <Mirroring CheckBox="ON" gda="" UUID="67e16ad8-f140-4976-b9f7-e9c4f9e7005b"><tag Value="-1"/><Name Value=""/><SupportEntities><Entities><Model>Model_1.gda</Model><Body>'''+str(outputBodies).replace("'",'"').strip('()')+'''</Body></Entities></SupportEntities><PlanePoints><RegionObject><Plane><Finite Value="0"/><Pt1 Value="0.000000000000 mm,24.548256049718 mm,6.451743950282 mm,"/><Pt2 Value="0.000000000000 mm,24.548256049718 mm,-6.451743950282 mm,"/><Pt3 Value="0.000000000000 mm,37.451743950282 mm,-6.451743950282 mm,"/><Pt4 Value="0.000000000000 mm,37.451743950282 mm,6.451743950282 mm,"/></Plane></RegionObject></PlanePoints><Copy Value="1"/><WithLBC Value="0"/><LocalCoordinateSystem Value=""/><Output/></Mirroring>''';
simlab.execute(Mirroring);

在这里插入图片描述

3、旋转面得到体SLBodies

Revolve=''' <FEARevolve gda="" UUID="AF802B42-CD64-41a1-B964-52906F7AA888"><tag Value="-1"/><Name Value="FEARevolve1"/><DirNode Value=""/><FirstEntity><Entities><Model>Model_1.gda</Model><Face></Face></Entities></FirstEntity><SecondEntity Value=""/><RowIndex Value=""/><LimitFaces Value=""/><Model Value=""/><layers Value="'''+str(Nk)+'''"/><elementtype Value="0"/><revolvetype Value="0"/><GuideEdge Value=""/><Output/><RegionObject><Axis><Pt1 Value="-4.000000000000 mm,0.000000000000 mm,0.000000000000 mm,"/><Pt2 Value="4.000000000000 mm,0.000000000000 mm,0.000000000000 mm,"/><Angle Value="360.000000000000 deg"/></Axis></RegionObject></FEARevolve>''';
simlab.execute(Revolve);

在这里插入图片描述

4、simlab.getBodiesWithSubString（）得到SLBodies名称

outputBodies=simlab.getBodiesWithSubString("Model_1.gda", ["SLBody*"])

在这里插入图片描述

5、MergeBodies

MergeBodies=''' <BodyMerge gda="" UUID="FA9128EE-5E6C-49af-BADF-4016E5622020"><tag Value="-1"/><Name Value="BodyMerge2"/><SupportEntities><Entities><Model>Model_1.gda</Model><Body>'''+str(outputBodies).replace("'",'"').strip('()')+'''</Body></Entities></SupportEntities><Delete_Shared_Faces Value="0"/><Output_Body_Name Value=""/><RedoFlag Value=""/><Output/></BodyMerge>''';
simlab.execute(MergeBodies);

在这里插入图片描述

6、平移拷贝轴瓦

outputBodies=simlab.getBodiesWithSubString("Model_1.gda", ["SLBody*"])

在这里插入图片描述

Translation=''' <Translation CheckBox="ON" UUID="26b085e8-d9a3-4baf-95f2-9f3eb75ac454"><tag Value="-1"/><Name Value=""/><SupportEntities><Entities><Model>Model_1.gda</Model><Body>'''+str(outputBodies).replace("'",'"').strip('()')+'''</Body></Entities></SupportEntities><RegionObject><TranslationDirection><DirectionVector Value="1.000000000000 mm,0.000000000000 mm,0.000000000000 mm,"/></TranslationDirection></RegionObject><UnitVector Value="0"/><Magnitude Value="1"/><DistanceVector Value="'''+str(L)+''' mm"/><CoorId Value="0"/><Copy Value="1"/><WithLBC Value="0"/><No.Of.Copies Value="'''+str(Ni-1)+'''"/><LocalCoordinateSystem Value=""/><Output/></Translation>''';
simlab.execute(Translation);

在这里插入图片描述

7、根据节点坐标重定义节点号

SetNodeID=''' <SetNodeId UUID="28706164-a6f5-4544-b4a9-c052c4ccc60f"><tag Value="-1"/><Type Value="1" SetNodeIdType="2"/><RenumberPickedNodes><StartNodeId Value="0"/><Offset Value="0"/><CheckDuplicates Value="0"/><Nodes/><SetRBEId Value="0"/></RenumberPickedNodes><RenumberNodesByLocation><Model Value="Model_1.gda" Given="1"/><Tolerance Value="0.001 mm"/><CheckDuplicates Value="0"/><ImportOrExport Option="1" Value="./nodes.csv"/><Nodes Value="" ModelIds="" EntityTypes=""/><NewNodeIds/></RenumberNodesByLocation><RenumberPistonNodes><StartNodeId/><Faces/><StartNodes/><IgnoredNodes Value=""/><ModifyNodes Value=""/><Symmetry Value="0"/><RadialOffset Value="0"/><AxialOffset Value="0"/><Skip Value="0"/><Mirror Value="0"/><CreateFace Value="0"/><PistonTolerance Value="0.0"/><PlaneRegionPoints/></RenumberPistonNodes><RenumberBearingNodes><StartNodeId Value="0"/><Offset Value="0"/><StartNode/><AxialNode/><CircularNode/><Faces/></RenumberBearingNodes><RenumberPinBoss><StartNode/><ArcNode/><StartNodeId/><RadialOffset/><AxialOffset/><Faces/><IncludeMidNodes Value=""/></RenumberPinBoss><RenumberCrankRbeNumbering><JournalRBE Value="0"/><JrnlStartId Value="0"/><IncrAcrossJrnl Value="0"/><IncrWithinJrnl Value="0"/><PinRBE Value="0"/><PinStartId Value="0"/><IncrAcrossPin Value="0"/><IncrWithinPin Value="0"/></RenumberCrankRbeNumbering><RenumberEdgeNodes><Edges/><StartNode/><NextNode/><StartNodeId Value="1"/><Offset Value="1"/></RenumberEdgeNodes><Output/></SetNodeId>''';
simlab.execute(SetNodeID);

### 三级目录