以下是關(guān)于用戶自定義材料用戶子程序userpl的翻譯，這是我根據(jù)ansys的理論手冊(cè)以及這個(gè)程序后面的主體部分理解，花費(fèi)一個(gè)星期研究并翻譯的，以后會(huì)把程序部分的翻譯拿上來(lái)。關(guān)于upfs的資料很少，我個(gè)人研究勢(shì)單力薄，希望這方面有興趣的同志一同研究，謝謝大家：

其中有些部分我還不是很明白，希望大家討論并糾正。
源程序在.....ANSYScustomUserintel下面

*deck,userpl parallel user dc
subroutine userpl (elem,intpt,mat,ncomp,kfirst,kfsteq,e,nu,dens,
x prop,d,ktform,timval,timinc,tem,dtem,toffst,flu,dflu,epel,eppl,
x statev,usvr,epeq,plwork,sigepl,sigrat,depeq,dt)
c前面是用到的一些子程序
c *** primary function: allow users to write their own plasticity laws.
c this logic is accessed with tb,user.
c the below demonstration logic is the same as using
c tb,bkin, without adaptive descent (nropt,,,off).
c Other plasticity rules may require internal
c iterations and/or the more general definition of
c plasticity theory, discussed in the Theory
c Manual.
c**函數(shù)的主要功能：如許用戶自己定義塑性準(zhǔn)則，接受命令tb，user中的變
c 量。下面的程序是使用tb，bink的示例。其他的如塑性增長(zhǎng)需要內(nèi)部的迭代
c 過(guò)程進(jìn)行處理。更多的塑性理論方面的定義請(qǐng)參見(jiàn)理論手冊(cè)。
c *** secondary function: demonstrate the use of user-written plasticity laws
c in this routine:
c a. update the nonlinear strain history
c b. compute the material tangent matrix if requested
c **函數(shù)的次要功能：在子程序中，示范了用戶所寫(xiě)的塑性準(zhǔn)則
c a. 根新非線性應(yīng)變的歷史變化。
c b. 如果有需要就計(jì)算材料的切線應(yīng)力應(yīng)變矩陣。
c *** Notice - This file contains ANSYS Confidential information ***
c
c
c *** ansys(r) copyright(c) 2000
c *** ansys, inc.
c
c input arguments:
c 輸入要求
c variable (type,size,intent) description
c 變量（類型，大小，輸入還是輸出） 變量描述
c
c elem (int,sc,in) - element number (label)
c 整形數(shù)，輸入 單元編號(hào)
c intpt (int,sc,in) - element integration point number
c 整形數(shù)，輸入 單元節(jié)點(diǎn)編號(hào)
c mat (int,sc,in) - material reference number
c 整形數(shù)，輸入 材料的相關(guān)編號(hào)
c ncomp (int,sc,in) - no. of stress/strain components (1,4 or 6)
c 1 - x
c 4 - x,y,z,xy
c 6 - x,y,z,xy,yz,xz
c 整形數(shù)，輸入 應(yīng)力和應(yīng)變組成部分的編號(hào)
c kfirst (int,sc,in) - 1 if first time through, 0 otherwise
c (useful for initializing state variables
c to a non-zero value)
c 整形數(shù)，輸入 如果第一次通過(guò)為1，否則為0
c 對(duì)初始化狀態(tài)變量通常是一個(gè)非零值
c kfsteq (int,sc,in) - 1 if first equilibrium iteration of a
c substep, 0 otherwise
c 整形數(shù)，輸入 一個(gè)子步第一次平衡迭代為1，否則為0
c e (dp,sc,in) - average elastic modulus
c 雙精度數(shù)，輸入 彈性模量
c nu (dp,sc,in) - average poisson ratio
c 雙精度數(shù)，輸入 泊松比
c dens (dp,sc,in) - current material density (mass/volume)
c 雙精度數(shù)，輸入 當(dāng)前材料的密度（質(zhì)量/體積）
c prop - linear material property array
c (dp,ar(9),in) (ex,ey,ez, gxy,gyz,gxz, nuxy,nuyz,nuxz)
c (dp,ar(1),in) if ncomp=1 (ex)
c 雙精度數(shù)一維組，輸入 線性材料特性數(shù)組
c d (dp,ar(ncomp,ncomp),in)- elastic stress-strain matrix
c 雙精度數(shù)二維組，輸入 彈性應(yīng)力應(yīng)變矩陣
c ktform (int,sc,in) - request key for tangent matrix formation
c (=1, form tangent .ne.1, do not form)
c 整形數(shù)，輸入 針對(duì)切線矩陣信息所要求的關(guān)鍵字
c （=1就形成切線矩陣，不等于1則不形成）
c
c timval (dp,sc,in) - current time value
c 雙精度數(shù)，輸入 當(dāng)前的時(shí)間值
c timinc (dp,sc,in) - time increment over this substep
c 雙精度數(shù)，輸入 對(duì)每一個(gè)子步過(guò)后的時(shí)間增長(zhǎng)值
c tem (dp,sc,in) - temperature at the end of this substep
c 雙精度數(shù)，輸入 在這個(gè)子步結(jié)束時(shí)候的溫度
c dtem (dp,sc,in) - temperature increment over this substep
c 雙精度數(shù)，輸入 這個(gè)子步結(jié)束的時(shí)候的溫度增長(zhǎng)值
c toffst (dp,sc,in) - temperature offset from absolute zero
c 雙精度數(shù)，輸入 從絕對(duì)0值的溫度偏移量
c flu (dp,sc,in) - fluence at the end of this substep
c 雙精度數(shù)，輸入 在這個(gè)子步結(jié)束時(shí)的影響
c dflu (dp,sc,in) - fluence increment over this substep
c 雙精度數(shù)，輸入 這個(gè)子步結(jié)束后的影響變化
c epel (dp,ar(ncomp),inout)- modified total strain (trial strain)
c epel = eptot - eppl - eptherm - ...
c if a large strain analysis, epel is
c rotation neutralized and is the hencky
c (i.e. log) strain
c 雙精度一維數(shù)組，輸入輸出 也叫trial strain，更改總的應(yīng)變值
c epel＝ 整個(gè)應(yīng)變-上一步塑性應(yīng)變-溫度應(yīng)變.....
c 如果是大變形分析，這里不是很明白：）
c eppl (dp,ar(ncomp),inout)- plastic strain from previous substep
c 雙精度一維數(shù)組，輸入輸出 從前一個(gè)子步等到的塑性應(yīng)變
c statev (dp,ar(ncomp,6),inout)- state variables from previous substep
c 雙精度二維數(shù)組，輸入輸出 從前一個(gè)子步得到的狀態(tài)變量
c usvr (dp,ar(nuval,nintp),inout)- additional state variables from
c previous equilibrium iteration (saved
c if the nsvr command is used)
c 雙精度二維數(shù)組，輸入輸出 從上一步平衡迭代得來(lái)的其他的狀態(tài)變量
c 如果用到這個(gè)命令則存儲(chǔ)
c
c epeq (dp,sc,inout) - effective plastic strain from prev substep
c 雙精度數(shù)，輸入輸出 從上一個(gè)子步得來(lái)的等效塑性應(yīng)變
c plwork (dp,sc,inout) - accumulated plastic work from prev substep
c 雙精度數(shù)，輸入輸出 從上一個(gè)子步得來(lái)的累計(jì)的塑性功
c
c output arguments:
c 輸出要求
c variable (type,size,intent) description
c 變量（類型，大小，輸入還是輸出） 變量描述
c epel (dp,ar(ncomp),inout)- elastic strain
c 雙精度一維數(shù)組，輸入輸出 彈性應(yīng)變
c eppl (dp,ar(ncomp),inout)- updated plastic strain
c 雙精度一維數(shù)組，輸入輸出 更新的塑性應(yīng)變
c
c statev (dp,ar(ncomp,6),inout)- updated state variables
c 雙精度二維數(shù)組，輸入輸出 更新的狀態(tài)變量
c usvr (dp,ar(nuval,nintp),inout)- updated additional state variables
c 雙精度一維數(shù)組，輸入輸出 更新的額外狀態(tài)變量
c
c epeq (dp,sc,inout) - updated effective plastic strain
c 雙精度數(shù)，輸入輸出 更新的等效塑性應(yīng)變
c plwork (dp,sc,inout) - updated accumulated plastic work
c 雙精度數(shù)，輸入輸出 根新的累計(jì)塑性功
c
c sigepl (dp,sc,out) - stress value on stress-strain curve at epeq
c 雙精度數(shù)，輸出 在更新后的等效塑性應(yīng)變后應(yīng)力應(yīng)變曲線上的應(yīng)力值
c sigrat (dp,sc,out) - ratio of trial stress to yield stress
c 雙精度數(shù)，輸出 trial stress到屈服面上的比率
c depeq (dp,sc,out) - increment in plastic strain (equivalent)
c (used for auto time stepping - time step
c is reduced if it exceeds .05)
c 雙精度數(shù)，輸出 等效塑性應(yīng)變的增長(zhǎng)
c （利用自動(dòng)時(shí)間步－如果它超過(guò)0.05時(shí)間步減少）
c
c dt (dp,ar(ncomp,ncomp),out)- material tangent modulus
c 雙精度二維數(shù)組，輸出 材料的切線模量
c
c fortran parameters (to be defined by the user):
c fortran的參數(shù)（由用戶自己定義的）
c variable (type) description
c 變量（類型） 描述
c numinp (int) - number of data items in the user-defined
c data table (tbdat commands)
c 整形 用戶通過(guò)tbdat命令定義數(shù)據(jù)的編號(hào)
c nuval (int) - number of additional state variables per
c integration point
c 整形 每個(gè)積分點(diǎn)額外狀態(tài)變量的編號(hào)
c nintp (int) - maximum number of integration points of
c an element to be used with this routine
c (14 is the maximum)
c 整形 這個(gè)程序使用的一個(gè)單元積分點(diǎn)的最大編號(hào)
c note: nuval x nintp = nstv(on nsvr command); cannot exceed 840!
c
c internal variables:
c 內(nèi)部變量
c variable (type,size) description
c 變量（類型，大小） 描述
c b (dp,ar(6,6)) - 2nd derivative of the yield function
c 雙精度二維數(shù)組 屈服函數(shù)的二階倒數(shù)
c c (dp,ar(6,12)) - part of deff = (i + dlamb*d*b)-1*d; b = (m - n*nt)/sigy
c 雙精度二維數(shù)組 是deff的一部分（這里有點(diǎn)模?/a>??
c con (dp,sc) - temporary variable
c 雙精度數(shù) 臨時(shí)變量
c deppl (dp,ar(6)) - plastic strain increment
c 雙精度一維數(shù)組 塑性應(yīng)變?cè)鲩L(zhǎng)
c dfds (dp,ar(6)) - derivative of the yield function (normal)
c 雙精度一維數(shù)組 屈服函數(shù)的偏倒數(shù)
c dlamb (dp,sc) - plastic multiplier
c 雙精度數(shù) 塑性因子λ。
c ep (dp,ar(6)) - shifted strain
c 雙精度一維數(shù)組 轉(zhuǎn)換應(yīng)變
c epshfo (dp,ar(6)) - initial shift strain (center of the yield surf)
c 雙精度一維數(shù)組 初始轉(zhuǎn)換應(yīng)變（屈服面中心）
c epshft (dp,ar(6)) - shift strain (center of the yield surface)
c 雙精度一維數(shù)組 轉(zhuǎn)換應(yīng)變（屈服面中心）
c et (dp,sc) - tangent modulus (stress/total strain)
c 雙精度數(shù) 切線模量（應(yīng)力/總應(yīng)變）－－不知道這個(gè)為什么呀：）
c h (dp,sc) - plastic tangent modulus (stress/plastic strain)
c 雙精度數(shù) 塑性切線模量（應(yīng)力/塑性應(yīng)變）
c n2 (int,sc) - ncomp squared, ncomp*ncomp
c 整形數(shù) ncomp的平方
c seqtr (dp,sc) - equivalent (von mises) trial stress
c 整形數(shù) 等效trialstress（mises屈服準(zhǔn)則）
c sigtr (dp,ar(6)) - trial stress
c 雙精度一維數(shù)組 trial stress
c sigy (dp,sc) - yield stress
c 雙精度數(shù) 屈服應(yīng)力
c vect (dp,ar(6)) - temporary vector
c 雙精度一維數(shù)組 臨時(shí)的向量or標(biāo)量

這里是變量的說(shuō)明，有了這些大家可以很快理解子程序，所剩下的就是對(duì)理論和算法的討論了。希望大家積極和我交流意見(jiàn)。：）