ansys輸出變量的云圖的案例
Abaqus_Python 后處理結點數據提取、場變量云圖輸出 ￥1.99
今天體會了Python腳本操作ODB文件進行數據的輸出、處理等方面的快速高效。
首先介紹一下背景，同事的項目，計算case較多。每個case需要輸出幾個路徑上的部分節點值，而且需要針對不同的case、不同的幀、不同的視角截圖、保存云圖，如果人工出圖和提取節點數據，工作量會非常大，費時會遠遠超出計算所用時間；并且，同事還要求我根據case名稱、組件名稱（模型有多個組件組成，且網格是在part上畫的，所以同一個節點編號可能對應多個節點）等保存成不同的文件。
結點溫度數據沒有保存，因此沒有輸出。下圖是自動輸出的png格式云圖。
比較倉促，代碼也比價粗糙?；仡^做一個詳細的介紹。
主要代碼見下。
展開 ANSYS中的自動化參數研究，自動建模/分網/多參數求解/自動輸出云圖/自動輸出所需結果
模型求解
/SOLU
FLST,2,1,5,ORDE,1
FITEM,2,5
DA,P51X,ALL,
FLST,2,1,5,ORDE,1
FITEM,2,6
SFA,P51X,1,PRES,-1000
EQSLV,PCG,1E-6
solve
后處理自動輸出應力云圖，自動保存所需數據
模型求解后，通過/ANUM、/TSPEC、/TLAB命令定義輸出圖上的注釋，如下圖所示。
/post1/
ANUM ,0,1,-0.59026,-0.7 !注釋位置，注釋字體設置
/TSPEC, 15, 1.200, 1, 0, 0
/TLABEL,-0.947,-0.7,Cylinder Radius =%holrad(count,1)%
將最大主應力云圖輸出在屏幕上，并自動保存為.JPEG格式圖片。
plnsol,s,1
/show,jpeg
plnsol,s,1
/show,close
/show,term
最大主應力通過*get命令提取并賦值給maxs1,儲存在holrad數組的第二列中。*enddo是匹配初始的*do命令。
NSORT,S,1
*GET,maxs1,SORT, ,MAX
holrad(count,2)=maxs1
*enddo
最后，使用*cfopen命令創建holard.txt文件，使用*vwrite寫入數據，holard.txt第一行寫入'Radius'和'Stress'標簽，這里4X指的在'Radius'和'Stress'之間有4個空格。然后將holrad(1,1),holrad(1,2)寫入，holrad(1,1)對應的孔直徑參數，holrad(1,2)對應存儲的最大主應力值。
展開 Abaqus后處理-云圖變量含義（部分）
結果中分量說明：
S11、S22、S33指各軸正應力；
S12指作用于XZ平面(與“2”，即Y軸垂直的平面)內，沿1方向剪應力；
S13指作用于XY平面內，沿1方向剪應力；
S23指作用于XY平面內，沿2方向剪應力。
若為柱坐標，S12、S13、S23分別指：由徑向向環向的剪應力、由徑向向軸向的剪應力、由環向向軸向的剪應力。
Abaqus后處理插件——云圖場變量查詢器 ￥50
基于Python對Abaqus進行后處理的二次開發，形成了場變量查詢插件，寫論文寫報告的實用小工具，話不多說直接上圖。
其實Abaqus后處理模塊本身具有場變量查詢功能（Tools-Query-Probe values），可以滿足使用需求，但其顯示效果非常不友好，以節點（單元）編號+場變量值的形式顯示，理工科氣息十足。不得不說，Abaqus在操作體驗方面而言和Workbench還是有一些差距。
PS：常規的二次開發通過fieldvariable.getsubset方式查詢到的場變量為節點直接插值的計算結果，即節點上的值是從相鄰單元插值后未經過平均化的值，是離散的，而結果顯示云圖默認是按75%平均化處理后的結果，更接近實際。
本插件查詢到的場變量為平均化之后的結果，和云圖的顯示的值一致。
展開 
ABAQUS能量平衡輸出變量
Total energy output quantities
ALLAE
“Artificial” strain energy associated with constraints used to remove singular modes (such as hourglass control), and with constraints used to make the drill rotation follow the in-plane rotation of the shell elements.
ALLCD
Energy dissipated by creep, swelling, and viscoelasticity.
ALLEE
Electrostatic energy.
ALLFD
Total energy dissipated through frictional effects. (Available only for the whole model.)
ALLIE
Total strain energy. (ALLIE = ALLSE + ALLPD + ALLCD + ALLAE + ALLQB + ALLEE + ALLDMD.)
ALLJD
Electrical energy dissipated due to flow of electrical current.
ALLKE
Kinetic energy.
ALLKL
Loss of kinetic energy at impact. (Available only for the whole model.)
ALLPD
Energy dissipated by rate-independent and rate-dependent
展開 ABAQUS中隱式和顯式的節點和單元的輸出變量解析
The output variables listed below are available in Abaqus/Explicit.
Mechanical analysis–nodal quantities
CFORCE
Field: yes History: no .fil: no
Contact normal force (CNORMF) and frictional shear force (CSHEARF).
CDISP
Field: yes History: no .fil: no
Contact opening (COPEN) and accumulated tangential motions (CSLIP1, CSLIP2, and
CSLIPEQ) for general contact analyses.
CEDGEACTIVE
Field: yes History: no .fil: no
Status of contact edges for general contact analyses (active as primary, active as
secondary and deactive).
CFRICWORK
Field: yes History: no .fil: no
Contact frictional work for general contact analyses.
CNAREA
Field: yes History: no .fil: no
Contact nodal area for each node with active contact forces in general
展開 【公益帖】VUMAT 中間變量輸出到指定文件中的方法
背景： 在進行用戶材料本構用戶子程序調試過程中，我們需要輸出一些中間變量，跟蹤其演化規律，用于debug， 在umat中我們往往用 write（6，*） write（7，*）寫入到相應的msg或者dat文件中，但是在VUMAT中這種方法卻不可以，下面提供一種VUMAT輸出變量的方法
需要在vumat子程序中添加如下語句
open(1,file='絕對路徑\output.txt')
write(1,*) 變量
這樣會在指定路徑下生成一個名為output.txt，希望對大家有用
ABAQUS歷史輸出中，各能量變量（ALLAE、ALLIE、ETOTAL等）意義
ABAQUS中，對于很多動態問題，尤其像高速沖擊模擬中，對結果評價很重要的一點就是要保證模型能量守恒，這就涉及到ABAQUS歷史輸出中各能量變量的意義，下面最各簡單整理：
ALLAE：人工偽應變能，六面體、殼網格中沙漏發生情況指標
ALLCD：蠕變、膨脹以及粘彈性消耗的能量
ALLFD：摩擦消耗的能量
ALLIE：總的內能，ALLIE=ALLSE+ALLPD+ALLCD+ALLAE+ALLQB+ALLEE+ALLDMD
ALLKE：動能
ALLKL：碰撞消耗的能量
ALLVD：粘性消耗的能量
ALLDMD：破壞消耗的能量
ETOTAL：所有能量的總和
NOTE：ALLAE（偽應變能）的理解：偽應變能就是控制沙漏變形所消耗的主要能量。如果偽應變能過高，說明過多的應變能被用來控制沙漏變形了。一般通過比較偽應變能和其他內部能量的值來判斷偽應變能是否過高，以及判斷過高的偽應變能的來源。
展開 ParaView---FLAC3D和3DEC輸出變量到VTK文件
1 引言
在《免費二維有限元分析軟件ADONIS計算步驟》中曾經提及ADONIS的計算結果可以輸出為VTK文件，事實上3DEC和FLAC3D的計算結果也可以輸出為VTK文件。這個筆記簡要討論了輸出步驟。
2 VTK命令
3DEC和FLAC3D使用block vtk和zone vtk命令輸出vtk文件，可以輸出的變量有：
(1) displacements: 網格點的位移矢量
(2) zone-extra: 單元索引i的額外變量
(3) gridpoint-extra: 網格點索引i的額外變量
(4) forces: 網格點不平衡力的信息，包括最后的局部應力比【FLAC/Slope的強度降低過程(Strength Reduction Procedure in FLAC/Slope)】
(5) groups: slot中單元設置的組
(6) model-fluid: 流體屬性
(7) model-mechanical: 力學屬性
(8) model-thermal: 熱屬性
(9) pore-pressure: 網格點的孔隙壓力值
(10) range: 設置range, 用來控制哪些單元包括在vtk文件中
(11) saturation: 網格點飽和度的值
(11) stresses: 每個單元的平均應力。不輸出內部四面體的應力。
(12) temperatures: 網格點的溫度值
(13) velocities: 網格點的速度矢量
(14) filename: 設置vtk的文件名
默認情況下，只輸出應力(stresses)，組(groups)和位移(displacements). vtk的文件格式是XML語言。
展開 【LSDYNA操作小技巧七】LSDYNA中具體材料歷史變量的定義與輸出問題 ￥1.99
LSDYNA作為專門針對材料非線性、大變形問題而開發的大型顯示動力學仿真軟件，其擁有相當齊全的材料庫，針對其中材料的本構關系都有非常明確的軟件設置與輸出步驟，方便用戶對本構的快速設置甚至二次開發。本帖以典型的材料大變形問題磨粒切削加工引起的加工損傷為案例講解損傷歷史變量的定義與輸出步驟。
在我們沒有更好的文獻參考時，官方仿真軟件的幫助文檔是一個最佳的選擇。以下提供在LSDYNA中具體查閱歷史變量定義的鏈接如下，
History Variables for Certain Material Models — Welcome to the LS-DYNA support site.
基于此，必須首先定義目標歷史變量如圖1所示。具體通過，在關鍵字 *DATABASE_EXTENT_BINARY中定義NEIPH 或 NEIPS 實現。這里NEIPH 或 NEIPS 只能填入具體的數字加以定義材料具體哪些的歷史變量。
圖1
具體以此為例：為輸出磨粒切削加工引起的工件損傷，首先找到對應的工件材料編號110（工件為K9玻璃，采用JH-2本構），可知損傷因定義為2如圖1所示，故在NEIPH輸入數字2即可，若為了查看其他材料歷史變量，對應輸出1-4即可實現。
圖2
完成定義后，完成模型建模輸出為k文件，在LSPP中進行損傷云圖的輸出。圖3為定義的損傷在后處理中的查看方式，圖4為輸出的損傷云圖，紅色粒子表示損傷。
圖3
圖4
以上，基于此案例，可在LSDYNA軟件中定義并輸出所給材料的歷史變量。
展開 【ANSYS文件輸出】ANSYS怎么批量文件輸出？
這期是ANSYS不定期更新，這期主要講一下ANSYS怎么批量文件輸出，還有批量輸出前的一些預處理操作。</span></p><p><strong>01</strong></p><p><strong>怎么去除右上角的時間?</strong></p><p>/PLOPTS,DATE,0 !都不顯示</p><p>/PLOPTS,DATE,1 !只顯示日期</p><p>/PLOPTS,DATE,2 !顯示日期和時間</p><p><strong>02</strong></p><p><strong>怎么去除右上角的LOGO？</strong></p><p>/plopts,logo,on !以logo形式</p><p>/plopts,logo,off !以文本形式</p><p> 用上面兩種形式，logo還是會存在的，還有一種方法就是在界面，鼠標箭頭放在logo上，右鍵點擊data，REPLOT一下，這樣logo就會不再出現。</p><p><strong>03</strong></p><p><strong>怎么把背景變成白色？</strong></p><p>!
展開 
LS-DYNA 復合材料層合板每層應力輸出方法和損傷云圖顯示！ ￥60
一個多層的shell 在球狀物體的沖擊下，共定義20層，現在想輸出：
（1）20層每層上的應力
（2）復合材料損傷云圖
基本的準備工作
（1）定義復合材料層合板（K文件中采用*Part_composite進行
*Part_composite
（2）在*DATABASE_EXTENT_BINARY選項中要設定 MAXINT=20 用來存儲20層的應力等變量
（3）查看損傷 歷史變量有如下6個變量，如果要輸出損傷變量NEIPS=6
提取某單元厚度方向應力的方法：
其實只是顯示方法的問題，采用如下方法就可以輸出全部20個層的數據
history----》int pt--->選中單元——》plot
損傷云圖用var6進行標定，如下圖所示
本案例的k文件 以及查看損傷云圖的方法見付費文檔
展開 ANSYS中的變量
ANSYS中的變量總的來說分為兩大類：一類是標量參數；一類是數組參數。其意義與C語言中的參數和數組類似。標量參數是指單個的變量，而數組則是由一系列具有相同意義的數據組成。
一：標量參數
標量參數分為兩種：數值型的和字符型的。
1：標量參數的定義或修改：可以使用*SET命令或直接寫出賦值表達式。修改時，只需改變數據即可。
例如：
*SET,LENGTH,10 !定義了一個名為LENGTH的數值型變量，并且值為10
*SET,COMMENT,'DISP' !定義了一個名為COMMENT的字符型變量，并且賦值為DISP
以上兩句等價于：
LENGTH=10
COMMENT='DISP'
2：標量參數的刪除：同樣可以使用*SET命令或直接寫出賦值表達式或使用*DEL命令。將需要刪除的變量賦空值，注意不是數字0或空格。
例如：我們將上面定義的兩個參數刪除。語句如下：
*SET,LENGTH,
*SET,COMMENT,
或：
LENGTH=
COMMENT=
又或：
*DEL,LENGTH
*DEL,COMMENT
3：標量參數的GUI操作菜單。GUI路徑為：Utility Menu->Parameters->Scalar Parameters。標量參數的定義，刪除和修改都可以在這個窗口中完成。
二：數組參數
數組參數按維數可分為：一維數組，二維數組和三維數組。
按存儲的數據類型可分為：
1）一般數組參數，也成為ARRAY Parameter。是一種默認的數據類型，也就是說，如果用戶不聲明數據類型，系統就會自動使用該數據類型。該數組的行、列、面的索引均為從1開始的連續整數，元素值可以為整數或實數。
2）字符數組參數，又名CHAR Array Parameter。元素值為不多于8個的字符或數字組成。其行，列，面的索引值也是從1開始的連續整數。
展開 ANSYS APDL如何求變量的最大值
在其他語言中求最大值非常容易，比如有三個變量分別是MXS6，MXS7和MXS8，要求他們的最大值賦予MaxS，用到的函數往往只是一個函數MaxS=max(MXS6,MXS7,MXS8)。但是在ANSYS Parametric Design Language Guide第三章第8節 Parametric Functions部分卻找不到這樣的函數。
無奈自力更生吧。
MaxS=MXS6
*if,MXS7,GT,MaxS,THEN
MaxS=MXS7
*ENDIF
*if,MXS8,GT,MaxS,THEN
MaxS=MXS8
*ENDIF
這是三個變量，如果有更多變量方法也是一樣的，變量太多就用個循環，如果需要留言給我，我給大家把代碼寫出來。
展開 【ANSYS經驗技巧】- 學會使用變量數組函數（轉載）
5.直到你為函數中所有的狀態的所有變量提供賦值，才能保存為表格式參數使用。
注意：在ansys分析中用函數加載必須的兩個步驟：
1利用函數編輯器創建任意方程或函數
2在利用函數加載器加載函數，并以函數定義表參數
轉自：三維網。 作者：4kpolo
原帖鏈接：http://www.3dportal.cn/discuz/viewthread.php?tid=746468&extra=page%3D1%26amp%3Bfilter%3Ddigest
