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

本次分享的是基于Modified Newton–Raphson的非線性有限元編程小案例，將會(huì)介紹牛頓法的優(yōu)缺點(diǎn)，以及如何一步步修正，最后結(jié)合彈簧單元的例子進(jìn)行Modified Newton–Raphson數(shù)值實(shí)現(xiàn)，望大家喜歡。

Abaqus/Standard應(yīng)用Newton-Raphson算法解決非線性問題，木木本期就為同學(xué)們“盡可能”全面講解該算法，從Abaqus內(nèi)部算法到數(shù)學(xué)問題中的非線性方程Newton-Raphson算法理論，最后結(jié)合具體非線性方程給出相應(yīng)的代碼，如此一來，更加生動(dòng)地演繹Newton-Raphson迭代過程。

有關(guān)Newton–Raphson的算法細(xì)節(jié)，可參考數(shù)值分析教材或者木木前幾期的推文，本質(zhì)就是一個(gè)切線逼近的過程，本次推文的目的就是帶著大家將該算法應(yīng)用到具體有限元非線性求解案例中。

主程序：Newton迭代法計(jì)算方程的根 use NEWTON OPEN(UNIT=11,FILE='FOUT1510.TXT') OPEN(UNIT=12,FILE='IM_RESULT1510.TXT') CALL SOLVE(X,ITER) WRITE(11,46)X,ITER46 FORMAT(T5,'Newton迭代法計(jì)算方程的根',//,&

? 收斂速度：Newton-IHTs方法收斂最快。? 運(yùn)行時(shí)間：Newton-IHTs方法效率顯著優(yōu)于SD和IHTs方法。

Newton-Raphson update d_ga = d_ga - f_gamma_delta / Df_gamma_d_d_gamma END DO !

非線性有限元一直是有限元中較為困難的一部分，在非線性有限元中我們經(jīng)常碰到諸如Newton-Raphson迭代法，切線剛度陣等概念，今天我們就單的介紹一下非線性吧。

分析和求解控制 使用初始Newton-Raphson方法進(jìn)行非線性靜態(tài)分析。載荷極限（結(jié)構(gòu)能夠承受的最大載荷）通過全局Newton-Raphson解的發(fā)散來確定。 在單個(gè)子步中計(jì)算原位應(yīng)力狀態(tài)。 使用0.05的初始時(shí)間步長計(jì)算載荷極限。為了精確捕捉結(jié)構(gòu)完整性損失，時(shí)間增量可以減少到0.001。

# 求解器constraints Plain;         # 約束處理numberer RCM;             # 自由度編號(hào)test EnergyIncr 1.0e-6 10;  # 收斂測試algorithm Newton

全局剛度矩陣是不對稱的，因此需要非對稱Newton-Raphson求解器選項(xiàng)（NROPT，UNSYM）。

牛頓-迭代硬閾值（Newton-IHTs）算法：是IHTs的性能升級(jí)版本，在IHTs的迭代邏輯基礎(chǔ)上，額外引入Hessian矩陣的逆矩陣參與中間結(jié)果計(jì)算——這一改進(jìn)能顯著加快迭代收斂速度，進(jìn)一步縮短非線性光刻優(yōu)化的運(yùn)算周期。

示例包括 Newton's Method 和 Quasi-Newton Methods。這些方法可以比 Gradient Descent 更快地收斂，但需要更多的計(jì)算并且可能不太穩(wěn)定。 Gradient Descent 是一種迭代優(yōu)化算法，用于查找函數(shù)的最小值。

圖7 仿真模型和仿真結(jié)果對比 表5 計(jì)算速率比對 2、自適應(yīng)時(shí)間步進(jìn)算法 基于多級(jí)時(shí)間步長控制算法，使用Newton迭代次數(shù)作為非線性的指示，減少線性部分的時(shí)間步數(shù)，增加弱/中等非線性段的時(shí)間步數(shù)。從而減少總的時(shí)間步數(shù)，降低Newton迭代次數(shù)、縮短仿真時(shí)間。

F(r)形式較復(fù)雜，很難求得準(zhǔn)確的解析解，因此可用Newton-Rapson迭代法求出數(shù)值解，鑒于Newton-Rapson迭代法對初值的要求較高，可先采用二分法大致確定初值范圍后再進(jìn)行迭代計(jì)算，迭代方程如下：式中：rk+1為第k+1個(gè)參數(shù)r;rk為第k個(gè)參數(shù)r;F(rk)為關(guān)于rk的函數(shù)；F'(rk)為關(guān)于rk的函數(shù)的導(dǎo)函數(shù)。

介紹弧長法之前，必須了解Newton-Raphson法的載荷控制和位移控制： 如下圖的位移-載荷曲線，如果使用載荷控制，只能夠達(dá)到Fcr。如果使用位移控制，有可能會(huì)跳過不穩(wěn)定點(diǎn)，但是必須要知道是什么位移，在復(fù)雜載荷下，一般不知道位移狀態(tài)。

Solve for Smises and deqpl using Newton's method deqpl=0.d0 Et=E*xn*(1.d0+E*eqplas/Sy)**(xn-1) do kewton=1,newton rhs=Smises-(3.d0*eg)*deqpl-Sf deqpl=deqpl+rhs/((3.d0*eg)+

在模型中，點(diǎn)擊在線靜平衡分析對話框進(jìn)行靜平衡計(jì)算，計(jì)算時(shí)間相比Newton法明顯變長，得到的仿真結(jié)果如下，說明模型已經(jīng)基本處于平衡狀態(tài)。得到平衡狀態(tài)下的模型如圖所示。4.

在input文件里，坐標(biāo)表示為，定義節(jié)點(diǎn)組合與邊界條件為，為了讓模型收斂性更好，采用quasi-newton 求解器。

Phase field fracture modelling using quasi-Newton methods and a new adaptive step scheme. Theoretical and applied fracture mechanics. 2020;107:102446.