CN109840387A - 一种主副簧式复合材料板弹簧复合刚度计算方法 - Google Patents

Abstract

本发明涉及一种主副簧式复合材料板弹簧复合刚度计算方法。该计算方法首先根据复合材料力学计算复合材料主簧和副簧各横截面的弯曲刚度，通过公式推导给出了复合材料主簧和副簧各横截面弯曲刚度的计算公式。然后，基于有限差分理论，给出了复合材料主副簧共同工作后复合刚度的计算步骤，并推导了主副簧式复合材料板弹簧复合刚度的计算公式。本发明提出的计算方法不但为主副簧式复合材料板弹簧复合刚度计算问题提供了理论依据，而且计算精度高，适用于具有任意截面形状的主副簧式复合材料板弹簧复合刚度计算问题。此外，该计算方法也适宜编程计算，计算速度远高于有限元分析等传统计算方法，能显著缩短复合材料板簧的研发周期，降低研发成本。

Description

一种主副簧式复合材料板弹簧复合刚度计算方法

【技术领域】

本发明涉及一种主副簧式复合材料板弹簧复合刚度计算方法，其应用于汽车悬架，属于板弹簧技术领域。

【背景技术】

板弹簧是汽车悬架中广泛应用的弹性元件，其重量占汽车非簧载质量的10％-20％。复合材料板弹簧是一种采用纤维增强树脂基复合材料的板弹簧。在功能相同的前提下，用复合材料板弹簧替换钢板弹簧可减重50％以上。此外，复合材料板弹簧也具有比模量、比强度高，疲劳性能好，可设计性强，耐腐蚀，易于加工，安全冗余高等显著的优点。因此，复合材料板弹簧的综合性能明显优于钢板弹簧，是汽车轻量化领域的研究热点。目前，国内外已对单片复合材料板弹簧的接头可靠性，结构设计，铺层设计，刚度计算，性能优化及性能测试方面进行了深入研究。然而，汽车主要分为满载和空载这两种状态，对应要求板簧提供数值差别较大的两级刚度，但现有的单片复合材料板簧在固化成型后刚度不可调，只能达到某种折衷下的较优减振效果，无法使悬架刚度最优。因此，现有的单片复合材料板簧无法满足日益增长的汽车性能需求。而采用与主副簧式钢板弹簧结构类似的主副簧式复合材料板弹簧，不但能起到显著的轻量化效果，而且能兼顾汽车性能需求，是实现复合材料板簧变刚度特性的有效途径之一。

目前，对多片复合材料板弹簧的研究主要集中于有限元仿真、接头设计及性能测试方面，现有设计理论仍然无法满足主簧-副簧式复合材料板弹簧的应用需要。刚度特性是主副簧式复合材料板弹簧的核心性能参数，不但直接决定了汽车的操纵稳定性及平顺性，也是其区别于单片复合材料板簧的重要结构特征，进而决定了其应用价值。准确的刚度计算，不但能有效预测样件刚度，而且能缩短产品开发的周期，进而压缩复合材料板弹簧的产品开发成本。然而，由于复合材料具有各向异性，复合材料板弹簧也具有变截面几何特征，且主副簧共同工作后的接触边界有非线性因素，这都导致主副簧式复合材料板弹簧的复合刚度计算变得非常棘手和复杂。

因此，为解决上述技术问题，确有必要提供一种创新的主副簧式复合材料板弹簧复合刚度计算方法，以克服现有技术中的所述缺陷。

【发明内容】

为解决上述问题，本发明的目的在于提供一种主副簧式复合材料板弹簧复合刚度计算方法。

为实现上述目的，本发明采取的技术方案为：一种主副簧式复合材料板弹簧复合刚度计算方法，其包括如下步骤：

1)，计算复合材料主簧和副簧各横截面的弯曲刚度；

2)，计算复合材料主副簧共同工作后的复合刚度；具体为：假设主簧和副簧在主簧D点接触，作用力及反作用力分别为F_D和F_D′，主副簧总成在主簧接头中心O点受到竖直向下的载荷F_O，则复合材料主副簧共同工作后的复合刚度通过以下方法来计算：

2-1，假设只有主簧受到力F_O的作用，计算D点的挠度及O点的挠度；

2-2，假设主簧D点单独受到作力F_D作用，计算末端O点的挠度；

2-3，对副簧GI段使用有限差分法，求得截面I的挠度；

2-4，求出F_D的具体数值。

2-5，求出主副簧总成后半段的复合刚度；

2-6，计算出主副簧总成的复合刚度。

本发明的主副簧式复合材料板弹簧复合刚度计算方法进一步为：所述步骤1)具体为：对复合材料板弹簧的任意一小段，设复合材料板弹簧的纵向为x轴，厚度方向为z轴，且向上为正，该小段的弯曲中性轴为y轴；假定复合材料层合结构在弯曲变形时满足平面假设，即变形前相距dx的两个横截平面，弯曲变形后仍保持平面，仅相对旋转了一个角度dθ，

根据应变的定义，第i层单层沿x轴的应变为：

对复合材料板弹簧任意截面，设其第k层与该截面的几何中面的距离为

式(2)中，λ是中性层位移系数，λ＝d/t，d为中性层与几何中面之间的距离向量；t为层合结构在该截面的厚度；δ为单层厚度；

对复合材料层合结构的任意一个横截面，沿x轴方向的力可表示为：

式(3)中，A^t为截面受拉区域面积，A^c为截面受压区域面积；

不考虑各单层之间的层间应力，且假设各单层的宽度相等，沿z轴对层合结构分层积分，将式(3)化为：

其中，b为层合结构的宽度，n₁为截面受拉区域层数；n₂为截面受压区域层数，为截面中第k层沿x轴的应力；为截面受拉区域第k₁层沿x轴的应力；为截面受压区域第k₂层沿x轴的应力；

因弯曲半径ρ不随单层距中性层的距离而变化，且在各单层内模量保持不变，故将式(4)化为：

因复合材料层合结构受纯弯矩的作用，故有：

F_N＝0 (6)

进而有：

联立式(2)和式(7)，即可解得中性层位移系数λ；

对于复合材料板弹簧任意一个纵向截面，其受到绕y轴的弯矩可表示为：

设截面抗弯刚度为K，则：

即有：

上式即为复合材料主簧和副簧变形部分各横截面弯曲刚度的计算公式。

本发明的主副簧式复合材料板弹簧复合刚度计算方法进一步为：所述步骤2-1)具体为：把BE段分成m+1段，则形成m个新截面，加上固定端面BG，共有(m+1)个截面；设每段长度为a，则有：

对BE段，根据有限差分理论，可得到由(m+1)个方程组成的方程组：

其中M_i＝-F_Ox＝-F_O(L_r-S_r-ia)(i＝0,1,2,...,m)；

在i＝0处，即BG截面的转角和挠度均为0，则由转角表达式得f_-1＝f₁又f₀＝0，所以由方程组的第一个方程可得代入第二个方程可得f₂＝a²(M₀/K₀+M₁/K₁)，依次代入下去可得到第i个截面处的挠度表达式为：

当i取(m+1)时即可求得截面E的挠度

再根据式(13)计算出f_m-1，然后把f_E和f_m-1代入式计算出θ_m；

则截面E的转角：

则F_O单独作用下主簧末端的挠度为：

f_O＝f_O1+f_E+f_Eθ＝f_E+θ_E×l₁ (16)

其中，f_O1为假设E处固定，O端作用载荷F_O时，O端的挠度；由于EO段被金属接头覆盖，视为刚性体，因此f_O1＝0

再求F_O单独作用下主簧截面D的挠度：由S_r+ia＝L_r-l₁-l₂得当时，对应i₁截面的挠度即为截面D的挠度，

则截面D的挠度为：

其中i₁为计算数值取整。

本发明的主副簧式复合材料板弹簧复合刚度计算方法进一步为：所述步骤2-3)具体为：使用有限差分法，将BD段分成(n+1)段，就有n个新截面，则有：

加上固定端面，共有(n+1)个截面，这样就可以得到由(n+1)个方程组成的方程组：

其中M_i′＝F_D(L_r-S_r-l₁-l₂-ia′)(i＝0,1,2,...,n)；

在i＝0处，即B端面处转角和挠度为0，则由转角表达式得f_-1＝f₁又f₀＝0，所以由方程组的第一个方程可以得到：代入第二个方程可以得到f₂′＝a′²(M₀′/K₀′+M₁′/K₁′)，依次代入下去可以得到第i截面处的挠度表达式为

当i取(n+1)时即可求得主簧D处只受到力F_D作用时截面D的挠度

由式(20)计算出f_n′-₁，把f_D′和f_n′-₁代入式计算出截面D的转角θ′_D：

当主簧受到力F_D作用时，DO段无外部集中载荷作用，且自身重力可忽略不计，因此当主簧只在D处受到力F_D作用时，其末端O点的挠度为：

f_O′＝f_D′+f_D′_θ＝f_D′+θ′_D×(l₁+l₂) (22)

其中，f_D′_θ是D处弯曲转角使O端产生的位移。

本发明的主副簧式复合材料板弹簧复合刚度计算方法进一步为：所述步骤2-3)具体为：把GI段分成(p+1)段，则形成p个新截面，加上固定端面BG，共有(p+1)个截面；设每段长度为a″，则有：

对GI段，根据有限差分理论，可得到由(p+1)个方程组成的方程组：

其中M_i″＝-F_D′x′＝-F_D(L_r-S_r-l₁-l₂-ia″)(i＝0,1,2,...,p)

同理，可得到第i个截面处的挠度表达式为

当i取(p+1)时即可求得截面I的挠度

本发明的主副簧式复合材料板弹簧复合刚度计算方法进一步为：所述步骤2-4)具体为：在主簧与副簧共同工作后，主簧D点的挠度为副簧末端对D点的集中力F_D与F_O的共同作用产生的，即主簧D点的挠度为f_D″＝f_D+f_D′；而主簧和副簧在接触点D处的挠度相同，因此f_D″＝f_I；根据该式可求出

其中，F_O为已知量；m，n，p为程序中取的各截面数；K_i，K_i′和K_i″分别是各板簧簧身的截面弯曲刚度。

本发明的主副簧式复合材料板弹簧复合刚度计算方法进一步为：所述步骤2-5)具体为：主簧O点的挠度f_O″是集中力F_D与F_O的共同作用产生的，即

f_O″＝f_O+f_O′＝(f_E+θ_E×l₁)+[f_D′+θ′_D×(l₁+l₂)] (28)

其中，f_E及θ_E与F_D无关，f_D′及θ′_D与F_D相关；将求出的F_D代入上式即可求出主簧和副簧共同工作后主簧O点的挠度f_O″，则主副簧总成后半段的复合刚度为：

本发明的主副簧式复合材料板弹簧复合刚度计算方法还为：所述步骤2-5)具体为：通过下式计算出主副簧总成的复合刚度：

其中，K_f和K_r分别为总成前半段和后半段的复合刚度，L_f和L_r分别为总成前后半段半长。

与现有技术相比，本发明具有如下有益效果：

(1)本发明的主副簧式复合材料板弹簧复合刚度计算方法不但为主副簧式复合材料板弹簧复合刚度计算问题提供了理论依据，而且其正确性已经通过了仿真及试验验证，计算误差在3％以内，因此具有较高的计算精度，完全满足工程应用的要求。

(2)由于该计算方法引入了微积分思想，因此适用于具有任意截面形状的主副簧式复合材料板弹簧复合刚度的计算，具有普遍性适用性。

(3)由于该计算方法适宜编程计算，因此计算速度远高于有限元分析等传统计算方法，这将显著缩短复合材料板簧的研发周期，降低研发成本。

【附图说明】

图1为复合材料层合板发生纯弯曲变形前后的状态示意图。

图2为某主副簧式复合材料板弹簧后半段的结构及相关几何尺寸示意图。

【具体实施方式】

本发明为一种主副簧式复合材料板弹簧复合刚度计算方法，其包括如下步骤：

1.计算复合材料主簧和副簧各横截面的弯曲刚度

复合材料主簧和副簧本质上是具有一定结构特征的层合结构。复合材料层合结构在弯曲变形时，由于各单层的拉伸模量E_t与压缩模量E_c不等，将导致弯曲中性层在层合结构厚度方向上偏离几何中面一定距离，这对复合材料板弹簧刚度计算精度有很大影响。因此，应首先确定复合材料主簧和副簧中的弯曲中性层的位置。对复合材料板弹簧的任意一小段，设复合材料板弹簧的纵向为x轴，厚度方向为z轴，且向上为正，该小段的弯曲中性轴为y轴。当复合材料板弹簧发生垂向变形时，该小段近似受到纯弯矩作用。假定复合材料层合结构在弯曲变形时满足平面假设，即变形前相距dx的两个横截平面，弯曲变形后仍保持平面，仅相对旋转了一个角度dθ，如图1所示。

在图1中，根据应变的定义，第i层单层沿x轴的应变为：

根据式(1)，复合材料层合结构各单层沿x轴的应变大小与单层距中性层的距离z_i成正比。

对复合材料板弹簧任意截面，设其第k层与该截面的几何中面的距离为

式(2)中，λ是中性层位移系数，λ＝d/t，d为中性层与几何中面之间的距离向量；t为层合结构在该截面的厚度；δ为单层厚度。

对于复合材料板簧这一纵向维度远大于横向维度的部件，各单层横向应变ε_y及剪切应变γ_xy相对于纵向应变ε_x来说可忽略不计，因此忽略纵向应变与其它各向的泊松耦合效应。因复合材料的拉伸模量E_t与压缩模量E_c不等，因此对截面内受拉区域和受压区域分别积分，

对复合材料层合结构的任意一个横截面，沿x轴方向的力可表示为：

式(3)中，A^t为截面受拉区域面积，A^c为截面受压区域面积，且有

不考虑各单层之间的层间应力，且假设各单层的宽度相等，沿z轴对层合结构分层积分，将式(3)化为：

其中，b为层合结构的宽度，n₁为截面受拉区域层数；n₂为截面受压区域层数，为截面中第k层沿x轴的应力；为截面受拉区域第k₁层沿x轴的应力；为截面受压区域第k₂层沿x轴的应力。

因弯曲半径ρ不随单层距中性层的距离而变化，且在各单层内模量保持不变，故将式(4)化为：

因复合材料层合结构受纯弯矩的作用，故有：

F_N＝0 (6)

进而有：

联立式(2)和式(7)，即可解得中性层位移系数λ。

对于复合材料板弹簧任意一个纵向截面，其受到绕y轴的弯矩可表示为：

设截面抗弯刚度为K，则：

即有：

上式即为复合材料主簧和副簧变形部分各横截面弯曲刚度的计算公式，将在后面的计算过程中调用。

2.计算复合材料主副簧共同工作后的复合刚度

在完成复合材料主副簧横截面弯曲刚度的计算后，主副簧各横截面的弯曲刚度K_i已知。以复合材料主副簧后半段为例，主副簧的结构及相关参数如图2所示。板簧宽度特征既可以是等宽的，也可以是任意形式的。其中，主簧和副簧在主簧D点接触，作用力及反作用力分别为F_D和F_D′，主副簧总成在主簧接头中心O点受到竖直向下的载荷F_O。则复合材料主副簧共同工作后的复合刚度通过以下方法来计算：

(2-1)假设只有主簧受到力F_O的作用(即忽略F_D对主簧变形的影响)，计算D点的挠度及O点的挠度

在该状态下，AB段实际上受到悬架骑马螺栓的固定，其挠度及转角忽略不计，因此将BG截面视为主副簧的固定端。把BE段(包括抛物线段BC及等厚段CE)分成m+1段，则形成m个新截面，加上固定端面BG，共有(m+1)个截面。设每段长度为a，则有：

对BE段，根据有限差分理论，可得到由(m+1)个方程组成的方程组：

其中M_i＝-F_Ox＝-F_O(L_r-S_r-ia)(i＝0,1,2,...,m)(注意此式的x原点为O)

在i＝0处，即BG截面的转角和挠度均为0，则由转角表达式得f_-1＝f₁又f₀＝0，所以由方程组的第一个方程可得代入第二个方程可得f₂＝a²(M₀/K₀+M₁/K₁)，依次代入下去可得到第i个截面处的挠度表达式为：

当i取(m+1)时即可求得截面E的挠度

再根据式(13)计算出f_m-1，然后把f_E和f_m-1代入式计算出θ_m。则截面E的转角：

则F_O单独作用下主簧末端的挠度为：

f_O＝f_O1+f_E+f_Eθ＝f_E+θ_E×l₁ (16)

其中，f_O1为假设E处固定，O端作用载荷F_O时，O端的挠度。由于EO段被金属接头覆盖，视为刚性体，因此f_O1＝0。f_Eθ是E处的弯曲转角使O端产生的位移。

则可根据式(16)求出f_O。下面求F_O单独作用下主簧截面D的挠度：由于m趋于无穷大，因此可将主簧截面D的挠度看作取在D附近最近截面对应的挠度。由S_r+ia＝L_r-l₁-l₂得当时，对应i₁(将求得的i数值取整)截面的挠度即为截面D的挠度，则截面D的挠度为：

其中i₁为计算数值取整(最接近的整数)。

(2-2)假设主簧D点单独受到作力F_D作用，则对BD段，再次使用有限差分法，将BD段(包括BC段变厚度段和CD段等厚段)分成(n+1)段，就有n个新截面，则有：

加上固定端面，共有(n+1)个截面，这样就可以得到由(n+1)个方程组成的方程组：

其中M_i′＝F_D(L_r-S_r-l₁-l₂-ia′)(i＝0,1,2,...,n)。

在i＝0处，即B端面处转角和挠度为0，则由转角表达式得f_-1＝f₁又f₀＝0，所以由方程组的第一个方程可以得到：代入第二个方程可以得到f₂′＝a′²(M₀′/K₀′+M₁′/K₁′)，依次代入下去可以得到第i截面处的挠度表达式为

当i取(n+1)时即可求得主簧D处只受到力F_D作用时截面D的挠度

由式(20)计算出f_n′-₁，把f_D′和f_n′-₁代入式计算出截面D的转角θ′_D：

当主簧受到力F_D作用时，DO段无外部集中载荷作用，且自身重力可忽略不计，因此当主簧只在D处受到力F_D作用时，其末端O点的挠度为：

f_O′＝f_D′+f_D′_θ＝f_D′+θ′_D×(l₁+l₂) (22)

其中，f_D′_θ是D处弯曲转角使O端产生的位移。

(2-3)对副簧GI段再次使用有限差分法。把GI段(包括抛物线段GH及等厚段HI，I为副簧片段降噪垫块的对称中心)分成(p+1)段，则形成p个新截面，加上固定端面BG，共有(p+1)个截面。设每段长度为a″，则有：

对GI段，根据有限差分理论，可得到由(p+1)个方程组成的方程组：

其中M_i″＝-F_D′x′＝-F_D(L_r-S_r-l₁-l₂-ia″)(i＝0,1,2,...,p)

同理，可得到第i个截面处的挠度表达式为

当i取(p+1)时即可求得截面I的挠度

(2-4)实际上在主簧与副簧共同工作后，主簧D点的挠度为副簧末端对D点的集中力F_D与F_O的共同作用产生的，即主簧D点的挠度为f_D″＝f_D+f_D′。而主簧和副簧在接触点D处的挠度相同，因此f_D″＝f_I。根据该式可求出

其中，各几何参数均为已知量，F_O为已知量；m，n，p为程序中取的各截面数，数值越大计算精度越高；K_i，K_i′和K_i″分别是各计算问题中板簧簧身的截面弯曲刚度，在板簧整套参数已知的情况下，这些弯曲刚度值是已知量。则可通过上式求出F_D的具体数值。

(2-5)实际上在主簧和副簧共同工作后，主簧O点的挠度f_O″是集中力F_D与F_O的共同作用产生的，即

f_O″＝f_O+f_O′＝(f_E+θ_E×l₁)+[f_D′+θ′_D×(l₁+l₂)] (28)

其中，f_E及θ_E与F_D无关，f_D′及θ′_D与F_D相关。将求出的F_D代入上式即可求出主簧和副簧共同工作后主簧O点的挠度f_O″，则主副簧总成后半段的复合刚度为：

(2-6)同样地，可通过上述理论计算出总成前半段主副簧复合刚度，通过下式计算出主副簧总成的复合刚度：

其中，K_f和K_r分别为总成前半段和后半段的复合刚度，L_f和L_r分别为总成前后半段半长。

因此，在主副簧式复合材料板弹簧的几何尺寸参数及铺层参数已知的情况下，可通过上述计算方法编程计算出主副簧式复合材料板弹簧复合刚度。

实施例

如图2为主副簧式复合材料板弹簧后半段的结构及相关几何尺寸示意图。在本例中，各参数的具体值为：

主簧：L_r＝715mm，板簧宽度b＝70mm，H＝31mm，l_r＝687mm，l₃＝50mm，l₂＝58mm，l₁＝92mm，h＝17mm,抛物线段BC厚度的变化规律：

副簧：H′＝31mm，l₅＝20mm，副簧总长l₄＝100mm，h′＝13.5mm，抛物线段GH厚度的变化规律：

主副簧式复合材料板弹簧前半段的结构及相关几何尺寸与后半段基本相同，不同的参数包括：主簧前半段长度为700mm，副簧前半段长度为：560mm。

根据本发明给出的计算公式，在MATLAB软件中编写相应的计算程序，将上述已知参数输入计算程序中，并设置合适的截面数，提交MATLAB软件计算。计算得到该主副簧式复合材料板弹簧总成的复合刚度计算值为130.54N/mm,与实际刚度之间的误差为2.6％。

以上的具体实施方式仅为本创作的典型实施例，并不用以限制本创作，凡在本创作的精神及原则之内所做的任何修改、等同替换、改进等，均应包含在本创作的保护范围之内。

Claims (8)

1.一种主副簧式复合材料板弹簧复合刚度计算方法，其特征在于：包括如下步骤：

1)，计算计算复合材料主簧和副簧各横截面的弯曲刚度；

2)，计算复合材料主副簧共同工作后的复合刚度；具体为：假设主簧和副簧在主簧D点接触，作用力及反作用力分别为F_D和F′_D，主副簧总成在主簧接头中心O点受到竖直向下的载荷F_O，则复合材料主副簧共同工作后的复合刚度通过以下方法来计算：

2-1，假设只有主簧受到力F_O的作用，计算D点的挠度及O点的挠度；

2-2，假设主簧D点单独受到作力F_D作用，计算末端O点的挠度；

2-3，对副簧GI段使用有限差分法，求得截面I的挠度；

2-4，求出F_D的具体数值。

2-5，求出主副簧总成后半段的复合刚度；

2-6，计算出主副簧总成的复合刚度。

2.如权利要求1所述的主副簧式复合材料板弹簧复合刚度计算方法，其特征在于：所述步骤1)具体为：对复合材料板弹簧的任意一小段，设复合材料板弹簧的纵向为x轴，厚度方向为z轴，且向上为正，该小段的弯曲中性轴为y轴；假定复合材料层合结构在弯曲变形时满足平面假设，即变形前相距dx的两个横截平面，弯曲变形后仍保持平面，仅相对旋转了一个角度dθ，

根据应变的定义，第i层单层沿x轴的应变为：

对复合材料板弹簧任意截面，设其第k层与该截面的几何中面的距离为

式(2)中，λ是中性层位移系数，λ＝d/t，d为中性层与几何中面之间的距离向量；t为层合结构在该截面的厚度；δ为单层厚度；

对复合材料层合结构的任意一个横截面，沿x轴方向的力可表示为：

式(3)中，A^t为截面受拉区域面积，A^c为截面受压区域面积；

不考虑各单层之间的层间应力，且假设各单层的宽度相等，沿z轴对层合结构分层积分，将式(3)化为：

其中，b为层合结构的宽度，n₁为截面受拉区域层数；n₂为截面受压区域层数，为截面中第k层沿x轴的应力；为截面受拉区域第k₁层沿x轴的应力；为截面受压区域第k₂层沿x轴的应力；

因弯曲半径ρ不随单层距中性层的距离而变化，且在各单层内模量保持不变，故将式(4)化为：

因复合材料层合结构受纯弯矩的作用，故有：

F_N＝0 (6)

进而有：

联立式(2)和式(7)，即可解得中性层位移系数λ；

对于复合材料板弹簧任意一个纵向截面，其受到绕y轴的弯矩可表示为：

设截面抗弯刚度为K，则：

即有：

上式即为复合材料主簧和副簧变形部分各横截面弯曲刚度的计算公式。

3.如权利要求1所述的主副簧式复合材料板弹簧复合刚度计算方法，其特征在于：所述步骤2-1)具体为：把BE段分成m+1段，则形成m个新截面，加上固定端面BG，共有(m+1)个截面；设每段长度为a，则有：

对BE段，根据有限差分理论，可得到由(m+1)个方程组成的方程组：

其中M_i＝-F_Ox＝-F_O(L_r-S_r-ia)(i＝0,1,2,...,m)；

在i＝0处，即BG截面的转角和挠度均为0，则由转角表达式得f_-1＝f₁又f₀＝0，所以由方程组的第一个方程可得代入第二个方程可得f₂＝a²(M₀/K₀+M₁/K₁)，依次代入下去可得到第i个截面处的挠度表达式为：

当i取(m+1)时即可求得截面E的挠度

再根据式(13)计算出f_m-1，然后把f_E和f_m-1代入式计算出θ_m；

则截面E的转角：

则F_O单独作用下主簧末端的挠度为：

f_O＝f_O1+f_E+f_Eθ＝f_E+θ_E×l₁ (16)

其中，f_O1为假设E处固定，O端作用载荷F_O时，O端的挠度；由于EO段被金属接头覆盖，视为刚性体，因此f_O1＝0

再求F_O单独作用下主簧截面D的挠度：由S_r+ia＝L_r-l₁-l₂得当时，对应i₁截面的挠度即为截面D的挠度，则截面D的挠度为：

其中i₁为计算数值取整。

4.如权利要求3所述的主副簧式复合材料板弹簧复合刚度计算方法，其特征在于：所述步骤2-3)具体为：使用有限差分法，将BD段分成(n+1)段，就有n个新截面，则有：

加上固定端面，共有(n+1)个截面，这样就可以得到由(n+1)个方程组成的方程组：

其中M′_i＝F_D(L_r-S_r-l₁-l₂-ia′)(i＝0,1,2,...,n)；

在i＝0处，即B端面处转角和挠度为0，则由转角表达式得f_-1＝f₁又f₀＝0，所以由方程组的第一个方程可以得到：代入第二个方程可以得到f′₂＝a′²(M′₀/K′₀+M′₁/K′₁)，依次代入下去可以得到第i截面处的挠度表达式为

当i取(n+1)时即可求得主簧D处只受到力F_D作用时截面D的挠度

由式(20)计算出f′_n-1，把f′_D和f′_n-1代入式计算出截面D的转角θ′_D：

当主簧受到力F_D作用时，DO段无外部集中载荷作用，且自身重力可忽略不计，因此当主簧只在D处受到力F_D作用时，其末端O点的挠度为：

f′_O＝f′_D+f′_Dθ＝f′_D+θ′_D×(l₁+l₂) (22)

其中，f′_Dθ是D处弯曲转角使O端产生的位移。

5.如权利要求4所述的主副簧式复合材料板弹簧复合刚度计算方法，其特征在于：所述步骤2-3)具体为：把GI段分成(p+1)段，则形成p个新截面，加上固定端面BG，共有(p+1)个截面；设每段长度为a″，则有：

对GI段，根据有限差分理论，可得到由(p+1)个方程组成的方程组：

其中M″_i＝-F′_Dx′＝-F_D(L_r-S_r-l₁-l₂-ia″)(i＝0,1,2,...,p)

同理，可得到第i个截面处的挠度表达式为

当i取(p+1)时即可求得截面I的挠度

6.如权利要求5所述的主副簧式复合材料板弹簧复合刚度计算方法，其特征在于：所述步骤2-4)具体为：在主簧与副簧共同工作后，主簧D点的挠度为副簧末端对D点的集中力F_D与F_O的共同作用产生的，即主簧D点的挠度为f″_D＝f_D+f′_D；而主簧和副簧在接触点D处的挠度相同，因此f″_D＝f_I；根据该式可求出

其中，F_O为已知量；m，n，p为程序中取的各截面数；K_i，K′_i和K″_i分别是各板簧簧身的截面弯曲刚度。

7.如权利要求6所述的主副簧式复合材料板弹簧复合刚度计算方法，其特征在于：所述步骤2-5)具体为：主簧O点的挠度f″_O是集中力F_D与F_O的共同作用产生的，即

f″_O＝f_O+f′_O＝(f_E+θ_E×l₁)+[f′_D+θ′_D×(l₁+l₂)] (28)

其中，f_E及θ_E与F_D无关，f′_D及θ′_D与F_D相关；将求出的F_D代入上式即可求出主簧和副簧共同工作后主簧O点的挠度f″_O，则主副簧总成后半段的复合刚度为：

8.如权利要求7所述的主副簧式复合材料板弹簧复合刚度计算方法，其特征在于：所述步骤2-5)具体为：通过下式计算出主副簧总成的复合刚度：

其中，K_f和K_r分别为总成前半段和后半段的复合刚度，L_f和L_r分别为总成前后半段半长。