车轨耦合动力学外部激励——车轮扁疤

最新推荐文章于 2024-06-06 11:05:42 发布

cofisher

最新推荐文章于 2024-06-06 11:05:42 发布

阅读量1.1k

点赞数

分类专栏： Matlab 车轨耦合动力学文章标签：动力学车轨耦合

本文链接：https://blog.csdn.net/qq_36758914/article/details/124865483

版权

Matlab 同时被 2 个专栏收录

8 篇文章 24 订阅

订阅专栏

车轨耦合动力学

4 篇文章 5 订阅

订阅专栏

本文详细介绍了如何在Matlab中使用轮径变化模拟车轮扁疤效应，通过《车辆－轨道耦合动力学》提供的公式，展示了如何在动力学模型中集成这种影响。关键步骤包括在循环中计算车轮半径变化和轮对响应的调整。参考资料丰富，覆盖扁疤对车辆振动的影响及计算方法。

摘要由CSDN通过智能技术生成

问题说明

扁疤会对车辆系统的振动产生影响，一般采用轮径变化来模拟扁疤激励。

车轮扁疤模型

本文采用车轮轮径变化来描述扁疤模型，如下图所示，在车轮滚动一周内，车轮圆周的半径变化为
$\Delta r(x)\approx h-\frac{1}{2R}(x-L/2)^2$
其中， $0\leq x\leq L$ ， $h\approx \frac{L^2}{8R}$ ， $\Delta r(x)$ 为车轮圆周半径变化值， $h$ 为扁疤深度。
在这里插入图片描述

Matlab代码说明

根据科学出版社出版的《车辆－轨道耦合动力学（第四版）》建立完常规的车轨耦合动力学模型后，只需在时间遍历循环中的最开头添加以下语句即可。

if 0<=mod(vcar*t, 2*pi*R_all) && mod(vcar*t, 2*pi*R_all)<=L
    delta_r = h - ((mod(vcar*t, 2*pi*R_all)-L/2)^2)/2/R_all;
else
    delta_r = 0;
end
R = R_all - delta_r;

其中，vcar表示车辆运行速度，R_all表示车轮圆周半径。
另外，在计算轮轨激励时也要考虑车轮直径的变化：

if z(i, 33)-Zr1-delta_r > 0
    NLz(1) = ((z(i, 33)-Zr1-delta_r)/G)^1.5; NRz(1) = ((z(i, 33)-Zr1-delta_r)/G)^1.5;
else
    NLz(1) = 0; NRz(1) = 0;
end
if z(i, 43)-Zr2-delta_r > 0
    NLz(2) = ((z(i, 43)-Zr2-delta_r)/G)^1.5; NRz(2) = ((z(i, 43)-Zr2-delta_r)/G)^1.5;
else
    NLz(2) = 0; NRz(2) = 0;
end
if z(i, 53)-Zr3-delta_r > 0
    NLz(3) = ((z(i, 53)-Zr3-delta_r)/G)^1.5; NRz(3) = ((z(i, 53)-Zr3-delta_r)/G)^1.5;
else
    NLz(3) = 0; NRz(3) = 0;
end
if z(i, 63)-Zr4-delta_r > 0
    NLz(4) = ((z(i, 63)-Zr4-delta_r)/G)^1.5; NRz(4) = ((z(i, 63)-Zr4-delta_r)/G)^1.5;
else
    NLz(4) = 0; NRz(4) = 0;
end