《计算机在材料科学与工程中的应用实验设计与指导 教学课件 ppt 作者 叶卫平 位错应力场模拟》由会员分享,可在线阅读,更多相关《计算机在材料科学与工程中的应用实验设计与指导 教学课件 ppt 作者 叶卫平 位错应力场模拟(7页珍藏版)》请在金锄头文库上搜索。
1、MATLAB在位错应力场的计算机模拟,利用弹性力学理论得刃位错的应力场,化为极坐标形式,设x=r cos,y=r sin,带入化简得:,以 xx为例,模拟计算前对公式进行处理,设:,讨论在与定值r上的所有点的应力情况。设m=A/r(r,A都为常数)方程可简化为:,编写M程序,theta=0:pi/180:2*pi; r=sin(theta).*(2+cos(2*theta); x1=r.*cos(theta); y1=r.*sin(theta); figure(1) plot(x1,y1, r); hold on; x2=r.*cos(theta); y2=-r.*sin(theta); pl
2、ot(x2,y2); hold off,编写M程序如下: theta=0:pi/180:2*pi; r=sin(theta).*cos(2*theta); x1=r.*cos(theta); y1=r.*sin(theta); figure(1) plot(x1,y1,r); hold on; x2=r.*cos(theta); y2=-r.*sin(theta); plot(x2,y2); axis(-pi/3 pi/3 -1 1) hold off,化简yy得: yy =msincos2,化简xy得:xy=coscos2,编写M程序如下: theta=0:pi/180:2*pi; r=cos(theta).*cos(2*theta); x1=r.*cos(theta); y1=r.*sin(theta); plot(x1,y1); hold on; x2=-r.*cos(theta); y2=r.*sin(theta); plot(x2,y2,r); axis(-pi/3 pi/3 -1 1) hold off,