UMAT全過程技術(shù)篇

UMAT全過程-"技術(shù)篇"[寫在前面:

這篇文章是UMAT全過程-感想篇的姊妹篇,容許要給大家寫的一篇帖子,同時也是為了記錄自己的學習過程,與大家分享!首先指出,俺的"技術(shù)篇"--是加了引號的,因為確實稱不上有多么大的技術(shù)含量,還望大家莫笑偶!只不過一是跟那個感想篇形成一個對照,同時主要內(nèi)容為自己編子程序過程中涉及的技術(shù)邊邊上的小問題的一些解決方法,供仿友們參考!偶不是謙虛,也不是一個低調(diào)的人,大家謝謝和支持的話,我先行謝過啦!更希望大家能提出質(zhì)疑或者別的更好的方法,大家相互交流,共同進步!]----------------------------------------------------------------------------*轉(zhuǎn)*入*正*題*第一局部:相關知識[特別聲明,這局部來自于華中科技大學楊曼娟同學的碩士學位論文,在此對作者表示感謝!--大家可以去知網(wǎng)下載]----------------------------------------------------------------------------1.ABAQUS中材料非線性問題的處理

ABAQUS中材料非線性問題用Newton-Raphson法來求解。首先將載荷分為假設干個微小增量，結(jié)構(gòu)受到一個微小增量△P。ABAQUS用與初始結(jié)構(gòu)位移相對應的初始剛度矩陣K0和荷載增量△P計算出結(jié)構(gòu)的在這一步增量后的位移修正Ca、修正后的位移值Ua和相應的新的剛度矩陣Ka。ABAQUS用新的剛度矩陣計算結(jié)構(gòu)的內(nèi)力Ia，荷載P和Ia的差值為迭代的剩余力Ra，即Ra=P-Ia。如果Ra在模型內(nèi)的每個自由度上的值都為零，如圖2-2中的a點，那么結(jié)構(gòu)處于平衡狀態(tài)。但在非線性問題中，通常Ra是不可能為零，ABAQUS為此設置了一個剩余力容差。如果Ra小于這個數(shù)字，ABAQUS就認為結(jié)構(gòu)的內(nèi)外力是平衡的。一般這個缺省值取為平均內(nèi)力的0.5％(如圖2-2)。

在Ra滿足小于剩余力容差的條件后，ABAQUS還要檢查位移修正系數(shù)Ca是否遠遠小于結(jié)構(gòu)的位移增量△Ua=Ua-Uo.如果大于位移增量的1％，ABAQUS將自動繼續(xù)進行迭代。第二次迭代時，ABAQUS用剛度矩陣Ka和第一次迭代的剩余力Ra計算，得到一個新的位移修正值。同理，第二次迭代中ABAQUS用新的剩余力、新的位移系數(shù)，繼續(xù)驗證收斂性，直到收斂為止〔如圖2-3〕。因此非線性問題的計算量要比線性問題的計算量大得多，計算的時候需要分配更多的內(nèi)存和磁盤空間。

ABAQUS會自動調(diào)整荷載增量的大小，用戶只需合理的設置一個第一次迭代的增量值即可，如果用戶不給出第一次迭代的增量值，ABAQUS會在第一次迭代的過程中把設置的所有荷載都加到結(jié)構(gòu)上去，然后不斷的試算自動調(diào)整第一次迭代的增量大小。尋找一個增量值的收斂解迭代的次數(shù)取決于系統(tǒng)非線性的程度。ABAQUS中缺省的設置是，如果在一個增量值的作用下，迭代16次以后結(jié)果仍不收斂，ABAQUS將放棄這個增量值，把增量取為前一次增量值的25％再進行計算。如果ABAQUS在迭代次數(shù)小于5次的情況下就找到了收斂解，ABAQUS就自動把增量值增加50％再進行計算，得到更為合理的增量值。2.用戶子程序接口

ABAQUS不僅提供標準的有限元分析程序，而且具有良好的開放性，可利用它

提供的用戶子程序接口生成非標準的分析程序來滿足用戶的需要，在實際工程中得

到廣泛應用。ABAQUS允許用戶通過子程序以代碼的形式來擴展主程序的功能，并

給用戶提供了強大而又靈活的用戶子程序接口和應用程序接口〔UtilityRoutine〕，

ABAQUS6.4一共有42個用戶子程序接口，15個應用程序接口，用戶可以定義包括

邊界條件、荷載條件、接觸條件、材料特性以及利用用戶子程序和其它應用軟件進

行數(shù)值交換等等。這些用戶子程序接口使得用戶解決一些問題時有很大的靈活性，

同時大大的擴充了ABAQUS的功能。例如：

通過用戶定義單元接口，用戶自定義的任何類型的線性或非線性單元都可以被引入模型中，對于線性單元剛度矩陣和質(zhì)量矩陣可以直接確定。同時，用戶子程序也可被用來定義這些單元的線性和非線性特性。通過用戶材料子程序接口，用戶可定義任何補充的材料模型，不但任意數(shù)量的

材料常數(shù)都可以作為資料被讀取，而且ABAQUS對于任何數(shù)量的與解相關的狀態(tài)變

量在每一材料計數(shù)點都提供了存儲功能，以便在這些子程序中應用。3.用戶子程序和主程序的結(jié)合

ABAQUS的用戶子程序是根據(jù)ABAQUS提供的相應接口，按照Fortran語法，用戶自己編寫的代碼。它是一個獨立的程序單元，可以獨立的被存儲和編譯，也能被其它程序單元引用，因此，利用它可帶回大量數(shù)據(jù)供引用程序使用，也可以用它來完成各種特殊的功能。它的一般結(jié)構(gòu)形式是：SUBROUTINES(x1,x2,……,xn)

INCLUDE‘ABA_PARAM.INC’〔用于ABAQUS/Standard用戶子程序中〕

ORINCLUDE‘VABA_PARAM.INC’)〔用于ABAQUS/Explicit用戶子程序中〕

……

RETURN

END

x1，x2，……，xn是ABAQUS提供的用戶子程序的接口參數(shù)，有些參數(shù)是ABAQUS傳到用戶子程序中的，例如SUBROUTINEDLOAD中的KSTEP、KINC、COORDS；有些是需要用戶自己定義的，例如F。文件aba_param.inc和vaba_param.inc隨著ABAQUS軟件的安裝而包含在操作系統(tǒng)中，它們含有重要的參數(shù)，幫助ABAQUS主求解程序?qū)τ脩糇映绦蜻M行編譯和鏈接。當控制遇到RETURN語句時便返回到引用程序單元中去，END語句是用戶子程序結(jié)束的標志。

在一個算例中，用戶可以用到多個用戶子程序，但必須把它們放在一個以.for為

擴展名的文件中。運行帶有用戶子程序的算例同時有兩種方法：一是在CAE中運行，

在EDITJOB菜單中的GENERAL子菜單的USERSUBROUTINEFILE對話框中選擇用戶子程序所在的文件即可；另外是在ABAQUS.COMMAND中運行，語法如下：

abaqusjob=job-nameuser={source-file|object-file}編制用戶子程序時應注意:(1)用戶子程序相互之間不能調(diào)用，但可以調(diào)用用戶自己編寫的Fortran子程序

和ABAQUS應用程序。ABAQUS應用程序必須由用戶子程序調(diào)用。當用戶編寫

Fortran子程序時，建議子程序名以K開頭，以免和ABAQUS內(nèi)部程序沖突。

(2)當用戶在用戶子程序中利用OPEN翻開外部文件時，要注意以下兩點：一是設備號的選擇是有限制的，只能取15～18和大于100的設備號，其余的都已被ABAQUS占用；二是用戶需提供外部文件的絕對路徑而不是相對路徑。

(3)對于不同的用戶子程序ABAQUS調(diào)用的時間是不同的，有的是在每個STEP

的開始，有的是STEP的結(jié)尾，有的是在每個INCREMENT的開始等等。當ABAQUS

調(diào)用用戶子程序時，都會把當前的STEP和INCREMENT利用用戶子程序的兩個實

參KSTEP和KINC傳給用戶子程序，用戶可把它們輸出到外部文件中，這樣就可清

楚的知道ABAQUS何時調(diào)用該用戶子程序。為保證用戶子程序的正確執(zhí)行，子程序的書寫必須遵循ABAQUS的相關規(guī)定.4.用戶材料子程序UMAT接口的原理

用戶材料子程序〔User-definedMaterialMechanicalBehavior，簡稱UMAT〕是

ABAQUS提供應用戶定義自己的材料屬性的Fortran程序接口，使用戶能使用ABAQUS材料庫中沒有定義的材料模型。用戶材料子程序UMAT通過與ABAQUS主求解程序的接口實現(xiàn)與ABAQUS的資料交流.在輸入文件中，使用關鍵詞“*USERMATERIAL〞表示定義用戶材料屬性。[更詳細相關知識,請看帖子ABAQUS-UMAT-自學知識整理貼[已經(jīng)初步完成,不斷完善更新,請跟帖討論]]

由于UMAT子程序在單元的積分點上調(diào)用，增量步開始時，主程序路徑將通過

UMAT的接口進入UMAT，單元當前積分點必要變量的初始值將隨之傳遞給UMAT

的相應變量。在UMAT結(jié)束時，變量的更新值將通過接口返回主程序。

ABAQUS調(diào)用UMAT過程及次數(shù)如圖2-4所示。5.UMAT子程序流程

UMAT子程序采用Fortran語言編制，可以包括以下幾個局部：子程序定義語句、ABAQUS定義的參數(shù)說明、用戶定義的局部變量說明、用戶編制的程序主體、子程序返回和結(jié)束語句。

主要求解過程：每一個增量加載步開始時，ABAQUS主程序在單元的積分點上調(diào)用UMAT子程序，并傳入應變增量、時間步長及荷載增量，同時也傳入當前狀態(tài)的應力、應變及其它與求解過程相關的變量；UMAT子程序根據(jù)本構(gòu)方程求解應力增量并更新應力及其它相關的變量，提供Jacobian矩陣給ABAQUS主程序以形成整體剛度矩陣；主程序結(jié)合當前荷載增量求解位移增量，繼而進行平衡校核；如果不滿足指定的誤差，ABAQUS將進行迭代直到認為收斂，然后進行下一增量步的求解。

ABAQUS中應力σ采用Cauchy應力張量描述，剪應變分量εij按照工程剪應變

的定義存儲。在求解的過程中，需記錄保存與求解過程相關的變量，這可以儲存在UMAT子程序的STATEV數(shù)組中，并且可以定義儲存的數(shù)量。對本文所述的本構(gòu)方程，定義了13個與求解過程相關的變量，即6個彈性應變，6個塑性應變，1個塑性參數(shù)增量。非非的ABAQUS用戶子程序[zz]轉(zhuǎn)自當用到某個用戶子程序時，用戶所關心的主要有兩方面：一是ABAQUS提供的用戶子程序的接口參數(shù)。有些參數(shù)是ABAQUS傳到用戶子程序中的，例如SUBROUTINEDLOAD中的KSTEP，KINC，COORDS；有些是需要用戶自己定義的，例如F。二是ABAQUS何時調(diào)用該用戶子程序，對于不同的用戶子程序ABAQUS調(diào)用的時間是不同的。有些是在每個STEP的開始，有的是STEP結(jié)尾，有的是在每個INCREMENT的開始等等。當ABAQUS調(diào)用用戶子程序是，都會把當前的STEP和INCREMENT利用用戶子程序的兩個實參KSTEP和KINC傳給用戶子程序，用戶可編個小程序把它們輸出到外部文件中，這樣對ABAQUS何時調(diào)用該用戶子程序就會有更深的了解?！沧映绦蛑泻苤匾木褪且烙蒩baqus提供的那些參量的意義，如下〕首先介紹幾個子程序：一．SUBROUTINEDLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,COORDS,JLTYP,SNAME)參數(shù)：

1．F為用戶定義的是每個積分點所作用的荷載的大?。?/p>

2．KSTEP,KINC為ABAQUS傳到用戶子程序當前的STEP和INCREMENT值；

3．TIME(1),TIME(2)為當前STEPTIME和INCREMENTTIME的值；

4．NOEL,NPT為積分點所在單元的編號和積分點的編號；

5．COORDS為當前積分點的坐標；

6．除F外，所有參數(shù)的值都是ABAQUS傳到用戶子程序中的。功能：

1．荷載可以被定義為積分點坐標、時間、單元編號和單元節(jié)點編號的函數(shù)。

2．用戶可以從其他程序的結(jié)果文件中進行相關操作來定義積分點F的大小。例1：這個例子在每個積分點施加的荷載不僅是坐標的函數(shù)，而且是隨STEP變化而變化的。

SUBROUTINEDLOAD(P,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,COORDS,

1JLTYP,SNAME)

INCLUDE'ABA_PARAM.INC'C

DIMENSIONTIME(2),COORDS(3)

CHARACTER*80SNAME

PARAMETER(PLOAD=100.E4)

IF(KSTEP.EQ.1)THEN！當STEP=1時的荷載大小

P=PLOAD

ELSEIF(KSTEP.EQ.2)THEN！當STEP＝2時的荷載大小

P=COORDS(1)*PLOAD！施加在積分點的荷載P是坐標的函數(shù)

ELSEIF(KSTEP.EQ.3)THEN！當STEP=3時的荷載大小

P=COORDS(1)**2*PLOAD

ELSEIF(KSTEP.EQ.4)THEN！當STEP=4時的荷載大小

P=COORDS(1)**3*PLOAD

ELSEIF(KSTEP.EQ.5)THEN！當STEP=5時的荷載大小

P=COORDS(1)**4*PLOAD

ENDIF

RETURN

END

UMAT子程序具有強大的功能，使用UMAT子程序：

(1)

可以定義材料的本構(gòu)關系，使用ABAQUS材料庫中沒有包含的材料進行計算，擴充程序功能。

(2)幾乎可以用于力學行為分析的任何分析過程，幾乎可以把用戶材料屬性賦予ABAQUS中的任何單元；(3)必須在UMAT中提供材料本構(gòu)模型的雅可比〔Jacobian〕矩陣，即應力增量對應變增量的變化率。

(4)可以和用戶子程序“USDFLD〞聯(lián)合使用，通過“USDFLD〞重新定義單元每一物質(zhì)點上傳遞到UMAT中場變量的數(shù)值。由于主程序與UMAT之間存在數(shù)據(jù)傳遞，甚至共用一些變量，因此必須遵守有關UMAT的書寫格式，UMAT中常用的變量在文件開頭予以定義，通常格式為：

SUBROUTINEUMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,

1RPL,DDSDDT,DRPLDE,DRPLDT,

2STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME,

3NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT,

4CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC)

INCLUDE'ABA_PARAM.INC'

CHARACTER*80CMNAME

DIMENSIONSTRESS(NTENS),STATEV(NSTATV),

1DDSDDE(NTENS,NTENS),DDSDDT(NTENS),DRPLDE(NTENS),

2STRAN(NTENS),DSTRAN(NTENS),TIME(2),PREDEF(1),DPRED(1),

3PROPS(NPROPS),COORDS(3),DROT(3,3),DFGRD0(3,3),DFGRD1(3,3)

usercodingtodefineDDSDDE,STRESS,STATEV,SSE,SPD,SCD

and,ifnecessary,RPL,DDSDDT,DRPLDE,DRPLDT,PNEWDT

RETURN

END

UMAT中的應力矩陣、應變矩陣以及矩陣DDSDDE，DDSDDT，DRPLDE等，都是直接分量存儲在前，剪切分量存儲在后。直接分量有NDI個，剪切分量有NSHR個。各分量之間的順序根據(jù)單元自由度的不同有一些差異，所以編寫UMAT時要考慮到所使用單元的類別。下面對UMAT中用到的一些變量進行說明：

DDSDDE〔NTENS，NTENS〕是一個NTENS維的方陣，稱作雅可比矩陣，，是應力的增量，是應變的增量，DDSDDE〔I，J〕表示增量步結(jié)束時第J個應變分量的改變引起的第I個應力分量的變化。通常雅可比是一個對稱矩陣，除非在“*USERMATERIAL〞語句中參加了“UNSYMM〞參數(shù)。

STRESS〔NTENS〕應力張量矩陣，對應NDI個直接分量和NSHR個剪切分量。在增量步的開始，應力張量矩陣中的數(shù)值通過UMAT和主程序之間的接口傳遞到UMAT中，在增量步的結(jié)束UMAT將對應力張量矩陣更新。對于包含剛體轉(zhuǎn)動的有限應變問題，一個增量步調(diào)用UMAT之前就已經(jīng)對應力張量的進行了剛體轉(zhuǎn)動，因此在UMAT中只需處理應力張量的共旋局部。UMAT中應力張量的度量為柯西〔真實〕應力。

STATEV〔NSTATEV〕用于存儲狀態(tài)變量的矩陣，在增量步開始時將數(shù)值傳遞到UMAT中。也可在子程序USDFLD或UEXPAN中先更新數(shù)據(jù)，然后增量步開始時將更新后的數(shù)據(jù)傳遞到UMAT中。在增量步的結(jié)束必須更新狀態(tài)變量矩陣中的數(shù)據(jù)。和應力張量矩陣不同的是：對于有限應變問題，除了材料本構(gòu)行為引起的數(shù)據(jù)更新以外，狀態(tài)變量矩陣中的任何矢量或者張量都必須通過旋轉(zhuǎn)來考慮材料的剛體運動。狀態(tài)變量矩陣的維數(shù)，等于關鍵字“*DEPVAR〞定義的數(shù)值。狀態(tài)變量矩陣的維數(shù)通過ABAQUS輸入文件中的關鍵字“*DEPVAR〞定義，關鍵字下面數(shù)據(jù)行的數(shù)值即為狀態(tài)變量矩陣的維數(shù)。材料常數(shù)的個數(shù)，等于關鍵字“*USERMATERIAL〞中“CONSTANTS〞常數(shù)設定的值。

PROPS〔NPROPS〕材料常數(shù)矩陣，矩陣中元素的數(shù)值對應于關鍵字“*USERMATERIAL〞下面的數(shù)據(jù)行。

SSE，SPD，SCD分別定義每一增量步的彈性應變能，塑性耗散和蠕變耗散。它們對計算結(jié)果沒有影響，僅僅作為能量輸出。其他變量：

STRAN(NTENS)：應變矩陣；DSTRAN(NTENS)：應變增量矩陣；

DTIME：增量步的時間增量；

NDI：直接應力分量的個數(shù)；

NSHR：剪切應力分量的個數(shù)；

NTENS：總應力分量的個數(shù)，NTENSNDINSHR=+。使用UMAT時需要注意單元的沙漏控制剛度和橫向剪切剛度。通常減縮積分單元的沙漏控制剛度和板、殼、梁單元的橫向剪切剛度是通過材料屬性中的彈性性質(zhì)定義的。這些剛度基于材料初始剪切模量的值，通常在材料定義中通過“*ELASTIC〞選項定義。但是使用UMAT的時候，ABAQUS對程序輸入文件進行預處理的時候得不到剪切模量的數(shù)值。所以這時候用戶必須使用“*HOURGLASSSTIFFNESS〞選項來定義具有沙漏模式的單元的沙漏控制剛度，使用“*TRANSVERSESHEARSTIFFNESS〞選項來定義板、殼、梁單元的橫向剪切剛度。幾個關于子程序的問題及相應解答

Q:本人在用umat作本構(gòu)模型時，

*static,

1,500,0.000001,0.1此時要求的增量步很多，即每次增量要很小，

*static

1,500時，在彈性向塑性過度時，出現(xiàn)錯誤，增量過大，出現(xiàn)尖點.?

A:YOUCANTRYASFOLLOWS:

*STEP,EXTRAPOLATION=NO,INC=2000000

*STATIC

0.001,500.0,0.00001,0.1。

Q:在abaqus中，如果采用umat,利用自己的本構(gòu)，如何讓abaqus明白這種材料的彈塑性應變，也就是說，如何讓程序返回彈性應變與塑性應變，好在output中輸出，我曾想用最笨地方法，在uvarm中定義輸出，利用getvrm獲取材料點的值，但無法獲取增量應力，材料常數(shù)等，研究了幫助中的例子，umatmst3.inp,umatmst3.for,他采用misesJ2流動理論，我在outputhistory顯示他已進入塑性狀態(tài)，但他的PE仍然為0?。?

A:用uvar()勉強成功。

Q:偶在umat中調(diào)用求主應力函數(shù)

CALLSPRINC(STRESS,PS,LSTR,NDI,NSHR)后，存儲主應力得數(shù)組PS中各個主應力排列順序是什么？

PS1>PS2>PS3?

PS1<PS2<PS3?

PS1>PS3>PS2?

A:第二個。個人覺得：

umat實現(xiàn)自己的本構(gòu)沒有固定的方法，對于不同的本構(gòu)有可能必須采用不同的方法。這要靠自己不斷地摸索。有可能一種方法對于簡單加載問題還行，但有可能對于復雜問題并不收斂。最重要一點，就是umat中采用的算法必須consistent.再就是ddsdde必須正確，〔如果采用back_Euler方法等一些算法，ddsdde錯誤有時不影響結(jié)果〔對于簡單加載問題沒有影響，能收斂，〕，但對于復雜問題不收斂。

uptonow，你這個算法對于Mises,hill,J2，J2d等一類的屈服函數(shù)是正確的，但具體的本構(gòu)還要靈活運用，這我也正學習，正在摸索。有時，umat需要很強的有限元根底，并且對采用的本構(gòu)要很熟悉，不要在一顆樹上吊死才好。首先要確認自己的umat沒有錯誤，如果沒有，但就是不收斂〔在不斷減小加載步長的情況下，當然最好對步長不敏感，特別是對于粘彈性，粘塑性，內(nèi)變量一類的材料，有的本構(gòu)取決于背應力的計算〕。那就應該考慮換一種算法。一點體會，請大家探討。

Q:abaqus-uamt的老問題，缺少'ABA_PARAM.INC'文件?

A:在cvf6.5調(diào)試時，顯示缺少'ABA_PARAM.INC'文件！這個沒有任何關系的，這個錯誤將在ABAQUS調(diào)用UMAT的時候自動會找到，僅僅有這個錯誤將沒有任何影響的。也就是說，ABAQUS中調(diào)用的時候，實際并不存在這個錯誤。FT，忘了說一句了，你把或拷到你的程序工作空間后，應該將把或的后綴.sp或.dp去掉，即將ABA_PARAM.INC.dp或改名為ABA_PARAM.INC。呵呵，他的意思是在VisualFortran中調(diào)試其子程序，我覺得這是一個好方法，我當時也是這麼辦的，畢竟在ABAQUS中調(diào)試是非常麻煩的，只有當你的UMAT沒有語法或者明顯的邏輯錯誤，你在ABAQUS中調(diào)試才能事半功倍。

Q：uvarm可以輸出到哪里？.odb可以么？另外那個strav？？就是自己定義用于umat的那個數(shù)組里的數(shù)可以輸出到.odb里么？

A：在umat中，statev是不能用在output中的，statev只是作為一個解的狀態(tài)變量，說來慚愧，我是在umat中定義peeq,(peeq)的求解一般在彈塑性力學書上有〔等效塑性應變〕，用write()寫入一個臨時文件，〔〔切記：這個文件unit號不要與abaqus中的重合，因為他有一些系統(tǒng)默認的文件號，〕〕然后在uvarm中讀取，以uvarm輸出，因為uvarm可以以odb的形式輸出，支持output,field,output,history

Q：在本版看了一個一維固結(jié)的例子，其中含有用戶子程序，如下

SUBROUTINEUFIELD(FIELD,KFIELD,NSECPT,KSTEP,KINC,TIME,NODE,COORDS,TEMP,DTEMP)

INCLUDE'ABA_PARAM.INC'

DIMENSIONFIELD(NSECPT),TIME(2),COORDS(3),TEMP(NSECPT),DTEMP(NSECPT)

KFIELD=1

FIELD(1)=COORDS(2)

RETURN

END我在VF6.5中進行調(diào)試，提示找不到ABA_PARAM.INC請問大俠這如何解決。還有我因為是初次接觸用戶子程序，我查閱了本版所有的相關貼子，都講的不太詳細，我將問題總結(jié)一下，大俠們能不能詳細的講解一下，

1子程序格式(程序后綴是.f;.f90;.for;.obj??)

2CAE中如何調(diào)用，command下如何調(diào)用？

3假設有多個子程序同時存在，如何處理

4我對VF不是很熟，是否可以用VC，C＋＋編寫子程序？

A:假設要在vf中調(diào)試，那么應該根據(jù)需要把SITE文件夾中的ABA_PARAM_DP.INC〔雙精度〕或ABA_PARAM_SP.INC〔單精度〕拷到相應的位置，并改名為ABA_PARAM.INC即可。

1。我試過，.for格是應該是不可以的，至少6.2和6.3版本應該是不行，其他的沒用過，沒有發(fā)言權(quán)。在Abaqus中，運行abaqusj=jobnameuser=username時，默認的用戶子程序后綴名是.for〔.f,.f90應該都不行的，手冊上也有講過〕，只有在username.for文件沒有找到的情況下，才會去搜索username.