abaqus合并曲面的案例
【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文件。
展開 【網格優化功能：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曲面方程參數化建模 ￥79.9
abaqus曲面方程參數化建模，需要用到pyhon腳本參數化建模，可以在曲面上拉伸厚度。直接輸入x,y的范圍，厚度建模。以下例子曲面方程為
該腳本可以輸入方程，給定區間建模。
【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中基于圓柱坐標系設置關于坐標函數的表面力（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的反方向。
展開 