资源描述：
《材料力学上机作业》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、络嚼廈力掘fc荐HarbinInstituteofTechnology材料力学上机作业课程名称:设计题目:材料力学应力状态分析院系:班级:姓名:学号:指导教师:设计时间:哈尔滨工业人学1.本程序用Matlab编程，程序如下:clearall;g=input(「如果求二向应力，输入2,如果求三向应力，输入3*);whileg==3;clearallsx=input(1sx(Mpa)=1);sy=input(1sy(Mpa)=1);sz=input(1sz(Mpa)=1);txy=input(1txy(Mpa)=!);tyz=input(!tyz(Mpa)=!);txz=input(
2、1txz(Mpa)=1);il=sx+sy+sz;i2=sx*sy-txy*txy+sy*sz-tyz*tyz+sx*sz-txz;i31=[sxztxy,txz;txy,sy,tyz;txzztyzAsz];i3=det(i31);a=[l,-ilzi2,-i3]；s=roots(a);sl=max(s);s3=min(s);fork=l:1:3if((s(k)~=sl)&&(s(k)~=s3))s2=s(k);endendg=i;tmax=(sl-s3)/2;rl3=(sl-s3)/2;rl2=(sl-s2)/2;r23=(s2-s3)/2;alphal=0:pi/100:2
3、*pi;[0,2*pi]R=rl3;%半径xl=R*cos(alphal)+s3+rl3;yl=R*sin(alphal);plot(xlzylf1*1)zgrid,holdonaxisequalv=axis;line([v(1),v(2)],[0,0]);line([0z0]z[v(3)zv(4)]);alpha2=0:pi/100:2*pi;2*pi]R=rl2;%半径x2=R*cos(alpha2)+s2+rl2;y2=R*sin(alpha2);plot(x2,y2r1*1)Agrid,holdonalpha3=0:pi/100:2*pi;度[0,2*pi]R=r23;%
4、半径x3=R*cos(alpha3)+s3+r23;y3=R*sin(alpha3);plot(x3,y3,'k-grid,holdonfprintf('主应力1=%8.5fMpa',si);fprintf('^JAc^/2=%8.5fMpa'zs2);fprintf('上应力3=%8.5fMpa',s3);fprintf{'最人剪应力=%8.5fMpa*,tmax);holdoffendwhileg==2;sx=input(1sx(Mpa)=1);sy=input(1sy(Mpa)=1);txy=input('txy(Mpa)=!);a=linspace(0,p
5、i,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,sx,txy,1b-1);axisequal;v=axis;line([v(l),v(2)],[0,0]);line([0,0],[v(3),v(4)])hold,plot(sa,0,1x1)smax=max(sigma),smin=min(sigma),tmax=max(tau);asigma=((atan((2*txy)/(sx-sy)))/2)/pi*180
6、;fprintf(1力最人值=%8.5fMpa1,smax);fprintf(，主应力最小值=%8.5fMpa1,smin);fprintf(1切应力最大值=%8.5fMpa1,tmax);fprintf(T主方向角=%8・5f度和％8・5f度*rasigma,asigma+90);h=input「如果不求应力，输入0,如果要再求应力，输入1*);whileh==l;a=input('纟合出斜截面方向角a=(弧度)*)sigma=sa+sd*cos(2*3)-txy*sin(2*a)tau=sd*sin(2*3)+txy*cos(2*a)plot(sigm“tau
7、,1or1)h=input(*如果不继续求丿'V力，输入0,如果还要求应力，输入1end,holdoffg=o；end2•程序在Matlab+运行时的界面：（以三向应力为例）CommandWindow①NewtoMATLAB?WatchthisVideo,seeDemos,orreadGettinaStarted.X»ymgli2fx如果求二向应力，输入2,如果求三向应力，输入3
8、应力圆示意图:XdNewtoMATLAB?WatchthisVideo,seeDemos,orrea
