Hi,“Radioss乾貨”技術專題又和大家見面啦~ 關於彈簧的相關問題,之前我們給大家分享過一期“彈簧單元失效”,今天我們就來講講彈簧的內力和彎矩計算。
彈簧在Altair Radioss™ 中有以下三種不同的模型,如下圖所示依次為:純彈簧模型、純阻尼器模型以及彈簧和阻尼器串聯模型。

不同彈簧模型有不同的彈簧內力和彎矩的計算,以及不同的時間步長的計算,並且這三種彈簧模型都可以描述彈簧的線性和非線性的力學性能,比如,彈簧單元屬性TYPE4,8 ,12,13,25 都是這樣的類型。

首先,線性部分力和彎矩如下計算:

以TYPE13 彈簧單元中的i=1 軸向拉壓為例,線性階段需要輸入彈簧剛度係數和阻尼係數即可。
然而,非線性部分的力和彎矩計算較為複雜,公式如下:
▇ i =1, 2, 3
▇ i =4, 5, 6
力矩的計算和內力的計算非常相似,下面我們以內力的計算為例,將看上去有些複雜的計算分解一下,它大致可以分為:
□ 靜態載荷下的力和位移f 曲線;
□ 動態載荷下的考慮伸長率的力和位移關係
- 使用對數函數形式(形式一)描述
- 使用與f 曲線相關的g 係數曲線輸入(形式二)
- 使用與f 曲線無關的獨立的h 曲線輸入(形式三)
- 考慮動態載荷下的線性阻尼作用

Radioss
靜態載荷下的力和位移關係
由輸入曲線fct_ID 1i 來描述靜態載荷下的力和位移關係(以下簡稱f曲線)。
注意這個輸入的力和位移的f曲線通常在曲線最後有一個非常高的斜率以便在彈簧受壓時(尤其受到比較大的壓力時)不會壓潰。如果輸入曲線的斜率不高那麼Radioss 會給出下面的警告信息:
WARNING ID: 506
** WARNING IN SPRING PROPERTY
** WARNING IN SPRING PROPERTY SET ID=XXX
STIFFNESS VALUE 100 IS NOT CONSISTENT WITH THE MAXIMUM SLOPE (4550)
OF THE YIELD FUNCTION ID=X
THE STIFFNESS VALUE IS CHANGED TO 1000
所以非線性靜態載荷階段需要用戶輸入f曲線(fct_ID 1i)及用於f曲線x軸係數A i和y軸係數Ascale1即可。以TYPE13彈簧單元中的i=1軸向拉壓為例,如下圖所示紅色標註的參數需要輸入,這裡用於描述彈簧線性的剛度K i或者阻尼C i(根據選擇的彈簧模型不同)是一定需要輸入的。
Radioss
動態載荷下的力和位移關係
動態載荷下的力和位移關係,實際上是需要考慮不同的伸長率下力和位移的關係,在TYPE4,8,12,13,25中基本可以分為以下三部分,分別描述或組合描述:
▇ 使用對數函數形式(形式一)

由於使用對數形式的函數這種方式考慮動態效應時仍然需要f曲線。所以對數函數實際上是描述力的係數和伸長率的關係(如下圖所示)。

使用參數輸入
彈簧在伸長率(elongation velocity)超過D i後,彈簧因伸長率而內力額外增加。
是力的係數,最終的力是這個係數乘以準靜態下的f曲線來描述,並且參數B i是這個對數曲線的比例因子,用於調節曲線中的力的係數的大小。
所以使用這種形式考慮彈簧動態效應下的力,需要用戶除了輸入靜態載荷時需要的f曲線(fct_ID 1i)及用於f曲線x軸係數A i 和y軸係數Ascale1以外,還需要輸入係數B i,D i 用於描述對數函數。以TYPE13彈簧單元中的i=1軸向拉壓為例,如下圖所示紅色標註的參數需要輸入。
▇ 使用與f曲線相關的g係數曲線輸入(形式二)

這裡輸入的曲線(輸入在fct_ID 2i中的係數曲線,以下簡稱g曲線)是彈簧各個伸長率下的力的比例係數曲線,也就是說最終的對於某一時刻彈簧某一伸長率下的力和位移的實際曲線是f曲線和g曲線的乘積f · g。同樣參數E i 是這個比例係數曲線的比例因子,用於調節g曲線。

使用g 曲線
所以使用這種形式考慮彈簧動態效應下的力,同樣需要用戶除了輸入靜態載荷時需要的f曲線(fct_ID 1i)及用於f曲線x軸係數A i和y軸係數Ascale1以外,還需要在輸入g曲線(fct_ID 2i)和g曲線的x軸係數F i 和y軸的係數E i。以TYPE13彈簧單元中的i=1軸向拉壓為例,如下圖所示紅色標註的參數需要輸入。
▇ 使用與f曲線無關的獨立的h曲線輸入(形式三)
這裡輸入的曲線(輸入在fct_ID 4i中的係數曲線,以下簡稱h曲線)是一個獨立的伸長率和力的關係的曲線。它和上面的g曲線的區別在於:
- 對於某一伸長率下,通過g 曲線(形式二)得到相應伸長率下得力的比例係數,然後這個比例係數乘以f 曲線才是最終的力和位移的曲線。
- 同樣對於某一伸長率下,通過h 曲線(形式三)相應伸長率下的力的值,然後這個力的值直接加到f 曲線上去即是最終的力和位移的曲線。

使用h 曲線
參數Hscale i是這個h曲線的比例因子,用於調節曲線中的力的大小。
所以使用這種形式考慮彈簧動態效應下的力,需要在輸入h曲線(fct_ID 4i)和h曲線的x軸係數Fi和y軸的係數Hscale i 。以TYPE13彈簧單元中的i=1軸向拉壓為例,如下圖所示紅色標註的參數需要輸入。
使用這種形式既可以不考慮靜態載荷,只輸入h曲線,也可以同時考慮靜態載荷,再輸入需要的f曲線(fct_ID 1i)。

阻尼部分使用線性阻尼來描述動態載荷下的內力計算。
注意
關於Altair
Altair(納斯達克股票代碼:ALTR)是一家全球技術公司,在產品開發、高性能計算(HPC)和數據分析等領域提供軟件和雲解決方案。
自1985年成立以來一直致力於為企業開髮用於仿真分析、優化、數據分析、信息可視化、流程自動化、雲計算以及物聯網領域的高端技術。
公司總部位於美國密歇根州,服務於11000多家全球企業,應用行業包括汽車、消費電子、航空航天、能源、機車車輛、造船、國防軍工、金融、零售等。
欲了解更多信息,歡迎訪問: