gear_blankgeometryparameters.asv

用于分析螺旋锥齿轮接触分析的MATLAB源码

                                                     %-程序地位：子 程 序
                                                 %-程序功能：计算右旋大轮的轮坯几何参数
                                                 %-编程时间：2008年8月26日
                                                 %-编 程 者：樊红卫 谷霁红
%***********************************1.定义M文件主函数的输出参数、函数名、输入参数****************************************%
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        小轮齿数
                   %Z_2        大轮齿数
                   %m_t        端面模数
                   %beta_m0    节锥母线中点处的螺旋角
                   %bigsigema0 两轮轴线的夹角
                   %cosi_k     工作齿高系数
                   %cosi_t     全齿高系数
                   %f0         齿高系数
                   %c          顶隙系数
                   %K          刀盘直径系数
                   %theta_f20  大轮齿根角
%**********************************2.读入基本参数,对每个参数赋值,并将角度的单位从度转换到弧度******************************%                   
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.计算轮坯几何尺寸（求派生几何参数）*********************************************%
fid=fopen('E:\项目\编程\gear\result_GearBlank.txt','W'); %事先创建一个存储大轮齿坯计算结果的文本文件
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); %小轮的节锥角
biggama_2=bigsigma-biggama_1; %大轮的节锥角
cosi_h=0.39*(1-(Z_1*cos(biggama_2))/(Z_2*cos(biggama_1)));%轮齿法面内齿高方向的变位系数(大轮取为负，小轮取为正）
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)); %大轮节锥顶点到根锥顶点的距离
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]; %创建一个存放输出参数的列数组
fprintf(fid, '%12.8f\n',array_1); %将计算结果存放到已创建好的文本文件中(小数点后保留8位，输出数据的总宽度不超过12位）
fclose(fid) %关闭存储文件
%***********************************************************结 束*********************************************************%
end