大家用AI办公了吗

mayzhang1024

我尝试用WorkBuddy自己写NX后处理，在一问一答间自动生成后处理文件，不知成功不成功，大家感觉怎么样

1. """
   
2. AC双摆头五轴后处理运动学验证脚本
   
3. 用法: python verify_kinematics.py
   
4. """
   
5. import math
   
6. 
   
7. PI = math.pi
   
8. RAD = 180.0 / PI
   
9. 
   
10. 
    
11. def calc_angles(ax, ay, az):
    
12.     """刀轴矢量(I,J,K) -> A角, C角"""
    
13.     # A角: 刀轴与-Z的夹角
    
14.     # az=-1 -> A=0(向下), az=0 -> A=90(水平)
    
15.     a_rad = math.acos(-az)
    
16.     a_ang = a_rad * RAD
    
17. 
    
18.     # C角: atan2(-ax, ay)
    
19.     if abs(a_ang) < 0.5:
    
20.         c_ang = 0.0  # 奇异区冻结
    
21.     else:
    
22.         c_ang = math.atan2(-ax, ay) * RAD
    
23. 
    
24.     return a_ang, c_ang, a_rad
    
25. 
    
26. 
    
27. def calc_machine_pos(Xp, Yp, Zp, L, a_ang, c_ang, a_rad):
    
28.     """编程坐标 -> 机床坐标"""
    
29.     c_rad = c_ang / RAD
    
30.     sinA = math.sin(a_rad)
    
31.     cosA = math.cos(a_rad)
    
32.     sinC = math.sin(c_rad)
    
33.     cosC = math.cos(c_rad)
    
34. 
    
35.     Xm = Xp + L * sinA * sinC
    
36.     Ym = Yp - L * sinA * cosC
    
37.     Zm = Zp + L * cosA
    
38.     return Xm, Ym, Zm
    
39. 
    
40. 
    
41. def verify(name, ax, ay, az, Xp, Yp, Zp, L,
    
42.            exp_X, exp_Y, exp_Z, exp_A, exp_C, tol=0.01):
    
43.     """验证单个测试用例"""
    
44.     a_ang, c_ang, a_rad = calc_angles(ax, ay, az)
    
45.     Xm, Ym, Zm = calc_machine_pos(Xp, Yp, Zp, L, a_ang, c_ang, a_rad)
    
46. 
    
47.     ok = (abs(Xm - exp_X) < tol and abs(Ym - exp_Y) < tol and
    
48.           abs(Zm - exp_Z) < tol and abs(a_ang - exp_A) < tol and
    
49.           abs(c_ang - exp_C) < tol)
    
50. 
    
51.     status = "PASS" if ok else "FAIL"
    
52.     print(f"\n{'='*60}")
    
53.     print(f"  测试: {name}  [{status}]")
    
54.     print(f"{'='*60}")
    
55.     print(f"  刀轴矢量:  I={ax:.4f}  J={ay:.4f}  K={az:.4f}")
    
56.     print(f"  编程坐标:  X={Xp:.3f}  Y={Yp:.3f}  Z={Zp:.3f}")
    
57.     print(f"  摆长 L:    {L:.3f}")
    
58.     print(f"  ---")
    
59.     print(f"  A角:  计算={a_ang:.4f}°  期望={exp_A:.4f}°  {'OK' if abs(a_ang-exp_A)<tol else 'MISMATCH'}")
    
60.     print(f"  C角:  计算={c_ang:.4f}°  期望={exp_C:.4f}°  {'OK' if abs(c_ang-exp_C)<tol else 'MISMATCH'}")
    
61.     print(f"  X轴:  计算={Xm:.4f}    期望={exp_X:.4f}    {'OK' if abs(Xm-exp_X)<tol else 'MISMATCH'}")
    
62.     print(f"  Y轴:  计算={Ym:.4f}    期望={exp_Y:.4f}    {'OK' if abs(Ym-exp_Y)<tol else 'MISMATCH'}")
    
63.     print(f"  Z轴:  计算={Zm:.4f}    期望={exp_Z:.4f}    {'OK' if abs(Zm-exp_Z)<tol else 'MISMATCH'}")
    
64.     return ok
    
65. 
    
66. 
    
67. def main():
    
68.     print("=" * 60)
    
69.     print("  AC双摆头五轴运动学验证 (非RTCP, 摆心距=0)")
    
70.     print("=" * 60)
    
71. 
    
72.     L = 150.0
    
73.     Xp, Yp, Zp = 100.0, 100.0, -50.0
    
74.     all_pass = True
    
75. 
    
76.     # 测试1: A=0°, C=0°, 刀轴垂直向下
    
77.     all_pass &= verify(
    
78.         "A=0 C=0 刀轴向下",
    
79.         ax=0, ay=0, az=-1,
    
80.         Xp=Xp, Yp=Yp, Zp=Zp, L=L,
    
81.         exp_A=0, exp_C=0,
    
82.         exp_X=100, exp_Y=100, exp_Z=100  # Z = -50 + 150*1 = 100
    
83.     )
    
84. 
    
85.     # 测试2: A=90°, C=0°, 刀轴指向+Y
    
86.     all_pass &= verify(
    
87.         "A=90 C=0 刀轴指+Y",
    
88.         ax=0, ay=1, az=0,
    
89.         Xp=Xp, Yp=Yp, Zp=Zp, L=L,
    
90.         exp_A=90, exp_C=0,
    
91.         exp_X=100, exp_Y=-50, exp_Z=-50  # Y=100-150=-50, Z=-50+0=-50
    
92.     )
    
93. 
    
94.     # 测试3: A=90°, C=90°, 刀轴指向-X
    
95.     all_pass &= verify(
    
96.         "A=90 C=90 刀轴指-X",
    
97.         ax=-1, ay=0, az=0,
    
98.         Xp=Xp, Yp=Yp, Zp=Zp, L=L,
    
99.         exp_A=90, exp_C=90,
    
100.         exp_X=250, exp_Y=100, exp_Z=-50  # C=90: sin(90)=1, X=100+150*1*1=250  # X=100+150=175
     
101.     )
     
102. 
     
103.     # 测试4: A=90°, C=180°, 刀轴指向-Y
     
104.     all_pass &= verify(
     
105.         "A=90 C=180 刀轴指-Y",
     
106.         ax=0, ay=-1, az=0,
     
107.         Xp=Xp, Yp=Yp, Zp=Zp, L=L,
     
108.         exp_A=90, exp_C=180,
     
109.         exp_X=100, exp_Y=250, exp_Z=-50  # Y=100+150=250
     
110.     )
     
111. 
     
112.     # 测试5: A=90°, C=-90°, 刀轴指向+X
     
113.     # C=-90: sin(-90)=-1, X=100+150*1*(-1)=100-150=-50
     
114.     all_pass &= verify(
     
115.         "A=90 C=-90 刀轴指+X",
     
116.         ax=1, ay=0, az=0,
     
117.         Xp=Xp, Yp=Yp, Zp=Zp, L=L,
     
118.         exp_A=90, exp_C=-90,
     
119.         exp_X=-50, exp_Y=100, exp_Z=-50
     
120.     )
     
121. 
     
122.     # 测试6: A=45°, C=0°
     
123.     a45 = math.radians(45)
     
124.     all_pass &= verify(
     
125.         "A=45 C=0",
     
126.         ax=0, ay=math.sin(a45), az=-math.cos(a45),
     
127.         Xp=Xp, Yp=Yp, Zp=Zp, L=L,
     
128.         exp_A=45, exp_C=0,
     
129.         exp_X=100, exp_Y=100 - L*math.sin(a45),
     
130.         exp_Z=-50 + L*math.cos(a45)
     
131.     )
     
132. 
     
133.     # 测试7: A=45°, C=45°
     
134.     all_pass &= verify(
     
135.         "A=45 C=45",
     
136.         ax=-math.sin(a45)*math.sin(a45),
     
137.         ay=math.cos(a45)*math.sin(a45),
     
138.         az=-math.cos(a45),
     
139.         Xp=Xp, Yp=Yp, Zp=Zp, L=L,
     
140.         exp_A=45, exp_C=45,
     
141.         exp_X=100 + L*math.sin(a45)*math.sin(a45),
     
142.         exp_Y=100 - L*math.sin(a45)*math.cos(a45),
     
143.         exp_Z=-50 + L*math.cos(a45)
     
144.     )
     
145. 
     
146.     # 测试8: A=110°(过顶), C=0°
     
147.     a110 = math.radians(110)
     
148.     all_pass &= verify(
     
149.         "A=110 C=0 过顶区域",
     
150.         ax=0, ay=math.sin(a110), az=-math.cos(a110),
     
151.         Xp=Xp, Yp=Yp, Zp=Zp, L=L,
     
152.         exp_A=110, exp_C=0,
     
153.         exp_X=100, exp_Y=100 - L*math.sin(a110),
     
154.         exp_Z=-50 + L*math.cos(a110)  # cos(110)<0, Z会减小
     
155.     )
     
156. 
     
157.     # 测试9: 刀尖不动验证 (不同A角, 同一编程点)
     
158.     print("\n" + "=" * 60)
     
159.     print("  测试9: 刀尖定点不动验证")
     
160.     print("=" * 60)
     
161.     print("  编程点固定 (100, 100, -50), L=150")
     
162.     print("  不同刀轴方向, 后处理输出机床坐标应使刀尖不动\n")
     
163. 
     
164.     angles = [0, 15, 30, 45, 60, 75, 90, 105, 110]
     
165.     for a_deg in angles:
     
166.         a_r = math.radians(a_deg)
     
167.         ax_t = 0
     
168.         ay_t = math.sin(a_r)
     
169.         az_t = -math.cos(a_r)
     
170.         a_out, c_out, _ = calc_angles(ax_t, ay_t, az_t)
     
171.         Xm, Ym, Zm = calc_machine_pos(Xp, Yp, Zp, L, a_out, c_out, math.radians(a_out))
     
172.         # 反算刀尖位置: Xm - L*sin(A)*sin(C), Ym + L*sin(A)*cos(C), Zm - L*cos(A)
     
173.         tip_x = Xm - L * math.sin(math.radians(a_out)) * math.sin(math.radians(c_out))
     
174.         tip_y = Ym + L * math.sin(math.radians(a_out)) * math.cos(math.radians(c_out))
     
175.         tip_z = Zm - L * math.cos(math.radians(a_out))
     
176.         ok = abs(tip_x - Xp) < 0.01 and abs(tip_y - Yp) < 0.01 and abs(tip_z - Zp) < 0.01
     
177.         all_pass &= ok
     
178.         print(f"  A={a_deg:6.1f}°  C={c_out:8.3f}°  "
     
179.               f"机床X={Xm:9.3f}  Y={Ym:9.3f}  Z={Zm:9.3f}  "
     
180.               f"刀尖反算=({tip_x:.3f}, {tip_y:.3f}, {tip_z:.3f})  "
     
181.               f"{'OK' if ok else 'FAIL'}")
     
182. 
     
183.     print("\n" + "=" * 60)
     
184.     if all_pass:
     
185.         print("  >>> 所有测试通过! 运动学公式正确 <<<")
     
186.     else:
     
187.         print("  >>> 部分测试失败, 请检查公式 <<<")
     
188.     print("=" * 60)
     
189. 
     
190. 
     
191. if __name__ == "__main__":
     
192.     main()
     

复制代码

白小灰

豆包做文件算嘛

mayzhang1024

白小灰 发表于 2026-8-10 09:41
豆包做文件算嘛

当然算啊，能轻松不少。

北心弦

大佬牛逼9999

Hap2y

大佬nb，666

黄帝啊

没用过嗷嗷嗷

默默Q

  没有过！