哈工大材力上机报告应力状态分析.pdf

材料力学上机报告 ——应力状态分析 班号： 学号： 姓名： 哈尔滨工业大学 MATLAB 源代码 clear all ; flag=input( 计算双向应力请输入 2 计算三向应力请输入 3\n ); if flag==2 sx=input( 请输入 x 向正应力 = ); sy=input( 请输入 y 向正应力 = ); txy=input( 请输入 xy 平面切应力 = ); a=linspace(0,pi,37); sa=(sx+sy)/2; sd=(sx-sy)/2; sigma=sa+sd*cos(2*a)+txy*sin(2*a); tau=-sd*sin(2*a)+txy*cos(2*a); plot(sigma,tau, b- ),grid axis equal ; v=axis; line([v(1),v(2)],[0,0]); line([0,0],[v(3),v(4)]) hold,plot(sa,0, x ) smax=max(sigma);smin=min(sigma);tmax=max(tau); asigma=((atan((2*txy)/(sx-sy)))/2)/pi*180; fprintf( σmax=%8.5fMpa\n ,smax); fprintf( σmin=%8.5fMpa\n ,smin); fprintf( τmax=%8.5fMpa\n ,tmax); fprintf( α 0=%8.5f 度 和 %8.5f 度 \n ,asigma,asigma+90); a=input( 请输入方位角 α/ °= ); a=a*pi/180; sigma=sa+sd*cos(2*a)+txy*sin(2*a); tau=sd*sin(2*a)-txy*cos(2*a); plot(sigma,tau, or ) fprintf( σ=%8.5fMpa\n ,sigma); fprintf( τ=%8.5fMpa\n ,tau); end if flag==3 del1=5; sx=input( 请输入 x 向正应力 : ); sy=input( 请输入 y 向正应力 : ); sz=input( 请输入 z 向正应力 : ); txy=input( 请输入 xy 平面切应力： ); tyz=input( 请输入 yz 平面切应力： ); tzx=input( 请输入 zx 平面切应力： ); I1=sx+sy+sz; tmp1=[sx,txy txy,sy]; tmp2=[sy,tyz tyz,sz]; tmp3=[sz,tzx tzx,sx]; I2=det(tmp1)+det(tmp2)+det(tmp3); tmp4=[sx,txy,tzx