【正文】 雜多自由度系統(tǒng)固有頻率具有參考意義，并為初步分析各構(gòu)件的振動(dòng)情況以及解耦分析系統(tǒng)響應(yīng)奠定了基礎(chǔ)。clc。b(4)=。 fplot(yx,[0,3])。)。clc。b(4)=。 fplot(yx,[0,3])。)。為求式（413）的解，假定 （417）其中，與為常量。亮的固有頻率可由式（415）計(jì)算，即 （423）函數(shù)稱為梁的固有振型函數(shù)，為振動(dòng)的固有頻率。假設(shè)梁的撓度為 （428）其中，為第階固有振型函數(shù)或滿足微分方程的特征函數(shù)：， （429）為對(duì)應(yīng)的廣義坐標(biāo)。a=10。w=100。 w1(i)=(2*f0/(ro*A*1))*sin(n*pi*a/1)*sin(n*pi*x/1)*sin(w*t(i))/(wn^2w^2)。endfor i=1:1001 t(i)=3*(i1)/1000。)。)。subplot(39。ylabel(39。x=20,n=1,239。)。)。title(39。在計(jì)算系統(tǒng)固有頻率與主振型后可以采用MATLAB的方法把抽象的理論圖像數(shù)據(jù)化分析，相對(duì)于以前所采用的算法語(yǔ)言（Fortran、C語(yǔ)言等），本文采用MATLAB語(yǔ)言，不僅相對(duì)語(yǔ)句少，可讀性強(qiáng)，還可以利用MATLAB的繪圖功能對(duì)結(jié)果進(jìn)行直觀地分析。 致 謝四年的本科大學(xué)學(xué)習(xí)生活即將結(jié)束，在完成畢業(yè)論文之際，深感畢業(yè)課題的順利完成是與許多人對(duì)筆者的幫助和關(guān)心分不開的，在此寫下感激之情。對(duì)此，筆者對(duì)方老師表示最誠(chéng)摯的敬意和最衷心的感謝！在課題的實(shí)驗(yàn)研究過(guò)程中，得到了力學(xué)與工程學(xué)院各位領(lǐng)導(dǎo)和老師的大力指導(dǎo)和熱情幫助。衷心感謝張智慧老師對(duì)筆者的大力幫助！幾年來(lái)，父母在生活上和事業(yè)上給予了筆者莫大的支持和鼓勵(lì)，付出了辛勤的勞動(dòng)，謹(jǐn)以此文獻(xiàn)給二老。再一次向曾經(jīng)培養(yǎng)、教育、關(guān)心和幫助過(guò)筆者的領(lǐng)導(dǎo)、前輩、老師和朋友們致以最誠(chéng)摯、最衷心的感謝！參考文獻(xiàn)[1] 許本文,[M].北京:機(jī)械工業(yè)出版社,1988.[2] 方治華,賈宏玉,[J],包頭鋼鐵學(xué)院學(xué)報(bào),2004,12.[3] 陳予怒,非線性振動(dòng),天津科學(xué)技術(shù)出版社[M],1982.[4] 傅志方,華宏星,模態(tài)分析理論與應(yīng)用[M],上海交通大學(xué)出版社,2002.[5] 周傳榮,趙淳生,機(jī)械振動(dòng)參數(shù)識(shí)別及其應(yīng)用[M],科學(xué)出版社,1989.[6] 席平原,張海濤,機(jī)械動(dòng)態(tài)系統(tǒng)可視化仿真及應(yīng)用[J],煤礦機(jī)械,2004,25（8）:3l33.[7] Edward ,高會(huì)生,李新葉,胡智奇等譯,Matlab原理與工程應(yīng)用[J],電子工業(yè)出版社,2002.[8] 張廷中,MATLAB使用指南[M],北京:科學(xué)技術(shù)文獻(xiàn)出版社,1998.[9] 張志涌等,精通MATLAB[M],北京:北京航空航天大學(xué)出版社,2000.[10] Hanselman D,Littlefield [M].張 航（譯）[M]北京:清華大學(xué)出版社,2002.[11] 夏永源,張阿舟機(jī)械振動(dòng)問題的計(jì)算機(jī)解法[M],北京國(guó)防工業(yè)出版社,1993.[12] 劉延柱等,振動(dòng)力學(xué)[M],北京:高等教學(xué)出版社,1998.[13] 孟立凡,鄭賓傳感器原理及技術(shù)[M],北京:國(guó)防工業(yè)出版社,2005.[14] 付巍,基于Madah的單自由度振動(dòng)系統(tǒng)的數(shù)學(xué)仿真實(shí)驗(yàn)[J],機(jī)械工程與自動(dòng)化,2005(6):234.[15] 湯姆遜,達(dá)利,振動(dòng)力學(xué)[M],清華大學(xué)出版社,2005.[16] 哈爾濱工業(yè)大學(xué)理論力學(xué)教研室,理論力學(xué)（下）[M],北京:高等教學(xué)出版社,1997[17] 程偉,李敏,工程振動(dòng)基礎(chǔ)[:北京航空航天大學(xué)出版社,2004.[18] 蔣偉,機(jī)械動(dòng)力學(xué)分析[M],北京:中國(guó)傳媒大學(xué)出版社,2004.[19] (美)Singiresu S. Rao,著 李欣業(yè),張明路,編譯[M],北京:清華大學(xué)出版社,2009.[20] 鄭兆昌,機(jī)械振動(dòng)[M],北京:機(jī)械工業(yè)出版社,1980.[21] 于長(zhǎng)官,現(xiàn)代控制理論[M].哈爾濱:哈爾濱工業(yè)大學(xué)出版設(shè),2004.[22] 師漢民,機(jī)械振動(dòng)系統(tǒng)——[M],華中科技大學(xué)出版社,2004.[23] 程耀東,機(jī)械振動(dòng)學(xué)[J],杭州:浙江大學(xué)出版社,~66.附錄A譯文第三章 簡(jiǎn)諧激勵(lì)振動(dòng)間歇激勵(lì)是工程系統(tǒng)中經(jīng)常遇見的。我們將首先考慮受簡(jiǎn)諧力激勵(lì)的具有粘性阻尼的單自由度系統(tǒng)，： （） 這個(gè)方程的解分為兩部分，一部分是所謂余函數(shù)，即齊次方程的解，另一部分是特殊積分。把方程（）代入微分方程（），就可以找到上述方程中的振幅和相位角。于是得到了振幅和相位角的無(wú)因次表示式： （） （）這些公式表明，無(wú)因次的振幅及相位角僅僅是頻率比和阻率的函數(shù)，在靠近共振的頻率范圍內(nèi)，阻率對(duì)振幅和相位角有很大的影響。這時(shí)慣性力比較大，他與彈簧力相平衡。這里考慮限制在垂直方向移動(dòng)的彈簧質(zhì)量系統(tǒng)，收到旋轉(zhuǎn)部分失衡質(zhì)量的激勵(lì)，它的偏心距為并以角速度旋轉(zhuǎn)。這種力根據(jù)它們?cè)谵D(zhuǎn)子上的分布而產(chǎn)生靜失衡和動(dòng)失衡。 具有靜失衡的系統(tǒng)動(dòng)失衡 當(dāng)不平衡質(zhì)量不只在一個(gè)平面內(nèi)時(shí)，形成了一個(gè)力和一個(gè)偶，這屬于動(dòng)失衡。一般地說(shuō)，長(zhǎng)轉(zhuǎn)子，如電動(dòng)機(jī)轉(zhuǎn)子和汽車發(fā)動(dòng)機(jī)的曲線，可以看成是由一系列個(gè)圓盤組成的，每一個(gè)都可能有某種失衡。從本質(zhì)上看，平衡機(jī)是由支撐座構(gòu)成的，支承座用彈簧安裝，以便通過(guò)運(yùn)動(dòng)來(lái)檢測(cè)不平衡力。圓盤被支撐在受彈簧約束的軸承上。 進(jìn)行薄圓盤平衡試驗(yàn)接著將實(shí)驗(yàn)重量加到輪子的任一點(diǎn)，用同樣的速度重復(fù)上述過(guò)程。旋曲的定義為由彎曲軸線和支撐中心線組成的平面的轉(zhuǎn)動(dòng)。軸的旋曲是一個(gè)難以捉摸的課題。在這一節(jié)中我們將考慮同步旋曲這一最簡(jiǎn)單情況，在這中情況下，旋曲速度等于軸的轉(zhuǎn)速。用和代表軸心的的位置，質(zhì)量中心的坐標(biāo)是()和（）。我們令代表支撐點(diǎn)的諧位移，而質(zhì)量的坐標(biāo)是相對(duì)于固定坐標(biāo)的。（）和（）的圖解附錄B 中文題目：基于MATLAB的振動(dòng)系統(tǒng)編程分析外文題目：ANALYSIS OF VIBRATION SYSTEM BASED ON MATLAB ROGRAMMING畢業(yè)設(shè)計(jì)（論文）共 63 頁(yè)（其中：外文文獻(xiàn)及譯文27頁(yè)） 圖紙共0張 完成日期 2011年6月 答辯日期 2011年6月 。把這些公式代入方程（），得；或 （）因而振幅比的絕對(duì)值是： ()為了求相位角，我們使公式（）兩邊的實(shí)數(shù)部分和虛數(shù)部分各自相等便可確定和，然后取兩者的比值便得相位角的計(jì)算公式： （）確定穩(wěn)態(tài)振幅和相位的方程式（）和（）。當(dāng)旋轉(zhuǎn)速度等于臨界速度（或軸的橫向振動(dòng)的固有頻率）時(shí)，出現(xiàn)共振的條件，這時(shí)的振幅只是由于阻尼才受到限制。軸承中心線穿過(guò)盤平面的點(diǎn)，軸中性的變位為。軸旋曲的一般運(yùn)動(dòng)的研究超出本課程的范圍。軸的旋曲可以發(fā)生在軸的轉(zhuǎn)向的同向或反向。矢量差ab就是試驗(yàn)重量w1的數(shù)值增加到w0（oa/ab），矢量ab與矢量oa大小相等，方向相反，這時(shí)因?yàn)閄1為零，因此輪便被平衡。裝在軸承上的加速度計(jì)及閉光測(cè)頻儀可以做這些檢測(cè)。 有動(dòng)失衡的系統(tǒng) 轉(zhuǎn)子平衡機(jī)盡管薄圓盤可以用靜力法平衡，但同樣可以用動(dòng)平衡法。這種轉(zhuǎn)子必須用快速旋轉(zhuǎn)法（Spin）。例如，我們考慮一個(gè)有兩個(gè)圓盤的軸（），如果兩個(gè)不平衡質(zhì)量相等且相對(duì)于，則轉(zhuǎn)子將相對(duì)軸心處于穩(wěn)定平衡。這種失衡可以用靜力試驗(yàn)法來(lái)找；將輪軸安放在一對(duì)水平軌道上當(dāng)輪子比較重的點(diǎn)滾到輪軸最下方一點(diǎn)時(shí)輪子就停止了。當(dāng)轉(zhuǎn)速增加到大大超過(guò)共振頻率時(shí)，求系統(tǒng)的阻率。共振時(shí)的振幅可以根據(jù)公式（）或公式（），： （）當(dāng)時(shí)，接近，外力幾乎全花費(fèi)在克服慣性力上。對(duì)于的情況，慣性力和阻尼力都小，這使得相位角也小，外加力的值接近等于彈簧力。從圖很容易看出： （） （）現(xiàn)在我們把公式（）和（）寫成無(wú)因次的形式，這樣就能夠把結(jié)果用簡(jiǎn)明的圖像表示。上述方程的特解是頻率與激勵(lì)頻率相同的穩(wěn)態(tài)振動(dòng)。盡管純簡(jiǎn)單間歇激勵(lì)似乎比其他周期激勵(lì)要較少遇見，但為了更好的理解在更一般的激勵(lì)下系統(tǒng)的相應(yīng)，明白在間歇激勵(lì)下系統(tǒng)的特性是必要的。感謝幫助過(guò)我的同學(xué)們，是你們讓我的大學(xué)生活豐富多彩。對(duì)理論力學(xué)的所有老師和工程力學(xué)系的所有老師在理論研究過(guò)程中提供的幫助表示感謝。導(dǎo)師淵博的知識(shí)、嚴(yán)謹(jǐn)?shù)闹螌W(xué)作風(fēng)、精益求精的鉆研精神和忘我工作的精神都值得筆者不斷的學(xué)習(xí)。5結(jié)論本文以振動(dòng)力學(xué)為基礎(chǔ)對(duì)于抽象的振動(dòng)系統(tǒng)進(jìn)行了分析和數(shù)學(xué)模擬，再引用了MATLAB編程軟件對(duì)單自由度的自由振動(dòng)數(shù)學(xué)模型進(jìn)行圖像數(shù)據(jù)模擬，使振動(dòng)數(shù)據(jù)更加明顯，對(duì)粘性阻尼自由振動(dòng)的圖像實(shí)現(xiàn)充分體現(xiàn)了MATLAB的優(yōu)點(diǎn)是真的數(shù)據(jù)更加明顯清晰，對(duì)于單自由的系統(tǒng)的MATLAB實(shí)現(xiàn)了振幅放大因子的三個(gè)有關(guān)圖像，而多自由度的振動(dòng)系統(tǒng)中MATLAB的應(yīng)用時(shí)主振型的計(jì)算更加簡(jiǎn)單速度，連續(xù)體的振動(dòng)同時(shí)也應(yīng)用了MATLAB的公式圖像實(shí)現(xiàn)使公式更加明顯。)。w(x,t)39。xlabel(39。subplot(39。)。)。x=20,n=139。ylabel(39。endsubplot(39。for i=1:1001 t(i)=3*(i1)/1000。wn=(n^2)*。l=40。運(yùn)用杜哈美積分，式（431）的解可以表示為 （434）算例作圖表示 （435）%x=20。常見梁的邊界條件 1)自由端：彎矩和剪力分別為零，即 （424） 2)簡(jiǎn)支端：撓度和彎矩分別為零，即 （425） 3)固定端：撓度和轉(zhuǎn)角分別為零，即 （425）4)梁