Simulink：车辆部件温度模型-建模思路

1 博客内容

      Simulink搭建模型的一个细节，用到的时候记不清。这里文字记录，方便同行的人参考。模型搭建的原理比较简单：输入功率=速差*扭矩/9550、流体润滑油散热功率=C₁m₁△T₁、部件材料吸收热量=C₂m₂△T₂，适合入门阶段了解。
0. 模型总体概览
      模型由信号输入、过程计算、输出三部分组成。

2 信号输入

• 使用Subsystem新建输入、过程和输出模块。

• 使用Source→Signal Builder导入表格文件。



• 表格数据如下：

Time	T_OilPan	P_Clu1	P_Clu2	N_Input	N_Clu1	N_Clu2
0	90	4.5	0	1500	1000	0
1	90	4.5	0	1500	1100	0
1.5	90	4.5	0	1500	1200	0
2	90	4.5	0	1500	1300	0
2.5	90	4.5	0	1500	1400	0
3	90	4.5	0	1500	1500	0
3.5	90	4.5	0	1500	1500	0
4	90	4.5	0	1500	1500	0
4.5	90	4.5	0	1500	1500	0
5	90	4.5	0	1500	1500	0
5.5	90	4.5	0	1500	1500	0
6	90	4.5	0	1500	1500	0
6.5	90	4.5	0	1500	1500	0
7	90	4.5	0	1500	1500	0
7.5	90	4.5	0	1500	1500	0
8	90	4.5	0	1500	1500	0
8.5	90	4.5	0	1500	1500	0
9	90	4.5	0	1500	1500	0
9.5	90	4.5	0	1500	1500	0
10	90	4.5	0	1500	1500	0
10.5	90	4.5	0	1500	1000	0
11	90	4.5	0	1500	1100	0
11.5	90	4.5	0	1500	1200	0
12	90	4.5	0	1500	1300	0
12.5	90	4.5	0	1500	1400	0
13	90	4.5	0	1500	1500	0
13.5	90	4.5	0	1500	1500	0
14	90	4.5	0	1500	1500	0
14.5	90	4.5	0	1500	1500	0
15	90	4.5	0	1500	1500	0
15.5	90	4.5	0	1500	1500	0
16	90	4.5	0	1500	1500	0
16.5	90	4.5	0	1500	1500	0
17	90	4.5	0	1500	1500	0
17.5	90	4.5	0	1500	1500	0
18	90	4.5	0	1500	1500	0
18.5	90	4.5	0	1500	1500	0
19	90	4.5	0	1500	1500	0
19.5	90	4.5	0	1500	1500	0
20	90	4.5	0	1500	1500	0
20.5	90	4.5	0	1500	1000	0
21	90	4.5	0	1500	1100	0
21.5	90	4.5	0	1500	1200	0
22	90	4.5	0	1500	1300	0
22.5	90	4.5	0	1500	1400	0
23	90	4.5	0	1500	1500	0
23.5	90	4.5	0	1500	1500	0
24	90	4.5	0	1500	1500	0
24.5	90	4.5	0	1500	1500	0
25	90	4.5	0	1500	1500	0
25.5	90	4.5	0	1500	1500	0
26	90	4.5	0	1500	1500	0
26.5	90	4.5	0	1500	1500	0
27	90	4.5	0	1500	1500	0
27.5	90	4.5	0	1500	1500	0
28	90	4.5	0	1500	1500	0
28.5	90	4.5	0	1500	1500	0
29	90	4.5	0	1500	1500	0
29.5	90	4.5	0	1500	1500	0
30	90	4.5	0	1500	1500	0
30.5	90	4.5	0	1500	1500	0
31	90	4.5	0	1500	1500	0
31.5	90	4.5	0	1500	1500	0
32	90	4.5	0	1500	1500	0
32.5	90	4.5	0	1500	1500	0
33	90	4.5	0	1500	1500	0
33.5	90	4.5	0	1500	1500	0
34	90	4.5	0	1500	1500	0
34.5	90	4.5	0	1500	1500	0
35	90	4.5	0	1500	1500	0
35.5	90	4.5	0	1500	1500	0
36	90	4.5	0	1500	1500	0
36.5	90	4.5	0	1500	1500	0
37	90	4.5	0	1500	1500	0
37.5	90	4.5	0	1500	1500	0
38	90	4.5	0	1500	1500	0
38.5	90	4.5	0	1500	1500	0
39	90	4.5	0	1500	1500	0
39.5	90	4.5	0	1500	1500	0
40	90	4.5	0	1500	1500	0

• 使用goto、from、Gain、Outport库建立输入：Signal_Input。

• 使用constant模块增加两个部件Kp输入，由8个信号构成Input。

• 新建脚本Parameters_Script.m，内容如下：

%单位转换
bar2Pa = 1.0e+5;
Pa2bar = 1.0e-5;
Pa2kPa = 1.0e-3;
kW2W = 1.0e+3;
W2kW = 1.0e-3;
cc2L = 1.0e-3;

%限值
OilTemperature_Lowlimit = -40;
OilTemperature_Uplimit = 150;

P_Clu_Lowlimit = 0;
P_Clu_Uplimit = 1000000;

N_Spd_Lowlimit=0;
N_Spd_Uplimit=6000;

%Friction_table
 Friction_3d_map = zeros(6,13,5);
 Friction_3d_map(:,:,1) = [ 0 59 84 109 134 159 184 209 234 259 284 309 334;
0 61 86 111 136 161 186 211 236 261 286 311 336;
0 63 88 113 138 163 188 213 238 263 288 313 338;
0 69 94 119 144 169 194 219 244 269 294 319 344;
0 79 104 129 154 179 204 229 254 279 304 329 354;
0 101 126 151 176 201 226 251 276 301 326 351 376];
 Friction_3d_map(:,:,2) = [ 0 44 69 94 119 144 169 194 219 244 269 294 319;
0 46 71 96 121 146 171 196 221 246 271 296 321;
0 48 73 98 123 148 173 198 223 248 273 298 323;
0 54 79 104 129 154 179 204 229 254 279 304 329;
0 64 89 114 139 164 189 214 239 264 289 314 339;
0 86 111 136 161 186 211 236 261 286 311 336 361];
 Friction_3d_map(:,:,3) = [ 0 34 59 84 109 134 159 184 209 234 259 284 309;
0 36 61 86 111 136 161 186 211 236 261 286 311;
0 38 63 88 113 138 163 188 213 238 263 288 313;
0 44 69 94 119 144 169 194 219 244 269 294 319;
0 54 79 104 129 154 179 204 229 254 279 304 329;
0 76 101 126 151 176 201 226 251 276 301 326 351];
 Friction_3d_map(:,:,4) = [ 0 28 53 78 103 128 153 178 203 228 253 278 303;
0 30 55 80 105 130 155 180 205 230 255 280 305;
0 32 57 82 107 132 157 182 207 232 257 282 307;
0 38 63 88 113 138 163 188 213 238 263 288 313;
0 48 73 98 123 148 173 198 223 248 273 298 323;
0 70 95 120 145 170 195 220 245 270 295 320 345];
 Friction_3d_map(:,:,5) = [ 0 25 50 75 100 125 150 175 200 225 250 275 300 ;
0 27 52 77 102 127 152 177 202 227 252 277 302;
0 29 54 79 104 129 154 179 204 229 254 279 304; 
0 35 60 85 110 135 160 185 210 235 260 285 310;
0 45 70 95 120 145 170 195 220 245 270 295 320;
0 67 92 117 142 167 192 217 242 267 292 317 342];

Friction_3d_x =[ 0 50 100 150 200 250 300 350 400 450 500 550 600];
Friction_3d_y =[ 0 100 200 500 1000 2000];
Friction_3d_z =[ -40,0,40,100,140];

• 将脚本与模型关联。

3 过程计算

      过程计算包含：输入信号限值过滤，部件1和2的滑摩扭矩计算，部件最大流量、请求流量和实际流量计算、部件内材料温度和流体油出温度计算。

• 信号过滤如下：

• 计算输入功率。使用nD-Lookup，通过转速差、压力（实际压力-Kp点）和油温三个输入得到部件滑摩扭矩传。使用公式：功率=扭矩×速差/9550，计算输入功率。

• 部件最大输入流量（系统输入流量-泄漏量）。其中系统输入流量：转速×排量×效率，温度对泄漏量影响：油底壳温度×系数，主线压力（最小7.5bar）对泄漏量影响：（主线压力：部件需求压力+1.5bar）×系数。

• 部件实际流量。通过输入功率查表得出的需求流量和系统提供的最大流量对比，得出部件实际流量。

• 部件材料温度和流体油出油温度计算。
        （1）流体油出油温度计算，设定输入功率70%由流体油携带，通过Q_1=C1_m1_△T1~计算油温增加量。
        （2）部件材料温度计算，设定输入功率30%用于材料温升(需要使用积分和积分清零)，通过Q₂=C₂m₂△T₂计算钢片温度增加量。

4 信号输出

• 通过Scope（直接观测）和To Workspace（放入工作区）观察输出信号。

5 运行结果

• 设定步长，运行模型。

• 作图，观察变量。
  
        大部分数据主观臆定，部分数据参考维普论文。