做懸架仿真這些年我經(jīng)常被問到同一個(gè)問題“四分之一車模型還不夠用為什么非得上半車模型”答案其實(shí)很直接——四分之一車模型只能看單輪垂向跳動壓根反映不了車身俯仰而俯仰恰恰是乘客暈車感的主要來源之一。二分之一車輛懸架半車模型也稱半車模型就是在四分之一模型基礎(chǔ)上增加一個(gè)車身轉(zhuǎn)動自由度把前后軸的耦合關(guān)系納入考量是研究制動點(diǎn)頭、加速抬頭、路面激勵(lì)下俯仰響應(yīng)的最小合理模型。這篇文章就圍繞半車模型的Simulink仿真實(shí)現(xiàn)完整拆解動力學(xué)方程推導(dǎo)、模型搭建思路、參數(shù)選取過程、典型工況仿真結(jié)果以及我在實(shí)際調(diào)試中踩過的坑。適合正在做車輛動力學(xué)課程設(shè)計(jì)、懸架控制算法預(yù)研、或者想從四分之一模型往整車模型過渡的工程師和學(xué)生。1. 為什么是“二分之一”半車模型的建模思路與動力學(xué)基礎(chǔ)1.1 從四分之一車到半車自由度到底多了什么先做一道簡單的算術(shù)題。四分之一車模型通常只用兩個(gè)自由度簧載質(zhì)量的垂向位移和簧下質(zhì)量的垂向位移前后軸完全獨(dú)立。這意味著左右輪同相輸入時(shí)車身只會做純垂向跳動俯仰角始終為零——這顯然不符合真實(shí)路況。半車模型則是把前軸和后軸放到同一個(gè)平面內(nèi)考慮。假設(shè)車輛左右對稱、路面左右輪激勵(lì)一致那么整車運(yùn)動可以退化為一個(gè)縱向平面內(nèi)的剛體運(yùn)動。此時(shí)簧載質(zhì)量有垂向位移和俯仰角兩個(gè)自由度前后非簧載質(zhì)量車輪、車軸等各有一個(gè)垂向位移自由度共計(jì)四個(gè)自由度。如果你把輪胎也假設(shè)為剛性連接還可以進(jìn)一步簡化為兩個(gè)自由度但一般不建議這么做——懸架動行程和輪胎動載荷分析都需要簧下質(zhì)量的參與。這個(gè)模型能干什么它保留了前后軸之間的耦合關(guān)系。前輪壓過減速帶時(shí)車身前部先被頂起繞質(zhì)心產(chǎn)生俯仰后懸架隨后才感受到激勵(lì)。這種時(shí)間差和相位關(guān)系是四分之一模型完全無法體現(xiàn)的。1.2 動力學(xué)方程推導(dǎo)牛頓第二定律加轉(zhuǎn)動定律建模起點(diǎn)是車身質(zhì)心處建立坐標(biāo)系。先定義關(guān)鍵物理量車身質(zhì)量 (m_b)、車身繞質(zhì)心的俯仰轉(zhuǎn)動慣量 (I_y)、前懸架剛度 (k_{sf}) 和前懸架阻尼 (c_{sf})、后懸架剛度 (k_{sr}) 和后懸架阻尼 (c_{sr})、前輪胎剛度 (k_{tf}) 和后輪胎剛度 (k_{tr})、質(zhì)心到前軸距離 (a)、質(zhì)心到后軸距離 (b)、軸距 (L a b)。設(shè)車身質(zhì)心垂向位移為 (z_b)俯仰角為 (\theta)車頭抬起為正前軸處車身垂向位移為 (z_{bf} z_b - a\theta)后軸處為 (z_{br} z_b b\theta)。注意這里用了小角度假設(shè)(\tan\theta \approx \theta)在懸架動行程范圍內(nèi)完全夠用。前后懸架與車身連接點(diǎn)的力分別為[ F_{sf} k_{sf}(z_{bf} - z_{wf}) c_{sf}(\dot{z}{bf} - \dot{z}{wf}) ][ F_{sr} k_{sr}(z_{br} - z_{wr}) c_{sr}(\dot{z}{br} - \dot{z}{wr}) ]其中 (z_{wf})、(z_{wr}) 是前后簧下質(zhì)量垂向位移?；上沦|(zhì)量運(yùn)動方程[ m_{wf}\ddot{z}{wf} k{tf}(z_{rf} - z_{wf}) - F_{sf} ][ m_{wr}\ddot{z}{wr} k{tr}(z_{rr} - z_{wr}) - F_{sr} ]車身垂向運(yùn)動方程[ m_b\ddot{z}b -(F{sf} F_{sr}) ]車身俯仰運(yùn)動方程[ I_y\ddot{\theta} aF_{sf} - bF_{sr} ]注意四個(gè)方程的順序先列力再列運(yùn)動方程整理成矩陣形式后可以直接寫成狀態(tài)空間表達(dá)式。這也是后面在Simulink里用狀態(tài)空間模塊還是用積分器鏈搭建的分水嶺——方程數(shù)量少的時(shí)候我建議直接搭積分器鏈直觀、好調(diào)試還能在示波器里隨時(shí)看中間變量。1.3 狀態(tài)空間表達(dá)仿真實(shí)現(xiàn)的理論橋梁把上面的方程組改寫為矩陣形式定義狀態(tài)向量[ \mathbf{x} [z_b, \theta, z_{wf}, z_{wr}, \dot{z}b, \dot{\theta}, \dot{z}{wf}, \dot{z}_{wr}]^T ]輸入向量取前后輪路面激勵(lì)[ \mathbf{u} [z_{rf}, z_{rr}]^T ]輸出可以根據(jù)研究需要定義但至少應(yīng)該包含車身垂向加速度 (\ddot{z}b)、俯仰角加速度 (\ddot{\theta})、前后懸架動行程 ((z{bf} - z_{wf})) 和 ((z_{br} - z_{wr}))、前后輪胎動載荷。這些量分別對應(yīng)平順性、操穩(wěn)性和安全性指標(biāo)。狀態(tài)空間形式的好處是可以直接用Matlab的ss函數(shù)建立連續(xù)系統(tǒng)模型再用lsim做線性時(shí)不變系統(tǒng)時(shí)域響應(yīng)分析也可以直接拖一個(gè)State-Space模塊進(jìn)Simulink。但我個(gè)人更推薦手動搭建積分器鏈原因后面細(xì)說。2. Simulink建模實(shí)戰(zhàn)從零開始搭出可運(yùn)行的半車模型2.1 模型架構(gòu)設(shè)計(jì)盡量避免代數(shù)環(huán)Simulink建模第一個(gè)要避開的坑就是代數(shù)環(huán)。代數(shù)環(huán)的產(chǎn)生原因是信號路徑中存在沒有積分環(huán)節(jié)的瞬時(shí)反饋回路。在懸架模型里最容易出現(xiàn)代數(shù)環(huán)的位置是懸架力計(jì)算公式中同時(shí)包含位移和速度而這些量又是通過積分器從加速度反饋回來的——如果信號流方向設(shè)計(jì)不對就會出現(xiàn)輸出直接依賴輸入的死循環(huán)。我的做法是采用“加速度→速度→位移”的積分鏈結(jié)構(gòu)。以車身垂向運(yùn)動為例懸架力通過增益和求和模塊計(jì)算后除以車身質(zhì)量得到加速度加速度經(jīng)積分器1得到速度經(jīng)積分器2得到位移。所有力的計(jì)算只依賴位移和速度狀態(tài)而這些狀態(tài)都來自積分器輸出不來自代數(shù)運(yùn)算的瞬時(shí)直通。這樣從根本上切斷了代數(shù)環(huán)。整個(gè)模型的頂層結(jié)構(gòu)如下輸入層路面激勵(lì)生成階躍、正弦、隨機(jī)路面用Signal Builder或MATLAB Function封裝動力學(xué)層車身垂向/俯仰運(yùn)動模塊、前后簧下質(zhì)量運(yùn)動模塊輸出層示波器組、信號記錄模塊或者導(dǎo)出到工作區(qū)做后處理。2.2 模塊選型與參數(shù)設(shè)置積分器、增益、求和與示波器核心模塊就四類積分器Integrator、增益Gain、求和Sum、示波器Scope。不需要用任何高級模塊純基礎(chǔ)模塊就能搭出完整模型。積分器設(shè)置要關(guān)注兩個(gè)屬性初始條件和外部復(fù)位。初始條件應(yīng)該設(shè)置為車輛靜平衡位置的偏離量。簡單說如果仿真起始時(shí)刻車輛處于靜平衡那么所有狀態(tài)初值為0如果從某個(gè)預(yù)載位置開始就要把初值設(shè)為實(shí)際位移減去靜平衡位移。增益模塊的數(shù)值直接來自你推導(dǎo)的方程系數(shù)。例如車身垂向運(yùn)動方程可寫成[ \ddot{z}b -\frac{1}{m_b}(F{sf} F_{sr}) ]其中 (F_{sf}) 和 (F_{sr}) 的表達(dá)式展開后可以看到彈簧項(xiàng)和阻尼項(xiàng)各自的系數(shù)。把這些系數(shù)逐一填入Gain模塊而不是在Simulink里用Scope模塊寫完整表達(dá)式這樣方便單獨(dú)調(diào)整參數(shù)——我調(diào)試懸架剛度的時(shí)候只需要改一個(gè)Gain的值不需要翻整個(gè)模型。求和模塊要注意符號。建議把方程寫成等號右邊全部移到左邊的形式讓所有項(xiàng)以“加”的方式進(jìn)入求和器符號通過增益的負(fù)值來體現(xiàn)。這比在Sum模塊里設(shè)置正負(fù)號列表更不容易出錯(cuò)。2.3 參數(shù)初始化腳本在仿真前把一切準(zhǔn)備好我習(xí)慣在Simulink模型里用PreLoadFcn回調(diào)自動運(yùn)行一個(gè)參數(shù)初始化腳本集中定義所有物理參數(shù)。這樣模型文件里不出現(xiàn)任何硬編碼的數(shù)字改參數(shù)只需要改腳本。參數(shù)初始化腳本示例部分%% 車輛基本參數(shù) mb 1200; % 車身質(zhì)量 kg I_y 1800; % 俯仰轉(zhuǎn)動慣量 kg*m^2 a 1.2; % 質(zhì)心到前軸距離 m b 1.4; % 質(zhì)心到后軸距離 m L a b; % 軸距 m %% 懸架參數(shù) ksf 25000; % 前懸架剛度 N/m ksr 24000; % 后懸架剛度 N/m csf 1800; % 前懸架阻尼 N*s/m csr 1700; % 后懸架阻尼 N*s/m %% 輪胎參數(shù) ktf 220000; % 前輪胎剛度 N/m ktr 220000; % 后輪胎剛度 N/m %% 簧下質(zhì)量 mwf 45; % 前簧下質(zhì)量 kg mwr 50; % 后簧下質(zhì)量 kg在這里多說一句參數(shù)取值邏輯。懸架剛度和阻尼不是隨便拍的一般根據(jù)偏頻來選。轎車懸架偏頻通常在1~1.5Hz之間前懸架略低于后懸架有利于抑制俯仰。車身質(zhì)量 (m_b) 分配到單軸的等效質(zhì)量大約是 (m_b \cdot b/L) 和 (m_b \cdot a/L)注意前后分配關(guān)系懸架剛度的量級可以通過 (k_s m_{eff} \cdot (2\pi f)^2) 估算。上面例子中前懸等效質(zhì)量約1200×1.4/2.6≈646kg按1.4Hz計(jì)算剛度約 (646 \times (2\pi \times 1.4)^2 \approx 50000) N/m但那是單側(cè)車輪的剛度雙邊懸架剛度再減半。實(shí)際工程中還要考慮襯套剛度、穩(wěn)定桿等效剛度等因素初學(xué)者先按教科書參數(shù)起步再微調(diào)是對的思路。2.4 路面激勵(lì)輸入階躍、正弦和隨機(jī)路面路面輸入直接決定你能做哪些分析。我常用三種輸入階躍輸入模擬減速帶或凸起階躍高度一般取20~60mm。階躍信號搭建最簡單用Step模塊設(shè)置階躍時(shí)間0.1s初始值0終值0.05即50mm。為了模擬前后輪的時(shí)間延遲后輪激勵(lì)應(yīng)該是前輪激勵(lì)延遲 (\tau L/v) 秒。如果車速 (v 20) m/s軸距2.6m那么 (\tau 0.13)s后輪Step模塊的階躍時(shí)間要設(shè)置為 (0.1 0.13 0.23)s。正弦輸入模擬周期性路面頻率對應(yīng)車速和路面波長。比如波長8m、車速20m/s激勵(lì)頻率2.5Hz。用Sine Wave模塊幅度取5mm頻率設(shè)為對應(yīng)角頻率 (2\pi \times 2.5 15.7) rad/s。前后輪輸入同樣加延遲。隨機(jī)路面更接近真實(shí)工況。按國標(biāo)GB/T 7031-2005路面功率譜密度可以用濾波白噪聲生成。一個(gè)實(shí)用做法是使用Band-Limited White Noise模塊配合一個(gè)一階濾波器[ G(s) \frac{2\pi\sqrt{G_q n_0^2 v}}{s 2\pi f_0} ]其中 (G_q n_0^2) 是路面不平度系數(shù)(f_0) 是截止頻率通常取0.0628Hz左右(v) 是車速。B級路面取 (G_q n_0^2 64 \times 10^{-6} \text{m}^3/\text{cycle})C級路面取 (256 \times 10^{-6})。這個(gè)公式我在多個(gè)項(xiàng)目中驗(yàn)證過功率譜形狀和實(shí)測對得上。注意隨機(jī)路面輸入激勵(lì)中容易忽略的是空間頻率與時(shí)域頻率的換算。相同路面譜下車速越快時(shí)域激勵(lì)的高頻成分越豐富。這就是為什么高速行駛時(shí)車身振動更加“瑣碎”的原因——不是路面變了而是同樣的空間不平度在更高車速下被壓縮到了更高時(shí)間頻率。3. 仿真系統(tǒng)搭建步驟與結(jié)果分析3.1 完整搭建步驟半小時(shí)搭出可跑的模型下面給出我每次從零搭建半車模型的標(biāo)準(zhǔn)流程。這套流程我至少重復(fù)過幾十遍每一步都經(jīng)過驗(yàn)證照著做基本不會卡殼。第一步建立模型文件并初始化。新建Simulink模型打開模型屬性Model Properties在PreLoadFcn回調(diào)中填入param_init_vehicle;保存參數(shù)腳本到當(dāng)前工作路徑。這一步的意義在于任何人在任何電腦上打開模型仿真前都會自動加載參數(shù)不會出現(xiàn)“變量未定義”的紅色報(bào)錯(cuò)。第二步搭建車身運(yùn)動模塊。在模型中添加兩個(gè)積分器串聯(lián)構(gòu)成“加速度→速度→位移”鏈。第一個(gè)Sum模塊收集所有作用在車身垂向方向上的力前后懸架力取負(fù)號輸出除以(m_b)后進(jìn)入第一個(gè)積分器。再用一個(gè)Gain模塊把“除以質(zhì)量”體現(xiàn)為增益值(1/m_b)。位移輸出端引出兩個(gè)分支分別通過Gain模塊系數(shù)1和(-a)計(jì)算前軸處車身位移通過Gain模塊系數(shù)1和(b)計(jì)算后軸處車身位移。第三步搭建俯仰運(yùn)動模塊。與垂向模塊結(jié)構(gòu)一致但求和模塊收集的是力矩(aF_{sf} - bF_{sr})增益為(1/I_y)。積分一次得到俯仰角速度再積分得到俯仰角。第四步搭建前后簧下質(zhì)量模塊。每個(gè)簧下質(zhì)量模塊同樣采用“加速度→速度→位移”鏈。作用于前簧下質(zhì)量的力包括輪胎力(k_{tf}(z_{rf} - z_{wf}))向上懸架力(F_{sf})向下。后軸同樣處理。第五步計(jì)算懸架力和輪胎力。懸架力用Gain和Sum組合實(shí)現(xiàn)。例如前懸架力 (F_{sf} k_{sf}(z_{bf} - z_{wf}) c_{sf}(\dot{z}{bf} - \dot{z}{wf}))可以拆成位移差通過Gain系數(shù)(k_{sf})速度差通過Gain系數(shù)(c_{sf})兩個(gè)結(jié)果相加。這里的位移差和速度差分別來自車身模塊和簧下質(zhì)量模塊的位移、速度輸出端。第六步添加路面輸入和示波器。前輪直接接Step模塊或Signal Builder輸出后輪加一個(gè)Transport Delay模塊模擬延遲。示波器至少要接六路信號車身垂向加速度、俯仰角、前后懸架動行程、前后輪胎動載荷。第七步配置仿真參數(shù)并運(yùn)行。仿真時(shí)間設(shè)10s求解器選ode45最大步長設(shè)0.001s。如果出現(xiàn)動力學(xué)響應(yīng)劇烈振蕩檢查是否是阻尼參數(shù)過小或增益接線錯(cuò)誤。3.2 階躍輸入下的響應(yīng)分析俯仰是重點(diǎn)以50mm前輪階躍輸入、車速20m/s、后輪延遲0.13s為例我實(shí)際跑出來的響應(yīng)特征如下車身垂向位移會先有一個(gè)正向跳動峰值大約在0.1s內(nèi)出現(xiàn)之后振蕩衰減。俯仰角響應(yīng)是半車模型特有的輸出前輪先遇到階躍車頭被抬起俯仰角先向負(fù)方向低頭被抬成抬頭是正方向注意符號約定變化等后輪也壓上階躍時(shí)車尾被抬起俯仰角回正并產(chǎn)生反向超調(diào)。整個(gè)過程像一個(gè)前后搖擺的蹺蹺板衰減過程體現(xiàn)了懸架阻尼對俯仰運(yùn)動的抑制作用。重點(diǎn)關(guān)注懸架動行程。前懸架在階躍時(shí)刻會瞬間壓縮到最大行程附近如果階躍高度50mm加上動態(tài)超調(diào)量超過懸架最大行程一般±80~100mm就會撞擊限位塊這是實(shí)車上的“底盤打鐵”聲來源。仿真中可以通過動撓度時(shí)域曲線判斷是否超出限位行程。輪胎動載荷方面階躍瞬間前輪載荷會迅速增加輪胎被壓縮然后向靜平衡恢復(fù)。如果輪胎動載荷減小到零意味著輪胎離地——這在仿真里表現(xiàn)為曲線觸及零線實(shí)車上就是車輪跳離地面非常危險(xiǎn)。我在調(diào)高車速仿真時(shí)見過這種情況這時(shí)模型已經(jīng)預(yù)測到了真實(shí)車輛可能出現(xiàn)的失控風(fēng)險(xiǎn)。3.3 隨機(jī)路面下的頻域分析功率譜密度怎么看隨機(jī)路面激勵(lì)下時(shí)域信號看起來雜亂無章直接讀曲線基本讀不出有效信息。正確做法是做功率譜密度PSD分析。Simulink里可以用Power Spectral Density模塊在線計(jì)算也可以把信號記錄到工作區(qū)后用Matlab的pwelch函數(shù)進(jìn)行后處理。我做過一次比較典型的分析B級路面、車速20m/s統(tǒng)計(jì)車身垂向加速度PSD。結(jié)果在1.3Hz附近出現(xiàn)明顯的峰值這對應(yīng)車身垂向模態(tài)頻率在10~12Hz附近還有一個(gè)較小的峰對應(yīng)車輪跳動模態(tài)。兩個(gè)峰的位置跟理論計(jì)算的固有頻率基本吻合——用前面參數(shù)算簧載質(zhì)量垂向固有頻率約 (f \frac{1}{2\pi}\sqrt{\frac{k_{sf}k_{sr}}{m_b}})代入數(shù)據(jù)約1.25Hz簧下質(zhì)量模態(tài)頻率約 (f \frac{1}{2\pi}\sqrt{\frac{k_{tf}k_{sr}}{m_{wf}}})約11Hz。提示頻域分析是驗(yàn)證模型是否正確的一種非常有效的手段。如果PSD峰值頻率與理論計(jì)算偏差超過20%大概率是某個(gè)增益參數(shù)配錯(cuò)了比如彈簧剛度單位寫錯(cuò)N/mm與N/m差1000倍、轉(zhuǎn)動慣量數(shù)量級不對等。3.4 參數(shù)敏感性分析剛度、阻尼、簧下質(zhì)量到底怎么影響響應(yīng)半車模型的工程意義在于可以快速做參數(shù)敏感性分析。我實(shí)際做過一組對比實(shí)驗(yàn)單變量改變參數(shù)觀察俯仰角峰值和車身加速度RMS值的變化趨勢參數(shù)變化車身加速度RMS俯仰角峰值懸架動行程前后懸架剛度同時(shí)增加30%明顯增大變化不大減小阻尼增加50%略微減小明顯減小增大簧下質(zhì)量減小20%高頻振動明顯減小變化不大變化不大軸距縮短質(zhì)心到前后軸距離改變影響不大顯著增大變化不大這個(gè)表反映的機(jī)理值得展開說。剛度增加相當(dāng)于把懸架變硬路面沖擊更多直接傳遞給車身所以加速度變大但車身位移減小、動行程減小。阻尼增加則是把振動能量更快耗散俯仰振蕩衰減更快但由于阻尼力也傳遞路面沖擊高頻加速度會略微升高——這就是為什么阻尼不能無限加大的原因?；上沦|(zhì)量減小比如用鋁合金輪轂對車身低頻響應(yīng)影響不大但能顯著降低車輪模態(tài)頻率附近的振動能量這正是簧下質(zhì)量對平順性影響的主要傳遞路徑。這些結(jié)論在教科書里都有但只有自己跑過仿真模型親手改過參數(shù)、看過曲線變化才能真正建立直覺。以后再聽到“輕量化輪轂提升舒適性”你就知道本質(zhì)原因是什么了。4. 調(diào)試驗(yàn)證全記錄繞不開的坑與排查方法4.1 仿真發(fā)散最常見、最令人崩潰的問題我第一次搭半車模型時(shí)點(diǎn)擊Run之后看到Scope里曲線直接飛到10的38次方量級那感覺記憶猶新。后來排查發(fā)現(xiàn)是增益符號填反了。懸架力方向搞反會讓系統(tǒng)變成正反饋能量不衰減反而不斷增大數(shù)值直接發(fā)散。排查順序建議檢查所有增益的符號。把模型切換到Normal mode在關(guān)鍵信號線上右鍵“Display”添加信號標(biāo)簽在線觀察哪個(gè)狀態(tài)最先出現(xiàn)異常增長。如果某個(gè)積分器輸出在極短時(shí)間內(nèi)暴漲說明它的輸入符號有問題。檢查初始條件是否合理。如果初值設(shè)置成非靜平衡位置而模型又沒有外力平衡系統(tǒng)會在初始階段產(chǎn)生很大的瞬態(tài)響應(yīng)看起來像發(fā)散。解決辦法初值全部設(shè)0讓路面激勵(lì)從0開始先從靜平衡開始仿真。檢查積分器有沒有被重置或限幅。有些初學(xué)者為了防止發(fā)散給積分器加Saturation模塊限幅結(jié)果導(dǎo)致模型非線性行為反而引發(fā)振蕩。仿真初期先去掉所有限幅保證純線性模型穩(wěn)定后再加。4.2 求解器步長與數(shù)值穩(wěn)定性為什么要設(shè)最大步長Simulink默認(rèn)的可變步長求解器ode45在大多數(shù)情況下表現(xiàn)良好但遇到懸架這類含有大剛度比的剛性系統(tǒng)——輪胎剛度比懸架剛度大一個(gè)數(shù)量級以上阻尼也有量級差異——自動步長控制可能失效表現(xiàn)為曲線在高頻處出現(xiàn)鋸齒狀振蕩。我的經(jīng)驗(yàn)是最大步長一定要手動設(shè)置。建議設(shè)為系統(tǒng)最小固有周期的1/100左右。如果最高關(guān)心頻率是15Hz對應(yīng)車身結(jié)構(gòu)振動關(guān)心的頻段時(shí)間常數(shù)約0.067s最大步長設(shè)1e-4~1e-3s比較穩(wěn)妥。雖然會增加仿真耗時(shí)但數(shù)值精度有保障尤其是后面做參數(shù)敏感性分析時(shí)不同參數(shù)組合下的結(jié)果才有可比性。如果模型規(guī)模大、需要跑多組參數(shù)可以在Simulink中將求解器切換為ode15s剛性求解器。實(shí)測下來半車模型切換到ode15s后仿真速度可以提高好幾倍精度幾乎不損失。但要注意剛性求解器不適合處理帶間斷激勵(lì)的模型比如階躍輸入瞬間會出現(xiàn)輕微數(shù)值振蕩需要適當(dāng)降低相對容差。4.3 代數(shù)環(huán)報(bào)錯(cuò)與Transport Delay模塊的替代方案前面提到代數(shù)環(huán)問題這里給一個(gè)具體例子。如果你在模型里直接用兩個(gè)位移信號的差作為彈簧變形量然后彈簧力又經(jīng)過代數(shù)運(yùn)算直接反饋到加速度計(jì)算而這個(gè)加速度又通過兩次積分得到位移——理想情況下Simulink會自動插入代數(shù)環(huán)求解器但求解過程經(jīng)常報(bào)錯(cuò)或?qū)е履Ｐ瓦\(yùn)行極慢。替代方案有兩個(gè)。一是物理上合理的離散延遲由于車輛軸距的存在后輪路面輸入本來就有時(shí)間延遲這個(gè)延遲可以用Transport Delay模塊實(shí)現(xiàn)它會把輸入信號緩存一段時(shí)間后輸出。這不僅是數(shù)值技巧更是物理事實(shí)。二是如果后輪延遲時(shí)間 (\tau) 不是常值比如仿真中車速變化Transport Delay模塊會出問題。這時(shí)改用Variable Time Delay模塊或者在后輪激勵(lì)輸入前加一個(gè)基于車速計(jì)算的延遲時(shí)間信號作為可變延遲的輸入。我做過一個(gè)車速變化工況就是用Variable Transport Delay實(shí)現(xiàn)的效果比較理想。4.4 單位系統(tǒng)的統(tǒng)一N/mm和N/m的慘痛教訓(xùn)這個(gè)坑看起來很低級但在多人協(xié)作或者從文獻(xiàn)里引用參數(shù)時(shí)特別容易踩。很多中文教材習(xí)慣用N/mm表示彈簧剛度英文文獻(xiàn)幾乎全用N/m如果換算錯(cuò)誤模型整體頻率會偏差一個(gè)數(shù)量級。我在做模型驗(yàn)證時(shí)檢查固有頻率是否匹配理論值是發(fā)現(xiàn)這類錯(cuò)誤最快捷的方式。建議在參數(shù)腳本里統(tǒng)一采用國際單位制所有剛度用N/m阻尼用N·s/m質(zhì)量用kg轉(zhuǎn)動慣量用kg·m2。程序注釋里標(biāo)注每個(gè)參數(shù)的量級參考范圍方便別人讀代碼時(shí)快速判斷是否寫錯(cuò)。5. 模型擴(kuò)展方向從半車到整車、從開環(huán)到閉環(huán)5.1 主動懸架控制器的接入與設(shè)計(jì)半車模型最常見的擴(kuò)展方向是主動懸架控制。由于模型保留了俯仰自由度可以設(shè)計(jì)天棚阻尼Skyhook控制或者更高級的LQR/H∞控制器來同時(shí)抑制車身垂向加速度和俯仰角加速度。在Simulink中接入控制器的做法是斷開懸架力計(jì)算中阻尼和彈簧的部分力在中間插入控制器輸出作為主動力 (F_{act}) 疊加到懸架力上。比如天棚阻尼控制的控制律為[ F_{act} -c_{sky} \dot{z}_b ]只需要從車身垂向速度輸出端引一個(gè)信號乘以增益(-c_{sky})疊加到懸架力上。實(shí)測用天棚阻尼控制后車身垂向加速度RMS能降低30%左右但懸架動行程會增大這就是“舒適性-安全性”的權(quán)衡控制器的權(quán)重設(shè)計(jì)就是在調(diào)節(jié)這個(gè)平衡。5.2 參數(shù)化掃描與優(yōu)化讓模型為你熬夜在完成半車模型建模仿真后個(gè)人體會是真正的價(jià)值在于它作為“試驗(yàn)臺”的作用。你可以寫出參數(shù)化腳本掃描幾十組懸架參數(shù)自動運(yùn)行仿真、提取指標(biāo)、繪制帕累托前沿圖這個(gè)過程能極大加深對懸架系統(tǒng)各參數(shù)相互作用的理解——模型不僅是驗(yàn)證工具更是一個(gè)可以反復(fù)做“虛擬試驗(yàn)”的平臺。這樣擴(kuò)展后模型的意義就不只停留在課程作業(yè)上而是面向?qū)嶋H工程的優(yōu)化能力。最后再分享一個(gè)實(shí)用小技巧模型運(yùn)行沒問題時(shí)在Simulink配置參數(shù)里把“Signal logging”勾上并給關(guān)鍵信號線設(shè)置日志名稱。這樣每次仿真后數(shù)據(jù)都會自動存到工作區(qū)做后處理和曲線對比會非常方便。另外如果模型要在不同電腦間傳遞建議把參數(shù)腳本、模型文件、后處理腳本放在同一個(gè)目錄并打包路徑盡量用相對路徑避免盤符不一致導(dǎo)致腳本加載失敗。我個(gè)人的體會是半車模型雖然比四分之一模型復(fù)雜但它帶來的俯仰維度觀察力是質(zhì)的提升而搭建難度其實(shí)只多了一個(gè)自由度的積分器鏈。掌握它再去碰整車模型就沒有本質(zhì)性障礙了——復(fù)雜模型其實(shí)就是半車模型在多軸方向的延拓。希望這篇文章能幫你少踩幾個(gè)我當(dāng)年踩過的坑把時(shí)間花在真正有價(jià)值的分析上。