?作者簡介熱愛科研的Matlab仿真開發(fā)者擅長數(shù)學建模、數(shù)據(jù)處理、建模仿真、程序設計、完整代碼獲取、論文復現(xiàn)及科研仿真。 往期回顧關注個人主頁Matlab科研工作室 關注我領取海量matlab電子書和數(shù)學建模資料個人信條格物致知,完整Matlab代碼獲取及仿真咨詢內(nèi)容私信。 內(nèi)容介紹多自由度受迫振動系統(tǒng)廣泛存在于機械結構、建筑框架、車輛懸架、精密儀器隔振平臺等工程場景中其在外激勵作用下的穩(wěn)態(tài)響應幅值與共振特性直接決定了結構的安全裕度與工作精度。與單自由度系統(tǒng)僅存在單一共振峰不同多自由度系統(tǒng)的響應曲線會出現(xiàn)多個共振峰且各階模態(tài)的貢獻隨激勵頻率、阻尼分布、參數(shù)偏差而劇烈變化單一工況分析難以反映系統(tǒng)的真實動力學行為。本研究以典型兩自由度與三自由度受迫振動系統(tǒng)為對象建立質(zhì)量-彈簧-阻尼矩陣形式的動力學方程采用模態(tài)疊加法與新mark-β逐步積分法求解穩(wěn)態(tài)與瞬態(tài)響應圍繞激勵頻率、阻尼比、剛度比、質(zhì)量比四類關鍵參數(shù)開展系統(tǒng)性參數(shù)掃描繪制幅頻響應族曲線與峰值分布圖譜定量揭示共振峰位置遷移、峰值幅值演化、模態(tài)耦合強度變化與共振點分裂的規(guī)律。研究結果表明在雙自由度系統(tǒng)中當兩階固有頻率比接近1時會出現(xiàn)拍現(xiàn)象與響應放大當激勵頻率跨越兩階固有頻率之間的反共振點時主質(zhì)量響應會出現(xiàn)深度谷值可用于設計動力吸振器阻尼比對第一共振峰峰值的抑制效率明顯高于第二共振峰當阻尼比從0.02增大至0.20時第一共振峰峰值下降約78%而第二共振峰峰值僅下降約52%剛度參數(shù)掃描顯示當?shù)诙偠扰c第一剛度之比在0.8~1.5區(qū)間時兩階模態(tài)耦合最強共振峰密集且響應峰值最高。本研究形成的參數(shù)掃描流程與圖譜可為多自由度系統(tǒng)的減振設計、參數(shù)靈敏度分析提供量化依據(jù)。?關鍵詞?多自由度系統(tǒng)受迫振動參數(shù)掃描幅頻響應共振峰模態(tài)耦合一、引言工程結構在絕大多數(shù)情況下都不能簡化為單自由度系統(tǒng)。建筑框架的層間振動、車輛懸架的多體耦合、精密設備的減振平臺、航空發(fā)動機的轉(zhuǎn)子-支承系統(tǒng)其動力學行為都需要用多自由度模型描述。多自由度受迫振動系統(tǒng)的核心特征在于其響應特性由系統(tǒng)的固有頻率、振型與阻尼分布共同決定而這三者都是系統(tǒng)物理參數(shù)的函數(shù)。在實際工程中系統(tǒng)參數(shù)往往存在制造誤差、溫度漂移、磨損老化、負載變化等帶來的不確定性。一個在設計工況下表現(xiàn)良好的結構在參數(shù)發(fā)生小幅偏移后可能出現(xiàn)共振峰位置移動、峰值幅值急劇放大甚至出現(xiàn)新的共振點從而引發(fā)疲勞破壞或精度喪失。因此僅分析標稱參數(shù)下的響應是不夠的必須開展參數(shù)掃描全面掌握響應隨參數(shù)變化的演化規(guī)律識別敏感參數(shù)區(qū)間與危險工況。參數(shù)掃描研究具有三個層面的價值其一明確共振峰隨參數(shù)變化的遷移規(guī)律為避振設計提供頻率規(guī)劃依據(jù)其二識別參數(shù)敏感區(qū)間為加工精度分配與公差設計提供量化要求其三揭示反共振現(xiàn)象與模態(tài)耦合規(guī)律為動力吸振器等被動減振裝置的設計提供理論支撐。本報告以典型雙自由度與三自由度受迫振動系統(tǒng)為對象建立完整的動力學模型與求解方法圍繞激勵頻率、阻尼比、剛度比、質(zhì)量比四類參數(shù)開展系統(tǒng)掃描繪制響應圖譜并量化分析各參數(shù)的靈敏度形成可供工程設計參考的結論。二、多自由度受迫振動系統(tǒng)的建模與求解2.1 動力學方程建立?? 運行結果 部分代碼% modal_analysis.m%{Calculates natural (undamped) frequencies, mode shapes and decay times of each rack level- Natural frequencies: undamped resonant frequencies of the rack.- Mode shapes: relative angular displacements (of the displacementwaveforms) of each level, for each natural frequency.- eg: L1 -0.2 rad, L2 0 rad, L 0.2 rad means L1 and L3 are moving in opposite directions while L2 is stationary.- Horizontally arranged vectors, where each column corresponds to aspecific natural frequency.- Decay time τ: time taken for natural response to decay to e^-1 oforiginal magnitude%}function [phi, omega, tau] modal_analysis(input_params)[phi, omega2] eig(input_params.K, input_params.M);omega diag(sqrt(omega2));% calculate time constant tau (time taken for ~63.2% decay)tau 1 ./ (input_params.p_run.zeta * omega);end 參考文獻往期回顧掃掃下方二維碼