CN106812851B - 非等偏频型三级渐变刚度板簧最大限位挠度的仿真验算法 - Google Patents

非等偏频型三级渐变刚度板簧最大限位挠度的仿真验算法 Download PDF

Info

Publication number
CN106812851B
CN106812851B CN201710023283.5A CN201710023283A CN106812851B CN 106812851 B CN106812851 B CN 106812851B CN 201710023283 A CN201710023283 A CN 201710023283A CN 106812851 B CN106812851 B CN 106812851B
Authority
CN
China
Prior art keywords
spring
level
simulation calculation
leaf spring
main
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201710023283.5A
Other languages
English (en)
Other versions
CN106812851A (zh
Inventor
周长城
于曰伟
朱召辉
汪晓
杨腾飞
王凤娟
邵明磊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong University of Technology
Original Assignee
Shandong University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shandong University of Technology filed Critical Shandong University of Technology
Priority to CN201710023283.5A priority Critical patent/CN106812851B/zh
Publication of CN106812851A publication Critical patent/CN106812851A/zh
Application granted granted Critical
Publication of CN106812851B publication Critical patent/CN106812851B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F16—ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
    • F16F—SPRINGS; SHOCK-ABSORBERS; MEANS FOR DAMPING VIBRATION
    • F16F3/00—Spring units consisting of several springs, e.g. for obtaining a desired spring characteristic
    • F16F3/02—Spring units consisting of several springs, e.g. for obtaining a desired spring characteristic with springs made of steel or of other material having low internal friction
    • F16F3/023—Spring units consisting of several springs, e.g. for obtaining a desired spring characteristic with springs made of steel or of other material having low internal friction composed only of leaf springs
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F16—ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
    • F16F—SPRINGS; SHOCK-ABSORBERS; MEANS FOR DAMPING VIBRATION
    • F16F1/00—Springs
    • F16F1/02—Springs made of steel or other material having low internal friction; Wound, torsion, leaf, cup, ring or the like springs, the material of the spring not being relevant
    • F16F1/18—Leaf springs
    • F16F1/185—Leaf springs characterised by shape or design of individual leaves
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F16—ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
    • F16F—SPRINGS; SHOCK-ABSORBERS; MEANS FOR DAMPING VIBRATION
    • F16F1/00—Springs
    • F16F1/02—Springs made of steel or other material having low internal friction; Wound, torsion, leaf, cup, ring or the like springs, the material of the spring not being relevant
    • F16F1/18—Leaf springs
    • F16F1/26—Attachments or mountings
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F16—ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
    • F16F—SPRINGS; SHOCK-ABSORBERS; MEANS FOR DAMPING VIBRATION
    • F16F2238/00—Type of springs or dampers
    • F16F2238/02—Springs
    • F16F2238/022—Springs leaf-like, e.g. of thin, planar-like metal
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06—Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Vehicle Body Suspensions (AREA)
  • Springs (AREA)

Abstract

本发明涉及非等偏频型三级渐变刚度板簧最大限位挠度的仿真验算法，属于车辆悬架钢板弹簧技术领域。本发明可根据各片主簧和各级副簧的结构参数，弹性模量，最大许用应力，骑马螺栓夹紧距，主簧及其与各级副簧的夹紧刚度，初始切线弧高，在最大许用载荷和接触载荷及三级渐变夹紧刚度仿真计算的基础上，对非等偏频型三级渐变刚度板簧的最大限位挠度进行仿真验算。通过样机加载应力试验可知，本发明所提供的非等偏频型三级渐变刚度板簧最大限位挠度的仿真验算法是正确的。利用该方法可得到准确可靠的最大限位挠度仿真验算值，提高产品的设计水平、质量、可靠性和使用寿命及车辆行驶安全性；同时，降低设计及试验费用，加快产品开发速度。

Description

非等偏频型三级渐变刚度板簧最大限位挠度的仿真验算法
技术领域
本发明涉及车辆悬架板簧，特别是非等偏频型三级渐变刚度板簧最大限位挠度的仿真验算法。
背景技术
为了满足在不同载荷下的车辆行驶平顺性，可将原一级渐变刚度板簧的主簧和副簧分别拆分为两级，即采用三级渐变刚度板簧；同时，为了满足主簧的应力强度，通常通过主簧和三级副簧初始切线弧高及三级渐变间隙，使三级副簧适当提前承担载荷，从而降低主簧的应力，即采用非等偏频型三级渐变刚度板簧悬架，其中，为了提高非等偏频型三级渐变刚度板簧在冲击载荷下的可靠性及使用寿命，以最大许用应力及最大许用载荷所对应的最大挠度作用板簧的最大限位挠度，并且依据最大限位挠度设置一限位装置，在冲击载荷下对板簧其保护作用，防止板簧因受冲击而断裂，提高板簧的可靠性和使用寿命及车辆行驶安全性。对于给定设计结构的非等偏频型三级渐变刚度板簧，限位装置是否真正在冲击载荷下对板簧起保护作用，必须对最大限位挠度进行仿真验算。然而，由于非等偏频型三级渐变刚度板簧的挠度计算非常复杂，并且受最大许用载荷和接触载荷仿真计算问题的制约，先前国内外一直未给出非等偏频型三级渐变刚度板簧最大限位挠度的仿真验算法，不能满足非等偏频型三级渐变刚度板簧设计及CAD软件开发要求。随着车辆行驶速度及其对平顺性要求的不断提高，对渐变刚度板簧悬架提出了更高要求，因此，必须建立一种精确、可靠的非等偏频型三级渐变刚度板簧最大限位挠度的仿真验算法，为非等偏频型三级渐变刚度板簧的最大限位挠度进行仿真验算提供可靠的技术方法，满足车辆行业快速发展、车辆行驶平顺性及对非等偏频型三级渐变刚度板簧的设计要求，提高产品的设计水平、质量、可靠性和使用寿命及车辆行驶安全性；同时，降低设计及试验费用，加快产品开发速度。
发明内容
针对上述现有技术中存在的缺陷，本发明所要解决的技术问题是提供一种简便、可靠的非等偏频型三级渐变刚度板簧最大限位挠度的仿真验算法，其仿真验算流程如图1所示。三级渐变刚度板簧的一半对称结构如图2所示，是由主簧1、第一级副簧2和第二级副簧3和第三级副簧4所组成的，三级渐变刚度板簧的总跨度的一半等于主簧首片的一半作用长度 L1T，骑马螺栓夹紧距的一半为L0，钢板弹簧的宽度为b，弹性模量为E，最大许用应力[σ]。其中，主簧片数n片，主簧各片的厚度为hi，一半作用长度为LiT，一半夹紧长度 Li＝L1iT-L0/2，i＝1,2,…,n。第一级副簧片数为n1，第一级副簧各片的厚度为hA1j，一半作用长度为LA1jT，一半夹紧长度LA1j＝LA1jT-L0/2，j＝1,2,…,n1。第二级副簧片数为n2，第二级副簧各片的厚度为hA2k，一半作用长度LA2kT，一半夹紧长度LA2k＝LA2kT-L0/2，k＝1,2,…,n2。第三级副簧片数为n3，第三级副簧各片的厚度为hA3l，一半作用长度LA3lT，一半夹紧长度 LA3l＝LA3lT-L0/2，l＝1,2,…,n3。通过主簧和各级副簧的初始切线弧高，在主簧1的末片下表面与第一级副簧2的首片上表面之间设置有第一级渐变间隙δMA1；第一级副簧2的末片下表面与第二级副簧3的首片上表面之间设置有第二级渐变间隙δA12；第二级副簧3的末片下表面与第三级副簧4的首片上表面之间设置有第三级渐变间隙δA23，以满足渐变刚度板簧的接触载荷、渐变刚度、应力强度、悬架偏频及车辆行驶平顺性和安全性的设计要求。依据最大限位挠度所设计的限位装置，在冲击载荷下对板簧起保护作用，提高板簧可靠性和使用寿命及车辆行驶安全性。根据主簧各片和各级副簧的结构参数，弹性模量，最大许用应力，骑马螺栓夹紧距，主簧夹紧刚度及主簧与各级副簧的复合夹紧刚度，初始切线弧高，在最大许用载荷、接触载荷、渐变夹紧刚度仿真计算的基础上，对非等偏频型三级渐变刚度板簧的最大限位挠度进行仿真验算。
为解决上述技术问题，本发明所提供的非等偏频型三级渐变刚度板簧最大限位挠度的仿真验算法，其特征在于采用以下仿真验算步骤：
(1)主簧及其与各级副簧的根部重叠部分等效厚度和主簧最大厚度板簧的厚度的确定：
A步骤：主簧及其与各级副簧的根部重叠部分等效厚度的确定
根据主簧片数n，主簧各片的厚度hi，i＝1,2,...,n；第一级副簧片数n1，第一级副簧各片的厚度hA1j，j＝1,2,...,n1；第二级副簧片数n2,第二级副簧各片的厚度hA2k，k＝1,2,...,n2；第三级副簧片数n3,第三级副簧各片的厚度hA3l，l＝1,2,...,n3；对主簧的根部重叠部分的等效厚度 hMe、主簧与各级副簧的根部重叠部分的等效厚度hMA1e、hMA2e、hMA3e分别进行计算，即：
B步骤：主簧的最大厚度板簧的厚度hmax的确定
根据主簧片数n，主簧各片的厚度hi，i＝1,2,...,n，确定主簧的最大厚度板簧的厚度hmax，即
hmax＝max(hi)；
(2)非等偏频型三级渐变刚度板簧的主簧和各级副簧的曲率半径的计算：
I步骤：第一级主簧末片下表面初始曲率半径RM0b计算
根据主簧片数n，主簧各片的厚度hi，i＝1,2,…,n；主簧首片的一半夹紧长度L1，主簧的初始切线弧高HgM0，对主簧末片下表面初始曲率半径RM0b进行计算，即
II步骤：第一级副簧首片上表面初始曲率半径RA10a计算
根据第一级副簧首片的一半夹紧长度LA11，第一级副簧的初始切线弧高HgA10，对第一级副簧末片上表面初始曲率半径RA10a进行计算，即
III步骤：第一级副簧末片下表面初始曲率半径RA10b计算
根据第一级副簧片数n1，第一级副簧各片的厚度hA1j，j＝1,2,…,n1；II步骤中计算得到的 RA10a，对第一级副簧末片下表面初始曲率半径RA10b进行计算，即
IV步骤：第二级副簧首片上表面初始曲率半径RA20a的计算
根据第二级副簧首片的一半夹紧长度LA21，第二级副簧的初始切线弧高设计值HgA20，对第二级副簧首片上表面初始曲率半径RA20a进行计算，即
V步骤：第二级副簧首片下表面初始曲率半径RA20b的计算
很据第二级副簧片数n2，第二级副簧各片的厚度hA2k，k＝1,2,…,n2，及IV步骤所确定的 RA20a，对第二级副簧首片下表面初始曲率半径RA20b进行计算，即
VI步骤：第三级副簧首片上表面初始曲率半径RA30a的计算
根据第三级副簧首片的一半夹紧长度LA31，第三级副簧的初始切线弧高HgA30，对第三级副簧首片上表面初始曲率半径RA30a进行计算，即
(3)非等偏频型三级渐变刚度板簧的各次接触载荷的仿真计算：
a步骤：第1次开始接触载荷Pk1的仿真计算
根据非等偏频型三级渐变刚度板簧的宽度b，弹性模量E；主簧首片的一半夹紧长度L1，步骤(1)中确定的hMe，步骤(2)中计算得到的RM0b和RA10a，对第1次开始接触载荷Pk1进行仿真计算，即
b步骤：第2次开始接触载荷Pk2的仿真计算
根据非等偏频型三级渐变刚度板簧的宽度b，弹性模量E；主簧首片的一半夹紧长度L1，步骤(1)中确定的hMA1e，步骤(2)中计算得到的RA10b和RA20a，及a步骤中仿真计算得到的Pk1，对第2次开始接触载荷Pk2进行仿真计算，即
c步骤：第3次开始接触载荷Pk3的仿真计算
根据非等偏频型三级渐变刚度板簧的宽度b，弹性模量E；主簧首片的一半夹紧长度L1，步骤(1)中确定的hMA2e，步骤(2)中计算得到的RA20b和RA30a，及b步骤中仿真计算得到的Pk2，对第3次开始接触载荷Pk3进行仿真计算，即
d步骤：第3次完全接触载荷Pw3的仿真计算
根据b步骤中仿真计算得到的Pk2，c步骤中仿真计算得到的Pk3，对非等偏频型三级渐变刚度板簧的第3次完全接触载荷Pw3进行仿真计算，即
(4)非等偏频型三级渐变刚度板簧的最大许用载荷Pmax的仿真计算：
根据非等偏频型三级渐变刚度板簧的宽度b，最大许用应力[σ]；主簧首片的一半夹紧长度 L1，步骤(1)中所确定的hMe、hMA1e、hMA2e、hMA3e和hmax，步骤(3)中仿真计算得到的Pk1、Pk2和Pk3，对非等偏频型三级渐变刚度板簧的最大许用载荷Pmax进行仿真计算，即
(5)非等偏频型三级渐变刚度钢板弹簧的渐变夹紧刚度的仿真计算：
i步骤：第一级渐变夹紧刚度Kkwp1的仿真计算：
根据主簧夹紧刚度KM，主簧与第一级副簧的复合夹紧刚度KMA1；步骤(3)中仿真计算得到的 Pk1和Pk2，对载荷P在[Pk1,Pk2]范围内的第一级渐变夹紧刚度KkwP1进行仿真计算，即
ii步骤：第二级渐变夹紧刚度KkwP2的仿真计算
根据主簧与第一级副簧的复合夹紧刚度KMA1，主簧与第一级副簧和第二级副簧的复合夹紧刚度KMA2，步骤(3)中仿真计算得到的Pk2和Pk3，对载荷P在[Pk2,Pk3]范围内的第二级渐变夹紧刚度KkwP2进行仿真计算，即
iii步骤：第三级渐变夹紧刚度KkwP3的仿真计算
根据主簧与第一级副簧和第二级副簧的复合夹紧刚度KMA2，主副簧的总复合夹紧刚度 KMA3；步骤(3)中仿真计算得到的Pk3和Pw3，对载荷P在[Pk3,Pw3]范围内的第三级渐变夹紧刚度KkwP3进行仿真计算，即
(6)非等偏频型三级渐变刚度板簧的最大限位挠度fMmax的仿真验算：
根据主簧夹紧刚度KM，主副簧的总复合夹紧刚度KMA3；步骤(3)中仿真计算得到的Pk1， Pk2，Pk3和Pw3，步骤(4)中仿真计算得到的Pmax，步骤(5)中仿真计算得到的KkwP1、KkwP2和KkwP3，对非等偏频型三级渐变刚度板簧的最大限位挠度fMmax进行仿真验算，即
本发明比现有技术具有的优点
由于非等偏频型三级渐变刚度板簧的挠度计算非常复杂，并且受最大许用载荷和接触载荷仿真计算问题的制约，先前国内外一直未给出非等偏频型三级渐变刚度板簧最大限位挠度的仿真验算法，不能满足非等偏频型三级渐变刚度板簧设计及CAD软件开发要求。本发明可根据主簧各片和各级副簧的结构参数，弹性模量，最大许用应力，骑马螺栓夹紧距，主簧及其与各级副簧的夹紧刚度，初始切线弧高，在最大许用载荷和接触载荷及三级渐变夹紧刚度仿真计算的基础上，对非等偏频型三级渐变刚度板簧的最大限位挠度进行仿真验算。通过样机加载挠度试验可知，本发明所提供的非等偏频型三级渐变刚度板簧最大限位挠度的仿真验算法是正确的，为最大限位挠度仿真验算提供了可靠的技术方法。利用该方法可得到准确可靠的最大限位挠度仿真验算值，确保限位装置在冲击载荷下对板簧起板簧作用，可提高产品的设计水平、质量、可靠性和使用寿命及车辆行驶安全性；同时，降低设计及试验费用，加快产品开发速度。
附图说明
为了更好地理解本发明，下面结合附图做进一步的说明。
图1是非等偏频型三级渐变刚度板簧最大限位挠度的仿真验算流程图；
图2是非等偏频型三级渐变刚度板簧的一半对称结构示意图；
图3是实施例的三级渐变刚度钢板弹簧的夹紧刚度KP随载荷P的变化曲线。
具体实施方式
下面通过实施例对本发明作进一步详细说明。
实施例：某非等偏频型三级渐变刚度板簧的宽度b＝63mm，骑马螺栓夹紧距的一半L0＝50mm，弹性模量E＝200GPa，最大许用应力[σ]＝800MPa。主副簧的总片数N＝5，其中，主簧片数n＝2,主簧各片的厚度h1＝h2＝8mm；首片主簧的一半作用长度为L1T＝525mm，一半夹紧长度为L1＝L1T-L0/2＝500mm。第一级副簧的片数n1＝1，厚度hA11＝8mm。第二级副簧的片数n2＝1,厚度hA21＝13mm。第三级副簧的片数n3＝1，厚度hA31＝13mm。主簧夹紧刚度 KM＝51.43N/mm，主簧与第一级副簧的复合夹紧刚度KMA1＝75.4N/mm；主簧与第一级和第二级副簧的复合夹紧刚度KMA2＝144.5N/mm，主副簧的总复合夹紧刚度KMA3＝172.9N/mm。主簧的初始切线弧高HgM0＝102.3mm，第一级副簧的初始切线弧高HgA10＝18.8mm，第二级副簧的初始切线弧高HgA20＝6mm，第三级副簧的初始切线弧高HgA30＝1.6mm。根据主簧各片和各级副簧的结构参数，弹性模量，最大许用应力，骑马螺栓夹紧距，主簧夹紧刚度及其与各级副簧的复合夹紧刚度，各级板簧的初始切线弧高，对该非等偏频型三级渐变刚度板簧的最大限位挠度进行仿真验算。
本发明实例所提供的非等偏频型三级渐变刚度板簧最大限位挠度的仿真验算法，其仿真验算流程如图1所示，具体仿真验算步骤如下：
(1)主簧及其与各级副簧的根部重叠部分等效厚度和主簧最大厚度板簧的厚度的确定：
A步骤：主簧及其与各级副簧的根部重叠部分等效厚度的确定
根据主簧片数n＝2,主簧各片的厚度h1＝h2＝8mm；第一级副簧片数n1＝1，厚度hA11＝8mm；第二级副簧片数n2＝1,厚度hA21＝13mm；第三级副簧片数n3＝1,厚度hA31＝13mm；对主簧的根部重叠部分的等效厚度hMe、主簧与各级副簧的根部重叠部分的等效厚度hMA1e、hMA2e、hMA3e进行确定，即：
B步骤：主簧的最大厚度板簧的厚度的确定
根据主簧片数n＝2，主簧各片的厚度h1＝h2＝8mm，确定主簧的最大厚度板簧的厚度hmax，即
hmax＝max(h1,h2)＝8mm；
(2)非等偏频型三级渐变刚度板簧的主簧和各级副簧的曲率半径的计算：
I步骤：第一级主簧末片下表面初始曲率半径RM0b计算
根据主簧片数n＝2，主簧各片的厚度hi＝8mm，i＝1,2,...,n；主簧首片的一半夹紧长度 L1＝500mm，主簧的初始切线弧高HgM0＝102.3mm，对主簧末片下表面初始曲率半径RM0b进行计算，即
II步骤：第一级副簧首片上表面初始曲率半径RA10a计算
根据第一级副簧首片的一半夹紧长度LA11＝325mm，第一级副簧的初始切线弧高HgA10＝18.8mm，对第一级副簧末片上表面初始曲率半径RA10a进行计算，即
III步骤：第一级副簧末片下表面初始曲率半径RA10b计算
根据第一级副簧片数n1＝1，厚度hA11＝8mm；II步骤中计算得到的RA10a＝2818.6mm，对第一级副簧末片下表面初始曲率半径RA10b进行计算，即
RA10b＝RA10a+hA11＝2826.6mm；
IV步骤：第二级副簧首片上表面初始曲率半径RA20a的计算
根据第二级副簧首片的一半夹紧长度LA21＝225mm，第二级副簧的初始切线弧高HgA20＝6mm，对第二级副簧首片上表面初始曲率半径RA20a进行计算，即
V步骤：第二级副簧首片下表面初始曲率半径RA20b的计算
很据第二级副簧片数n2＝1，厚度hA21＝13mm，及IV步骤所确定的RA20a＝4221.8mm，对第二级副簧首片下表面初始曲率半径RA20b进行计算，即
RA20b＝RA20a+hA21＝4234.8mm；
VI步骤：第三级副簧首片上表面初始曲率半径RA30a的计算
根据第三级副簧首片的一半夹紧长度LA31＝125mm，第三级副簧的初始切线弧高HgA30＝1.6mm，对第三级副簧首片上表面初始曲率半径RA30a进行计算，即
(3)非等偏频型三级渐变刚度板簧的各次接触载荷的仿真计算：
a步骤：第1次开始接触载荷Pk1的仿真计算
根据非等偏频型三级渐变刚度板簧的宽度b＝63mm，弹性模量E＝200Gpa；主簧首片的一半夹紧长度L1＝500mm，步骤(1)中确定的hMe＝10.1mm，步骤(2)中计算得到的RM0b＝1289mm 和RA10a＝2818.6mm，对第1次开始接触载荷Pk1进行仿真计算，即
b步骤：第2次开始接触载荷Pk2的仿真计算
根据非等偏频型三级渐变刚度板簧的宽度b＝63mm，弹性模量E＝200GPa；首片主簧的一半夹紧长度L1＝500mm，步骤(1)中确定的hMA1e＝11.5mm，步骤(2)中计算得到的 RA10b＝2826.6mm和RA20a＝4221.8mm，及a步骤中仿真计算得到的Pk1＝1810N，对第2次开始接触载荷Pk2进行仿真计算，即
c步骤：第3次开始接触载荷Pk3的仿真计算
根据非等偏频型三级渐变刚度板簧的宽度b＝63mm，弹性模量E＝200GPa；主簧首片的一半夹紧长度L1＝500mm，步骤(1)中确定的hMA2e＝15.5mm，步骤(2)中所确定的RA20b＝4234.8mm 和RA30a＝4883.6mm，及b步骤中验算得到的Pk2＝2565N，对第3次开始接触载荷Pk3进行仿真计算，即
d步骤：第3次完全接触载荷Pw3的仿真计算
根据b步骤中仿真计算得到的Pk2＝2565N，c步骤中仿真计算得到的Pk3＝3057N，对该非等偏频型三级渐变刚度板簧的第3次完全接触载荷Pw3进行仿真计算，即
(4)非等偏频型三级渐变刚度板簧的最大许用载荷Pmax的仿真计算：
根据非等偏频型三级渐变刚度板簧的宽度b＝63mm，最大许用应力[σ]＝800MPa；主簧首片的一半夹紧长度L1＝500mm，步骤(1)中所确定的hMe＝10.1mm、hMA1e＝11.5mm、 hMA2e＝15.5mm、hMA3e＝18.1mm和hmax＝8mm，步骤(3)中仿真计算得到的Pk1＝1810N， Pk2＝2565N和Pk3＝3057N，对该非等偏频型三级渐变刚度板簧的最大许用载荷Pmax进行计算，即
(5)非等偏频型三级渐变刚度钢板弹簧的渐变夹紧刚度的仿真计算：
i步骤：第一级渐变夹紧刚度Kkwp1的仿真计算
根据主簧夹紧刚度KM＝51.43N/mm，主簧与第一级副簧的复合夹紧刚度KMA1＝75.4N/mm；步骤(3)中仿真计算得到的Pk1＝1810N和Pk2＝2565N，对载荷P在[Pk1,Pk2]范围内的第一级渐变夹紧刚度KkwP1进行仿真计算，即
ii步骤：第二级渐变夹紧刚度KkwP2的仿真计算
根据主簧与第一级副簧的复合夹紧刚度KMA1＝75.4N/mm，主簧与第一级和第二级副簧的复合夹紧刚度KMA2＝144.5N/mm，步骤(3)中仿真计算得到的Pk2＝2565N和Pk3＝3057N，对载荷 P在[Pk2,Pk3]范围内的第二级渐变夹紧刚度KkwP2进行仿真计算，即
iii步骤：第三级渐变夹紧刚度KkwP3的仿真计算
根据主簧与第一级和第二级副簧的复合夹紧刚度KMA2＝144.5N/mm，主副簧的总复合夹紧刚度KMA3＝172.9N/mm；步骤(3)中仿真计算得到的Pk3＝3057N和Pw3＝3643N，对载荷P在[Pk3, Pw3]范围内的第三级渐变夹紧刚度KkwP3进行仿真计算，即
利用Matlab计算程序，仿真计算得到的该三级渐变刚度钢板弹簧的夹紧刚度KP随载荷 P的变化曲线，如图3所示，其中，当载荷P<Pk1＝1810N时，渐变夹紧刚度 KP＝KM1＝51.4N/mm；当载荷P＝Pk2＝2565N时，KP＝KM2＝75.4N/mm,当载荷P＝Pk3＝3057N 时，KP＝KMA1＝144.5N/mm，当载荷P>PW3＝3643N时，KP＝KMA＝172.9N/mm。
(6)非等偏频型三级渐变刚度板簧的最大限位挠度fMmax的仿真验算:
根据主簧夹紧刚度KM＝51.4N/mm，主副簧的总复合夹紧刚度KMA3＝172.9N/mm；步骤(3)中仿真计算得到的Pk1＝1810N，Pk2＝2565N，Pk3＝3057N和Pw3＝3643N，步骤(4)中仿真计算得到的Pmax＝13339N，步骤(5)中仿真计算得到的KkwP1、KkwP2和KkwP3，对该非等偏频型三级渐变刚度板簧的最大限位挠度fMmax进行仿真验算，即
通过样机加载挠度试验可知，本发明所提供的非等偏频型三级渐变刚度板簧最大限位挠度的仿真验算法是正确的，为最大限位挠度仿真验算提供了可靠的技术方法。利用该方法可得到准确可靠的最大限位挠度仿真验算值，确保限位装置在冲击载荷下对板簧起保护作用，可提高产品的设计水平、质量、可靠性和使用寿命及车辆行驶安全性；同时，降低设计及试验费用，加快产品开发速度。

Claims (1)

1.非等偏频型三级渐变刚度板簧最大限位挠度的仿真验算法，其中，各片板簧为以中心穿装孔对称的结构，安装夹紧距的一半为骑马螺栓夹紧距的一半；板簧是由主簧和三级副簧构成，通过主簧和各级副簧的初始切线弧高及三级渐变间隙，确保满足板簧接触载荷、渐变刚度和主簧应力强度的设计要求，即非等偏频型三级渐变刚度板簧；依据最大限位挠度设计值设置一限位装置，在冲击载荷下对板簧起保护作用，提高板簧可靠性和使用寿命及车辆行驶安全性；根据各片板簧的结构参数，弹性模量，最大许用应力，骑马螺栓夹紧距，主簧夹紧刚度及主簧与各级副簧的复合夹紧刚度，初始切线弧高，在最大许用载荷、接触载荷、渐变夹紧刚度仿真计算的基础上，对非等偏频型三级渐变刚度板簧的最大限位挠度进行仿真验算，具体仿真验算步骤如下：
(1)主簧及其与各级副簧的根部重叠部分等效厚度和主簧最大厚度板簧的厚度的确定：
A步骤：主簧及其与各级副簧的根部重叠部分等效厚度的确定
根据主簧片数n，主簧各片的厚度hi，i＝1,2,...,n；第一级副簧片数n1，第一级副簧各片的厚度hA1j，j＝1,2,...,n1；第二级副簧片数n2,第二级副簧各片的厚度hA2k，k＝1,2,...,n2；第三级副簧片数n3,第三级副簧各片的厚度hA3l，l＝1,2,...,n3；对主簧的根部重叠部分的等效厚度hMe、主簧与各级副簧的根部重叠部分的等效厚度hMA1e、hMA2e、hMA3e分别进行计算，即：
B步骤：主簧的最大厚度板簧的厚度hmax的确定
根据主簧片数n，主簧各片的厚度hi，i＝1,2,...,n，确定主簧的最大厚度板簧的厚度hmax，即
hmax＝max(hi)；
(2)非等偏频型三级渐变刚度板簧的主簧和各级副簧的曲率半径的计算：
I步骤：第一级主簧末片下表面初始曲率半径RM0b计算
根据主簧片数n，主簧各片的厚度hi，i＝1,2,…,n；主簧首片的一半夹紧长度L1，主簧的初始切线弧高HgM0，对主簧末片下表面初始曲率半径RM0b进行计算，即
II步骤：第一级副簧首片上表面初始曲率半径RA10a计算
根据第一级副簧首片的一半夹紧长度LA11，第一级副簧的初始切线弧高HgA10，对第一级副簧首 片上表面初始曲率半径RA10a进行计算，即
III步骤：第一级副簧末片下表面初始曲率半径RA10b计算
根据第一级副簧片数n1，第一级副簧各片的厚度hA1j，j＝1,2,…,n1；II步骤中计算得到的RA10a，对第一级副簧末片下表面初始曲率半径RA10b进行计算，即
IV步骤：第二级副簧首片上表面初始曲率半径RA20a的计算
根据第二级副簧首片的一半夹紧长度LA21，第二级副簧的初始切线弧高设计值HgA20，对第二级副簧首片上表面初始曲率半径RA20a进行计算，即
V步骤：第二级副簧首片下表面初始曲率半径RA20b的计算
很据第二级副簧片数n2，第二级副簧各片的厚度hA2k，k＝1,2,…,n2，及IV步骤所确定的RA20a，对第二级副簧首片下表面初始曲率半径RA20b进行计算，即
VI步骤：第三级副簧首片上表面初始曲率半径RA30a的计算
根据第三级副簧首片的一半夹紧长度LA31，第三级副簧的初始切线弧高HgA30，对第三级副簧首片上表面初始曲率半径RA30a进行计算，即
(3)非等偏频型三级渐变刚度板簧的各次接触载荷的仿真计算：
a步骤：第1次开始接触载荷Pk1的仿真计算
根据非等偏频型三级渐变刚度板簧的宽度b，弹性模量E；主簧首片的一半夹紧长度L1，步骤(1)中确定的hMe，步骤(2)中计算得到的RM0b和RA10a，对第1次开始接触载荷Pk1进行仿真计算，即
b步骤：第2次开始接触载荷Pk2的仿真计算
根据非等偏频型三级渐变刚度板簧的宽度b，弹性模量E；主簧首片的一半夹紧长度L1，步骤(1)中确定的hMA1e，步骤(2)中计算得到的RA10b和RA20a，及a步骤中仿真计算得到的Pk1，对第2次开始接触载荷Pk2进行仿真计算，即
c步骤：第3次开始接触载荷Pk3的仿真计算
根据非等偏频型三级渐变刚度板簧的宽度b，弹性模量E；主簧首片的一半夹紧长度L1，步骤(1)中确定的hMA2e，步骤(2)中计算得到的RA20b和RA30a，及b步骤中仿真计算得到的Pk2，对第3次开始接触载荷Pk3进行仿真计算，即
d步骤：第3次完全接触载荷Pw3的仿真计算
根据b步骤中仿真计算得到的Pk2，c步骤中仿真计算得到的Pk3，对非等偏频型三级渐变刚度板簧的第3次完全接触载荷Pw3进行仿真计算，即
(4)非等偏频型三级渐变刚度板簧的最大许用载荷Pmax的仿真计算：
根据非等偏频型三级渐变刚度板簧的宽度b，最大许用应力[σ]；主簧首片的一半夹紧长度L1，步骤(1)中所确定的hMe、hMA1e、hMA2e、hMA3e和hmax，步骤(3)中仿真计算得到的Pk1、Pk2和Pk3，对非等偏频型三级渐变刚度板簧的最大许用载荷Pmax进行仿真计算，即
(5)非等偏频型三级渐变刚度钢板弹簧的渐变夹紧刚度的仿真计算：
i步骤：第一级渐变夹紧刚度Kkwp1的仿真计算：
根据主簧夹紧刚度KM，主簧与第一级副簧的复合夹紧刚度KMA1；步骤(3)中仿真计算得到的Pk1和Pk2，对载荷P在[Pk1,Pk2]范围内的第一级渐变夹紧刚度KkwP1进行仿真计算，即
ii步骤：第二级渐变夹紧刚度KkwP2的仿真计算
根据主簧与第一级副簧的复合夹紧刚度KMA1，主簧与第一级副簧和第二级副簧的复合夹紧刚度KMA2，步骤(3)中仿真计算得到的Pk2和Pk3，对载荷P在[Pk2,Pk3]范围内的第二级渐变夹紧刚度KkwP2进行仿真计算，即
iii步骤：第三级渐变夹紧刚度KkwP3的仿真计算
根据主簧与第一级副簧和第二级副簧的复合夹紧刚度KMA2，主副簧的总复合夹紧刚度KMA3；步骤(3)中仿真计算得到的Pk3和Pw3，对载荷P在[Pk3,Pw3]范围内的第三级渐变夹紧刚度KkwP3进行仿真计算，即
(6)非等偏频型三级渐变刚度板簧的最大限位挠度fMmax的仿真验算：
根据主簧夹紧刚度KM，主副簧的总复合夹紧刚度KMA3；步骤(3)中仿真计算得到的Pk1，Pk2，Pk3和Pw3，步骤(4)中仿真计算得到的Pmax，步骤(5)中仿真计算得到的KkwP1、KkwP2和KkwP3，对非等偏频型三级渐变刚度板簧的最大限位挠度fMmax进行仿真验算，即
CN201710023283.5A 2017-01-12 2017-01-12 非等偏频型三级渐变刚度板簧最大限位挠度的仿真验算法 Expired - Fee Related CN106812851B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710023283.5A CN106812851B (zh) 2017-01-12 2017-01-12 非等偏频型三级渐变刚度板簧最大限位挠度的仿真验算法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710023283.5A CN106812851B (zh) 2017-01-12 2017-01-12 非等偏频型三级渐变刚度板簧最大限位挠度的仿真验算法

Publications (2)

Publication Number Publication Date
CN106812851A CN106812851A (zh) 2017-06-09
CN106812851B true CN106812851B (zh) 2019-02-26

Family

ID=59109779

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710023283.5A Expired - Fee Related CN106812851B (zh) 2017-01-12 2017-01-12 非等偏频型三级渐变刚度板簧最大限位挠度的仿真验算法

Country Status (1)

Country Link
CN (1) CN106812851B (zh)

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102094921A (zh) * 2009-12-09 2011-06-15 湖南易通汽车配件科技发展有限公司 渐变钢度钢板弹簧
CN202022134U (zh) * 2011-01-12 2011-11-02 三一重工股份有限公司 可变刚度少片弹簧总成及混凝土搅拌运输车
CN102734364A (zh) * 2012-07-17 2012-10-17 山东理工大学 汽车钢板弹簧弧高与曲面形状的解析设计方法
CN105673753B (zh) * 2016-03-25 2018-06-08 山东恒日悬架弹簧有限公司 端部和根部加强型少片变截面钢板弹簧的设计方法

Also Published As

Publication number Publication date
CN106812851A (zh) 2017-06-09

Similar Documents

Publication Publication Date Title
CN106812851B (zh) 非等偏频型三级渐变刚度板簧最大限位挠度的仿真验算法
CN106812849B (zh) 非等偏频型三级渐变刚度板簧的接触载荷的验算方法
CN106777804B (zh) 基于偏频仿真的三级渐变刚度板簧接触载荷的调整设计法
CN106802996A (zh) 两级副簧式非等偏频型渐变刚度板簧接触载荷的验算方法
CN106704429B (zh) 非等偏频型三级渐变刚度板簧的最大限位挠度的设计方法
CN106599525B (zh) 非等偏频型三级渐变刚度板簧悬架偏频特性的仿真计算法
CN106682357B (zh) 高强度三级渐变刚度板簧悬架系统偏频特性的仿真计算方法
CN106802998B (zh) 非等偏频型三级渐变刚度板簧夹紧刚度特性的仿真计算法
CN106763390B (zh) 非等偏频型三级渐变刚度板簧的挠度特性的仿真计算法
CN106812846B (zh) 基于偏频仿真的两级主簧式渐变刚度板簧接触载荷调整设计法
CN106548003B (zh) 非等偏频型三级渐变刚度板簧根部最大应力的仿真计算法
CN106763384B (zh) 两级副簧式非等偏频型渐变刚度板簧切线弧高的设计方法
CN106777793B (zh) 两级副簧式非等偏频型渐变刚度板簧刚度特性的计算方法
CN106599524B (zh) 非等偏频型三级渐变刚度板簧的初始切线弧高的设计方法
CN106777806B (zh) 高强度等偏频三级渐变刚度板簧接触载荷的验算方法
CN106763385B (zh) 两级主簧式非等偏频型渐变刚度板簧限位挠度的设计方法
CN106678224B (zh) 等渐变偏频两级渐变刚度板簧最大限位挠度的仿真验算法
CN106812845B (zh) 两级主簧式非等偏频型渐变刚度板簧限位挠度的仿真验算法
CN106548002B (zh) 非等偏频型三级渐变刚度板簧的渐变间隙的设计方法
CN106812847B (zh) 非等偏频型三级渐变刚度板簧的强度校核法
CN106763381B (zh) 高强度一级渐变刚度板簧最大限位挠度的设计方法
CN106763387B (zh) 高强度三级渐变刚度板簧最大限位挠度的仿真验算法
CN106682359B (zh) 两级副簧式非等偏频型渐变刚度板簧主簧挠度的计算方法
CN106682360B (zh) 高强度两级渐变刚度主副簧的最大应力特性的仿真计算法
CN106650174B (zh) 高强度三级渐变刚度板簧的各次接触载荷的仿真计算方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190226

Termination date: 20210112