實用子程序SPRINC

在UMAT子程序k">ABAQUS中使用UMAT子程序時有時會使用到最大主應力進行計算。通過查閱幫助文檔，ABAQUS實用子程序SPRINC可以在UMAT中計算最大主應力和最大主應變，SPRIND可以計算最大主應力和最大主應變的方向。

 

下面是ABAQUS幫助文檔關于實用子程序SPRINC的介紹：

SPRINC (calculate principal values)

 

Interface

     CALL SPRINC(S,PS,LSTR,NDI,NSHR)

Variables to be provided to the utility routine

S

Stress or strain tensor.

LSTR

An identifier. LSTR=1 indicates that S contains stresses; LSTR=2 indicates that S contains strains.

NDI

Number of direct components.

NSHR

Number of shear components.

Variables returned from the utility routine

PS(I), I=1,2,3

The three principal values.

 

 

調(diào)用的時候直接在UMAT子程序中輸入：CALL SPRINC(S,PS,LSTR,NDI,NSHR)

需要提供給該實用子程序的變量包括：

實用子程序返回UMAT的變量

實例：

CALL SPRINC(S,PS,LSTR,NDI,NSHR) 可以直接寫成

CALL SPRINC(STRESS,PS,1,3,3)或CALL SPRINC(STRESS,PS,1, NDI,NSHR)

這行代碼的意思是將UMAT中的應力分量數(shù)組STRESS提供給實用子程序SPRINC計算；1表示LSTR=1，表明提供給實用子程序SPRINC的為應力分量數(shù)組；后面的兩個“3”分別表示正應力分量個數(shù)為3和切應力分量個數(shù)為3，也可以不用改變直接使用NDI和NSHR，因為在三維模型的UMAT中NDI和NSHR已經(jīng)被賦予了值為3。PS數(shù)組儲存返回的計算得到的主應力，這個需要自己在umat程序中定義數(shù)組PS.

下面是關于返回主應力順序的驗證：

     CALL SPRINC(STRESS,PS,1,3,3)

     STATEV(4)=PS(1)

     STATEV(5)=PS(2)

     STATEV(6)=PS(3)

需要注意以下點：

1.PS數(shù)組一般需要在程序中自定義數(shù)組變量；

2.S數(shù)組一般也需要自定義數(shù)組變量，且需要為S數(shù)組賦予應力分量或應變分量的值；

3.如果SPRINC或SPRIND返回的主應力為零，檢查輸入的變量或數(shù)組是否未定義或者未賦值，其次檢查輸出數(shù)組是否未定義。