钢环分离锥式无级变速器弹流润滑膜厚的计算
流体润滑(弹流润滑)
几种摩擦的界限常以膜厚比来大致估计：
hmin R R
2 q1 2 q2
式中：hmin——最小公称油膜厚度，m
Rq1 ——接触表面轮廓的均方根偏差，m
Rq2 ——接触表面轮廓的均方根偏差，m
≤0.4,干摩擦 ≤1，边界摩擦； ＝1～3，混合摩擦； ＞3，流体摩擦
不同润滑状态下的摩擦因数
总结：粘压效应和弹性变形效应有利于提高润滑膜的承载能力。
前面所讨论的流体动压润滑理论及计算，是假定两个润滑表面相对运 动时仍保持完全的刚性，未产生弹性变形，这在低副接触时是正确的。但是， 对于高副接触，如齿轮、滚动轴承等，其比压很大，运用流体动压润滑理论 就不再合适了。 低副 （面接触） 润滑表面刚性 流体动压润滑理论 高副 （点、线接触） 润滑表面弹性变形 弹性流体动压润滑理论 定义一：相对运动表面的弹性变形与流体动压作用都对润滑油的润滑性 能起着重要作用的一种润滑状态。 定义二：弹性流体动力润滑是指流体进入在两个相互运动的固体摩擦接 触表面后，受到接触表面产生的巨大接触压力而发生的性状改变，以分 割固体摩擦接触表面，减少摩擦。 从广义上说：凡是表面弹性变形量和最小油膜厚度处在同一量级上的 润滑问题，都属于弹流问题。 弹性流体动压润滑理论是流体动压润滑理论的重要发展，可以说弹性流 体动力润滑是流体膜润滑的一种特殊形式。它主要研究名义上是点、线 接触的摩擦副润滑问题（如齿轮副、滚动轴承等）。
两个弹性圆柱的接触，可等效于一当量圆柱和一刚性平面的接触问题，因此 在弹流润滑的研究中，可以将接触区视为平面。
等效圆柱
刚性平面
为了分析弹性流体动压润滑机理，首先观察一下对偶表面干接触时的情况。 如图所示模型是弹性圆柱体与一刚性平面干接触的情况。在载荷作用下，弹性 圆柱体发生弹性变形，使线接触变成了小面积接触，载荷所造成的接触压力常 称为赫兹压力，其分布情况是在接触区域内成抛物线形分布，中间的压力最高 而至边缘降低为零。
钢环分离锥式无级变速器设计
目录第一章绪论 (1)1.1无级变速器的介绍 (1)1.2摩擦式无级变速器 (1)1.3摩擦式无级变速器运动原理 (1)1.4钢环分离锥式无级变速器的优点 (3)1.5本次课题设计任务 (3)第二章钢环分离锥式无级变速器设计方案 (4)2.1钢环分离锥式无级变速器简图 (4)2.2传动零件尺寸 (4)2.3钢环分离锥式无级变速器受力分析 (5)2.4强度验算 (7)2.4.1恒功率传动情况时 (8)2.4.2恒扭矩传动情况时 (10)2.4.3钢环强度效验计算 (11)第三章钢环分离锥式无级变速器的计算 (13)3.1计算锥轮的尺寸和参数 (13)3.2钢环的设计 (14)3.3轴系零件设计 (14)3.4调速操纵机构设计 (16)3.4.1确定齿轮的参数 (16)3.4.2确定齿条的参数 (17)3.4.3计算螺杆 (16)3.5设计箱体 (18)第四章强度校核 (19)4.1刚换强度验算 (19)4.2校验轴的强度 (20)设计总结 (22)致谢 (23)参考文献 (24)附录：英文文献翻译及原文 (25)摘要钢环分离锥锥轮无级变速器是机械摩擦式的一种变速器，它以钢环为中间原件，以改变主、从动锥轮的工作半径来实现无级变速。
它能实现对称变速而且无需再设加压装，结构简单，时常将这种变速器应用在传动系统的高速级。
首先查找变速器相关资料，了解其传动原理及设计要求和计算公式，选择材料。
通过已知给定参数先求出变速器主要零件钢环和主从锥轮的相关尺寸，再根据已算出的数据和配合关系选定其主要配合原件轴承型号，然后确定锥轮各段长度和大小。
再进行轴的设计，通过公式选取轴的最少直径，再结合与锥轮配合关系确定轴的各段长度及选取键和轴键等相关尺寸，根据设计手册选取有关尺寸的配合公差，选取设计调速操作机构，再由已知的零件尺寸和配合关系，根据设计手册确定箱体和端盖的基本尺寸，其后对轴和钢环进行强度校核，以确定尺寸是否满足要求。
面接触润滑测量系统油膜厚度自动计算方法
面接触润滑测量系统油膜厚度自动计算方法白清华;郭峰;田鹏晖【摘要】To overcome the low efficiency in artificial estimation of the film thickness in an optical slider-on-disc lubricating film tester,a robust algorithm was proposed to recognize interference intensity variation and its order change,which enables automatic measurement of film thickness.Theoretical and experimental results show that the multi-beam interference principle can be applied to the film thickness measurement in the slider-on-disc contact,and the interference intensity profile can be divided into regions Ⅰ and Ⅱ within one film thickness period.Freeman chain code and simple moving average filter technique were used to smooth the data.A normalization algorithm based on threshold value was proposed for determining in which region the data are,and the interference order can be obtained accurately.A template matching method was used to infer the shape of the start and the end areas.This approach was applied to process measured data in the experiments,and results were in good agreement with the oit film thickness calculated by the old method.%为了解决面接触润滑油膜光干涉测量系统人工计算油膜厚度效率低下的问题,提出光强数据的高效处理算法,实现油膜厚度的自动计算.通过理论和实验验证多光束干涉理论应用于滑块-玻璃盘面接触副的可行性.对整周期内的光强分布进行区域定义,应用Freeman方向链码和均值滤波对数据曲线进行分段光滑,设计基于阈值的归一算法,通过模板对比算法,得到光强曲线的始末端形态,自动计算出光强数据所属区域Ⅰ和Ⅱ的数量,实现干涉级次的准确计算;以该方法编写程序应用于实际测量,结果与手动计算的油膜厚度有很好的一致性.【期刊名称】《润滑与密封》【年(卷),期】2017(042)001【总页数】7页(P34-40)【关键词】光干涉;油膜厚度;Freeman方向链码;模板对比【作者】白清华;郭峰;田鹏晖【作者单位】青岛理工大学计算机工程学院山东青岛266520;青岛理工大学机械工程学院山东青岛266520;青岛理工大学机械工程学院山东青岛266520;青岛理工大学机械工程学院山东青岛266520【正文语种】中文【中图分类】TH117摩擦副表面之间形成一定厚度的流体润滑油膜可以减少或避免磨损，降低摩擦阻力。
滚动轴承点接触弹流油膜厚度及摩擦力矩的分析计算
滚动轴承的摩擦力矩问题涉及到弹性力学 、 接触力学 、摩擦 、润滑等学科 , 且各种因素相互影 响 ,相互作用 ,使摩擦力矩的精确分析计算不易进 行 。新的摩擦力矩计算公式可以定量地分析轴承 弹流油膜厚度与摩擦力矩的关系 , 也可以应用于 轴承的设计和工艺 ,提高轴承的性能 , 进而增强机 械设备的可靠性 ,增加设备利用率 。
(1
-
e ) - 0. 68k
(3)
材料参数 G′=αE′
(4)
E1′=
1 21 (-源自ν2 1E1+
1
-
ν2 2
)
E2
速度参数 U ′=ηE0′RUx
(5)
U
π =
(Dp2w
-
Dw2 co sβ)
120Dpw
| ne
-
ni
|
载荷参数W
′=
W
E ′R2x
(6)
W = Fa / Z sinβ
椭圆率 k = 1. 03 ( Ry ) 0. 64
(1. 河南科技大学 机电工程学院 ,河南 洛阳 471003; 2. 洛阳轴承研究所 ,河南 洛阳 471039)
摘要 :点接触弹流润滑的实际计算较为复杂 ,对点接触弹流油膜厚度及摩擦力矩的计算进行了详细分析并给出 了计算实例 。同时对弹流油膜厚度与轴承摩擦力矩的关系进行了探讨 ,提出了新的摩擦力矩经验计算公式 。 关键词 :滚动轴承 ;点接触 ;弹流润滑 ;油膜厚度 ;摩擦力矩 中图分类号 : TH133. 33; TH117 文献标志码 : B 文章编号 : 1000 - 3762 (2008) 04 - 0010 - 03
对试验所得的数据进行整理分析 , 提出了新 的摩擦力矩经验计算公式为
金属带式无级变速器钢带环弹性特性及其造成的能量损失
金属带式无级变速器钢带环弹性特性及其造成的能量损失张武1，2 刘凯1 崔亚辉1 原园1 赵桐1【摘要】摘要:基于弹性力学理论，建立了金属带式无级变速器钢带环伸长模型。
在圆弧段上，钢带环被简化为组合薄壁圆筒模型，在直线段上被简化为单向拉伸模型。
模型计算结果显示:随着传动比的增大，应力应变在主动轮上不断增大，在从动轮上先减小后增大;径向位移在主动轮上先减小后增大，在从动轮上不断增大;当传动比一定时，最内层钢带环的应力应变和径向位移最大，最外层最小;钢带环总的伸长量随着传动比的增大而增大，但当传动比为1时有所减小;钢带环线应变范围是0.04% ～0.08%。
功率损失分析表明:钢带环应变能功率损失范围为27.15～86.73W，当传动比为1时，总功率损失为34.85W。
【期刊名称】中国机械工程【年(卷),期】2012(023)021【总页数】7【关键词】关键词:钢带环;金属带;无级变速器;伸长量0 引言金属带式无级变速器(MB－CVT)作为最成功的无级变速器之一，具有节省燃料、操控简单、驾驶平顺舒适和低排放等优点。
最新研究显示:在相同的工作状态下，装备MB－CVT的汽车较装备自动挡的汽车节油7% ～15%［1］。
MB－CVT 属于摩擦传动系统，因此其传递动力的能力受到一定限制。
Kong等［1］以单个钢带环为基础建立了MBCVT动力学模型，研究了带轮与最内层钢带环以及钢带环之间的相互作用及每层钢带环在传递动力方面的负载分布状况，并且认为一层或多层钢带环过早断裂是变速器失效的一个重要原因。
Kuwabara等［2］在研究中忽略了钢带环之间的相互作用，并将其假设为一个整体。
Kim等［3］完全忽略MB－CVT各部件间的相互作用，且进一步认为两组钢带环和数以百计的金属块是一个连续体。
Akehurst等［4］假设所有的钢带环都具有相同的圆弧段，且每组钢带环是由数层钢带环组合而成，但是对于每层钢带环承载力矩的能力却没有研究。
弹流润滑数值计算方法-赵江
弹流润滑数值计算方法1.1等温线接触弹流润滑数值计算方法与程序等温点接触弹流润滑的基本方程包括：Reynolds 方程、膜厚方程、变形方程、粘压方程、密压方程和载荷平衡方程，主要形式如下Reynolds 方程：xρh 12μ=x p ηρh x 3d )d()d d (d d s （1） 膜厚方程)(+2R+=)(20x υx h x h （2）变形方程∫ex +d ln(2-=x c s x -s x p Eπx υ))()( （3）粘压方程{}])()[.(z p p ηηη000+1+1-679+ln exp = （4）密压方程)..ppρρ71+160+1=0（ （5）载荷平衡方程∫0=d -0ex x x x p ω)( （6）计算程序计算工况参数 节点数：N=30量纲一花起始点坐标：X0=-4.0 量纲一花终止点坐标：XE=1.5 载荷：W=1.0E5综合弹性模量E1=2.2E11 初始粘度：EDA0=0.05 圆柱半径：R=0.05 速度:US=1.5 主程序：FORGRAM LINEEHLCOMMON/COM1/ENDA,A1,A2,A3,Z,HM0,DH/COM2/EDA0/COM4/X0,XE/COM3/E1,PH,B,RDATA PAI,Z,P0/3.14159565,0.68,1.96E8/DATA N,X0,XE,W,E1,EDA0,R,US/130,-4.0,1.5,1.0E5,2.2E11,0.05,0.05,1.5/OPEN(8,FILE=’OUT.DAT’,STA TUS=’UNKNOWN’)W1=W/(E1*R)C W无量纲公式，PH=E1*SQRT(0.5*W1/PAI)C PH为最大hertz接触应力Ph，单位paA1=(ALOG(EDA0)+9.67)C 粘度公式计算时的参数A2=PH/P0C 粘度公式计算时的参数A3=0.59/(PH*1.E-9)C 密度公式计算时参数B=4.*R*PH/E1C B为HERTZ接触区半宽，单位mALFA=Z*A1/P0 C粘度系数G=ALFA*E1 C材料参数U=EDA0*US/(2.*E1*R)C 苏联人提出的入口区分析解中定义的CC1=SQRT(2.*U)AM=2.*PAI*(PH/E1)**2/CC1ENDA=3.*(PAI/AM)**2/8C ENDA量钢化reynolds中的栏目达HM0=1.6*(R/B)**2*G**0.6*U**0.7*W1**(-0.13)C hm0最小膜厚，公式为最小膜厚公式WRITE(*,*)N,X0,XE,W,E1,EDA0,R,USCALL SUBAK(N)CALL EHL(N)STOPENDSUBROUTINE EHL(H)C 弹流润滑子程序DIMENSION X(1100),P(1100),H(1100),RO(1100),EPS(1100),EDA(1100),V(1100)C X P压力、H膜厚、RO密度、EPS雷诺方程中的E、EDA粘度、V弹性变形COMMON/COM1/ENDA,A1,A2,A3,Z,HM0,DH/COM4/X0,XEC ENDA量钢化reynolds中的栏目达, hm0最小膜厚COMMON/COM3/E1,PH,B,RRC E1当量弹性模量E，PH为最大hertz接触应力ph，B为HERTZ接触区半宽，RRMK=1DX=(XE-X0)/(N-1.0)C DX等距节点间距DO 10 I=1,NX(I)=X0+(I-1)*DXIF(ABS(X(I)).GE,1.0)P(I)=0.0C 压力分布为抛物线形式,IF(ABS(X(I)).LT,1.0)C 压力分布为抛物线形式10 CONTINUECALL HREE(N,KK,DX,X,P,H,RO,EPS,EDA,V)CALL FZ(N,P,POLD)C 对应于框图中,调用HREE子程序计算个节点膜厚,密度ro,和粘度EDA，C 并调用子程序VI计算个节点的弹性变形C 调用子程序FZ将P(I)的值赋给POLD,POLD第i个节点的前一次迭代压力值14 KK=19CALL ITER(N,KK,DX,X,P,H,RO,EPS,EDA,V)C 调用ITER子程序进行雷诺方程迭代，重新计算hree压力分布下的节点膜厚/粘度/密度MK=MK+1CALL ERROP(N,P,POLD,ERP)C 调用ERROP子程序计算迭代前后的压力差并输出WRITE(*,*)’ERP=’,ERPIF(ERP.GT.1.E-5.AND.DH.GT.1.E-6)THENC .GT.大于IF(MK.GE.50)THENMK=1DH=0.5*DHENDIFGOTO14ENDIFC 判断ERP和DH是否满足条件IF(DH.LE.1.E-6)WRITE(*,*)’pressure are not convergent’H2=1.E3P2=0.0DO 106 I=1,NIF(H(I).LT.H2)H2=H(I)IF(P(I).GT.P2)P2=P(I)C 找出最小膜厚H2和最大压力P2,量纲一化并输出106 CONTINUEH3=H2*B*B/RRP3=P2*PHC 这个地方还不是很懂,可能是获得最小的值然后量纲计算出有单位的值110 FORMAT(6(1X,E12.6))C 格式化输出1X把输出位置向右跳过n个位置C E12.6以12个字符宽输出指数类型的浮点数，小数部分占6个字符宽120 CONTINUEWRITE(*,*)‘P2,H2,P3,H3=’,P2,H2,P3,H3CALL OUTHP(N,X,P,H)RETURNSUBROUTINE OUTHP(N,X,P,H)C 输出子程序DIMENSION X(N),P(N),H(N)DO 10 I=1,NWRITE(8,20)X(I),P(I),H(I)10 CONTINUE20 FORMA T(1X,6(E12,6,1X))RETURNENDSUBROUTINE HREE(N,DX,X,P,H,RO,EPS,EDA,V)DIMENSION X(N),P(N),H(N),RO(N),EPS(N),EDA(N),V(N)COMMON/COM1/ENDA,A1,A2,A3,Z,HM0,DH/COM2/EDA0/COMAK/AK(0:1100)DATA KK,PAI1,G0/0,0.318309886,1.570796325/IF(KK.NE.0)GOTO3H00=0.03 W1=0.0DO 4 I=1,N4 W1=W1+P(I)C3=(DX*W1)/G0DW=1.-C3C 不是很懂CALL VI(N,DX,P,V)HMIN=1.E3DO 30 I=1,NH0=0.5*X(I)*X(I)+V(I)C 和膜厚公式有关IF(H0.LT.HMIN) HMIN=H0H(I)=H030 CONTINUEIF(KK.NE.0)GOTO 32KK=1DH=0.005*HM0C 逐步减小法H00=-HMIN+HM032 IF(DW.LT.0.0)H00=H00+DHIF(DW.GT.0.0)H00=H00-DHDO 60 I=1,NH(I)=H00+H(I)EDA(I)=EXP(A1*(-1.+(1.+A2*P(I)**Z))C 就是粘压公式RO(I)=(A3+1.35*P(I))/(A3+P(I))EPS(I)=RO(I)*H(I)**3/(ENDA*EDA(I))C 求出了雷诺方程中的一副三拉，为后面雷诺方程的迭代计算做基础60 CONTINUERETURNSUBROUTINE ITER(N,KK,DX,X,P,H,RO,EPS,EDA,V)C 雷诺方程中的迭代DIMENSION X(N),P(N),H(N),RO(N),EPS(N),EDA(N),V(N)COMMON/COMAK/AK(0:1100)DATA PAI1/0.318309886/DO 100 K=1,KKD2=0.5*(EPS(1)+EPS(2))D3=0.5*(EPS(2)+EPS(3))DO 70 I=2, N-1D1=D2D2=D3IF(I.NE.N-1)D3=0.5*(EPS(I+1)+EPS(I+2))D8=RO(I)*AK(0)*PAI1D9=RO(I-1)*AK(1)*PAI1D10=1.0/(D1+D2+(D9-D8)*DX)D11=D1*P(I-1)+D2*P(1+1)D12=(RO(I)*H(I)-RO(I-1)*H(I-1)+(D8-D9)*P(I))*DXP(I)=(D11-D12)*D10IF(P(I).LT.0.0)P(I)=0.0C 计算雷诺差分方程个系数，计算节点压力，并保证节点压力不小于0(这里还要仔细考虑)70 CONTINUECALL HREE(N,DX,X,P,H,RO,EPS,EDA,V)C 调用HREE方程重新计算各节点的膜厚，100 CONTINUERETURNENDSUBROUTINE VI(N,DX,P,V)DIMENSION P(N),V(N)COMMON/COMAK/AK(0:1100)PAI1=0.318309886C=ALOG(DX)DO 10 I=1,NV(I)=0.0DO 10 J=1,NIJ=IABS(I-J)10 V(I)=V(I)+(AK(IJ)+C)*DX*P(J)DO I=1,NV(I)=-PAI1*V(I)ENDDORETURNENDSUBROUTINE SUBAK(MM)COMMON/COMAK/AK(0:1100)DO 10 I=1,MM10 AK(I)=(I+0.5)*(ALOG(ABS(I+0.5))-1.)-(I-0.5)*(ALOG(ABS(I-0.5))-1.)ENDSUBROUTINE FZ(N,P,POLD)DIMENSION P(N),POLD(N)DO 10 I=1,N10 POLD(I)=P(I)RETURNENDSUBROUTINE ERROP(N,P,POLD,ERP)DIMENSION P(N),POLD(N)SD=0.0SUM=0.0DO 10 I=1,NSD=SD+ABS(P(I)-POLD(I))10 SUM=SUM+P(I)ERP=SD/SUMRETURENENDFORGRAM LINEEHLCOMMON/COM1/ENDA,A1,A2,A3,Z,HM0,DH/COM2/EDA0COMMON/COM4/X0,XE/COM3/E1,PH,B,RDATA PAI,Z,P0/3.14159565,0.68,1.96E8/DATA N,X0,XE,W,E1,EDA0,R/130,-4.0,1.5,1.0E5,2.2E11,0.05,0.05/ DATA US/1.5/OPEN(8,FILE=’OUT.DAT’,STATUS=’UNKNOWN’)W1=W/(E1*R)C W无量纲公式，PH=E1*SQRT(0.5*W1/PAI)C PH为最大hertz接触应力Ph，单位paA1=(ALOG(EDA0)+9.67)C 粘度公式计算时的参数A2=PH/P0C 粘度公式计算时的参数A3=0.59/(PH*1.E-9)C 密度公式计算时参数B=4.*R*PH/E1C B为HERTZ接触区半宽，单位mALFA=Z*A1/P0 C粘度系数G=ALFA*E1 C材料参数U=EDA0*US/(2.*E1*R)C 苏联人提出的入口区分析解中定义的CC1=SQRT(2.*U)ENDA=3.*(PAI/AM)**2/8C ENDA量钢化reynolds中的栏目达HM0=1.6*(R/B)**2*G**0.6*U**0.7*W1**(-0.13)C hm0最小膜厚，公式为最小膜厚公式WRITE(*,*)N,X0,XE,W,E1,EDA0,R,USCALL SUBAK(N)CALL EHL(N)STOPENDSUBROUTINE EHL(H)C 弹流润滑子程序DIMENSION X(1100),P(1100),H(1100),RO(1100),EPS(1100),EDA(1100)DIMENSION V(1100)C X P压力、H膜厚、RO密度、EPS雷诺方程中的E、EDA粘度、V弹性变形COMMON/COM1/ENDA,A1,A2,A3,Z,HM0,DH/COM4/X0,XEC ENDA量钢化reynolds中的栏目达, hm0最小膜厚COMMON/COM3/E1,PH,B,RRC E1当量弹性模量E，PH为最大hertz接触应力ph， B为HERTZ接触区半宽，RRMK=1DX=(XE-X0)/(N-1.0)C DX等距节点间距DO 10 I=1,NX(I)=X0+(I-1)*DXC 定位X(I)的位置,对应于框图中的计算个节点X(I),和hertz压力值IF(ABS(X(I)).GE,1.0)P(I)=0.0C 压力分布为抛物线形式,IF(ABS(X(I)).LT,1.0)C 压力分布为抛物线形式10CONTINUECALL HREE(N,KK,DX,X,P,H,RO,EPS,EDA,V)CALL FZ(N,P,POLD)C 对应于框图中,调用HREE子程序计算个节点膜厚,密度ro,和粘度EDA，C 并调用子程序VI计算个节点的弹性变形C 调用子程序FZ将P(I)的值赋给POLD,POLD第i个节点的前一次迭代压力值14 KK=19CALL ITER(N,KK,DX,X,P,H,RO,EPS,EDA,V)C 调用ITER子程序进行雷诺方程迭代，重新计算hree压力分布下的节点膜厚/粘度/密度 MK=MK+1CALL ERROP(N,P,POLD,ERP)C 调用ERROP子程序计算迭代前后的压力差并输出WRITE(*,*)’ERP=’,ERPIF(ERP.GT.1.E-5.AND.DH.GT.1.E-6)THENC .GT.大于IF(MK.GE.50)THENMK=1ENDIFGOTO14ENDIFC 判断ERP和DH是否满足条件IF(DH.LE.1.E-6)WRITE(*,*)’pressure are not convergent’H2=1.E3P2=0.0DO 106 I=1,NIF(H(I).LT.H2)H2=H(I)IF(P(I).GT.P2)P2=P(I)C 找出最小膜厚H2和最大压力P2,量纲一化并输出106CONTINUEH3=H2*B*B/RRP3=P2*PHC 这个地方还不是很懂,可能是获得最小的值然后量纲计算出有单位的值110FORMAT(6(1X,E12.6))C 格式化输出1X把输出位置向右跳过n个位置C E12.6以12个字符宽输出指数类型的浮点数，小数部分占6个字符宽120CONTINUEWRITE(*,*)‘P2,H2,P3,H3=’,P2,H2,P3,H3CALL OUTHP(N,X,P,H)RETURNENDSUBROUTINE OUTHP(N,X,P,H)C 输出子程序DIMENSION X(N),P(N),H(N)DO 10 I=1,NWRITE(8,20)X(I),P(I),H(I)10CONTINUE20FORMAT(1X,6(E12,6,1X))RETURNENDSUBROUTINE HREE(N,DX,X,P,H,RO,EPS,EDA,V)DIMENSION X(N),P(N),H(N),RO(N),EPS(N),EDA(N),V(N)COMMON/COM1/ENDA,A1,A2,A3,Z,HM0,DH/COM2/EDA0/COMAK/AK(0:1100)DATA KK,PAI1,G0/0,0.318309886,1.570796325/IF(KK.NE.0)GOTO3H00=0.03 W1=0.0DO 4 I=1,N4 W1=W1+P(I)C3=(DX*W1)/G0DW=1.-C3C 不是很懂HMIN=1.E3DO 30 I=1,NH0=0.5*X(I)*X(I)+V(I)C 和膜厚公式有关IF(H0.LT.HMIN) HMIN=H0H(I)=H030CONTINUEIF(KK.NE.0)GOTO 32KK=1DH=0.005*HM0C 逐步减小法H00=-HMIN+HM032IF(DW.LT.0.0)H00=H00+DHIF(DW.GT.0.0)H00=H00-DHDO 60 I=1,NH(I)=H00+H(I)EDA(I)=EXP(A1*(-1.+(1.+A2*P(I)**Z))C 就是粘压公式RO(I)=(A3+1.35*P(I))/(A3+P(I))EPS(I)=RO(I)*H(I)**3/(ENDA*EDA(I))C 求出了雷诺方程中的一副三拉，为后面雷诺方程的迭代计算做基础60CONTINUERETURNENDSUBROUTINE ITER(N,KK,DX,X,P,H,RO,EPS,EDA,V)C 雷诺方程中的迭代DIMENSION X(N),P(N),H(N),RO(N),EPS(N),EDA(N),V(N)COMMON/COMAK/AK(0:1100)DATA PAI1/0.318309886/DO 100 K=1,KKD2=0.5*(EPS(1)+EPS(2))D3=0.5*(EPS(2)+EPS(3))DO 70 I=2, N-1D1=D2D2=D3IF(I.NE.N-1)D3=0.5*(EPS(I+1)+EPS(I+2))D8=RO(I)*AK(0)*PAI1D9=RO(I-1)*AK(1)*PAI1D10=1.0/(D1+D2+(D9-D8)*DX)D11=D1*P(I-1)+D2*P(1+1)D12=(RO(I)*H(I)-RO(I-1)*H(I-1)+(D8-D9)*P(I))*DXP(I)=(D11-D12)*D10IF(P(I).LT.0.0)P(I)=0.0C 计算雷诺差分方程个系数，计算节点压力，并保证节点压力不小于0(这里还要仔细考虑) 70CONTINUEC 调用HREE方程重新计算各节点的膜厚，100CONTINUERETURNENDSUBROUTINE VI(N,DX,P,V)DIMENSION P(N),V(N)COMMON/COMAK/AK(0:1100)PAI1=0.318309886C=ALOG(DX)DO 10 I=1,NV(I)=0.0DO 10 J=1,NIJ=IABS(I-J)10 V(I)=V(I)+(AK(IJ)+C)*DX*P(J)DO I=1,NV(I)=-PAI1*V(I)ENDDORETURNENDSUBROUTINE SUBAK(MM)COMMON/COMAK/AK(0:1100)DO 10 I=1,MM10 AK(I)=(I+0.5)*(ALOG(ABS(I+0.5))-1.)-(I-0.5)*(ALOG(ABS(I-0.5))-1.)RETURNENDSUBROUTINE FZ(N,P,POLD)DIMENSION P(N),POLD(N)DO 10 I=1,N10 POLD(I)=P(I)RETURNENDSUBROUTINE ERROP(N,P,POLD,ERP)DIMENSION P(N),POLD(N)SD=0.0SUM=0.0DO 10 I=1,NSD=SD+ABS(P(I)-POLD(I))10 SUM=SUM+P(I)ERP=SD/SUMRETURNENDPROGRAM MAINWRITE(*,*) 'HELLO'WRITE(*,*)1'HELLO'100WRITE(*,*)'HELLO'10STOPEND1.2等温点接触弹流润滑数值计算方法与程序Program pointehlCOMMON/COM1/ENDA,A1,A2,A3,Z,HM0,DHDATA PAI,Z/3.14159265,0.68/DATA N,PH,E1,EDA0,RX,US,X0,XE/33,0.8E9,2.21E11,0.05,0.02,1.0,-2.5,1.5/OPEN(4,FILE=’OUT.DAT’,STATUS=’UNKNOWN’)OPEN(8,FILE=’FILM.DAT’,STATUS=’UNKNOWN’)OPEN(10,FILE=’PRESSURE.DAT’, STATUS=’UNKNOWN’)MM=N-1A1=ALOG(EDA0)+9.67A2=5.1E-9*PHA3=0.59/(PH*1.E-9)U=EDA0*US/(2.*E1*RX)B=PAI*PH*RX/E1W0=2.*PAI*PH/(3.*E1)*(B/RX)**2ALFA=Z*5.1E-9*A1G=ALFA*E1HM0=3.6*(RX/B)**2*G**0.49*U**0.68*W0**(-0.073)ENDA=12.*U*(E1/PH)*(RX/B)**3WRITE(*,*)N,X0,XE,W0,PH,E1,EDA0,RX,USWRITE(4,*)N,X0,XE,W0,PH,E1,EDA0,RX,USWRITE(*,*)’ WAIT PLEASE’CALL SUBAK(MM)CALL EHL(N,X0,XE)STOPENDSUBROUTINE EHL(N,X0,XE)DIMENSION X(65),Y(65),H(4500),RO(4500),EPS(4500),EDA(4500),P(4500),POLD(4500),V(4500) COMMON/COM1/ENDA,A1,A2,A3,Z,HM0,DHDATA MK,G0/1,2.0943951/CALL INITI(N,DX,X0,XE,X,Y,P,POLD)KK=0CALL HREE(N,DX,KK,H00,G0,X,Y,H,RO,EPS,EDA,P,V)14 KK=15CALL ITER(N,KK,DX,H00,G0,X,Y,H,RO,EPS,EDA,P,V)MK=MK+1CALL ERP(N,ER,P,POLD)IF(ER.GT.1.E-5.AND,DH,GT.1.E-6)THENIF(MK.GE.20)THENMK=1DH=0.5*DHENDIFGOTO 14ENDIFIF(DH.LE.1.E-6)WRITE(*,*)’PRESSURE ARE NOT CONVERGENT’CALL OUTPUT(N,DX,X,Y,H,P)RETURNENDSUBROUTINE ERP(N,ER,P,POLD)DIMENSION P(N,N),POLD(N,N)ER=0.0SUM=0.0DO 10 I=1,NDO 10 J=0,NER=ER+ABS(P(I,J)-POLD(I,J))POLD(I,J)=P(I,J)SUM=SUM+P(I,J)10 CONTINUEER=ER/SUMRETURNENDSUBROUTINE INITI(N,DX,X0,XE,X,Y,P,POLD)DIMENSION X(N),Y(N),P(N,N),POLD(N,N)NN=(N+1)/2DX=(XE-X0)/(N-1.)Y0=-0.5*(XE-X0)DO 5 I=1,N,X(I)=X0+(I-1)*DXY(I)=Y0+(I-1)*DX5 CONTINUEDO I=1,ND=1.-X(I)*X(I)DO J=1,NC=D-Y(J)*Y(J)IF(C.LE.0.0)P(I,J)=0.0IF(C.GT.0.0)P(I,J)=SQRT(C)POLD(I,J)=P(I,J)EnddoEnddoReturnEndSubroutine hREE(n,dx,kk,h00,G0,X,Y,H,RO,EPS,EDA,P,V)DATA PAI,PAI1/3.14159265,0.2026423/NN=(N+1)/2CALL VI(N,DX,P,V)HMIN=1.E3DO 30 I=1,NDO 30 J=1,NRAD=X(I)*X(I)+Y(J)*Y(J)W1=0.5*RADH0=W1+V(I,J)IF(H0.LT.HMIN)HMIN=H030 H(I,J)=H0IF(KK.EQ.0)THENKK=1DH=0.01*HM0H00=-HMIN+HM0ENDIFW1=0.0DO 32 I=1,NDO 32 J=1,N32 W1=W1+P(I,J)W1=DX*DX*W1/G0DW=1.-W1IF(DW.LT.0.0)H00=H00+DHIF(DW.GT.0.0)H00=HOO-DHDO 60 I=1,NDO 60 J=1,NH(I,J)=H00+H(I,J)EDA1=EXP(A1*(-1.+(1.+A2*P(I,J))**Z)EDA(I,J)=EDA1RO(I,J)=(A3+1.34*P(I,J))/(A3+P(I,J))60 EPS(I,J)=RO(I,J)*H((I,J)**3/(ENDA*EDA1)DO 70 J=NN+1,NJJ=N-J+1DO 70 I=1,NH(I,J)=H(I,JJ)RO(I,J)=RO(I,JJ)EDA(I,J)=EDA(I,JJ)70 EPS(I,J)=EPS(I,JJ)RETURNENDSUBROUTINE ITER(N,KK,DX,H00,G0,X,Y,H,RO,EPS,EDA,P,V)DIMENSION X(N),Y(N),P(N,N),H(N,N),ro(n,n),EPS(N,N),EDA(N,N)V(N,N) COMMON/COMAK/AK(0:65,0:65)DATA KG1,PAI/0,3.14159265/IF(KG1.NE.0)GOTO2AK00=AK(0,0)AK10=AK(1,0)2 NN=(N+1)/2DO 100 K=1,KKDO 70 J=2,NNJ0=J-1J1=J+1D2=0.5*(EPS(1,J)+EPS(2,J))DO 70 I=2,N-1I0=I-1I1=I+1D1=D2D2=0.5*(EPS(I1,J)+EPS(I,J))D4=0.5*(EPS(I1,J0)+EPS(I,J))D5=0.5*(EPS(I,J1)+EPS(I,J))D8=2.0*RO(I,J)*AK00/PAI**2D9=2.0*RO(I0,J)*AK10/PAI**2D10=D1+D2+D4+D5+D8*DX-D9*DXD11=D1*P(I0,J)*+D2*P(I1,J)+D4*P(I,J0)+D5*P(I,J1)D12=(RO(I,J)*H(I,J)-D8*P(I,J)-RO(I0,J)*H(I0,J)+D9*P(I,J)*DXP(I,J)=(D11-D12)/D10IF(P(I,J).LT.0.0)P(I,J)=0.070 CONTINUEDO 80 J=1,NNJJ=N+1-JDO 80 I=1,N80 P(I,JJ)=P(I,J)CALL HREE(N,DX,KK,H00,G0,X,Y,H,RO,EPS,EDA,P,V)100 CONYINUERETURNENDSUBROUTINE VI(N,DX,P,V)DIMENSION P(N,N),V(N,N)COMMON/COMAK/AK(0:65,0:65)PAI1=0.0202643DO 40 I=1,NDO 40 J=1,NH0=0.0DO 30 K=1,NIK=IABS(I-K)DO 30 L=1,NJL=IABS(J-L)30 H0=H0+AK(IK,JL)*P(K,L)40 V(I,J)=H0*DX*PAI1RETURNSUBROUTINE SUBAK(MM)COMMON/COMAK/AK(0:65,0:65)S(X,Y)=X+SQRT(X**2+Y**2)DO 10 I=1,NXP=I+0.5XM=I-0.5DO 10 J=1,NYP=J+0.5YM=J-0.5A1=S(YP,XP)/S(YM,XP)A2=S(XM,YM)/S(XP,YM)A3=S(YM,XM)/S(YP,XM)A4=S(XP,YP)/S(XM,YP)AK(I,J)=XP*ALOG(A1)+YM*ALOG(A2)+XM*ALOG(A3)+YP*ALOG(A4) 10 AK(J,I)=AK(I,J)RETURNENDSUBROUTINE OUTPUT(N,DX,X,Y,H,P)SIMENSION X(N),Y(N),H(N,N),P(N,N)A=0.0WRITE(8,110)A,(Y(I),I=1,N)DO I=1,NWRITE(8,110)X(I),(H(I,J),J=1,N)ENDDOWRITE(10,110)A,(Y(I),I=1,N)DO I=1,NWRITE(10,110)X(I),(P(I,J),J=1,N)ENDDO110 FORMAT(66(E12.5,1X))RETURNEND。
圆锥滚子轴承挡边弹流润滑膜厚的一种计算方法
・26 ・ 平均卷吸速度
1. 2 无量纲方程的建立
洛 阳 工 学 院 学 报 2000
u=
1 ( u + uvb) 2 vi
(1)
平均卷吸速度确定后 , 便可建立模型中各种无量纲方程 , 方程中使用的无量纲参数如下 : η Rx h 0u w P H= , W = 2 , U = , P =
圆锥滚子轴承挡边弹流润滑膜厚的一种计算方法
邱 明1 ,王中宇1 ,郑培斌2
(1. 洛阳工学院 机电工程系 , 河南 洛阳 471039 ;2. 洛阳工学院 机设工程系 , 河南 洛阳 471039)
摘要 :以弹性流体动力润滑理论为基础 ,建立了圆锥轴承滚子端面与内圈挡边间膜厚的计算模型 。对现有膜 厚的计算方法 (直接迭代法) 进行了改进 ,提出了该模型的复合迭代法 。通过对圆锥轴承挡边膜厚的计算表 明 ,该方法的求解误差小 ,收敛速度快 ,是一种有效的求解方法 。 关键词 :滚动轴承 ; 弹性流体动压润滑 ; 油膜 ; 计算方法 ;建立模型 《中图法》 分类号 : TH133. 331 文献标识码 :A
1. 2. 5 载荷平衡方程 0 . 6 ×10 - 9 pH P 1 + 1 . 7 ×10 - 9 pH P 2 π 3
(5)
(4)
1. 2. 4 无量纲密度方程
κ P ( X , Y) d Xd Y =
Μ
(6)
2 求解方法
2. 1 直接迭代法的不足及改进
传统的直接迭代法是将 Reynolds 方程与弹性变形方程截然分开、 交替求解的方法。这样 ,由压力求膜 厚造成的膜厚误差和由膜厚求压力造成的压力误差互相影响 ,逐渐扩大 ,导致求解失败。研究认为 ,将 Reynolds 方程与弹性变形方程融为一体 ,使 Reynolds 方程等号右边只含未知压力 P , 而不含膜厚 h ,可减少 膜厚误差对压力的影响。具体实施过程如下 : 方程 (2) 采用三点中心差分法 ,并将方程 (3) (4) (5) 代入 ,经化简可得到下面的差分方程组 :
摩擦式机械无级变速器结构设计说明书.1
摩擦式机械无级变速器结构设计目录第一章绪论 (4)§1.1机械无级变速器的发展背景及现状 (4)§1.2机械无级变速器的特征和应用 (4)§1.3机械无级变速器的选用和润滑密封 (5)§1.4本文的主要内容及要求 (7)第二章摩擦无级变速器的机械特性加压装置和调速机构 (8)§2.1机械特性 (8)§2.2调速操纵机构 (9)§2.3加压装置 (9)第三章摩擦式无级变速器设计说明和计算过程 (9)§3.1摩擦机械无级变速器的工作原理 (9)§3.2摩擦无级变速器的特点.............................................. 错误!未定义书签。
§3.3锥轮的设计与计算...................................................... 错误!未定义书签。
§3.4钢环的设计与计算...................................................... 错误!未定义书签。
3.4.1、钢环尺寸和参数的确定........................... 错误!未定义书签。
3.4.2、强度验算 ...................................... 错误!未定义书签。
§3.5轴系的设计 ..................................................................... 错误!未定义书签。
§3.6轴的结构设计 ................................................................. 错误!未定义书签。
第四章主要零件的校核........................................... 错误!未定义书签。