新计算机在材料科学与工程中的应用实验设计与指导 教学课件 叶卫平 位错应力场模拟.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θ * * 在线教务辅导网： 教材其余课件及动画素材请查阅在线教务辅导网 QQ:349134187 或者直接输入下面地址： 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([-pi/3 pi/3 -1 1]) hold off 化简σyy得： σyy =msinθcos2θ 编写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