下载文档原格式
目录第一章绪论 (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，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。
摩擦式机械无级变速器结构设计目录第一章绪论 (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轴的结构设计 ................................................................. 错误!未定义书签。
第四章主要零件的校核........................................... 错误!未定义书签。
高副接触弹流润滑条件下的油膜厚度分析一高副接触弹流润滑条件下的油膜厚度分析1 弹流润滑条件下的油膜厚度公式1）线接触弹流润滑条件下的油膜厚度公式线接触弹流润滑油膜厚度公式选用Dowson-Higginson 提出的油膜厚度公式【１】，其最小油膜厚度公式为13.003.0'13.043.07.0054.0min)(65.2wE L R u h ηα= （1-1）式中，h min 为最小油膜厚度，m ；R 是综合曲率半径，21111R R R +=，其中R 1、R 2为两接触体在接触点处的曲率半径，m ；u 是接触点卷吸速度，221u u u +=，其中u 1、u 2为两接触体在接触点处的线速度，m/s ；η0是润滑油在大气压下的粘度，Pa ·s ，；α是粘压系数，m 2/N ；E '是综合弹性模量，)11(211222121'E E E μμ-+-=，其中，μ1、μ2为两接触体的泊松比，E 1、E 2为两接触体的弹性模量，Pa ；L 是接触区域轴向长度，m ；w 是滚动体承受的载荷，N 。
从最小油膜厚度公式可以推导出中心油膜厚度公式为13.003.0'13.043.07.0054.0)(53.3w E L R u h c ηα=（1-2）最小油膜厚度公式的无量纲形式为13.07.054.0min65.2W U G H =（1-3）式中，min H 为无量纲最小油膜厚度，R h H /min min =；G 为无量纲材料参数，'E G α=；U 为无量纲速度参数，RE uU '0η=；W 为无量纲载荷参数，RLE wW '=。
从最小油膜厚度公式可以推导出中心油膜厚度公式的无量纲形式为13.07.054.053.3WU G H c =（1-4） 2）点接触弹流润滑条件下的油膜厚度公式点接触弹流润滑油膜厚度公式选用Hamrock-Dowson 提出的油膜厚度公式【2】，其最小油膜厚度公式为)1()(63.368.0073.0117.0'493.049.068.00min k e w E R u h ----=αη （1-5）式中，min h 为最小油膜厚度，m ；R 是综合曲率半径，21111R R R +=，其中1R 、2R 为两接触体在接触点处的曲率半径，m ；u 是接触点卷吸速度，221u u u +=，其中1u 、2u 为两接触体在接触点处的速度，m/s ；η0是润滑油在大气压下的粘度，Pa ·s ；α是粘压系数，m 2/N ；E '是综合弹性模量，)11(211222121'E E E μμ-+-=，其中，1μ、2μ为两接触体的泊松比，1E 、2E 为两接触体的弹性模量，Pa ；w 是滚动体承受的载荷，N 。
一高副接触弹流润滑条件下的油膜厚度分析1 弹流润滑条件下的油膜厚度公式1）线接触弹流润滑条件下的油膜厚度公式线接触弹流润滑油膜厚度公式选用Dowson-Higginson 提出的油膜厚度公式【１】，其最小油膜厚度公式为13.003.0'13.043.07.0054.0min)(65.2wE L R u h ηα= （1-1） 式中，h min 为最小油膜厚度，m ；R 是综合曲率半径，21111R R R +=，其中R 1、R 2为两接触体在接触点处的曲率半径，m ；u 是接触点卷吸速度，221u u u +=，其中u 1、u 2为两接触体在接触点处的线速度，m/s ；η0是润滑油在大气压下的粘度，Pa ·s ，；α是粘压系数，m 2/N ；E '是综合弹性模量，)11(211222121'E E E μμ-+-=，其中，μ1、μ2为两接触体的泊松比，E 1、E 2为两接触体的弹性模量，Pa ；L 是接触区域轴向长度，m ；w 是滚动体承受的载荷，N 。
从最小油膜厚度公式可以推导出中心油膜厚度公式为13.003.0'13.043.07.0054.0)(53.3w E L R u h c ηα=（1-2） 最小油膜厚度公式的无量纲形式为13.07.054.0min65.2W U G H =（1-3） 式中，min H 为无量纲最小油膜厚度，R h H /min min =；G 为无量纲材料参数，'E G α=；U 为无量纲速度参数，RE uU '0η=；W 为无量纲载荷参数，RLE wW '=。
从最小油膜厚度公式可以推导出中心油膜厚度公式的无量纲形式为13.07.054.053.3WU G H c =（1-4） 2）点接触弹流润滑条件下的油膜厚度公式点接触弹流润滑油膜厚度公式选用Hamrock-Dowson 提出的油膜厚度公式【2】，其最小油膜厚度公式为)1()(63.368.0073.0117.0'493.049.068.00min k e w E R u h ----=αη （1-5）式中，min h 为最小油膜厚度，m ；R 是综合曲率半径，21111R R R +=，其中1R 、2R 为两接触体在接触点处的曲率半径，m ；u 是接触点卷吸速度，221u u u +=，其中1u 、2u 为两接触体在接触点处的速度，m/s ；η0是润滑油在大气压下的粘度，Pa ·s ；α是粘压系数，m 2/N ；E '是综合弹性模量，)11(211222121'E E E μμ-+-=，其中，1μ、2μ为两接触体的泊松比，1E 、2E 为两接触体的弹性模量，Pa ；w 是滚动体承受的载荷，N 。
