首页 / 专利库 / 轴承类 / 导轨 / 一种滚柱直线导轨副精度可靠性的优化设计方法

一种滚柱直线导轨精度可靠性的优化设计方法

阅读:97发布:2023-01-21

专利汇可以提供一种滚柱直线导轨精度可靠性的优化设计方法专利检索,专利查询,专利分析的服务。并且本 发明 提出了一种滚柱直线 导轨 副 精度 可靠性的优化设计方法,旨在实现对滚柱 直线导轨 副精度可靠性的定量预测,并提高精度可靠性。实现步骤为:建立滚柱直线导轨副精度损失模型δV(S);推导精度可靠性极限状态函数δX(X);建立精度可靠度数学模型RFM、均值灵敏度模型和方差灵敏度模型;选取随机变量xi和随机变量xj;建立多目标稳健优化模型的尺寸约束和可靠度约束;以随机变量xi均值灵敏度最小和随机变量xj方差灵敏度最小为目标,建立滚柱直线导轨副多目标稳健优化模型F(xi,xj);运用最短理想点法对滚柱直线导轨副多目标稳健优化模型进行寻优,得到随机变量xi和随机变量xj的最优解。,下面是一种滚柱直线导轨精度可靠性的优化设计方法专利的具体信息内容。

1.一种滚柱直线导轨精度可靠性的优化设计方法,其特征在于包括如下步骤:
(1)建立滚柱直线导轨副精度损失模型δV(S):
(1.1)计算滚柱直线导轨副的弹性变形量δV1;
(1.2)建立滚柱直线导轨副的滚柱与滚道接触面的磨损预测模型δV2(S);
(1.3)根据滚柱直线导轨副的弹性变形量δV1和滚柱直线导轨副的滚柱与滚道接触面的磨损预测模型δV2(S),建立滚柱直线导轨副精度损失模型δV(S);
(2)推导滚柱直线导轨副精度可靠性极限状态函数δX(X):利用滚柱直线导轨副的精度指标U和滚柱直线导轨副精度损失模型δV的差值,计算滚柱直线导轨副精度可靠性极限状态函数δX(X),δX(X)=U-δV(S),其中,X表示基本随机变量向量,可以从工作参数、材料参数和几何参数中选取;
(3)建立滚柱直线导轨副精度可靠度数学模型RFM:
(3.1)通过Taylor公式,对滚柱直线导轨副精度可靠性极限状态函数δX(X)展开,得到精度可靠性极限状态函数δX(X)保留至二次项的Taylor展开式;
(3.2)将各随机变量的均值分别代入精度可靠性极限状态函数δX(X)的Taylor展开式中,得到极限状态函数δX(X)的均值表达式μg;
(3.3)分别计算Taylor展开式中一阶偏导向量的2次Kronecker幂、3次Kronecker幂和4次Kronecker幂,同时将2次Kronecker幂与基本随机变量向量X的方差向量相乘,3次Kronecker幂和三阶矩向量相乘,4次Kronecker幂与基本随机变量向量X的四阶矩向量相乘,得到极限状态函数δX(X)的方差表达式 三阶矩表达式θg和四阶矩表达式ηg;
(3.4)将极限状态函数δX(X)的均值表达式、方差表达式、三阶矩表达式和四阶矩表达式代入基于HOMST的单失效模式可靠度计算公式,得到可靠性指标表达式βFM;
(3.5)将可靠性指标表达式βFM代入RFM=φ(βFM),得到滚柱直线导轨副精度可靠度数学模型RFM;
(4)将滚柱直线导轨副精度可靠度数学模型RFM分别代入四阶矩均值灵敏度公式和四阶矩方差灵敏度公式,得到基本随机变量向量X的均值灵敏度模型 和基本随机变量向量X的方差灵敏度模型
(5)将滚柱直线导轨副的工作参数、材料参数和各随机变量的均值和方差均代入均值灵敏度模型,计算滚柱直线导轨副每一个跑合位置点的各随机变量的均值灵敏度,得到各随机变量的均值灵敏度曲线,同时,将滚柱直线导轨副的工作参数、材料参数和各随机变量的均值和方差均代入方差灵敏度模型,计算滚柱直线导轨副每一个跑合位置点的各随机变量的方差灵敏度,得到各随机变量的方差灵敏度曲线;
(6)从各随机变量的均值灵敏度曲线中选取偏离X轴最远的曲线所对应的随机变量xi,同时,从各随机变量的方差灵敏度曲线中选取偏离X轴最远的曲线所对应的随机变量xj;
(7)分别选取随机变量xi和随机变量xj的上限和下限,得到滚柱直线导轨副多目标稳健优化模型的尺寸约束:
xi1≤xi≤xi2
xj1≤xj≤xj2
其中,xi1和xi2分别为随机变量xi的下限和上限,xj1和xj2分别为随机变量xj的下限和上限;
(8)将任务给定的目标可靠度R0代入可靠度约束函数FR=μg-φ-1(R0)·σg,得到滚柱直线导轨副多目标稳健优化模型的可靠度约束:
μg-φ-1(R0)·σg≥0;
(9)以随机变量xi的均值对滚柱直线导轨副的可靠度影响和以随机变量xj的方差对滚柱直线导轨副的可靠度影响同时最小为目标,建立滚柱直线导轨副多目标稳健优化模型F(xi,xj);
(10)分别计算随机变量xi和随机变量xj的最优解:
(10.1)在尺寸约束和可靠度约束的范围内,利用非线性规划的寻优算法,对滚柱直线导轨副多目标稳健优化模型F(xi,xj)中的每个单目标函数分别进行寻优,得到两个单目标函数相应的理想点 和
(10.2)判断两个单目标函数相应的理想点 和 是否相同,若是,则该理想
点值就是随机变量xi和随机变量xj的最优解,否则执行步骤(10.3);
(10.3)计算滚柱直线导轨副多目标稳健优化模型F(xi,xj)与理想点集合 的距离模,得到评价函数U(xi,xj):
(10.4)利用非线性规划的寻优算法,对评价函数U(xi,xj)进行寻优,得到该评价函数U(xi,xj)的最优解得到随机变量xi和随机变量xj的最优解;
(11)将滚柱直线导轨副的工作参数、材料参数和各随机变量的均值和方差均代入滚柱直线导轨副精度可靠度数学模型,计算滚柱直线导轨副每一个跑合位置点的精度可靠度,得到精度可靠度曲线;
(12)获取优化后滚柱直线导轨副的可靠度曲线、随机变量xi的均值灵敏度曲线以及方差灵敏度曲线和随机变量xj的均值灵敏度曲线以及方差灵敏度曲线:
(12.1)将滚柱直线导轨副的工作参数、材料参数、随机变量xi的最优解、随机变量xj的最优解和其他各随机变量的均值和方差,代入滚柱直线导轨副的可靠度模型,计算滚柱直线导轨副每一个跑合位置点的精度可靠度,得到优化后滚柱直线导轨副的可靠度曲线;
(12.2)将滚柱直线导轨副的工作参数、材料参数、随机变量xi的最优解、随机变量xj的最优解和其他各随机变量的均值和方差代入均值灵敏度模型,分别计算滚柱直线导轨副每一个跑合位置点的随机变量xi和随机变量xj的均值灵敏度,得到优化后滚柱直线导轨副的随机变量xi的均值灵敏度曲线和随机变量xj的均值灵敏度曲线;
(12.3)将滚柱直线导轨副的工作参数、材料参数、随机变量xi的最优解、随机变量xj的最优解和其他各随机变量的均值和方差代入方差灵敏度模型,分别计算滚柱直线导轨副每一个跑合位置点的随机变量xi和随机变量xj的方差灵敏度,得到优化后滚柱直线导轨副的随机变量xi的方差灵敏度曲线和随机变量xj的方差灵敏度曲线;
(13)将优化后滚柱直线导轨副的可靠度曲线与步骤(11)的精度可靠度曲线、优化后滚柱直线导轨副的随机变量xi的均值灵敏度曲线和步骤(5)的随机变量xi的均值灵敏度曲线、优化后滚柱直线导轨副的随机变量xj的均值灵敏度曲线和步骤(5)的随机变量xj的均值灵敏度曲线、优化后滚柱直线导轨副的随机变量xi的方差灵敏度曲线和步骤(5)随机变量xi的方差灵敏度曲线以及优化后的滚柱直线导轨副的随机变量xj的方差灵敏度曲线和步骤(5)随机变量xj的方差灵敏度曲线分别进行对比,并判断优化后滚柱直线导轨副的各曲线是否均缓于对比曲线,若是,则完成滚柱直线导轨副精度可靠性的优化设计,否则,采用随机变量xi的最优解,分别对随机变量xi上限和下限进行调整,同时采用随机变量xj的最优解分别对随机变量xj的上限和下限进行调整,并执行步骤(7)~步骤(13)。
2.根据权利要求1所述的一种滚柱直线导轨副精度可靠性的优化设计方法,其特征在于,步骤(2)中所述的滚柱直线导轨副精度可靠性极限状态函数,其表达式为:
其中,ν1和ν2分别为滚柱和滚道的泊松比,E1和E2分别为滚柱和滚道的弹性模量、F为外加载荷,Z为承载滚柱数,K为粘着磨损系数,S为跑合里程以及H为材料硬度,Da为滚柱直径,le为滚柱有效长度,α为滚柱与滚道间接触,F0为预紧
3.根据权利要求1所述的一种滚柱直线导轨副精度可靠性的优化设计方法,其特征在于,步骤(3.1)中所述精度可靠性极限状态函数δX(X)保留至二次项的Taylor展开式,具体为:
其中,μX为各随机变量的均值, 为Taylor展开式中一阶偏导向量, 为Taylor展开式中二阶偏导向量。
4.根据权利要求1所述的滚柱直线导轨副精度可靠性的优化设计方法,其特征在于,步骤(3.4)中所述基于HOMST的单失效模式可靠度计算公式,其表达式为:
其中,α3g为滚柱直线导轨副精度可靠性极限状态函数的偏态系数表达式,其表达式为α4g为滚柱直线导轨副精度可靠性极限状态函数的峰值系数表达式,其表达式为 βSM为滚柱直线导轨副精度可靠性极限状态函数的二阶矩可靠性指标。
5.根据权利要求1所述的滚柱直线导轨副精度可靠性的优化设计方法,其特征在于,步骤(4)中所述四阶矩均值灵敏度公式和四阶矩方差灵敏度公式,其表达式分别为:
其中,
6.根据权利要求1所述的滚柱直线导轨副精度可靠性的优化设计方法,其特征在于,步骤(9)中所述建立滚柱直线导轨副多目标稳健优化模型F(xi,xj),其建立步骤为:
(9.1)根据随机变量xi的均值灵敏度模型,建立随机变量xi的均值对滚柱直线导轨副的可靠度影响目标函数f1(xi):
(9.2)根据随机变量xj的方差灵敏度模型,建立随机变量xj的方差对滚柱直线导轨副的可靠度影响目标函数f2(xj):
(9.3)以随机变量xi的均值对滚柱直线导轨副的可靠度影响目标函数f1(xi)和随机变量xj的方差对滚柱直线导轨副的可靠度影响目标函数f2(xj)同时最小为目标,建立滚柱直线导轨副多目标稳健优化模型F(xi,xj):
7.根据权利要求1所述的滚柱直线导轨副精度可靠性的优化设计方法,其特征在于,步骤(10.3)中所述评价函数U(xi,xj),其表达式为:
其中 为滚柱直线导轨副多目标稳健优化模型F(xi,xj)与理想点集

说明书全文

一种滚柱直线导轨精度可靠性的优化设计方法

技术领域

[0001] 本发明属于数控机床技术领域,涉及一种滚柱直线导轨副精度可靠性的优化设计方法,可用于滚柱直线导轨副精度可靠性的定量预测和优化改进。

背景技术

[0002] 根据运动学原理,所谓导轨就是将运动构件约束到只有一个自由度的装置。这一个自由度可以是直线运动或者是回转运动,而做直线运动的导轨被称为直线导轨。滚动直线导轨是一种承受固定、引导移动并减少其摩擦的一种装置,用于直线往复运动场合,拥有比直线轴承更高的额定负载,同时可以承担一定的扭矩,可在高负载的情况下实现高精度的直线运动。按结构特点和摩擦特性,滚动直线导轨又包括滑动导轨、静压导轨和滚动导轨等。其中滚动直线导轨是由导轨、滑滚动体保持架等组成的,根据滚动体的不同,滚动直线导轨副又包括滚柱直线导轨副、滚针直线导轨副和滚珠直线导轨副。滚柱直线导轨副相较于滚珠直线导轨副,其刚性和承载能更高。由于滚柱直线导轨副具有定位精度高、动静摩擦系数小、可维护性好等众多优点,其已经广泛用作机床的关键导向承载领域。
[0003] 近几年来,机床正不断向着高速、高精度和长寿命的方向发展,并对滚柱直线导轨副这一关键功能部件的性能和可靠性提出了更高的要求。滚柱直线导轨副可靠性是指导轨副在规定的条件和时间内完成规定功能的能力。规定的条件有导轨副在设备上的安装形式、跑和速度、加载情况和工作环境等,规定的时间即是导轨副在使用过程中的跑合里程或者跑合总时间;规定的功能是指在指滚动直线导轨副在跑合的过程中性能良好,精度、噪声和振动在可接受范围内,而且没有出现卡死、剥落和点蚀等情况出现。在滚柱直线导轨副的使用过程中,噪音和振动对滚柱直线导轨副带来的直接影响也可以由滚柱直线导轨副的精度保持性体现。精度是滚柱直线导轨副性能的重要性能参数,精度的高低直接影响机床等设备的加工或运行精度。滚柱直线导轨副的精度保持性是指滚柱直线导轨副在工作过程中保持原有精度指标的能力。精度保持性低的滚柱直线导轨副在使用一定时间后,导轨、滑块及滚动体之间发生接触变形和磨损,影响机械传动精度,进而影响设备运行或加工精度。精度可靠性是指精度保持稳定的概率值。所以对滚柱直线导轨副的精度可靠性的研究与优化设计进行深入研究是十分必要的。
[0004] 从目前的公开的资料来看,相关领域的现有的滚柱直线导轨副的精度可靠性优化设计方法方面的研究较为匮乏,部分学者通过ANSYS仿真获得试验数据,并应用相应可靠性理论对试验数据进行分析的方法对滑动导轨的磨损可靠性进行研究,但是由于滚柱直线导轨副与滑动导轨副的结构及磨损机理不大相同,且该方法对接触面表面其他因素考虑不够全面,不能完整描述整个精度损失过程,只对精度损失初始阶段的跑合磨损适用,无法满足缺少试验数据情况下的可靠性研究。而在对于滚柱直线导轨副的优化设计中,目前主要是通过对滚柱直线导轨副进行仿真实验,通过分析实验结果对滚柱直线导轨副进行结构优化,从而提高滚柱直线导轨副的刚度。例如,宋现春等人在其发表的论文“高速滚柱直线导轨副结构优化设计及其性能测试”(《制造技术与机床》2015,(12):141-145)中公开了一种高速滚柱直线导轨副结构优化设计方法,该方法利用多体动力学理论结合数值计算分析,通过对滚柱直线导轨副的运动仿真分析,对滚柱直线导轨副中滚柱、滑块、整体结构布局及其主要结构参数进行优化设计,改善其受力状况,提高其结构刚度。但是该种方法并未充分考虑实际工况中相关参数的分散性和不确定性,且滚柱直线导轨副作为关键功能零部件,其精度保持性维持时间的长短对于提高数控机床的可靠性是至关重要的。

发明内容

[0005] 本发明的目的在于克服上述现有技术存在的缺陷,提供了一种滚柱直线导轨副精度可靠性的优化设计方法,旨在实现对滚柱直线导轨副精度可靠性的定量预测,并提高滚柱直线导轨副精度可靠性。
[0006] 为实现上述目的,本发明采取的技术方案包括如下步骤:
[0007] (1)建立滚柱直线导轨副精度损失模型δV(S):
[0008] (1.1)计算滚柱直线导轨副的弹性变形量δV1;
[0009] (1.2)建立滚柱直线导轨副的滚柱与滚道接触面的磨损预测模型δV2(S);
[0010] (1.3)根据滚柱直线导轨副的弹性变形量δV1和滚柱直线导轨副的滚柱与滚道接触面的磨损预测模型δV2(S),建立滚柱直线导轨副精度损失模型δV(S);
[0011] (2)推导滚柱直线导轨副精度可靠性极限状态函数δX(X):利用滚柱直线导轨副的精度指标U和滚柱直线导轨副精度损失模型δV的差值,计算滚柱直线导轨副精度可靠性极限状态函数δX(X),δX(X)=U-δV(S),其中,X表示基本随机变量向量,可以从工作参数、材料参数和几何参数中选取;
[0012] (3)建立滚柱直线导轨副精度可靠度数学模型RFM:
[0013] (3.1)通过Taylor公式,对滚柱直线导轨副精度可靠性极限状态函数δX(X)展开,得到精度可靠性极限状态函数δX(X)保留至二次项的Taylor展开式;
[0014] (3.2)将各随机变量的均值分别代入精度可靠性极限状态函数δX(X)的Taylor展开式中,得到极限状态函数δX(X)的均值表达式μg;
[0015] (3.3)分别计算Taylor展开式中一阶偏导向量的2次Kronecker幂、3次Kronecker幂和4次Kronecker幂,同时将2次Kronecker幂与基本随机变量向量X的方差向量相乘,3次Kronecker幂和三阶矩向量相乘,4次Kronecker幂与基本随机变量向量X的四阶矩向量相乘,得到极限状态函数δX(X)的方差表达式 三阶矩表达式θg和四阶矩表达式ηg;
[0016] (3.4)将极限状态函数δX(X)的均值表达式、方差表达式、三阶矩表达式和四阶矩表达式代入基于HOMST的单失效模式可靠度计算公式,得到可靠性指标表达式βFM;
[0017] (3.5)将可靠性指标表达式βFM代入RFM=φ(βFM),得到滚柱直线导轨副精度可靠度数学模型RFM;
[0018] (4)将滚柱直线导轨副精度可靠度数学模型RFM分别代入四阶矩均值灵敏度公式和四阶矩方差灵敏度公式,得到基本随机变量向量X的均值灵敏度模型 和基本随机变量向量X的方差灵敏度模型
[0019] (5)将滚柱直线导轨副的工作参数、材料参数和各随机变量的均值和方差均代入均值灵敏度模型,计算滚柱直线导轨副每一个跑合位置点的各随机变量的均值灵敏度,得到各随机变量的均值灵敏度曲线,同时,将滚柱直线导轨副的工作参数、材料参数和各随机变量的均值和方差均代入方差灵敏度模型,计算滚柱直线导轨副每一个跑合位置点的各随机变量的方差灵敏度,得到各随机变量的方差灵敏度曲线;
[0020] (6)从各随机变量的均值灵敏度曲线中选取偏离X轴最远的曲线所对应的随机变量xi,同时,从各随机变量的方差灵敏度曲线中选取偏离X轴最远的曲线所对应的随机变量xj;
[0021] (7)分别选取随机变量xi和随机变量xj的上限和下限,得到滚柱直线导轨副多目标稳健优化模型的尺寸约束:
[0022] xi1≤xi≤xi2
[0023] xj1≤xj≤xj2
[0024] 其中,xi1和xi2分别为随机变量xi的下限和上限,xj1和xj2分别为随机变量xj的下限和上限;
[0025] (8)将任务给定的目标可靠度R0代入可靠度约束函数FR=μg-φ-1(R0)·σg,得到滚柱直线导轨副多目标稳健优化模型的可靠度约束:
[0026] μg-φ-1(R0)·σg≥0;
[0027] (9)以随机变量xi的均值对滚柱直线导轨副的可靠度影响和以随机变量xj的方差对滚柱直线导轨副的可靠度影响同时最小为目标,建立滚柱直线导轨副多目标稳健优化模型F(xi,xj);
[0028] (10)分别计算随机变量xi和随机变量xj的最优解:
[0029] (10.1)在尺寸约束和可靠度约束的范围内,利用非线性规划的寻优算法,对滚柱直线导轨副多目标稳健优化模型F(xi,xj)中的每个单目标函数分别进行寻优,得到两个单目标函数相应的理想点 和
[0030] (10.2)判断两个单目标函数相应的理想点 和 是否相同,若是,则该理想点值就是随机变量xi和随机变量xj的最优解,否则执行步骤(10.3);
[0031] (10.3)计算滚柱直线导轨副多目标稳健优化模型F(xi,xj)与理想点集合的距离模,得到评价函数U(xi,xj):
[0032] (10.4)利用非线性规划的寻优算法,对评价函数U(xi,xj)进行寻优,得到该评价函数U(xi,xj)的最优解得到随机变量xi和随机变量xj的最优解;
[0033] (11)将滚柱直线导轨副的工作参数、材料参数和各随机变量的均值和方差均代入滚柱直线导轨副精度可靠度数学模型,计算滚柱直线导轨副每一个跑合位置点的精度可靠度,得到精度可靠度曲线;
[0034] (12)获取优化后滚柱直线导轨副的可靠度曲线、随机变量xi的均值灵敏度曲线以及方差灵敏度曲线和随机变量xj的均值灵敏度曲线以及方差灵敏度曲线:
[0035] (12.1)将滚柱直线导轨副的工作参数、材料参数、随机变量xi的最优解、随机变量xj的最优解和其他各随机变量的均值和方差,代入滚柱直线导轨副的可靠度模型,计算滚柱直线导轨副每一个跑合位置点的精度可靠度,得到优化后滚柱直线导轨副的可靠度曲线;
[0036] (12.2)将滚柱直线导轨副的工作参数、材料参数、随机变量xi的最优解、随机变量xj的最优解和其他各随机变量的均值和方差代入均值灵敏度模型,分别计算滚柱直线导轨副每一个跑合位置点的随机变量xi和随机变量xj的均值灵敏度,得到优化后滚柱直线导轨副的随机变量xi的均值灵敏度曲线和随机变量xj的均值灵敏度曲线;
[0037] (12.3)将滚柱直线导轨副的工作参数、材料参数、随机变量xi的最优解、随机变量xj的最优解和其他各随机变量的均值和方差代入方差灵敏度模型,分别计算滚柱直线导轨副每一个跑合位置点的随机变量xi和随机变量xj的方差灵敏度,得到优化后滚柱直线导轨副的随机变量xi的方差灵敏度曲线和随机变量xj的方差灵敏度曲线;
[0038] (13)将优化后滚柱直线导轨副的可靠度曲线与步骤(11)的精度可靠度曲线、优化后滚柱直线导轨副的随机变量xi的均值灵敏度曲线和步骤(5)的随机变量xi的均值灵敏度曲线、优化后滚柱直线导轨副的随机变量xj的均值灵敏度曲线和步骤(5)的随机变量xj的均值灵敏度曲线、优化后滚柱直线导轨副的随机变量xi的方差灵敏度曲线和步骤(5)随机变量xi的方差灵敏度曲线以及优化后的滚柱直线导轨副的随机变量xj的方差灵敏度曲线和步骤(5)随机变量xj的方差灵敏度曲线分别进行对比,并判断优化后滚柱直线导轨副的各曲线是否均缓于对比曲线,若是,则完成滚柱直线导轨副精度可靠性的优化设计,否则,采用随机变量xi的最优解,分别对随机变量xi上限和下限进行调整,同时采用随机变量xj的最优解分别对随机变量xj的上限和下限进行调整,并执行步骤(7)~步骤(13)。
[0039] 本发明与现有技术相比,具有如下优点:
[0040] 1、本发明由于采用滚柱直线导轨副精度损失模型δV(S),从弹性变形和磨损两方面充分考虑滚柱直线导轨副与轨道之间的精度损失过程,并运用涉及的假设较少,结果误差较小的四阶矩法对滚柱直线导轨副精度损失进行可靠性分析,建立了滚柱直线导轨副精度可靠度数学模型RFM,实现了滚柱直线导轨副精度可靠性的定量预测。
[0041] 2、本发明充分考虑滚柱直线导轨副设计中客观存在的不确定性因素,将可靠性设计理论、稳健优化设计理论和滚柱直线导轨副精度损失模型相结合,建立了一种多目标稳健优化模型,因而能在获得满足设计要求的最优解的同时,保证该最优解的可靠性和稳健性。附图说明
[0042] 图1为本发明的实现流程图
[0043] 图2为本发明实施例的滚柱直线导轨副受力分析图;
[0044] 图3为本发明实施例的滚柱直线导轨副各随机变量均值灵敏度随跑合里程增加的变化曲线图;
[0045] 图4为本发明实施例的滚柱直线导轨副各随机变量方差灵敏度随跑合里程增加的变化曲线图;
[0046] 图5为本发明实施例的滚柱直线导轨副精度可靠度RFM随跑合里程增加的变化曲线图;
[0047] 图6为本发明实施例的滚柱直线导轨副精度可靠度RFM与优化后滚柱直线导轨副精度可靠度 的对比曲线图;
[0048] 图7为本发明实施例的滚柱直线导轨副随机变量xα的均值灵敏度与优化后滚柱直线导轨副的随机变量xα的均值灵敏度对比曲线图;
[0049] 图8为本发明实施例的滚柱直线导轨副随机变量xα的方差灵敏度与优化后滚柱直线导轨副的随机变量xα的方差灵敏度对比曲线图。

具体实施方式

[0050] 以下结合附图和具体实施例,对本发明作进一步详细说明。
[0051] 参照图1、一种滚柱直线导轨副精度可靠性的优化设计方法,包括如下步骤:
[0052] 步骤1、建立滚柱直线导轨副精度损失模型δV(S):建立导轨副精度损失模型时要从弹性变形和磨损两方面考虑,建立精度损失模型δV(S)的步骤如下:
[0053] 步骤1.1)计算滚柱直线导轨副的弹性变形量δV1:
[0054] 设定滚柱直线导轨副受到垂直向下加载在滑块上的作用载荷F,每个滚柱承受的接触面法向载荷为Qn,接触面法线方向与垂直方向的夹为α,Z为导轨上每排滚柱的数目,滚柱直线导轨副所受预紧力F0。其受力情况如附图2所示。根据受力平衡可得:
[0055]
[0056] 采用Palmgren经验公式,综合考虑法向接触载荷Qn和法向接触变形δn联系以及预紧力F0和初始变形为δ0联系,计算得到滚柱直线导轨副的弹性变形量δV1为:
[0057]
[0058] 其中,ν1和ν2分别为滚柱和滚道的泊松比,E1和E2分别为滚柱和滚道的弹性模量
[0059] 步骤1.2)建立滚柱直线导轨副的滚柱与滚道接触面的磨损预测模型δV2(S):根据Archard粘着磨损理论,可得到滚柱的体积磨损量WV计算公式:
[0060]
[0061] 其中,K为粘着磨损系数,S为运行里程,H为较软材料的硬度,弹性蠕滑率[0062] 根据Hertz弹性接触理论,可得到单列滚柱的轮廓接触面积Ap为:
[0063]
[0064] 其中载荷W=Qn+F0/2Z=F/(2Zcosα)+F0/2Z,R为滚柱半径;
[0065] 所以滚柱直线导轨副的滚柱与滚道接触面的磨损预测模型δV2(S)为:
[0066]
[0067] 步骤1.3)根据滚柱直线导轨副的弹性变形量δV1和滚柱直线导轨副的滚柱与滚道接触面的磨损预测模型δV2(S),建立滚柱直线导轨副精度损失模型δV(S):
[0068]
[0069] 其中,δV(S)为滑块对导轨底面基准位移偏差的总变动量;
[0070] 步骤2、推导滚柱直线导轨副精度可靠性极限状态函数δX(X):利用滚柱直线导轨副的精度指标U和滚柱直线导轨副精度损失模型δV的差值,计算滚柱直线导轨副精度可靠性极限状态函数δX(X),δX(X)=U-δV(S),其中,X表示基本随机变量向量,选取外加载荷F、材料硬度H、滚柱直径Da、滚柱有效直径le、接触角α以及预紧力F0。由于各随机变量均服从正态分布,故查阅资料可知各随机变量的偏态系数为:CsX=(0,0,0,0,0,0)T,各随机变量的峰度系数为:CkX=(3,3,3,3,3,3)T;
[0071] 步骤3、建立滚柱直线导轨副精度可靠度数学模型RFM:
[0072] 步骤3.1)通过Taylor公式,对滚柱直线导轨副精度可靠性极限状态函数δX(X)展开,得到精度可靠性极限状态函数δX(X)保留至二次项的Taylor展开式:
[0073]
[0074] 步骤3.2)将各随机变量的均值分别代入精度可靠性极限状态函数δX(X)的Taylor展开式中,且各随机变量相对独立,得到极限状态函数δX(X)的均值表达式μg:
[0075]
[0076] 步骤3.3)分别计算Taylor展开式中一阶偏导向量的2次Kronecker幂、3次Kronecker幂和4次Kronecker幂,同时将2次Kronecker幂与基本随机变量向量X的方差向量相乘,3次Kronecker幂和三阶矩向量相乘,4次Kronecker幂与基本随机变量向量X的四阶矩向量相乘,得到极限状态函数δX(X)的方差表达式 三阶矩表达式θg和四阶矩表达式ηg分别为:
[0077]
[0078]
[0079]
[0080] 步骤3.4)将极限状态函数δX(X)的均值表达式、方差表达式、三阶矩表达式和四阶矩表达式代入基于HOMST的单失效模式可靠度计算公式,得到可靠性指标表达式βFM,其表达式为:
[0081]
[0082] 其中,α3g为滚柱直线导轨副精度可靠性极限状态函数的偏态系数表达式,其表达式为 α4g为滚柱直线导轨副精度可靠性极限状态函数的峰值系数表达式,其表达式为 βSM为滚柱直线导轨副精度可靠性极限状态函数的二阶矩可靠性指标;
[0083] 步骤3.5)将可靠性指标表达式βFM代入RFM=φ(βFM),得到滚柱直线导轨副精度可靠度数学模型RFM;
[0084] 步骤4、将滚柱直线导轨副精度可靠度数学模型RFM分别代入四阶矩均值灵敏度公式和四阶矩方差灵敏度公式,四阶矩均值灵敏度公式和四阶矩方差灵敏度公式,其表达式分别为:
[0085]
[0086]
[0087] 其中,
[0088]
[0089]
[0090]
[0091]
[0092]
[0093] 得到基本随机变量向量X的均值灵敏度模型 和基本随机变量向量X的方差灵敏度模型
[0094] 步骤5、采用计算机应用软件MATLAB分别对各基本随机变量向量X的均值灵敏度模型 和各基本随机变量向量X的方差灵敏度模型 进行编程,并以表1所示的滚柱直线导轨副的主要参数进行计算,可分别得到本发明实施例的滚柱直线导轨副各随机变量均值灵敏度随跑合里程增加的变化曲线图和本发明实施例的滚柱直线导轨副各随机变量方差灵敏度随跑合里程增加的变化曲线图,分别见附图3和附图4:
[0095] 表1滚柱直线导轨副导轨主要参数表
[0096]
[0097] 步骤6、由于外加载荷F、预紧力F0、和材料硬度H具有不确定性,故应从随机变量中的设计参数接触角α、滚柱直径Da和有效长度le的中选取需要优化的随机变量。从本发明实施例的滚柱直线导轨副各随机变量均值灵敏度随跑合里程增加的变化曲线图中可以看出,接触角α的均值灵敏度随跑合里程增加的向负方向增大,即接触角α对滚柱直线导轨副的精度可靠性呈负影响,且在运动过程影响越来越大。并且从本发明实施例的滚柱直线导轨副各随机变量方差灵敏度随跑合里程增加的变化曲线图中可以看出,接触角α的方差灵敏度随跑合里程增加的向负方向先增大后减小,但接触角α对滚柱直线导轨副的精度可靠性呈负影响,且在运动过程较其他随机变量的影响是最大的,综上所述,选取接触角α作为需要优化的随机变量;
[0098] 步骤7、选取随机变量xα的上限和下限,得到滚柱直线导轨副多目标稳健优化模型的尺寸约束:
[0099] 40≤xα≤50
[0100] 步骤8、将任务给定的目标可靠度R0=0.9和跑合里程S=618km代入可靠度约束函数FR=μg-φ-1(R0)·σg,得到滚柱直线导轨副多目标稳健优化模型的可靠度约束:
[0101] μg-φ-1(0.9)·σg≥0;
[0102] 步骤9、以随机变量xα的均值对滚柱直线导轨副的可靠度影响和以随机变量xα的方差对滚柱直线导轨副的可靠度影响同时最小为目标,建立滚柱直线导轨副多目标稳健优化模型F(xα):
[0103]
[0104] 步骤10、计算随机变量xα最优解:
[0105] 步骤10.1)在尺寸约束和可靠度约束的范围内,分别利用MATLAB优化工具箱中的fmincon函数求解,设定接触角α的初值α0为45,对滚柱直线导轨副多目标稳健优化模型F(xα)中的每个单目标函数分别进行寻优,得到两个单目标函数相应的理想点分别为和
[0106] 步骤10.2)由于两个单目标函数相应的理想点 和 不同,则计算滚柱直线导轨副多目标稳健优化模型F(xα)与理想点集合 的距离模,得到评价函数U(xα):
[0107]
[0108] 步骤10.3)利用MATLAB优化工具箱中的fmincon函数求解,对评价函数U(xα)进行寻优,得到该评价函数U(xα)的最优解得到随机变量xα最优解为40.0034;
[0109] 步骤11、采用计算机应用软件MATLAB对滚柱直线导轨副精度可靠度数学模型RFM进行编程,通过表1滚柱直线导轨副主要参数的均值和方差,得到滚柱直线导轨副运动精度可靠度RFM随跑合里程增加的变化曲线,如图5所示,本实例的滚柱直线导轨副在运动初期精度可靠性较高,当跑和里程达到400km时,精度可靠性开始下降,当跑和里程达到800km时,精度可靠度达到0.5;
[0110] 步骤12、获取优化后滚柱直线导轨副的可靠度曲线、随机变量xα的均值灵敏度曲线以及方差灵敏度曲线:
[0111] 步骤12.1)将滚柱直线导轨副的工作参数、材料参数、随机变量xα的最优解其他各随机变量的均值和方差,代入滚柱直线导轨副的可靠度模型,采用计算机应用软件MATLAB,计算滚柱直线导轨副每一个跑合位置点的精度可靠度,得到优化后滚柱直线导轨副的精度可靠度 变化曲线;
[0112] 步骤12.2)将滚柱直线导轨副的工作参数、材料参数、随机变量xα的最优解和其他各随机变量的均值和方差代入均值灵敏度模型,采用计算机应用软件MATLAB,分别计算滚柱直线导轨副每一个跑合位置点的随机变量xα,得到优化后滚柱直线导轨副的随机变量xα的均值灵敏度曲线;
[0113] 步骤12.3)将滚柱直线导轨副的工作参数、材料参数、随机变量xα的最优解和其他各随机变量的均值和方差代入方差灵敏度模型,采用计算机应用软件MATLAB,分别计算滚柱直线导轨副每一个跑合位置点的随机变量xα,得到优化后滚柱直线导轨副的随机变量xα的方差灵敏度曲线;
[0114] 步骤13、根据附图(6)所示,将优化后滚柱直线导轨副的可靠度曲线与步骤(11)的精度可靠度曲线进行对比,发现优化后滚柱直线导轨副的可靠度曲线下降明显晚于步骤(11)的精度可靠度曲线,趋势较为缓和,且相同跑合里程时的精度可靠性明显提升;根据附图(7)所示,将滚柱直线导轨副的随机变量xα的均值灵敏度曲线和步骤(5)的随机变量xα的均值灵敏度曲线进行对比,发现优化后滚柱直线导轨副的随机变量xα的均值灵敏度曲线下降明显晚于步骤(5)的随机变量xα的均值灵敏度曲线,且趋势较为缓和,接触角α的负影响减小,使得滚柱直线导轨副精度可靠性不易于失效;根据附图(8)所示,将优化后滚柱直线导轨副的随机变量xα的方差灵敏度曲线和步骤(5)随机变量xα的方差灵敏度曲线进行对比,优化后滚柱直线导轨副的随机变量xα的方差灵敏度曲线开始下降明显晚于步骤(5)的随机变量xα的方差灵敏度曲线,且趋势较为缓和,接触角α的负影响减小,使得滚柱直线导轨副精度可靠性不易于失效。随机变量xα的最优解可以有效提高精度可靠性,降低接触角α的均值灵敏度和方差灵敏度,完成滚柱直线导轨副的优化设计。
相关专利内容
标题 发布/更新时间 阅读量
导轨电车 2020-05-13 68
一种导轨 2020-05-13 573
一种电梯方向导引对重导轨架 2020-05-12 782
一种直线导轨端面专用倒角装置 2020-05-12 726
一种抗压电梯导轨 2020-05-11 541
升降机导轨的制造方法及调整设备 2020-05-12 658
一种直线导轨端面专用倒角装置 2020-05-13 373
导轨 2020-05-11 961
导轨 2020-05-11 98
直线导轨 2020-05-12 77
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈