摩尔库伦强度准则
clc;clear;
close all;
m=15;
n=100;
phi = 40*pi/180;
p=linspace(0,200,m);
%用p、θ表示q
for i = 1:m
theta = linspace(0,2*pi,n);
cos1=cos(theta); %第一卦限转换公式中σ1中的余弦值
cos2=cos(theta-2*pi/3);%第一卦限转换公式中σ2中的余弦值
cos3=cos(theta+2*pi/3);%第一卦限转换公式中σ3中的余弦值
for k = 1:n
if theta(k)<=pi/3 %σ1≥σ2≥σ3第一卦限
cosmax(k)=cos1(k); %最大主应力转换公式中的余弦值
cosmin(k)=cos3(k); %最小主应力转换公式中的余弦值
elseif (theta(k)<=2*pi/3)&&(theta(k)>pi/3) %σ2≥σ1≥σ3第二卦限
cosmax(k)=cos2(k); %最大主应力转换公式中的余弦值
cosmin(k)=cos3(k); %最小主应力转换公式中的余弦值
elseif (theta(k)<=3*pi/3)&&(theta(k)>2*pi/3) %σ2≥σ3≥σ1第三卦限
cosmax(k)=cos2(k); %最大主应力转换公式中的余弦值
cosmin(k)=cos1(k); %最小主应力转换公式中的余弦值
elseif (theta(k)<=4*pi/3)&&(theta(k)>3*pi/3) %σ3≥σ2≥σ1第四卦限
cosmax(k)=cos3(k); %最大主应力转换公式中的余弦值
cosmin(k)=cos1(k); %最小主应力转换公式中的余弦值
elseif (theta(k)<=5*pi/3)&&(theta(k)>4*pi/3) %σ3≥σ1≥σ2第五卦限
cosmax(k)=cos3(k); %最大主应力转换公式中的余弦值
cosmin(k)=cos2(k); %最小主应力转换公式中的余弦值
elseif (theta(k)<=6*pi/3)&&(theta(k)>5*pi/3)%σ1≥σ3≥σ2第六卦限
cosmax(k)=cos1(k); %最大主应力转换公式中的余弦值
cosmin(k)=cos2(k); %最小主应力转换公式中的余弦值
end
end
q=3*p(i)*ones(1,n)*sin(phi)./(cosmax*(1-sin(phi))-cosmin*(1+sin(phi)));%引入莫尔库伦强度准则
s1(i,:)=p(i)*ones(1,n)+2/3.*q.*cos1;
s2(i,:)=p(i)*ones(1,n)+2/3.*q.*cos2;
s3(i,:)=p(i)*ones(1,n)+2/3.*q.*cos3;
end
surf(s1,s2,s3)
view(120,30)
页:
[1]