male實驗報告調研報告WORD版

1、本文格式為Word版，下載可任意編輯male實驗報告調研報告WORD版 e male 試驗報告調研報告精編 D WORD 版 機電系統(tǒng)設計和仿真 姓名 ： 陳俊英 學號 ： SY1317315 指導老師 ：李慶波 儀器科學與光電工程學院 4 2021 年 年 6 6 月 月 2 22 日 Maple- - 滑塊擺 試驗 . 1. 試驗要求 滑塊擺由一置于光滑桿上的質量為 m 的滑塊 A、一質量為 M 的小球 B 和長度為 L ，質量不計的剛性桿鉸接而成，不計各處摩擦，以過 A 點的水平面為零勢能面，通過Lagrange 方程建立系統(tǒng)的運動方程，利用 Maple 軟件畫出： 1） 滑塊 A 的位

2、移 x 隨時間 t 的變化曲線 2） 角度 隨時間 t 的變化曲線 3） 滑塊擺的運動動畫 圖 1 滑塊擺示意圖 . 2. 試驗說明 設定初始條件為：m=1Kg，M=1Kg，g=9.8，L=2m (0) = 0rad， x(0) = 0m， (0) = -1.3rad/s,，x(0) = 1m/s 如下定義的拉格朗日方程： 其中:q：x(t)和 (t)的自由度 D：由于摩擦而消耗的能量 F q ：由自由度 q 產生的力 E c 和 E p ：系統(tǒng)的動能和勢能 系統(tǒng)有兩個自由度，以 x 和 j 為廣義坐標，以過 A 點的水平面為零勢能面，系統(tǒng)的動能和勢能分別為 系統(tǒng)的 Lagrange 方程為：

3、 計算出諸導數(shù)： 帶入 Lagrange 方程，得到系統(tǒng)的運動微分方程： . 3. 試驗步驟 步驟 1：安裝 Maple 軟件； 步驟 2：依據(jù)單擺試驗的試驗說明和程序例程學習 maple 的編程語言； 步驟 3：依據(jù)滑塊擺的試驗說明，根據(jù)試驗要求編寫 maple 仿真程序； 步驟 4：得到試驗結果。 4. 試驗結果 1）滑塊 A 的位移 x 隨時間 t 的變化曲線 圖 2 滑塊的位移隨時間的變化曲線 2）角度 隨時間 t 的變化曲線 圖 3 角度隨時間的變化曲線 3）滑塊擺的運動動畫 圖 4 滑塊擺運動動畫 附錄程序： restart; with (DEtools): # # 調用微分方程工

4、具包 with (plots): # # 調用圖形包 with (plottools): # # 調用圖形工具包 g:=9.8;L:=2;M:=1;m:=1; equ1:=diff(x(t),t$2)*(m+M)+M*L*diff(theta(t),t$2)*cos(theta(t)- -M*L*diff(theta(t),t)2*sin(theta(t)=0; equ2 2 := L*diff(theta(t),t$2)+diff(x(t),t$2)*cos(theta(t)+g*sin(theta(t)=0; val:=array(1.100); for i to 100 do vali:

5、=i/10 end do; S:=dsolve(equ1,equ2,theta(0)=0,D(theta)(0)=- -1.3,x(0)=0,D(x)(0)=1,x(t),theta( (t t) ) ,method=rkf45,type=numeric,output=val); eval(s); # # 繪制位移隨時間變化曲線圖 DEplot(equ1,equ2,x(t),theta(t),t=0.10,theta(0)=0,D(theta)(0)=- -1.3,x(0)=0,D(x)(0)=1,stepsize=0.1,linecol or=t/2); # # 繪制相位隨時間變化曲線圖 D

6、Eplot(equ1,equ2,theta(t),x(t),t=0.10,theta(0)=0,D(theta)(0)=- -1.3,x(0)=0,D(x)(0)=1,stepsize=0.1,linecolor=t/2); # # 繪制滑塊擺運動動畫 for t to 100 do x1(t):=S2, 1t, 4:x(t):=x1(t)+L*sin(S2, 1t, 2):y(t):=- - L*cos(S2, 1t, 2): end do: for t to 100 do: wall:=curve(- - 4.0,0,0.5,0,color=black, linestyle=solid, thickness=2): fixp:=rectangle(x1(t)- - 0.1,0.1, x1(t)+0.1,- - 0.1,color=turquoise): pendulum:=disk(x(t),y(t),0.10,color=turquoise): lineW_P:=curve(x1(t),0,x(t),y(t ),color=red): G1t:=display(pendulum,lineW_P,wall,fixp): end do: d