楔横轧阶梯轴成形过程的有限元仿真与分析
|
摘要
楔横轧经过几十年的发展，其实践经验技术已达到很高水平。但由于楔横轧成形机理十分复杂，加之影响因素众多，理论上某些方面还有待于进一步的探讨和研究。而有限元模拟技术有助于解决这些问题，但由于有限元模拟技术在楔横轧成形模拟中的应用正处于起步阶段，其研究面较窄，人们对于阶梯轴的楔横轧成形仍采用传统的实践实验法进行研究。因此，进行楔横轧阶梯轴成形过程的有限元仿真与分析，对于楔横轧的理论发展和技术推广有着重大的意义。
本文采用有限元分析软件ANSYS／LS-DYNA，首次成功模拟了铝合金AL1100阶梯轴的楔横轧轧制全过程。并根据研究的需要，本文设计、分析了四个模型。借助于ANSYS／LS-DYNA的通用后处理器，得到了轧制过程中的应变分布云图，对轧件横截面、纵截面上的切向应变、径向应变、轴向应变、及等效应变的分布特征作了详细探讨，并给出了合理的解释。借用ANSYS／LS-DYNA的时间历程后处理器，输出了等效应变、轧制力的时间历程数据，经绘图、对比分析，总结出了轧制过程中等效应变、轧制力的变化特征，详细探讨了成形角α、展宽角β、断面收缩率ψ三工艺参数对等效应变、轧制力的影响。对进一步揭示楔横轧制过程中的轧件变形特征、应变特征、轧制力特征，指导冷楔横轧的实际生产具有一定的参考价值。
With the development in past years, practice and experiment technology of cross
wedge rolling has been to a high level. Yet the deforming mechanism of cross wedge rolling is quite complicated, with multitudinous factors affecting, there are still questions in theory to need discussing and researching further. Fortunately, technology of finite element simulation can help to resolving these problems. However, application of the technology used in simulating deformation of cross wedge rolling is now in initial stage, studying range of which is limited. Up to now, technicians still adopt conventional method of practice and experiment to study shaft deformation by cross wedge rolling. So simulating and analyzing shaft deformation process by cross wedge rolling has very important sense for theory developing and technology popularizing.
This article successfully simulated the whole deformation process of shaft by cross wedge rolling, where material of workpiece is AL1100, with the finite element software ANSYS/LS-DYNA, totally designed and analyzed four element models for studying. Using the general post processor of ANSYS/LS-DYNA, the article showed pictures of strain distribution of workpiece during deformation, and analyzed strain distribution characters of cross section and longitudinal section, with reasonable and detail describing. In virtue of the time history post processor of ANSYS/LS-DYNA, the article gained data of strains and rolling forces of time history, by plotting and contrasting, summarized varying characters of strain and rolling forces in rolling process, further more, studied the effect on strain and rolling forces of these three technics factors, forming angle, stretching angle and section reduction. To sum up, this article has definite reference value for further discovering characters of workpiece deformation, strain distribution and rolling forces in cross wedge rolling, instructing practice and production.
本文采用有限元分析软件ANSYS／LS-DYNA，首次成功模拟了铝合金AL1100阶梯轴的楔横轧轧制全过程。并根据研究的需要，本文设计、分析了四个模型。借助于ANSYS／LS-DYNA的通用后处理器，得到了轧制过程中的应变分布云图，对轧件横截面、纵截面上的切向应变、径向应变、轴向应变、及等效应变的分布特征作了详细探讨，并给出了合理的解释。借用ANSYS／LS-DYNA的时间历程后处理器，输出了等效应变、轧制力的时间历程数据，经绘图、对比分析，总结出了轧制过程中等效应变、轧制力的变化特征，详细探讨了成形角α、展宽角β、断面收缩率ψ三工艺参数对等效应变、轧制力的影响。对进一步揭示楔横轧制过程中的轧件变形特征、应变特征、轧制力特征，指导冷楔横轧的实际生产具有一定的参考价值。
With the development in past years, practice and experiment technology of cross
wedge rolling has been to a high level. Yet the deforming mechanism of cross wedge rolling is quite complicated, with multitudinous factors affecting, there are still questions in theory to need discussing and researching further. Fortunately, technology of finite element simulation can help to resolving these problems. However, application of the technology used in simulating deformation of cross wedge rolling is now in initial stage, studying range of which is limited. Up to now, technicians still adopt conventional method of practice and experiment to study shaft deformation by cross wedge rolling. So simulating and analyzing shaft deformation process by cross wedge rolling has very important sense for theory developing and technology popularizing.
This article successfully simulated the whole deformation process of shaft by cross wedge rolling, where material of workpiece is AL1100, with the finite element software ANSYS/LS-DYNA, totally designed and analyzed four element models for studying. Using the general post processor of ANSYS/LS-DYNA, the article showed pictures of strain distribution of workpiece during deformation, and analyzed strain distribution characters of cross section and longitudinal section, with reasonable and detail describing. In virtue of the time history post processor of ANSYS/LS-DYNA, the article gained data of strains and rolling forces of time history, by plotting and contrasting, summarized varying characters of strain and rolling forces in rolling process, further more, studied the effect on strain and rolling forces of these three technics factors, forming angle, stretching angle and section reduction. To sum up, this article has definite reference value for further discovering characters of workpiece deformation, strain distribution and rolling forces in cross wedge rolling, instructing practice and production.
引文
[1] 胡正寰．楔横轧技术的现状与展望[J]．锻压技术，1995(5)：25-27
[2] 王宝雨，胡正寰．板式楔横轧轧件成形过程几何模拟[J]．钢铁1999(10)：27-31
[3] 宋玉泉．连续局部塑性成形的发展前景[J]．中国机械工程，2000，11(1-2)：65-67
[4] Yaomion Dong, Kaveh A. Tagavi and Michael R. Lovell.Analysis of interracial slip in cross-wedge rolling: a numerical and phenomenological investigation[J].Journal of Materials Processing Technology, 2000, 97(1-3):44-53
[5] 胡正寰，张康生，王宝雨．楔横轧理论与应用[M]．北京：冶金工业出版社，1996：13-25
[6] 胡正寰，许协和，沙德元．斜轧与楔横轧—原理、工艺及设备[M]．北京：冶金工业出版社，1981：14-17
[7] 曲建国．小变形楔横轧件加工可能性及其质量控制的研究[D]．长春：吉林工业大学经济管理学院，1984
[8] 徐春国，王景梁，任广升．楔横轧成形直角台阶精确截齐曲线的计算[J]．塑性工程学报，1999，6(2)：93-98
[9] 张康生，胡正寰．楔横轧精确直角台阶轧齐曲线[J]．锻压技术，1996(6)：27-33
[10] 张康生，胡正寰．楔横轧轧齐基本理论研究[J]．锻压技术，1996(4)：30-34
[11] 张康生，胡正寰．直角密集台阶轧齐曲线[J]．塑性工程学报，1996，3(1)：44-49
[12] 束学道，张康生，胡正寰，杨翠苹．楔横轧力能参数影响因素分析[J]．重型机械，2002(4)：29-33
[13] Yaomin Dong, Michael Lovell and Kaveh Tagavi .Analysis of interfacial slip in cross-wedge rolling: an experimentally verified finite-element model[J]. Journal of Materials Processing Technology, 1998(80-81):273-281
[14] Qiang,Li,Michael R.Lovell,William Slaughter,Kaveh Tagavi. Investigation of
the morphology of internal defects in cross wedge rolling[J]. Journal of Materials Processing Technology, 2002(125-126):248-257
[15] 马振海，杨翠苹，胡正寰．楔横轧轧件端头凹心影响因素地研究[J]．锻压技术，2002(1)：29-31
[16] 梁继才，任广升等．空心件楔横轧旋转条件的分析[J]．农业机械学报，1995，26(1)：97-101
[17] 梁继才，任广升等．空心件楔横轧压扁失稳的实验研究[J]．吉林工业大学学报，1994，24(4)：78-81
[18] 崔波，梁继才，张昕，刘国平．空心件楔横轧成型过程的应力应变分析[J]．农业机械学报，2001，32(1)：102-104
[19] 颜世公，李玉京．冷楔横轧的实验研究[J]．锻压技术，1997，5：31-34
[20] 权修华，董定福．冷楔横轧工艺实验与研究[J]．锻压机械，1996(4)：24-26
[21] 董定福，权修华，曹诗倬．冷楔横轧中几个工艺问题的探讨[J]．锻压技术，1996(6)：34-36
[22] 白志斌，李东平等．汽车刹车蹄片轴的偏心精密成形[J]．吉林工业大学自然科学学报，2000，30(4)：21-23
[23] 任广升．小型拖拉机阶梯轴类锻坯楔横轧成形[J]．农业机械学报，1993，24(2)：77-80
[24] 吕增亭．楔横轧轴类零件的矫正[J]．机械工程师，1998，12：23-23
[25] 卫原平，彭颖红，阮雪榆．金属塑性成形过程的计算机模拟系统[J]．上海交通大学学报，1996，30(3)：132-137
[26] 洪深泽，胡积威，陈建学．金属成形过程的刚塑性有限元数值模拟[J]．合肥工业大学学报(自然科学版)，1995，18(1)：22-27
[27] 周思柱，X．Q．Li．金属成形的三维有限元模拟系统的发展与展望[J]．锻压技术，1998(4)：46-50
[28] 刘光宇，虞跃生，万贤毅．计算机模拟技术在金属塑性成形中的应用[J]．汽车科技，2000(6)：34-37
[29] 鹿守理，赵辉，张鹏．金属塑性加工的计算机模拟[J]．轧钢，1997(4)：54-57
[30] 谢配良，盛伟，游庆．金属塑性成形过程的计算机模拟[J]．锻压机
械．1997(4)：45-48
[31] 单德彬，吕炎，王真．金属体积成形过程三维刚塑性有限元模拟技术的研究[J]．塑性工程学报，1997，4(3)：98-102
[32] Lee N K,Yoon J H and Yang D Y. Finite Element Analysis of Large Deformation by Automatic Renoding as a Weak Remeshing Technique[J]. Int.J.Mech.Sci, 1992, 34(4):255-273
[33] Dyduch M,Habraken A M and Cescotto S.Automatic Adaptive Remeshing for Numerical Simulation of Metal Forming[J]. Computer Methods in Applied Mechanics and Engineering,1992(101):283-298
[34] 王克智，杜庆柏，马世宁．我国塑性加工摩擦学研究展望[J]．中国表面工程，1995(1)：33-37
[35] 张士宏．塑性加工技术的科学化与中国塑性加工技术的发展[J]．世界科技研究与发展，2001(5)：5-11
[36] 王纪武，刘庄，王本一．体积成形的三维有限变形弹塑性有限元法的研究[J]，塑性工程学报，2003，7(1)：41-44
[37] 胡忠．塑性有限元模拟技术的最新进展[J]．塑性工程学报，1994，1(3)：3-12
[38] 寇淑清，博沛富，杨慎华，胡平，李运兴．三维大变形弹塑性有限元分析与边界条件的处理[J]．塑性工程学报，1998，5(4)：65-70
[39] 闫洪，鲍乐，王美艳，包忠诩．弹塑性有限元法在金属塑性加工中的应用[J]．模具技术，2000(5)：12-14
[40] 孟凡中．弹塑性有限变形理论和有限元方法[M]．北京：清华大学出版社，1985：57-211
[41] 王祖城，汪家才．弹性和塑性理论及有限单元法[M]．北京：冶金工业出版社，1983：1-54
[42] 刘涛，杨凤鹏．精通ANSYS[M]．北京：清华大学出版社，2002：1-170
[43] 博嘉科技．ANSYS融会与贯通[M]．北京：中国水利水电出版社，2002：126-165
[44] 马振海，胡正寰，杨翠萍，束学道．楔横轧展宽段的变形特征与应力应变分析[J]．北京科技大学学报，2002，24(3)：309-312
[45] Yaomin Dong, Kaveh A. Tagavi, Michael R. Lovell and Zhi Deng. Analysis of stress in cross wedge rolling with application to failure[J]. International Journal of Mechanical Sciences, 2000(42): 1233-1253
[46] 黄忠耀，李冬梅，王景先．Pro/ENGINEER 2001基础训练教程[M]．北京：清华大学出版社，2002：368-390
[47] 孙江宏，黄小龙．Pro/ENGINEER 2001中文版入门与提高[M]．北京：清华大学出版社，2003：50-202
[48] 何玉安．楔横轧机模具CAD[J]．模具工业，1998(11)：13-16
[49] 冯建春，张巍，胡正寰．楔横轧模具智能CAD的基本构想[J]．塑性工程学报，1995，2(1)：61-65
[50] 杜凤山，刘庆玲．基于特征建模技术的楔横轧模具CAD系统研究[J]．燕山大学学报，2003，27(1)：8-13
[51] 应富强，张更超，潘孝勇．三维有限元模拟技术在金属塑性成形中的应用[J]．锻压装备与制造技术，2003(5)：10-13
[52] 陈军，王玉国，王跃先，董万鹏，阮雪榆．金属塑性成形有限元模拟中的六面体单元再划分技术研究进展[J]．塑性工程学报，2002，9(3)：1-4
[53] 刘建生，陈慧琴，郭晓霞．金属塑性加工有限元模拟技术与应用[M]．北京：冶金工业出版社，2003：1-12
[54] 刘建生，王立元，原向阳．金属塑性成形的三维刚塑性有限元模拟技术研究[J]．太原重型机械学院学报，1999：20(4)：283-288
[55] 杨伟军，包忠诩，柳和生，江雄心．金属塑性成形有限元分析中的网格生成与重划[J]．金属成形工艺，1998，16(4)：38-40
[56] 张建华，叶尚辉．有限元网格自动生成典型方法及发展方向[J]．计算机辅助设计与制造，1996，2：28-31
[57] 陈如欣，胡忠民．塑性有限元法及其在金属成形中的应用[M]．重庆：重庆大学出版社，1989：180-182
[58] 杨翠苹．楔横轧展宽过程的三维有限元分析[D]．北京科技大学硕士学位论文，2001
[59] Zb. Pater, W. Weroski, J. Kazanecki and A. Gontarz. Study of the process stability of cross wedge rolling[J]. Journal of Materials Processing Technology,
1999(92-93)：458-462
[60] 汪建敏．二辊楔横轧工艺参数的研究[J]．模具工业，1997(12)：35-40
[61] 梁继才，付沛富，李义．用楔横轧工艺成形阶梯轴类件的防止拉细和加工界限研究[J]．吉林工业大学自然科学学报，2000，30(4)：18-20
[62] 徐春国，任广升．楔横轧工艺两次轧制时断面缩减率的分配[J]．塑性工程学报，1999，6(4)：102-105
[63] 王全先．板式楔横轧力能参数的理论计算[J]．锻压技术，1994，2：32-35
[64] 张庆勇，白志斌，刘化民，任广升．楔横轧变形载荷的上限分析[J]．吉林工业大学学报，1995，1：53-61
[65] 束学道，胡正寰．大型楔横轧机轧制力和轧制力矩有限元算法[J]．轧钢，2003，20(6)：4-6
[66] 王红梅．楔橫轧展宽过程三维热—力耦合有限元仿真[D]．北京科技大学硕士学位论文，2002：48-52
[2] 王宝雨，胡正寰．板式楔横轧轧件成形过程几何模拟[J]．钢铁1999(10)：27-31
[3] 宋玉泉．连续局部塑性成形的发展前景[J]．中国机械工程，2000，11(1-2)：65-67
[4] Yaomion Dong, Kaveh A. Tagavi and Michael R. Lovell.Analysis of interracial slip in cross-wedge rolling: a numerical and phenomenological investigation[J].Journal of Materials Processing Technology, 2000, 97(1-3):44-53
[5] 胡正寰，张康生，王宝雨．楔横轧理论与应用[M]．北京：冶金工业出版社，1996：13-25
[6] 胡正寰，许协和，沙德元．斜轧与楔横轧—原理、工艺及设备[M]．北京：冶金工业出版社，1981：14-17
[7] 曲建国．小变形楔横轧件加工可能性及其质量控制的研究[D]．长春：吉林工业大学经济管理学院，1984
[8] 徐春国，王景梁，任广升．楔横轧成形直角台阶精确截齐曲线的计算[J]．塑性工程学报，1999，6(2)：93-98
[9] 张康生，胡正寰．楔横轧精确直角台阶轧齐曲线[J]．锻压技术，1996(6)：27-33
[10] 张康生，胡正寰．楔横轧轧齐基本理论研究[J]．锻压技术，1996(4)：30-34
[11] 张康生，胡正寰．直角密集台阶轧齐曲线[J]．塑性工程学报，1996，3(1)：44-49
[12] 束学道，张康生，胡正寰，杨翠苹．楔横轧力能参数影响因素分析[J]．重型机械，2002(4)：29-33
[13] Yaomin Dong, Michael Lovell and Kaveh Tagavi .Analysis of interfacial slip in cross-wedge rolling: an experimentally verified finite-element model[J]. Journal of Materials Processing Technology, 1998(80-81):273-281
[14] Qiang,Li,Michael R.Lovell,William Slaughter,Kaveh Tagavi. Investigation of
the morphology of internal defects in cross wedge rolling[J]. Journal of Materials Processing Technology, 2002(125-126):248-257
[15] 马振海，杨翠苹，胡正寰．楔横轧轧件端头凹心影响因素地研究[J]．锻压技术，2002(1)：29-31
[16] 梁继才，任广升等．空心件楔横轧旋转条件的分析[J]．农业机械学报，1995，26(1)：97-101
[17] 梁继才，任广升等．空心件楔横轧压扁失稳的实验研究[J]．吉林工业大学学报，1994，24(4)：78-81
[18] 崔波，梁继才，张昕，刘国平．空心件楔横轧成型过程的应力应变分析[J]．农业机械学报，2001，32(1)：102-104
[19] 颜世公，李玉京．冷楔横轧的实验研究[J]．锻压技术，1997，5：31-34
[20] 权修华，董定福．冷楔横轧工艺实验与研究[J]．锻压机械，1996(4)：24-26
[21] 董定福，权修华，曹诗倬．冷楔横轧中几个工艺问题的探讨[J]．锻压技术，1996(6)：34-36
[22] 白志斌，李东平等．汽车刹车蹄片轴的偏心精密成形[J]．吉林工业大学自然科学学报，2000，30(4)：21-23
[23] 任广升．小型拖拉机阶梯轴类锻坯楔横轧成形[J]．农业机械学报，1993，24(2)：77-80
[24] 吕增亭．楔横轧轴类零件的矫正[J]．机械工程师，1998，12：23-23
[25] 卫原平，彭颖红，阮雪榆．金属塑性成形过程的计算机模拟系统[J]．上海交通大学学报，1996，30(3)：132-137
[26] 洪深泽，胡积威，陈建学．金属成形过程的刚塑性有限元数值模拟[J]．合肥工业大学学报(自然科学版)，1995，18(1)：22-27
[27] 周思柱，X．Q．Li．金属成形的三维有限元模拟系统的发展与展望[J]．锻压技术，1998(4)：46-50
[28] 刘光宇，虞跃生，万贤毅．计算机模拟技术在金属塑性成形中的应用[J]．汽车科技，2000(6)：34-37
[29] 鹿守理，赵辉，张鹏．金属塑性加工的计算机模拟[J]．轧钢，1997(4)：54-57
[30] 谢配良，盛伟，游庆．金属塑性成形过程的计算机模拟[J]．锻压机
械．1997(4)：45-48
[31] 单德彬，吕炎，王真．金属体积成形过程三维刚塑性有限元模拟技术的研究[J]．塑性工程学报，1997，4(3)：98-102
[32] Lee N K,Yoon J H and Yang D Y. Finite Element Analysis of Large Deformation by Automatic Renoding as a Weak Remeshing Technique[J]. Int.J.Mech.Sci, 1992, 34(4):255-273
[33] Dyduch M,Habraken A M and Cescotto S.Automatic Adaptive Remeshing for Numerical Simulation of Metal Forming[J]. Computer Methods in Applied Mechanics and Engineering,1992(101):283-298
[34] 王克智，杜庆柏，马世宁．我国塑性加工摩擦学研究展望[J]．中国表面工程，1995(1)：33-37
[35] 张士宏．塑性加工技术的科学化与中国塑性加工技术的发展[J]．世界科技研究与发展，2001(5)：5-11
[36] 王纪武，刘庄，王本一．体积成形的三维有限变形弹塑性有限元法的研究[J]，塑性工程学报，2003，7(1)：41-44
[37] 胡忠．塑性有限元模拟技术的最新进展[J]．塑性工程学报，1994，1(3)：3-12
[38] 寇淑清，博沛富，杨慎华，胡平，李运兴．三维大变形弹塑性有限元分析与边界条件的处理[J]．塑性工程学报，1998，5(4)：65-70
[39] 闫洪，鲍乐，王美艳，包忠诩．弹塑性有限元法在金属塑性加工中的应用[J]．模具技术，2000(5)：12-14
[40] 孟凡中．弹塑性有限变形理论和有限元方法[M]．北京：清华大学出版社，1985：57-211
[41] 王祖城，汪家才．弹性和塑性理论及有限单元法[M]．北京：冶金工业出版社，1983：1-54
[42] 刘涛，杨凤鹏．精通ANSYS[M]．北京：清华大学出版社，2002：1-170
[43] 博嘉科技．ANSYS融会与贯通[M]．北京：中国水利水电出版社，2002：126-165
[44] 马振海，胡正寰，杨翠萍，束学道．楔横轧展宽段的变形特征与应力应变分析[J]．北京科技大学学报，2002，24(3)：309-312
[45] Yaomin Dong, Kaveh A. Tagavi, Michael R. Lovell and Zhi Deng. Analysis of stress in cross wedge rolling with application to failure[J]. International Journal of Mechanical Sciences, 2000(42): 1233-1253
[46] 黄忠耀，李冬梅，王景先．Pro/ENGINEER 2001基础训练教程[M]．北京：清华大学出版社，2002：368-390
[47] 孙江宏，黄小龙．Pro/ENGINEER 2001中文版入门与提高[M]．北京：清华大学出版社，2003：50-202
[48] 何玉安．楔横轧机模具CAD[J]．模具工业，1998(11)：13-16
[49] 冯建春，张巍，胡正寰．楔横轧模具智能CAD的基本构想[J]．塑性工程学报，1995，2(1)：61-65
[50] 杜凤山，刘庆玲．基于特征建模技术的楔横轧模具CAD系统研究[J]．燕山大学学报，2003，27(1)：8-13
[51] 应富强，张更超，潘孝勇．三维有限元模拟技术在金属塑性成形中的应用[J]．锻压装备与制造技术，2003(5)：10-13
[52] 陈军，王玉国，王跃先，董万鹏，阮雪榆．金属塑性成形有限元模拟中的六面体单元再划分技术研究进展[J]．塑性工程学报，2002，9(3)：1-4
[53] 刘建生，陈慧琴，郭晓霞．金属塑性加工有限元模拟技术与应用[M]．北京：冶金工业出版社，2003：1-12
[54] 刘建生，王立元，原向阳．金属塑性成形的三维刚塑性有限元模拟技术研究[J]．太原重型机械学院学报，1999：20(4)：283-288
[55] 杨伟军，包忠诩，柳和生，江雄心．金属塑性成形有限元分析中的网格生成与重划[J]．金属成形工艺，1998，16(4)：38-40
[56] 张建华，叶尚辉．有限元网格自动生成典型方法及发展方向[J]．计算机辅助设计与制造，1996，2：28-31
[57] 陈如欣，胡忠民．塑性有限元法及其在金属成形中的应用[M]．重庆：重庆大学出版社，1989：180-182
[58] 杨翠苹．楔横轧展宽过程的三维有限元分析[D]．北京科技大学硕士学位论文，2001
[59] Zb. Pater, W. Weroski, J. Kazanecki and A. Gontarz. Study of the process stability of cross wedge rolling[J]. Journal of Materials Processing Technology,
1999(92-93)：458-462
[60] 汪建敏．二辊楔横轧工艺参数的研究[J]．模具工业，1997(12)：35-40
[61] 梁继才，付沛富，李义．用楔横轧工艺成形阶梯轴类件的防止拉细和加工界限研究[J]．吉林工业大学自然科学学报，2000，30(4)：18-20
[62] 徐春国，任广升．楔横轧工艺两次轧制时断面缩减率的分配[J]．塑性工程学报，1999，6(4)：102-105
[63] 王全先．板式楔横轧力能参数的理论计算[J]．锻压技术，1994，2：32-35
[64] 张庆勇，白志斌，刘化民，任广升．楔横轧变形载荷的上限分析[J]．吉林工业大学学报，1995，1：53-61
[65] 束学道，胡正寰．大型楔横轧机轧制力和轧制力矩有限元算法[J]．轧钢，2003，20(6)：4-6
[66] 王红梅．楔橫轧展宽过程三维热—力耦合有限元仿真[D]．北京科技大学硕士学位论文，2002：48-52