橡膠模型的穩定性
動畫中顯示的橡膠構件是使用不同的橡膠材料模型進行分析的。當使用穩定性指數為負數的橡膠材料模型時這個分析是不收斂的。橡膠材料的穩定特性可以通過ADINA8.5（8.5.2或更高版本）中新的穩定性視圖功能顯示出來。
關鍵詞：橡膠材料，ADINA，穩定性，Mooney-Rivlin，Sussman-Bathe

橡膠材料模型

圖1和2顯示的是對某種橡膠材料做單軸和雙軸拉伸時的實驗數據。這些圖畫出的是工程應力-應變關係曲線，這是描述橡膠類材料常用的關係曲線。圖中也顯示了兩種擬合得很好的橡膠材料：9項的Mooney-Rivlin材料和Sussman-Bathe材料。這兩種材料的擬合結果與實驗數據非常吻合。

圖3顯示的是根據單軸拉/壓得到的數據，用真實應力-對數應變來表示。注意這條曲線包含了受拉區和受壓區。

圖1 工程應力-應變，單軸拉伸

圖2 工程應力-應變，雙軸拉伸

圖3 真實應力-對數應變。這條曲線是在圖1和2顯示的單軸和雙軸拉伸曲線的基礎上得到的

橡膠構件的分析

上邊的兩種材料模型都用在了橡膠構件分析中。橡膠構件在四個方向上受相同大小的位移荷載，如圖4所示。

圖4 橡膠構件的初始和*終狀態

圖5顯示了不同的材料模型得到的力和變形的曲線。使用9項的Mooney-Rivlin材料模型，當位移大於1.4時無法得到收斂解。但是使用Sussman-Bathe材料模型，在很大位移荷載條件下仍然能得到收斂解。

圖5 力和變形的曲線

穩定性視圖

這個全新的穩定性視圖功能提供了一個查看橡膠構件收斂性的窗口。圖6和7顯示的是兩種材料模型的穩定性視圖。

圖6 由材料數據得到的穩定性曲線-Sussman-Bathe材料模型（Material1表示隻定義了一種材料）

圖7 由材料數據得到的穩定性曲線-Mooney-Rivlin材料模型

下麵是穩定性視圖的一些想法。考慮一個均勻的橡膠片受單軸拉伸作用。對每一個應變等級，增加的剛度陣（相應於變化的位移荷載產生的變化的力）是確定的。計算增加的剛度陣的特征值，並且將*小的特征值作為穩定性指數。如果穩定性指數大於零，則材料穩定（對於施加變化的力），否則材料不穩定。對於純剪力和雙軸拉伸執行相同的過程。

穩定性視圖顯示Sussman-Bathe材料模型對三種模式的變形都是穩定的，但是9項的Mooney-Rivlin材料模型在雙軸拉伸作用下真實應變大於0.4時就不穩定了。由於橡膠材料是雙軸拉伸，所以用9項Mooney-Rivlin材料在很小的荷載/變形下不收斂也就不奇怪了。

討論

我們希望材料模型有一個非常好的特點，如果實驗得到的應力應變數據符合一種非常穩定的材料，則這種材料模型也應該是穩定的。顯然，在這個例子裡，9項的Mooney-Rivlin材料模型冇有這個特點。

當然，在Mooney-Rivlin材料模型中為了使穩定性指數為正數，可以定義不同的材料常數，但是這樣材料模型就不能很好得擬合實驗數據了。