12 / 2 页

楼主: liaoma010101

三维应力状态下主应力平面如何确定？

[复制链接]

字体大小: 正常放大

hundunwei

3 主题	4 听众	89 积分

升级 88.42%

TA的每日心情

	郁闷 2012-6-9 15:23

签到天数: 16 天

[LV.4]偶尔看看III

群组: Matlab讨论组

群组: 2012第三期美赛培训

11^#

发表于 2012-2-8 19:23 |只看该作者

|招呼Ta 关注Ta

三维的类似，画三向应力圆

我也说一句

发表

使用道具举报

hundunwei

3 主题	4 听众	89 积分

升级 88.42%

TA的每日心情

	郁闷 2012-6-9 15:23

签到天数: 16 天

[LV.4]偶尔看看III

群组: Matlab讨论组

群组: 2012第三期美赛培训

12^#

发表于 2012-2-8 19:32 |只看该作者 |招呼Ta 关注Ta

这有一个三向应力圆例子
ox=30,oy=60,oz=20,Txy=Tyz=Txz=0
C:\Users\Administrator\Desktop\捕获
再想了解你就看书吧
张少实《新编材料力学》等

我也说一句

发表

使用道具举报

hundunwei

3 主题	4 听众	89 积分

升级 88.42%

TA的每日心情

	郁闷 2012-6-9 15:23

签到天数: 16 天

[LV.4]偶尔看看III

群组: Matlab讨论组

群组: 2012第三期美赛培训

13^#

发表于 2012-2-9 15:05 |只看该作者 |招呼Ta 关注Ta

诶？图片显示不了？

我也说一句

发表

使用道具举报

hundunwei

3 主题	4 听众	89 积分

升级 88.42%

TA的每日心情

	郁闷 2012-6-9 15:23

签到天数: 16 天

[LV.4]偶尔看看III

群组: Matlab讨论组

群组: 2012第三期美赛培训

14^#

发表于 2012-2-9 15:08 |只看该作者 |招呼Ta 关注Ta

给个程序吧，matlab的，看程序，想方法
%三向应力状态
ox=input('x向应力');
oy=input('y向应力');
oz=input('z向应力');
Txy=input('xy向切应力');
Tyz=input('yz向切应力');
Tzx=input('zx向切应力');
A=[ox -Txy;Txy oy];
B=[oy -Tyz;Tyz oz];
C=[oz -Tzx;Tzx ox];
D=[ox -Txy Tzx;Txy oy -Tyz;-Tzx Tyz oz];
I1=ox+oy+oz;
I2=det(A)+det(B)+det(C);
I3=det(D);
E=roots([1,-I1,I2,-I3]);
F=sort(E);
a1=F(3)
a2=F(2)
a3=F(1)
Tmax=(a1-a3)/2
w=0:0.01*pi:2*pi;
x1=(a1+a2)/2+(a2-a1)/2*cos(w);
y1=(a2-a1)/2*sin(w);
x2=(a1+a3)/2+(a3-a1)/2*cos(w);
y2=(a3-a1)/2*sin(w);
x3=(a3+a2)/2+(a3-a2)/2*cos(w);
y3=(a3-a2)/2*sin(w);
v=-30:0.01:50;
plot(x1,y1,x2,y2,x3,y3,v,0,0,v);
title('应力圆')