一种机床主轴的流固热耦合数值仿真方法.pdf

《一种机床主轴的流固热耦合数值仿真方法.pdf》由会员分享，可在线阅读，更多相关《一种机床主轴的流固热耦合数值仿真方法.pdf（9页完整版）》请在专利查询网上搜索。

1、10申请公布号CN102063549A43申请公布日20110518CN102063549ACN102063549A21申请号201110002309022申请日20110107G06F17/5020060171申请人西安交通大学地址710049陕西省西安市咸宁西路28号申请人西安瑞特快速制造工程研究有限公司72发明人李宝童洪军张进华杜飞郭俊康74专利代理机构西安通大专利代理有限责任公司61200代理人陆万寿54发明名称一种机床主轴的流固热耦合数值仿真方法57摘要本发明公开了一种机床主轴的流固热耦合数值仿真方法，建立主轴冷却液的模型，将主轴部件中的轴承用相同尺寸的等效圆环代替，进行模型简化；计。

2、算主轴温度场相关外部边界条件；导入网格划分软件ANSYSICEMCFD中，进行流体和固体网格划分；导入到CFX软件中进行流固热耦合计算；将简化后的几何模型导入WORKBENCH中稳态温度场模块，将冷却液模型忽略，并重新划分网格；根据温度场边界条件，对主轴重新进行温度场计算；各个轴承的径向和轴向刚度计算等效圆环的径向和轴向的弹性模量，并在WORKBENCH中修改各个等效圆环的弹性模量；将稳态温度场计算得到的温度导入WORKBENCH稳态静力分析模块，计算热变形。51INTCL19中华人民共和国国家知识产权局12发明专利申请权利要求书2页说明书5页附图1页CN102063559A1/2页21一种机。

3、床主轴的流固热耦合数值仿真方法，其特征在于，包括如下步骤1在几何建模软件PRO/E中建立主轴冷却液的模型，将主轴部件中的轴承用相同尺寸的等效圆环代替，并对主轴几何模型进行模型简化；2根据下面公式，计算主轴温度场相关外部边界条件；HF1047104NM1其中HF为轴承发热量，N为轴承转速，M为轴承摩擦力矩；P1P2其中P1为传递到主轴锥面的切削热量，P为主轴切削功率，为切削热传递到主轴锥面处的效率；其中为滑枕内部空气强迫对流换热系数，NU为努谢尔特数，为冷却液热传导系数，DE为特征尺寸；15W/M2K4其中1为滑枕外部与空气接触处的自然对流换热系数；3将简化后的模型导入网格划分软件ANSYSIC。

4、EMCFD中，进行流体和固体网格划分；4根据步骤2计算所得边界条件，将步骤3划分好的网格模型导入到CFX软件中进行流固热耦合计算；5将步骤1简化后的几何模型导入WORKBENCH中稳态温度场模块，将冷却液模型忽略，并重新划分网格；6将步骤4CFX中计算得到的冷却液流道的温度导入WORKBENCH稳态温度场模块，根据步骤2计算所得温度场边界条件，对主轴重新进行温度场计算；7根据各个轴承的径向和轴向刚度计算等效圆环的径向和轴向的弹性模量，并在WORKBENCH中修改各个等效圆环的弹性模量；8将步骤6稳态温度场计算得到的温度导入WORKBENCH稳态静力分析模块，计算热变形。2根据权利要求1所述仿真。

5、方法，其特征在于，所述的几何模型中包括冷却液模型，利用与轴承尺寸相同的圆环代替主轴部件中的轴承。3根据权利要求1所述仿真方法，其特征在于，所述的主轴部件几何模型的简化方法为1将小于1MM的台阶删除；将螺纹孔全部删除；2将主轴部件内部的拉刀机构删除，在热变形计算的有限元模型中利用相同质量的质点代替；3将主轴部件中的电机删除，并在热变形计算的有限元模型中利用相同质量的质点代替；4保留主轴支架零件，保证主轴支架与主轴连接处的接触面积不变。4根据权利要求1所述的仿真方法，其特征在于，所述的在ANSYSICEMCFD中的网格划分方法为划分网格时必须对流体定义边界层，边界层厚度控制在1MM2MM之间。权利。

6、要求书CN102063549ACN102063559A2/2页35根据权利要求1所述的仿真方法，其特征在于，所述的CFX软件中的流固热耦合计算流体的边界条件包括，流体入口速度、入口温度、出口压强，在CFX计算中将主轴支座做绝热处理。6根据权利要求1所述的仿真方法，其特征在于，所述的网格重新划分方法为选择MULTIZONE法进行网格划分。7根据权利要求1所述的仿真方法，其特征在于，所述的等效圆环弹性模量确定方法为1径向弹性模量的确定在WORKBENCH结构分析模块中，将圆环外圈固定约束，然后在圆环内圈施加一个径向的面力；当圆环内圈径向变形的平均值与相应轴承刚度对应的变形值相同时，此时圆环弹性模量。

7、即为所求径向弹性模量；2轴向弹性模量的确定在WORKBENCH结构分析模块中，在圆环一侧施加轴向压力，在另一侧约束轴向自由度，在圆环外圈约束住除轴向以外的自由度；轴向压力施加位置为对应轴承内圈的侧面；轴向约束处的为之为轴承外圈的侧面；当内圈轴向变形平均值与轴承刚度值对应的变形值相同时，此时圆环弹性模量即为所求轴向弹性模量。8根据权利要求1所述的仿真方法，其特征在于，所述的热变形计算中主轴轴承与主轴轴套之间接触选为BONDED，其余轴承与轴套的接触均选为NOSEPARATION；轴承隔套与轴或者轴套的接触均选为NOSEPARATION，其与旁边轴承的接触均选为BONDED。权利要求书CN1020。

8、63549ACN102063559A1/5页4一种机床主轴的流固热耦合数值仿真方法技术领域0001本发明涉及一种计算机辅助分析工具领域的方法，具体是一种机床主轴的流固热耦合数值仿真方法。背景技术0002机床主轴发热变形不但会严重影响加工精度，而且会影响轴承寿命。在主轴设计阶段利用有限元软件对其进行温度场及热变形分析，对比不同情况下的发热及热变形，并优化主轴冷却参数和结构对于主轴的设计定型具有重要意义。0003目前尚没有利用有限元软件对主轴热变形进行准确计算的方法。目前的主轴有限元计算方法主要有00041二维计算方法，即认为主轴模型为轴对称的，将其简化为二维模型进行温度场及热变形计算。对于第一种。

9、方法，实际主轴并不是轴对称的，所以会存在较大的误差。00052单个零件的三维计算方法，目前主轴的三维计算方法主要是对于单个零件如单个主轴零件、单个主轴箱零件进行三维刚度、模态计算，或将几个零件粘为一个零件进行计算。该方法则无法预知整个主轴组建的温度场与热变形情况，一般不用来计算温度场及热变形。00063考虑多个零件的三维计算方法，主要计算主轴部件的刚度、模态、温度场、热变形等。该方法首先将零件模型进行装配，并对装配体进行整体求解，其结果的精度可大为提高。所以这种是计算温度场及热变形的主要方法。0007在主轴温度场计算中，主轴冷却液换热系数对于温度场分析精确与否至关重要，但目前上述第一种及第三种。

10、方法均为利用简化的理论公式进行估算，但在工程实践中往往很难全部满足公式所要求的各种假设条件，所以较常用的方法是先根据公式求得基本值，然后利用实验进行修正，以得到经验参数。该方法过于复杂，效率低。0008在热变形计算中，轴承的简化方法对结果的影响至关重要，目前轴承一般采用弹簧单元代替，或者采用圆环代替，前者一般不能用于热变形；后者一般在计算中采用各项同性材料，而此时则不能反映主轴所受轴向力对于主轴变形的影响，计算精度差。0009另外在一般的主轴热变形计算中，一般直接对主轴进行位移约束，没有考虑主轴支架滑枕座的影响，而实际上滑枕座对主轴变形影响巨大。发明内容0010针对上述现有技术存在的缺陷或不足。

11、，本发明的目的在于提供一种主轴部件温度场及热变形的精确计算方法。0011本发明的技术方案是这样解决的0012包括如下步骤00131在几何建模软件PRO/E中建立主轴冷却液的模型，将主轴部件中的轴承用相同尺寸的等效圆环代替，并对主轴几何模型进行模型简化；说明书CN102063549ACN102063559A2/5页500142根据下面公式，计算主轴温度场相关外部边界条件；0015HF10471O4NM10016其中HF为轴承发热量，N为轴承转速，M为轴承摩擦力矩；0017P1P20018其中P1为传递到主轴锥面的切削热量，P为主轴切削功率，为切削热传递到主轴锥面处的效率；00190020其中为滑。

12、枕内部空气强迫对流换热系数，NU为努谢尔特数，为冷却液热传导系数，DE为特征尺寸；002115W/M2K40022其中1为滑枕外部与空气接触处的自然对流换热系数；00233将简化后的模型导入网格划分软件ANSYSICEMCFD中，进行流体和固体网格划分；00244根据步骤2计算所得边界条件，将步骤3划分好的网格模型导入到CFX软件中进行流固热耦合计算；00255将步骤1简化后的几何模型导入WORKBENCH中稳态温度场模块，将冷却液模型忽略，并重新划分网格；00266将步骤4CFX中计算得到的冷却液流道的温度导入WORKBENCH稳态温度场模块，根据步骤2计算所得温度场边界条件，对主轴重新进行。

13、温度场计算；00277根据各个轴承的径向和轴向刚度计算等效圆环的径向和轴向的弹性模量，并在WORKBENCH中修改各个等效圆环的弹性模量；00288将步骤6稳态温度场计算得到的温度导入WORKBENCH稳态静力分析模块，计算热变形。0029所述的几何模型中包括冷却液模型，利用与轴承尺寸相同的圆环代替主轴部件中的轴承。0030所述的主轴部件几何模型的简化方法为00311将小于1MM的台阶删除；将螺纹孔全部删除；00322将主轴部件内部的拉刀机构删除，在热变形计算的有限元模型中利用相同质量的质点代替；00333将主轴部件中的电机删除，并在热变形计算的有限元模型中利用相同质量的质点代替；00344保。

14、留主轴支架零件，保证主轴支架与主轴连接处的接触面积不变。0035所述的在ANSYSICEMCFD中的网格划分方法为划分网格时必须对流体定义边界层，边界层厚度控制在1MM2MM之间。0036所述的CFX软件中的流固热耦合计算流体的边界条件包括，流体入口速度、入口温度、出口压强，在CFX计算中将主轴支座做绝热处理。0037所述的网格重新划分方法为选择MULTIZONE法进行网格划分。0038所述的等效圆环弹性模量确定方法为说明书CN102063549ACN102063559A3/5页600391径向弹性模量的确定在WORKBENCH结构分析模块中，将圆环外圈固定约束，然后在圆环内圈施加一个径向的面。

15、力；当圆环内圈径向变形的平均值与相应轴承刚度对应的变形值相同时，此时圆环弹性模量即为所求径向弹性模量；00402轴向弹性模量的确定在WORKBENCH结构分析模块中，在圆环一侧施加轴向压力，在另一侧约束轴向自由度，在圆环外圈约束住除轴向以外的自由度；轴向压力施加位置为对应轴承内圈的侧面；轴向约束处的为之为轴承外圈的侧面；当内圈轴向变形平均值与轴承刚度值对应的变形值相同时，此时圆环弹性模量即为所求轴向弹性模量。0041所述的热变形计算中主轴轴承与主轴轴套之间接触选为BONDED，其余轴承与轴套的接触均选为NOSEPARATION；轴承隔套与轴或者轴套的接触均选为NOSEPARATION，其与旁边。

16、轴承的接触均选为BONDED。0042本发明的技术效果00431、解决了一般机床主轴稳态温度场仿真中不能精确考虑冷却液影响的问题，主轴温度场的数值仿真更加精确，避免了传统方法利用实验来修正冷却液壁面换热系数带来的效率低下问题，使在设计阶段对温度场进行准确的估计成为可能。00442、在热变形计算中，更加准确方便的模拟了轴承的刚度特性，使仿真结果更加准确。00453、在热变形计算中，更加准确的考虑了主轴部件内部各零件之间的接触情况。00464、在热变形计算中，更加准确的考虑了主轴支座滑枕座对于主轴变形的影响。附图说明0047图1是圆环轴向等效弹性模量计算边界条件位置0048其中，A施加轴向压力位置。

17、；0049B施加轴向约束位置；0050C施加径向约束位置；0051下面结合附图对本发明的内容作进一步详细说明。具体实施方式00521、几何建模005311利用几何建模软件PRO/E对原来滑枕模型做必要的简化00541删除小于1MM的阶梯；00552删除螺纹孔；00563删除滑枕后端进给电机及同步带部分；00574删除拉刀机构；00585保留后尾筒上面薄钢板；00596保留简化后的滑枕座模型，滑枕座模型仅保留其筋板及与滑枕接触处，其余细节特征全部简化掉。006012利用PRO/E软件在原来滑枕模型的基础上，根据冷却液流道的尺寸建立冷却液的几何模型。006113利用PRO/E软件，利用与轴承相同尺。

18、寸的圆环代替滑枕中的轴承。说明书CN102063549ACN102063559A4/5页700622、计算温度场分析边界条件0063为进行共轭传热计算及稳态温度场计算，考虑的温度场边界条件有轴承发热、切削热、滑枕内部旋转部件与空气的强制对流换热、滑枕外部与空气的自然对流换热、滑枕零件与主轴箱、滑枕座接触处的温度，其值按照下面选取0064HF1047104NM10065其中HF为轴承发热量，N为轴承转速，M为轴承摩擦力矩；0066P1P20067其中P1为传递到主轴锥面的切削热量，P为主轴切削功率，为切削热传递到主轴锥面处的效率；00680069其中为滑枕内部空气强迫对流换热系数，NU为努谢尔特。

19、数，为冷却液热传导系数，DE为特征尺寸；007015W/M2K40071其中1为滑枕外部与空气接触处的自然对流换热系数。0072滑枕与主轴箱、滑枕座等其他部件接触处的温度取35。00733、计算等效圆环弹性模量0074等效圆环采用各项异性材料。根据滑枕所选用轴承刚度，确定与其相对应的等效圆环的弹性模量。其径向及轴向弹性模量分别根据轴承径向及轴向刚度进行计算。且对于角接触轴承组，每个轴承的刚度按照轴承组的平均值计算，具体按下面方法进行计算00751径向弹性模量的确定在WORKBENCH结构分析模块中，将圆环外圈固定约束，然后在圆环内圈施加一个径向的面力；当圆环内圈径向变形的平均值与相应轴承刚度对。

20、应的变形值相同时，此时圆环弹性模量即为所求径向弹性模量。00762轴向弹性模量的确定在WRKBENCH结构分析模块中，在圆环一侧施加轴向压力，在另一侧约束轴向自由度，在圆环外圈约束住除轴向以外的自由度；轴向压力施加位置为对应轴承内圈的侧面；轴向约束处的为之为轴承外圈的侧面；当内圈轴向变形平均值与轴承刚度值对应的变形值相同时，此时圆环弹性模量即为所求轴向弹性模量。00774、ANSYSICEMCFD网格划分0078为保证流体边界处在径向有足够密的网格，在划分网格时须对流体定义边界层，此时定义为2MM，网格划分完成后，将其保存为GTM格式。00795、CFX共轭传热计算0080在CFX软件中，导入。

21、GTM格式的网格文件。计算前定义流体计算域、固体计算域。流体计算域的边界条件为流体入口速度及温度、出口压力。固体计算域的边界条件为上述步骤2所述的温度、热流密度、对流换热边界条件；在轴承的固体计算域中设置计算子域，定义发热体边界条件。在各个计算域之间设置计算域交界面并定义相关接触面热阻。其中滑枕座计算域按照绝热处理。0081由于冷却液为油基冷却液，根据入口流速及冷却液粘度，流体模型采用层流模型进行计算，计算后得到冷却液流道壁面的温度及换热系数。00826、网格重新划分说明书CN102063549ACN102063559A5/5页80083将PRO/E中几何模型导入WORKBENCH中的稳态温度。

22、场模块，将冷却液模型忽略，并利用划分网格。为同时保证温度场计算和热变形计算的计算精度，必须对于温度梯度大和主要受力部位的网格大小进行控制，网格划分方法采用MULTIZONE。00847、稳态温度场的计算0085导入CFX计算所得冷却液流道壁面温度及简化的滑枕座外表面温度，其余边界条件与上述第二步中相同，重新计算温度场。00868、WORKBENCH中进行热变形计算0087跟据第7步计算结果设置轴承等效圆环的轴向及径向弹性模量。0088将上述第7步计算所得整体温度场导入结构分析中。0089进行接触设置轴承与轴之间均为BONDED；主轴轴承与主轴轴套之间设置为BONDED；其余轴承与轴套的接触均为。

23、NOSEPARATION；轴承隔套与轴承间的接触均为BONDED；轴承隔套与轴/轴套的接触均为NOSEPARATION；运动结合面包括镗轴与钢套之间的接触面等有相互运动的结合面接触设置均为NOSEPARATION；镗轴中锁紧螺母与轴均为BONDED。滑枕座与滑枕间的接触为增广拉格朗日接触。0090设置滑枕力边界条件和位移边界条件进行计算。0091本发明包含有如下步骤在几何建模软件PRO/E中建立主轴冷却液的模型，将主轴部件中的轴承用相同尺寸的等效圆环代替，并对主轴几何模型进行模型简化；根据公式14，计算主轴温度场相关外部边界条件；将简化后的模型导入网格划分软件ANSYSICEMCFD中，进行流。

24、体和固体网格划分；根据上述计算所得温度场边界条件，将划分好的网格模型导入到CFX软件中进行流固热耦合计算；将简化后的几何模型导入WORKBENCH中稳态温度场模块，将冷却液模型忽略，并重新划分网格；将CFX中计算得到的冷却液流道的温度导入WORKBENCH稳态温度场模块，根据温度场边界条件，对主轴重新进行温度场计算；根据各个轴承的径向和轴向刚度计算等效圆环的径向和轴向的弹性模量，并在WORKBENCH中修改各个等效圆环的弹性模量；将稳态温度场计算得到的温度导入WORKBENCH稳态静力分析模块，计算热变形。说明书CN102063549ACN102063559A1/1页9图1说明书附图CN102063549A。