一种基于有限差分法的重型龙门铣床横梁重力变形预测方法
【技术领域】
[0001] 本发明涉及重型龙门铣床横梁重力变形预测方法，特别涉及一种基于有限差分法 的重型龙门铣床横梁重力变形预测方法。
【背景技术】
[0002] 重型数控机床作为加工母机广泛应用于国防、航空航天、能源、船舶、冶金等重点 领域[1]，其精度的好坏基于有限差分法的重型龙门铣床横梁重力变形曲线计算方法映了 以一个国家制造业的水平。由于重型龙门铣床横梁自身大尺寸、大跨距等结构因素，在自身 重力作用下，会造成一定程度的变形，且该重力引起的变形误差无法忽略。
[0003] 横梁作为重型龙门铣床的核心部件，垂直刀架移动对工作台面的平行度（G5项精 度）是其最重要的精度指标。通过对横梁承载面加工反变形曲线进行补偿，能有效提高机 床的G5项精度。但由于铸造过程的不可控性，重型机床的结构件无可避免地存在夹砂、气 孔等各种缺陷，导致横梁材料属性、尺寸等不一致，使目前横梁反变形计算采用的有限元方 法计算准确性仅能达到40%~50%，横梁需经过多次实验校核，反复拆装修配才能满足精 度要求，成本较高且非常耗时。因此，准确计算横梁重力变形曲线，能减少横梁拆装次数，降 低成本。图2为重型龙门铣床横梁俯视图。
[0004] 张雁亭（《双柱立式车床横梁弹性变形的补偿方法》）通过近似计算得到双柱立式 车床横梁的弹性变形曲线，提出获得合理的导轨几何形状所需采用的预变形方法，提高了 机床的精度。该方法计算过程过于繁琐，且由于采用近似简化模型，计算精度较差。郭铁能 (《重载大跨距横梁承载曲线分析与实验研宄》）等利用ANSYS对重型龙门铣床进行有限元 分析，得到横梁上25个等间距工作位置的变形量，绘制得到横梁的承载曲线，通过实验表 明在预测横梁承载曲线时需要增加7%~16%的预估量。该方法仅通过对比有限元分析与 实验结果给出加工承载面时的预估量缺乏理论支撑，泛用性较差。王道明（《龙门机床横梁 的仿形加工》）等提出一种基于有限元分析，同时结合实际检测的方法得出横梁反变形加工 曲线，降低了成本，提高了装配效率。该方法并未全面考虑材料属性的不均一性，仅横梁自 重变形曲线由实验得到，而外力作用曲线由有限元仿真得到。
[0005] 综上所述，理论计算方法过程过于繁琐，计算精度较差，但能通过公式中的横梁材 料属性反映出横梁实际变形情况；使用有限元分析方便快速，计算精度较高，但前处理过程 仅能定义部件整体的材料属性，无法考虑实际材料的不均一性，不符合实际情况，导致计算 结果与实际变形值相差很大。
【发明内容】
[0006] 本发明的目的是为了解决现有的有限元分析计算方法无法在实际材料属性不均 一的情况下准确计算横梁重力变形曲线，导致计算结果与实际变形值相差大的问题，而提 出一种基于有限差分法的重型双柱立车横梁重力变形预测方法。
[0007] 上述的发明目的是通过以下技术方案实现的：
[0008] 步骤一、通过在横梁与立柱装配处放置垫铁，模拟实际装配条件设计重型机床横 梁自重变形实验，得到横梁自重变形曲线；
[0009] 步骤二、利用材料力学理论，根据横梁在自重作用下的受力情况将横梁简化为横 梁自重变形模型和横梁扭转变形模型；
[0010] 步骤三、对步骤二得到的横梁自重变形模型离散化后结合有限差分法建立横梁重 力变形离散化模型：
[0011] 步骤四、结合步骤一所述重型机床横梁自重变形曲线和步骤三所述横梁重力变形 离散化模型，计算横梁离散微段的当量抗弯刚度；其中，当量抗弯刚度表征横梁的材质属 性；
[0012] 步骤五、通过有限元方法模拟重型机床横梁的实际装配条件，将横梁与垂直刀架 装配后计算横梁有限元重力变形曲线；其中，横梁有限元重力变形包括横梁弯曲变形和扭 转变形；
[0013] 步骤六、利用横梁重力变形有限元仿真分离方法将有限元重力变形曲线分离得到 垂直刀架装配后横梁弯曲变形有限元曲线和扭转变形有限元曲线；
[0014] 步骤七、利用步骤四计算得到的当量抗弯刚度，基于有限差分法对步骤六计算得 到垂直刀架装配后横梁弯曲变形及横梁扭转变形的有限元曲线进行校正，得到最终的横梁 重力变形曲线；即完成了一种基于有限差分法的重型龙门铣床横梁重力变形预测方法。
[0015] 发明效果
[0016] 由于理论计算方法与有限元计算方法具有很好的互补性，因此为解决由于横梁材 质、制造工艺等因素导致有限元分析结果不准确的问题，本发明基于有限差分法提出了一 种结合材料力学、自重变形实验与有限元方法的横梁重力变形曲线计算方法。
[0017] 由于重力导致的弯曲变形和扭转变形对龙门铣床的精度均具有较大的影响，需要 同时考虑其对龙门铣床精度的影响。实现基于有限差分法的龙门铣床横梁重力变形曲线计 算方法的主要流程如下：通过模拟实际装配条件设计重型机床横梁自重变形实验，得到考 虑材质不均一性的横梁自重变形曲线。利用材料力学理论将横梁简化为简支梁弯曲变形力 学模型和固支梁扭转变形力学模型，将该模型离散成微段，结合有限差分法建立重型机床 横梁的重力变形离散化模型。结合横梁自重变形实验、重力变形离散化模型计算横梁各离 散段的当量抗弯刚度，表征横梁的实际材料属性。通过有限元方法模拟横梁的实际装配条 件，计算垂直刀架装配后横梁的有限元重力变形曲线，并利用横梁弯曲变形和扭转变形的 有限元仿真数据分离方法得到弯曲变形及扭转变形的仿真数据。应用计算得到的当量抗弯 刚度，基于有限差分法针对弯曲变形和扭转变形校正有限元方法计算的重力变形曲线，得 到准确的横梁重力变形曲线。将横梁重力变形仿真曲线、实际变形曲线及基于有限差分法 的重力变形曲线计算方法校正结果进行对比，如图15所示。验证该计算方法的正确性。
[0018] 从图16可以看出校正后的曲线相比有限元仿真曲线更贴近实际的横梁变形情 况。经计算，原有限元计算结果与实际横梁变形的平均误差率为26. 86%，而基于有限差分 法的横梁重力变形计算方法得到的横梁Z向变形与实际横梁变形的平均误差率为8. 37%， 主要加工区域的误差值最大为0.0564mm。证明了基于有限差分法和扭转计算的有限元结果 校正方法的正确性。
【附图说明】
[0019] 图1为【具体实施方式】一提出的横梁重力反变形曲线计算流程图；
[0020] 图2为【具体实施方式】四提出的重型龙门铣床横梁俯视图；
[0021] 图3为【具体实施方式】二提出的建立的重型龙门铣床横梁俯视图坐标系示意图；
[0022] 图4为实施例提出的横梁约束条件示意图；其中，A、B和C表示进行有限元分析时 添加约束类型的编号；
[0023] 图5为实施例提出的横梁载荷定义不意图；
[0024] 图6为【具体实施方式】四提出的横梁重力载荷弯曲计算模型简图；其中，L为简支梁 支点之间长度的一半，2L为简支梁支点之间的长度；Ll为两端矩形梁的长度；L2为中段矩 形梁的长度的一半，2L2为中段矩形梁的长度；a为外伸梁的长度；ql为两端矩形梁截面的 重力载荷集度值；qll为中段矩形梁截面的重力载荷集度值;A和B为两个支点位置的名称 符号，易于表述；
[0025] 图7为【具体实施方式】四提出的横梁刀架重力载荷扭转计算模型简图；L为固支梁 支点之间长度的一半，2L为固支梁支点之间的长度，与简支梁力学模型中的L相同；S表示 的是刀架移动到该位置处的名称符号，易于表述，如刀架位置S ;s为刀架在X轴方向上距离 坐标系原点O的距离；
[0026] 图8为【具体实施方式】五提出的横梁离散化模型示意图；
[0027] 图9为实施例提出的重力作用下的Z向变形示意图；
[0028] 图10为实施例提出的横梁刀尖点Z向重力变形仿真曲线示意图；
[0029] 图11为【具体实施方式】七提出的横梁肋板截面扭心有限元计算示意图；
[0030] 图12为【具体实施方式】七提出的横梁无肋板截面扭心有限元计算示意图；
[0031] 图13为【具体实施方式】七提出的横梁扭心观察点设置示意图；
[0032] 图14为实施例