公路车桥耦合振动程序（考虑路面不平整度），适用于研究生参考 有两套模型，一套纯ansys，一套是ansys与matlab联合的 模型一： 基于MATLAB与ANSYS软件联合开发的通用性车桥耦合程序，基础车辆模型为双轴四自由度1/2车辆模型，桥梁为拱桥模型，可自由更改车辆模型与桥梁模型参数，十分简便。 程序运行基础步骤如下： 1、更改MATLAB与ANSYS软件运行目录于同一路径下（同一文件夹下） 2、修改文件中车辆参数（车体质量、车轮质量、各弹簧刚度、各阻尼等），修改桥梁模型（桥梁模型需为ANSYS模型） 3、建立加载节点way.txt文件，放入运行文件夹中 4、 建立路面不平整度 5、建立result.txt，内容为桥梁模型中需要提取计算结果的节点编号，此文本输入格式为：第一行为需要提取数据的节点的总数，从第二行开始依次为各节点的编号，节点编号的总数应与第一个数据相同。 6、运行主程序 模型二： 此套模型纯基于ansys软件进行模拟车桥耦合，将原matlab中对于车辆计算部分转移到ansys中，避免了手算车辆运动方程错误。 程序运行基础步骤如下： 1、选择所需车辆模型，输入自身车辆参数 2、修改自身桥梁模型 3、建立路面不平整度 4、运行程序 两套系统均可导出主要结果有：1、车辆四自由度的位移、速度、加速度时程曲线；2、车辆行驶过程中，桥梁各点位移、速度、加速度时程曲线；3、桥梁ANSYS有限元各种计算结果。 下图为参考某篇文献跨中位移时程曲线与软件运行结果对比

最近在研究车桥耦合振动的问题，发现了一个挺有意思的东西，就是利用ANSYS和MATLAB联合开发的车桥耦合振动程序。这个程序不仅能模拟车辆在桥梁上的动态响应，还能考虑路面不平整度的影响，对研究生来说是个不错的参考工具。

两套模型，各有千秋

首先，这个程序有两套模型，一套是纯ANSYS的，另一套是ANSYS与MATLAB联合开发的。我先简单介绍一下它们的特点。

模型一：MATLAB + ANSYS 联合开发

这套模型的核心是基于MATLAB和ANSYS的联合开发，车辆模型采用的是双轴四自由度的1/2车辆模型，桥梁模型则是拱桥模型。它的优点是灵活性很高，可以自由修改车辆和桥梁的参数，非常适合用来做参数化研究。

运行这套程序的步骤如下：

1. 设置运行路径：把MATLAB和ANSYS的运行目录设置到同一个文件夹下，这样程序运行起来会更顺畅。
2. 修改参数：在文件中修改车辆参数（比如车体质量、车轮质量、弹簧刚度、阻尼等）和桥梁模型（桥梁模型需要是ANSYS模型）。
3. 加载节点文件：创建一个way.txt文件，里面写上加载节点的信息，然后把这个文件放到运行文件夹里。
4. 生成路面不平整度：这个步骤有点关键，因为路面不平整度直接影响车辆的动态响应。
5. 设置结果提取节点：创建一个result.txt文件，里面写上需要提取计算结果的桥梁节点编号。
6. 运行主程序：最后运行主程序，就可以得到车辆和桥梁的动力学响应了。

模型二：纯 ANSYS 模拟

这套模型完全基于ANSYS软件，把原来MATLAB中车辆计算的部分转移到了ANSYS中，这样就避免了手算车辆运动方程的麻烦。

运行步骤相对简单：

1. 选择车辆模型：根据自己的需求选择车辆模型，并输入车辆参数。
2. 修改桥梁模型：同样需要修改桥梁模型，以适应不同的研究需求。
3. 生成路面不平整度：和模型一一样，路面不平整度是关键因素。
4. 运行程序：直接运行程序，得到结果。

程序输出结果

两套模型的输出结果都挺全面的，主要有以下几类：

1. 车辆动力学响应：包括车辆四自由度的位移、速度、加速度时程曲线。
2. 桥梁动力学响应：车辆行驶过程中，桥梁各点的位移、速度、加速度时程曲线。
3. 桥梁有限元结果：ANSYS计算的桥梁有限元结果，比如应力、应变等。

代码分析

这里简单分析一下模型一的代码，主要是MATLAB部分。代码的大致结构如下：

% 读取车辆参数
vehicle_params = readVehicleParams('vehicle_config.txt');

% 读取桥梁模型
bridge_model = readBridgeModel('bridge_model.txt');

% 生成路面不平整度
road_roughness = generateRoadRoughness(vehicle_params, bridge_model);

% 设置加载节点
load_nodes = readLoadNodes('way.txt');

% 设置结果提取节点
result_nodes = readResultNodes('result.txt');

% 运行主程序
[vehicle_response, bridge_response] = runCouplingSimulation(vehicle_params, bridge_model, road_roughness, load_nodes, result_nodes);

这段代码的核心是runCouplingSimulation函数，它负责耦合振动的计算。代码逻辑清晰，参数化设计做得很好，方便用户修改参数。

实际应用案例

假设我们要研究不同路面不平整度对桥梁振动的影响，可以按照以下步骤操作：

1. 设置路面不平整度：比如生成一个随机的路面不平整度时程。
2. 运行程序：分别用两套模型运行程序，得到车辆和桥梁的响应。
3. 对比结果：比如对比两套模型的桥梁跨中位移时程曲线。

下图是一个参考文献中的跨中位移时程曲线与程序运行结果的对比：

公路车桥耦合振动程序（考虑路面不平整度），适用于研究生参考 有两套模型，一套纯ansys，一套是ansys与matlab联合的 模型一： 基于MATLAB与ANSYS软件联合开发的通用性车桥耦合程序，基础车辆模型为双轴四自由度1/2车辆模型，桥梁为拱桥模型，可自由更改车辆模型与桥梁模型参数，十分简便。 程序运行基础步骤如下： 1、更改MATLAB与ANSYS软件运行目录于同一路径下（同一文件夹下） 2、修改文件中车辆参数（车体质量、车轮质量、各弹簧刚度、各阻尼等），修改桥梁模型（桥梁模型需为ANSYS模型） 3、建立加载节点way.txt文件，放入运行文件夹中 4、 建立路面不平整度 5、建立result.txt，内容为桥梁模型中需要提取计算结果的节点编号，此文本输入格式为：第一行为需要提取数据的节点的总数，从第二行开始依次为各节点的编号，节点编号的总数应与第一个数据相同。 6、运行主程序 模型二： 此套模型纯基于ansys软件进行模拟车桥耦合，将原matlab中对于车辆计算部分转移到ansys中，避免了手算车辆运动方程错误。 程序运行基础步骤如下： 1、选择所需车辆模型，输入自身车辆参数 2、修改自身桥梁模型 3、建立路面不平整度 4、运行程序 两套系统均可导出主要结果有：1、车辆四自由度的位移、速度、加速度时程曲线；2、车辆行驶过程中，桥梁各点位移、速度、加速度时程曲线；3、桥梁ANSYS有限元各种计算结果。 下图为参考某篇文献跨中位移时程曲线与软件运行结果对比

!对比图

从图中可以看出，程序的计算结果与文献中的结果吻合得不错，说明程序的准确性是有保障的。

总结

总的来说，这套车桥耦合振动程序无论是从功能还是灵活性上，都挺适合研究生用来做相关研究的。两套模型各有特点，可以根据自己的需求选择适合的方案。如果想更深入地研究车辆动力学，模型一可能更适合；如果更关注桥梁的响应，模型二可能更方便。

希望这篇博文能帮到正在研究车桥耦合振动的同学们！如果有任何问题，欢迎留言讨论~