不知火舞的被虐|伊人天伊人天天综合网|博洛尼亚天气|任你懆这里只有精品4|久久美日韩精品久久|掌中之物漫画免费阅读观看|0丨d老妇

今日推文主要分享一個非線性方程的牛頓-拉弗森迭代解法，借助Fortran語言，講述Fortran編程時需要注意的地方。理論及在Abaqus中的實現過程已在上幾期推文基于Abaqus的Newton-Raphson算法中說明，本次主要說明Fortran編程時需要注意的地方，本文代碼主要參考：《Fortran程序設計權威指南》。

非線性靜力分析需要求解非線性平衡方程，程序采用全牛頓-拉弗森或修正全牛頓-拉弗森迭代法。許多問題涉及歷程相關的響應，因此通常通過一系列的增量步得到求解結果，在每個增量中通過迭代得到平衡。有時必須保持增量較小(意味著轉動和應變增量必須小),以確保對歷程相關影響的正確模擬，但最常見的是，增量步大小的選擇是計算效率的問題，如果增量步太大則需要更多迭代。

別忘了思考牛頓-拉夫森方法在這種情況下扮演什么角色

這是一個完美的模擬示例，用牛頓-拉夫森方法來展示幾何形狀、材料和接觸非線性。 橡膠靴形密封件在許多工業應用中用于保護柔性接頭 在兩個物體之間。在汽車行業中，橡膠防塵罩密封件持續覆蓋著 驅動軸上的速度接頭，用于保護其免受外部因素（如灰塵）的影響，潮濕、泥濘等環境。 這些橡膠靴的設計旨在適應這些環境關節的最大可能擺動角度，以及補償軸長變化。

需要使用牛頓-拉普森方法來計算位移。在撓曲較大的情況下，需要仔細設置材料屬性、尺寸和時間增量。如果膜沒有任何殘余應力或預應力，則垂直方向的剛度值非常小。由于模型應該使用牛頓-拉普森迭代，這會導致在第一步時計算非常不穩定。膜在所有邊緣處都受到所有自由度的約束。毫無疑問，膜應包含足夠細的網格以計算最大撓曲。因此，這些都是計算膜撓曲時需要克服的挑戰。

以下是一些額外的建議：1）繪制剩余力：牛頓-拉夫森剩余力的高值通常表明導致不收斂的特定接觸對。2）在接觸區域細化網格：這將使接觸壓力分布在更多的單元上，并增加接觸點的數量。相對較少的接觸點可能會導致非常高的接觸應力，從而導致單元過度變形和收斂困難。對于非線性材料，這尤其成問題。

而非線性分析中，剛度矩陣隨計算過程變化，需要通過牛頓-拉夫遜法等算法進行多次迭代，計算量呈幾何倍數增長。 2?? 幾何非線性 (Geometric Nonlinearity) 當結構發生“大位移”、“大轉動”或“大應變”時，初始構型發生顯著改變（如釣魚竿受力）。此時，必須開啟大變形開關，以修正剛度矩陣對構型變化的響應。

在ANSYS里還是牛頓－拉普森法和弧長法。牛頓－拉普森法是常用的方法，收斂速度較快，但也和結構特點和步長有關。弧長法常被某些人推崇備至，它能算出力加載和位移加載下的響應峰值和下降響應曲線。但也發現：在峰值點，弧長法仍可能失效，甚至在非線性計算的線性階段，它也可能會無法收斂。 為此，盡量不要從開始即激活弧長法，還是讓程序自己激活為好（否則出現莫名其妙的問題）。

上述為非線性等式，我們須將其正切剛性線性化，并迭代求解。線性化后的平衡系統可以以下公式表示：為了進行迭代計算，我們采用牛頓-拉弗森方法──解決非線性數學問題的最著名的方法。此分析會持續收斂，直到殘余力小于收斂標準時即完成解答。幾何擾動(Imperfection)模型應用有時在數值分析過程中，結構分析中不容易出現非線性情況。

上述為非線性等式，我們須將其正切剛性線性化，并迭代求解。線性化后的平衡系統可以以下公式表示： 為了進行迭代計算，我們采用牛頓-拉弗森方法──解決非線性數學問題的最著名的方法。此分析會持續收斂，直到殘余力小于收斂標準時即完成解答。 幾何擾動(Imperfection)模型應用 有時在數值分析過程中，結構分析中不容易出現非線性情況。

程序中塑性迭代部分并非主流的牛頓-拉夫遜和梯度下降方法，但經過驗證能夠適用于該模型，介意請勿拍。附件內容：1. inp算例模型(低速沖擊工況，1/4模型，層間使用cohesive element) 2. 子程序 3 .使用方法 4.參考論文名稱首先介紹該子程序的使用方法與效果1. 在ABAQUS中建立三維復合材料模型，這里建立一個簡單的方塊。