基于MATLAB的壓床機構(gòu)運動學(xué)與動力學(xué)分析

潘宇 陸鳳儀 岳一領(lǐng) 肖欣一

摘 要：文章建立了壓床機構(gòu)運動分析和動力學(xué)分析模型，并運用MATLAB對其進(jìn)行分析與仿真。該方法可為其他連桿分析提供借鑒。

關(guān)鍵詞：壓床機構(gòu)；MATLAB；運動學(xué)；動力學(xué)

平面機構(gòu)運動學(xué)與動力學(xué)分析方法有圖解法和解析法。圖解法雖然具有形象直觀的特點，但是從現(xiàn)代科技和工業(yè)發(fā)展的要求來看，它不僅精度較低，費時較多，而且不便于把機構(gòu)分析的問題和機構(gòu)綜合問題聯(lián)系起來。解析法正好能克服上述缺點，隨著計算機技術(shù)的發(fā)展和普及，其應(yīng)用將是越來越廣泛[1]。在機械原理課程設(shè)計中，因壓床機構(gòu)具有典型的代表性，故以此作為設(shè)計題目。本文以壓床機構(gòu)為研究對象，建立其運動學(xué)和動力學(xué)的分析模型，并利用工程中常用且適宜平面機構(gòu)分析和仿真的Matlab語言，并編制了集分析和仿真于一體的分析軟件。

1 壓床機構(gòu)運動分析模型

1.1 位置分析

壓床機構(gòu)機構(gòu)運動簡圖見圖1。

式中 Li-各構(gòu)件的長度，i=1，2，3，4。

？茲i-各構(gòu)件矢量的方位角，i=1，2，3，4。

1.2 速度分析

1.3 加速度分析

2 壓床機構(gòu)動力分析模型

壓床機構(gòu)各運動副反力寫為FRij，且FRij=FRijx+FRijy。同一運動副中，作用在不同構(gòu)件上的兩個力，大小相等，方向相反，即FRij=-FRji。各運動副處的力對質(zhì)心的力矩計算方法與理論力學(xué)中的相同，方向逆時針為正，順時針為負(fù)。限于篇幅慣性力和慣性力矩的計算略。僅給出運動副反力和平衡力矩矩陣形式的平衡方程：

AFR=D

式中：A-為系數(shù)矩陣，F(xiàn)R-為未知力列陣，D-為已知力列陣。其中：

3 基于MATLAB 分析及仿真的實現(xiàn)

基于上述分析，運用MATLAB 語言編寫實現(xiàn)壓床運動學(xué)、動力學(xué)分析及仿真程序。其主要使用的功能函數(shù)和命令有：程序運行前，clc命令用來清除命令窗口的內(nèi)容，clear命令可以用來清除工作空間的內(nèi)容，防止誤操作和數(shù)據(jù)混亂。運用控制流for-end循環(huán)結(jié)構(gòu)，按照迭代的方法，以曲柄作為循環(huán)控制變量實現(xiàn)對整個運動周期，按照運動規(guī)律進(jìn)行逐點的參數(shù)計算。為了循環(huán)結(jié)構(gòu)的運行速度更快，可以再循環(huán)結(jié)構(gòu)執(zhí)行之前采用全零數(shù)組命令zeros（n）預(yù)先定義數(shù)組的大小。

利用本文編制的程序，按表1和2中給出的參數(shù)進(jìn)行壓床機構(gòu)的運動學(xué)、動力學(xué)運算與仿真。圖2為運動學(xué)分析仿真曲線。圖3為動力學(xué)分析仿真曲線。

4 結(jié)束語

本文基于MATLAB，完成了壓床機構(gòu)運動學(xué)分析和動力學(xué)分析，并實現(xiàn)了運動仿真。運算程序簡潔明了、易于實現(xiàn)?？蔀槠渌愋推矫孢B桿機構(gòu)運動學(xué)和動力學(xué)提供較好的借鑒。

參考文獻(xiàn)

[1]陸鳳儀.機械原理課程設(shè)計[M].（2版）北京：機械工業(yè)出版社，2011.

[2]孫桓，陳作模，葛文杰.機械原理[M].（7版）北京：高等教育出版社，2006.

[3]郭仁生.機械工程設(shè)計分析和MATLAB應(yīng)用[M].北京：機械工業(yè)出版社，2005.

[4]李濱城，徐超.機械原理MATLAB輔助分析 [M].北京：化學(xué)工業(yè)出版社，2011.

作者簡介：潘宇（1989，5-），機自101202班本科生，研究方向：機械制造及自動化。

陸鳳儀（1958-），女，江蘇常州人，太原科技大學(xué)教授，碩導(dǎo)，省級教學(xué)名師，機械設(shè)計類課程建設(shè)負(fù)責(zé)人、機械設(shè)計基礎(chǔ)實驗教學(xué)中心主任。