%-程序地位：子 程 序
                                                 %-程序功能：計算右旋大輪的輪坯幾何參數(shù)
                                                 %-編程時間：2008年8月26日
                                                 %-編 程 者：樊紅衛(wèi) 谷霽紅
%***********************************1.定義M文件主函數(shù)的輸出參數(shù)、函數(shù)名、輸入?yún)?shù)****************************************%
function [h_k,h_t,d_2,biggama_1,biggama_2,cosi_h,L_e,b,A_m,smallgama_2,D_c2,h_f2,Z_R2]=Gear_BlankGeometryParameters()
global Z_1  Z_2  m_t  beta_m0  bigsigma0  cosi_k  cosi_t  f0  c  K  theta_f20
% 定義（全局）變量：%Z_1        小輪齒數(shù)
                   %Z_2        大輪齒數(shù)
                   %m_t        端面模數(shù)
                   %beta_m0    節(jié)錐母線中點處的螺旋角
                   %bigsigema0 兩輪軸線的夾角
                   %cosi_k     工作齒高系數(shù)
                   %cosi_t     全齒高系數(shù)
                   %f0         齒高系數(shù)
                   %c          頂隙系數(shù)
                   %K          刀盤直徑系數(shù)
                   %theta_f20  大輪齒根角
%**********************************2.讀入基本參數(shù),對每個參數(shù)賦值,并將角度的單位從度轉(zhuǎn)換到弧度******************************%                   
load input_gear_blankdata_fhw.txt 

Z_1=input_gear_blankdata_fhw(1);     
Z_2=input_gear_blankdata_fhw(2);
m_t=input_gear_blankdata_fhw(3);
beta_m0=input_gear_blankdata_fhw(4);
beta_m=beta_m0*pi/180;
bigsigma0=input_gear_blankdata_fhw(5);
bigsigma=bigsigma0*pi/180;
cosi_k=input_gear_blankdata_fhw(6);
cosi_t=input_gear_blankdata_fhw(7);
f0=input_gear_blankdata_fhw(8);
c=input_gear_blankdata_fhw(9);
K=input_gear_blankdata_fhw(10);
theta_f20=input_gear_blankdata_fhw(11);
theta_f2=theta_f20*pi/180;
%******************************************3.計算輪坯幾何尺寸（求派生幾何參數(shù)）*********************************************%
fid=fopen('D:\1.齒輪研究\4.項目\2.編程\gear\result_GearBlank.txt','W'); %事先創(chuàng)建一個存儲大輪齒坯計算結(jié)果的文本文件

h_k=cosi_k*m_t; %工作齒高
h_t=cosi_t*m_t; %全齒高
d_2=Z_2*m_t; %大輪的分度圓直徑
biggama_1=atan(Z_1/Z_2); %小輪的節(jié)錐角
biggama_2=bigsigma-biggama_1; %大輪的節(jié)錐角
cosi_h=0.39*(1-(Z_1*cos(biggama_2))/(Z_2*cos(biggama_1)));%輪齒法面內(nèi)齒高方向的變位系數(shù)(大輪取為負，小輪取為正）
L_e=d_2/(2*sin(biggama_2));%外錐距
b=min(10*m_t,0.3*L_e); %齒寬
A_m=L_e-0.5*b; %中點錐距
smallgama_2=biggama_2-((h_t-((0.5*h_k)-(cosi_h*m_t)))/L_e); %大輪的根錐角（機床根角、大輪安裝角）
D_c2=sqrt((2*(K^2)*(L_e^2))-(A_m^2)*(2-(sin(beta_m))^2))+A_m*sin(beta_m); %大輪加工刀盤的名義直徑                         
h_f2=(f0+c+cosi_h)*m_t;%大輪齒根高
Z_R2=(L_e-h_f2/tan(theta_f2))*(sin(biggama_2)/tan(biggama_2-theta_f2)-cos(biggama_2)); %大輪節(jié)錐頂點到根錐頂點的距離

array_1=[h_k;h_t;d_2;biggama_1;biggama_2;cosi_h;L_e;b;A_m;smallgama_2;D_c2;h_f2;Z_R2]; %創(chuàng)建一個存放輸出參數(shù)的列數(shù)組
fprintf(fid, '%12.8f\n',array_1); %將計算結(jié)果存放到已創(chuàng)建好的文本文件中(小數(shù)點后保留8位，輸出數(shù)據(jù)的總寬度不超過12位）
fclose(fid) %關(guān)閉存儲文件
%***********************************************************結(jié) 束*********************************************************%
end