提高少片簧承载能力的设计方法研究
|
摘要
少片簧通常用于轻微型车的悬架系统中，是汽车中重要的高负荷安全部件，主要用来传递垂直力，并缓和由不平路面引起的冲击和振动，作为传递力和力矩连接装置，传递性能的好坏直接影响到操纵稳定性、制动安全性、驾驶舒适性和动力加速性等汽车动力学性能。只有配备性能优良的少片簧，才会得到性能优良的汽车。在使用钢板弹簧之前，计算它的刚度和应力是非常重要的。通常可以使用共同曲率法、集中载荷法和有限元法。
首先，本文综述了少片钢板弹簧的研究现状；采用结构力学方法，对不同的少片钢板弹簧的设计理念进行介绍，并推导相应力学参数的计算公式；并在哈密顿辛对偶变量体系下，利用精细积分法和正交试验法，分析片数、弧高、截面尺寸对刚度的影响，并选出了优水平、主次因素和最优搭配。其中正交试验中所用的数据由精细积分法算出的；
其次利用界带法分析接触对少片簧刚度的影响；基于界带理论，把接触赋予厚度，接触对板簧刚度的影响看成是对于板簧刚度的摄动，是一个小量，由共同曲率法知道，接触和板簧在片端的挠度一样，看成是和板簧的并联，接触产生的刚度是钢板弹簧刚度的摄动小量，板簧系统的柔度估计偏大，计算中引入了一个负的量把估计偏大的部分给减掉后。分析了接触长度对少片簧刚度的影响，接触长度增加时，少片簧的刚度变化不明显。接触厚度不同时，接触对少片簧的影响是随着接触厚度的增大，板簧的刚度降低。不同激励频率下接触对少片簧刚度的影响是，激励频率从10-4Hz增加到2Hz时，接触对少片簧的刚度影响不明显;当从两赫兹增加到十的四次方赫兹时，接触随着激励频率的增加而增加。
再次，利用有限元分析软件ABAQUS对少片钢板弹簧进行有限元分析，计算了少片钢板弹簧的刚度、应力特性曲线。用PRO/E，对板簧建模，用HYPEMESH对其进行了前处理，用ABAQUS对变截面少片钢板弹簧刚度应力进行分析。
最后，利用正交实验法分析,在7因素2水平情况交互作用和无交互作用条件下，分析不同长度、厚度组合等不同截面尺寸情况下的刚度特性。在正交试验中，利用ABAQUS软件对不同情况进行计算。利用ABAQUS软件计算了21种交互作用后发现无交互作用的不是最优搭配。
Less leaf spring is usually used for minor car's suspension system and is an importantautomotive safety components under the high load conditions. It's main role is to pass thevertical force, to ease the shock and vibration caused by the uneven road, as connecteddevices for the transmission of the force and torque, and the transfer properties is directlyimpact on the vehicle dynamics performance such an the handling stability, the brakingsafety, the driving comfort and the driving force of acceleration. If we want to get a goodperformance car, only with the excellent performance of the small leaf spring. Before theproduction of leaf springs, it is very important to calculate its stiffness and stress. We canusually use a common curvature, concentrated load method and finite element method.
First, this paper reviews the status of a small piece of leaf spring; using the method ofstructural mechanics to introduce small piece of leaf spring design concepts and derive theformula of the corresponding mechanical parameters; under the Hamiltonian symplecticdual variables system using the precise integration method and orthogonal test to analyzethe impact of the number of pieces, the arc height and the section size for the stiffness andelect a superior level of primary and secondary factors and the optimal mix. And the data ofthe orthogonal experiment is taken by the precise integration method.
Then using the inter-belt method analyzes the exposure for the less leaf springstiffness's impact; basing on the inter-belt theory, giving the contact thickness, consideringthe contact's impact of leaf spring stiffness as a leaf spring stiffness perturbation, which is asmall amount and known by the common curvature method, contact and the leaf springdeflection in the film side is the same and consider it parallel to the leaf springs, thestiffness caused by the contact is a small amount perturbation for the spring stiffness, theestimate of the flexibility of the leaf spring system is too large, so while calculationintroducing a negative amount to lose the estimated large part. Analyzing the impact of thecontact length on the leaf spring stiffness finds that while the contact length increases, theleaf spring stiffness did not change significantly. While the contact thickness is different, the thicker the contact thickness the smaller the leaf spring stiffness. Under different excitationfrequency the impact of the contact of the leaf spring stiffness is that while the change ofexcitation frequency is from10-4to2Hz, the contact has no obvious influence for thespring stiffness; but while the excitation frequency is from2to10Hz, contact increaseswith the excitation frequency increasing..
Again using the finite element analysis software ABAQUS to do the finite elementanalysis, calculating the stiffness of the leaf spring and drawing the stress curve. The paperuses the PRO/E software to model the leaf spring, the HYPERMESH software to do thepre-treatment and uses the ABAQUA software to analyze the impact of the variablecross-section of the stiffness of the leaf spring.
Finally under the seven factors and two levels of interaction and no interactionconditions using the orthogonal experiment analyze stiffness characteristics under differentsectional dimensions of length, thickness combination. In orthogonal experiment using theABAQUS software calculates the different situations. While using the ABAQUS softwareto calculate the21kinds of interaction，it is found that the no interaction is not the optimalmatch.
首先，本文综述了少片钢板弹簧的研究现状；采用结构力学方法，对不同的少片钢板弹簧的设计理念进行介绍，并推导相应力学参数的计算公式；并在哈密顿辛对偶变量体系下，利用精细积分法和正交试验法，分析片数、弧高、截面尺寸对刚度的影响，并选出了优水平、主次因素和最优搭配。其中正交试验中所用的数据由精细积分法算出的；
其次利用界带法分析接触对少片簧刚度的影响；基于界带理论，把接触赋予厚度，接触对板簧刚度的影响看成是对于板簧刚度的摄动，是一个小量，由共同曲率法知道，接触和板簧在片端的挠度一样，看成是和板簧的并联，接触产生的刚度是钢板弹簧刚度的摄动小量，板簧系统的柔度估计偏大，计算中引入了一个负的量把估计偏大的部分给减掉后。分析了接触长度对少片簧刚度的影响，接触长度增加时，少片簧的刚度变化不明显。接触厚度不同时，接触对少片簧的影响是随着接触厚度的增大，板簧的刚度降低。不同激励频率下接触对少片簧刚度的影响是，激励频率从10-4Hz增加到2Hz时，接触对少片簧的刚度影响不明显;当从两赫兹增加到十的四次方赫兹时，接触随着激励频率的增加而增加。
再次，利用有限元分析软件ABAQUS对少片钢板弹簧进行有限元分析，计算了少片钢板弹簧的刚度、应力特性曲线。用PRO/E，对板簧建模，用HYPEMESH对其进行了前处理，用ABAQUS对变截面少片钢板弹簧刚度应力进行分析。
最后，利用正交实验法分析,在7因素2水平情况交互作用和无交互作用条件下，分析不同长度、厚度组合等不同截面尺寸情况下的刚度特性。在正交试验中，利用ABAQUS软件对不同情况进行计算。利用ABAQUS软件计算了21种交互作用后发现无交互作用的不是最优搭配。
Less leaf spring is usually used for minor car's suspension system and is an importantautomotive safety components under the high load conditions. It's main role is to pass thevertical force, to ease the shock and vibration caused by the uneven road, as connecteddevices for the transmission of the force and torque, and the transfer properties is directlyimpact on the vehicle dynamics performance such an the handling stability, the brakingsafety, the driving comfort and the driving force of acceleration. If we want to get a goodperformance car, only with the excellent performance of the small leaf spring. Before theproduction of leaf springs, it is very important to calculate its stiffness and stress. We canusually use a common curvature, concentrated load method and finite element method.
First, this paper reviews the status of a small piece of leaf spring; using the method ofstructural mechanics to introduce small piece of leaf spring design concepts and derive theformula of the corresponding mechanical parameters; under the Hamiltonian symplecticdual variables system using the precise integration method and orthogonal test to analyzethe impact of the number of pieces, the arc height and the section size for the stiffness andelect a superior level of primary and secondary factors and the optimal mix. And the data ofthe orthogonal experiment is taken by the precise integration method.
Then using the inter-belt method analyzes the exposure for the less leaf springstiffness's impact; basing on the inter-belt theory, giving the contact thickness, consideringthe contact's impact of leaf spring stiffness as a leaf spring stiffness perturbation, which is asmall amount and known by the common curvature method, contact and the leaf springdeflection in the film side is the same and consider it parallel to the leaf springs, thestiffness caused by the contact is a small amount perturbation for the spring stiffness, theestimate of the flexibility of the leaf spring system is too large, so while calculationintroducing a negative amount to lose the estimated large part. Analyzing the impact of thecontact length on the leaf spring stiffness finds that while the contact length increases, theleaf spring stiffness did not change significantly. While the contact thickness is different, the thicker the contact thickness the smaller the leaf spring stiffness. Under different excitationfrequency the impact of the contact of the leaf spring stiffness is that while the change ofexcitation frequency is from10-4to2Hz, the contact has no obvious influence for thespring stiffness; but while the excitation frequency is from2to10Hz, contact increaseswith the excitation frequency increasing..
Again using the finite element analysis software ABAQUS to do the finite elementanalysis, calculating the stiffness of the leaf spring and drawing the stress curve. The paperuses the PRO/E software to model the leaf spring, the HYPERMESH software to do thepre-treatment and uses the ABAQUA software to analyze the impact of the variablecross-section of the stiffness of the leaf spring.
Finally under the seven factors and two levels of interaction and no interactionconditions using the orthogonal experiment analyze stiffness characteristics under differentsectional dimensions of length, thickness combination. In orthogonal experiment using theABAQUS software calculates the different situations. While using the ABAQUS softwareto calculate the21kinds of interaction，it is found that the no interaction is not the optimalmatch.
引文
[1][苏]А.М.Горелик.蒋立盛译.少片簧的研究[J].苏.АП.1981年第1期:17-23
[2]马文伦.重型载货车少片变截面钢板弹簧的研究[D].哈尔滨工程大学工学硕士论文.2009:1-7
[3]郑银环.汽车钢板弹簧计算模型研究[D].武汉理工大学工学博士论文.2005:5-9
[4]张宁.汽车钢板弹簧设计[D].吉林大学工学硕士论文.2007:30-40
[5]步一鸣.车用少片变截面弹簧.当代汽车，1989，（3）：16-25
[6]李建辉.钢板弹簧迟滞特性非线性有限元分析[D].吉林大学本科生论文.2010
[7]刘帅.钢板弹簧理论设计方法及可视化[D].吉林大学本科生论文.2010
[8]陈家瑞.汽车构造（下册）[M].人民交通出版社,1976
[9]王望予.汽车设计[M].机械工业出版社,2000
[10]张洪欣.汽车设计[M].机械工业出版社,1996
[11]刘维信.汽车设计[M].清华大学出版社,2001
[12]小田柿浩三.汽车设计[M].机械工业出版社,1978
[13]肖军.汽车钢板弹簧的应用及其发展趋势[J].城市车辆,2007
[14]刘凤景.多片钢板弹簧参数化有限元分析与优化系统研究[D].山东科技大学硕士论文.2009:1-12
[15]李未.复合材料板簧的优化设计及有限元分析[D].哈尔滨工业大学硕士论文.2006:1-18
[17]凌荣江.汽车板簧的建模及非线性分析[D].南京航空航天大学硕士论文.2008:1-15
[18]房月根.高速叉车悬架系统研究及分析[D].长安大学硕士论文.2008:1-21
[19]李轶石.基于有限元的某多片钢板弹簧性能仿真研究[D].湖南大学硕士论文2007:1-8
[20]石湘.渐变刚度钢板弹簧的模拟仿真[D].华中科技大学硕士论文.2007:1-9
[21]雷昌浩.汽车钢板弹簧计算的一种新模型及其静动态特性研究[D].广西工学院硕士论文.2010:5-8
[22]何彬.汽车钢板弹簧自组织配置设计方法的研究[D].武汉理工大学博士论文.2009:4-8
[23]蒋立盛，方坚.少片弹簧的计算及其尺寸参数的选择[J].汽车技术.1984:23-25
[24]张俊海等.车辆发动机气门弹簧的设计方法探讨[J].拖拉机与农用运输车.2007:10-15
[25]周站福等.少片钢板弹簧在实际应用中的开发与设计[J].汽车技术.2009:27-32
[26]秦东晨等.汽车钢板弹簧多体模型建立的一种方法[J].武汉理工大学学报.2007:19-24
[27] H.P.帕尔希洛夫斯基.汽车的钢板弹簧[M].人民交通出版社,1959
[28] H.H.乌斯潘斯基.汽车悬架设计[M].人民交通出版社,1980
[29]杨宗孟.汽车钢板弹簧的设计与计算[J].汽车技术,1979:12-28
[30]郭孔辉.钢板弹簧刚度计算的主片分析法[J].汽车工程,1984:22-28
[31]林金木.有限单元法解钢板弹簧[J].汽车工程,1984:20-36
[32]林金木.应用余能原理求解渐变刚度钢板弹簧[J].汽车工程,1996:301-308
[33]林金木.渐变刚度钢板弹簧的有限元计算[J].农业机械学报,1996:155-160
[34]林金木，柏伟武.钢板弹簧设计计算[J].汽车工程,1995:169-177
[35]林国成，王黎钦，秦满.变厚截面板簧的有限元分析[J].弹簧工程,1996:7-9
[36]李金伯，周中坚.车辆前钢板弹簧的有限元计算模拟[J].上海汽车,1998:29-30
[37]姚层林.汽车钢板弹簧的有限元分析及程序化实现[J].汽车科技,1999:32-33
[38]刘巧伶林逸，张义民等.汽车多片钢板弹簧的可靠性设计[J].中国公路学报,1998:109-112
[39]王其东等.基于多体理论的汽车悬架系统分析、设计与控制研究[J].合肥工业大学学报（自然科学版）：1999,22（6）：11-15
[40]胡玉梅，邓兆祥，王欣，宾洋.汽车后悬架的非线性有限元分析[J].重庆大学学报(自然科学版),2003:38-41
[41]丁能根等.钢板弹簧迟滞特性的有限元分析[J].汽车工程2003,25(1):12-14
[42] Zahavi E. Analysis of a contact problem in leaf springs[J].MechanicsResearch Communications,1992,19(1):21-27
[43] Zahavi, E.Contact. Problems with friction in machine design[J].ComputersStructures,1993,4(17):591-594
[44] Li Q.A contact finite element algorinthm for the multi-leaf spring ofvehicle suspension systems[J].Journal of AutomobileEngineering,2004,218(3):305-314
[45]李建辉等.多片钢板弹簧减摩措施研究[J].哈尔滨工业大学学报，2011:35-37
[46]姚征.辛体系算法在波的传播与振动问题中的应用[D].大连理工大学博士论文.2007:14-25
[47]韩爱红等.基础隔振结构地震响应时程分析新方法[J].振动与冲击.2007:36-39
[48]胡启平等.铁摩辛柯梁弯曲问题的精细积分方法[J].工业建筑.2007:35-37
[49]陈畅亚.改性粉煤灰对城市污泥的调理及机理分析[D].湖南大学硕士论文.2009:4-11
[50]测试用例设计技术之五-正交实验法之正交表的由来[J/OL]《网络（http://www.test-edu.net/newsshow-23-246.html）》
[51]正交实验设计-MBA智库百科[J/OL].《网络（http://wiki.mbalib.com/wiki/%e6%ad%a3%e4%ba%a4%e8%af%95%e9%aa%8c%e8%ae%be%e8%ae%a1）》
[52]李建文.多因素致某地下结构破坏成因的数值模拟.燕山大学硕士论文.2011:15-28
[53]郭洪强等.汽车钢板弹簧的接触有限元分析[J].湖北汽车工业学院学报2006,20(2):11-3
[54]于安和等.钢板弹簧刚度特性的有限元分析[J].汽车技术2007(2):26-29
[55]丁华等.钢板弹簧刚度计算影响因数分析[J].拖拉机与农用运输车,2008,35(6):56-57
[56]郑银环等.钢板弹簧力学特性的非线性有限元分析[J].机械科学与技术,2008,27(7):914-916
[57]王其东等.汽车钢板弹簧多体动力学建模及动特性仿真研究[J].合肥工业大学学报（自然科学版）,1999,22(6):35-37
[58]郑银环等.汽车钢板弹簧多柔体建模及仿真研究[J].湖北工业大学学报,2007,22(4):35-37
[59]雷昌浩等.汽车多片钢板弹簧刚度分析的一种新模型及计算方法[J].拖拉机与农用运输车2010,37(1):52-54
[60]徐劲力.汽车钢板弹簧设计计算分析及软件开发[J].湖北工学院学报2003,18(5):27-29
[61]李未.复合材料板簧的优化设计及有限元分析哈尔滨[D].哈尔滨工业大学出版社
[62]佟刚等.利用Matlab优化设计少片变截面钢板弹簧.沈阳航空工业学院学报,2007,17(34):11-13
[63]郑银环等.利ANSY进行变截面板簧的优化分析[J].现代机械,2005(2):30-31
[64]王勖成.有限单元法[M].清华大学出版社,2003
[65]郭小明，赵惠麟.工程结构接触问题的研究及进展[J].东南大学学报,2003,577-582
[66]周仲坚卢耀祖.机械与汽车结构的有限元分析[M].同济大学出版社,1996
[67]王勖成，邵敏.有限单元法基本原理和数值方法[M].清华大学出版社,1998
[68]赵腾伦. ABAQUS6.6在机械工程中的应用[M].水利水电出版社,2007
[69]韩志仁，陶华，黄赟.机械设计中有限元分析的几个关键问题[J].机械设计与制造,2004:58-60
[70]熊珍兵，罗会信.基于HyperMesh的有限元前处理技术[J].排灌机械,2006
[71]于开平，周传月，谭惠丰等. HyperMesh从入门到精通[M].科学出版社,2005
[72]张洪武等.界带分析的基本理论和计算方法[J].计算力学学报.2006:19-24
[73]钟万勰等.界带分析[C].中国力学学会学术大会2005论文摘要集.20050801
[74]姚伟岸，钟万勰.辛弹性力学[D].高等教育出版社.2002.4
[75]钟万勰等.辛破茧辛拓展辛层次[M].大连理工大学出版社,2011
[76]钟万勰.力功能量与辛数学[M].大连理工大学出版社,2005
[77]钟万勰.应用力学对偶体系[M].科学出版社,2002
[78]任露泉.实验优化技术[M].机械工业出版,1987
[79]钟万勰.应用力学的辛数学方法[M].高等教育出版社,2006
[80]冯康，秦孟兆.Hamilton体系的辛计算格式[M].浙江科技出版社,2004
[81]高强，钟万勰. Hamilton系统的保辛-守恒算法[J].动力学与控制学报,2009,7（3）：193-197
[82]钟万勰，张洪武，武承伟.参变量变分原理及其应用[M].科学出版社,1997
[83]张洪武.参变量变分原理与材料和结构力学分析[M].科学出版社,2010
[84]钟万勰，姚征.椭圆函数的精细积分算法[J].应用力学进展,科学出版社，2004:106-111
[85]张洪武，姚征，钟万勰.界带分析的基本理论和计算方法[J].计算力学学报,2006，23:257-263
[86]姚征，张洪武，王晋宝，钟万勰.基于界带模型的炭纳米管声子谱的辛分析.固体力学学报[J].2008,29（1）：13-22
[87] И.Пархиловский.Автомобильные листовые рессоры [M]. Машиностроение,1978
[88]张启江.变截面钢板弹簧梯形叶片的斜面位置分析[J].汽车技术，2004（2）:15-17
[2]马文伦.重型载货车少片变截面钢板弹簧的研究[D].哈尔滨工程大学工学硕士论文.2009:1-7
[3]郑银环.汽车钢板弹簧计算模型研究[D].武汉理工大学工学博士论文.2005:5-9
[4]张宁.汽车钢板弹簧设计[D].吉林大学工学硕士论文.2007:30-40
[5]步一鸣.车用少片变截面弹簧.当代汽车，1989，（3）：16-25
[6]李建辉.钢板弹簧迟滞特性非线性有限元分析[D].吉林大学本科生论文.2010
[7]刘帅.钢板弹簧理论设计方法及可视化[D].吉林大学本科生论文.2010
[8]陈家瑞.汽车构造（下册）[M].人民交通出版社,1976
[9]王望予.汽车设计[M].机械工业出版社,2000
[10]张洪欣.汽车设计[M].机械工业出版社,1996
[11]刘维信.汽车设计[M].清华大学出版社,2001
[12]小田柿浩三.汽车设计[M].机械工业出版社,1978
[13]肖军.汽车钢板弹簧的应用及其发展趋势[J].城市车辆,2007
[14]刘凤景.多片钢板弹簧参数化有限元分析与优化系统研究[D].山东科技大学硕士论文.2009:1-12
[15]李未.复合材料板簧的优化设计及有限元分析[D].哈尔滨工业大学硕士论文.2006:1-18
[17]凌荣江.汽车板簧的建模及非线性分析[D].南京航空航天大学硕士论文.2008:1-15
[18]房月根.高速叉车悬架系统研究及分析[D].长安大学硕士论文.2008:1-21
[19]李轶石.基于有限元的某多片钢板弹簧性能仿真研究[D].湖南大学硕士论文2007:1-8
[20]石湘.渐变刚度钢板弹簧的模拟仿真[D].华中科技大学硕士论文.2007:1-9
[21]雷昌浩.汽车钢板弹簧计算的一种新模型及其静动态特性研究[D].广西工学院硕士论文.2010:5-8
[22]何彬.汽车钢板弹簧自组织配置设计方法的研究[D].武汉理工大学博士论文.2009:4-8
[23]蒋立盛，方坚.少片弹簧的计算及其尺寸参数的选择[J].汽车技术.1984:23-25
[24]张俊海等.车辆发动机气门弹簧的设计方法探讨[J].拖拉机与农用运输车.2007:10-15
[25]周站福等.少片钢板弹簧在实际应用中的开发与设计[J].汽车技术.2009:27-32
[26]秦东晨等.汽车钢板弹簧多体模型建立的一种方法[J].武汉理工大学学报.2007:19-24
[27] H.P.帕尔希洛夫斯基.汽车的钢板弹簧[M].人民交通出版社,1959
[28] H.H.乌斯潘斯基.汽车悬架设计[M].人民交通出版社,1980
[29]杨宗孟.汽车钢板弹簧的设计与计算[J].汽车技术,1979:12-28
[30]郭孔辉.钢板弹簧刚度计算的主片分析法[J].汽车工程,1984:22-28
[31]林金木.有限单元法解钢板弹簧[J].汽车工程,1984:20-36
[32]林金木.应用余能原理求解渐变刚度钢板弹簧[J].汽车工程,1996:301-308
[33]林金木.渐变刚度钢板弹簧的有限元计算[J].农业机械学报,1996:155-160
[34]林金木，柏伟武.钢板弹簧设计计算[J].汽车工程,1995:169-177
[35]林国成，王黎钦，秦满.变厚截面板簧的有限元分析[J].弹簧工程,1996:7-9
[36]李金伯，周中坚.车辆前钢板弹簧的有限元计算模拟[J].上海汽车,1998:29-30
[37]姚层林.汽车钢板弹簧的有限元分析及程序化实现[J].汽车科技,1999:32-33
[38]刘巧伶林逸，张义民等.汽车多片钢板弹簧的可靠性设计[J].中国公路学报,1998:109-112
[39]王其东等.基于多体理论的汽车悬架系统分析、设计与控制研究[J].合肥工业大学学报（自然科学版）：1999,22（6）：11-15
[40]胡玉梅，邓兆祥，王欣，宾洋.汽车后悬架的非线性有限元分析[J].重庆大学学报(自然科学版),2003:38-41
[41]丁能根等.钢板弹簧迟滞特性的有限元分析[J].汽车工程2003,25(1):12-14
[42] Zahavi E. Analysis of a contact problem in leaf springs[J].MechanicsResearch Communications,1992,19(1):21-27
[43] Zahavi, E.Contact. Problems with friction in machine design[J].ComputersStructures,1993,4(17):591-594
[44] Li Q.A contact finite element algorinthm for the multi-leaf spring ofvehicle suspension systems[J].Journal of AutomobileEngineering,2004,218(3):305-314
[45]李建辉等.多片钢板弹簧减摩措施研究[J].哈尔滨工业大学学报，2011:35-37
[46]姚征.辛体系算法在波的传播与振动问题中的应用[D].大连理工大学博士论文.2007:14-25
[47]韩爱红等.基础隔振结构地震响应时程分析新方法[J].振动与冲击.2007:36-39
[48]胡启平等.铁摩辛柯梁弯曲问题的精细积分方法[J].工业建筑.2007:35-37
[49]陈畅亚.改性粉煤灰对城市污泥的调理及机理分析[D].湖南大学硕士论文.2009:4-11
[50]测试用例设计技术之五-正交实验法之正交表的由来[J/OL]《网络（http://www.test-edu.net/newsshow-23-246.html）》
[51]正交实验设计-MBA智库百科[J/OL].《网络（http://wiki.mbalib.com/wiki/%e6%ad%a3%e4%ba%a4%e8%af%95%e9%aa%8c%e8%ae%be%e8%ae%a1）》
[52]李建文.多因素致某地下结构破坏成因的数值模拟.燕山大学硕士论文.2011:15-28
[53]郭洪强等.汽车钢板弹簧的接触有限元分析[J].湖北汽车工业学院学报2006,20(2):11-3
[54]于安和等.钢板弹簧刚度特性的有限元分析[J].汽车技术2007(2):26-29
[55]丁华等.钢板弹簧刚度计算影响因数分析[J].拖拉机与农用运输车,2008,35(6):56-57
[56]郑银环等.钢板弹簧力学特性的非线性有限元分析[J].机械科学与技术,2008,27(7):914-916
[57]王其东等.汽车钢板弹簧多体动力学建模及动特性仿真研究[J].合肥工业大学学报（自然科学版）,1999,22(6):35-37
[58]郑银环等.汽车钢板弹簧多柔体建模及仿真研究[J].湖北工业大学学报,2007,22(4):35-37
[59]雷昌浩等.汽车多片钢板弹簧刚度分析的一种新模型及计算方法[J].拖拉机与农用运输车2010,37(1):52-54
[60]徐劲力.汽车钢板弹簧设计计算分析及软件开发[J].湖北工学院学报2003,18(5):27-29
[61]李未.复合材料板簧的优化设计及有限元分析哈尔滨[D].哈尔滨工业大学出版社
[62]佟刚等.利用Matlab优化设计少片变截面钢板弹簧.沈阳航空工业学院学报,2007,17(34):11-13
[63]郑银环等.利ANSY进行变截面板簧的优化分析[J].现代机械,2005(2):30-31
[64]王勖成.有限单元法[M].清华大学出版社,2003
[65]郭小明，赵惠麟.工程结构接触问题的研究及进展[J].东南大学学报,2003,577-582
[66]周仲坚卢耀祖.机械与汽车结构的有限元分析[M].同济大学出版社,1996
[67]王勖成，邵敏.有限单元法基本原理和数值方法[M].清华大学出版社,1998
[68]赵腾伦. ABAQUS6.6在机械工程中的应用[M].水利水电出版社,2007
[69]韩志仁，陶华，黄赟.机械设计中有限元分析的几个关键问题[J].机械设计与制造,2004:58-60
[70]熊珍兵，罗会信.基于HyperMesh的有限元前处理技术[J].排灌机械,2006
[71]于开平，周传月，谭惠丰等. HyperMesh从入门到精通[M].科学出版社,2005
[72]张洪武等.界带分析的基本理论和计算方法[J].计算力学学报.2006:19-24
[73]钟万勰等.界带分析[C].中国力学学会学术大会2005论文摘要集.20050801
[74]姚伟岸，钟万勰.辛弹性力学[D].高等教育出版社.2002.4
[75]钟万勰等.辛破茧辛拓展辛层次[M].大连理工大学出版社,2011
[76]钟万勰.力功能量与辛数学[M].大连理工大学出版社,2005
[77]钟万勰.应用力学对偶体系[M].科学出版社,2002
[78]任露泉.实验优化技术[M].机械工业出版,1987
[79]钟万勰.应用力学的辛数学方法[M].高等教育出版社,2006
[80]冯康，秦孟兆.Hamilton体系的辛计算格式[M].浙江科技出版社,2004
[81]高强，钟万勰. Hamilton系统的保辛-守恒算法[J].动力学与控制学报,2009,7（3）：193-197
[82]钟万勰，张洪武，武承伟.参变量变分原理及其应用[M].科学出版社,1997
[83]张洪武.参变量变分原理与材料和结构力学分析[M].科学出版社,2010
[84]钟万勰，姚征.椭圆函数的精细积分算法[J].应用力学进展,科学出版社，2004:106-111
[85]张洪武，姚征，钟万勰.界带分析的基本理论和计算方法[J].计算力学学报,2006，23:257-263
[86]姚征，张洪武，王晋宝，钟万勰.基于界带模型的炭纳米管声子谱的辛分析.固体力学学报[J].2008,29（1）：13-22
[87] И.Пархиловский.Автомобильные листовые рессоры [M]. Машиностроение,1978
[88]张启江.变截面钢板弹簧梯形叶片的斜面位置分析[J].汽车技术，2004（2）:15-17