葉片曲率對少片變截面鋼板彈簧力學特性的影響*

楊林 于曰偉 周長城 鄭偉

（山東理工大學，淄博 255049）

1 前言

鋼板彈簧作為彈性元件兼導向機構，目前廣泛應用于乘用車和商用車。隨著節能環保、汽車輕量化等相關政策的實施，少片簧以其質量輕、噪聲小、性能好、壽命長等優點，日益受到重視。為滿足不同工況下的應力約束條件、剛度設計要求和加工工藝需要，少片簧通常加工成拋物線型、斜線型、根部加強型、端部加強型、兩端加強型等不同結構型式。針對不同結構型式的少片簧，研究人員建立了相應的少片變截面鋼板彈簧剛度、應力計算公式[1-5]。為了使分析結果更加精確，研究人員開始考慮材料[6-10]、片間摩擦[11-14]、安裝約束[15-16]等因素的影響，并利用有限元法對少片變截面鋼板彈簧的力學特性進行了分析[17-18]，為少片鋼板彈簧的前期設計和強度評價提供了科學依據。然而，這些研究均未考慮葉片曲率對少片簧力學特性的影響，不可避免地使設計結果與實際情況存在較大的偏差。

在汽車行駛過程中，少片鋼板彈簧除了承受車輛自身及載物質量外，還要承受因路面不平引起的沖擊載荷，其主要失效方式是疲勞斷裂和產生永久性塑性變形[19]。為在一定程度上減小板簧應力、延長板簧壽命，使各片板簧壽命接近，通常在自由狀態下將各片板簧加工成不同的曲率，裝配壓緊前各片間保持一定自由間隙，裝配后各片緊密貼合，呈現自然向上彎曲的狀態。但在少片鋼板彈簧的設計、分析計算中，往往忽略葉片曲率的影響，將少片簧視為平直的簧片疊加，導致了實際應用時少片簧的性能和壽命難以與設計值相吻合。

為了在現有基礎上進一步提高少片變截面鋼板彈簧的分析、設計精度，本文通過考慮葉片曲率對少片簧的影響，建立符合少片簧實際結構特征的有限元模型，對不同曲率的少片變截面鋼板彈簧的力學特性進行分析，以期為少片鋼板彈簧的設計和疲勞壽命預測提供有效的理論參考和技術借鑒。

2 少片變截面鋼板彈簧建模與仿真

2.1 結構特征

通常，少片變截面鋼板彈簧由1～4 片中間厚兩端薄、等寬、等長、曲率半徑不同的彈簧鋼片疊合而成，其中，少片簧的一半對稱結構根據厚度特征劃分為3個部分，即根部平直段、變厚度段和端部平直段。為保持等強度特性，變厚度段厚度可按拋物線型或線性變化。板簧中部安裝中心螺栓以保證裝配時各片的相對位置，形成一根彈性梁縱置于車橋與車架之間，如圖1 所示。

圖1 少片變截面鋼板彈簧

2.2 仿真模型的建立

根據少片變截面鋼板彈簧的實際結構尺寸，考慮葉片曲率的影響，利用SolidWorks軟件建立相應的少片簧仿真模型。模型簡化原則為：

a.考慮到少片簧沿片長方向的對稱性，在確保不影響結果精度的前提下，建立鋼板彈簧的1/2模型進行仿真。

b.假設板簧在折彎成一定曲率后，各位置的厚度不變。

c.假設第n片板簧與第(n+1)片板簧之間的墊片厚度與第n片板簧端部平直段厚度之和等于第n片板簧根部平直段厚度。

2.2.1 上表面曲線的確定

基于少片簧結構以及設計、安裝要求的特殊性，折彎加工前其葉片上表面為水平直線段，折彎成一定曲率后，除根部平直段被騎馬螺栓固定仍為直線段外，其余部分變為該片板簧對應曲率圓的一段圓弧，幾何關系如圖2所示。

圖2 少片簧幾何關系

圖2中，O0、R分別為板簧上表面曲率圓的圓心和半徑，h1、h2分別為板簧端部和根部平直段厚度，直線段OA、直線段O′A′分別為1/2 板簧模型的根部平直段上、下表面曲線，、曲線A′B′分別為模型的變厚度段上、下表面曲線，分別為模型的端部平直段上、下表面曲線，H為弧高，即板簧根部平直段上表面與卷耳下部切線間的垂直距離，θ為對應的圓心角。

由圖2可得：

式中，LAC為的長度。

由于各片板簧被壓緊且由螺栓固定，根據模型簡化原則c，其余各片上表面曲率半徑可根據第1片確定。以3片簧為例，設第1片、第2片、第3片葉片上表面曲率半徑分別為R1、R2、R3，根部平直段厚度分別為h12、h22、h32，則

2.2.2 下表面曲線的確定

折彎加工后，對應的少片簧的下表面曲線由直線段O′A′、曲線A′B′、組成。其中是以O0為圓心，(R+h1)為半徑的一段圓弧。以A為原點，根據簡化原則b，曲線A′B′上某點(x1,y1)與曲率圓圓心連線交于點(x0,y0)，點(x0,y0)與點(x1,y1)的連線長度為hx，即折彎加工前點(x1,y1)對應的葉片厚度。利用角度θ建立(x0,y0)與(x1,y1)的關系：

式中，θx為點(x0,y0)與圓心O0連線與垂線O0A的夾角。

運用MATLAB 進行求解，得出曲線A′B′上一系列點，通過對這些點進行8 次曲線擬合，可得到曲線A′B′的表達式。以某片板簧為例，該板簧結構參數如表1所示。

表1 結構參數 mm

所得到的曲線擬合結果如圖3 所示。其中曲線A′B′的表達式為：

圖3 MATLAB曲線擬合

2.2.3 三維模型的建立

將得到的曲線擬合表達式輸入SolidWorks軟件中，得到少片簧三維模型。

2.3 有限元仿真過程設置

2.3.1 定義單元屬性和材料屬性

定義單元類型為SOLID186，SOLID186 是一個高階3維20節點固體結構單元，每個節點有3個分別沿x、y、z方向平移的自由度，支持塑性、超彈性、蠕變、應力鋼化、大變形和大應變，三維實體接觸單元采用TARGE170及CONTAC174來模擬片與片之間的接觸和力的傳遞。鋼板彈簧的材料為60Si2Mn，彈性模量為206 GPa，泊松比為0.29，密度為7.74×103kg/m3，屈服強度為1.176 GPa，抗拉強度為1.274 GPa。

2.3.2 少片簧的網格劃分

為提高網格質量，對模型進行合理分塊，網格劃分主要采用六面體單元，根據板簧結構參數設置合適的單元尺寸來保證計算精度，一般設置為5 mm。

2.3.3 定義邊界條件

在兩葉片根部平直段、端部平直段分別定義1個接觸對來模擬葉片間接觸和力的傳遞，每個接觸對都由1個接觸面和1 個目標面組成，定義凸曲面為接觸面，凹曲面為目標面，接觸類型定義為摩擦接觸，按簧片在良好的潤滑狀況下計算，摩擦因數設置為0.2，接觸剛度系數為0.001。由于各葉片之間采用中心螺栓連接且在鋼板彈簧中部用U形騎馬螺栓夾緊在橋殼或橋體上，故在各葉片根部施加寬度方向上的固定約束，各葉片的端部為自由端，允許其自由變形。

2.3.4 加載-卸載和求解

集中載荷加載在首片少片簧端部的節線上，為了獲得不同載荷下少片簧的變形，以分步加載的方式施加載荷，加載過程分解為7 個載荷步，逐級加載至最大載荷5 280 N。卸載過程同樣分為7 個載荷步，逐級卸載至空載。由于少片簧在受力過程中產生大變形，故打開求解控制中的幾何大變形選項。因板簧間的接觸的特點為不穿透、自由的分離和帶摩擦的滑移，故關閉修剪接觸（Trim Contact）并選擇接觸算法為罰函數法（Pure Penalty），以提高計算準確性。添加總位移分布云圖（Total Deformation）和應力分布云圖（Equivalent Stress）分析求解。

3 模型試驗驗證

以某少片斜線型變截面鋼板彈簧為例，對所提出的建模方法的正確性進行試驗驗證。該少片簧的結構參數如表2所示。

表2 結構參數

3.1 ANSYS仿真

根據所建立的仿真模型搭建方法，構建的少片斜線型變截面鋼板彈簧仿真模型如圖4所示。

圖4 仿真模型

利用ANSYS有限元分析軟件對該少片簧進行靜力學特性仿真，通過解算器輸出檢查，確定該有限元模型共計52 422個節點、12 103個單元。

3.2 試驗結果與仿真結果對比

根據GB/T 19844—2018所規定的試驗方法，對表2所示少片簧進行力學特性試驗，試驗測試設備為TPW-100 型鋼板彈簧靜剛度試驗機，試驗平臺如圖5 所示。試驗中，逐級對該少片簧進行加載和卸載，測試少片簧的變形和應力等參數。通過有限元軟件對曲率k為0和0.8 的板簧模型進行加載-卸載仿真，所得到的試驗數據與仿真數據對比結果如圖6所示，力學特性對比結果如表3所示。

圖5 試驗平臺和測試過程

圖6 試驗與仿真對比結果

表3 力學特性對比結果

由圖6和表3可知，加載至最大載荷時，仿真板簧1與仿真板簧2 的變形與試驗結果的相對偏差分別為7.8%和1.0%，對加載過程的變形曲線進行線性擬合，得到剛度數據，仿真板簧1與仿真板簧2的剛度與試驗值的相對偏差分別為7.987%和0.403%。結果表明，仿真板簧2與試驗板簧剛度吻合度較高，即建模時考慮曲率能夠較好地反映少片簧的剛度特性。產生相對偏差的主要原因在于：建模時忽略了實際板簧在由平直狀態折彎成一定曲率的過程中因受拉力、壓力而造成的葉片厚度變化。加載至最大載荷時，仿真板簧1 與仿真板簧2的最大應力均出現在板簧根部，與試驗結果的相對偏差分別為9.94%和7.83%，由于仿真時沒有考慮裝配所產生的應力，導致應力仿真結果小于試驗結果，為準確模擬該曲率板簧在裝配和受壓后的應力情況，對裝配過程進行仿真，在考慮裝配、加載、卸載3 個階段后，得出的應力對比結果如圖7所示。

圖7 板簧根部應力對比

加載至最大載荷時，仿真板簧2的應力與試驗結果的相對偏差為3.38%，通過上述分析，證明了該模型是準確可靠的。

4 葉片曲率對少片變截面鋼板彈簧力學特性的影響

以3種不同結構參數的少片變截面鋼板彈簧為例，其中少片簧1為斜線型鋼板彈簧，少片簧2與少片簧3為拋物線型主副鋼板彈簧，利用所建立的建模方法搭建不同曲率的少片簧仿真模型，為方便探究曲率的影響，設葉片向上彎曲時曲率k>0，葉片向下彎曲時曲率k<0，葉片平直，即不考慮曲率時k=0，曲率的正負僅表示葉片彎曲方向，曲率大小以|k|表示。對于曲率k≠0的主副簧復合式鋼板彈簧，以其副簧端點與首片主簧上表面對應點間距離為主副簧間隙。進行葉片曲率對少片簧力學特性的影響規律的分析，少片簧結構參數如表4所示。

表4 少片簧結構參數

4.1 葉片曲率對剛度特性的影響

圖8所示為仿真得到的葉片曲率對表4所示少片變截面鋼板彈簧剛度特性影響的曲線。

圖8 剛度特性結果對比

由圖8a可以看出：k<0時，相同端部作用力下，隨著|k|的增大，板簧的變形減??；當曲率k≥0時，相同端部作用力下，隨著|k|的增加，板簧的變形增加，但增加程度較曲率k<0時小。由圖8b、圖8c可以看出：k<0時，相同端部作用力下，隨著曲率的增大，板簧的變形減??；當k>0時，隨著曲率的增加，板簧的變形增加；板簧2 與板簧3是主副簧復合式鋼板彈簧，k=0時，主副簧間隙一定，無法等同于其他曲率板簧的主副簧間隙，逐級加載時，其副簧早于其他非零曲率副簧進入接觸，所以一定載荷范圍內，其剛度較大，變形較小。

為了進一步分析葉片曲率對剛度特性的影響，圖9給出了仿真所得的葉片曲率對表4 所示少片變截面鋼板彈簧剛度影響的曲線，其中板簧1的剛度為加載至最大力下的剛度，板簧2與板簧3的剛度為主副簧完全接觸后的復合剛度。

圖9 不同曲率下的少片簧剛度

由圖9 可以看出：葉片向下彎曲時，隨著曲率的增大，少片簧的剛度增大，且增大趨勢明顯；葉片向上彎曲時，隨著曲率的增大，少片簧剛度逐漸減小，減小趨勢較為平緩。整體看來，圖中曲線近似二次曲線。出現上述情況的原因可能是：葉片曲率變化導致端部作用力的作用線與板簧中心垂線的水平距離發生變化，從而導致了彎矩、片間接觸摩擦的變化，且葉片向上彎曲時，施加載荷有使葉片曲率減小的趨勢，葉片向下彎曲時，施加載荷有使葉片曲率增大的趨勢。這種減小或增大趨勢減緩或加劇了上述變化的程度，這種變化反映在剛度上使剛度曲線在不同曲率下近似呈二次曲線變化。通常，在板簧初始狀態下，葉片呈現向上彎曲的形狀。仿真結果顯示，隨著葉片曲率的增大，板簧抵抗彈性變形的能力降低，在懸架垂直載荷一定時，剛度減小導致懸架垂直變形增大，增加了撞擊限位塊的概率。對于主副簧復合式鋼板彈簧，剛度降低會造成副簧過早進入接觸狀態，導致板簧不能滿足實際工作要求。

為了進一步分析葉片曲率對剛度的影響程度，表5給出了仿真所得的不同非零曲率下的少片簧剛度對于不考慮曲率時的變化量。從表5中可以看出，隨著曲率的增大，變化量逐漸增大，即曲率對剛度的影響程度逐漸增大。通常，少片簧初始狀態下的曲率在0.5～0.8 m-1范圍內，為了便于工程應用，同時簡化求解過程，可以在不考慮板簧曲率進行計算分析的基礎上，以曲率-剛度修正系數γ進行修正，γ的取值范圍為0.06～0.08。

表5 剛度變化量

4.2 葉片曲率對應力特性的影響

圖10 給出了仿真所得的葉片曲率對表4 所示少片變截面鋼板彈簧的第2 片主簧同一節線距離板簧端部不同位置處的應力曲線。

圖10 板簧不同位置處的應力曲線

由圖10a 可以看出：當k=0 時，在端部力的作用下，板簧1 某一節線上的應力由端部平直段到變厚度段迅速增加，而后增加趨勢趨于平緩，并在根部平直段出現最大應力；當板簧1 的葉片曲率改變時，同一節線上的應力整體變化趨勢與曲率k=0時一致，曲率的存在并沒有改變板簧1的應力分布，僅改變了應力的大小。由圖10b 可以看出：當k=0 時，在端部力的作用下，板簧2 某一節線上的應力先是沿長度方向迅速增加，而后在拋物線段出現等應力現象，最后在根部平直段出現最大應力；當板簧2 的葉片曲率改變時，同一節線上的應力分布發生了變化，如曲率k>0 時，在主簧對應副簧端點位置處出現了較大的應力。由圖10c可以看出：當k=0時，在端部力的作用下，應力由端部平直段到變厚度段迅速增加，在主簧對應副簧端點位置處出現了較大的應力，下降后趨于不變，而后又再增加并在板簧根部平直段出現了應力最大值；當板簧3的k<0時，隨著曲率的增大，主簧對應副簧端點位置處的應力逐漸減小，在此位置不再出現較大的應力分布。整體來看，當板簧長度、厚度不變時，曲率的存在可能會改變其應力的大小或分布情況。

為了進一步分析葉片曲率對應力特性的影響，圖11 給出了仿真所得的葉片曲率對表4 所示少片變截面鋼板彈簧的根部最大應力影響的曲線。

圖11 不同曲率下的少片簧根部最大應力

由圖11可以看出，在葉片曲率的影響下，少片簧的應力呈非線性變化。k>0 時，應力隨曲率增加而增加，k<0時，應力隨曲率的增加而減小。出現上述情況的原因可能是：葉片曲率變化導致端部作用力的作用線與板簧中心垂線的水平距離發生變化，從而導致了彎矩、片間接觸摩擦的變化，這種變化反映在應力上使應力曲線在不同曲率下呈非線性。圖11 中板簧2 與板簧3 的應力曲線在k=0 處出現尖點的原因是：板簧2 與板簧3 是主副簧復合式鋼板彈簧，k=0時，主副簧間隙一定，無法等同于其他曲率的板簧。

為了進一步分析曲率對應力的影響程度，表6給出了仿真所得的不同非零曲率下的少片簧應力對于不考慮曲率時的變化量。從表6中可以看出，隨著曲率的增大，變化量逐漸增大，即曲率對應力的影響程度逐漸增大。為了便于工程應用同時簡化求解過程，可以在不考慮板簧曲率進行計算分析的基礎上，以曲率-應力修正系數λ進行修正，λ的取值范圍為0.02～0.04。

表6 應力變化量

5 結束語

本文通過考慮曲率對少片變截面鋼板彈簧的影響，建立了其上、下表面曲線的構造方法，并利用臺架試驗驗證了該方法的正確性；通過有限元分析，對不同曲率的少片簧的力學特性進行了探究，得到以下結論：

a.少片變截面鋼板彈簧剛度及應力隨葉片曲率變化的一般規律為：葉片向上彎曲時，隨著曲率的增大，少片簧的剛度逐漸減小，應力逐漸增大；葉片向下彎曲時，隨著曲率的增大，少片簧的剛度逐漸增大，應力逐漸減小。

b.為了便于工程應用，同時簡化求解過程，可以以曲率-剛度修正系數γ與曲率-應力修正系數λ對不考慮葉片曲率的板簧模型的剛度和應力進行修正，從而提高設計精度，優化板簧設計，延長板簧壽命，保證少片簧更好地滿足實際工作要求。