CN106777790A - 两级主簧式非偏频式渐变刚度板簧刚度特性的仿真计算法 - Google Patents

Abstract

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

Description

两级主簧式非偏频式渐变刚度板簧刚度特性的仿真计算法

技术领域

本发明涉及车辆悬架钢板弹簧，特别是两级主簧式非偏频式渐变刚度板簧刚度特性的仿真计算法。

背景技术

为了进一步提高车辆在半载情况下的行驶平顺性，可采用两级主簧式渐变刚度板簧，即将原一级渐变刚度板簧的主簧拆分为两级主簧；同时，为了确保主簧的应力强度，通常通过第一级主簧、第二级主簧和副簧的初始切线弧高及两级渐变间隙，使第二级主簧和副簧适当提前承担载荷，即给次接触载荷适当提前，从而降低第一级主簧的应力，即采用两级主簧式非等偏频型渐变刚度板簧悬架，其中，板簧的夹紧刚度特性不仅与各片板簧的结构有关，而且与各次接触载荷有关，并且影响悬架偏频特性及车辆行驶平顺性和安全性。然而，由于受两级主簧式非等偏频渐变刚度板簧的各级板簧夹紧刚度和渐变刚度计算及接触载荷的仿真计算的制约，先前一直未能给出两级主簧式非偏频式渐变刚度板簧刚度特性的仿真计算法，因此，不能满足车辆行业快速发展及悬架弹簧现代化CAD设计要求。随着车辆行驶速度及其对平顺性要求的不断提高，对渐变刚度板簧悬架提出了更高要求，因此，必须建立一种精确、可靠的两级主簧式非偏频式渐变刚度板簧刚度特性的仿真计算法，为两级主簧式非等偏频型渐变刚度板簧的特性仿真及CAD软件开发奠定可靠的技术基础，满足车辆行业快速发展、车辆行驶平顺性及对渐变刚度板簧的设计要求，提高两级主簧式非等偏频型渐变刚度板簧的设计水平、产品质量和性能、悬架可靠性及车辆行驶平顺性；同时，降低设计及试验费用，加快产品开发速度。

发明内容

针对上述现有技术中存在的缺陷，本发明所要解决的技术问题是提供一种简便、可靠的两级主簧式非偏频式渐变刚度板簧刚度特性的仿真计算法，仿真验算流程如图1所示。两级主簧式非等偏频型渐变刚度板簧的一半对称结构如图2所示，是由第一级主簧1、第二级主簧2和副簧3组成。采用两级主簧，并通过第一级主簧1、第二级主簧2和副簧的初始切线弧高H_gM10、H_gM20和H_gA0，在第一级主簧1与第二级主簧2和第二级主簧2与副簧3之间设有两级渐变间隙δ_M12和δ_MA，以提高半载情况下的车辆行驶平顺性。为了确保满足第一级主簧1应力强度设计要求，第二级主簧2和副簧3适当提前承担载荷，悬架渐变载荷偏频不相等，即将板簧设计为非等偏频型渐变刚度板簧。板簧的一半跨度等于第一级主簧首片的一半作用长度L_11T，骑马螺栓夹紧距的一半为L₀，宽度为b，弹性模量为E。第一级主簧1的片数为n₁，各片一级主簧的厚度为h_1i，一半作用长度为L_1iT，一半夹紧长度L_1i＝L_i＝L_iT-L₀/2，i＝1,2,…,n₁。第二级主簧2的片数为n₂，各片二级主簧的厚度为h_2j，一半作用长度为L_2jT，一半夹紧长度L_2j＝L_n1+j＝L_2jT-L₀/2，j＝1,2,…,n₂。第一级和第二级主簧的片数之和n＝n₁+n₂。副簧3的片数为m，各片副簧的厚度为h_Ak，一半作用长度为L_AkT，一半夹紧长度L_Ak＝L_n+k＝L_AkT-L₀/2，k＝1,2,…,m。主副簧的总片数N＝n₁+n₂+m。板簧的夹紧刚度特性不仅与各片板簧的结构有关，而且与各次接触载荷有关，并且影响悬架偏频特性及车辆行驶平顺性和安全性。根据各片第一级和第二级主簧及副簧的结构参数，弹性模量，骑马螺栓夹紧距，初始切线弧高，在各次接触载荷、夹紧刚度和渐变刚度的仿真计算的基础上，对两级主簧式非等偏频型渐变刚度板簧在不同载荷下的夹紧刚度特性进行仿真计算。

为解决上述技术问题，本发明所提供的两级主簧式非偏频式渐变刚度板簧刚度特性的仿真计算法，其特征在于采用以下仿真计算步骤：

(1)两级主簧式非等偏频渐变刚度板簧的不同片数重叠段的等效厚度进行计算：

根据第一级主簧片数n₁,第一级主簧各片的厚度h_1i，i＝1,2,…,n₁；第二级主簧片数n₂，第二级主簧各片的厚度h_2j，j＝1,2,…,n₂；副簧片数m，副簧各片的厚度h_Ak，k＝1,2,…,m；第一级与第二级主簧的片数之和n＝n₁+n₂，主副簧的总片数N＝n₁+n₂+m，对两级主簧式非等偏频渐变刚度板簧的不同片数l重叠段的等效厚度h_le进行计算，l＝1,2,…,N,即

其中，第一级主簧的根部重叠部分的等效厚度第一级主簧与第二级主簧的根部重叠部分的等效厚度h_M2e＝h_ne，主副簧的根部重叠部分的总等效厚度h_MAe＝h_Ne；

(2)两级主簧式非等偏频渐变刚度板簧的各级板簧初始曲率半径的计算：

I步骤：第一级主簧末片下表面初始曲率半径R_M10b计算

根据第一级主簧片数n₁，第一级主簧各片的厚度h_1i，i＝1,2,…,n₁；第一级主簧首片的一半夹紧长度L₁₁，第一级主簧的初始切线弧高H_gM10，对第一级主簧末片下表面初始曲率半径R_M10b进行计算，即

II步骤：第二级主簧首片上表面初始曲率半径R_M20a计算

根据第二级主簧首片的一半夹紧长度L₂₁，第二级主簧的初始切线弧高H_gM20，对第二级主簧末片上表面初始曲率半径R_M20a进行计算，即

III步骤：第二级主簧末片下表面初始曲率半径R_M20b计算

根据第二级主簧片数n₂，第二级主簧各片的厚度h_2j，j＝1,2,…,n₂；II步骤中计算得到的R_M20a，对第二级主簧末片下表面初始曲率半径R_M20b进行计算，即

IV步骤：副簧首片上表面初始曲率半径R_A0a计算

根据副簧首片的一半夹紧长度L_A1，副簧的初始切线弧高H_gA0，对副簧末片上表面初始曲率半径R_A0a进行计算，即

(3)两级主簧式非等偏频型渐变刚度板簧的各次接触载荷P_k1、P_k2、P_w2的仿真计算：

A步骤：第1次开始接触载荷P_k1的仿真计算

根据两级主簧式非等偏频型渐变刚度板簧的宽度b，弹性模量E；第一级主簧片数n₁，第一级主簧各片的厚度h_1i，i＝1,2,…,n₁，第一级主簧首片的一半夹紧跨长度L₁₁，步骤(1)中计算得到的h_M1e，步骤(2)中计算得到的R_M10b和R_M20a，对第1次开始接触载荷P_k1进行仿真计算，即

B步骤：第2次开始接触载荷P_k2的仿真计算

根据两级主簧式非等偏频型渐变刚度板簧的宽度b，弹性模量E；第一级主簧首片的一半夹紧跨长度L₁₁；第二级主簧片数n₂，第二级主簧各片的厚度h_2j，j＝1,2,…,n₂；步骤(1)中计算得到的h_M2e，步骤(2)中计算得到的R_M20b和R_A0a，及A步骤中仿真计算得到的P_k1，对第2次开始接触载荷P_k2进行仿真计算，即

C步骤：第2次完全接触载荷P_w2的仿真计算

根据A步骤中仿真计算得到的P_k1，B步骤中仿真计算得到的P_k2，对第2次完全接触载荷P_w2进行验算，即

(4)两级主簧的渐变刚度板簧的各级夹紧刚度的仿真计算：

a步骤：第一级主簧夹紧刚度K_M1的仿真计算：

根据两级主簧式非等偏频型渐变刚度板簧的宽度b，弹性模量E；第一级主簧片数n₁,第一级主簧各片主簧的一半夹紧长度L_1i＝L_i，i＝1,2,…,n₁，及步骤(1)中计算得到的h_le，l＝i＝1,2,…,n₁，对第一级主簧夹紧刚度K_M1进行仿真计算，即

b步骤：第一级主簧与第二级主簧的复合夹紧刚度K_M2的仿真计算：

根据两级主簧式非等偏频型渐变刚度板簧的宽度b，弹性模量E；第一级主簧片数n₁，第一级主簧各片的一半夹紧长度L_1i＝L_i，i＝1,2,…,n₁；第二级主簧片数n₂,第二级主簧各片的一半夹紧长度L_2j＝L_n1+j，j＝1,2,…,n₂；第一级主簧和第二主簧的总片数n＝n₁+n₂，及步骤(1)中计算得到的h_le，l＝1,2,…,n，对第一级主簧与第二级主簧的复合夹紧刚度K_M2进行仿真计算，即

c步骤：主副簧总复合夹紧刚度K_MA的仿真计算：

根据两级主簧式非等偏频型渐变刚度板簧的宽度b，弹性模量E；第一级主簧的片数n₁，第一级主簧各片的一半夹紧长度L_1i＝L_i，i＝1,2,…,n₁；第二级主簧的片数n₂，第二级主簧各片的一半夹紧长度L_2j＝L_n1+j，j＝1,2,…,n₂；副簧片数m；副簧各片的一半夹紧长度分别为L_A1＝L_n+k，k＝1,2,…,m，主副簧的总片数N＝n₁+n₂+m；及步骤(1)中计算得到的h_le，l＝1,2,…,N，对主副簧的总夹紧复合刚度K_MA进行仿真计算，即，即

(5)两级主簧式非等偏频型渐变刚度板簧的渐变夹紧刚度的仿真计算：

i步骤：第一级渐变夹紧刚度K_kwP1的仿真计算

根据步骤(3)中仿真计算得到的P_k1和P_k2，步骤(4)中仿真计算得到的K_M1和K_M2；对载荷P在[P_k1,P_k2]范围内的第一级渐变夹紧刚度K_kwP1进行仿真计算，即

ii步骤：第二级渐变夹紧刚度K_kwP2的仿真计算

根据步骤(3)中仿真计算得到的P_k2和P_w2，步骤(4)中仿真计算得到的K_M2和K_MA；对载荷P在[P_k2,P_w2]范围内的第二级渐变夹紧刚度K_kwP2进行仿真计算，即

(6)两级主簧式非等偏频型渐变刚度板簧的夹紧刚度特性的仿真计算：

根据步骤(3)中仿真计算得到的P_k1、P_k2和P_w2；步骤(4)中仿真计算得到的K_M1、K_M2和K_MA，及步骤(5)骤中仿真计算得到的K_kwP1和K_kwP2，对两级主簧式非等偏频型渐变刚度板簧在不同载荷下的夹紧刚度特性进行仿真计算，即

本发明比现有技术具有的优点

由于受两级主簧式非等偏频渐变刚度板簧的各级板簧夹紧刚度和渐变刚度计算及接触载荷的仿真计算的制约，先前一直未能给出两级主簧式非偏频式渐变刚度板簧刚度特性的仿真计算法，因此，不能满足车辆行业快速发展及悬架弹簧现代化CAD设计要求。本发明可根据各片第一级和第二级主簧及副簧的结构参数，弹性模量，骑马螺栓夹紧距，初始切线弧高，及额定载荷，在各次接触载荷、夹紧刚度和渐变刚度的仿真计算的基础上，对两级主簧式非等偏频型渐变刚度板簧在不同载荷下的夹紧刚度特性进行仿真计算。通过样机的ANSYS仿真和加载挠度及刚度试验测试可知，本发明所提供的两级主簧式非偏频式渐变刚度板簧刚度特性的仿真计算法是正确的，利用该方法可得到可靠的在不同载荷下的夹紧刚度仿真计算值，为两级主簧式非偏频式渐变刚度板簧的夹紧刚度特性仿真提供了可靠的技术方法。利用该方法可提高产品设计水平、可靠性及车辆行驶安全性；同时，降低产品设计及试验费用，加快产品开发速度。

附图说明

为了更好地理解本发明，下面结合附图做进一步的说明。

图1是两级主簧式非偏频式渐变刚度板簧刚度特性的仿真计算流程图；

图2是两级主簧式非等偏频渐变刚度板簧的一半对称结构示意图；

图3是实施例的两级主簧式非等偏频型渐变刚度板簧的夹紧刚度随载荷的变化特性曲线。

具体实施方案

下面通过实施例对本发明作进一步详细说明。

实施例：某两级主簧式非等偏频渐变刚度板簧的宽度b＝63mm，骑马螺栓夹紧距的一半L₀＝50mm，弹性模量E＝200GPa，最大许用应力[σ]＝800MPa。第一级主簧片数n₁＝2，第一级主簧各片的厚度h₁₁＝h₁₂＝8mm，第一级主簧各片的一半作用长度L_11T＝525mm，L_12T＝450mm，一半夹紧长度L₁＝500mm，L₂＝425mm。第二级主簧片数n₂＝1，厚度h₂₁＝8mm，第二级主簧首片的一半作用长度L_21T＝350mm，一半夹紧长度L₂₁＝L₃＝L_21T-L₀/2＝325mm。第一级主簧与第二级主簧的片数之和n＝n₁+n₂＝3。副簧片数m＝2，副簧各片的厚度h_A1＝h_A2＝13mm；副簧各片的一半作用长度L_A1T＝250mm，L_A2T＝150mm；一半夹紧长度分别为L_A1＝L₄＝L_A1t-L₀/2＝225mm，L_A2＝L₅＝L_A2t-L₀/2＝125mm。第一级主簧的初始切线弧高设计值H_gM10＝103.7mm，第二级主簧的初始切线弧高H_gM20＝18.8mm，副簧的初始切线弧高H_gA0＝6mm。根据各片第一级和第二级主簧及副簧的结构参数，弹性模量，骑马螺栓夹紧距，初始切线弧高，在各次接触载荷、夹紧刚度和渐变刚度的仿真计算的基础上，对该两级主簧式非等偏频型渐变刚度板簧在不同载荷下的夹紧刚度特性进行仿真计算。

本发明实例所提供的两级主簧式非偏频式渐变刚度板簧刚度特性的仿真计算法，其仿真验算流程如图1所示，具体仿真验算步骤如下：

(1)两级主簧式非等偏频渐变刚度板簧的不同片数重叠段的等效厚度进行计算：

根据第一级主簧片数n₁＝2,各片厚度h₁₁＝h₁₂＝8mm；第二级主簧片数n₂＝1，厚度h₂₁＝8mm；副簧片数m＝2,各片副簧厚度h_A1＝h_A2＝13mm；第一级和第二级主簧的片数之和n＝n₁+n₂＝3，板簧的总片数N＝n+m＝5，其中，对两级主簧式非等偏频渐变刚度板簧的不同片数l重叠段的等效厚度h_le进行计算，l＝1,2,…,N,即

h_1e＝h₁＝8.0mm；

其中，第一级主簧的根部重叠部分的等效厚度h_M1e＝h_2e＝10.1mm；第一级主簧与第二级主簧的根部重叠部分的等效厚度h_M2e＝h_3e＝11.5mm；主副簧的根部重叠部分的等效厚度h_MAe＝h_5e＝18.1mm。

(2)两级主簧式非等偏频渐变刚度板簧的各级板簧初始曲率半径的计算：

I步骤：第一级主簧末片下表面初始曲率半径R_M10b计算

根据第一级主簧片数n₁＝2，各片主簧的厚度h_1i＝8mm，i＝1,2,…,n₁；第一级主簧首片的一半夹紧长度L₁₁＝500mm，第一级主簧的初始切线弧高H_gM10＝103.7mm，对第一级主簧末片下表面的曲率半径R_M10b进行计算，即

II步骤：第二级主簧首片上表面初始曲率半径R_M20a计算

根据第二级主簧首片的一半夹紧长度L₂₁＝325mm，第二级主簧的初始切线弧高H_gM20＝18.8mm，对第二级主簧末片上表面初始曲率半径R_M20a进行计算，即

III步骤：第二级主簧末片下表面初始曲率半径R_M20b计算

根据第二级主簧片数n₂＝1，厚度h₂₁＝8mm；II步骤中计算得到的R_M20a＝2818.6mm，对第二级主簧末片下表面初始曲率半径R_M20b进行计算，即

R_M20b＝R_M20a+h₂₁＝2826.6mm；

IV步骤：副簧首片上表面初始曲率半径R_A0a计算

根据副簧首片的一半夹紧长度L_A1＝225mm，副簧的初始切线弧高H_gA0＝6mm，对副簧末片上表面初始曲率半径R_A0a进行计算，即

(3)两级主簧式非等偏频型渐变刚度板簧的各次接触载荷P_k1、P_k2、P_w2的仿真计算：

A步骤：第1次开始接触载荷P_k1的仿真计算

根据两级主簧式非等偏频型渐变刚度板簧的宽度b＝63mm，弹性模量E＝200GPa，第一级主簧首片的一半夹紧跨长度L₁₁＝500mm，步骤(1)中计算得到的h_M1e＝10.1mm，步骤(2)中计算得到的R_M10b＝1273.3mm和R_M20a＝2818.6mm，对第1次开始接触载荷P_k1进行仿真计算，即

B步骤：第2次开始接触载荷P_k2的仿真计算

根据两级主簧式非等偏频型渐变刚度板簧的宽度b＝63mm，弹性模量E＝200GPa，第一级主簧首片的一半夹紧跨长度L₁₁＝500mm；第二级主簧片数n₂＝1，厚度h₂₁＝8mm；步骤(1)中计算得到的h_M2e＝11.5mm；步骤(2)中计算得到的R_M20b＝2826.6mm和R_A0a＝4221.8mm，及A步骤中仿真计算得到的P_k1＝1851N，对第2次开始接触载荷P_k2进行仿真计算，即

C步骤：第2次完全接触载荷P_w2的仿真计算

根据A步骤中仿真计算得到的P_k1＝1851N，B步骤中仿真计算得到的P_k2＝2606N，对第2次完全接触载荷P_w2进行仿真计算，即

(4)两级主簧的渐变刚度板簧的各级夹紧刚度的仿真计算：

a步骤：第一级主簧夹紧刚度K_M1的仿真计算：

根据渐变刚度板簧的宽度b＝63mm，弹性模量E＝200GPa；第一级主簧片数n₁＝2,其中，各片主簧的一半夹紧长度L₁＝500mm，L₂＝425mm，及步骤(1)中计算得到的h_1e＝8.0mm，h_2e＝10.1mm，对第一级主簧夹紧刚度K_M1进行仿真计算，即

b步骤：第一级主簧与第二级主簧夹紧复合刚度K_M2的仿真计算：

根据渐变刚度板簧的宽度b＝63mm，弹性模量E＝200GPa；第一级主簧片数n₁＝2,第一级主簧各片的一半夹紧长度L₁＝500mm，L₂＝425mm；第二级主簧片数n₂＝1,一半夹紧长度L₂₁＝L₃＝325mm；第一级主簧与第二主簧的片数之和n＝3，及步骤(1)中计算得到的h_1e＝8.0mm，h_2e＝10.1mm，h_3e＝11.5mm，对第一级主簧与第二级主簧的夹紧复合刚度K_M2进行仿真计算，即

c步骤：主副簧总复合夹紧刚度K_MA的仿真计算：

根据两级主簧式非等偏频型渐变刚度板簧的宽度b＝63mm，弹性模量E＝200GPa，第一级主簧的片数n₁＝2，第一级主簧各片的一半夹紧长度L₁＝500mm，L₂＝425mm；第二级主簧的片数n₂＝1，一半夹紧长度L₂₁＝L₃＝325mm；副簧片数m＝2，副簧各片的一半夹紧长度分别为L_A1＝L₄＝225mm，L_A2＝L₅＝125mm；主副簧的总片数N＝5，其中，及步骤(1)中计算得到的h_1e＝8.0mm，h_2e＝10.1mm，h_3e＝11.5mm，h_4e＝15.5mm，h_5e＝18.1mm，对主副簧的总夹紧复合刚度K_MA进行仿真计算，即，即

(5)两级主簧式非等偏频型渐变刚度板簧的渐变夹紧刚度的仿真计算：

i步骤：第一级渐变夹紧刚度K_kwP1的仿真计算

根据步骤(3)中仿真计算得到的P_k1＝1851N和P_k2＝2602N，步骤(4)中仿真计算得到的K_M1＝51.43N/mm和K_M2＝75.4N/mm；对载荷P在[P_k1,P_k2]范围内的第一级渐变夹紧刚度K_kwP1进行仿真计算，即

ii步骤：二级渐变夹紧刚度K_kwP2的仿真计算

根据步骤(3)中仿真计算得到的P_k2＝2602N和P_w2＝3667N，步骤(4)中仿真计算得到的K_M2＝75.4N/mm和K_MA＝172.9N/mm；对载荷P在[P_k2,P_w2]范围内的第二级渐变夹紧刚度K_kwP2进行仿真计算，即

(6)两级主簧式非等偏频型渐变刚度板簧的夹紧刚度特性的仿真计算：

根据步骤(3)中仿真计算得到的P_k1＝1851N、P_k2＝2606N和P_w2＝3667N；步骤(4)中仿真计算得到的K_M1＝51.43N/mm、K_M2＝75.4N/mm和K_MA＝172.9N/mm，及步骤(5)骤中仿真计算得到的K_kwP1和K_kwP2，对该两级主簧式非等偏频型渐变刚度板簧在不同载荷下的夹紧刚度特性进行仿真计算，即

利用Matlab计算程序，仿真计算得到的该两级主簧式非等偏频型渐变刚度板簧的夹紧刚度随载荷的变化特性曲线，如图3所示，其中，在P_k1、P_k2、P_w2载荷下的夹紧刚度分别为K_k1＝51.43N/mm，K_k2＝75.4N/mm，K_w2＝172.9N/mm，在渐变接触过程中，夹紧刚度K随载荷P而变化。

通过样机的ANSYS仿真和加载挠度及刚度试验测试可知，本发明所提供的两级主簧式非偏频式渐变刚度板簧刚度特性的仿真计算法是正确的，利用该方法可得到可靠的在不同载荷下的夹紧刚度仿真计算值，为两级主簧式非偏频式渐变刚度板簧的夹紧刚度特性仿真提供了可靠的技术方法。利用该方法可提高产品设计水平、可靠性及车辆行驶安全性；同时，降低产品设计及试验费用，加快产品开发速度。

Claims (1)

1.两级主簧式非偏频式渐变刚度板簧刚度特性的仿真计算法，其中，各片板簧为以中心穿装孔对称的结构，安装夹紧距的一半为骑马螺栓夹紧距的一半；将原一级渐变刚度板簧的主簧拆分设计为两级主簧，通过两级主簧和副簧的初始切线弧高及两级渐变间隙，提高半载情况下的车辆行驶平顺性；同时，为了确保满足第一级主簧应力强度设计要求，第二级主簧和副簧适当提前承担载荷，悬架在渐变载荷下的偏频不相等，即两级主簧式非等偏频型渐变刚度板簧；根据各片板簧的结构参数，弹性模量，骑马螺栓夹紧距，初始切线弧高，在各次接触载荷、夹紧刚度和渐变刚度的仿真计算的基础上，对两级主簧式非等偏频型渐变刚度板簧在不同载荷下的夹紧刚度特性进行仿真计算，具体仿真计算步骤如下：

(1)两级主簧式非等偏频渐变刚度板簧的不同片数重叠段的等效厚度进行计算：

根据第一级主簧片数n₁,第一级主簧各片的厚度h_1i，i＝1,2,…,n₁；第二级主簧片数n₂，第二级主簧各片的厚度h_2j，j＝1,2,…,n₂；副簧片数m，副簧各片的厚度h_Ak，k＝1,2,…,m；第一级与第二级主簧的片数之和n＝n₁+n₂，主副簧的总片数N＝n₁+n₂+m，对两级主簧式非等偏频渐变刚度板簧的不同片数l重叠段的等效厚度h_le进行计算，l＝1,2,…,N,即

$h_{le}=\left\{\begin{array}{cc}\sqrt[3]{\underset{i=1}{\overset{l}{\Σ}}{h}_{1i}^3},& 1\≤l\≤n_1\\ \sqrt[3]{\underset{i=1}{\overset{n_1}{\Σ}}{h}_{1i}^3+\underset{j=1}{\overset{l-n_1}{\Σ}}{h}_{2j}^3},& 1+n_1\≤l\≤n\\ \sqrt[3]{\underset{i=1}{\overset{n_1}{\Σ}}{h}_{1i}^3+\underset{j=1}{\overset{n_2}{\Σ}}{h}_{2j}^3+\underset{k=1}{\overset{l-n}{\Σ}}{h}_{Ak}^3},& 1+n\≤l\≤N\end{array}\right.;$

其中，第一级主簧的根部重叠部分的等效厚度第一级主簧与第二级主簧的根部重叠部分的等效厚度h_M2e＝h_ne，主副簧的根部重叠部分的总等效厚度h_MAe＝h_Ne；

(2)两级主簧式非等偏频渐变刚度板簧的各级板簧初始曲率半径的计算：

I步骤：第一级主簧末片下表面初始曲率半径R_M10b计算

根据第一级主簧片数n₁，第一级主簧各片的厚度h_1i，i＝1,2,…,n₁；第一级主簧首片的一半夹紧长度L₁₁，第一级主簧的初始切线弧高H_gM10，对第一级主簧末片下表面初始曲率半径R_M10b进行计算，即

$R_{M10b}=\frac{L_{11}^2+H_{gM10}^2}{2H_{gM10}}+\underset{i=1}{\overset{n_1}{\Σ}}{h}_{1i};$

II步骤：第二级主簧首片上表面初始曲率半径R_M20a计算

根据第二级主簧首片的一半夹紧长度L₂₁，第二级主簧的初始切线弧高H_gM20，对第二级主簧末片上表面初始曲率半径R_M20a进行计算，即

$R_{M20a}=\frac{L_{21}^2+H_{gM20}^2}{2H_{gM20}};$

III步骤：第二级主簧末片下表面初始曲率半径R_M20b计算

根据第二级主簧片数n₂，第二级主簧各片的厚度h_2j，j＝1,2,…,n₂；II步骤中计算得到的R_M20a，对第二级主簧末片下表面初始曲率半径R_M20b进行计算，即

$R_{M20b}=R_{M20a}+\underset{j=1}{\overset{n_2}{\Σ}}{h}_{2j};$

IV步骤：副簧首片上表面初始曲率半径R_A0a计算

根据副簧首片的一半夹紧长度L_A1，副簧的初始切线弧高H_gA0，对副簧末片上表面初始曲率半径R_A0a进行计算，即

$R_{A0a}=\frac{L_{A1}^2+H_{gA0}^2}{2H_{gA0}};$

(3)两级主簧式非等偏频型渐变刚度板簧的各次接触载荷P_k1、P_k2、P_w2的仿真计算：

A步骤：第1次开始接触载荷P_k1的仿真计算

根据两级主簧式非等偏频型渐变刚度板簧的宽度b，弹性模量E；第一级主簧片数n₁，第一级主簧各片的厚度h_1i，i＝1,2,…,n₁，第一级主簧首片的一半夹紧跨长度L₁₁，步骤(1)中计算得到的h_M1e，步骤(2)中计算得到的R_M10b和R_M20a，对第1次开始接触载荷P_k1进行仿真计算，即

$P_{k1}=\frac{{\mathrm{Ebh}}_{M1e}^3(R_{M20a}-R_{M10b})}{6L_{11}{R}_{M20b}{R}_{M10a}};$

B步骤：第2次开始接触载荷P_k2的仿真计算

根据两级主簧式非等偏频型渐变刚度板簧的宽度b，弹性模量E；第一级主簧首片的一半夹紧跨长度L₁₁；第二级主簧片数n₂，第二级主簧各片的厚度h_2j，j＝1,2,…,n₂；步骤(1)中计算得到的h_M2e，步骤(2)中计算得到的R_M20b和R_A0a，及A步骤中仿真计算得到的P_k1，对第2次开始接触载荷P_k2进行仿真计算，即

$P_{k2}=P_{k1}+\frac{{\mathrm{Ebh}}_{M2e}^3(R_{A0a}-R_{M20b})}{6L_{11}{R}_{M20b}{R}_{A0a}};$

C步骤：第2次完全接触载荷P_w2的仿真计算

根据A步骤中仿真计算得到的P_k1，B步骤中仿真计算得到的P_k2，对第2次完全接触载荷P_w2进行验算，即

$P_{w2}=\frac{P_{k2}^2}{P_{k1}};$

(4)两级主簧的渐变刚度板簧的各级夹紧刚度的仿真计算：

a步骤：第一级主簧夹紧刚度K_M1的仿真计算：

根据两级主簧式非等偏频型渐变刚度板簧的宽度b，弹性模量E；第一级主簧片数n₁,第一级主簧各片主簧的一半夹紧长度L_1i＝L_i，i＝1,2,…,n₁，及步骤(1)中计算得到的h_le，l＝i＝1,2,…,n₁，对第一级主簧夹紧刚度K_M1进行仿真计算，即

$K_{M1}=\frac{bE}{2\[\frac{{(L_1-L_2)}^3}{h_{1e}^3}+\underset{l=2}{\overset{n_1-1}{\Σ}}\frac{{(L_1-L_{l+1})}^3-{(L_1-L_l)}^3}{h_{le}^3}+\frac{L_1^3-{(L_1-L_{n_1})}^3}{h_{n_{1}e}^3}\]}$

b步骤：第一级主簧与第二级主簧的复合夹紧刚度K_M2的仿真计算：

根据两级主簧式非等偏频型渐变刚度板簧的宽度b，弹性模量E；第一级主簧片数n₁，第一级主簧各片的一半夹紧长度L_1i＝L_i，i＝1,2,…,n₁；第二级主簧片数n₂,第二级主簧各片的一半夹紧长度L_2j＝L_n1+j，j＝1,2,…,n₂；第一级主簧和第二主簧的总片数n＝n₁+n₂，及步骤(1)中计算得到的h_le，l＝1,2,…,n，对第一级主簧与第二级主簧的复合夹紧刚度K_M2进行仿真计算，即

$K_{M2}=\frac{bE}{2\[\frac{{(L_1-L_2)}^3}{h_{1e}^3}+\underset{l=2}{\overset{n-1}{\Σ}}\frac{{(L_1-L_{l+1})}^3-{(L_1-L_l)}^3}{h_{le}^3}+\frac{L_1^3-{(L_1-L_n)}^3}{h_{ne}^3}\]};$

c步骤：主副簧总复合夹紧刚度K_MA的仿真计算：

根据两级主簧式非等偏频型渐变刚度板簧的宽度b，弹性模量E；第一级主簧的片数n₁，第一级主簧各片的一半夹紧长度L_1i＝L_i，i＝1,2,…,n₁；第二级主簧的片数n₂，第二级主簧各片的一半夹紧长度L_2j＝L_n1+j，j＝1,2,…,n₂；副簧片数m；副簧各片的一半夹紧长度分别为L_A1＝L_n+k，k＝1,2,…,m，主副簧的总片数N＝n₁+n₂+m；及步骤(1)中计算得到的h_le，l＝1,2,…,N，对主副簧的总夹紧复合刚度K_MA进行仿真计算，即，即

$K_{MA}=\frac{bE}{2\[\frac{{(L_1-L_2)}^3}{h_{1e}^3}+\underset{l=2}{\overset{N-1}{\Σ}}\frac{{(L_1-L_{l+1})}^3-{(L_1-L_l)}^3}{h_{le}^3}+\frac{L_1^3-{(L_1-L_N)}^3}{h_{Ne}^3}\]};$

(5)两级主簧式非等偏频型渐变刚度板簧的渐变夹紧刚度的仿真计算：

i步骤：第一级渐变夹紧刚度K_kwP1的仿真计算

根据步骤(3)中仿真计算得到的P_k1和P_k2，步骤(4)中仿真计算得到的K_M1和K_M2；对载荷P在[P_k1,P_k2]范围内的第一级渐变夹紧刚度K_kwP1进行仿真计算，即

$K_{kwP1}=\frac{P}{P_{k1}}{K}_{M1}+\frac{P-P_{k1}}{P_{k2}-P_{k1}}(K_{M2}-\frac{P_{k2}}{P_{k1}}{K}_{M1}),P\∈\[P_{k1},P_{k2}\];$

ii步骤：第二级渐变夹紧刚度K_kwP2的仿真计算

根据步骤(3)中仿真计算得到的P_k2和P_w2，步骤(4)中仿真计算得到的K_M2和K_MA；对载荷P在[P_k2,P_w2]范围内的第二级渐变夹紧刚度K_kwP2进行仿真计算，即

$K_{kwP2}=\frac{P}{P_{k2}}{K}_{M2}+\frac{P-P_{k2}}{P_{w2}-P_{k2}}(K_{MA}-\frac{P_{w2}}{P_{k2}}{K}_{M2}),P\∈\[P_{k2},P_{w2}\];$

(6)两级主簧式非等偏频型渐变刚度板簧的夹紧刚度特性的仿真计算：

根据步骤(3)中仿真计算得到的P_k1、P_k2和P_w2；步骤(4)中仿真计算得到的K_M1、K_M2和K_MA，及步骤(5)骤中仿真计算得到的K_kwP1和K_kwP2，对两级主簧式非等偏频型渐变刚度板簧在不同载荷下的夹紧刚度特性进行仿真计算，即

$K=\left\{\begin{array}{cc}{K}_{M1},& 0\&le;P<P_{k1}\\ \frac{P}{P_{k1}}{K}_{M1}+\frac{P-P_{k1}}{P_{k2}-P_{k1}}(K_{M2}-\frac{P_{k2}}{P_{k1}}{K}_{M1}),& P_{k1}\&le;P<P_{k2}\\ \frac{P}{P_{k2}}{K}_{M2}+\frac{P-P_{k2}}{P_{k2}-P_{k1}}(K_{MA}-\frac{P_{k2}}{P_{k2}}{K}_{M2}),& P_{k2}\&le;P<P_{w2}\\ K_{MA},& P_{w2}\&le;P\end{array}\right..$