CN106777789A - 两级主簧式非等偏频型渐变刚度板簧接触载荷的仿真计算法 - 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_1iT-L₀/2，i＝1,2,…,n₁。第二级主簧2的片数为n₂，第二级主簧各片的厚度为h_2j，一半作用长度为L_2jT，一半夹紧长度L_2j＝L_2jT-L₀/2，j＝1,2,…,n₂。副簧3的片数为m，各片副簧的厚度为h_Ak，一半作用长度为L_AkT，一半夹紧长度L_Ak＝L_AkT-L₀/2，k＝1,2,…,m。根据各片第一级、第二级主簧和副簧的结构参数，弹性模量，初始切线弧高，在等效厚度、渐变刚度计算及开始接触时板簧挠度仿真计算的基础上，对两级主簧式非等偏频型渐变刚度板簧的各次接触载荷进行仿真计算。

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

(1)两级主簧式非等偏频型渐变刚度板簧的两级主簧和副簧的曲率半径计算：

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进行计算，即

(2)两级主簧式非等偏频型渐变刚度板簧的第1次开始接触载荷P_k1的仿真计算：

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

式中，h_M1e为第一级主簧的根部重叠部分等效厚度，

(3)两级主簧式非等偏频型渐变刚度板簧的第2次开始接触载荷P_k2的仿真计算：

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

式中，h_M2e为第一级主簧和第二级主簧的根部重叠部分等效厚度，

(4)第2次开始接触时第一级主簧首片的挠度f_Mk2仿真计算：

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

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

b步骤：第2次开始接触时的第一级主簧首片挠度f_Mk2的仿真计算

根据第一级主簧的夹紧刚度K_M1；步骤(2)中仿真计算得到的P_k1，步骤(3)中仿真计算得到的P_k2，及a步骤中仿真计算得到的K_kwP1，对第2次开始接触时第一级主簧首片挠度f_Mk2进行仿真计算，即

(5)第一级主簧和第二级主簧叠加的单片等效厚度h_M12E的计算：

根据第一级主簧片数n₁，第一级主簧各片的厚度h_1i，第一级主簧各片的一半夹紧长度L_1i，

i＝1,2,…,n₁；第二级主簧片数n₂，第二级主簧各片的厚度h_2j，第二级主簧各片的一半夹紧长度L_2j，j＝1,2,…,n₂；对第一级主簧和第二级主簧叠加之后的单片等效厚度h_M12E计算，即

式中，h_1iE为第一级主簧各片的等效厚度，其中，η_1i-1为第一级主簧的第i片板簧相对于首片的缩短率，

h_2jE为第二级主簧各片等效于第一级主簧首片的等效厚度其中η_2j-1为第二级主簧的第j片板簧相对于第一级主簧首片板簧的缩短率，

(6)两级主簧式非等偏频型渐变刚度板簧的第2次完全接触载荷P_w2的仿真计算：

A步骤：第2次开始接触时的第一级主簧切线弧高H_gMk2的仿真计算

根据第一级主簧的初始切线弧高H_gM10，及步骤(4)的步骤中仿真计算得到的f_Mk2，对第2次开始接触时的第一级主簧切线弧高H_gMk2进行仿真计算，即

H_gMk2＝H_gM10-f_Mk2；

B步骤：第2次开始接触时的第二级主簧末片下表面曲率半径R_Mk2b的仿真计算

根据第一级主簧片数n₁，第一级主簧首片的一半夹紧长度L₁₁，第一级主簧各片的厚度h_1i，i＝1,2,…,n₁；第二级主簧的片数n₂，第二级主簧各片的厚度h_2j，j＝1,2,…,n₂；A步骤中计算得到的H_gMk2，对第2次开始接触时的第二级主簧末片下表面曲率半径R_Mk2b进行仿真计算，即

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

根据两级主簧式非等偏频型渐变刚度板簧的宽度b＝63mm，弹性模量E＝200GPa；第一级主簧首片的一半夹紧长度L₁₁；步骤(1)中仿真计算得到的R_A0a，步骤(3)中仿真计算得到的P_k2，步骤(5)中计算得到的h_M12E，及B步骤中仿真计算得到的R_Mk2b，对第2次完全接触载荷P_w2进行仿真计算，即

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

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

附图说明

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

图1是两级主簧式非等偏频型渐变刚度板簧接触载荷的仿真计算流程图；

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

具体实施方案

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

实施例：某两级主簧式非等偏频渐变刚度板簧的宽度b＝63mm，骑马螺栓夹紧距的一半L₀＝50mm，弹性模量E＝200GPa。第一级主簧片数n₁＝2，第一级主簧各片的厚度h₁₁＝h₁₂＝8mm，第一级主簧各片的一半作用长度分别为L_11T＝525mm，L_12T＝450mm；一半夹紧长度分别为L₁₁＝L_11T-L₀/2＝500mm，L₁₂＝L_12T-L₀/2＝425mm。第二级主簧片数n₂＝1，厚度h₂₁＝8mm，一半作用长度L_21T＝350mm，一半夹紧长度L₂₁＝L_21T-L₀/2＝325mm。副簧片数m＝2，副簧各片的厚度h_A1＝h_A2＝13mm；副簧首片的一半作用长度L_A1T＝250mm，一半夹紧长度为L_A1＝L_A1T-L₀/2＝225mm。第一级主簧的初始切线弧高设计值H_gM10＝103.7mm，第二级主簧的初始切线弧高H_gM20＝18.8mm，副簧的初始切线弧高H_gA0＝6mm。根据各片板簧的结构参数，弹性模量，初始切线弧高设计值，对该两级主簧式非等偏频渐变刚度板簧的各次接触载荷进行仿真计算，确保板簧初始切线弧高满足接触载荷的设计要求。

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

(1)两级主簧式非等偏频型渐变刚度板簧的两级主簧和副簧的曲率半径计算：

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

根据第一级主簧片数n₁＝2，第一级主簧各片的厚度h₁₁＝h₁₂＝8mm，第一级主簧首片的一半夹紧长度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进行计算，即

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

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

(2)两级主簧式非等偏频型渐变刚度板簧的第1次开始接触载荷P_k1的仿真计算：

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

式中，h_M1e为第一级主簧的根部重叠部分等效厚度，

(3)两级主簧式非等偏频型渐变刚度板簧的第2次开始接触载荷P_k2的仿真计算：

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

式中，h_M2e为第一级主簧和第二级主簧的根部重叠部分等效厚度，

(4)第2次开始接触时第一级主簧首片的挠度f_Mk2仿真计算：

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

根据第一级主簧的夹紧刚度K_M1＝51.43N/mm；第一级与第二级主簧的复合夹紧刚度K_M2＝75.4N/mm，步骤(2)中仿真计算得到的P_k1＝1851.3N，步骤(3)中仿真计算得到的P_k2＝2606.2N，对载荷P在[P_k1,P_k2]范围内的第一级渐变夹紧刚度K_kwP1进行仿真计算，即

b步骤：第2次开始接触时的第一级主簧首片挠度f_Mk2的仿真计算

根据第一级主簧的夹紧刚度K_M1＝51.43N/mm；步骤(2)中仿真计算得到的P_k1＝1851.3N，步骤(3)中仿真计算得到的P_k2＝2606.2N，及a步骤中仿真计算得到的K_kwP，对第2次开始接触时第一级主簧首片挠度f_Mk2进行仿真计算，即

(5)第一级主簧和第二级主簧的叠加单片等效厚度h_M12E的计算：

根据第一级主簧片数n₁＝2，第一级主簧各片的厚度h₁₁＝h₁₂＝8mm，第一级主簧各片的一半夹紧长度L₁₁＝500mm，L₁₂＝425mm；第二级主簧片数n₂＝1，厚度h₂₁＝8mm，一半夹紧长度L₂₁＝325mm；对第一级主簧和第二级主簧叠加之后的单片等效厚度h_M12E计算，即

式中，h_11E和h_12E分别第一级主簧各片的等效厚度，h_11E＝h₁₁＝8mm，其中，η_12-1为第一级主簧的第2片相对于首片的缩短率，

h_21E为第二级主簧等效于第一级主簧首片的等效厚度其中η_21-1为第二级主簧的一片板簧，相对于第一级主簧首片板簧的缩短率，(6)两级主簧式非等偏频型渐变刚度板簧的第2次完全接触载荷P_w2的仿真计算：

A步骤：第2次开始接触时的第一级主簧切线弧高H_gMk2的仿真计算

根据第一级主簧的初始切线弧高H_gM10＝103.7mm，及b步骤中仿真计算得到的f_Mk2＝48mm，对第2次开始接触时的第一级主簧的切线弧高H_gMk2进行仿真计算，即

H_gMk2＝H_gM10-f_Mk2＝55.7mm；

B步骤：第2次开始接触时的第二级主簧末片下表面曲率半径R_Mk2b的仿真计算

根据第一级主簧片数n₁＝2，第一级主簧首片的一半夹紧长度L₁₁＝500mm，第一级主簧各片的厚度h₁₁＝h₁₂＝8mm；第二级主簧的片数n₂＝1，厚度h₂₁＝8mm，A步骤中计算得到的H_gMk2＝55.7mm，对第2次开始接触时的第二级主簧末片下表面曲率半径R_Mk2b进行仿真计算，即

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

根据两级主簧的渐变刚度钢板弹簧的宽度b＝63mm，弹性模量E＝200GPa；第一级主簧首片的一半夹紧长度L₁₁＝500mm；步骤(1)中仿真计算得到的R_A0a＝4221.8mm，步骤(3)中仿真计算得到的P_k2＝2606.2N，步骤(5)中计算得到的h_M12E＝11.42mm，及B步骤中仿真计算得到的R_Mk2b＝2297.7mm，对第2次完全接触载荷P_w2进行仿真计算，即

通过接触载荷的仿真计算值与设计要求值，比较可知：该两级主簧的渐变刚度钢板弹簧接触载荷的仿真计算值与设计要求值相吻合，说明该板簧的设计结构及初始切线弧高设计值是可靠的，能满足板簧接触载荷设计要求。

通过样机加载挠度试验测试可知，本发明所提供的两级主簧式非等偏频型渐变刚度板簧接触载荷的仿真计算法是正确的，为两级主簧式非等偏频型渐变刚度板簧的特性验证及CAD软件开发奠定了可靠的技术基础。利用该方法可得到可靠的接触载荷仿真计算值，确保板簧的设计结构满足接触载荷、渐变刚度及应力强度的设计要求，提高产品设计水平、质量和性能及车辆行驶平顺性；同时，降低设计及试验费用，加快产品开发速度。

Claims (1)

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

(1)两级主簧式非等偏频型渐变刚度板簧的两级主簧和副簧的曲率半径计算：

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}};$

(2)两级主簧式非等偏频型渐变刚度板簧的第1次开始接触载荷P_k1的仿真计算：

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

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

式中，h_M1e为第一级主簧的根部重叠部分等效厚度，

(3)两级主簧式非等偏频型渐变刚度板簧的第2次开始接触载荷P_k2的仿真计算：

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

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

式中，h_M2e为第一级主簧和第二级主簧的根部重叠部分等效厚度，

(4)第2次开始接触时第一级主簧首片的挠度f_Mk2仿真计算：

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

根据第一级主簧的夹紧刚度K_M1；第一级与第二级主簧的复合夹紧刚度K_M2，步骤(2)中仿真计算得到的P_k1，步骤(3)中仿真计算得到的P_k2，对载荷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}\];$

b步骤：第2次开始接触时的第一级主簧首片挠度f_Mk2的仿真计算

根据第一级主簧的夹紧刚度K_M1；步骤(2)中仿真计算得到的P_k1，步骤(3)中仿真计算得到的P_k2，及a步骤中仿真计算得到的K_kwP1，对第2次开始接触时第一级主簧首片挠度f_Mk2进行仿真计算，即

$f_{Mk2}=\frac{P_{k1}}{K_{M1}}+{\∫}_{P_{k1}}^{P_{k2}}\frac{dP}{K_{kwP1}};$

(5)第一级主簧和第二级主簧叠加的单片等效厚度h_M12E的计算：

根据第一级主簧片数n₁，第一级主簧各片的厚度h_1i，第一级主簧各片的一半夹紧长度L_1i，i＝1,2,…,n₁；第二级主簧片数n₂，第二级主簧各片的厚度h_2j，第二级主簧各片的一半夹紧长度L_2j，j＝1,2,…,n₂；对第一级主簧和第二级主簧叠加之后的单片等效厚度h_M12E计算，即

$h_{M12E}=\sqrt[3]{\underset{i=1}{\overset{n_1}{\Σ}}{h}_{1iE}^3+\underset{j=1}{\overset{n_2}{\Σ}}{h}_{2jE}^3};$

式中，h_1iE为第一级主簧各片的等效厚度，其中，η_1i-1为第一级主簧的第i片板簧相对于首片的缩短率，

h_2jE为第二级主簧各片等效于第一级主簧首片的等效厚度其中η_2j-1为第二级主簧的第j片板簧相对于第一级主簧首片板簧的缩短率，j＝1,2,…,n₂；

(6)两级主簧式非等偏频型渐变刚度板簧的第2次完全接触载荷P_w2的仿真计算：

A步骤：第2次开始接触时的第一级主簧切线弧高H_gMk2的仿真计算

根据第一级主簧的初始切线弧高H_gM10，及步骤(4)的步骤中仿真计算得到的f_Mk2，对第2次开始接触时的第一级主簧切线弧高H_gMk2进行仿真计算，即

H_gMk2＝H_gM10-f_Mk2；

B步骤：第2次开始接触时的第二级主簧末片下表面曲率半径R_Mk2b的仿真计算

根据第一级主簧片数n₁，第一级主簧首片的一半夹紧长度L₁₁，第一级主簧各片的厚度h_1i，i＝1,2,…,n₁；第二级主簧的片数n₂，第二级主簧各片的厚度h_2j，j＝1,2,…,n₂；A步骤中计算得到的H_gMk2，对第2次开始接触时的第二级主簧末片下表面曲率半径R_Mk2b进行仿真计算，即

$R_{Mk2b}=\frac{L_{11}^2+{H^2}_{gMk2}}{2H_{gMk2}}+\underset{i=1}{\overset{n_1}{\mathrm{\Σ}}}{h}_{1i}+\underset{j=1}{\overset{n_2}{\Σ}}{h}_{2j};$

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

根据两级主簧式非等偏频型渐变刚度板簧的宽度b＝63mm，弹性模量E＝200GPa；第一级主簧首片的一半夹紧长度L₁₁；步骤(1)中仿真计算得到的R_A0a，步骤(3)中仿真计算得到的P_k2，步骤(5)中计算得到的h_M12E，及B步骤中仿真计算得到的R_Mk2b，对第2次完全接触载荷P_w2进行仿真计算，即

$P_{w2}=P_{k2}+\frac{{\mathrm{Ebh}}_{M12E}^3}{6L_{11}}(\frac{1}{R_{Mk2b}}-\frac{1}{R_{A0}}).$