下面將貼出我用python寫得一維彈簧單元的直接剛度法：

例子計算：

結果如下：

可以看到，輸出結果和書上的答案一致。

下面進行ABAQUS模擬：

添加彈簧單元

添加邊界條件

顯示單元編號、節點編號如下所示，紅色表示單元編號，黃色是節點編號

ABAQUS計算結果如下：

首先是變形圖前后對比

位移云圖如下所示，基本和直接剛度法計算的結果一致

整體剛度矩陣如下所示，因為ABAQUS彈簧單元是三維的，每一個節點有3個自由度，15x15，原味的剛度矩陣如下

我們把剛度矩陣轉化為一維的，方便和編程的結果對比

從結果可以看出，ABAQUS的整體剛度矩陣和直接剛度法計算出來的整體剛度矩陣有些差異，如圖標紅的所示。

那么在整體剛度矩陣上為什么ABAQUS會和直接剛度法的整體剛度矩陣有差異呢，到底ABAQUS的整體剛度矩陣對不對呢？答案將在下一期揭曉。歡迎大家積極討論。