滚动轴承动力学模型附上程序和网上的paper 程序百分百为博主自研并且花费了较大精力，故可以保质保量，可以对照程序和文章学习建模，以便考虑新的因素，故对轴承动力学小白十分友好 后支持程序 刚性保持架模型：综合考虑滚动体与保持架的相互作用关系，滚动体与内外圈的接触力和摩擦力，阻尼作用，滚动体离心力，得到了内圈质心轨迹，保持架转速，保持架打滑率，滚动体与保持架接触力，滚动体打滑率，滚动体公转、自转、径向加速度等动力学响应 柔性保持架模型：在刚性模型基础上根据论文内容进行了模型建立 可以额外输出保持架相邻质量块间的弹簧作用力等（此论文未提现） 单独刚性保持架模型的提前沟通好，联系别付款我改价 需要整个paper文件夹的即可

程序概述

本程序包提供了一套完整的滚动轴承动力学仿真工具，包含刚性保持架和柔性保持架两种模型。程序基于MATLAB环境开发，采用四阶龙格-库塔数值积分方法求解轴承系统的动力学响应。

核心算法模块

1. 数值求解器 (RK_4.m)

• 实现四阶龙格-库塔定步长数值积分算法
• 输入参数：微分方程函数句柄、时间范围、初始条件、步长
• 输出：时间序列和状态变量矩阵

2. 摩擦系数计算 (u.m)

提供两种摩擦系数模型：

% 原始文献模型（注释状态）
% function u=u(deltaV)
%  if abs(deltaV) < 0.05
%      u=0.4*(abs(deltaV));
%  else
%      u=0.002;
%  end
% end

% 当前采用的指数衰减模型
function u=u(deltaV)
  A=-7.5e-3; B=1.99e-2; C=1.6; D=7.5e-3;
  if deltaV==0
     u=0;
  else
     u=((A+B.*abs(deltaV))*exp(-C*abs(deltaV))+D)*(deltaV/abs(deltaV));
  end
end

刚性保持架模型

模型文件

• 动力学函数：Rigidcage6304.m
• 主程序：mainRigidcage_6304.m

物理参数设定

基于深沟球轴承6304规格：

• 内圈质量：0.4kg
• 滚动体数目：7个
• 滚动体半径：4.76mm
• 节圆半径：18mm
• 径向游隙：5μm
• 接触刚度：内圈3.19e9 N/m^1.5，外圈3.38e9 N/m^1.5
• 系统阻尼：300 N·s/m

建模特点

1. 保持架假设：视为刚性整体，具有转动惯量Jc=3e-5 kg·m²
2. 接触力学：
   - 滚动体与内外圈：Hertz接触理论，指数n=1.5
   - 滚动体与保持架：线性刚度Krc=1e8 N/m，摩擦系数uc=0.002
3. 自由度配置：
   - 内圈：径向平移(x,y) + 转动
   - 滚动体：径向位移 + 公转 + 自转
   - 保持架：整体转动

输出响应

程序生成15个图形输出，包括：

• 保持架转速时域图
• 滚动体与保持架碰撞力
• 内圈接触力（时域和角度域）
• 内圈X/Y向加速度
• 内圈质心轨迹
• 保持架打滑率
• 滚动体打滑率
• 滚动体运动状态（公转、自转、径向加速度）
• 内圈位移时程

柔性保持架模型

模型文件

• 动力学函数：flexcage6304.m
• 主程序：mainflexcage_6304.m

扩展功能

在刚性模型基础上增加：

1. 保持架离散化：将保持架分解为多个质量块
2. 单元间耦合：引入弹簧刚度Kcc=5.2083e7 N/m连接相邻质量块
3. 附加自由度：每个保持架离散单元具有独立的转动自由度

新增输出

• 保持架相邻质量块间弹簧作用力
• 单个保持架离散单元的运动响应
• 保持架系统内部相互作用力

后处理与可视化程序 (huizhi.m)

功能特性

1. 数据加载：从matlab.mat文件读取仿真结果
2. 多视图展示：
   - 时域接触力曲线
   - 角度域接触力分布
   - 单周期接触力细节
   - 极坐标载荷分布图

1. 极坐标绘图特色：

% 承载区与非承载区分段绘制
theta1=linspace(0,a,L1); % 承载区角度范围
theta2=linspace(a,2*pi-a,L2); % 非承载区
theta3=linspace(2*pi-a,2*pi,L3); % 承载区

% 极坐标样式设置
rlim([0 800]) % 半径范围
thetaticks(0:45:360) % 角度刻度
ThetaZeroLocation = 'bottom' % 0度位置

程序执行流程

1. 参数初始化：设定轴承几何、材料、工况参数
2. 数值积分：调用RK_4求解动力学方程
3. 后处理计算：提取接触力、速度、加速度等响应量
4. 结果可视化：生成多维度图形展示动力学行为

技术特点

1. 物理完整性：涵盖接触力学、摩擦、离心效应、阻尼耗散
2. 数值稳定性：采用鲁棒的龙格-库塔算法
3. 输出丰富性：提供全面的动力学响应指标
4. 模型对比性：刚性/柔性保持架模型便于性能对比分析

本程序包为滚动轴承动态特性研究提供了从基础到高级的完整仿真框架，适用于轴承设计优化、故障诊断和性能评估等工程应用场景。

滚动轴承动力学模型附上程序和网上的paper 程序百分百为博主自研并且花费了较大精力，故可以保质保量，可以对照程序和文章学习建模，以便考虑新的因素，故对轴承动力学小白十分友好 后支持程序 刚性保持架模型：综合考虑滚动体与保持架的相互作用关系，滚动体与内外圈的接触力和摩擦力，阻尼作用，滚动体离心力，得到了内圈质心轨迹，保持架转速，保持架打滑率，滚动体与保持架接触力，滚动体打滑率，滚动体公转、自转、径向加速度等动力学响应 柔性保持架模型：在刚性模型基础上根据论文内容进行了模型建立 可以额外输出保持架相邻质量块间的弹簧作用力等（此论文未提现） 单独刚性保持架模型的提前沟通好，联系别付款我改价 需要整个paper文件夹的即可