非線性-弧長法-源代碼

    在之前的帖子《非線性行為初識》中，我們通過簡單的彈簧桿件結構介紹了非線性問題，牛頓代法和弧長法。回顧一下該問題：如圖所示，中間節點作用一個F的力，會產生一個位移v，

由靜力平衡關系可得到

對于弧長法，我們已經介紹過其原理和迭代格式，下面為其具體的算法流程。

采取弧長法求解上面的問題，取如下約束方程：

該函數為一個圓，這更清晰的說明了弧長法的含義，下圖為k=0時的載荷位移曲線，除了極值點處有一些不足(代碼未加弧長控制)，弧長法得到了完整地載荷位移曲線。

python版源代碼如下，歡迎討論