單元質(zhì)心的案例
Abaqus二次開發(fā)-部件體素化(依據(jù)體素?cái)?shù)目) ￥20
結(jié)果如圖所示：
對比：
插件說明
體素生成技術(shù)說明：
該插件首先根據(jù)部件的尺寸大小，構(gòu)建出可以剛好包裹部件的網(wǎng)格體，再遍歷每一個(gè)單元，根據(jù)單元的質(zhì)心是否包絡(luò)再部件實(shí)體中，來決定單元的取舍。
本次以一個(gè)球體來演示單元質(zhì)心的包絡(luò)情況。如圖所示：
情況一：不包絡(luò)
情況二：包絡(luò)
在情況一中，單元質(zhì)心在球體之外，未包絡(luò)，但是單元卻與部件接觸，這種情況，單元會(huì)判斷為None(也就是舍去),在情況二中，質(zhì)心和單元都包絡(luò)在球體中，這種情況，單元會(huì)判斷為保留。這里就列舉這兩種情況。詳細(xì)關(guān)系見下表：
質(zhì)心，單元與部件的位置關(guān)系對單元取舍的影響
其中，質(zhì)心與單元都與部件進(jìn)行接觸時(shí)（也就是說質(zhì)心剛好在部件表面時(shí)），單元也進(jìn)行保留。
這里的包絡(luò)是指質(zhì)點(diǎn)或單元剛好被部件包裹的情況（也就是說其剛好在部件體內(nèi)，不超出表面）。接觸是指與部件產(chǎn)生交集。
最后，當(dāng)遍歷了每一個(gè)單元后，就可以完成部件體素的轉(zhuǎn)換。
注意：
由于要遍歷每一個(gè)單元，所以當(dāng)單元總數(shù)較多時(shí)，比較耗時(shí)。
XYZ值越大，網(wǎng)格越細(xì)致，單元數(shù)愈多，耗時(shí)越長。
遍歷單元數(shù)可由XYZ的乘積簡單計(jì)算得出。
操作對象：單一part，不適用于殼模型。
體素大小不按尺寸定義。
體素大小：由全局坐標(biāo)系下三個(gè)方向的最大體素?cái)?shù)和包圍part的box的大小決定。
體素宏觀尺寸與原部件一致。
問題排除：
由于引用abaqus內(nèi)核函數(shù)來判斷單元存留，當(dāng)遍歷單元對不保留單元進(jìn)行判斷時(shí)會(huì)發(fā)生警告。已嘗試引用warnings模塊最高等級對警告進(jìn)行抑制，但是，毫無效果。所以該部分已在源代碼中刪除。對該問題進(jìn)行保留，所以，在進(jìn)行比較細(xì)致的體素轉(zhuǎn)化時(shí)，收到警告是正常現(xiàn)象，非bug。
展開 Abaqus二次開發(fā)-部件體素化（依據(jù)體素大小） ￥20
插件說明
體素生成技術(shù)說明：
（與上一篇插件生成技術(shù)大致相同）
該插件首先根據(jù)部件的尺寸大小，構(gòu)建出可以包裹部件的最小網(wǎng)格體，再遍歷每一個(gè)單元，根據(jù)單元的質(zhì)心是否包絡(luò)再部件實(shí)體中，來決定單元的取舍。
本次以一個(gè)球體來演示單元質(zhì)心的包絡(luò)情況。如圖所示：
情況一：不包絡(luò)
情況二：包絡(luò)
在情況一中，單元質(zhì)心在球體之外，未包絡(luò)，但是單元卻與部件接觸，這種情況，單元會(huì)判斷為None(也就是舍去),在情況二中，質(zhì)心和單元都包絡(luò)在球體中，這種情況，單元會(huì)判斷為保留。這里就列舉這兩種情況。詳細(xì)關(guān)系見下表：
質(zhì)心，單元與部件的位置關(guān)系對單元取舍的影響
其中，質(zhì)心與單元都與部件進(jìn)行接觸時(shí)（也就是說質(zhì)心剛好在部件表面時(shí)），單元也進(jìn)行保留。
這里的包絡(luò)是指質(zhì)點(diǎn)或單元剛好被部件包裹的情況（也就是說其剛好在部件體內(nèi)，不超出表面）。接觸是指與部件產(chǎn)生交集。
最后，當(dāng)遍歷了每一個(gè)單元后，就可以完成部件體素的轉(zhuǎn)換。
這里由于是按體素大小進(jìn)行生成，所以不可避免與原部件宏觀尺寸產(chǎn)生偏差，3個(gè)方向的偏差至多為設(shè)定的單個(gè)體素相應(yīng)的長寬高的數(shù)值。體素生成的基準(zhǔn)點(diǎn)為剛好包裹原部件的box的中心。此插件也將偏差進(jìn)行輸出，偏差為宏觀尺寸在全局坐標(biāo)下三個(gè)方向的偏差。
注意：
由于要遍歷每一個(gè)單元，所以當(dāng)單元總數(shù)較多時(shí)，比較耗時(shí)。
體素越小，網(wǎng)格越細(xì)致，單元數(shù)愈多，耗時(shí)越長，但是誤差也越小（這種情況不包括公因子，如果體素塊按公因子設(shè)定，就與原部件宏觀尺寸無偏差）向原部件逼近。
遍歷單元數(shù)可由剛好包裹原部件的box的體積與體素塊的體積之商粗略計(jì)算得出。
操作對象：單一part，不適用于殼模型。
體素大小按尺寸定義。
展開 有限體積方法處理非結(jié)構(gòu)網(wǎng)格-多邊形和多面體的幾何的處理
上一篇文章中已經(jīng)講了通過離散求解的流程，目的就是先構(gòu)造出每個(gè)網(wǎng)格的代數(shù)方程，代數(shù)方程的具體形式如下，碰到的問題就是需要求解三個(gè)系數(shù)項(xiàng)bc，af，ac
但是這三個(gè)系數(shù)項(xiàng)就是采用有限體積方法離散獲取的，每個(gè)都有自己的特點(diǎn)， 有的需要當(dāng)前網(wǎng)格的值或者這個(gè)值的梯度，有的需要與當(dāng)前網(wǎng)格相鄰的單元的值。而且前面也列出了這三個(gè)公式的具體表達(dá)式，如下
前面我們已經(jīng)知道求解這三項(xiàng)只需要知道gDiff和Tf還有▽?三個(gè)就可以了這三個(gè)就可以了，而且這三相都是和網(wǎng)格相關(guān)的，而且都需要長篇大論的寫，所以我們就一個(gè)一個(gè)的處理，首先是這個(gè)gDiff吧，先給出一個(gè)網(wǎng)格的圖片，圖片來自互聯(lián)網(wǎng)哈。
一個(gè)典型的非結(jié)構(gòu)網(wǎng)格
這里把公式先寫出來
Ef是表面向量的類正交擴(kuò)散分量，dCF是兩個(gè)單元的質(zhì)心距離，e是兩個(gè)網(wǎng)格單元質(zhì)心向量的單元向量，Sf是面的法向向量,已經(jīng)知道的是網(wǎng)格中所有的點(diǎn)的坐標(biāo)都。
首先來完成單元的質(zhì)心計(jì)算吧，當(dāng)然質(zhì)心的計(jì)算其實(shí)也是一個(gè)比較漫長的步驟，首先是把每個(gè)面的面心計(jì)算出來，對于三角形來說每個(gè)面的面心直接加權(quán)就可以，但是對于多邊形來說還需要轉(zhuǎn)化為三角形再處理。
1.大致流程是這樣的，通過面上的頂點(diǎn)坐標(biāo)加和然后除以頂點(diǎn)的數(shù)目可以計(jì)算到幾何中心
2.將幾何中心與邊上的兩個(gè)點(diǎn)相連得到n邊形對應(yīng)的n個(gè)三角形
3.將這n個(gè)三角形的面矢量和面積分別計(jì)算出來，還有三角形的幾何中心
4.將這n個(gè)三角形的幾何中心和對應(yīng)的面積加權(quán)就可以得到面的面心，三角形面積矢量加和除以2就是這個(gè)面的面積矢量Sf。
為了省去你的麻煩，我把函數(shù)已經(jīng)定義好了，當(dāng)然你也可以自己編寫，最好以類的形式寫出來，計(jì)算結(jié)果也存為類里面的屬性，可以減少計(jì)算量。
展開 COMSOL 軟件 5.4 版本新功能： 通過材料激活仿真制造工藝
在實(shí)踐中，這意味著如果單個(gè)網(wǎng)格單元擁有多個(gè)積分點(diǎn)，那么它可以實(shí)現(xiàn)部分激活。
對網(wǎng)格單元中的各個(gè)積分點(diǎn)進(jìn)行激活。
案例 2：逐單元激活
現(xiàn)在考慮另一種情況：使整個(gè)網(wǎng)格單元激活，但并非基于單個(gè)積分點(diǎn)。為此，你需要修改激活表達(dá)式，改成在每個(gè)網(wǎng)格單元中對每個(gè)積分點(diǎn)進(jìn)行同等計(jì)算。這一操作可以利用質(zhì)心 運(yùn)算符來完成。我們修改了上一個(gè)案例中的激活表達(dá)式，如下圖所示。現(xiàn)在，我們在網(wǎng)格單元質(zhì)心上計(jì)算 Y 坐標(biāo)，這意味著對于網(wǎng)格單元中的所有積分點(diǎn)，激活表達(dá)式將獲得相同的值。
逐單元材料激活的激活表達(dá)式。
在下圖中的網(wǎng)格單元內(nèi)部，我們對單元質(zhì)心運(yùn)算了激活表達(dá)式，所以四個(gè)積分點(diǎn)均處于被激活狀態(tài)。
通過使用質(zhì)心運(yùn)算符對網(wǎng)格單元中的所有積分點(diǎn)進(jìn)行激活。
案例 3：使用之前的激活狀態(tài)
假設(shè)你需要模擬激光熔覆工藝，在工藝過程中，填充材料隨著時(shí)間的推移熔化并沉積。激光束位置隨時(shí)決定材料的沉積位置。工藝過程開始后，激光束的整個(gè)軌跡限定了之前被激活材料的區(qū)域。（關(guān)于模擬激光束運(yùn)動(dòng)的詳細(xì)信息，請閱讀介紹移動(dòng)載荷與約束建模的博客文章。）借助變量 wasactive，我們不必以通過數(shù)學(xué)方式描述激光束軌跡。這類情況的激活表達(dá)式可以表述為：
(描述激光束當(dāng)前位置的邏輯表達(dá)式）|| solid.wasactive
該表達(dá)式規(guī)定，如果“描述激光束當(dāng)前位置的邏輯表達(dá)式”為真，或者 如果材料在仿真的任何之前時(shí)間（或參數(shù)步驟）曾處于被激活狀態(tài)，則材料處于激活狀態(tài)。如果激活表達(dá)式?jīng)]有使用 wasactive 變量，那么激光束經(jīng)過后，材料將變成失活狀態(tài)，這很可能與預(yù)期效果相反。
結(jié)果可視化
假設(shè)你模擬過了材料隨時(shí)間沉積的時(shí)間依賴性工藝，那么僅顯示域中被激活部分的結(jié)果這一功能或許會(huì)吸引你的興趣。如下圖所示，在過濾器 節(jié)點(diǎn)下的包含邏輯表達(dá)式 中添加變量isactive，即可實(shí)現(xiàn)這項(xiàng)操作。
展開 
abaqus中計(jì)材料參數(shù)不為常數(shù)的兩種處理方法-Part2
UMAT子程序具有強(qiáng)大的功能，使用UMAT子程序：
(1)可以定義材料的本構(gòu)關(guān)系，對ABAQUS 材料庫進(jìn)行擴(kuò)充；
(2)幾乎可以將任何本構(gòu)關(guān)系運(yùn)用到ABAQUS中的任何單元中；
但是，用戶必須在UMAT中提供材料本構(gòu)模型的雅可比（Jacobian）矩陣，即應(yīng)力增量對應(yīng)變增量的變化率，而實(shí)際上，這也是Umat子程序最主要的任務(wù)。
承前Part1所述，由于假定彈性模量與坐標(biāo)相關(guān)，使用Umat子程序的首要任務(wù)就是找到彈性模量E和坐標(biāo)之間的關(guān)系；之后，才是根據(jù)坐標(biāo)得到彈性模量和泊松比之后對DDEDDS的編寫。這里，依據(jù)Abaqus的幫助文檔，給出了Umat子程序中一些變量的定義：
圖1 Abaqus提供的Umat函數(shù)接口
圖2 部分Umat變量定義
詳細(xì)的參數(shù)介紹大家可以查看abaqus幫助。不難發(fā)現(xiàn)其中并不包含有單元或者節(jié)點(diǎn)的坐標(biāo)信息，但是卻包含有單元的單元編號信息(NOEL)。而通過單元編號又可以關(guān)聯(lián)起來單元的坐標(biāo)。
圖3 線彈性材料的DDEDDS
圖3中EMOD=PROPS(1)*A(NOEL)/0.125 就是根據(jù)單元NOEL的坐標(biāo)確定其彈性模量的語句，不難看出這個(gè)彈性模量是INP文件中給定的值得A(NOEL)/0.125倍，其中A數(shù)組為自定義的數(shù)組，包含有單元的坐標(biāo)值數(shù)組。
5 算例
對于如圖5所示的4階魔方結(jié)構(gòu)，x軸左端面約束x方向自由度，x軸右端面加載位移。材料彈性模量由約束面到加載面分別為10、30、50、70GPa，而其對應(yīng)的單元質(zhì)心x軸坐標(biāo)分別為0.125、0.375、0.625、0.875。圖5中所示云圖為單元應(yīng)變E11。
展開 Abaqus利用region讀取結(jié)果數(shù)據(jù)
單元節(jié)點(diǎn)（未平均）
CENTROID ！單元質(zhì)心點(diǎn)
Abaqus利用region讀取結(jié)果數(shù)據(jù)
單元節(jié)點(diǎn)（未平均）
CENTROID ！單元質(zhì)心點(diǎn)
轉(zhuǎn)自公眾號——ABAQUS大世界
旨在分享，若侵即刪.
基于深度學(xué)習(xí)的電阻抗、電磁與電容層析成像方法研究
所有目標(biāo)區(qū)域通過判斷 FEM 網(wǎng)格單元質(zhì)心是否落入給定幾何范圍內(nèi)構(gòu)造掩碼，并統(tǒng)一賦予電導(dǎo)率擾動(dòng)值（如 +1）。
2.3 正向電壓模擬與樣本存儲
對每一組 img_target 與其對應(yīng)的 img_bg 執(zhí)行前向求解，得到擾動(dòng)前后邊界電壓分布 v_bg 與 v_target。最終每個(gè)樣本包含：
img_target：擾動(dòng)電導(dǎo)率圖像（有限元模型格式）；
v_bg：背景電壓（無目標(biāo)）；
v_target：目標(biāo)擾動(dòng)后的電壓；
可視化圖像：模型可視化圖 .jpg；
存儲格式：所有樣本統(tǒng)一保存為 dataset_XXX.mat 文件，保存在 dataset 文件夾中。
數(shù)據(jù)集共構(gòu)建上萬組樣本，確保具備豐富的幾何結(jié)構(gòu)、多目標(biāo)組合以及隨機(jī)擾動(dòng)條件，適用于神經(jīng)網(wǎng)絡(luò)訓(xùn)練與泛化能力評估。
圖2 數(shù)據(jù)集構(gòu)建過程
神經(jīng)網(wǎng)絡(luò)（圖像回歸分析）
3.1 網(wǎng)絡(luò)輸入構(gòu)建：電壓圖像化
原始電壓數(shù)據(jù)為來自前向正問題仿真的邊界測量（208 通道），通過下述步驟轉(zhuǎn)換為網(wǎng)絡(luò)可處理的二維圖像：
歸一化差分電壓：
圖像重構(gòu)格式：將 208 通道映射為13×16的二維電壓圖像。
尺寸標(biāo)準(zhǔn)化：所有圖像 resize 為 224×224，并復(fù)制為 3 通道 RGB 格式，以適配 CNN 輸入。
展開 【ANSYS 17】更快、更好的計(jì)算流體動(dòng)
– 尤其有助于具有本地多面體網(wǎng)格和/或高度拉伸單元的案例。
? 代數(shù)多網(wǎng)格求解器現(xiàn)在自動(dòng)對線性系統(tǒng)進(jìn)行重新排序。
? 在使用基于壓力的求解器時(shí)可使用彎曲面梯度校正選項(xiàng)。這對具有十分不平整面或單元質(zhì)心位于單元外的六面體或多面體單元非常有用。
高性能計(jì)算（HPC）提供經(jīng)過驗(yàn)證的可擴(kuò)展能力，最大可達(dá)12.9萬個(gè)內(nèi)核
ANSYS Fluent不斷打破之前的仿真世界記錄，現(xiàn)已擴(kuò)展到129,000個(gè)計(jì)算內(nèi)核，且運(yùn)行效率高達(dá)90%。因此，即使是最大的仿真模型也可提供快速及時(shí)的結(jié)果。由于對分區(qū)、通信和AMG求解器的進(jìn)一步優(yōu)化，F(xiàn)luent求解器已經(jīng)展現(xiàn)出前所未有的擴(kuò)展性能與效率，可充分滿足并行處理和模型復(fù)雜性的要求。此外，針對滑動(dòng)界面交叉和負(fù)載平衡的優(yōu)化算法，可為MDM和滑動(dòng)網(wǎng)格案例帶來改善的擴(kuò)展能力。比如剎車片、發(fā)動(dòng)機(jī)曲軸箱潤滑模型和內(nèi)燃機(jī)（ICE）案例等。例如，160萬單元的ICE案例在384個(gè)內(nèi)核上并行效率提高到55%。最后同樣重要的是，案例/數(shù)據(jù)IO時(shí)間通過HDF5格式得到顯著縮短。
在劣質(zhì)網(wǎng)格上，彎曲面梯度校正的結(jié)果與優(yōu)質(zhì)網(wǎng)格的結(jié)果更加貼近
這些HPC改進(jìn)并非只針對極端案例。通過優(yōu)化通信和去除妨礙案例讀取和并行構(gòu)建的瓶頸，較小的案例和HPC集群（尤其是2000個(gè)核以上的集群）將大獲裨益。包括分區(qū)和鄰域創(chuàng)建在內(nèi)的HPC設(shè)置可加快30倍。由于所有讀/寫階段的輸入/輸出時(shí)間縮短，HDFIO得以增強(qiáng)。寫入時(shí)間也得到最顯著的改善。
滑移與變形網(wǎng)格
Fluent中的一些滑移與變形網(wǎng)格案例獲得了顯著提速，包括有額外物理場出現(xiàn)的情況。具體改進(jìn)包括鄰域優(yōu)化、滑動(dòng)界面優(yōu)化和并行求解器優(yōu)化等。例如，曲軸箱潤滑仿真加快85%，氣缸內(nèi)燃燒仿真加快55%。
加快收斂
針對耦合求解器，F(xiàn)luent現(xiàn)在默認(rèn)使用保守的網(wǎng)格粗化。
展開 NO.20 聚能射流成型（SPH）
該方法是將有限元網(wǎng)格單元轉(zhuǎn)換為 SPH 粒子，新生成的 SPH 粒子坐標(biāo)位于原有網(wǎng)格單元的質(zhì)心，粒子的質(zhì)量與網(wǎng)格單元的質(zhì)量相同，其占據(jù)的體積為原有網(wǎng)格單元的體積。（2）先分別建立炸藥和藥型罩的空心外殼，然后在其中填充等間距的SPH粒子。本文采用第二種方法構(gòu)建有限元模型。
有限元模型
均勻排布的炸藥和藥型罩粒子
開展三種工況模擬，通過修改相關(guān)參數(shù)，不斷優(yōu)化SPH聚能射流的成型形態(tài)。
（1）FORM=1
射流出現(xiàn)斷裂、粒子飛散
（2）FORM=0
射流成型較好，沒有粒子飛散，但頭部粒子較分散
射流速度分布
（3）FORM=0，考慮人工粘性
射流 成型良好 ，頭部粒子未分散
射流速度分布
射流成型過程
未經(jīng)許可，不得私自轉(zhuǎn)發(fā)
展開 準(zhǔn)確性、收斂性和網(wǎng)格質(zhì)量
他提出了三種類型的有效性檢查：
類型 1 檢查單元格是否具有相互不相交的正體積和面。這里又是“這個(gè)網(wǎng)格壞了嗎？”的一個(gè)例子。問題。
類型 2 檢查內(nèi)部單元面是否與另一個(gè)內(nèi)部面唯一匹配，以及邊界單元面是否位于被劃分網(wǎng)格的對象的幾何模型上。
類型 3 檢查幾何模型的每個(gè)表面是否完全被邊界單元面覆蓋，幾何的每個(gè)硬邊緣是否被邊界單元面的邊緣覆蓋，以及邊界面面積的總和是否與實(shí)際幾何表面積匹配。
圖 1：一個(gè)簡單的演示，說明從單元幾何角度（右）來看，較差的網(wǎng)格如何比具有“完美”單元（左）的網(wǎng)格產(chǎn)生更低的離散化誤差。來自參考文獻(xiàn) 1c。
弗吉尼亞理工大學(xué)的 Christopher Roy 教授展示了一個(gè)反直覺的例子（至少從先驗(yàn)指標(biāo)的角度來看），即 2D Burger 方程在自適應(yīng)網(wǎng)格（具有廣泛變化的傾斜、縱橫比和其他指標(biāo)的單元格）上的解與完美正方形網(wǎng)格上的解相比，離散化誤差要小得多。僅從這個(gè)例子就可以清楚地看出，僅基于單元格幾何形狀的指標(biāo)并不是網(wǎng)格質(zhì)量的良好指標(biāo)，因?yàn)樗c解決方案的準(zhǔn)確性有關(guān)。
求解器的角度
研討會(huì)很幸運(yùn)有幾位流動(dòng)求解器開發(fā)人員參與，他們分享了有關(guān)求解器如何受網(wǎng)格質(zhì)量影響的詳細(xì)信息。所有這些的共同點(diǎn)是收斂性和穩(wěn)定性受網(wǎng)格質(zhì)量的影響比求解精度更直接。
差價(jià)合約++
Metacomp Technologies 的 Vinit Gupta 將單元偏度和單元大小變化列為結(jié)構(gòu)化網(wǎng)格需要注意的兩個(gè)質(zhì)量問題。特別是，在梯度較低的遠(yuǎn)場中跨塊邊界的網(wǎng)格細(xì)化對收斂具有強(qiáng)烈的負(fù)面影響。對于非結(jié)構(gòu)化和混合網(wǎng)格，邊界層中的各向異性四面體以及從棱柱到邊界層外四面體的過渡也可能存在問題。
Gupta 還指出了與度量計(jì)算相關(guān)的兩個(gè)問題。依賴于將細(xì)胞分解成四組的細(xì)胞體積計(jì)算不是唯一的，并且取決于分解的方式。
展開 
ABAQUS中材料參數(shù)不為常數(shù)的兩種處理辦法——Part1
大家知道，abaqus材料中材料是在*section對截面屬性定義的時(shí)候?qū)?em>單元組進(jìn)行賦值的。并且GUI中也沒有提供變參數(shù)的接口(這里順路提一下，載荷是可以通過函數(shù)給出變載荷的，GUI中有接口)。這就需要我們采取一些特殊的措施。本教程給出2種行之有效的方法，它們分別是：(1)編輯inp文件和(2)使用Umat。
2方法
方法
要求
難度
編輯inp
會(huì)一門編程語言，能夠完成對文本文件的讀取和寫入
4星
使用Umat
會(huì)使用abaqus的Umat用戶子程序，計(jì)算文件中不需要使用其它的用戶子程序
3星(如果對abaqus用戶子程序不熟悉，則為5星)
上表中給出兩種方法的一些基本要求，編輯inp文件的方法要求更低，它通過直接編輯inp文件的方法來解決問題。而使用Umat的方法要求更高，需要用戶對子程序有一定的了解，并且也需要會(huì)一門編程語言。
同樣是根據(jù)坐標(biāo)確定材料參數(shù)，使用Umat時(shí)，可以直接通過abaqus用戶子程序自帶的子函數(shù)直接獲得單元的節(jié)點(diǎn)坐標(biāo)，進(jìn)而通過坐標(biāo)確定材料參數(shù)；而使用編輯inp的方法則需要首先通過編寫程序?qū)δ阋呀?jīng)生成的inp文件(包含節(jié)點(diǎn)、單元信息)進(jìn)行讀取，之后對每一個(gè)單元進(jìn)行材料賦值。從這方面講，對熟悉Umat的用戶而言，使用Umat解決這個(gè)問題更為方便。
3編輯inp方法的簡要實(shí)施過程
以Fortran+VS為例，在編輯inp文件時(shí)，涉及到了單元節(jié)點(diǎn)信息的讀取、計(jì)算單元(質(zhì)心)坐標(biāo)、計(jì)算彈性模量E、和輸出這幾個(gè)子程序。如圖1所示，prereadnew為讀取節(jié)點(diǎn)、單元信息子程序；processmat為通過單元坐標(biāo)計(jì)算彈性模量子程序；writeabaqus為輸出子程序。
展開 設(shè)計(jì)仿真 | 新版本Dytran 用戶子程序調(diào)用及案例分享
案例說明
在液體晃蕩分析過程中，我們通常需要將液體部分的質(zhì)心進(jìn)行輸出，以查看其在晃蕩過程中液體對結(jié)構(gòu)狀態(tài)的影響，要實(shí)現(xiàn)該功能需要通過其EXOUT子程序接口將所有Euler單元的質(zhì)量及位置坐標(biāo)進(jìn)行處理得到整體質(zhì)心位置并輸出。
將Dytran安裝目錄下/usr-subrtns/文件夾內(nèi)的dytran_usersub.f子程序源文件作為模板，編輯該文件下的EXOUT子程序，遍歷每個(gè)歐拉單元的質(zhì)量以及質(zhì)心，通過數(shù)據(jù)處理后，實(shí)現(xiàn)歐拉單元內(nèi)液體部分質(zhì)心輸出的功能。
在Dytran Explorer中，在提交作業(yè)時(shí)，將UDS路徑指向編輯好的Fortran源代碼，點(diǎn)擊RUN按鈕后即可進(jìn)行作業(yè)提交。
程序會(huì)首先將源代碼編譯為動(dòng)態(tài)鏈接庫，并在目錄下生成dytran_usersub.dll文件，該文件可以代替子程序源文件，在沒有編譯器的環(huán)境下進(jìn)行子程序的提交。當(dāng)采用動(dòng)態(tài)鏈接庫進(jìn)行作業(yè)提交時(shí)，只能采用CMD形式，具體命令如下：
文章小結(jié)
Dytran UDS接口讓用戶方便的實(shí)現(xiàn)求解器的二次開發(fā)，大大擴(kuò)充了Dytran的功能，有效提高了Dytran在工程上的通用性。
展開 2019.07.29 關(guān)于VUMAT的一些參數(shù)
分析步開始的時(shí)間為totalTime-stepTime
dt：時(shí)間增量的大小
cmname：用戶指定的材料名（左對齊）避免與內(nèi)部材料名沖突，不要用“ABQ_”命名
coordMp（nblock，*）：材料點(diǎn)坐標(biāo) 殼單元中性面的材料點(diǎn)，梁和管單元的質(zhì)心
charLength（nblock）：特征元素長度，它要么是基于幾何平均值的默認(rèn)值，要么是用戶子程序VUCHARLENGTH中定義的用戶自定義特征元素長度。一階元素的默認(rèn)值是跨元素的一根特征線的長度（一階四邊形 面積開二次方，一階六面體 體積開三次方）;對于二階元素，它是一般長度的一半。對于梁、管道和桁架，默認(rèn)值是沿著元素軸的特征長度。對于膜和殼，它是參考表面的一個(gè)特征長度。對于軸對稱單元，它是r-z平面上的一個(gè)特征長度。對于chohesive單元等于其本構(gòu)的厚度。
props（nprops）：用戶提供的材料特性
density（nblock）：Current density at the material points in the midstep configuration.This value may be inaccurate in problems where the volumetric strain increment is very small. If an accurate value of the density is required in such cases, the analysis should be run in double precision. This value of the density is not affected by mass scaling.
展開 Abaqus2023新功能-1
隨機(jī)響應(yīng)分析性能和功能增強(qiáng)
產(chǎn)品：Abaqus/Standard
在隨機(jī)響應(yīng)分析程序中，Abaqus/Standard節(jié)點(diǎn)和單元輸出的性能得到了顯著提高。此功能在Abaqus 2022 FD04 （FP.2232）版本中首次提供。
將隨機(jī)響應(yīng)分析的節(jié)點(diǎn)輸出階段加快2–500倍
將典型隨機(jī)響應(yīng)分析的求解時(shí)間縮短2–20倍
允許在隨機(jī)響應(yīng)過程中使用更多特征模態(tài)，并獲得大量節(jié)點(diǎn)的輸出
支持歷史記錄、字段和打印輸出
在隨機(jī)響應(yīng)分析程序中，在單元輸出和計(jì)算輸出變量MISES 和RMISES 方面，性能得到了顯著提高。此功能在 Abaqus 2022 FD02 （FP.2214）版本中首次提供。可以請求計(jì)算積分節(jié)點(diǎn)、單元質(zhì)心、單元節(jié)點(diǎn)處的馮米塞斯應(yīng)力和RMS馮米塞斯應(yīng)力，以及節(jié)點(diǎn)處的平均值；馮米塞斯和RMS馮米塞斯應(yīng)力的計(jì)算值保存在ODB和SIM數(shù)據(jù)庫中。
2. 固有頻率提取中的殘差模態(tài)增強(qiáng)
產(chǎn)品：Abaqus/Standard
固有頻率提取過程中，提供了多項(xiàng)殘差模態(tài)增強(qiáng)功能。此功能在Abaqus 2022 FD04 （FP.2232）版本中首次提供。
可以使用Abaqus/Standard中其他線性動(dòng)態(tài)過程
支持的任何荷載類型（分布式荷載、子結(jié)構(gòu)荷載等），在AMS固有
頻率提取過程中生成殘差模態(tài)。
可以在固有頻率提取分析中定義實(shí)載荷和虛載荷，以生成殘差模態(tài)。
可以在純聲學(xué)模型的固有頻率提取分析中定義殘差模態(tài)。
3.
展開 