• python二次开发Solidworks：画砂轮

    先根据输入参数计算出绘制砂轮需要的数据，然后绘制草图，完全标注后生成旋转体，具体代码如下：

    1. import sympy as sy
    2. import numpy as np
    3. import matplotlib.pyplot as plt
    4. '''
    5. x1,y1为第一条直线端点坐标(-10,0),theta_l1为角度,取5*np.pi/3
    6. l1为第一条直线的长度，取值2
    7. x2,y2为第二条直线端点坐标(10,0),theta_l2为角度，取4*np.pi/3
    8. x3,y3为第二条直线端点坐标，由其参数方程确定
    9. l2为第二条直线的长度，取值3
    10. r1为第一条圆弧的半径，取值3
    11. r2为第二条圆弧的半径，待求解
    12. l为直线方程的参数
    13. x,y为未知数
    14. d为
    15. theta为圆方程的参数
    16. '''
    17. #输入参数
    18. Line1_p1=[-10,0]
    19. Line1_angle=5*np.pi/3
    20. Line1_length=5
    21. Line2_p1=[10,0]
    22. Line2_angle=4*np.pi/3
    23. Line2_length=3
    24. Arc1_r=8
    25. Arc1_angle=np.pi/6
    26. #定义符号变量
    27. x1,y1,x2,y2,l1,l2,theta_l1,theta_l2=sy.symbols('x1,y1,x2,y2,l1,l2,theta_l1,theta_l2')
    28. r1,theta_arc1_start,theta_arc1_end,xr1_center,yr1_center=sy.symbols('r1,theta_arc1_start,theta_arc1_end,xr1_center,yr1_center')
    29. r2,theta_arc2,theta_arc2_start=sy.symbols('r2,theta_arc2,theta_arc2_start')
    30. l,theta,x,y,d,da=sy.symbols('l,theta,x,y,d,da')
    31. #points存储重要的点坐标
    32. points=[]
    33. #存直线1的起点坐标，第1点
    34. points.append(Line1_p1)
    35. #存直线2的起点坐标，第2点
    36. points.append(Line2_p1)
    37. ​
    38. #第一条直线的参数方程，x1=Line1_p1[0],y1=Line1_p1[1],theta_l1=Line1_angle,l为参数
    39. xl1=x1+l*sy.cos(theta_l1)
    40. yl1=y1+l*sy.sin(theta_l1)
    41. ​
    42. ​
    43. fig=plt.figure(figsize=(5,5))
    44. plt.rcParams['xtick.direction'] = 'in'  # 将x周的刻度线方向设置向内
    45. plt.rcParams['ytick.direction'] = 'in'  # 将y轴的刻度方向设置向内
    46. #生成直线1的横坐标，l1=Line1_length
    47. a=np.arange(0, Line1_length, 0.01)
    48. #计算直线1数据点并绘图
    49. h0L1=xl1.subs([(x1,Line1_p1[0]),(theta_l1,Line1_angle)])
    50. h1L1=yl1.subs([(y1,Line1_p1[1]),(theta_l1,Line1_angle)])
    51. data_l1_x=[]
    52. data_l1_y=[]
    53. for i in a:
    54.     data_l1_x.append(h0L1.subs([(l,i)]))
    55.     data_l1_y.append(h1L1.subs([(l,i)]))
    56. plt.plot(data_l1_x,data_l1_y,c='blue')
    57. ​
    58. #存直线1的终点坐标，也就是圆弧1的起点坐标，第3点
    59. points.append([data_l1_x[-1],data_l1_y[-1]])
    60. ​
    61. #圆弧1初始角度
    62. theta_arc1_start=theta_l1-np.pi/2
    63. start=theta_arc1_start.subs([(theta_l1,Line1_angle)])
    64. #圆弧1终止角度,圆心角取np.pi/3
    65. theta_arc1_end=theta_arc1_start+theta
    66. end=theta_arc1_end.subs([(theta_arc1_start,start),(theta,Arc1_angle)])
    67. #圆弧1的圆心,r1=3,l1=2
    68. xr1_center=x1+l1*sy.cos(theta_l1)+r1*sy.cos(theta_l1+np.pi/2)
    69. yr1_center=y1+l1*sy.sin(theta_l1)+r1*sy.sin(theta_l1+np.pi/2)
    70. ​
    71. #存圆弧1的圆心点坐标，第4点
    72. o1x=xr1_center.subs([(x1,Line1_p1[0]),(theta_l1,Line1_angle),(l1,Line1_length),(r1,Arc1_r)])
    73. o1y=yr1_center.subs([(y1,Line1_p1[1]),(theta_l1,Line1_angle),(l1,Line1_length),(r1,Arc1_r)])
    74. points.append([o1x,o1y])
    75. ​
    76. #第一条圆弧上任意一点的坐标
    77. xr1=xr1_center+r1*sy.cos(theta)
    78. yr1=yr1_center+r1*sy.sin(theta)
    79. ​
    80. ​
    81. #计算圆弧1坐标并作图
    82. b=np.arange(start,end, 0.01)
    83. h0R1=xr1.subs([(x1,Line1_p1[0]),(theta_l1,Line1_angle),(r1,3),(l1,Line1_length)])
    84. h1R1=yr1.subs([(y1,Line1_p1[1]),(theta_l1,Line1_angle),(r1,3),(l1,Line1_length)])
    85. data_r1_x=[]
    86. data_r1_y=[]
    87. for i in b:
    88.     data_r1_x.append(h0R1.subs([(theta, i)]))
    89.     data_r1_y.append(h1R1.subs([(theta, i)]))
    90. plt.plot(data_r1_x,data_r1_y,c='red')
    91. ​
    92. ​
    93. #存圆弧1终止点,，第5点
    94. points.append([data_r1_x[-1],data_r1_y[-1]])
    95. ​
    96. ​
    97. ​
    98. #第二条直线的参数方程，x2=10,y2=0,theta_l2=4*np.pi/3,l为参数
    99. xl2=x2+l*sy.cos(theta_l2)
    100. yl2=y2+l*sy.sin(theta_l2)
    101. #生成直线1的横坐标，l2=3
    102. d=np.arange(0, 3, 0.01)
    103. #计算直线2数据点并绘图
    104. h0L2=xl2.subs([(x2,Line2_p1[0]),(theta_l2,Line2_angle)])
    105. h1L2=yl2.subs([(y2,Line2_p1[1]),(theta_l2,Line2_angle)])
    106. data_l2_x=[]
    107. data_l2_y=[]
    108. for i in d:
    109.     data_l2_x.append(h0L2.subs([(l,i)]))
    110.     data_l2_y.append(h1L2.subs([(l,i)]))
    111. plt.plot(data_l2_x,data_l2_y,c='blue')
    112. ​
    113. ​
    114. #存直线2终点,，第6点
    115. points.append([data_l2_x[-1],data_l2_y[-1]])
    116. ​
    117. #直线2终点处法线方程
    118. x_ln=data_l2_x[-1]+l*sy.cos(theta_l2-np.pi/2)
    119. y_ln=data_l2_y[-1]+l*sy.sin(theta_l2-np.pi/2)
    120. ​
    121. lh0=x_ln.subs([(x2, Line2_p1[0]),(theta_l2,Line2_angle),(l2, Line2_length)])
    122. lh1=y_ln.subs([(y2,Line2_p1[1]),(theta_l2,Line2_angle),(l2, Line2_length)])
    123. ​
    124. #中垂线方程:不包括斜率不存在的情况
    125. x3=data_r1_x[-1]
    126. y3=data_r1_y[-1]
    127. x4=data_l2_x[-1]
    128. y4=data_l2_y[-1]
    129. flz=(x4-x3)*(x-(x4+x3)/2)+(y4-y3)*(y-(y4+y3)/2)
    130. ​
    131. #求解第二条圆弧的圆心及半径
    132. '''
    133. 刚开始把exp=[x-lh0,y-lh1,flz]写成exp=[x-lh0,y-lh0,flz]
    134. 导致圆弧2圆心及半径求解错误，圆弧2始终与其它段接不上
    135. 花了较长时间才找到错误
    136. '''
    137. exp=[x-lh0,y-lh1,flz]
    138. m=sy.solve(exp,[x,y,l])
    139. xr2_center=m[x]
    140. yr2_center=m[y]
    141. r2c=m[l]
    142. #存圆弧2圆心,，第7点
    143. points.append([xr2_center,yr2_center])
    144. ​
    145. #求解第二条圆弧的圆心角
    146. fj=r2c*sy.sin(da)-sy.sqrt((x3-x4)**2+(y3-y4)**2)/2
    147. ​
    148. theta_arc2=sy.solveset(fj,da,domain=sy.Interval(0,np.pi/2)).args[0]
    149. ​
    150. ​
    151. #第二条圆弧的初始角度
    152. theta_arc2_start=theta_l2+np.pi/2
    153. #第二条圆弧的终止角度
    154. theta_arc2_end=theta_arc2_start-2*theta_arc2
    155. ​
    156. ​
    157. ​
    158. #第二条圆弧上任意一点的坐标
    159. xr2=xr2_center+r2c*sy.cos(theta)
    160. yr2=yr2_center+r2c*sy.sin(theta)
    161. ​
    162. start2=theta_arc2_start.subs([(theta_l2,Line2_angle)])
    163. end2=theta_arc2_end.subs([(theta_l2,Line2_angle)])
    164. p=np.arange(end2,start2, 0.01)
    165. data_r2_x=[]
    166. data_r2_y=[]
    167. for i in p:
    168.     data_r2_x.append(xr2.subs([(theta, i)]))
    169.     data_r2_y.append(yr2.subs([(theta, i)]))
    170. plt.plot(data_r2_x,data_r2_y,c='red')
    171. ​
    172. clist = ['blue', 'red', 'green', 'black', 'darkgreen', 'lime', 'gold', 'purple', 'green', 'cyan', 'salmon', 'grey',
    173.              'mediumvioletred', 'darkkhaki', 'gray', 'darkcyan', 'violet', 'powderblue']
    174. markerlst = ['o', '*', 'x', '>', '<', '^', 'D', 'd', '1', '2', '3', '4', 'o', '*', 'x', '>', '<', '^', 'D', 'd',
    175.                  '1', '2', '3', '4', 'o', '*', 'x', '>', '<', '^', 'D', 'd', '1', '2', '3', '4']
    176. plt.xticks(np.linspace(-10,10,11))
    177. plt.yticks(np.linspace(-10, 0.0, 11))
    178. plt.grid(True, linestyle='--', alpha=0.5)
    179. for i in range(len(points)):
    180.     plt.scatter(points[i][0],points[i][1], c=clist[i], marker=markerlst[i],label=f"P{i}")
    181. font1 = {'family': 'SimHei', 'weight': 'normal', 'size': 16}
    182. # 图例展示位置，数字代表第几象限
    183. plt.legend(loc=4,prop=font1,ncol=2,framealpha=0.5)
    184. plt.tight_layout()
    185. plt.show()
    186. ​
    187. ​
    188. import win32com.client as win32
    189. import pythoncom
    190. ​
    191. swApp = win32.Dispatch('sldworks.application')
    192. swApp.Visible = True
    193. Nothing = win32.VARIANT(pythoncom.VT_DISPATCH, None)
    194. Part = swApp.NewDocument(r"C:\ProgramData\SolidWorks\SOLIDWORKS 2016\templates\gb_part.prtdot", 0, 0, 0)
    195. # Part = swApp.ActiveDoc
    196. boolstatus = Part.Extension.SelectByID2("前视基准面", "PLANE", 0, 0, 0, False, 0, Nothing, 0)
    197. Part.SketchManager.InsertSketch(True)
    198. Part.ClearSelection2(True)
    199. Part.SketchManager.CreateLine(points[0][0],points[0][1],0,points[2][0],points[2][1],0)
    200. Part.ClearSelection2(True)
    201. Part.SketchManager.CreateArc(points[3][0],points[3][1],0,points[2][0],points[2][1],0,points[4][0],points[4][1],0,1)
    202. Part.ClearSelection2(True)
    203. Part.SketchManager.CreateArc(points[6][0],points[6][1],0,points[4][0],points[4][1],0,points[5][0],points[5][1],0,1)
    204. Part.ClearSelection2(True)
    205. Part.SketchManager.CreateLine(points[5][0],points[5][1], 0,points[1][0],points[1][1],0)
    206. Part.ClearSelection2(True)
    207. Part.SketchManager.CreateCenterLine(0,5,0,0,-10,0)
    208. Part.ClearSelection2(True)
    209. Part.SketchManager.CreateCenterLine(-12,0,0,12,0,0)
    210. Part.ClearSelection2(True)
    211. boolstatus = Part.Extension.SelectByID2("Point1", "SKETCHPOINT", -10, 0, 0, False, 0, Nothing, 0)
    212. boolstatus = Part.Extension.SelectByID2("Line3", "SKETCHSEGMENT", 0.157744366042774, 0, -1.91530567862208, True, 0, Nothing, 0)
    213. myDisplayDim = Part.AddDimension2(-5,4,0)
    214. Part.ClearSelection2(True)
    215. Part.ClearSelection2(True)
    216. boolstatus = Part.Extension.SelectByID2("Point1", "SKETCHPOINT", -10, 0, 0, False, 0, Nothing, 0)
    217. boolstatus = Part.Extension.SelectByID2("Line4", "SKETCHSEGMENT", -8.77718041256473, -1.03520491789499E-02, 0, True, 0, Nothing, 0)
    218. Part.SketchAddConstraints('sgCOINCIDENT')
    219. Part.ClearSelection2(True)
    220. boolstatus = Part.Extension.SelectByID2("Line1", "SKETCHSEGMENT", -8.6434692451422, -2.28344189536183, 0, False, 0, Nothing, 0)
    221. myDisplayDim = Part.AddDimension2(-8.24233574287464, -0.919587987652101, 0)
    222. Part.ClearSelection2(True)
    223. boolstatus = Part.Extension.SelectByID2("Line4", "SKETCHSEGMENT", -6.58431726683536, 1.63901843055543E-02, 0, False, 0, Nothing, 0)
    224. boolstatus = Part.Extension.SelectByID2("Line1", "SKETCHSEGMENT", -8.88414934650275, -1.88230839309426, 0, True, 0, Nothing, 0)
    225. myDisplayDim = Part.AddDimension2(-12.4943508669108, 1.03259505671672, 0)
    226. Part.ClearSelection2(True)
    227. boolstatus = Part.Extension.SelectByID2("Point8", "SKETCHPOINT", 10, 0, 0, False, 0, Nothing, 0)
    228. boolstatus = Part.Extension.SelectByID2("Line3", "SKETCHSEGMENT", 3.26156827918567E-03, 1.89203245972462, 0, True, 0, Nothing, 0)
    229. myDisplayDim = Part.AddDimension2(3.97469657891235, 3.95129505783071, 0)
    230. Part.ClearSelection2(True)
    231. Part.ClearSelection2(True)
    232. boolstatus = Part.Extension.SelectByID2("Point8", "SKETCHPOINT", 10, 0, 0, False, 0, Nothing, 0)
    233. boolstatus = Part.Extension.SelectByID2("Line4", "SKETCHSEGMENT", 8.09322177512453, -9.36850455919629E-02, 0, True, 0, Nothing, 0)
    234. Part.SketchAddConstraints('sgCOINCIDENT')
    235. Part.ClearSelection2(True)
    236. boolstatus = Part.Extension.SelectByID2("Line2", "SKETCHSEGMENT", 9.08608052778282, -1.45426926219777, 0, False, 0, Nothing, 0)
    237. myDisplayDim = Part.AddDimension2(8.7183550638353, -1.08654379825025, 0)
    238. Part.ClearSelection2(True)
    239. boolstatus = Part.Extension.SelectByID2("Line4", "SKETCHSEGMENT", 7.13713556886099, -9.36850455919629E-02, 0, False, 0, Nothing, 0)
    240. boolstatus = Part.Extension.SelectByID2("Line2", "SKETCHSEGMENT", 9.19639816696708, -1.41749671580302, 0, True, 0, Nothing, 0)
    241. myDisplayDim = Part.AddDimension2(11.9175666001787, 1.15658153182959, 0)
    242. Part.ClearSelection2(True)
    243. boolstatus = Part.Extension.SelectByID2("Point5", "SKETCHPOINT", -6.41126394835776, -5.42278825446442, 0, False, 0, Nothing, 0)
    244. boolstatus = Part.Extension.SelectByID2("Line4", "SKETCHSEGMENT", -6.55131639307265, -0.127252746499565, 0, True, 0, Nothing, 0)
    245. myDisplayDim = Part.AddDimension2(-5.24079660871201, -2.44874493593841, 0)
    246. Part.ClearSelection2(True)
    247. Part.ClearSelection2(True)
    248. boolstatus = Part.Extension.SelectByID2("Point5", "SKETCHPOINT", -6.41126394835776, -5.42278825446442, 0, False, 0, Nothing, 0)
    249. boolstatus = Part.Extension.SelectByID2("Line3", "SKETCHSEGMENT", 0.113612795961455, -5.48166215117303, 0, True, 0, Nothing, 0)
    250. myDisplayDim = Part.AddDimension2(-2.20787939347739, -9.63788203871677, 0)
    251. Part.ClearSelection2(True)
    252. boolstatus = Part.Extension.SelectByID2("Arc2", "SKETCHSEGMENT", 2.50999183022091, -7.35383327168823, 0, False, 0, Nothing, 0)
    253. myDisplayDim = Part.AddDimension2(6.36666433848222, -10.6862978662053, 0)
    254. Part.ClearSelection2(True)
    255. Part.SetPickMode
    256. Part.ClearSelection2(True)
    257. boolstatus = Part.Extension.SelectByID2("Point3", "SKETCHPOINT", -6.60104030516506, -5.58829422721332, 0, False, 0, Nothing, 0)
    258. boolstatus = Part.Extension.SelectByID2("Point5", "SKETCHPOINT", -6.41126394835776, -5.42278825446442, 0, True, 0, Nothing, 0)
    259. Part.SketchAddConstraints("sgMERGEPOINTS")
    260. Part.ClearSelection2(True)
    261. boolstatus = Part.Extension.SelectByID2("Arc1", "SKETCHSEGMENT", -7.0380808844066, -4.77023712537726, 0, False, 0, Nothing, 0)
    262. boolstatus = Part.Extension.SelectByID2("Line1", "SKETCHSEGMENT", -7.8992795998436, -3.45971734101662, 0, True, 0, Nothing, 0)
    263. Part.SketchAddConstraints("sgTANGENT")
    264. Part.ClearSelection2(True)
    265. Part.SketchManager.CreateLine(-10, 0, 0, -10, 40, 0)
    266. Part.SketchManager.CreateLine(-10, 40, 0, 10, 40, 0)
    267. Part.SketchManager.CreateLine(10, 40, 0, 10, 0, 0)
    268. Part.Extension.SelectByID2("Line7", "SKETCHSEGMENT", 0, 0, 0, False, 0, Nothing, 0)
    269. boolstatus = Part.DeSelectByID("Line7", "SKETCHSEGMENT", 0, 0, 0)
    270. boolstatus = Part.Extension.SelectByID2("Line5", "SKETCHSEGMENT", -9.79560003363511, 21.2523700341158, 0, False, 0, Nothing, 0)
    271. myDisplayDim = Part.AddDimension2(-28.8824799708768, 18.5256729002241, 0)
    272. Part.ClearSelection2(True)
    273. Part.ShowNamedView2("*上下二等角轴测", 8)
    274. boolstatus = Part.Extension.SelectByID2("Line6@草图1", "EXTSKETCHSEGMENT", -93.6539010025325, 25.9402592619527, 0, False, 16, Nothing, 0)
    275. myFeature = Part.FeatureManager.FeatureRevolve2(True, True, False, False, False, False, 0, 0, 6.2831853071796, 0, False, False, 0.01, 0.01, 0, 0, 0, True, True, True)
    276. Part.SelectionManager.EnableContourSelection = False
    277. Part.ViewZoomtofit2()

    问题：

    绘制草图时出现了一点小问题，以相同的一个端点绘制的两端圆弧的交点却不重合，原因可能是圆心或半径的计算误差所导致的。

  • 相关阅读:
    2023年中国泌尿内窥镜市场规模、竞争格局及细分产品规模分析[图]
    Axios在vue项目中的基本用法
    ML.NET相关资源整理
    jsp+ajax+json练习（表单数据检测）
    vlan实验
    解决swiper组件autoplay报错问题
    Pandas数据分析17——pandas数据清洗(缺失值、重复值处理)
    OpenCV官方教程中文版 —— 模板匹配
    framework源码读后感
    在Java中使用io实现TCP协议编程
  • 原文地址：https://blog.csdn.net/T20151470/article/details/133913354