发明名称 基于混合维模型的精密偶件配合间隙泄漏特性预测方法
摘要 本发明公开了一种基于混合维模型的精密偶件配合间隙中流体泄漏特性预测方法。根据精密偶件的尺寸参数和尺寸偏差使用样条曲面模拟精密偶件配合面的整数维曲面分量,然后,通过归一化分数维表面高度场特征函数模拟精密偶件配合面的分数维分量,并计算分数维分量的高度修正系数,使其满足特定的形状误差,建立满足尺寸偏差和形状误差的精密偶件配合面混合维模型,最后,通过计算两精密偶件混合维表面的间隙截面积,预测精密偶件配合间隙的流体泄漏特性。本发明采用规则整数维曲面与不规则分数维曲面相叠加的精密偶件配合面的非理想几何模型计算精密偶件的配合间隙,准确地预测了尺寸偏差和形状误差综合作用下精密偶件配合间隙中流体泄漏特性。
申请公布号 CN102722656A 申请公布日期 2012.10.10
申请号 CN201210192812.1 申请日期 2012.06.08
申请人 浙江大学 发明人 刘振宇;裘辿;谭建荣
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 林怀禹
主权项 1.一种基于混合维模型的精密偶件配合间隙中流体泄漏特性预测方法,其特征在于包括以下步骤:(1)采用标准样条曲面模拟精密偶件规则的整数维外形特征,根据零件的尺寸参数和尺寸偏差建立精密偶件的整数维分量p<sub>id</sub>(u,w);(2)基于归一化分数维表面高度场特征函数在(u,w)变量空间中模拟精密偶件的分数维分量h<sub>fd</sub>(u,w),其生成方法为:<maths num="0001"><![CDATA[<math><mrow><msub><mi>h</mi><mi>fd</mi></msub><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>w</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><msub><mi>n</mi><mi>max</mi></msub></munderover><msup><mi>&gamma;</mi><mrow><mrow><mo>(</mo><mi>D</mi><mo>-</mo><mn>3</mn><mo>)</mo></mrow><mi>n</mi></mrow></msup><mo>&CenterDot;</mo><mo>{</mo><mi>cos</mi><msub><mi>&phi;</mi><mrow><mi>m</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>-</mo><mi>cos</mi><mo>[</mo><mfrac><mrow><mn>2</mn><mi>&pi;</mi><msup><mi>&gamma;</mi><mi>n</mi></msup><msqrt><msup><mi>x</mi><mn>2</mn></msup><mo>+</mo><msup><mi>y</mi><mn>2</mn></msup></msqrt></mrow><mi>L</mi></mfrac><mi>cos</mi><mo>[</mo><msup><mi>tan</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mfrac><mi>y</mi><mi>x</mi></mfrac><mo>)</mo></mrow><mo>-</mo><mfrac><mi>&pi;m</mi><mi>M</mi></mfrac><mo>]</mo><mo>+</mo><msub><mi>&phi;</mi><mrow><mi>m</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>]</mo><mo>}</mo><mi>u</mi><mo>,</mo><mi>w</mi><mo>&Element;</mo><mo>[</mo><mn>0,1</mn><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,x和y是整数维分量上任意一点的坐标值,计算结果h<sub>fd</sub>(u,w)为该点在分数维粗糙表面上的高度值;D(2<D<3)为分数维粗糙表面的分数维维数;γ(γ>1)为尺度参数,表征了粗糙表面的频谱密度;M表示构造表面时叠加轮廓峰的数量;φ<sub>m,n</sub>为一随机相位,在[0,2π]范围内分布;L为整数维分量的单方向最大长度;n为累加的次数,最大累加次数为<img file="FDA00001742073000012.GIF" wi="358" he="135" />方括号为取整运算,n<sub>0</sub>为整数维分量的单方向最大采样点数;(3)对步骤(2)生成的归一化分数维表面高度场进行几何结构评定,如果当前分数维分量的几何结构参数F<sub>p</sub>不满足形状误差要求,则需要对其进行高度修正,计算其高度修正系数K,否则,令高度修正系数K为1;设精密偶件的形状误差为<img file="FDA00001742073000013.GIF" wi="90" he="55" />则高度修正系数的计算方法如下式:<maths num="0002"><![CDATA[<math><mrow><mi>K</mi><mo>=</mo><mfrac><msubsup><mi>F</mi><mi>p</mi><mi>m</mi></msubsup><msub><mi>F</mi><mi>p</mi></msub></mfrac><mo>=</mo><mfrac><msubsup><mi>F</mi><mi>p</mi><mi>m</mi></msubsup><mrow><msub><mi>h</mi><mi>fd</mi></msub><msub><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>w</mi><mo>)</mo></mrow><mi>max</mi></msub><mo>-</mo><msub><mi>h</mi><mi>fd</mi></msub><msub><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>w</mi><mo>)</mo></mrow><mi>min</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>(4)以高度修正系数K为偏移幅值、以整数维表面上每点的单位法向量V<sub>hd</sub>为偏移方向,对整数维分量p<sub>id</sub>上的每个点进行偏移,偏移基数为该点在归一化分数维表面高度场h<sub>fd</sub>中的高度值,得到精密偶件表面的混合维几何模型为:p<sub>hd</sub>(u,w)=p<sub>id</sub>(u,w)+V<sub>hd</sub>Kh<sub>fd</sub>(u,w)  u,w∈[0,1]  (3)(5)使用间隙面积来衡量精密偶件的配合精度,某一截面的有效截面积由该截面的轮廓线计算求得;计算精密偶件各个截面的间隙截面积,其中最小间隙截面积A<sub>p</sub>为:<img file="FDA00001742073000021.GIF" wi="1556" he="237" />其中,X<sub>内i,j</sub>和X<sub>外i,j</sub>分别代表内、外配合表面第i个截面中第j个点的X轴坐标值;Z<sub>内i,j</sub>和Z<sub>外i,j</sub>分别代表内、外配合表面第i个截面中第j个点的Z轴坐标值;(6)根据精密偶件的混合维表面模型特征参数和计算出的配合精度间隙截面积,推导出在尺寸偏差、配合精度、相对运动速度、油液的温度和偏心量因素影响下精密偶件配合间隙的“压差-剪切流”油液泄漏量:<maths num="0003"><![CDATA[<math><mrow><mi>Q</mi><mo>=</mo><mfrac><mrow><msubsup><mi>A</mi><mi>p</mi><mn>3</mn></msubsup><mi>&Delta;p</mi><mo>[</mo><mn>1</mn><mo>+</mo><mn>1.5</mn><msup><mrow><mo>(</mo><mfrac><mi>e</mi><mrow><msub><mi>&delta;</mi><mi>R</mi></msub><mo>-</mo><msub><mi>&delta;</mi><mi>r</mi></msub></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>]</mo></mrow><mrow><mn>48</mn><msup><mi>&pi;</mi><mn>2</mn></msup><mi>l</mi><msub><mi>&mu;</mi><mn>0</mn></msub><msup><mi>e</mi><mrow><mo>-</mo><mi>&lambda;</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow></msup><msup><mrow><mo>(</mo><mi>R</mi><mo>+</mo><mfrac><mrow><msub><mi>&delta;</mi><mi>R</mi></msub><mo>+</mo><msub><mi>&delta;</mi><mi>r</mi></msub></mrow><mn>2</mn></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>&PlusMinus;</mo><mfrac><mrow><mi>U</mi><msub><mi>A</mi><mi>p</mi></msub></mrow><mn>2</mn></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,Δp为精密偶件两端的压差,e为实际偏心量,l为环形缝隙长度,μ<sub>0</sub>为温度为t<sub>0</sub>时的油液动力粘度,λ为油液的粘温系数,t为当前的油液温度,R为名义尺寸半径,δ<sub>R</sub>和δ<sub>r</sub>分别为内径和外径的尺寸偏差,U为相对运动速度。
地址 310027 浙江省杭州市西湖区浙大路38号