Abaqus輸出矩陣的方法

2017-03-18 by:CAE仿真在線來源:互聯網

在比較深入的CAE分析過程中,我們可能會關注模型中某些部分的剛度分布情況,因此需要提取剛度矩陣參數,因為使用相對較少,目前的Abaqus前處理模塊CAE中還無法直接實現該設置,因此需要對其inp文件進行局部修改來定義。

Abaqusn輸出剛度矩陣主要分為兩大類,一類是輸出某一部分的單元的剛度矩陣,一類是輸出整體模型的剛度矩陣(可含單元、彈簧、質量、阻尼、力等)。

下面我們分別對這兩中情況分別討論:

1、輸出單元矩陣

輸出單元剛度矩陣的方式相對容易理解,當然你也可以從自帶的英文幫助中需找答案,你只需要在inp文件中添加下述關鍵字即可:

*File Format,ASCII
*Element Matrix Output,Elset=xxx,

File Name=xxx,Frequency=1,Output File=User Defined,Stiffness=Yes

其中第一行的目的是將軟件默認生成的二進制文件轉化為ASCII格式文件,以方便直接查看和讀取避免亂碼。Elset后面的xxx對應的是之前已經定義好的單元集的名稱,也就是你想要獲取的那部分單元。File Name后面的xxx是abaqus計算后所提取出的單元矩陣所保存的文件名稱,后綴為.mtx ,也就是說最終在工作目錄中會生成一個xxx.mtx的文件,可以通過記事打開的剛度矩陣文件。

2、輸出總體矩陣

對于整體矩陣參數的輸出,跟單元矩陣的關鍵字設置有細微的差別,我們通過軟件自帶的標準模板進行解釋:

**********************************************************

*HEADING
…
**
*STEP
Options to define the preloading history for the model.
*END STEP
********************************************************************
*STEP
*MATRIX GENERATE, STIFFNESS, MASS, VISCOUS DAMPING,STRUCTURAL DAMPING, LOAD
*MATRIX OUTPUT, STIFFNESS, MASS, VISCOUS DAMPING,STRUCTURAL DAMPING, LOAD, FORMAT=MATRIX INPUT
*BOUNDARY
Options to define the boundary conditions for the matrix generation step.
**
*LOAD CASE, NAME=LC1
Options to define the loading for the first load case.
*END LOAD CASE
**
*LOAD CASE, NAME=LC2
Options to define the loading for the second load case.
*END LOAD CASE
Any number of load cases can be defined.
*END STEP

************************************************************

如過不考慮輸出整體矩陣,正常情況下,所設置的參數在紅色之前即結束,計算后即可獲取應力應變變形等結果。但為了獲取整體矩陣參數,就需要補充紅色的關鍵字部分,實際上(個人理解而已),紅色添加的部分即是你關注的那個分析步計算后的整體矩陣,只需要把之前的分析步中的邊界條件和工況設置復制到最后,并更改分析部名稱和添加整體矩陣關鍵字到相應位置(其中的藍色部分)即可。

對于整體矩陣的輸出,包括剛度、質量、阻尼、力等,具體的需要輸出哪些可以自己的情況調整。比如你想要獲取剛度和力的整體矩陣,關鍵字的調整可以如下:

*STEP, NAME=STEPMLC1

*MATRIX GENERATE, STIFFNESS

*MATRIX OUTPUT, STIFFNESS,LOAD, FORMAT=MATRIX INPUT

*Boundary

_PickedSet5, 1, 1

_PickedSet5, 2, 2, -5.

_PickedSet5, 3, 3

_PickedSet5, 4, 4

_PickedSet5, 5, 5

_PickedSet5, 6, 6

*END STEP