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

HarbinInstituteofTechnology

材料力学上机报告

课程名称：材料力学

设计题目：二向应力状态分析

院系：XXXXXX

班级：XXXXXX

设计者：XXXXXX

学号：XXXXXX

设计时间：2013.06.18

哈尔滨工业大学

二向应力状态分析

一：课题要求

1.输入：任意一点的应力状态：（σ、σ、τ）；某截面方位角α

xyxy

2.输出：输入点的主应力（σ、σ、σ），方位角α斜截面上的应力σ、τ

123αα

及主方向角α。

0

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;

ifx10

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);gridon%绘制应力圆

以上程序为在matlab中使用GUI编程时的主代码，界面代码请见m文件。

四：运行过程、结果和应力圆

在matlab中打开m文件，按F5使程序运行，显示窗口如下：

左侧为输入窗口，中间为相应的主切应力和斜截面应力的输出窗口，右侧为二向

应力圆的输出窗口。现在在左侧输入一些值，显示如下：

点击Enter键，即可得到相应的值与应力圆，如下图所示：