riks的案例
加筋板屈曲riks法分析實例
算例來源如下圖所示：（對于這種沿著壓載方向剛度均勻的結構來說，不能直接進行riks分析，必須先進行擾動分析）
先進行線性屈曲特征值分析，inp文件請查看附件，僅計算第一階模態，模態的構型如下圖所示：
求得特征值為5.9655，施加的線壓力載荷為100，故屈曲線壓力載荷為596.55，總的屈曲力為32214，第一圖文獻1的給出的屈曲載荷為29160，可見誤差不大，可以接受。
在線性屈曲分析的基礎上，進行riks后屈曲分析，inp文件請查看附件，分別引入初始缺陷為shell厚度的1%、10%、50%、100%進行比較，所得到的整理后的載荷比例因子—位移如下圖所示：
可知：1）初始缺陷的引入對后屈曲行為有較大的影響；2）此壁板結構在初次屈曲失穩后仍可繼續加載（如1%情況），且未出現負剛度。
buckling.rar
BUCKLING-CAE.rar
展開 加筋板屈曲riks法分析實例
算例來源如下圖所示：（對于這種沿著壓載方向剛度均勻的結構來說，不能直接進行riks分析，必須先進行擾動分析）
先進行線性屈曲特征值分析，inp文件請查看附件，僅計算第一階模態，模態的構型如下圖所示：
求得特征值為5.9655，施加的線壓力載荷為100，故屈曲線壓力載荷為596.55，總的屈曲力為32214，第一圖文獻1的給出的屈曲載荷為29160，可見誤差不大，可以接受。
在線性屈曲分析的基礎上，進行riks后屈曲分析，inp文件請查看附件，分別引入初始缺陷為shell厚度的1%、10%、50%、100%進行比較，所得到的整理后的載荷比例因子—位移如下圖所示：
可知：1）初始缺陷的引入對后屈曲行為有較大的影響；2）此壁板結構在初次屈曲失穩后仍可繼續加載（如1%情況），且未出現負剛度。
額外之言：如果僅進行線性特征值分析，所有操作均可在CAE中完成；若想在其基礎上繼續進行后屈曲分析，則需要手動更改inp文件添加*NODE FILE語句，以生成后屈曲分析所需要的fil文件；后屈曲分析也需要在inp文件中手動添加*imperfection語句，詳細情況請參看分析手冊。所以本帖附件沒有給出CAE文件，僅給出inp文件。希望此小例子能給大家帶來幫助！
buckling.rar
BUCKLING-CAE.rar
展開 abaqus用riks仿真出的結果和implicit不中準靜態模擬結果一樣，哪個可信度更高？
abaqus用riks仿真出的結果和implicit不中準靜態模擬結果一樣，哪個可信度更高？
圓拱曲線彎曲分析
在Riks程序中，不能指定負載大小的實際值。相反，它們作為解決方案的一部分計算，因為“負載比例因子”乘以加載數據線上給出的負載大小。用戶規定的負載大小僅用于定義方向并估計步長的負載初始增量的大小。該初始負載增量是初始時間增量與時間段的比率以及加載選項中給出的負載量值的乘積。用戶可以通過指定最大負載比例因子或節點處的最大位移或兩者來終止Riks分析。當計算出這些限制中的任何一個的解點時，分析將停止。無論如何，或者如果兩個選項都不被使用，分析結束時超出步長的最大增量。
在諸如這些的直通式研究中，結構可以在完全捕捉之后承載增加的負載。因此，通過指定最大負載比例因子便于方便地進行分析。
對于夾緊的淺拱，初始咬合發生在約-1000（力/長度2單位）的壓力下。因此，-250（力/長度2單位）似乎是要應用的負載的第一增量的合理估計。因此，在1.0的時間段和-5000的壓力負荷（力/長度2單位）中規定了初始時間增量為0.05。該解決方案將在約-2000（力/長度2單位）的壓力下充分發展。因此，當負載比例系數超過0.4時，分析結束。
為了說明在幾個步驟中使用Riks，還包括第二步，其中將壓力從拱頂上取下，使其將朝向其初始配置。在Riks分析中的任何一點，實際負載由下式給出： p0在前一步驟結束時的負載 是當前步驟中規定的負載大小，并且 是負載比例系數。拱被卸載，使得在初始時間增量中，去除約0.15 的壓力。在1.0的時間段內使用0.05的初始時間增量，對該重新啟動的步驟規定了負載。此外，我們希望分析結束時，所有的負載被刪除，并且拱已經返回到其初始配置。因此，為拱的中心設定位移閾值0.0。當這個限制被越過時，分析結束。因為Abaqus必須在初始Riks步驟結束時接收負載大小以啟動下一步，所以Riks步驟之后的任何步驟只能作為上一步驟中的重新啟動作業完成。
展開 
第二章 abaqus分析收斂準則
（例如，* STEP, UNSYMM=YES）
10 對于全局不穩定問題，如全局屈曲、坍塌或突變，其中非線性不穩定區域是發生突變4的區域，平衡路徑從一個穩定點A轉到另一個新的穩定點B，可以使用Riks5。如果使用RIKS，在創建使用RIKS的附加步驟之前，在需要之前不要使用RIKS。必須注意，使用位移控制比RIK更有效。對于RIKS分析中的回溯，請在*STATIC，RIKS下指定最大弧長，如1.5。11 對于局部不穩定問題，使用自動穩定和監測阻尼能量。這不能與RIK一起使用，但可以與位移控制*STATIC, STABILIZE *ENERGY OUTPUT,*ENERGY PRINT或 *ENERGY FILE，以監視能量ELSD、6ESDDEN7和ALLSD。
12 在*PLASTIC材質定義的完全塑性區域中添加略微增大的坡度。
13 將混合單元用于高度不可壓縮單元（泊松比接近0.5）或各向異性超彈性公式（單元中的大剛度差異，如彎曲與軸向剛度）。
14松散的收斂條件（如果可能，請避免）。在為后續步驟使用默認參數之前，可能需要在需要接觸的初始小步驟中進行此操作。*controls, PARAMETERS=FIELD.。
2.5 接觸穩定工具
這里，提供了一些適用于靜態平衡的建議，這取決于接觸相互作用的故障排除。
1 創建一個非常小的初始步驟來啟動接觸。
2 用位移控制代替負荷控制。將所需的節點力和位移寫入.dat文件，然后利用x y data特性生成一個（x-y）荷載與位移數據文件，繪制在viewport中。
3 添加與總負載相比剛度較低的彈簧，以使接觸對產生一些阻力，直到建立接觸。如果彈簧力太大，則可以建立第二個步驟，以便在建立接觸后拆下彈簧。四。為了在沒有剛體運動的情況下獲得初始階段建立的接觸，應使用逼近參數。
展開 ABAQUS生死單元（Model change）入門
You can create a model change interaction in all
Abaqus/Standard steps except for static, Riks steps and linear perturbationsteps. For a detailed discussion, see Element and contact pair removal andreactivation.”
通過這段說明我們可以知道，生死單元的使用條件是：in all Abaqus/Standard steps except for static, Riks steps andlinear perturbation steps。
生死單元在CAE界面實現方法為： create interaction》選擇需要移除或者激活的分析步 step并選擇model change，如下圖中選擇step-3使用 model change。
點擊繼續之后，可選取相應的區域來施加model change以及相關選項。
下面這個案例簡要的演示了model change的使用，該案例為一個簡單的拉伸模型，僅為演示用法，并無特殊的工程意義，讀者在使用生死單元的時候可根據自己的實際工況加以靈活選用。本案例中在STEP-2中移除了長方體模型中間的model change區域，然后在step-3中激活。
模型inp文件見附件：
ModelChange-AXin.zip
展開 基于多點位移控制增量的網殼結構穩定性分析
最終依據式（1）定義的約束方程如下：
為便于對比，采用同樣模型進行弧長法計算，僅僅是取消約束方程，而代替以實際的荷載施加，并用static，riks類型的step進行同樣的計算。
（三）計算結果：
網殼的最終屈曲變形形狀如下圖：
荷載位移曲線：
從上述結果可以看出，結構在此種荷載模式下的最大基底反力約為0.9e6N。同樣對該結構，采用常規的弧長法（static，riks）進行分析，與采用多點位移控制的曲線對比如下圖：
從計算結果可以看出，采用本研究的多點位移控制和弧長法計算結果基本一致，表明了采用約束方程實現多點位移控制的非線性屈曲的準確性。
計算環境與計算時長：
計算環境：AMD3900X 12核CPU，內存16G
計算時采用USE multiprocessors設為12，multiprocessing model采用MPI
計算時長：多點位移控制技術和采用弧長法計算時長均為200s左右。
結論：采用約束方程+位移加載可以有效地實現特定荷載分布下的單層網殼結構非線性屈曲分析，從而改善荷載增量法無法得到下降段的缺陷，是一種堪比弧長法的高效靜力非線性計算方法。同時，該方法相對于弧長法來說，其可以在該step后添加其他step，從而實現網殼結構屈曲后繼續承載下的結構分析，是一種高效的計算方法。
展開 NX NASTRAN 非線性
NX NASTRAN可以確定屈曲和后屈曲屬性，對于屈曲問題，可以同時考慮材料和幾何非線性；對于后屈曲問題，NX NASTRAN提供三種弧長法的自適應混合使用，包括Crisfield 法、Riks法及改進的Riks法,大大提高分析效率。
?
SOL601隱式求解用于進行靜態和隱式動態高級非線性分析。適合求解低速動態問題，如金屬成型、低速碰撞分析、地震響應及生物醫學等問題。
?
SOL701顯式求解用于顯式動態高級非線性問題的分析，適合高速動態問題，如爆炸響應、高速碰撞等波的傳播問題。
SOL601與SOL701的輸入非常相似，很多情況下，用戶可以從一個分析類型的結果開始RESTAR另一個分析，滿足不同目的的數據恢復和求解，重啟動可在穩定區和非穩定區的任何一點開始。
展開 薄殼圓柱屈曲分析
今天說一說屈曲分析
大致變形動圖是醬紫的
1 關于屈曲現象的簡單說明
當我們用腳踩可樂易拉罐的時候，罐子往往顯示為下面這樣
這個就是屈曲最常見的現象
原本不應該發生變形的物件在受壓的狀態下會出現失穩，這就是屈曲
在abaqus里面屈曲分析往往要經過兩個步驟
線性屈曲（特征值分析）+非線性屈曲（靜強度分析）
也就是先進行buckling分析，得到大致的失穩模態以及失穩發生時的力值，這時是與加載力值大小無關的；然后進行riks非線性屈曲，加上大致的失穩力值進行分析。
2 操作步驟
a 建立模型
這里用到薄殼模型，建模步驟我就省略了
賦予材料的話，要給于塑形，各種參數見上圖
b 建立相關連接
先建立兩個參考點（0，0，0）和（0，0，500）
再將頂部的所有點、底部的所有點與相應的參考點進行耦合（剛性）
c 建立分析步
d 加載邊界條件
由于與力值大小無關，因此在頂端參考點給-1的力，底部參考點完全約束
e 修改inp活著關鍵字，將節點的位移量進行輸出
f 計算求解，確保你的結果里面有fil文件（網格自行劃分，種子數給的10）
注意到屈曲力大概在1.2e6這么大
g 復制model-1到model-2
h 刪除原來的分析步，建立riks分析步
i 加載邊界條件
這里的力值參考最大屈曲力1.2e6，但是實際上這個值需要嘗試，特別是加載imperfection關鍵字之后，力值往往會急劇見效，這里設置為1.2e5
j 加入自定義關鍵字imperfection
在step之前加入
這里的imperfection是將特征值屈曲時的模態按照組合進行初始化，更接近實際結果
k 計算
展開 ANSYS中弧長法的原理
關于弧長法的原理，推薦參考《非線性分析弧長法的讀書報告》、陸新征老師學生時代的作業：《基于預處理技術和弧長法的非線性方程通用求解子程序總結報告 》，以及Yusd的博文《弧長法（Riks Method）的基本原理》。喜歡編程的話，還可以參考他的另一篇文章《弧長法(Riks method)通用求解程序》。英文資料可閱讀蘇黎世聯邦理工學院結構工程研究所Prof. Dr. Eleni Chatzi的講稿：《The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems》
在這里只強調一下弧長法的一個獨特的優勢：
見下圖，當微小荷載增量可以引起顯著的位移突躍，則成為荷載控制的急速通過(Snap though)；當微小位移增量可以產生顯著的荷載突躍，則成為位移控制的急速返回(snap back)——這種復雜的加載路徑，單純荷載或位移控制已經難以奏效，而弧長法則可輕松處理此類問題。
轉自公眾號——ANSYS學習與應用
旨在分享，若侵即刪.
展開 Nastran 常用的非線性分析類型簡介
更普遍的是，靜態平衡狀態的響應不穩定時可以采用修正的Riks法。該方法用于比例加載即載荷的大小是由一個單一的標量參數控制。該方法通過控制每個增量中的載荷-位移曲線的路徑長度來獲得平衡解(而不是控制載荷或位移增量),從而使載荷量成為系統的未知量。即使在復雜、不穩定的情況下，該方法也能得到問題的解。Riks法不能用于接觸、傳熱、耦合或有強迫運動等情形。
3. 蠕變、粘塑性和粘彈性行為
靜態分析中隨時間變化的材料響應可能涉及蠕變和膨脹(一般發生在相當長的時間段),或屈服應力與速率有關(在相當快速的過程中，例如金屬加工問題中，這一點通常很重要)。屈服應力與速率有關，使用常規的靜力學分析必須引進一個合適的時間尺度才能使SOL 400正確處理粘塑性。利用向后差分算子對塑性應變進行積分，蠕變問題以及粘彈性模型，由蠕變求解程序分析(這是由包含一個非零的時間間隔的NLPARM卡片指定)。非線性蠕變問題通常通過非彈性應變的向前差分積分(“初始應變”法)有效地解決，因為該算子的數值穩定性極限通常足夠大，因而可以在不多的時間增量步中得到求解結果。線性粘彈性模型由一個簡單的、隱式的、無條件穩定的算子積分。在這種情況下，自動時間步長策略是由用戶指定的精度容差參數控制的。它限制了一個增量步中最大非彈性應變率變化量。
4.
展開 
[案例]薄壁彎管在內壓和彎矩作用下的彈塑性坍塌分析
---對應elbowcollapse_s8r5.inp和elbowcollapse_s8r5_fine.inp
步驟 5：設置分析步與輸出請求
創建分析步： 創建static Riks分析步，用于施加彎矩和壓力載荷。
控制參數：打開大變形NLgeom:On。由于涉及坍塌（極值點失穩），通常需要使用弧長法（Riks） 或設置非常小的初始增量步0.05來控制求解過程。
場輸出請求： 確保輸出應力（S）、應變（E）、位移（U）等。
增加輸出請求： 輸出Nout點集合的施加彎矩一端的反作用力矩（RM）和轉角（UR），用于繪制力矩-轉角曲線、橢圓變形等。
步驟 6：定義相互作用
綁定約束： 使用“Tie”約束將彎管段與直管段的端面完全連接在一起。若用S8R5單元則只有一個零件，不需要。
對稱邊界條件： 在對稱面上定義對稱邊界條件（XSYMM）。一端全約束。
步驟 7：施加載荷與邊界條件
固定端： 約束遠離彎管的直管段末端的全部自由度（ENCASTRE）。
加載端： 在另一個直管段的末端，創建一個參考點（RP），并將該端面的所有節點與RP進行運動耦合約束（Kinematic Coupling），以模擬剛性端蓋。
載荷施加：
內部壓力： 在分析步中，作為表面壓力載荷施加在所有管道的內表面上，3.45 MPa。
面內彎矩： 通過在加載端的參考點（RP）上施加轉角（Rotation） 來間接實現彎矩加載。分別施加正向和反向的轉角來模擬“張開”和“閉合”彎矩。
步驟 8：劃分網格
對殼體部件進行網格劃分。為確保精度，特別是在彎管曲率較大的區域，需要進行適當的網格細化。
案例中使用了S8R5（8節點四邊形殼）單元，并進行了網格收斂性研究。
展開 abaqus鋼梁四點彎曲實驗
buckle分析?riks分析，模擬荷載位移曲線沒有下降段，該怎么修改啊
abaqus問答精華
A:利用riks分析，記錄位移和反力，再在xy plot operation里面作圖。
*step,INC=100, NLGEOM
*static,RIKS
0.001,1.0,,0.04,,4,3,-0.080
riks指的是弧長法
也可以利用*NODE PRINT命令，把節點位移和支座反力輸出到.dat文件，用excel畫圖，很方便，修改也容易。
*node print只能在ABAQUS/standard中起作用，在ABAQUS/explicit中不能使用
Q:有個問題想請教大俠：
我在做一個雙材料界面的斜裂紋問題，模型見圖。裂紋兩側是不同材料。如果不定義局部坐標，輸出的裂紋附近的正應力和剪應力應該都是整體坐標系下的，顯然這時候的S22,S12并不滿足連續條件。
現在就是想定義局部坐標系，并且輸出局部坐標系定義下的S22'，S12'，不知道abaqus具體怎么樣實現，請指教。
看到過幫助里提到過orient子程序，不知道怎么用
A:在simwe上壇友的幫助下（當時只記了他的mail:asia_999@sina.com,忘了他的ID，很抱歉），解決了這個問題。
1。在打開ODB文件的時候，前面的read_only前面的勾去掉。
2。create coordinate system選擇一種方法建立局部坐標。
3。coordinate system manager選擇建立的csys-->move to odb
4。重新打開odb，main menu bar-->result-->options...-->transformation-->user-spcified-->choose your csys-->apply
Q:adaptive mesh是不是只能用在explicit里頭呀
A:如果想要在standard里面用，在cae里該如何操作？
展開 abaqus非線性屈曲odb轉為vtu文件（v2017） ￥200
<div contenteditable="false" width="100%">
maxlpfindex = 10 #8階，求解階數
</div><div contenteditable="false" width="100%">
meshtype = "UnstructuredGrid" #支持六面體網格
</div><div contenteditable="false" width="100%">
job_name = "Job-2" # odb名字
</div><div contenteditable="false" width="100%">
workdir = "E:/test" #從# 案例路徑
</div><div contenteditable="false" width="100%">
step_name = "Step-1"
</div><div contenteditable="false" width="100%">
#調用方法
</div><div contenteditable="false" width="100%">
Nonlinear_vtu(workdir, job_name, step_name, maxlpfindex,meshtype)
</div><div contenteditable="false" width="100%">1）適合非線性屈曲分析將odb轉為vtu ,STE類型 ：static Riks</div><div contenteditable="false" width="100%">2）雙擊 test.bat</div><p><br></p>
展開 