abaqus變形很小的案例
試驗鐵地板：精度高、變形小，測試臺面用得久
這個過程會徹和底消除鑄件內部的應力，確保平臺在長達數年的使用中不會發生絲毫變形，精度保持相當穩定。
核心二：精度標準——從“毫米”到“微米”的躍升
試驗鐵地板精度遠高于普通工業地板，它的精度是保證試驗數據可重復、可信賴的基礎。
平面度：這是比較核心的指標。普通平臺可能滿足于1mm/m的誤差，但試驗鐵地板通常要求達到 0級精度，即每平方米平面度誤差不超過 0.05毫米（約一根頭發絲的直徑）。
刮研工藝：要達到如此高的平面度，必和須依靠傳統的人工“刮研”手藝。工人們在平臺表面反復涂抹色漿，用刮刀一剔除高點，比較終使表面布滿均勻的接觸斑點（每平米超過25個點）。
水平度：安裝調試后，整個臺面的水平度誤差要求控制在 0.03mm/m 以內，確保設備運行不會產生傾斜偏差。
核心三：T型槽系統——靈活固定的關鍵
試驗鐵地板表面布滿平行的 T型槽，這是它功能性的體現。
作用：你可以利用T型槽和配套的螺栓，快速、牢固地固定發動機、測功機、振動臺等各種形狀各異的試驗設備，而無需在地面上打孔。
組成：一個完整的試驗鐵地板通常由平臺本體、T型槽、蓋板、地腳螺栓和專門的調平裝置（如墊鐵）組成。
試驗鐵地板的四種常見應用布局
熱模擬實驗室型：
主要應用場景：材料高溫、高壓試驗。
特點說明：側重耐熱性和熱穩定性，減少溫度變化帶來的變形。
單缸機實驗臺型：
主要應用場景：小型發動機或單體部件測試。
特點說明：平臺尺寸相對較小，但精度要求高，便于快速拆裝。
多缸機實驗臺型：
主要應用場景：大型多缸發動機、變速箱測試。
特點說明：平臺面積大，承重能力相當強，通常設計有復雜的油路和排水系統。
汽車疲勞振動試驗型：
主要應用場景：模擬車輛長期顛簸、振動環境。
展開 批量提取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`文件。
展開 小變形下的非線性彈性行為和其 UMAT實現——簡單的鄧肯張模型 ￥4.9
文檔介紹了非線性彈性行為的背景，鄧肯張模型的由來，和UMAT實現的代碼，展示如下：
【網格優化功能：Abaqus 極小曲面】之 極小曲面 III --免安裝綠色小程序2 ￥99
image_process=/format,webp/quality,q_40/resize,w_760" data-initial-src="https://img.jishulink.com/upload/202108/184d29023c7a4cf29da8947060e5ba28.jpg">
</div>
</div><p><br></p><p>前段時間跟大家分享了怎么用MATLAB 、 python或用自編的小程序 建立 3D打印用的極小曲面及將其輸出為stl格式的方法，具體請看：</p><ol><li><a href="https://www.yqgqt.org.cn/content/post/b9ec543f-74f1-4dda-add4-17c0deb4f303" rel="noopener noreferrer" target="_blank">Matlab生成極小曲面</a>，包括matlab腳本及生成為stl的腳本</li><li><a href="https://www.yqgqt.org.cn/content/post/1802096" rel="noopener noreferrer" target="_blank">python生成極小曲面</a>，包括python腳本、安裝包及生成stl的腳本</li><li><a href="https://www.yqgqt.org.cn/content/post/1812725" rel="noopener noreferrer" target="_blank">免安裝綠色小程序</a> ，是一個小程序，能夠生產帶 數學表達式的曲面；同時<strong>沿著法向偏移實現加厚</strong>；</li></ol><p>這些方法生成的曲面輸出的stl網格文件，一般網格質量較差；其實若做有限元仿真可以用其網格工具進行優化即可
展開 
【Abaqus 3D打印建模】之 極小曲面 III --免安裝綠色小程序 ￥79
<h1><strong>***已更新，請見【網格優化功能：Abaqus 極小曲面】之 極小曲面 III --免安裝綠色小程序2***</strong></h1><p>***已購買本算例的，可以聯系我免費更新***</p><h2 class="ql-align-center"><strong>0.前言</strong></h2><p>前段時間跟大家分享了怎么用MATLAB 和 python 建立 3D打印用的極小曲面及將其輸出為stl格式的方法，具體請看：</p><ol><li><a href="https://www.yqgqt.org.cn/content/post/b9ec543f-74f1-4dda-add4-17c0deb4f303" rel="noopener noreferrer" target="_blank">Matlab生成極小曲面</a>，包括matlab腳本及生成為stl的腳本</li><li><a href="https://www.yqgqt.org.cn/content/post/1802096" rel="noopener noreferrer" target="_blank">python生成極小曲面</a>，包括python腳本、安裝包及生成stl的腳本</li></ol><p>以上兩種方法基本上等效，不僅可以生產極小曲面，也能夠<strong>加厚</strong>或<strong>輸出為stl</strong>，只不過是所用軟件不一樣。但據網友反饋，存在寫不方便之處，如有些人需要安裝matlab，或python的依賴包安裝失敗等。
展開 【Abaqus 3D打印建模】之 極小曲面 I --Matlab生成極小曲面 ￥79
*cos(x)
[x,y,z]=meshgrid(-i*pi:0.05*pi:j*pi); %該區域中所有的 xyz坐標
v=f(x,y,z);c=isosurface(x,y,z,v,0) %找到v=0的等值面
h=patch(c);
isonormals(x,y,z,v,h);
set(h,'facecolor','g','edgecolor','b')
grid oncamlight;lighting gouraud;
</pre><p>其實只有前四行代碼是用來生成極小曲面的，后幾行是視覺效果之類的。</p><p>那么，對于做仿真的人來講，只有圖片肯定不行，至少得轉換成stl格式或obj格式或inp等 abaqus支持的格式，此時需要把matlab生成的頂點、面、連接順序等信息寫入為 stl或obj格式或inp（abaqus可支持obj導入）。
展開 【Abaqus 3D打印建模】之 極小曲面 II --python生成極小曲面 ￥79
***請注意，附件僅為現式和隱士極小曲面生成和輸出為stl的python代碼，并不包括加厚***
***有加厚的需求，請看其他算例。加厚后輸出的stl，也是空心封閉的殼，不是實心的實體***
0.算例
上一個帖子介紹了怎么用matlab建立極小曲面，詳情見Matlab創建極小曲面。
下面是個簡單的算例，在y方向壓縮極小曲面之Gyroid，幾何模型建立方法見下文，建立后陣列并有畫網格導入abaqus即可。
為了對比該極小曲面的應力水平，采用同樣的材料做了單軸壓縮，兩種情況對比如圖所示：
從圖中可以看到，如果僅去極小曲面上的一個點作為其應力應變，其曲線甚至比實心立方體還高，顯然這是不合理的。出現這種現象的主要原因是，此類細觀結構或變形不均勻時，不能取一個點代表整個模型，不然很可能會出現較大的誤差。對此問題，細觀力學有些方法，比如作用力反作用力法，體積平均法等，但也有人認為對于細觀結構，作用力反作用力法 不太合理。體積平均法的簡單表達式如下：
即模型中每一個單元的應力（應變）對單元體積積分后，除以模型整個體積。上述應力應變曲線也證實，采用該方法能夠得到較為真實的數據。
那么，怎么通過體積平均法獲得數據呢？ 在計算結束后，需要通過python腳本對數據進行處理，輸出中需要EVOL（單元體積）以便獲得各個單元的體積。
1.介紹
之前已經介紹過什么是極小曲面，同時根據表達式定義為 隱式 和顯式 極小曲面，主要區別是 隱式極小曲面一般只有一個方程，不容易將x,y,z獨立表達出來，如下圖所示：
今天介紹怎么用python生成上述極小曲面并輸出為stl文件。
展開 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應用中的小技巧 附ABAQUS field output 中frequency 不同選項
對于三維結構，只有模型區域滿足以下條件，才能被劃分為結構化網格：
①沒有孔洞、孤立的面、孤立的邊、孤立的點；
②面和邊上的弧度值應該小于90°；
③三維區域的所有面必須要保證可以運用二維結構化網格劃分方法；
④保證區域的每個頂點屬于三條邊；
⑤必須保證至少有四個面（如果包含虛擬拓撲，必須僅包含六條邊）；
⑥各面之間要盡可能地接近90°，如果面之間的角大于 150°，就應該對它進行分割；
⑦若三維區域不是立方體，每個面只能包含一個小面，若三維區域是立方體，每個面可以包含一些小面，但每個小面僅有四條邊，且面被劃分為規則的網格形狀。
>
>
>
>
三維結構的掃略網格劃分-掃略網格劃分技術
對于三維結構，只有模型區域滿足以下條件，才能被劃分為掃略網格：
①連接起始面和目標面的的每個面（稱為連接面）只能包含一個小面，且不能含有孤立的邊或點；
②目標面必須僅包含一個小面，且沒有孤立的邊或點；
③若起始面包含兩個及兩個以上的小面，則這些小面間的角度應該接近180°；④每個連接面應由四條邊組成，邊之間的角度應接近90°；
⑤每個連接面與起始面、目標面之間的角度應接近90°；
⑥如果旋轉體區域與旋轉軸相交，就不能使用掃略網格劃分技術；
⑦如果被劃分區域的一條或多條邊位于旋轉軸上，ABAQUS/CAE不能用六面體或楔形單元對該區域進行掃略網格劃分，而必須選擇 Hex-dominated形狀的單元；
⑧當掃略路徑是一條封閉的樣條曲線時，該樣條曲線必須被分割為兩段或更多。
展開 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 CEL(例11) 地震工況下的邊坡大變形模擬 ￥70
ABAQUS CEL(例11) 地震工況下的邊坡大變形模擬
一、建模技術
地震工況下邊坡可能失穩進而出現滑坡現象，為避免模擬滑坡時網格產生的畸變問題，采用耦合歐拉拉格朗日法(CEL)進行滑坡的大變形模擬；土體本構采用摩爾庫倫模型；采用模型底部小范圍內的周期性荷載模擬地震荷載。
二、模型及部分結果展示
圖1：藍色為邊坡；紅色為空氣層
圖2：網格的劃分
圖3：賦予模型初始應力
圖4：土體達到地應力平衡時的應力分布
圖5：土體底部的地震荷載施加區域
圖6：所施加的周期性荷載（地震荷載）
圖7：邊坡因地震荷載產生的位移
圖8：地震波產生的區域
展開 ABAQUS變形放大比例Deformation Scale Factor
The Deformation Scale Factor options are in the lower left corner of the Basic page，等比例的話就在Uniform value 這個地方填寫，如果填1，說明是按實際變形大小來plot
