基本信息

---

下载

   源码介绍

---

clc;
clear all;
req=[];

%%%%%%材料参数
%装药材料参数

CD=56/1000;  %装药内切圆半径,m

dc=0.5*pi/180;   %%%角度上的微元,δθ
mc=floor(pi/dc);  %角度上的微元数量

jj=75*pi/180;  %%%下底与腰的夹角


for i=1:1:mc;  %周向 
      c=i*dc;
      C=c*180/pi;
     if  c<=jj/2;   
         b1=CD/2/cos(c);
         b2=CD/2/cos(c);
     elseif  c>jj/2 && c<=jj;  
         b1=CD/2/cos(c);
         b2=CD/2/cos(jj-c);
     elseif c>jj && c<=pi/2-jj/2
         b1=CD/2/cos(c);
         b2=CD/2/cos(c-jj); 
     elseif c>pi/2-jj/2 && c<=pi/2 jj/2 
         b1=CD/2/cos(pi-jj-c);
         b2=CD/2/cos(c-jj);       
     elseif c>pi/2 jj/2  &&  c<=pi-jj 
           b1=CD/2/cos(pi-jj-c);
           b2=CD/2/cos(pi-c);
     elseif c>pi/2 jj/2 && c<=pi-jj/2; 
          b1=CD/2/cos(c-(pi-jj));
          b2=CD/2/cos(pi-c);
     elseif c>pi-jj/2;   
          b1=CD/2/cos(pi-c);
          b2=CD/2/cos(pi-c);
     end       
     
%   R(i)=b;
  Cr(i)=C;
X1(i)=b1*cos(c);
Y1(i)=b1*sin(c);
X2(i)=b2*cos(pi c);
Y2(i)=b2*sin(pi c);
end


figure(1);
plot(X1,Y1);
hold on
plot(X2,Y2);