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

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键，即可得到相应的值与应力圆，如下图所示：