Matlab找出最优钻压转速配合

function [A] = promote(d,tooth,land,KR,Af,M,lamda,C2,CH,CP,Cb,Cr,tt,n1,n2,n3,hf)
%根据钻机设备条件，确定转速的允许范围，求出不同钻压、转速配合时的钻进成本，
% 从中找出最低的最优钻压、转速配合；
%
%输入：  d = 钻头直径-根据所选钻头确定钻压影响系数Z1、Z2;
%        tooth = 齿型-输入1为铣齿钻头，2为镶齿钻头；
%        land = 适用底层类型-其中1~4为软地层，系列号为1；5~8为中地层，系列号为2，9~12为硬底层
%                  9~12为硬底层，系列号为3；用来确定转速影响系数a1、a2、C1；
%                  例：铣齿钻头适用于中硬地层的21型钻头就是land=5;
%                         适用地层      系列号      类型     land
%                           软            1          1        1
%                           软            1          2        2
%                           软            1          3        3
%                           软            1          4        4
%              
%                           中            2          1        5
%                           中            2          2        6
%                           中            2          3        7
%                           中            2          4        8
%        KR = 地层可钻性系数；
%        Af = 研磨性系数；
%        M = 门限钻压；
%        lamda = 转速指数；
%        Z1 = 钻压影响系数；
%        Z2 = 钻压影响系数；
%        a1 = 钻头类型决定的系数；
%        a2 = 钻头类型决定的系数；
%        C1 = 牙齿磨损减慢系数；
%        C2 = 牙齿磨损系数；
%        CH = 水力净化系数；
%        CP = 压差影响系数；
%        Cb = 钻头成本 元/只；
%        Cr = 钻机作业费 元/小时；
%        tt = 起下钻时间；
%        n1 = 转速；
%        n2 =                     
%        n3 = 
%        hf = 牙齿磨损量；
%  输出： A= [n;
%        a1*n+a2*n^3;
%             R;
%            R/F;
%            Wopt;
%             S;
%             tf;
%             J;
%             Hf;
%            Cpm]         表示不同转速时的最优钻压及其工作指标；
d=input('Enter the d=');

%%     确定钻压影响系数Z1、Z2;
X=[159 0.0198 5.5;
    171 0.0187 5.6;
    200 0.0167 5.94;
    220 0.0160 6.11;
    244 0.0148 6.38;
    251 0.0146 6.44;
    270 0.0139 6.68;
    311 0.0131 7.15;
    350 0.0124 7.56];
[i,j]=find(X==d);
x=X(i,:);
Z1=x(2);                 
Z2=x(3);

%%     确定转速影响系数a1、a2、C1；
tooth=input('Enter the tooth=');
if tooth==1
    land=input('Enter the land=');
    if land==1 || land==2
        a1=2.5;                        
        a2=1.088*10^(-4);
        C1=7;
    elseif land==3 || land==4
        a1=2.0;
        a2=0.870*10^(-4);
        C1=6;
    elseif land==5
        a1=1.5;
        a2=0.653*10^(-4);
        C1=5;
    elseif land==6 || land==7
        a1=1.2;
        a2=0.522*10^(-4);
        C1=4;
    elseif land==8
        a1=0.9;
        a2=0.392*10^(-4);
        C1=3;
    elseif land==9
        a1=0.65;
        a2=0.283*10^(-4);
        C1=2;
    else
        a1=0.5;
        a2=0.218*10^(-4);
        C1=2;
    end
else
    a1=0.5;
    a2=0.218*10^(-4);
    C1=2;
end

%%   输入地层各参数和作业标准
KR=input('Enter the KR=');
Af=input('Enter the Af=');
M=input('Enter the M=');
lamda=input('Enter the lamda=');
C2=input('Enter the C2=');
CH=input('Enter the CH=');
CP=input('Enter the Cp=');
Cb=input('Enter the Cb=');
Cr=input('Enter the Cr=');
tt=input('Enter the tt=');
n1=input('Enter the n1=');
n2=input('Enter the n2=');
n3=input('Enter the n3=');
hf=input('Enter the hf=');

%%   根据公式求出各参数并输出不同转速时的最优钻压及其工作指标
tE=Cb/Cr+tt;
E=(C1/C2)*hf-((C1-C2)/C2^2)*log(1+C2*hf);
F=hf+(C1/2)*hf^2;
n=[n1 n2 n3];
R=(tE*Af*(a1*n+a2*n.^3))./Z1;
Wopt=Z2/Z1+R/F-sqrt((R/F).*(R/F+Z2/Z1-M));
tf=((Z2-Z1*Wopt)./(Af*(a1*n+a2*n.^3))).*(hf+(C1*hf^2)/2);
J=CH*CP*KR*(Wopt-M).*n.^lamda;
S=(Af*(a1*n+a2*n.^3))./(Z2-Z1*Wopt);
Hf=(J./S).*E;
Cpm=(Cr./(J.*E)).*(tE*S+F);
A=[n1 n2 n3;(a1*n+a2*n.^3);R;R/F;Wopt;S;tf;J;Hf;Cpm];
end

作者：TianSL