time-vari-meshing-stiffness

斜齿轮时变啮合刚度计算matlab程序。原理见李瑰贤的《宽斜齿轮副啮合刚度计算及扭振特性的研究》。

function [a0amnphn] = dfseries（xx yy tol）
%     [a0amnphn] = dfseries（xx yy tol）
    %tol:total terms of （amnphn）
    %a0:mean term a0/2
    %amnphn:amplitude and phase of an*cos（wt）+bn*sin（wt）
    if nargin == 2
        tol=5;
    end
    xxn=length（xx）;
    ll=（xx（xxn）-xx（1））/2;
    a0=1/ll*trapz（xxyy）;
    an=zeros（tol1）;bn=an;
    for ii=1:1:tol
        an（ii）=1/ll*trapz（xxyy.*cos（ii*pi*xx/ll））;
        bn（ii）=1/ll*trapz（xxyy.*sin（ii*pi*xx/ll））;
    end
    amn=（an.^2+bn.^2）.^0.5.*sign（bn）;phn=atan（an./bn）;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         551  2015-03-25 12:44  dfseries.m
     文件        2454  2015-03-25 12:21  dynmeshk.m