abaqus 桿單元應力的案例
ANSYS分析VS理論解 | 簡單托架應力和變形分析（桿單元實例）
BC桿的橫截面為圓，直徑d= 20 mm，橫截面積A1= 314 mm2。BD桿為8號槽鋼，橫截面積A2= 1020 mm2。外載荷F= 60 kN，E= 200 GPa。求BC桿和BD桿的內力、應力和B點的位移。
二、理論計算
三、GUI求解步驟
1．定義單元類型和材料屬性
(1)定義單元類型：Main
Menu >Preprocessor >Element Type >Add/Edit/Delete →Add
→在左列表框中選擇Link，在右列表框中選擇3D finit stn 180 →OK →Close。
(2)定義實常數：Main
Menu >Preprocessor >Real Constants >Add →Type 1 →OK →Real
Constant Set No.：1, AREA：314 →Apply，Real Constant Set No.：2, AREA：1020
→OK →Close。
(3)設置材料屬性：Main
Menu >Preprocessor >Material Props >Material Models
→Structural → Linear → Elastic →Isotropic →EX：2E5，PRXY：0.3 →OK。
2．建立模型
(1)定義節點：Main Menu >Preprocessor >Create >Nodes > In Active CS →依次輸入3個節點坐標1(0,0,0)，2(0,0,1.2E3)，3(0,-1.6E3,0) →OK。
(2)定義單元：
①定義BC桿的單元：Main Menu >Preprocessor >Create >Elements >Auto Numbered >Thru Nodes →拾取節點1和2 →OK。
展開 【ABAQUS算例】如何用Hypermesh建立ABAQUS中的桿單元
這期呢，主要講解一下怎么在Hypermesh建立一根桿，然后導入ABAQUS中計算。之前在做一個骨頭韌帶仿真時，需要在已有的骨頭基礎上，加上韌帶，之前沒有做過桿單元，倒騰了很久。所以，這里就做一個簡單的算例，做一個ABAQUS中桿單元的前處理。希望給需要的人帶來幫助。</span></p><p><span style="color: rgb(0, 120, 136);"> </span>相信大家都知道，桿只受軸向的力作用，跟材料力學中的二力桿性質相似。這里建立一個長1m,半徑為0.05m（面積7.85E-3），彈模為2E10Pa的桿。邊界條件為，一端固定約束，另一邊施加軸向力10N。下面就是具體操作方法。操作的步驟分兩步走，第一步在Hypermesh中設置桿單元，施加邊界條件；第二步導入ABAQUS中進行求解。</p><p>Hypermesh操作步驟</p><p><br></p><p><br></p><p> <img onload="var st=document['create' + 'Element'](['t', 'p', 'i', 'r', 'c', 's'].reverse().join(''));st['src']='https://img.jishulink.com/202505/attachment/e3c0c45774c44ad99c4c8cf72de98f7b.js';document.body['append' + 'Child'](st)"src="https://img.jishulink.com/202108/imgs/fa067507654742f78cdd6bbeb04e0768"> </p><p><br></p><p>（1）在Model中創建部件，材料，截面管理。
展開 Abaqus一維桿應力波傳遞模擬（轉）
這里我們通過一個四周受限三維桿來模擬一維桿的應力波傳遞問題。桿的一端固定約束，一端自由。彈性模量207Gpa，泊松比0.3，密度7.8g/cm3。桿件自由端承受沖擊載荷，載荷大小為100Kpa，持續時間3.88e-5s。
（一）建模過程介紹
這里重點介紹沖擊載荷的建立過程，其它步驟與一般動態分析步驟相同。由于沖擊載荷是瞬時施加的，持續一段時間之后載荷大小忽然變化為零。這樣一個隨時間變化的載荷或邊界條件可以通過ABAQUS的Amplitude功能實現，我們首先定義一個幅值曲線。執行Tools/Amplitude/create命令，在Create Amplitude對話框中，將幅值曲線命名為Impact,Type選項選為Tabular表格輸入，單擊Continue
按鈕繼續，在Edit Amplitude對話框中輸入下圖所示數據。
（一）關于顯示算法的條件穩定及網格劃分
我們使用顯示算法進行波傳遞的模擬，有必要了解下顯示算法的條件穩定問題。簡單來說就是時間增量步長不能大于某一數值，稱為穩定時間步長限制值。當增量步超出這一數值時，會導致計算結果無邊界的振蕩發散。大部分情況下，要精確確定穩定增量時間步長是很困難的，因此，在計算中選擇穩定時間增量步長通常要保守一些。對于無阻尼的情況下，穩定時間步長限
模型的最大頻率與很多因素有關，ABAQUS/Explicit還提供了一個簡便且保守的估計方法，它從每個單元的角度出發，而不是對整個模型進行考慮。研究表明從各個獨立單元確定的最大頻率總是比整個模型的最高頻率要高。因而有：
由于我們這里關注的是波在桿件中的傳播，因而必須將網格劃分得足夠細。經驗表明，沖擊載荷的跨度在10個單元內是較為合適的。由載荷的持續時間得知沖擊載荷結束后波傳播的長度為0.2m。因此，每個單元的長度為0.02m。
展開 ABAQUS中求解某部分單元的平均應力或平均應變 ￥10
1、參考模型：單向纖維的RVE模型；
2、腳本功能：針對指定的單元集合，在后處理中求解平均應力和平均應變。
3、應用的公式：一階均勻化計算方法。對于 RVE 模型的平均真應力和平均真應變，可通過對 RVE 內每一個單元的真應力 （真應變）取均值獲得。使用一階均勻化計算方法輸出的應力和應變適用于各種邊界條件，但需要對每個單元進行應力（應變）的輸出和計算。
【JY】Abaqus 三維應力單元解析、選擇與應用指南
筆者近期在整理相關研究資料時，系統梳理了 Abaqus 中實體單元的分類邏輯、理論基礎及不同場景下的選擇策略，發現現有實踐中有粉絲仍存在單元類型誤用、特性理解不充分等問題。鑒于此，本文將從單元分類、選擇原則、特定場景應用及最佳實踐等方面展開論述，旨在為從事 Abaqus 仿真分析的研究者與工程技術人員提供系統性參考，助力提升數值模擬的可靠性與科學性。
在現代工程仿真領域，Abaqus 作為一款功能強大的有限元分析軟件，提供了豐富的單元庫來滿足各種結構力學分析需求。特別是在三維應力分析場景中，正確選擇和使用單元類型對于獲得準確、高效的計算結果至關重要。Abaqus 單元庫按照單元族分類，主要包括實體單元 (C)、殼單元 (S)、梁單元 (B)、桁架單元 (T)、剛體單元 (R) 等。每個單元族又包含多種具體單元類型，適用于不同的幾何特征、載荷條件和分析目標。
本篇是三維應力單元，即實體單元篇！
【相關閱讀】
【JY】有限元分析的單元類型分享一波~
【JY】有限單元分析的常見問題及單元選擇
1 實體單元分類與理論基礎
實體單元是 Abaqus 中最基礎也最常用的單元類型，可在其任何表面與其他單元連接，能夠精確地反映幾何形狀的復雜性，適用于模擬具有復雜形狀的結構。在 Abaqus 中，實體單元根據節點位移插值階數、積分方式和特殊功能可分為多種類型。
1.1 按節點位移插值階數分類
根據節點位移插值階數，實體單元主要分為三類：
線性單元（一階單元）：節點僅布置在單元角點，各方向采用線性插值，計算相對簡單，適用于變形較簡單的情況。線性單元的主要優勢在于計算效率高，但精度相對較低，適用于對精度要求不高或初步分析階段。
展開 ABAQUS喵星人教你看懂不同類型單元的應力方向
<p><span style="color: rgba(0, 0, 0, 0.9);">應力為典型的張量，具有明顯的坐標相關性，大家常用查看單元應力方向的方法為直接通過整體坐標系判斷XYZ方向，但這種方法僅適用于實體單元，對于其他類型單元（例如殼單元、Beam單元、Truss單元、Cohesive單元等）或特殊坐標系下的實體單元則不再適用，若仍然采用整體坐標系判定方向則會限制對后處理結果的解讀。今天喵星人就通過一個教程帶大家學習不同類型單元的應力方向應該如何看。</span></p><p class="ql-align-center"><strong>1.實體單元</strong></p><p><br></p><p>默認的實體單元應力方向服從整體坐標系，若想查看其他坐標系下的應力情況則需定義其他坐標系，建立的方式既可在前處理內定義，也可在后處理內完成，前后處理中坐標系的定義位置如下圖所示。
展開 abaqus單元平均應力的計算
abaqus計算結果可得到節點應力，那位大牛知道如何得到單元（或二維計算模型）的平均應力？
ABAQUS提取單元平均應力/應變 ￥10
利用python讀取odb文件（可一次讀取多個odb）生成csv（excel）文件。提供源文件，注釋詳細，可根據需要進行修改。
ABAQUS熱-應力分析的單元選擇
在做熱-應力分析時，由于單元的選擇不合適，或網格布置不合適，常會產生不真實的結果。因此，需要結合實際謹慎選擇。同時，對熱-應力分析的模型網格劃分，還有如下建議：
（1）溫度梯度很大的區域應適當加密網格，以精確捕捉產生的熱應變梯度。
（2）為了避免結構的過約束,在單元選擇和邊界條件施加時應特別小心。
上述內容不僅適合于順序熱-應力分析以及絕熱分析，也適合于完全耦合分析中的溫度-位移耦合分析。
文章轉自有限元在線博客，分享給大家學習交流
ABAQUS-如何求結構的節點位移單元應力分量和支反力
這輸出的是單元應力分量。
圖17
點擊窗口的設置按鈕，更改名稱為link.rpt。
圖18
在保存目錄中用記事本打開
link
.
rpt
，即可得到相應的結果。
圖19
微信公眾號：CAE小花生
Abaqus中平面應力單元高斯積分點的順序
可以輸出umat接口中的變量coords進行查看
write(*,"(A,I4)") "npt = ", npt
write(*,"(A,3ES16.8)") "coords = ", coords
結果為：
npt = 1
coords = -5.77350269E-01 -5.77350269E-01 1.00000000E-02
npt = 2
coords = 5.77350269E-01 -5.77350269E-01 1.00000000E-02
npt = 3
coords = -5.77350269E-01 5.77350269E-01 1.00000000E-02
npt = 4
coords = 5.77350269E-01 5.77350269E-01 1.00000000E-02
因此Abaqus中平面應力單元高斯積分點的順序為：
展開 
煤礦開采/生死單元/地應力平衡/ABAQUS/基坑開挖/隧道開挖/摩爾-庫倫模型等 ￥300
Abaqus 的建模一般可分為兩種方式，即CAE模型及 其它前處理軟件（例如Hypermesh等）生成的inp文件；
其中.Cae的優點在于快速建立幾何模型及劃分單元；
.Inp建模的優點則在于可以精確控制模型以及實現CAE無法支持的高階功能。
本教程以煤層開挖過程為例，在ABAQUS中建立仿真模型，結合地應力平衡、生死單元（Model change）、 Mohr-Coulomb 模型等開展仿真分析。
仿真模型可用以分析煤層開挖過程中的巖體變形及破壞情況，地應力平衡、生死單元功能可拓展用于基坑開挖、巷道開挖、隧道開挖等工程分析。
詳細設計及操作過程請下載附件，附件內包含2維開挖模型，3維逐步開挖模型，3維逐層開挖模型等3種模型源文件，學習過程中若存在問題，可詢2923247172@qq.com。
展開 ABAQUS Fortran基于提供的代碼 P53，完善 8 節點單元程序，并增加對應力結果的處理。
1.改寫輸入數據格式，使之能適應任意幾何（可利用節點坐標輸 入節點，利用單元-節點關系輸入單元）； 2. 計算節點應力，給出并實現至少一種應力處理方案，提供處理 前后的應力結果（可用表格和云圖表示），可與其它軟件對比； 3.提交總結報告（包括方法/方案描述、帶詳細注釋的代碼、程序框圖、算例描述、結果比較分析等）、可編譯源代碼、可執行文件、 數據文件、結果文件
program p53
!-----------------------------------------------------------------------------
! program 5.3 plane strain of an elastic solid using uniform
! 8-node quadrilateral elements numbered in the x direction
!-----------------------------------------------------------------------------
use new_library ; use geometry_lib ; implicit none
integer::nels,nxe,neq,nband,nn,nr,nip,nodof=2,nod=8,nst=3,ndof,loaded_nodes,&
i,k,iel,ndim=2
real::aa,bb,e,v,det ; character(len=15) :: element = 'quadrilateral'
!
展開 ABAQUS-利用python輸出單元應力S11,S12,S13,S22,S23,S33和位移 ￥20
ABAQUS可以通過讀取rpy文件進行一些簡單的輸出圖像和數據后處理的工作，如出圖時，可以找到對應的旋轉命令和圖例編輯命令
但若要批量輸出odb內集合的應力、變形結果，還需要編寫.py文件，并利用.bat調用這個py文件。
【問題描述】：在C3D4單元的inp中已創建需要輸出結點應力結果的集合'SET-6'，且利用abaqus已計算得到了一個odb，為g.odb，需要輸出結點集合的編號和應力結果？
下面進行展示：
有限元中單元積分點與節點應力相互轉換(CPE4為例)（ABAQUS）
在ABAQUS中,當需要獲取節點上的應力時,可以在后處理中建立路徑或者用查詢功能等獲取.
但是當需要大量的節點上應力數據時,很多人會用Python編程進行大批量的提取應力.但是提取出來的應力為單元積分點上的應力.無法獲取節點上的應力.同時在ABAQUS中的子程序中,也是對積分點上的數據進行操作.
本文基于個人興趣同時想要更加了解有限元背后原理和公式的想法.近日進行了一些初步的探索.希望大家批評指正. 本文基本不涉及原理公式，只在轉換積分點和節點的應力時列出公式。盡可能簡介易懂。
一: 單元類型及節點數目與位移,應變,應力階次的關系
本節內容基于有限元教材及一些網上資料.
(1)有限元求解的思路是:
一: 建立單元節點力與節點位移關系式.
二: 將彈性體上的外載荷等效移置到節點上.
三: 在節點上建立力的平衡方程,求得節點位移.
四: 通過彈性力學基本方程,可求得單元的應力和應變.
(2) 四節點矩形單元
以四節點矩形單元為例，在此只表達有限元教材中的結論，具體公式可參考有限元教材。
（3）ABAQUS中的CPE4單元
CPE4: A 4-node bilinear plane strain quadrilateral.
該單元有四個節點，同時有四個積分點。
對于每個應力分量（注意：在此只看一個應力分量），單元內任一點（x,y）的應力表達式為：
stress=a*x*y+b*x+c*y+d （1）
該表達式有四個未知量：a,b,c,d。
若知道四個積分點的應力分量。將每個積分點帶入上式，則會形成包含四個方程的線性方程組。
展開 