哈工大材力上机大作业讲义.doc

材料力学上机大作业 题目名称：应力状态分析 班号： 学号： 姓名： 学院： 机电工程学院 指导教师： 完成时间： 2015.6.20 使用软件： Matlab2014B 应力状态分析 解决问题方向 该程序可以解决的问题：① 平面应力状态，已知σx，σy，τxy和角度的情况下计算不同任意截面的应力分量，还可以计算平面应力状态主应力的大小和方向并能画出应力圆；② 空间应力状态，已知σx，σy，σz，τxy等的情况下计算主应力和最大切应力。 数学模型 平面应力状态任一斜截面上正应力切应力公式为： 平面应力状态主应力大小及方向主应力计算公式 主平面方位计算公式 3、最大切应力：45度角的平面方位 5、三向应力状态 特征方程： 3个不变量： 三、程序流程图 1.平面应力计算 （1）打开M程序之后，选择平面应力计算部分。提示输入数据σx，σy，τxy 。 （2）点击回车后计算。 （3）程序算出任意截面的应力分量以及三个主应力和最大切应力。 （4）启动绘图，可以得到相应的平面应力圆图像。 空间应力计算 （1）打开M程序之后，选择空间应力计算部分。提示输入数据σx，σy，σz，τxy ，τyz ，τxz。 （2）点击回车后计算。 （3）程序算出任意截面的应力分量以及三个主应力和最大切应力。 （4）启动绘图，可以得到相应的空间应力圆图像。 四、程序源代码 %材力大作业：应力状态分析 rum=input(选择应力状态方式(1或2)：);%1表示平面应力状态，2表示空间应力状态 if rum==1 %选择平面应力状态分析 %输入已知量，应力单位为MPa，转角单位为rad dnx=input(输入x轴方向正应力 cmx=); dny=input(输入y轴方向正应力 cmy=); txy=input(输入切应力 txy=); %求解主应力、主方向及最大剪应力并输出 dn1=(dnx+dny)/2+(((dnx-dny)/2)^2+txy^2)^(1/2); dn2=0; dn3=(dnx+dny)/2-(((dnx-dny)/2)^2+txy^2)^(1/2); tm=(dn1-dn3)/2; afc=(1/2)*atan(2*txy/(dnx-dny)); dnt=(dnx+dny)/2+(dnx-dny)*cos(2*afc)/2+txy*sin(2*afc); if dnt==dn1; verfc1=afc; verfc3=afc+pi/2; else verfc3=afc; verfc1=afc+pi/2; end display(主应力为：); display(dn1); display(dn2); display(dn3); display(主方向为：); display(verfc1); display(verfc3); display(最大切应力为：); display(tm); % 求解任意截面上的应力 afa=input(输入截面方位(以弧度表示) aerfa=); dnr=(dnx+dny)/2+(dnx-dny)*cos(2*afa)/2+txy*sin(2*afa); tr=-(dnx-dny)*sin(2*afa)/2+txy*cos(2*afa); display(截面处应力状况：); display(正应力：); display(dnr); display(切应力：); display(tr); %求解应力圆并作图 i=0; for theta=0:pi/200:pi dnt=(dnx+dny)/2+(dnx-dny)*cos(2*theta)/2+txy*sin(2*theta); t=-(dnx-dny)*sin(2*theta)/2+txy*cos(2*theta); i=i+1; C(i)=dnt;T(i)=t; plot(C,T),axis equal; title(应力圆);xlabel(正应力cgm/ MPa);ylabel(切应力t/MPa);grid on; end elsei