abaqus應(yīng)力驗(yàn)算的案例
【abaqus】個(gè)人筆記—應(yīng)力奇異&應(yīng)力平均&應(yīng)力集中
【abaqus】個(gè)人筆記—應(yīng)力奇異&應(yīng)力平均&應(yīng)力集中
ABAQUS熱應(yīng)力分析 附ABAQUS中初始地應(yīng)力的施加下載
軋輥與Cu層的熱傳導(dǎo)系數(shù)
下載地址：ABAQUS中初始地應(yīng)力的施加
Abaqus帶螺紋螺栓接觸應(yīng)力分析淺析 Abaqus帶螺紋螺栓接觸應(yīng)力分析淺析
目前的常規(guī)做法通常有兩種：1.簡化，用RBE2和beam梁來代替螺栓，這樣不能反映連接螺栓真實(shí)應(yīng)力，圖1為某結(jié)構(gòu)連接螺栓簡化的beam梁應(yīng)力云圖，沒有接觸應(yīng)力：
.直接做出來螺栓螺紋采用接觸分析，雖然得出的結(jié)果很精確，但這樣前處理工作量大（螺栓和螺紋用六面體網(wǎng)格建模）、計(jì)算量大（接觸收斂困難），如圖為某結(jié)構(gòu)帶螺紋螺栓和連接件模型（圖2）和計(jì)算得出的結(jié)果（圖3）：
圖3 計(jì)算結(jié)果
那么，有什么好辦法可以不用簡化帶螺紋螺栓，不用直接做出帶螺紋螺栓，又能得到足夠精確的結(jié)果？
運(yùn)用大型通用非線性有限元分析軟件Abaqus，只需要在接觸定義中設(shè)置跟實(shí)際螺紋形狀有關(guān)聯(lián)的參數(shù)，如牙角、螺距、螺栓小徑等，就可以模擬真實(shí)的連接螺栓接觸狀況。既可以得到足夠精確的分析結(jié)果，又節(jié)省了時(shí)間專注進(jìn)行其他的分析設(shè)置。如圖4，為連接螺栓接觸來定義帶螺紋螺栓：
圖4 連接螺栓接觸來定義帶螺紋螺栓
圖5為某結(jié)構(gòu)直徑10MM的帶螺紋的連接螺栓接觸壓力分布云圖：
圖5 某結(jié)構(gòu)直徑10mm帶螺紋的連接螺栓接觸壓力分部云圖
展開 ABAQUS-真實(shí)應(yīng)力和名義應(yīng)力轉(zhuǎn)化
ABAQUS-真實(shí)應(yīng)力和名義應(yīng)力轉(zhuǎn)化.doc
針對某袋除塵器整體進(jìn)行ABAQUS有限元分析，考慮九項(xiàng)載荷工況，分析設(shè)備靜應(yīng)力、熱應(yīng)力、變形及熱膨脹數(shù)值 ￥15
某袋除塵殼體結(jié)構(gòu)選型如下：
箱體板厚5mm
箱體角柱：角鋼L90*56*8
箱體加強(qiáng)筋：角鋼L90*56*6
花板厚6mm
花板下加強(qiáng)筋：橫向?yàn)楸怃?0*6，縱向?yàn)楸怃?00*6
箱體中間支撐管：鋼管Φ60*5
圖1 袋除塵殼體結(jié)構(gòu)示意圖
2、 建立模型
按照殼體結(jié)構(gòu)示意圖建立幾何模型如圖2所示。
圖2 建立幾何模型
三、約束條件及載荷
立柱底部約束如圖3所示。
圖3 立柱底部邊界約束
載荷：
（1）自重（軟件考慮）；
（2） 頂部載荷：檢修載（按400kg/m2）；
（3） 花板處載荷：濾袋、濾籠、濾袋積灰（積灰厚度按5mm）共3.06t；
（4） 灰斗積灰重：滿灰9.6t；
（5） 保溫載荷：按25kg/m2；
（6） 負(fù)壓11000Pa或正壓8000Pa兩種工況分別施加；
（7） 煙道及檢修平臺(tái)載荷：上煙道（出氣端）900kg，下煙道（進(jìn)氣端）
400kg，上中下三層檢修平臺(tái)檢修載荷均為400×2.85×3.25=3705kg。
注：此項(xiàng)載荷殼體和鋼支架各占一半。
（8） 灰斗卸灰口載荷（方向按照幾何模型坐標(biāo)系）：FX=4700N，F(xiàn)Y=3500N，F(xiàn)Z=-4700N，MX=3690N.m，MY=4800N.m，MZ=5540N.m。
（9） 頂部牛腿處檢修荷載：單個(gè)牛腿處載荷為1t，頂板為260×260，轉(zhuǎn)化為面壓添加，面壓為1×10×1000/260/260=0.148N/mm2。
下圖4所示為載荷添加圖示：
（a）負(fù)壓11000Pa （b）正壓8000Pa （c）花板處載荷
展開 ABAQUS中mises應(yīng)力云圖顯示的最大值還不到屈服應(yīng)力值為啥還有PEEQ值
ABAQUS中mises應(yīng)力云圖顯示的最大值還不到屈服應(yīng)力值為啥還有PEEQ值，PEEQ云圖有變形值
ABAQUS中求解某部分單元的平均應(yīng)力或平均應(yīng)變 ￥10
1、參考模型：單向纖維的RVE模型；
2、腳本功能：針對指定的單元集合，在后處理中求解平均應(yīng)力和平均應(yīng)變。
3、應(yīng)用的公式：一階均勻化計(jì)算方法。對于 RVE 模型的平均真應(yīng)力和平均真應(yīng)變，可通過對 RVE 內(nèi)每一個(gè)單元的真應(yīng)力 （真應(yīng)變）取均值獲得。使用一階均勻化計(jì)算方法輸出的應(yīng)力和應(yīng)變適用于各種邊界條件，但需要對每個(gè)單元進(jìn)行應(yīng)力（應(yīng)變）的輸出和計(jì)算。
通過Abaqus python腳本批量獲取節(jié)點(diǎn)的應(yīng)力 ￥25
背景
有限單元法計(jì)算單元積分點(diǎn)的應(yīng)力應(yīng)變值，而對于節(jié)點(diǎn)的應(yīng)力應(yīng)變值是通過外插得到的，Abaqus中云圖顯示的就是經(jīng)過插值和平均后的節(jié)點(diǎn)的值。通過工具欄的Query-Probe values可以查看單元或節(jié)點(diǎn)的應(yīng)力應(yīng)變等結(jié)果。
對于自動(dòng)化的后處理場景，通常需要自動(dòng)批量地獲取單元/節(jié)點(diǎn)的結(jié)果，通常都需要通過python腳本來實(shí)現(xiàn)。通過類似odb.steps['Step-1'].frames[-1].fieldOutputs['S']的場輸出可以比較方便地直接獲得單元的積分點(diǎn)應(yīng)力，但沒有直接的API可以獲取節(jié)點(diǎn)的應(yīng)力應(yīng)變等結(jié)果。
如果需要獲取部件表面節(jié)點(diǎn)應(yīng)力，可以通過創(chuàng)建路徑+XYData的方式實(shí)現(xiàn)，但想要獲得最大節(jié)點(diǎn)應(yīng)力，則該方式不便實(shí)現(xiàn)。
2. 通過python腳本獲取節(jié)點(diǎn)應(yīng)力結(jié)果
本文通過fieldOutput.getSub()函數(shù)獲取所有單元的節(jié)點(diǎn)結(jié)果，并對每一節(jié)點(diǎn)關(guān)聯(lián)的多個(gè)單元的節(jié)點(diǎn)值進(jìn)行平均后得到節(jié)點(diǎn)的結(jié)果。以下以某個(gè)簡單的odb結(jié)果進(jìn)行驗(yàn)證。
（1）批量獲得節(jié)點(diǎn)的mises應(yīng)力值
（2）批量獲得節(jié)點(diǎn)的X方向正應(yīng)力值
（3）批量獲得節(jié)點(diǎn)的最大主應(yīng)力值
（4）獲取節(jié)點(diǎn)的最大mises應(yīng)力及編號
3. 獲取節(jié)點(diǎn)應(yīng)變等結(jié)果
只需將腳本程序中的應(yīng)力場改為應(yīng)變成E等即可，此處不再演示。
以下為本文的python腳本代碼（代碼中作了必要的簡單注釋）。
展開 abaqus模擬圓孔結(jié)構(gòu)中應(yīng)力集中分析 ￥19.89
<p>1.問題描述 </p><p><br></p><figure style="text-align: center;" class="ql-align-center">
<figure class="figure-image" contenteditable="false" data-img="https://img.jishulink.com/202502/attachment/f337440f18204bb49994cddb8926c825.png" style="display: inline-block;">
<img src="https://img.jishulink.com/202502/attachment/f337440f18204bb49994cddb8926c825.png" data-mobile-src="https://img.jishulink.com/202502/attachment/f337440f18204bb49994cddb8926c825.png?image_process=/format,webp" data-pc-src="https://img.jishulink.com/202502/attachment/f337440f18204bb49994cddb8926c825.png?image_process=/format,webp" data-initial-src="https://img.jishulink.com/202502/attachment/f337440f18204bb49994cddb8926c825.png">
</figure>
</figure><figure style="text-align: center;" class="ql-align-center">
展開 ABAQUS 噴丸殘余應(yīng)力分析 ￥10
本案例適合哪些人學(xué)習(xí)：
1、學(xué)習(xí)型仿真工程師
2、理工科院校學(xué)生
3、對有限元分析感興趣的工程師
你會(huì)得到什么：
1、掌握噴丸三維模型的繪制
2、掌握顯示動(dòng)力學(xué)分析相關(guān)的材料參數(shù)設(shè)置
3、理解顯示動(dòng)力學(xué)分析步的建立
4、學(xué)習(xí)噴丸強(qiáng)化分析的相互關(guān)系的設(shè)置
5、了解顯示動(dòng)力學(xué)網(wǎng)格的劃分
6、學(xué)習(xí)結(jié)果后處理的查看與對比
案例介紹：
所使用軟件為ABAQUS2018.
本案例完整得提供了分析相關(guān)所有的分析文件。
?
ABAQUS實(shí)用子程序SPRINC提取主應(yīng)力
實(shí)用子程序SPRINC
在ABAQUS中使用UMAT子程序時(shí)有時(shí)會(huì)使用到最大主應(yīng)力進(jìn)行計(jì)算。通過查閱幫助文檔，ABAQUS實(shí)用子程序SPRINC可以在UMAT中計(jì)算最大主應(yīng)力和最大主應(yīng)變，SPRIND可以計(jì)算最大主應(yīng)力和最大主應(yīng)變的方向。
下面是ABAQUS幫助文檔關(guān)于實(shí)用子程序SPRINC的介紹：
SPRINC (calculate principal values)
Interface
CALL SPRINC(S,PS,LSTR,NDI,NSHR)
Variables to be provided to the utility routine
S
Stress or strain tensor.
LSTR
An identifier. LSTR=1 indicates that S contains stresses; LSTR=2 indicates that S contains strains.
NDI
Number of direct components.
NSHR
Number of shear components.
Variables returned from the utility routine
PS(I), I=1,2,3
The three principal values.
展開 
ABAQUS應(yīng)力場
求大神告知abaqus中考慮計(jì)算出的應(yīng)力場，是通過預(yù)定義場中的Initial-Mechanical-Stress，還是Step 1-Other-Field？
ABAQUS有效應(yīng)力
Us=AB*(Δσ1-Δσ3)/2
而對于孔壓系數(shù)A主要反應(yīng)剪切力對土體積變化影響，它在剪切過程中不是常數(shù)，與加載應(yīng)力應(yīng)變數(shù)率，加載卸載，排水條件和應(yīng)力歷史都有關(guān)系，一般實(shí)驗(yàn)可以測得。
孔隙水壓力增量u=Uc+Us=B*((Δσ1+Δσ3)/2+A*(Δσ1-Δσ3)/2)
③外力作用下飽和土土體強(qiáng)度
太沙基一維固結(jié)有助于理解有效應(yīng)力和孔隙水壓力的變化過程。
基本假定：
1.土層是均質(zhì)且完全飽和
2.土顆粒與水不可壓縮
3.水的滲出和土層壓縮只沿豎向發(fā)生
4.滲流符合達(dá)西定律且滲透系數(shù)保持不變
5.壓縮系數(shù)是常數(shù)
6.荷載均布，瞬時(shí)施加，總應(yīng)力不隨時(shí)間變化
太沙基一維固結(jié)模型：
有限元網(wǎng)格：
頂面排水，土體表面一次瞬時(shí)施加荷載。
分析結(jié)果：
①不同高度處孔壓隨時(shí)間的變化。
②固結(jié)度與時(shí)間因數(shù)理論值對比
理論解：
對比圖：
③超孔隙水壓力與有效應(yīng)力對比
轉(zhuǎn)自公眾號——abaqus有限元仿真
旨在分享，若侵即刪.
展開 ABAQUS地應(yīng)力平衡
對于涉及開挖、填充的動(dòng)態(tài)巖土工程問題，初始應(yīng)力場是正確模擬其施工過程的先決條件。初始應(yīng)力的加載必須滿足地應(yīng)力平衡，而地應(yīng)力平衡其實(shí)就是為了使數(shù)值模擬獲得一個(gè)存在初始應(yīng)力，而無初始應(yīng)變的狀態(tài)。當(dāng)考慮自重場是產(chǎn)生地應(yīng)力場的主要因素時(shí)，顯然重力是外力，初始應(yīng)力場是內(nèi)力，將提取出的內(nèi)力施加于數(shù)值模擬后再施加重力，是內(nèi)力和外力平衡，從而獲得較精確的沒有受到人為因素干擾情況下的數(shù)值模型的初始狀態(tài)。
在 ABAQUS 有限元軟件中，提供了 5 種定應(yīng)力平衡方法，這里我就比較常用的，使用比較普遍的一種進(jìn)行講解。用的例子還是上節(jié)ANSYS的例子。
接下來，小編就一個(gè)簡單的算例就行地應(yīng)力平衡。在一個(gè)均質(zhì)的地基上有一根立柱。
柱子的尺寸為1m×2m×10m，密度2400kg/m3，彈性模量2.8E8Pa（這里我為了方便觀察規(guī)律，降低這個(gè)彈性模量），泊松比0.167；均質(zhì)地基的尺寸為30m×30m×15m，密度為2000kg/m3，彈性模量3E7Pa，泊松比0.3。模型見圖1.
1、將模型導(dǎo)入ABAQUS中。
2、創(chuàng)建靜力荷載步
3、殺死柱子單元
4、施加約束（四周施加法向連桿約束，底部施加全約束）
5、施加自重荷載
6、提交計(jì)算
7、查看后處理結(jié)果
豎向位移
豎向應(yīng)力
8、應(yīng)力場輸出
9、在工作目中找到stress_out.rpt文件
10、經(jīng)過格式的轉(zhuǎn)化（文件的后面不要有空行）
11、在之前的inp文件中Step-1前面加上*Initial conditions,Type=stress,Input=stress_in.rpt這句命令，記得一定要把stress_in.rpt放在工作目錄里面。
展開 基于Abaqus的工字梁靜應(yīng)力分析
對于工字梁鋼架結(jié)構(gòu)，通常都需要對鋼梁的應(yīng)力及撓度進(jìn)行分析，本案例通過Abaqus有限元分析軟件對此材料力學(xué)的問題進(jìn)行了分析，相比利用材料力學(xué)知識的理論計(jì)算，不僅結(jié)果計(jì)算準(zhǔn)確簡潔還能將結(jié)果可視化呈現(xiàn)，這也是本文采用Abaqus仿真分析的一大優(yōu)勢。
2模型建立
本案例采用的工字型鋼梁的幾何尺寸如圖1所示，其中鋼梁的材料特性如表1所示。
圖1工字型鋼梁的結(jié)構(gòu)尺寸
表1工字型鋼梁的材料特性
彈性模量E（N/m2）
泊松比μ
屈服強(qiáng)度fy（N/m2）
2.1E11
0.3
3.45E8
3.有限元建模分析
3.1單位制
這里不論是否使用abaqus軟件，只要是使用有限元仿真程序，均需要首先明確建模分析使用的單位制，保證單位制的統(tǒng)一。單位制的統(tǒng)一根據(jù)牛頓定律（F=ma）進(jìn)行相互推導(dǎo)確認(rèn)，這里不再贅述。
3.2ABAQUS建模
啟動(dòng)ABAQUS/CAE，進(jìn)入分析界面，依照有限元分析一般流程步驟：前處理→求解→后處理。我們通過工具欄的module模塊選擇part模塊進(jìn)行工字梁幾何建模，按照圖1所示的幾何尺寸建立的模型如圖2所示。接著在property功能模塊定義工字梁的材料本構(gòu)和截面屬性，并將截面屬性賦予相應(yīng)的區(qū)域上。這里鋼梁的材料本構(gòu)選用理想的線性彈塑性本構(gòu)，即將切變模量（塑性應(yīng)變）設(shè)置為0完成材料本構(gòu)設(shè)置。應(yīng)當(dāng)注意要設(shè)置截面的屬性并給部件賦予上截面屬性。這是不同于ANSYS軟件分析流程的地方，ABAQUS不能把材料屬性直接賦予模型，而是先定義模型的截面，將材料屬性定義在截面上。之后定義裝配件并設(shè)置分析步。
展開 