计算机在材料科学与工程中的应用实验设计与指导 叶卫平 位错应力场模拟新.pptVIP

MATLAB在位错应力场的计算机模拟 利用弹性力学理论得刃位错的应力场 编写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 化简σxy得：σxy=cosθcos2θ * * An extra half plane of atom is inserted Slip  化为极坐标形式，设x=r cosθ，y=r sinθ，带入化简得： 以σ xx为例，模拟计算前对公式进行处理，设： 讨论在与定值r上的所有点的应力情况。设m=A/r（r，A都为常数）方程可简化为： 编写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([