摘要
摆线针轮行星传动具有传动比范围大、结构紧凑、可靠性高和寿命长等显著特点,因而获得了广泛的应用,研究也不断深入。该传动不仅广泛应用于通用传动领域,而且在微机械、机械人传动装置、精密机械传动、超小型传动、宇航设备,测量仪器、住宅智能化和高技术设备等方面有诱人的应用潜力。
在摆线针轮行星传动家族中,FA型摆线针轮行星传动变速器是一种新的传动装置。与一般的摆线针轮行星传动相比,它具有体积小、重量轻、传动比范围大、传动比的多样性、寿命长、刚度大、回转精度高、精度保持稳定、效率高、传动平稳等一系列的优点。该传动设计采用了许多先进的理念和技术,应用一种创新的结构形式,因此该种FA型针摆传动有效地克服了传统的结构不足,提高了传递的功率和容量,提高了传递的精度。
论文对此进行了系统的综述和分析,并提出了论文的主要研究内容。针对该传动的结构特点,论文从理论上和实际应用两方面研究并解决了以下几个问题:
1.在摆线针轮行星传动摆线轮和针齿齿面受力分析理论方面,提出了一种更加适用于工程实际的齿面有隙啮合受力分析的方法。该方法可以比较精确地计算在整个传递过程中针齿和摆线轮齿面接触力的大小和变化范围。该方法既有效地克服了按标准齿形进行理论受力分析计算误差大的缺点,又证明文献[1]中提出的理论只计算了整个传递过程中针齿和摆线轮齿面接触力的其中一个位置,是本文方法的一种特例。
2.通过对摆线轮和针齿齿面受力分析有限元计算,证明了本文提出的受力分析方法的正确性。两种受力分析计算结果在针齿的受力、同时接触针齿、接触应力的数值分布上均呈现了一致性,且两者最大接触应力的误差在5%左右。
3.在销孔式输出机构受力分析方面,首次提出了一种采用有隙啮合原理进行柱销与摆线轮上的柱销孔之间接触力的计算理论与方法,并推导出一套完整的计算公式,利用该理论比较准确地获得了柱销在整个转动过程中受力的变化区间,以及各个针齿的参与传递扭矩区间。计算结果表明了本文的计算方法更加接近于工程实际。
4.齿面修形理论与方法是国外保密的核心技术,本文在调研各种修形方法的基础上,基于齿面接触状态最佳的原则,首次提出了“反弓齿廓”概念,并给出了利用“负移距+正等距”简单组合修形方法获得反弓齿廓的条件,反弓齿廓初始间隙的计算方法,并采用优化理论获得了最佳的反弓齿廓所需要的“负移距+正等距”组合修形量。受力分析计算结果表明,最佳的反弓齿廓的接触状态可以达到最佳的受力状态。
5.在对FA型针摆传动几何回差和动态回差深入分析的基础上,针对高精传动的特点,本文提出了采用“负等距+正移距”的修形方法可以获得“弓背齿廓”的方法,而“弓背齿廓”的回差可以通过选择不同的修形量来满足。本文还提出了通过搜索计算获得满足许用回差的“弓背齿廓”修形量的计算方法。
6.在FA型摆线针轮行星传动变速器的研究方面,利用本文所提出的理论和研究成果,采用优化设计理论,对FA型针摆传动中的结构主要参数进行优化设计,编制了相关参数优化软件,并以FA45-59型号为例进行了产品设计,并由合作厂生产出样机。整机性能实际试验的结果表明了本文研究方法的正确性。
The cycloid drive has the essential advantages such as wide range of transmission ratio, compact structure, high reliability and long working life, so it gets broad application and its research goes deep constantly. The drive has been not only applied to traditional transmission field but also has alluring application in the aspects of micro machine, robot gear, precision machine transmission, super mini cycloid drive, astronaut equipment, measurement apparatus, tenement intelligence and high-tech equipment.In the family of the cycloid drive, the FA cycloid-pin wheel planetary drive reducer is a new transmission device. Compared with the common cycloid drive, it has a series of advantages such as small volume, lightweight, wide range of transmission, diversity transmission, long life time, high rigidity, high slewing precision, stable precision, high efficiency and stable transmission. The design of the drive adopts a lot of advanced theories and technology methods and uses a kind of innovation structure, so the FA drive has overcome the defect of the traditional structure, enhanced the transmission efficiency, volume, and the transmission precision.It has summarized and analyzed the FA cycloid-pin wheel planetary drive systematically and put forward the main research content in the paper. Aiming for the structure character, the paper has studied and solved the following problems in the two aspects of both the theory and the practical application.1. In the aspect of the force analysis of the cycloid and pin wheel of the cycloid-pin wheel planetary drive, a force analysis method of mesh with initial clearance is put forward, which is more accord with the actual engineering. The method can calculate the contact force and its change range between the pin gear and the cycloid during the whole transmission more precisely. The method has overcome efficiently the deficiency of the big calculation error in the force analysis of the standard gear profile, also has proved that the one position of the touch force between the gear and the cycloid during whole transmission which put forward by literature [1] is a special case of the paper.2. Though the finite element calculation of the force analysis of the cycloid and the pin gear surface, it has been proved that the method put forward by the paper is correct. The two force analysis calculation results are consistent in the pin gear force, the touch pin gear at the same time and the touch stress value distribution, and the biggest error of touch stress between them is less than 5 percent.3. In the force analysis of the pin-hole-output mechanism, based on the principle of mesh with initial clearance, a method has been put forward to calculate the touch stress between the pin and the pin-hole in the cycloid disk for the first time, and a set of calculation formulas has been deduced too. Using the theory the pin force variation range during whole transmission and the pin gear transmitting force area has been obtained correctly. The calculation result shows that the calculation method in this paper is more near to the actual engineering.
4. The teeth profile modification theory is the confidential core technology outside. On the base of investigating all kinds of the profile modification methods, and on the principle of optimum teeth touching surface, the concept of Inverse Arch-shaped Teeth Profile is been brought forward for the first time. It has also provided the method to get the Inverse Arch-shaped Teeth Profile with the combining modification method of 'plus equidistant + minus radial-moving' and the calculation method of the initial clearance of the Arch-shaped Teeth Profile. With the optimization theory the optimum inverse arch-shaped teeth profile is obtained. The calculation results about the various modification methods have showed that the touch state of the optimum inverse arch-shaped teeth profile can get the best stress.5. On the base of analyzing the geometry return difference and the dynamic return difference of the FA cycloid drive deeply, aiming for the character of the high precision drive, the Arch Back Teeth profile is obtained through the modification of the minus equidistant and plus radial-moving. While the return difference of the Arch Back Teeth profile can be satisfied by means of selecting different modification value. The calculation method of the modification value of the Arch Back Teeth profile of the allowable difference is also obtained by searching calculation.6. In the aspect of the FA cycloid-pin wheel planetary drive reducer, using the theory and research results put forward by the paper and adopting the optimized design theory, the main structure parameters in FA drive are optimized, the parameter optimization software is programmed, and taking the case of the type FA45-59 the product is designed and the sample machine is produced by the corporate factory. The practical test result of the whole sample function shows that the research methods of the paper are correct.7. According to the research results above, the serial software package of the FA cycloid-pin wheel planetary drive with three cycloids has been developed with the Visual C++ 6.0. The serial software concentrates the study result of the paper. The serial software is openness, and it can help users to design the main parameters of the FA cycloid-pin wheel planetary drive by providing the ordered transmission ratio. The software is able to provide all sets of technology information for cycloid drive manufactures and it can short the exploit cycle of the new products and economize the design costs.
引文
[1] 李力行.摆线针轮行星传动的齿形修正与受力分析.机械工程学报,1986(1):18~23
[2] CYCLOID FA & 1FA series,产品样本[R],日本住友重机械工业会社,1996.
[3] CYCLOID RV series,产品样本[R],日本住友重机械工业会社,1992.
[4] CYCLOID 90 series,产品样本[R],日本住友重机械工业会社,1992.
[5] CYCLOID 200 series,产品样本[R],日本住友重机械工业会社,1998.
[6] CYCLOID 4000 series,产品样本[R],日本住友重机械工业会社,1994
[7] CYCLOID 5000 series,产品样本[R],日本住友重机械工业会社,2000.
[8] CYCLOID 6000 series,产品样本[R],日本住友重机械工业会社,2001.
[9] 徐灏.机械设计手册.第3卷24篇,机械工业出版社,1992:24-92~24-136.
[10] 张国瑞,张展.行星传动技术.上海交通大学出版社,1989.
[11] 张展.实用机械传动设计手册.科学出版社,1994.
[12] 冯澄宙.摆线针轮行星行星传动.人民教育出版社,1981.
[13] 刘学厚,黎巨泉.行星传动设计.北京工业学院出版社,1988:148~194
[14] 王宏猷,石田武,日高照晃等.采用摆线齿轮的K-H-V行星齿轮装置回转传动误差的研究.日本机械学会论文集(C编),1994,60(578):286-293.
[15] 石田武,王宏猷,日高照晃等.采用摆线齿轮的K-H-V行星齿轮装置回转传动误差的研究.日本机械学会论文集(C编),1994,60(578):278~285.
[16] 日高照晃,王宏猷,石田武等.采用摆线齿轮的K-H-V行星齿轮装置回转传动误差的研究.日本机械学会论文集(C编),1994,60(570):281~289.
[17] 丰住滋,岩本信彦等.新系列摆线针轮行星传动减速器.日本住友重机械技报,1979,27(81):218-228
[18] B.H.库德洛甫采夫等著.江耕华,顾永寿译.齿轮减速器德结构与计算.上海:上海科学技术出版社.1982
[19] Birkholz, H. Drehzahlpruefung von Cycloid-Getrieben. IMW-Institutsmitteilung Nr. 2002(27): 67~70.
[20] Manfred Chmurawa, Antoni John. FEM in numerical analysis of stress and displacement distributions in planetary wheel of cycloidal gear. Naa2000, LNCS 1988.Springer-Verlag Berlin Heidelberg 2001: 772~779.
[21] Lehmann M. Berechnung und Messung der Krafte in einen Zykloiden-Kurvenschieben Getriebe. Dissertation. Technische Universtaet Muenchen, 1976.
[22] Lehmann M.摆线传动论文集.精密工程技术与测量技术,1980(7):1~88.
[23] S. J. Chung, H.Hein, T. Hirata, J. Mohr, T.Akashi. A micro cycloid-Gear system fabricated by multiexposure LIGA technique. Microsystem Technologies, 2000(6): 149~153
[24] Hongbin Xu, Quan Shi, Changhui Yang. Study on load carrying capacity of cycloid pinwheel planetary drive with parallel axis [A]. Proceedings of the international conference on mechanical transmissions [C]. Beijing: Mechnic Industry Publish House, 2001: 436~437.
[25] L. Vulkov, J. Wasniewski, P. Yalamov. FEM in Numerical Analysis of Stress and Displacement Distribution in Planetary Wheel of Cycloidal Gear. Manfred Chmurawa and Antoni John, NAA 2000, LNCS 1988, pp. 772~779, 2001.
[26] Chandrupatla T. R, Belegundu A. D. (1991) Introduction to FEM in engineering. Prentice Hall, London.
[27] Chmurawa M. Distribution of loads in cycloidal planetary gear.International Conference" Mechanics'99, Kaunas University, Lithuania, 1999
[28] Chmurawa M, John A, Kokot G. (1999) The influence of numerical model on distribution of loads and stress in cycloidal planetary gear. Proc. International Scientific Colloquium Cax Techniques, Bielefeld, Germany.
[29] Chmurawa M. (1999) Distribution of loads in cycloidal planetary gear. Proc. International Conference, Mechanics"99", Kaunas University, Lithuania.
[30] Chung S. J, Hein H, Hirata T etc. A micro cycloid-gear system fabricated by muti-exposure LIGA technique. The Third International Workshop on High Aspect Ration Microstructure Technology HARMST'99, 1999, 6(2000): 149~153.
[31] 万朝燕,兆文忠,李力行.二齿差摆线针轮减速器针齿壳内曲线参数优化.机械工程学报,2003(6):28~32.
[32] 万朝燕,李力行.摆线针轮行星传动的参数优化.大连铁道学院学报,1992(3):42~47.
[33] 关天民,赵岩.摆线针轮行星传动的齿面接触力的精确分析方法.面向21世纪的工程设计,人民交通出版社,1998:225~228.
[34] 赵岩,关天民.各种修形方式下摆线针轮行星传动的齿侧间隙.面向21世纪的工程设计,人民交通出版社,1998:221~224.
[35] 李力行,关天民,王子孚.大型摆线针轮行星传动的合理结构和齿形.机械工程学报,1988(3):24.28~32.
[36] 李力行,关天民,王子孚.摆线针轮行星传动的计算机辅助设计.大连铁道学院学报,1992(3):23~34.
[37] 关天民,李力行.摆线针轮行星传动减速器转臂轴承的可靠性计算.大连铁道学院学报,1992(1):93~95.
[38] 李力行,薛嘉庆.对摆线针轮行星减速器摆线齿形修正方式的分析方法.东北工学院学报,1981(1).
[39] 李力行,洪淳赫.摆线针轮行星传动中摆线轮齿形通用方程式的研究.大连铁道学院学报,1992(1):7~12.
[40] 刘银远,张致祥,李力行.摆线减速器计算机辅助绘图方法的研究.大连铁道学院学报,1992(1):57~61.
[41] 关天民.摆线针轮行星传动销孔式输出机构受力的准确计算方法.机械传动,2000,24(3):4~5~11.
[42] 关天民,范英,刘彬.修磨摆线轮时出现的三种现象的分析.组合机床与自动化加工技术,2000(6):47~49.
[43] 关天民,万朝燕.三片摆线轮新型针摆传动理论及其受力分析.大连铁道学院学报,1999(3):48~51.
[44] 赵岩,关天民.针摆传动转角组合修形下的受力分析新方法.机械设计与制造,1998(3):45~46.
[45] 关天民.摆线针轮行星传动减速器可靠性评价体系的初步研究.机械传动,2001(1):17~19~23.
[46] 关天民,孙英时.超小型摆线针轮行星传动及其受力分析.机械设计与制造,2001(3).
[47] Li Lixing Guan Tianmin etc. The computer aided design of cycloid drive. CHINESE JOURNAL OF MECHANICAL ENGINEERING, 1990, 3(2): 221~231.
[48] Li Lixing, Guan Tianmin etc. The new optimum tooth profile of the cycloid gear and the CAD of cycloid drive. Ninth world congress on the theory of machines and mechanisms (IFTOMM), August 29/September2, 1995, Millano, Italy: 355~359.
[49] Tianmin Guan, Xiaochun shi, The Accurate Force Analysis of the Pin-hole-output Mechanism in the Cycloid Drive. Proceedings of the International Conference on Mechanical Transmissions (ICMT'2001), April5~9, 2001, Chongqing, China: 645~647.
[50] 马英驹,孙英时,关天民.二齿差针摆行星传动中摆线轮等效代换齿廓的试验研究.大连铁道学院学报,1999(1):19~22~37.
[51] Li Lixing, Guan Tianmin etc. The rational construction and the tooth profile on cycloid disk of a large cycloid gearing. Chinese Journal of Mechanical Engineering, 2(1), 1989: 72~81
[52] 关天民.摆线针轮行星传动中修形所产生的回转误差计算与分析.组合机床与自动化加工技术,2001(10):15~18.
[53] 单丽君,关天民.摆线针轮行星传动动态理论回转误差计算与分析.机械传动,2002,26(2):29~32.
[54] 吴永宽,郑剑云等.摆线针轮行星传动的几何回差分析计算.大连铁道学院学报,1999,20(2):28~32.
[55] 关天民.针摆传动理论啮合效率公式的推导.机械设计与制造,2001(3):48~49.
[56] 林菁,张钰.摆线针轮传动输出转角滞后分析.东北大学学报(自然科学学报),1996,17(1):110~113.
[57] 姚文席.摆线针轮行星减速机的传动误差分析.机械传动,1998,22(4):14~18.
[58] 关天民,张东生.FA针摆传动回转精度分析.机械设计与制造,2004(3):90~91.
[59] Jifang Qu, Zijun An. The research of zero clearance cycloid ball transmission. Journal of Chinese Mechanical Engineering, 7(1), 1994: 17~22.
[60] Yang D. C. H, Blanche J. G. Design and application guidelines foe cycloid drives with machining tolerances. Mechanism and Machinery Theory, 1990, 25(5): 487~501.
[61] 李力行,何卫东等.机器人用高精度RV传动的研究.大连铁道学院学报,1999,20(2):1~11.
[62] 何卫东.机器人用高精度RV减速器中摆线轮的优化新齿形.机械工程学报,2000(3).
[63] 徐永贤,何卫东.王洪等.RV传动刚度计算方法.大连铁道学院学报,1999,20(2):18~23.
[64] 关天民,李力行.摆线针轮行星传动减速器转臂轴承的可靠性计算.大连铁道学院学报,1992(1):93~95.
[65] 关天民.摆线针轮行星传动中最佳修形量的确定方法.中国机械工程,2002,13(10):811~815.
[66] 关天民,雷蕾,孙英时.二齿差摆线针轮行星传动的受力分析.机械传动,2002,25(4):19~23.
[67] Guan Tianmin, Shi Xiaochun. Force Analysis considering Manufacture Error of the Pin-hole-output Mechanism in the Cycloid Drive. Chinese Journal of Mechanical Engineering (English Edition), 15(2), 2002: 142~144.
[68] 李力行.论摆线针轮行星传动新产品开发.大连铁道学院学报,1992(1):1~6.
[69] Y. Zhang, Z. Wu. Offset face gear driver tooth geometry and contact analysis. Transaction of the ASME, 1993, 119(3): 223~227.
[70] 雷蕾.三片摆线轮新针摆传动系列三维造型与样机研制.大连铁道学院硕士毕业论文.2002.
[71] 魏延刚,李力行,兆文忠.摆线针轮行星传动针齿壳的有限元分析.大连铁道学院学报,1992,13(1):53~56.
[72] 关天民,雷蕾,孙英时等.FA型摆线针轮行星传动装置的反求设计.中国机械工程,2003,14(1):68~72.
[73] 关天民,孙英时,于力牛.三片摆线轮减速机针齿与摆线轮啮合实际作用力分析.大连铁道学院学报,2002,23(2):20~23.
[74] Weidong He, Lixing Li, Xin Li, Linda Schmidt. A New Drive with High-Load Capacity and High Transmission Efficiency. Proceedings of DETC'00 ASME 2000 Design Engineering Technical Conferences Computers and Information in Engineering Conference, Baltimore, Maryland, September 10~13, 2000.
[75] 关天民,雷蕾.超小型摆线针轮行星传动减速器参数确定及其绘图软件的开发.机械传动,2002,26(4):46~49.
[76] 陈其廉,刘晓波,胡胜海.摆线针轮行星减速器的CAD设计.哈尔滨船舶工程学院学报,1994,15(3):53~59.
[77] Jianjun Zhou, Zichen Chen. Study and Performance Test of Full Complement Cycloidal Ball Reducer with Ceramic Balls as its Gear Teeth. CHINESE JOURNAL OF AERONAUTICS, 14(4), 2001 : 245~251.
[78] 曲继方,安子军.密珠摆线球齿传动研究[J].机械工程学报,1994,7(1):17~22.
[79] 安子军,曲志刚,张荣贤.摆线钢球传动的齿形综合研究[J].机械工程学报,1996,9(5):33~38.
[80] Jianjun Zhou, Zichen Chen. Performance of a Full Complement Ball cycloid Drive [C]. ICME' 2000 Shanghai, 2001: 436~437.
[81] Jianjun Zhou, Zichen Chen. Computerized Design of Cycloidal Gear Drive with Improved Gear Tooth Modification [C]. ICMT'2001, Chongqing 2001: 167~170.
[82] 周建军,陈子辰.摆线钢球行星传动接触疲劳强度研究 [A]. The 4st International Conference Frontiers of Design and Manufacturing [C], 2000: 515~518.
[83] Hidetsugu T, et al. Fundamental analysis of cycloid ball reducer(1st report)[J]. Journal of the Japan Society for Precision Engineering, 1988, 54(1): 67~74.
[84] 刘茂武,孟惠荣.重载摆线针轮传动的接触问题分析及三维有限元计算.机械传动,1994,18(1):11~16.
[85] 喻洪流,齐功相.重载摆线针轮行星传动针齿非赫兹弹性接触修形的研究.机械传动,2002,26(1):52~54~87.
[86] 裘建新.长幅外摆线行星传动学研究及齿廓圆弧优化替代.上海交通大学学报,2004,38(6):727~931.
[87] Litvin F. L. Pinhao Feng. Computerized design and generation of cycloidal gearing. Mech Mach Theory, 31(7), 1996: 891~910.
[88] 赵菊娣.圆弧摆线齿轮油泵的齿廓曲线参数优化设计.FLUID MACHINERY,2003,31(7):18~21~32.
[89] 喻洪流.摆线减速器轮齿啮合的数学模型及求解.机械设计,2004,21(3):14~16.
[90] Sundt E V. Cycloidal reduction gear [P]. America Patent: 3037400. 1962-06.
[91] 王厚宽,王卫荣.链条少齿差行星减速机的研究及其CAD[J].机械工程学报,1994,30(4):76~82.
[92] 戚厚军,孙涛,刘纪岩等.2K-V型摆线针轮减速器传动系统的动态特性分析.机械传动,2004,28(2):13~16.
[93] 芦新春,王明强.摆线针轮行星减速器的稳健优化设计.华东船舶工业学院学报,2004,18(5):61~65.
[94] 李海涛,魏文军.摆线齿锥齿轮齿面接触区的计算机辅助分析.中国农业大学学报,2004,9(5): 45~50.
[95] Lewicki. D. G, Ballarini. R. Effect of rim thickness on gear crack propagation path. Transaction of the ASME, 1997, 119(3): 48~52.
[96] Wang S. M, Morse I. E. Torsional response of a gear train system. ASME (B), 1972, 194(2): 583~594.
[97] Kahrman A, Singh R. Non-Linear Dynamics of a Geared Rotor-Bearing System with Mutiple Clearances. Journal of Sound and Vibration, 1991, 144(3).
[98] Tuplin W. A. Gear Load Capacity. Sir Isacc Pitman & Sans Ltd, 1962.
[99] Simon V. Load and stress distributions in spur and helical gears ASME Journal of Mechanisms Transmissions, and Automations in Design, 1988,110: 197~202.
[100] Schmidt G. R. Optimum tooth profile correction of helical gears. Third international Power Transmission and Gearing Conference, San Francisco, 1980: DET-110.
[101] 林清安.Pro/ENGINEER零件设计基础篇上、下,高级篇上,零件组合.北大宏博改编.北京大学出版社,2000.
[102] 李为真,伍北成,高国宏.最新Pro/Engineer2000i应用实例教程.冶金工业出版社,2000.
[103] 王福军,张志民,张师伟.AutoCAD 2000环境下C/Visual C++应用程序开发教程.北京希望电子出版社,2000.
[104] 刘良华,朱东海.AutoCAD 2000ARX开发技术.清华大学出版社,2000.
[105]Jon Bates,Tim Tompkins.实用Visual C++6.0教程:何健辉,董方鹏译.清华大学出版社,2000.
[106] 宋斌,陈玉亭.Visual C++6.0教程.北京希望电子出版社,2000:1~114.
[107] 陈火红.MSC.Marc接触分析培训教程.MSC.Software中国.2001.
[108] 陈火红.Marc有限元实例分析教程.机械工业出版社,2002.
[109] 黄炎.工程弹性力学.清华大学出版社.
[110] 张允真,曹富新.弹性力学及其有限元法.中国铁道出版社.
[111] 方新.I-DEAS产品三维设计指南.机械工业出版社.
[112] 刘清吉,张文奖,谢忠佑.I-DEAS进阶.科学出版社,2001.
[113] 王瑁成,邵敏.有限单元法基本原理和数值方法.清华大学出版社,1997.
[114] 钟万勰,张洪武,吴承伟.参变量变分原理及其在工程中的应用.科学出版社,1997.
[115] 吴昌华,张军.轮轨接触的有限元研究.大连铁道学院学报,1997(4):25~30.
[116] 佟维,吴昌华.SS_7型电力机车牵引齿轮系统的有限元三维接触分析.机车电传动,1999(6):11~13~23
[117] N. Ahmadi, L. M. Keeer and T. Mura. Non-Hertzian Contact Stress Analysis For An Elastic Half Space Normal and Sliding Contact. Int. J. Solids Structure, 19(4), 1993.
[118] M. Heydari and R. Gohar. The Influence of Axial Profile on Pressure Distribution in Radically Load Rollers. Journal Mechanical Engineering Science, 1979, 21 (6).
[119] J. W. Kannel. Comparison between Predicted and Measured Axial Pressure Distribution between Cylinders. ASME, 1984.
[120] Krishna P. Single and Burton Paul. Numerical Solution of Non-Hertzian Elastic Contact Problems. ASME, 1974.
[121] Litvin F. L. Gear Geometry and Applied Theory [M]. Prentice Hall Englewood Cliiffs, 1994: 22~98.
[122] Heyliger P. R, Reddy J. N. A Mixed Computational Algorithm for Plan Elastic Contact Problems [J]. Computers & Structures, 1987, 26(5): 621~634.
[123] Francavilla A, Zienkiewicz O. C. A Note on Numerical Computation of Elastic Contact Problem [J]. Int. J Num Meth Eng, 1975, 19(3): 813~924.
[124] Kyuichiro Washizu. Variational Methods in Elasticity and Plasticity [M]. Pergam on Press: 1982: 3~168.
[125] Zefeng Wen, Xuesong Jin, Weihuang Zhang. ELASTIC-PLASTIC FINITE ELEMENT ANALYSIS OF THREE-DEMENSIONAL CONTACT-IMPACT AT RAIL JOINT. CHINESE JOURNAL OF MECHANICAL ENGINEERING, Vol. 16, No. 4, 2003: 411~416.
[126] Bahattin KANBER, Ibrahim H. GUZELBEY, Ahmet ERKLIG. BOUNDARY ELEMENT ANALYSIS OF CONTACT PROBLEMS USING ARTIFICIAL BOUNDARY NODE APPROACH. ACTA MECHANICA SINICA, Vol. 19, No. 4, August 2003: 347~354.
[127] Karami G, Fenner RT. A two dimensional BEM method for thermo-elastic body forces contact problems. In: Boundary Element Ⅸ. Ⅴ. 01. 2: Stress Analysis Applications Uni of Stuttgart. Germany: Springer-Verlag, 1987.
[128] 杨生华.齿轮接触有限元分析.计算力学学报,2003,20(2):189~194.
[129] Barlam D, Zahavi E. The reliability of solutionsin in contact problems [J]. Comp & Struct, 1999, 70: 35~45.
[130] Zheng Pei. NEW METHOD TO RESEARCH MESHED GEARS AND NEW TOOTH PROFILE. The 11st International Conference on Geometry and Graphics, 1-5 August, Guangzhou, China, 3888~391.
[131] 杨琪,李乔.运用面向对象的Visual C++开发纯Win32桥梁软件.中南公路工程,1999,24(2):59~61.