理論力學(xué)仿真實(shí)驗(yàn)報(bào)告

理論力學(xué)仿真實(shí)驗(yàn)報(bào)告姓名班級(jí)學(xué)號(hào)成績(jī)理論力學(xué)仿真實(shí)驗(yàn)實(shí)驗(yàn)?zāi)康恼莆樟W(xué)基本原理，培養(yǎng)力學(xué)自學(xué)能力、力學(xué)建模能力（建立力學(xué)模型的能力和建立數(shù)學(xué)模型的能力）和力學(xué)分析能力（力學(xué)理論分析能力、力學(xué)實(shí)驗(yàn)分析能力、力學(xué)數(shù)值仿真能力）。實(shí)驗(yàn)內(nèi)容力系平衡問題、復(fù)擺運(yùn)動(dòng)分析研究實(shí)驗(yàn)原理在做力系平衡實(shí)驗(yàn)時(shí)根據(jù)梁受力平衡列出方程求解未知力，在復(fù)擺運(yùn)動(dòng)分析研究實(shí)驗(yàn)時(shí)根據(jù)微分方程列出微分方程組，將力學(xué)模型轉(zhuǎn)化為數(shù)學(xué)模型，代入MATLAB程序，會(huì)得出復(fù)擺的運(yùn)動(dòng)規(guī)律曲線。實(shí)驗(yàn)過程力系平衡問題⑴力學(xué)模型建立、描述⑵數(shù)學(xué)模型建立對(duì)左端列出平衡方程：對(duì)右端列出平衡方程：⑶數(shù)學(xué)模型求解仿真梁的左端列平衡方程：eq1='FAx-FD*cos(1/6*pi)=0';eq2='FAy+FBy+FD*sin(1/6*pi)-4*q=0';eq3='FBy*2-8*q-M+8*FD*sin(1/6*pi)=0';梁的右端列平衡方程：eq4='4*FD*sin(1/6*pi)-M=0';eq5='FCx-FD*cos(1/6*pi)=0';eq6='FCy+FD*sin(1/6*pi)=0';s=solve(eq1,eq2,eq3,eq4,eq5,eq6,...'FAx','FAy','FBy','FCx','FCy','FD');q=5000;M=20000;代入方程，F(xiàn)Ax=subs(s.FAx)FAy=subs(s.FAy)FBy=subs(s.FBy)FCx=subs(s.FCx)FCy=subs(s.FCy)FD=subs(s.FD)2、復(fù)擺運(yùn)動(dòng)分析研究⑴力學(xué)模型建立、描述⑵數(shù)學(xué)模型建立設(shè)m=5kgg=10N/kga=1.5mJ=15=ω=0列微分方程組：⑶數(shù)學(xué)模型求解仿真[]=[，]列微分方程組：>>[T,Y]=ode45('odefun',[01],[200]);>>plot(T,Y(:,1),'-',T,Y(:,2),'-.')實(shí)驗(yàn)結(jié)果分析討論1、力系平衡問題的MATLAN的求解結(jié)果為：FAx=8.6603e+003FAy=5000FBy=10000FCx=8.6603e+003FCy=-5000FD=10000故：2、復(fù)擺運(yùn)動(dòng)分析研究的MATLAB的求解結(jié)果為：故：>>[T,Y]=ode45('odefun',[