abaqus中輸出設置的案例
LS-DYNA中的操作及設置（三）（力的輸出，重力載荷，鉸鏈剛度）
本文翻譯自官方文檔，原文鏈接：
https://www.dynasupport.com/howtos/general
一、力的輸出（Force Analysis）
可以通過一下幾種方式獲得受力情況：
1.使用*DATABASE_CROSS_SECTION 和 *DATABASE_SECFORC可以獲得一個橫截面上的內力和內力矩。注意，在使用set選項設置橫截面時，必須提供用于定義橫截面路徑的節點集以及橫截面某一側的至少一個單元集。
2.使用*DATABASE_NODAL_FORCE_GROUP 和 *DATABASE_NODFOR可以獲得施加在節點或節點集上的外部力。
3.使用*DATABASE_BNDOUT可以獲得邊界反力。
4.使用*DATABASE_SPCFORC可以獲得由*BOUNDARY_SPC定義的邊界上的反力。
5.*DATABASE_RCFORC可以輸出接觸力合力。在單面接觸中，還需要使用*CONTACT_FORCE_TRANSDUCER_(option)來提取接觸力合力。
6.單個節點所受的力由*DATABASE_NCFORC輸出。必須在*CONTACT中設置SPR, MPR中至少一個輸出控制選項為1，以便向程序指定輸出側為接觸面主面側或從面側。
以上數據分別由程序輸出到ASCII文件SECFORC, NODFORC, BNDOUT, SPCFORC, RCFORC, 和 NCFORC中。
展開 ABAQUS中輸出Cohesive單元的斷裂形式
ABAQUS中其實已經提供了識別這兩種（拉伸和剪切）破壞形式的場輸出：MMIXDME、MMIXDMI。下面從定義、使用、結果演示等幾個方面進行介紹。
定義
幫助文檔中，二者定義如下：
MMIXDME，為損傷演化過程中混合斷裂模式的比例，定義為1-m1，位于單元積分點位置，單元未破壞時，其數值默認為-1；
MMIXDMI，為初始損傷時混合斷裂模式的比例，同樣定義為1-m1，位于單元積分點位置，單元未破壞時，其數值默認為-1；
后面看下m1的定義：
m1為Gn（一型張開斷裂能）與GT（三型斷裂能之和）的比值：
當m1為1，表示完全的拉伸破壞時，對應的MMIXDME和MMIXDMI數值為0；
當m1為0，表示完全的剪切破壞時，對應的MMIXDME和MMIXDMI數值為1；
通過這兩個參量，我們就可以判斷cohesive單元以哪種形式破壞為主。
使用方法
目前不支持GUI界面定義，Step模塊的場輸出無法找到MMIXDMI和MMIXDME。
展開 Abaqus中溫度輸出的規定 ￥10
Abaqus中溫度輸出的規定，詳細解釋了abaqus中實體單元，殼單元，梁單元溫度輸出的規定，并用實例進行了展示。
ABAQUS腳本將計算結果輸出到文本中 ￥1.99
ABAQUS的結果大多以字典存儲，之前的帖子講了怎么將結果輸出到Excel中，此處講述如何將多個數據輸出到文本文件中，方便用于其他目的。
附件中的結果文件是一個L型梁，輸出的數據包括單元ID，應變能、mises應力和應變等，如果不做額外處理輸出的結果很不整齊，看起來繁雜。為此在使用python腳本讀取結果數據時，在腳本中專門編寫了一個文本排列函數，如下所示：
def fb(character, identifier, length):
if identifier == 'f':
sentence = character + ' '*(length - len(character))
else:
sentence = ' '*(length - len(character)) + character
return sentence
輸出結果如下：
腳本在附件中，需要的朋友可以嘗試下載。
展開 
Abaqus動態分析中，如何快速查看整個響應過程中場輸出結果的最值 ￥9.9
<p>需求：動態分析（基于模態的瞬態動態響應分析、顯示動態分析等）中結果的響應也是一個動態的過程，不確定哪個時刻的結果是最大值或者最小值，或者說想知道整個響應過程中的最大值、最小值是多少。結果輸出中是不會直接輸出的，只能看到每幀場輸出中的最值，又不可能自己逐幀場輸出結果里去看，然后找到所有幀中的最值，那么Abaqus軟件內如何實現呢？</p><p><br></p><p><span style="background-color: rgb(255, 255, 255); color: rgb(25, 27, 31);">原創聲明：未經本人同意，禁止抄襲、二次創作及轉載！</span></p>
展開 abaqus 獲取odb場輸出中某一集合的數據
獲取odb中某一個集合的數據，有代碼實現和cae實現
代碼實現：主要的是一個getSubset的函數，可以參考abaqus手冊
https://help.3ds.com/2022/english/dssimulia_established/SIMACAECMDRefMap/simacmd-c-odbintroreadregpyc.htm?
ABAQUS如何輸出分析過程中每一時刻的最大溫度曲線
ABAQUS如何輸出分析過程中每一時刻的最大溫度曲線，或者輸出分析過程中每一時刻最大溫度值，如下圖所示
微信截圖_20180823102044.jpg
Abaqus二次開發讀取變形后的節點坐標并輸出到txt文件中。讀取模態頻率到txt 文件中。 ￥10
有時候在abaqus中，我們需要知道某一個集合中的節點變形后的坐標。以此為輸入來進行一些研究。這里我用一個自編函數，將節點集合變形后的坐標寫入到一個txt文件中。格式化的寫入文件，方便用此坐標來進行計算研究。
同樣有時候，需要將計算的模態頻率值提取出來。同樣用一個函數將模態頻率提取出來，放進txt文件中，方便后續研究。
讀取的節點結果如下圖所示：
ABAQUS中隱式和顯式的節點和單元的輸出變量解析
FSLIP
Field: yes History: no .fil: no
Length of contact slip path at secondary nodes during contact (FSLIPEQ) and in some
cases (see About contact pairs in Abaqus/Explicit) components of net contact slip in local
tangent directions (FSLIP1 and FSLIP2). These variables remain constant while a
secondary node is not in contact.
FSLIPR
Field: yes History: no .fil: no
Magnitude of contact slip rate at secondary nodes during contact (FSLIPR) and in some
cases (see About contact pairs in Abaqus/Explicit) components of contact slip rate in local
tangent directions (FSLIPR1 and FSLIPR2). These variables are set to zero while a
secondary node is not in contact.
BONDSTAT
Field: no History: yes .fil: no
Spot weld bond status.
展開 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（偽應變能）的理解：偽應變能就是控制沙漏變形所消耗的主要能量。如果偽應變能過高，說明過多的應變能被用來控制沙漏變形了。一般通過比較偽應變能和其他內部能量的值來判斷偽應變能是否過高，以及判斷過高的偽應變能的來源。
展開 abaqus中阻尼的設置
Abaqus阻尼設置方式 abaqus的阻尼分為兩類，與速度成比例的粘性阻尼；和與位移成比例的結構阻尼（在頻域分析中采用） abaqus引入阻尼的3中途徑： 材料和單元的阻尼 整體阻尼，包括粘性阻尼，瑞利阻尼，結構阻尼 模態阻尼，只能用于模態分析 在ABAQUS中阻尼可以應用在下面的動力學分析中： 非線性問題直接積分求解（顯式分析或者隱式分析） 直接法或子空間法穩態動力學分析 模態動力學分析（線性） 4. Abaqus阻尼設置 - 具體操作 針對模態動力學分析，在ABAQUS/Standard中可定義幾種不同類型的阻尼：直接模態阻尼（DirectModal Damping），瑞利阻尼（RayleighDamping），復合模態阻尼（Composite Modal Damping）和結構阻尼（StructureDamping）。 ABAQUS動力學分析中用*Modal Damping選項來定義阻尼。 以下內容是以在step分析步內定義阻尼的舉例，每階模態可以定義不同量值的阻尼，但其實也可以在Material分析步設置阻尼。 4.1 直接模態阻尼： 采用直接模態阻尼可以定義對應于每階模態的阻尼比。其典型的取值范圍是在臨界阻尼的1%~10%之間。直接模態阻尼允許用戶精確定義系統的每階模態的阻尼。 舉例：設置前10階振型的阻尼定義為4%的臨界模態阻尼，11~20階振型的阻尼為5%的臨界阻尼 界面操作：在分析步驟內定義直接模態阻尼，如下圖所示，激活直接模態阻尼選項（DirectModal），并在數據行內輸入數據。 4.2 瑞利阻尼 - 常用，需重點掌握 在瑞利阻尼中，假設阻尼矩陣可表示為質量矩陣和剛度矩陣的線性組合，即： 在Material分析步設置阻尼 ABAQUS中通過設置alpha和beta來求解瑞利阻尼，具體如下圖。
展開 
讀取ABAQUS結果文件中的單元應變能并輸出至excel文件的處理腳本 ￥1
腳本內容如題。
ABAQUS中接觸設置小問題總結
對于法向作用，ABAQUS中接觸壓力和間隙的默認關系是“硬接觸”（HARD CONTACT）。硬接觸：接觸面之間能夠傳遞的接觸壓力大小不受限制；當接觸壓力變為零或負值時，兩個接觸分離，并且去掉相應節點上的接觸約束。
對于切向作用，ABAQUS中常用的摩擦模型為庫侖摩擦，即使用摩擦系數來表示接觸之間的摩擦特性。默認的摩擦系數為零，即無摩擦。庫侖摩擦的計算公式為
τ=μ x p
τ是臨界切向力，μ是摩擦系數，p 是法向接觸壓強（CPRESS）。在切向力達到臨界切應力之前，摩擦面之間不會發生相對滑動。
主面/從面選擇一般性要求：
（1）小面為從面；?
（2）選擇剛度較大的面；
（3）相同剛度（不要僅從彈性模量角度，要從剛度角度，比如大E的殼比小E的塊要軟）的選擇單元粗的作為主面；
（4）主面不能是由節點構成的面，并且必須是連續的。如果是有限滑移，主面在發生接觸的部位必須是光滑的，不能有尖角。
有限滑移和小滑移
有限滑移（FINITE SLIDING） 兩個接觸面之間可以有任意的相對滑動和裝懂，這是定義接觸時的默認特性。在有限滑移的過程中，ABAQUS/STANDARD需要實時地判定從面節點和主面的哪一部分發生接觸，因此計算代價較大。一般在滑移量大于單元尺寸或者需要精確分析時，采用有限滑移。
小滑移（SMALL SLIDING） 兩個接觸面之間只有很小的相對滑動（小于單元尺寸的20%）。對于小滑移的接觸時，ABAQUS/STANDARD在分析的開始就確定了從面節點和主面的哪一部分發生接觸，在整個分析過程中這種接觸關系不會再發生變化。
展開 用ANSA設置ABAQUS中的過盈接觸
ANSA作為一款功能強大的前處理軟件，為ABAQUS高度非線性的過盈接觸提供了全面的支持。
2.ANSA中設置ABAQUS過盈接觸示例
過盈接觸屬于緊配合的一種，在工程中應用相對廣泛。在孔與軸，鍵槽與鍵等結構之間，都可能使用了過盈接觸。ABAQUS具有三種定義過盈接觸的方式，ANSA均提供了對應的支持。本文以上圖中的兩個環形結構為例，簡要介紹使用ANSA為ABAQUS設置過盈接觸的方法。
（1） 幾何過盈
在初始的裝配過程中，使結構出現適當的干涉（即過盈量），然后為過盈接觸的兩個面定義一般的接觸關系即可。
（2） 關鍵詞CLEARANCE
模型網格劃分完成后，打開ABAQUS DECK面板，在AUXILIARIES菜單下選擇CONTACT打開CONTACT菜單，單擊右鍵或直接在菜單底部選擇New->CONTACT，創建新的接觸對。
打開接觸對定義窗口，選擇接觸的主、從面ID，以及接觸屬性ID。可以在ANSA相應的輸入域中鍵入“?“，以打開相應輸入項的編輯定義窗口，對相應的項進行創建、編輯或確認。
這種方式只適用于小滑移這種點的相對滑動方式。在接觸類型下選擇點-面接觸或面-面接觸的接觸類型，相應的設置方法與一般接觸設置相同。此外，需要將接觸容差打開，并對其設定恰當的值。考慮到過盈接觸的特性，接觸容差值應為小于零的負值。
（3） 關鍵詞CONTACT INTERFERENCE
與（2）相同，在模型網格劃分完成后，打開ABAQUS DECK面板，在AUXILIARIES菜單下選擇CONTACT打開CONTACT菜單，單擊右鍵或直接在菜單底部選擇New->CONTACT，創建新的接觸對。
展開 質量點在Abaqus中的設置
簡單算例一個，看下圖所示：
上圖是一個10×10×10 的立方體，密度為1，施加X方向的初速度1，在interaction模塊下檢查其結構體的相關質量速度，可得到相關的質量參數， 1000
理論可得出1/2*m*V*V ，得到相關的動能為500 在后處理中可以查詢其相關的動力學信息。
加入質量為1000的質量塊：所有的操作都在intercation模塊下的完成
1、加入參考點RF
2、在參考點上加入質量點
Tools-special-Create-Point mass/intertia
3、將新加入的質量點 Coupling到一個相應的面之上 ，否則該質量點懸空在！
一定要記得將在初始速度中考慮該參考點，重新計算，可得到 動能 1000，結果沒有問題。
Note：在Abaqus CAE下如果通過quere Mass property，是不考慮質量點的，所以查詢的時候，顯示的仍然只有500，切記！ 要檢查質量，請開打status 文件，會有質量的提示。
展開 