abaqus變形系數的案例
批量提取Abaqus的節點坐標(初始坐標、指定Step下的變形量、變形后節點坐標) ￥40
<h2>摘要</h2><p>本文介紹如何使用Python腳本二次開發來批量提取ABAQUS輸出數據庫（ODB）文件中指定Step下的Set節點集變形量。通過詳細的步驟說明、代碼示例和圖片展示，您將學會如何使用該腳本，自動化輸出CSV文件包含（Node Label；Step Name、Increment、Step Time，U1，U2）。</p><p>如果還需要按Increment提取每個增量下的變形后的節點坐標的話，在提取變形量的基礎上，與初始坐標進行簡單的計算就可以求得坐標。 （備注：該代碼只提取了x，y方向的變形量）</p><h2>1. 問題描述</h2><p>在工程仿真和分析領域，提取ABAQUS輸出數據庫（ODB）文件中的節點集變形量是一項常見任務。然而，手動提取這些數據是一項繁瑣且容易出錯的工作。因此，需要一種自動化的方法來批量提取指定步驟下按節點集組織的變形量數據。</p><h2>2. 實例展示</h2><p>假設我們有一個名為`example.odb`的ODB文件，其中包含名為`Step-x`的步驟和名為`Set-x`的節點集。運行以上代碼后，腳本會自動將該步驟下節點集的變形量提取出來，并保存為`NodalDisplacement.csv`文件。
展開 橋梁橫向分部系數ABAQUS模擬
橋梁橫向分部系數計算方法有：杠桿法、剛性橫梁法、修正剛性橫梁法，鉸接板梁法、剛接板法和比擬正交法，其中剛性橫梁法用的較多，且重慶交院王老師編制了專門的計算程序，我采用ABAQUS模擬T梁，橫膈板采用剛性梁，用3D空間模擬，效果不錯，請大家鑒賞
axa.rar
橋梁橫向分部系數ABAQUS模擬
橋梁橫向分部系數計算方法有：杠桿法、剛性橫梁法、修正剛性橫梁法，鉸接板梁法、剛接板法和比擬正交法，其中剛性橫梁法用的較多，且重慶交院王老師編制了專門的計算程序，我采用ABAQUS模擬T梁，橫膈板采用剛性梁，用3D空間模擬，效果不錯，請大家鑒賞！
axa.rar
Abaqus安全系數的實現——UVARM子程序 ￥2
Abaqus安全系數的實現——UVARM子程序
在進行結構件強度校核有限元分析時，分析結果中安全系數的顯示能夠輔助對結構件進行優化設計。然而Abaqus軟件后處理器不能直接顯示安全系數，比較簡便的方式為調用UVARM子程序來實現。
UVARM子程序的代碼模板如下：
SUBROUTINE UVARM(UVAR,DIRECT,T,TIME,DTIME,CMNAME,ORNAME,
1 NUVARM,NOEL,NPT,LAYER,KSPT,KSTEP,KINC,NDI,NSHR,COORD,
2 JMAC,JMATYP,MATLAYO,LACCFLA)
INCLUDE 'ABA_PARAM.INC'
C
CHARACTER*80 CMNAME,ORNAME
CHARACTER*3 FLGRAY(15)
DIMENSION UVAR(NUVARM),DIRECT(3,3),T(3,3),TIME(2)
DIMENSION ARRAY(15),JARRAY(15),JMAC(*),JMATYP(*),COORD(*)
C The dimensions of the variables FLGRAY, ARRAY and JARRAY
C must be set equal to or greater than 15.
user coding to define UVAR
RETURN
END
該子程序中要在高亮顯示部分定義UVAR(NUVARM)的代碼，以便在后處理器中顯示安全系數場變量。
展開 
Abaqus 邊坡強度折減法計算安全系數 ￥5
Abaqus強度折減法計算邊坡的安全系數是采用設置場變量的方法，在分析計算過程中，逐步折減土體強度參數，當土體強度參數折減到很小時候，土體塑性區貫通，模型由于塑性變形過大無法計算下去的時候，這時候的場變量數值即為安全系數。
以某加固工程公路邊坡為原型，邊坡土體為黃土狀粉土，邊坡高度為11m，其主要物理力學性質列于表2.1，其中支護采用錨桿支護。
2.1 土體物理性質
土層名稱
厚度
m
重度
γ（kN/m3）
彈性模量
Mpa
泊松比
ν
粘聚力
c(kPa)
內摩擦角
φ（?）
黃土粉狀土
>30
16.8
12
0.3
15
22
定義場變量的地方為材料參數的第三列，number of filed variables設置為1，
定義兩個分析步，第一個分析步是重力場平衡，自重應力場，第二個分析步中進行折減。從菜單欄model—>edit kerwords進入到編輯關鍵詞界面，在第一個分析步開始之前添加如下關鍵詞，
在第二個分析步中添加如下關鍵詞。
其他建模步驟沒有特別需要注意的地方，完成這些就可以提交計算了，計算的得到邊坡安全系數為1.8，塑性區如下圖。
塑性區貫通
展開 ABAQUS-靜動態接觸分析中常用材料之間的摩擦系數
ABAQUS-靜動態接觸分析中常用材料之間的摩擦系數
ABAQUS-靜動態接觸分析中常用材料之間的摩擦系數.doc
Abaqus中利用fric_coef子程序定義摩擦系數與速度的關系
背景介紹
輪胎-路面摩擦模型在道路工程中應用十分普遍， 下圖為一典型的路面輪胎模型：
輪胎在路面正常前行時主要進行滾動，在有的文獻中（附件），摩擦系數描述為與滑動速度呈一定的函數關系，如下圖所示：
那么該如何在abaqus中實現摩擦系數隨滑動速度變化呢？答案就是fric_coef子程序
二。fric_coef子程序介紹
該子程序定義接觸面的摩擦屬性，其標準格式如下所示：
參考上述標準文件格式和幫助文檔的案例，即可完成摩擦系數的定義。
文獻中給出的結果如下：
附件文獻：
30-Al-Qadi、汪浩的胎路摩擦論文.pdf
最后，大家有相關需求可以關注“320科技工作室”的微信公眾號，更多干貨等你來撩~~
展開 Abaqus+PyQt+Python平面變形歐拉角計算
下面以簡單例子介紹平面變形、指向歐拉角的計算，包括絕對歐拉角、相對歐拉角。
1 簡化模型
下面的六面體為表面殼模型，下面由三段梁支持，三段梁分別沿X、Y、Z軸向。六個面的厚度不同，在上側3個面施加不同的壓力，如下左圖所示。位移云圖如下右圖所示。
2 計算要求
計算六面體上面3個面的變形歐拉角，包括3個面的絕對歐拉角，平面2、3相對與平面1的相對歐拉角。平面1、2、3如下圖所示。
3 數據處理
使用平面節點坐標、位移數據計算平面變形歐拉角。可以使用Python腳本輸出平面節點編號、節點坐標（X、Y、Z）、節點位移（U1、U2、U3），如下圖所示。下圖為平面1的10個工況的數據文件，打開的文本文件中7列數據為節點編號、坐標、位移。
三個平面10個工況的節點數據文件如下圖所示。每個文件中包含一個工況一個平面的節點編號、坐標、位移數據。
4 絕對歐拉角計算
使用PyQt+Python開發了一個簡單的小軟件，計算絕對歐拉角、相對歐拉角。
首先計算各平面的絕對歐拉角。
計算平面1的10個工況的絕對歐拉角。
平面1變形的絕對歐拉角計算結果如下圖所示。
伴隨絕對歐拉角計算結果，軟件同時寫出了平面變形前后的坐標系數據，如下圖。每行18個數據，每3個數據為一個坐標軸向量，變形前后2個坐標系，6個坐標軸，18個數據。
5 相對歐拉角計算
利用計算絕對歐拉角時得到的坐標系文件，計算平面變形相對歐拉角，如下圖所示，計算平面2相當于平面1、平面3相對與平面1的相對歐拉角。
計算結果如下圖所示。
6 小結
上述軟件用的算法申請了發明專利，軟件申請了軟著。
展開 Abaqus模擬橡膠大變形
Abaqus為用戶提供了多種本構關系來模擬超彈性材料，這種材料具有高度非線性，當Abaqus進行模擬時假設這種材料是具有彈性、各向同性，并且同時考慮幾何非線性效應。與材料的剪切柔度相比，對于大多數類似橡膠的固體材料，其可壓縮性非常小，當分析對象為平面應力問題、殼、薄膜、梁、桁架、或者鋼筋等，這個問題不值得關注。但是對于固體、平面應變或者軸對稱問題卻不能忽略。對此，Abaqus/Standard提供了雜交單元來模擬超彈性材料中完全的不可壓縮行為。
橡膠材料力學性能的描述方法主要為兩類：一類是認為橡膠為連續介質的現象學描述；另一類是基于熱力學統計的方法。基于連續介質力學的本構模型主要有Polynomial、Reduce Polynomial、Ogden模型等，其中Mooney-Rivlin模型是 Polynomial的特殊形式，Neo-Hookean 模型是Reduce Polynomial的特殊形式。基于熱力學統計主要有Arruda-Boyce和Van der Waals等本構模型。本文利用Abaqus模擬大變形的橡膠，具體步驟如下。
1、在Abaqus/CAE Sketch模塊中作出模型草圖，如圖1所示，然后在Part模塊中分別建立Push、Rubber、Base三個部件。其中Push為解析剛體，Base為離散剛體。
圖1 草圖
2、在Property模塊中定義橡膠的屬性，采用Mooney-Rivlin模型，參數如圖2所示，然后賦給Rubber部件。
圖2 橡膠參數設置
3、裝配，定義分析步，采用默認的場輸出和歷史輸出。為了保證剛開始能夠較容易收斂，設置分析步初始增量步為0.01，打開幾何非線性。
圖3 分析步定義
4、定義接觸對：Push下表面和橡膠表面，Base上表面和橡膠表面。
展開 ABAQUS CEL(例11) 地震工況下的邊坡大變形模擬 ￥70
ABAQUS CEL(例11) 地震工況下的邊坡大變形模擬
一、建模技術
地震工況下邊坡可能失穩進而出現滑坡現象，為避免模擬滑坡時網格產生的畸變問題，采用耦合歐拉拉格朗日法(CEL)進行滑坡的大變形模擬；土體本構采用摩爾庫倫模型；采用模型底部小范圍內的周期性荷載模擬地震荷載。
二、模型及部分結果展示
圖1：藍色為邊坡；紅色為空氣層
圖2：網格的劃分
圖3：賦予模型初始應力
圖4：土體達到地應力平衡時的應力分布
圖5：土體底部的地震荷載施加區域
圖6：所施加的周期性荷載（地震荷載）
圖7：邊坡因地震荷載產生的位移
圖8：地震波產生的區域
展開 ABAQUS 單向拉伸大變形模擬
靜態模擬一種軟材料POE的單向拉伸，拉伸應變希望到300%，但是總是在100%就失敗了。不知道哪里出了問題，有沒有高手幫幫忙。
HyperMesh導入Abaqus的變形后模型
HyperMesh導入Abaqus的變形后模型
經過Abaqus的操作可以將Abaqus計算后的變形單元和節點輸出為inp文件，這個inp文件可以通過HyperMesh讀取。
打開HyperMesh后設置user files為Abaqus，如圖1所示。
圖1
之后通過import Solve Deck導入Abaqus生成的inp文件，如圖2所示。
圖2
之后點擊import即可在窗口顯示導入的變形模型，如圖3所示。
圖3
在此基礎上，可以對網格進行編輯重劃分，修改結束后同樣可以導出修改后的模型inp文件，如圖4所示。
圖4
此時的inp文件記錄了在HyperMesh中修改后的模型。
再之后，可以通過Abaqus導入Hm生成的inp文件，如圖5所示。
圖5
導入后的模型在Abaqus中的顯式如圖6所示。
圖6
在HyperMesh編輯三維單元網格不是那么方便，建議的思路是先生成實體再進行網格編輯。
Abaqus6.12以后的版本據說可以由孤立的網格生成實體模型，這個還是比較方便的，生成實體模型后在Abaqus里面就可以進行網格重劃分。
展開 Abaqus薄板彎曲變形分析實
ABAQUS提供了業內領先的接觸建模能力，接觸中各種表面間的各類摩擦性質可以建立相應的模型模擬，來符合不同接觸行為的要求。
本文采用Abaqus/Standard求解器，進行薄板彎曲變形分析，用以簡單展示ABAQUS接觸建模及其分析功能。
1、 計算模型
如圖1所示，懸臂梁左端受剛性模具固定，右端受移動模具下壓產生變形。
2、 有限元模型
建立有限元模型，創建穩態分析步，分析薄板和剛性表面間的接觸，平板使用實體平面應變單元CPE4I, 該單元沿板厚方向只需要一個單元即可以準確模擬彎曲行為。剛性表面以解析剛性面模擬。
3、 接觸建立
ABAQUS中，接觸的一般需要三個步驟。
首先定義接觸表面。剛性表面一般作為接觸對的主面，本例中將剛性模具的面定義為主面，薄板面為從面。
進而定義接觸對。選擇發生接觸的主從面定義為接觸對。
最后定義接觸屬性。包括接觸類型，以及摩擦系數等相關接觸參數。本例選擇無摩擦的光滑接觸屬性。
本案例共包括三個接觸對，分別為三個剛性模具與薄板之間的接觸。
完成接觸設定后，對模型設定相關邊界條件：上下模具完全固定，沖頭向下移動60mm。薄板左端固定。
在此邊界條件下，沖頭向下移動時，薄板上的三個接觸對發生作用，使得薄板右端發生彎曲。
4、 接觸輸出
接觸設定中，對于多有表面的接觸信息，可以設定接觸應力、接觸位移等接觸輸出信息。
5、 分析結果
如圖所示，計算完成后薄板發生預想彎曲。案例設定了接觸應力輸出，接觸應力包括接觸壓力、摩擦剪切力的輸出，均可以在后處理中進行相應結果顯示。圖中所示云圖所示為接觸壓力云圖。
展開 ABAQUS模擬多道次變形的變量繼承方法
一、引言
使用ABAQUS進行多道次加工時，往往牽扯道次之間變量的繼承（如晶粒尺寸、累積損傷等），這對多道次變形模擬結果的準確性有較大的影響。本文以VUHARD子程序及簡單的熱壓縮模型為例，分享雙道次壓縮之間的晶粒尺寸的繼承方法。
ABAQUS變形放大比例Deformation Scale Factor
The Deformation Scale Factor options are in the lower left corner of the Basic page，等比例的話就在Uniform value 這個地方填寫，如果填1，說明是按實際變形大小來plot
