alpha=0:0.01:pi/2;
t1=286*sin(2*alpha)/2;
m1=286*(cos(alpha)-cos(alpha).^3);
x1=633+286*cos(2*alpha)/2;
[AX,H1,H2]=plotyy(x1,t1,x1,m1);
set(H2,'LineStyle','--');
set(AX,'XTick',[]);
set(AX,'YTick',[]);
hold on
t2=400*sin(2*alpha)/2;
m2=400*(cos(alpha)-cos(alpha).^3);
x2=846+400*cos(2*alpha)/2;
[AX,H1,H2]=plotyy(x2,t2,x2,m2);
hold off;
set(H2,'LineStyle','--')
set(AX,'Xlim',[0 1100]);
set(AX,'Ylim',[0 200]);
set(AX(1),'XTick',0:100:1100);
set(AX(2),'XTick',[]);
set(AX,'YTick',0:50:200);
set(get(AX(1),'Ylabel'),'String','剪应力');
set(get(AX(2),'Ylabel'),'String','有效力矩');
xlabel('主应力');
title('剪应力与有效力矩图');
