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); plot(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,