matlab課程設(shè)計(jì)word版(編輯修改稿)
2025-07-13 14:22 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 39。RungKutta法解 39。,2) hold on plot(x1,y1,39。bo39。,xr,yr,39。r*39。) 程序運(yùn)行結(jié)果： MyRunge_Kutta_x = 0 MyRunge_Kutta_y = 0 二、變成解決以下科學(xué)計(jì)算問(wèn)題 ： （一） [例 724] 材料力學(xué)復(fù)雜應(yīng)力狀態(tài)的分析 —— Moore 圓。 程序 說(shuō)明： 利用平面應(yīng)力狀態(tài)下斜截面應(yīng)力的一般公式，畫出任意平面應(yīng)力狀態(tài)下的應(yīng)力圓（ Moore 圓），求出相應(yīng)平面應(yīng)力狀態(tài)下的主應(yīng)力 （ 1? 、 3? ） ，并求出該應(yīng)力狀態(tài)下任意方位角 ? 的斜截面上 的 應(yīng)力 ? 、 ? 。 ?x?y?y?x???yx?yx?xy?xy? 程序流程圖： 程序代碼： clear。clc。 Sx=input(39。Sigma_x(MPa)=39。)。 %輸入該應(yīng)力狀態(tài)下的各應(yīng)力值 Sy=input(39。Sigma_y(MPa)=39。)。 Txy=input(39。Tau_xy(MPa)=39。)。 a=linspace(0,pi,100)。 %等分半圓周角 Sa=(Sx+Sy)/2。 Sd=(SxSy)/2。 Sigma=Sa+Sd*cos(2*a)Txy*sin(2*a)。 %應(yīng)力圓一般方程 Tau=Sd*sin(2*a)+Txy*cos(2*a)。 plot(Sigma,Tau,Sx,Txy,39。.r39。,Sy,Txy,39。.r39。)。 %畫出應(yīng)力圓，標(biāo)出該應(yīng)力狀態(tài)下各應(yīng)力參數(shù) line([Sx,Sy],[Txy,Txy])。 axis equal。 %控制各坐標(biāo)軸的分度使其相等 —— 使應(yīng)力圓變圓 開始 輸入待求應(yīng)力狀態(tài)的參數(shù) 畫出應(yīng)力圓 求某一方向角截面上的應(yīng)力 ? 輸入方向角 求出相應(yīng)正應(yīng)力、切應(yīng)力 是 否 得出該應(yīng)力狀態(tài)下的主應(yīng)力 求出主應(yīng)力平面方向角 結(jié)束 title(39。應(yīng)力圓 39。)。 xlabel(39。正應(yīng)力（ MPa） 39。)。 ylabel(39。剪應(yīng)力（ MPa） 39。)。 text(Sx,Txy,39。A39。)。 text(Sy,Txy,39。B39。)。 Smax=max(Sigma)。 Smin=min(Sigma)。 Tmax=max(Tau)。 Tmin=min(Tau)。 b=axis。
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號(hào)-1