| 1.3 主轴参数化有限元模型的建立
基于ANSYS的优化设计必须在结构、热态和模态等简单分析的基础上进行,不同的是需要建立参数化有限元模型。为此,笔者对需要优化的尺寸:2个前轴承之间的距离和前后轴承之间的跨距分别用参数B和L来标定。根据力学原理,把主轴简化为具有空心特征的梁,用BEAM23 单元来模拟,不同轴径的轴段取不同的实常数;把轴承简化为具有一定刚度的压缩弹簧,忽略其刚度的非线性因素和阻尼,用一个COMBIN14 弹簧一阻尼单元和MAss21结构点质量单元来组合模拟,弹簧的一端与主轴上的相应节点相联接,弹簧的另一端完全固接。每个弹簧的刚度按单列轴承的径向刚度并联计算,单个轴承静态预紧后的径向刚度Kr(N/um)按下式进行计算.

式中2为滚动体数目Db为滚动体直径,a为接触角,Fa0为轴向预紧力。
为保证仿真计算的精确性,将电机转子和过盈衬套以分布质量等效到主轴相应的节点上。
根据以上分析,即可建立该轴承一主轴系统的二维参数化有限元模型(见图2)。

1.4 生成优化文件
ANSYS的优化文件可以采用批处理和GUI2种方式生成。本文采用GUI方式生成优化文件。
在ANSYS的图形界面中建立图2所示的有限元模型,其中参数B和L的初值设为原设计数据,在主轴前端节点上施加径向载荷P=l000N,进行结构分析。分析结束后进入后处理,读取主轴前端的横向位移,设为参数Ul;读取主轴电机中部M点处的横向位移,设为参数UM。主轴部件最大静刚度问题也就是最小柔度问题,使主轴的刚度最大,即就是使轴端挠度最小,并且要保证电机中部的横向位移不能超过允许值。因此,结束结构分析进人优化设计模块后,设定B和L为设计变量,UM为状态变量,Ul为目标函数,并为各变量指定取值范围,然后选择优化算法和循环控制方式(本文选择一阶方法),生成优化文件。
1.5 优化计算
运行上一步生成的优化文件,得到表1所示的优化序列表。
由表1可以看出,在第7次迭代得到最优参数:B=20mm,L=412mm时,Min[Ul]= 4.5701um,即当2个前轴承靠紧布置,前后支承跨距缩短10mm时(与原设计值比较),主轴前端的挠度达到最小。此时,主轴的径向静刚度为1000/4.571=219N/mm,此值比原设计刚度(200N/mm)提高了9.5%。

2 一阶自振频率校核
优化迭代结束后ANSYS的当前数据库为最优参数对应的数据库(即第7次迭代结果),故可以从优化模块直接进人模态分析,设定新分析,选择模态提取方法为子空间法,求解并读取主轴第一阶横向自振频率为683 Hz。这一结果与原设计的一阶自振频率(7o6Hz)相比,降低了3%。由于本电主轴设计最高转速为15000r/min,即最高工作频率为250Hz,远比优化后的一阶自振频率低,因此主轴仍能有效避开共振区,在安全转速内工作。
3 结语
通过对主轴支承跨距进行优化,虽然使主轴的一阶自振频率略有降低,但却换来了主轴刚度的较大幅度提高。对于工作在安全转速内的球轴承高速电主轴,刚度是表征其加工能力和加工精度的重要指标,在设计中采用有效措施提高其量值,对电主轴加工性能的提升具有重要意义。 |