浏览量:0

基于FFT的深沟球轴承和角接触球轴承弹性变形快速计算方法

专利类型:发明专利 

语 言:中文 

申 请 号:CN201610203061.7 

申 请 日:20160331 

发 明 人:孟凡明刘越程志涛浦超 

申 请 人:重庆大学中国北方车辆研究所 

申请人地址:400044 重庆市沙坪坝区沙坪坝正街174号 

公 开 日:20180605 

公 开 号:CN105893327B 

代 理 人:赵荣之 

代理机构:北京同恒源知识产权代理有限公司 11275 

摘  要:本发明涉及一种基于FFT的深沟球轴承和角接触球轴承弹性变形快速计算方法,属于滚动轴承优化设计技术领域。该方法基于深沟球轴承和角接触球轴承的点接触应力分布和影响系数矩阵,将离散循环卷积定理和FFT技术应用于轴承弹性变形计算中,计算各滚动体与内圈或外圈之间的弹性变形量及各滚动体与内外圈总的弹性变形量。本发明提供的方法能够有效解决传统方法计算效率不高的问题,因为使用传统方法计算深沟球轴承和角接触球轴承弹性变形时需要进行复杂的四重循环计算,严重耗费计算机内存且计算量很大,而本方法只要正确运用离散循环卷积定理和FFT技术,便能准确、快速地获得计算区域内的轴承弹性变形。 

主 权 项:一种基于FFT的深沟球轴承和角接触球轴承弹性变形快速计算方法,其特征在于:包括以下步骤:S1:给定/修正轴承载荷分布范围参数ε;S2:查询载荷分布积分表并通过线性插值得到载荷分布积分Jr(ε);S3:由式δr=Pd/[2(1?2ε)]计算得到内、外圈在径向载荷方向上的相对位移δr,其中Pd表示径向游隙;S4:由式Frc=ZKnr?Pd/2)1.5Jr(ε)计算得到计算外载荷Frc,其中Z表示滚动体数目,Kn表示滚动体与内、外圈之间总的负荷?变形常数;S5:判断计算外载荷Frc与给定的径向载荷Fr之间的关系是否满足收敛条件,若|(Frc?Fr)/Fr|>ξ,则按照式εnew=εold(Frc/Fr)修正ε并重复步骤S1?步骤S4,其中,ξ取0.001,β在0.001与0.01之间取任意值;若|(Frc?Fr)/Fr|≤ξ,则依据式计算各滚球的载荷分布其中,ψi表示各滚动体的方位角,i=1,…,Z,执行下一步;S6:依据各滚球的载荷分布,运用Hertz接触理论,计算接触椭圆长、短半轴长a、b;S7:计算各滚球在内、外圈接触区上的接触应力分布p,将接触应力分布函数p(x,y)离散为(M+1)×(N+1)的矩阵p,然后补零使其成为2M×2N的应力矩阵p’;S8:以坐标原点为中心,在4倍的计算区域内求得影响系数矩阵C,维数为2M×2N;S9:对矩阵p’和C分别进行FFT变换,得到其FFT变换矩阵S10:的对应项相乘得到频域内的弹性变形矩阵作IFFT变换得到时域上的弹性变形矩阵D;S11:令Vi,j=DM+i,N+j,i=0,1,...,M,j=0,1,...,N,获得待求计算域内的弹性变形Vi,j,分别在各滚球与内、外圈接触区上进行步骤S6?步骤S11,得到各滚球与内、外圈间的弹性变形分布。 

关 键 词:角接触球轴承;深沟球轴承;弹性变形量;快速计算;循环卷积;轴承弹性;滚动体;滚动轴承;变形;影响系数矩阵;优化设计技术;计算机内存;计算区域;计算效率;应力分布;有效解决;点接触;计算量;内外圈;重循环;内圈;应用 

法律状态: 

IPC专利分类号:G06F17/14(2006.01)I