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

其中有些部分我還不是很明白，希望大家討論并糾正。
源程序在.....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**函數的主要功能：如許用戶自己定義塑性準則，接受命令tb，user中的變
c 量。下面的程序是使用tb，bink的示例。其他的如塑性增長需要內部的迭代
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 **函數的次要功能：在子程序中，示范了用戶所寫的塑性準則
c a. 根新非線性應變的歷史變化。
c b. 如果有需要就計算材料的切線應力應變矩陣。
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 整形數，輸入 單元編號
c intpt (int,sc,in) - element integration point number
c 整形數，輸入 單元節點編號
c mat (int,sc,in) - material reference number
c 整形數，輸入 材料的相關編號
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 整形數，輸入 應力和應變組成部分的編號
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 整形數，輸入 如果第一次通過為1，否則為0
c 對初始化狀態變量通常是一個非零值
c kfsteq (int,sc,in) - 1 if first equilibrium iteration of a
c substep, 0 otherwise
c 整形數，輸入 一個子步第一次平衡迭代為1，否則為0
c e (dp,sc,in) - average elastic modulus
c 雙精度數，輸入 彈性模量
c nu (dp,sc,in) - average poisson ratio
c 雙精度數，輸入 泊松比
c dens (dp,sc,in) - current material density (mass/volume)
c 雙精度數，輸入 當前材料的密度（質量/體積）
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 雙精度數一維組，輸入 線性材料特性數組
c d (dp,ar(ncomp,ncomp),in)- elastic stress-strain matrix
c 雙精度數二維組，輸入 彈性應力應變矩陣
c ktform (int,sc,in) - request key for tangent matrix formation
c (=1, form tangent .ne.1, do not form)
c 整形數，輸入 針對切線矩陣信息所要求的關鍵字
c （=1就形成切線矩陣，不等于1則不形成）
c
c timval (dp,sc,in) - current time value
c 雙精度數，輸入 當前的時間值
c timinc (dp,sc,in) - time increment over this substep
c 雙精度數，輸入 對每一個子步過后的時間增長值
c tem (dp,sc,in) - temperature at the end of this substep
c 雙精度數，輸入 在這個子步結束時候的溫度
c dtem (dp,sc,in) - temperature increment over this substep
c 雙精度數，輸入 這個子步結束的時候的溫度增長值
c toffst (dp,sc,in) - temperature offset from absolute zero
c 雙精度數，輸入 從絕對0值的溫度偏移量
c flu (dp,sc,in) - fluence at the end of this substep
c 雙精度數，輸入 在這個子步結束時的影響
c dflu (dp,sc,in) - fluence increment over this substep
c 雙精度數，輸入 這個子步結束后的影響變化
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 雙精度一維數組，輸入輸出 也叫trial strain，更改總的應變值
c epel＝ 整個應變-上一步塑性應變-溫度應變.....
c 如果是大變形分析，這里不是很明白：）
c eppl (dp,ar(ncomp),inout)- plastic strain from previous substep
c 雙精度一維數組，輸入輸出 從前一個子步等到的塑性應變
c statev (dp,ar(ncomp,6),inout)- state variables from previous substep
c 雙精度二維數組，輸入輸出 從前一個子步得到的狀態變量
c usvr (dp,ar(nuval,nintp),inout)- additional state variables from
c previous equilibrium iteration (saved
c if the nsvr command is used)
c 雙精度二維數組，輸入輸出 從上一步平衡迭代得來的其他的狀態變量
c 如果用到這個命令則存儲
c
c epeq (dp,sc,inout) - effective plastic strain from prev substep
c 雙精度數，輸入輸出 從上一個子步得來的等效塑性應變
c plwork (dp,sc,inout) - accumulated plastic work from prev substep
c 雙精度數，輸入輸出 從上一個子步得來的累計的塑性功
c
c output arguments:
c 輸出要求
c variable (type,size,intent) description
c 變量（類型，大小，輸入還是輸出） 變量描述
c epel (dp,ar(ncomp),inout)- elastic strain
c 雙精度一維數組，輸入輸出 彈性應變
c eppl (dp,ar(ncomp),inout)- updated plastic strain
c 雙精度一維數組，輸入輸出 更新的塑性應變
c
c statev (dp,ar(ncomp,6),inout)- updated state variables
c 雙精度二維數組，輸入輸出 更新的狀態變量
c usvr (dp,ar(nuval,nintp),inout)- updated additional state variables
c 雙精度一維數組，輸入輸出 更新的額外狀態變量
c
c epeq (dp,sc,inout) - updated effective plastic strain
c 雙精度數，輸入輸出 更新的等效塑性應變
c plwork (dp,sc,inout) - updated accumulated plastic work
c 雙精度數，輸入輸出 根新的累計塑性功
c
c sigepl (dp,sc,out) - stress value on stress-strain curve at epeq
c 雙精度數，輸出 在更新后的等效塑性應變后應力應變曲線上的應力值
c sigrat (dp,sc,out) - ratio of trial stress to yield stress
c 雙精度數，輸出 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 雙精度數，輸出 等效塑性應變的增長
c （利用自動時間步－如果它超過0.05時間步減少）
c
c dt (dp,ar(ncomp,ncomp),out)- material tangent modulus
c 雙精度二維數組，輸出 材料的切線模量
c
c fortran parameters (to be defined by the user):
c fortran的參數（由用戶自己定義的）
c variable (type) description
c 變量（類型） 描述
c numinp (int) - number of data items in the user-defined
c data table (tbdat commands)
c 整形 用戶通過tbdat命令定義數據的編號
c nuval (int) - number of additional state variables per
c integration point
c 整形 每個積分點額外狀態變量的編號
c nintp (int) - maximum number of integration points of
c an element to be used with this routine
c (14 is the maximum)
c 整形 這個程序使用的一個單元積分點的最大編號
c note: nuval x nintp = nstv(on nsvr command); cannot exceed 840!
c
c internal variables:
c 內部變量
c variable (type,size) description
c 變量（類型，大小） 描述
c b (dp,ar(6,6)) - 2nd derivative of the yield function
c 雙精度二維數組 屈服函數的二階倒數
c c (dp,ar(6,12)) - part of deff = (i + dlamb*d*b)-1*d; b = (m - n*nt)/sigy
c 雙精度二維數組 是deff的一部分（這里有點模?/a>??
c con (dp,sc) - temporary variable
c 雙精度數 臨時變量
c deppl (dp,ar(6)) - plastic strain increment
c 雙精度一維數組 塑性應變增長
c dfds (dp,ar(6)) - derivative of the yield function (normal)
c 雙精度一維數組 屈服函數的偏倒數
c dlamb (dp,sc) - plastic multiplier
c 雙精度數 塑性因子λ。
c ep (dp,ar(6)) - shifted strain
c 雙精度一維數組 轉換應變
c epshfo (dp,ar(6)) - initial shift strain (center of the yield surf)
c 雙精度一維數組 初始轉換應變（屈服面中心）
c epshft (dp,ar(6)) - shift strain (center of the yield surface)
c 雙精度一維數組 轉換應變（屈服面中心）
c et (dp,sc) - tangent modulus (stress/total strain)
c 雙精度數 切線模量（應力/總應變）－－不知道這個為什么呀：）
c h (dp,sc) - plastic tangent modulus (stress/plastic strain)
c 雙精度數 塑性切線模量（應力/塑性應變）
c n2 (int,sc) - ncomp squared, ncomp*ncomp
c 整形數 ncomp的平方
c seqtr (dp,sc) - equivalent (von mises) trial stress
c 整形數 等效trialstress（mises屈服準則）
c sigtr (dp,ar(6)) - trial stress
c 雙精度一維數組 trial stress
c sigy (dp,sc) - yield stress
c 雙精度數 屈服應力
c vect (dp,ar(6)) - temporary vector
c 雙精度一維數組 臨時的向量or標量

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