材料力学上机大作业(哈工大)

H a r b i n I n s t i t u t e o f T e c h n o l o g y
材料力学上机报告
课程名称：材料力学
设计题目：二向应力状态分析
院系：XXXXXX
班级：XXXXXX
设计者：XXXXXX
学号：XXXXXX
设计时间：2013.06.18
哈尔滨工业大学
二向应力状态分析
一：课题要求
1.输入：任意一点的应力状态：（σx、σy、τxy）；某截面方位角α
2.输出：输入点的主应力（σ1、σ2、σ3），方位角α斜截面上的应力σ
α、τα。

及主方向角α
3.画出应力圆示意图。

4.程序运行时为界面显示形式。

二：程序框图
三：所编程序
x=str2double(get(handles.edit1,'string'));
y=str2double(get(handles.edit2,'string'));
xy=str2double(get(handles.edit3,'string'));
M=str2double(get(handles.edit4,'string')); %将窗口输入值分别赋给x,y,xy,M b=sqrt((x/2-y/2)^2+xy^2);x1=(x+y)/2+b;x3=(x+y)/2-b;
x2=0;
if x1<0
x2=x1;
x1=0;
end
t=(x1-x3)/2;
M=M*pi/180;
b1=(x+y)/2+(x-y)*cos(2*M)/2-xy*sin(2*M);
b2=(x-y)*sin(2*M)/2+xy*cos(2*M);
b3=90*atan((-2*xy)/(x+y))/pi;%计算输出的主切应力大小、方向和截面上的应力并赋值set(handles.edit5,'string',x1);
set(handles.edit6,'string',x2);
set(handles.edit7,'string',x3);
set(handles.edit9,'string',t);
set(handles.edit10,'string',b3);
set(handles.edit11,'string',b1);
set(handles.edit12,'string',b2);%在输出窗口显示主切应力大小、方向和截面上应力
b4=sqrt(b.^2+t.^2);
v1=(x+y)/2-b4:0.001:(x+y)/2+b4;
b11=sqrt(b4.^2-(v1-(x+y)/2).^2);b12=-sqrt(b4.^2-(v1-(x+y)/2).^2);
%绘制应力圆上的点
axes(handles.axes1); %选择应力圆的输出地址
plot(v1,b11,v1,b12);grid on%绘制应力圆
以上程序为在matlab中使用GUI编程时的主代码，界面代码请见m文件。

四：运行过程、结果和应力圆
在matlab中打开m文件，按F5使程序运行，显示窗口如下：
左侧为输入窗口，中间为相应的主切应力和斜截面应力的输出窗口，右侧为二向
应力圆的输出窗口。

现在在左侧输入一些值，显示如下：
点击Enter键，即可得到相应的值与应力圆，如下图所示：。