哈工大材力上机 matlab,材力上机Matlab哈工大

《材力上机Matlab哈工大》由会员分享，可在线阅读，更多相关《材力上机Matlab哈工大(7页珍藏版)》请在金锄头文库上搜索。

1、一、解决题目由于之前只学过简单的 C语言， 对于本次程序设计的帮助不大， 所以重新学习了Matlab 语言，并对题目 7 进行处理题目 7 应力状态分析输入：1. 平面应力状态输入：( x、 y、 xy)；某截面方位角 2. 空间应力状态输入：( x、 y、 z、 xy、 yz、 zx)输出：1. 输出主应力( 1、 2、 3)2. 最大剪应力( max= 13=23-1 错误！未找到引用源。 )3. 如为平面应力状态则需输出方位角 斜截面上的应力 、 及主方向角 错误！未找到引用源。4. 画出应力圆示意图二、方向规定正面上应力以沿着坐标轴正方向者为正， 反之为负： 负面上应力以沿着坐标轴负方。

2、向者为正，反之为负。三、程序设计思路本程序分为两个部分， 分别对应平面应力分析与空间应力分析， 可由使用者进行选择。平面应力分析主要是主应力和任一截面位置上的应力状态，还有应力圆示意图。空间应力状态主要是主应力问题四、程序编写程序采用 Matlab 语言，所以对于一些特殊字符，无法识别。首先对变量进行定义：n 为用户选择进入的两个程序部分， n=1 时，进行平面应力状态分析， n=2 时，进行空间应力状态分析。x， y， z 分别对应三个方向的正应力。xy， yz， zx 分别对应切应力( yx， zy， xz 省略) 。j 为用户输入的某截面 角其他一些代码为程序内部执行代码，不再一一介绍。。

3、下面是主程序代码n=input(If plane,enter 1,If spatial,enter 2 n=); 用户选择并进入相应界面if n=1 平面应力分析x=input(x=); y=input(y=); xy=input(xy=); j=input(j=); 用户输入数据jj=j*pi/180; aa=(x+y)/2; bb=sqrt(x-y)/2)2+xy2); cc=(x-y)/2; a=aa+bb; b=aa-bb; if a=b a1=a a2=0 a3=b tmax=(a1-a3)/2 az=aa+cc*cos(2*jj)+xy*sin(2*jj) tz=-cc*sin(2。

4、*jj)+xy*cos(2*jj) alpha=0:pi/50:2*pi; X=sqrt(b2+xy2)*cos(alpha)+a; Y=sqrt(b2+xy2)*sin(alpha); plot(X,Y,-) else a1=b a2=0 a3=a tmax=(a1-a3)/2 az=aa+cc*cos(2*jj)+xy*sin(2*jj) tz=-cc*sin(2*jj)+xy*cos(2*jj) alpha=0:pi/50:2*pi; X=sqrt(b2+xy2)*cos(alpha)+a; Y=sqrt(b2+xy2)*sin(alpha); plot(X,Y,-) end 计算，绘图。

5、并结束程序else x=input(x=); y=input(y=); z=input(z=); xy=input(xy=); yz=input(yz=); zx=input(zx=); I1=x+y+z; I2=-xy*xy+x*y+y*z-yz*y+z*x-zx*zx; I3=-zx*y*zx-x*yz*yz-z*xy*xy+x*y*z+xy*yz*zx+zx*xy*yz; syms Z; solve(Z3-I1*Z2+I2*Z-I3) b1=eval(ans(3); b2=eval(ans(2); b3=eval(ans(1); b0=0; if b10 b0=b1; b1=b2; b2。

6、=b0; if b10 b0=b1; b1=b3; b3=b0; b0=b2; b3=b2; b3=b0; elseif b20 b0=b2; b2=b3; b3=b0; end else if b20 b0=b2; b2=b3; b3=b0; b3=b2; b2=b0; end end a1=b1 a2=b2 a3=b3 tmax=0.5*(a1-a3) t23=0.5*(a2-a3) t12=0.5*(a1-a2) alpha1=0:pi/50:2*pi; alpha2=0:pi/50:2*pi; alpha3=0:pi/50:2*pi; X1=tmax*cos(alpha1)+(a1+a。

7、3)/2; Y1=tmax*sin(alpha1); X2=t23*cos(alpha2)+(a2+a3)/2; Y2=t23*sin(alpha2); X3=t12*cos(alpha3)+(a1+a2)/2; Y3=t12*sin(alpha3); hold on plot(X1,Y1,k); plot(X2,Y2,k); plot(X3,Y3,k); end 五、实例展示下面我用书中的例题向大家展示我的程序例题 1: 材料力学教材 p33 2-7 ( a)平面应力已知 x=50, y=0, xy=20, =50 。选择平面应力分析 x=50 y=0 xy=20 =50 图片一输出结果图片。

8、二 图片三应力圆输出图片四例题 2：材料力学教材 p35 2-14 ( a)已知 x=60； y=30； z=-70；图片五图片六 图片七图片八六、程序缺陷1. 没有生成标准的 EXE文件，不能很好地应用2. 用户输入窗口不够人性化3. 应力圆输出并非圆形七、课程收获在这次上机编程作业中， 虽然花费了很长时间， 但是我感到收获很多。 在此过程中， 我学习了 Matlab 语言并对 Matlab 语言有了更加熟悉的运用， 这将成为我一生的财富， Matlab 语言也必将成为我日后工作的得力助手。在看到自己做出的最终结果以及它运行时的高效， 我不禁有了一丝成就感， 我为能有一次这样锻炼的机会而感到幸运， 希望以后还会更多这样的机会充实自己。 对于材料力学这门课程， 我想也会因为这次程序的编写而记忆的更加深刻， 我一定要争取多掌握一些这方面的知识，然后在工作中充分的应用进去！。