摩尔库伦强度准则

admin

[mw_shl_code=applescript,true]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)


[/mw_shl_code]