abaqus曲面映射的案例
關于DynaForm回彈補償曲面映射失敗的問題
最近做幾個例子，發現很多版本的dynaform不能正常映射曲面，這個是軟件本身的問題！！
比如DynaForm5.9.1系列就不行
目前比較正常的是DynaForm5.9.2.1
所以次功能失敗的可以考慮換一個版本的dynaform試試。或者用第三方的
Abaqus隨機材料映射網格插件：Random Material Mesh - AbyssFish ￥128
插件介紹
Random Material Mesh - AbyssFish 插件可在Abaqus軟件將材料隨機批量賦值給部件的網格單元。插件支持二維及三維部件的所有網格類型，可指定任意多種（實際材料種類＜10^7）不同材料。注意，插件僅提供空材料指定網格功能，并不能生成隨機的材料屬性，需要手動或借助其他插件設置材料行為參數。
模型展示
使用教程
模型以二維部件的軸壓試件為例，介紹隨機材料網格插件的使用方法，并與均質材料進行對比分析。
?打開Abaqus新建幾何模型，并劃分網格。
?在Abaqus的Plug-ins菜單下，找到AF_ RandomMaterialMesh，點擊打開插件。選擇需要編輯的模型及部件，并制定其材料種類的數量。這里制定10種不同的材料。
?切換到屬性模塊，可查看材料制定情況。
?打開材料管理器，依次對制定的材料進行參數設置。這里也可以采用復制的方式，將已有的材料屬性復制一份并稍作修改，注意復制后需要保證材料名稱與需要編輯的材料一致。本案例中材料（1~10）的彈性模量分別為（1~10）× 10^5。
?建立分析步，指定荷載等并提交分析。
?最終非均質材料與均質材料（E = 5.5×10^5）計算結果如下。
說明提醒
插件可運行在WindowsXP、7、8、10、11系統上，支持Abaqus6.14、Abaqus2017~2023及以上版本。
插件需要注冊，售價為單機許可的價格，購買后請聯系QQ：1135122921獲取許可證。
展開 【Abaqus 3D打印建模】之 極小曲面 I --Matlab生成極小曲面 ￥79
image_process=/format,webp/resize,w_760" data-initial-src="https://img.jishulink.com/upload/202105/256c0c045a704f8486cdf6812b97ba0d.png">
</div><h2 class="ql-align-center"><strong><em>斷面封閉的曲面</em></strong></h2><p><br></p><p><img src="https://img.jishulink.com/upload/202105/4404ab4e814e412fb0f4b8c893975b35.png" alt="minim3.png"></p><p>****上述封閉曲面為空心曲面，若要用實體，需用CAD軟件進行實體轉換 或 <strong>用abaqus的mesh模塊中tri to tet 填充為實體</strong></p><p><br></p><p>雖然上述stl文件導入后在abaqus中可以正常計算，但網格劃分有些不好。怎么有畫網格？大家可以提提高見。據我所知，由于此類曲線很<strong>難轉化為 幾何模型</strong>，只能通過網格優化軟件或 一些算法實現網格重新劃分。
展開 【Abaqus 3D打印建模】之 極小曲面 II --python生成極小曲面 ￥79
***請注意，附件僅為現式和隱士極小曲面生成和輸出為stl的python代碼，并不包括加厚***
***有加厚的需求，請看其他算例。加厚后輸出的stl，也是空心封閉的殼，不是實心的實體***
0.算例
上一個帖子介紹了怎么用matlab建立極小曲面，詳情見Matlab創建極小曲面。
下面是個簡單的算例，在y方向壓縮極小曲面之Gyroid，幾何模型建立方法見下文，建立后陣列并有畫網格導入abaqus即可。
為了對比該極小曲面的應力水平，采用同樣的材料做了單軸壓縮，兩種情況對比如圖所示：
從圖中可以看到，如果僅去極小曲面上的一個點作為其應力應變，其曲線甚至比實心立方體還高，顯然這是不合理的。出現這種現象的主要原因是，此類細觀結構或變形不均勻時，不能取一個點代表整個模型，不然很可能會出現較大的誤差。對此問題，細觀力學有些方法，比如作用力反作用力法，體積平均法等，但也有人認為對于細觀結構，作用力反作用力法 不太合理。體積平均法的簡單表達式如下：
即模型中每一個單元的應力（應變）對單元體積積分后，除以模型整個體積。上述應力應變曲線也證實，采用該方法能夠得到較為真實的數據。
那么，怎么通過體積平均法獲得數據呢？ 在計算結束后，需要通過python腳本對數據進行處理，輸出中需要EVOL（單元體積）以便獲得各個單元的體積。
1.介紹
之前已經介紹過什么是極小曲面，同時根據表達式定義為 隱式 和顯式 極小曲面，主要區別是 隱式極小曲面一般只有一個方程，不容易將x,y,z獨立表達出來，如下圖所示：
今天介紹怎么用python生成上述極小曲面并輸出為stl文件。
展開 
CFD結果映射到ABAQUS模型上
*Part, name=PART-1
*Node
1, 54.4015923, -5.89080048, 30.
2, 53.8370857, -11.9988546, 30.
3, -48.6207314, -3.89129758, 84.
4, -49.0763855, -4.09715223, 94.4679184
5, -53.4180222, -3.19716525, 94.5191727
6, -52.7271576, -2.91121149, 84.
7, -44.4379768, -4.5009594, 84.
8, -44.6520882, -4.47322035, 94.4163437
步驟三，參數映射：
ABAQUS中自帶了場定義模塊，通過粘貼坐標和對應變量值，自動實現插值映射。
如果有一些特殊情況，比如我的ABAQUS模型是帶厚度的實體單元，我想讓那個面的溫度映射到整個厚度方向單元上。
這就需要我們自己做插值算法。可以采用如下思路：
(1) 對于每個CAE模型的節點，計算它與所有CFD節點之間的距離。
(2) 查找到每個CAE節點最鄰近的CFD節點，將這個最鄰近CFD數據點的溫度插值到CAE中，這樣便可得到CAE溫度的插值數據。
(3) 最后使用插值結果，對inp進行自動修改，將溫度結果施加到每個節點上。
步驟四，渲染對照：
再開發一個可視化模塊，將原始CFD溫度分布和映射后面的CAE模型變量分布進行渲染。幫助我們確認是否有明顯的映射錯誤。
展開 【網格優化功能：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二次開發】骨料隨機填充及映射網格法
<h2 class="ql-align-center"><strong>0.出發點</strong></h2><p>用abaqus做細觀模擬，隨機填充顆粒/纖維是跳不過的坎，網上有很多插件或腳本，技術鄰也有不少大咖做了些插件，但個人而言更傾向于腳本，雖然插件使用起來方便，不過有時需要根據具體要求修改模型，此時腳本的靈活性就體現出來了。今天跟大家分享怎么通過python編寫隨機填充骨料的思路：</p><ol><li>二維模型，主要是將基體做partition，比較簡單；而做partition用的草圖就需要用python批量生成，操作跟普通CAE界面做partition類似。</li><li>三維模型直接用上述方法比較難，此時可以建立多個part，或將骨料作為一個part跟基體在assembly模塊進行裝配，并巧用Merge/cut功能即可。</li><li>網格映射法主要是通過將獲得骨料中的點坐標，把它映射到已畫好結構化網格的另一個part中；這種方法的有點在于施加<strong>周期性邊界條件</strong>很方便。</li></ol><h2 class="ql-align-center"><strong>1. 普通隨機骨料填充</strong></h2><p>下面是按照上述思路編寫的2D及3D骨料隨機填充模型，其中二維模型支持<strong>圓形、橢圓及多邊形</strong>（可定義邊數），指定粒徑和含量比較方便。
展開 ABAQUS二維混凝土細觀模型的數字化重建技術(二)圖像映射
ABAQUS二維混凝土細觀模型的數字化重建技術(一)幾何重構
https://www.yqgqt.org.cn/post/1990726
本篇介紹二維混凝土細觀模型在ABAQUS中數字化重建技術的第二種方法——基于ABAQUS背景網格的圖像映射方法。混凝土圖像前處理部分與第一種方案一致，這里不多做贅述，將處理完成的混凝土圖像通過ABAQUS Image To Part 2D插件進行導入，導入時縮放比例參數（Scaling）設置為0.3，可將分辨率為500×500 px的圖像建立尺寸為150×150 mm的試件模型（0.3 = 150/500）。模型建立后采用EasyCDP Mortar&ITZ插件設置混凝土損傷塑性材料參數，本案例不考慮骨料的損傷破壞。
設置分析步、載荷后建立作業，并在提交作業前采用ABAQUS CDED插件設置混凝土開裂。
提交作業完成模擬分析。
展開 原創#在abaqus中建立實際形貌的二維或準三維材料映射網格模型
</p><p>如果這方面需求</p><p>ABAQUS斷裂模擬收徒 ，保證快速學會各種ABAQUS斷裂模擬方法 1200/人（將享有各種插件以及程序，價值3000+、專門定制視頻、全程親自教學、各種模型調試及解答問題等等，傾囊相教）</p>
【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曲面方程參數化建模 ￥79.9
abaqus曲面方程參數化建模，需要用到pyhon腳本參數化建模，可以在曲面上拉伸厚度。直接輸入x,y的范圍，厚度建模。以下例子曲面方程為
該腳本可以輸入方程，給定區間建模。
在ABAQUS中基于圓柱坐標系設置關于坐標函數的表面力（keyword 曲面加載,圓柱坐標,面力）
例如下圖所示，受Y方向某拉力作用，各點應力狀態為：
在圓孔中心位置建立圓柱坐標系，該應力狀態在圓柱坐標系下的公式為：
在這種情況下反推物理量，需要對曲面施加基于圓柱坐標系的面力。
案例如下：在圓弧面基于圓柱坐標系施加等效于單向應力狀態的面力。
加載前先建立圓柱坐標系（注意R軸方向為0度位置，T軸方向為角度增大方向，示意圖見文后的加載圖）
具體設置方法為：Load>Create Load>Mechanical>surface traction
選中中間曲面后，先設置徑向力，按以下參數設置：
Distribution：應力分配，點擊后面的f(x)創建一個基于圓柱坐標系的表達式，Local system 要選擇圓柱坐標系，Th為角度變量。
Traction：選擇General，為一般力。
Vector：點擊選擇圖標后，依次選擇（0，0，0） （-1，0，0） ，坐標選擇建立的圓柱坐標系。
注：面力方向矢量是基于所選坐標系，（-1，0，0）就是沿圓柱坐標系下的R軸反向。
Magnitude：選擇應力大小為1。
然后在創建一個Load,設置切向力，如下圖所示，也是基于圓柱坐標系。
再創建一個Load，在整體坐標系下對兩側的平面施加Y方向的面力，大小為1，同時對后面的面施加全約束。
最后加載形式為下圖所示：
求解結果如下圖：
大部分位置應力在0.99~1.01之間，為單向應力狀態，加載方式正確。
本問題的關鍵是面力的方向問題，在選擇面力的方向矢量時，是基于所選坐標系。對于圓柱坐標系，切向力矢量為（0，-1，0）時，即力的方向只沿著theta的反方向。
展開 