clear
clf
%大主应力sig1，围压sig3
sig1=[114
229
362];
sig3=[100
200
300];
%计算莫尔圆
n=length(sig1);%莫尔圆个数
ci=(sig1+sig3)/2;%圆心
ri=(sig1-sig3)/2;%半径
rici=ri.*ci;
ci2=ci.*ci;
aci=sum(ci);
ari=sum(ri);
arici=sum(rici);
aci2=sum(ci2);
u=(n*arici-aci*ari)/(n*aci2-aci*aci);
v=(ari*aci2-aci*arici)/(n*aci2-aci*aci);
b=u/(1-u^2)^0.5;
a=v*(1+b^2)^0.5;
r=corrcoef(ci,ri);%相关系数
c=a;%粘聚力
fai=atan(b)*180/pi;%摩擦角
%绘制莫尔圆
k=1;
for i=1:n
alpha=0:pi/20:pi; %角度[0,2*pi]
R=(sig1(i)-sig3(i))/2; %半径
x(:,i)=R*cos(alpha)+(sig1(i)+sig3(i))/2;
y(:,i)=R*sin(alpha);
val1(:,k)=x(:,i);
k=k+1;
val1(:,k)=y(:,i);%莫尔圆数据
k=k+1;
plot(x(:,i),y(:,i),'-')%绘制莫尔圆
hold on
end
%绘制抗剪包络线
for j=1:1:sig1(end)
x1(j,:)=j;
y1(j,:)=b*j+a;
plot(x1,y1,'-')
end
val2=[x1,y1];%切线数据
%出图设置
set(gcf,'color','white','units','inches','position',[0 0 6 3.6]);%设置图片尺寸：6×4.8英寸
axis([0 sig1(end) 0 max(y1)]);%设置坐标范围
font_style='helvetica';%设置文字字体
font_size=8;%设置文字字号
set(gca,'FontSize',font_size,'FontName',font_style);%字体字号赋值
xlabel('\sigma / kPa','FontWeight','bold') ;
ylabel('\tau / kPa','FontWeight','bold') ;%设置横纵坐标
cstr=sprintf('c=%3.3f kPa\nf=%3.3f °\nR=%3.6f \n',c,fai,r(2));%显示c、f、r
hh=annotation('textbox',[0.3 0.65 0.1 0.1],'LineStyle','-','color',[0 0 0],'edgecolor',[1 1 1]);
set(hh, 'string',cstr, 'fontsize', font_size-1);
评论3