VUMAT中英文對(duì)照

1、25.3.4 VUMAT25.3.4 VUMATUser subroutine to define material behavior.定義材料本構(gòu)用戶子程序Product: ABAQUS/Explicit  Warning:  The use of this user subroutine generally requires considerable expertise. You are cautioned that the implementation of any realistic constitutive model require

2、s extensive development and testing. Initial testing on a single-element model with prescribed traction loading is strongly recommended.注意：用戶子程序的使用通常需要一定的專長。用戶需要知道執(zhí)行任何實(shí)際的本構(gòu)模型需要大量的試驗(yàn)數(shù)據(jù)。強(qiáng)烈建議用戶對(duì)用戶子程序進(jìn)行在指定拉力作用下單個(gè)單元的驗(yàn)證測試。The component ordering of the symmetric and nonsymmetric tensors for the three-dimen

3、sional case using C3D8R elements is different from the ordering specified in , and the ordering used in ABAQUS/Standard.C3D8R單元三維軸對(duì)稱及非軸對(duì)稱張量成分順序與及ABAQUS/Standard中指定的順序不同。References· “User-defined mechanical m· *USER MATERIALOverviewUser subroutine VUMAT: 用戶子程序VUMAT· is used to define t

4、he mechanical constitutive behavior of a material;· 用來定義材料的力學(xué)本構(gòu)關(guān)系· will be called for blocks of material calculation points for which the material is defined in a user subroutine ();· 可以被用戶子程序定義材料計(jì)算點(diǎn)調(diào)用· can use and update solution-dependent state variables;· 可以使用和更新結(jié)果依賴狀態(tài)變量&

5、#183; can use any field variables that are passed in;· 可以使用傳入的任何場變量· is described further in ; and· 在中詳細(xì)論述；· cannot be used in an adiabatic analysis.· 可以被用于絕熱分析Component ordering in tensors 張量組成順序The component ordering depends upon whether the tensor is symmetric or nonsymme

6、tric.張量組成順序取決于其是否為對(duì)稱或非對(duì)稱張量。Symmetric tensors 對(duì)稱張量For symmetric tensors such as the stress and strain tensors, there are ndir+nshr components, and the component order is given as a natural permutation of the indices of the tensor. The direct components are first and then the indirect components, begi

7、nning with the 12-component. For example, a stress tensor contains ndir direct stress components and nshr shear stress components, which are passed in as 對(duì)于如同應(yīng)力及應(yīng)變張量等的對(duì)稱張量，含有ndir+nshr分量，分量的序號(hào)按張量索引號(hào)的自然排序給出。首先是直接分量，然后是從12分量開始的間接分量。例如，包含ndir正應(yīng)力分量及nshr的剪應(yīng)力張量的應(yīng)力張量被按照下面的順序傳入Component2-D Case3-D Case12345&

8、#160;6 The shear strain components in user subroutine VUMAT are stored as tensor components and not as engineering components; this is different from user subroutine UMAT in ABAQUS/Standard, which uses engineering components. Nonsymmetric tensors 非對(duì)稱張量For nonsymmetric tensors there are ndir+2*n

9、shr components, and the component order is given as a natural permutation of the indices of the tensor. The direct components are first and then the indirect components, beginning with the 12-component. For example, the deformation gradient is passed as 對(duì)于非對(duì)稱張量含有ndir+2*nshr分量，分量的順序按照張量索引號(hào)的自然排序給出。首先是

10、直接分量，其次是從12分量開始的間接分量。例如，位移梯度按照下面的順序傳遞Component2-D Case3-D Case123456 7 8 9 Initial calculations and checks 最初計(jì)算和檢查In the data check phase of the analysis ABAQUS/Explicit calls user subroutine VUMAT with a set of fictitious strains and a totalTime and stepTime both equal to 0.0. T

11、his is done as a check on your constitutive relation and to calculate the equivalent initial material properties, based upon which the initial elastic wave speeds are computed.在ABAQUS/Explicit調(diào)用用戶子程序VUMAT分析的數(shù)據(jù)檢查階段,小應(yīng)變、總時(shí)間及時(shí)間步都為0。這作為對(duì)用戶本構(gòu)關(guān)系的一個(gè)檢查，基于計(jì)算得到的初始材料波速來計(jì)算等效初始材料屬性。Defining local orientations 定義

12、局部方向All stresses, strains, stretches, and state variables are in the orientation of the local material axes. These local material axes form a basis system in which stress and strain components are stored. This represents a corotational coordinate system in which the basis system rotates with the mat

13、erial. If a user-specified coordinate system () is used, it defines the local material axes in the undeformed configuration.所有的應(yīng)力、應(yīng)變、延伸及狀態(tài)變量均按局部材料軸的方向。這些局部材料軸形成一個(gè)應(yīng)力與應(yīng)變分量存儲(chǔ)的基本系統(tǒng)。即這個(gè)基本系是隨著材料聯(lián)合轉(zhuǎn)動(dòng)的坐標(biāo)系。如果使用用戶指定坐標(biāo)系，則它在無變形結(jié)構(gòu)中定義局部材料軸。Special considerations for various element types不同單元類型的特殊考慮The use of user

14、 subroutine VUMAT requires special consideration for various element types.用戶子程序VUMAT的使用需要對(duì)不同的單元類型進(jìn)行特殊的考慮。Shell and plane stress elements殼及平面應(yīng)力單元You must define the stresses and internal state variables. In the case of shell or plane stress elements, you must define strainInc(*,3), the thickness str

15、ain increment. The internal energies can be defined if desired. If they are not defined, the energy balance provided by ABAQUS/Explicit will not be meaningful.用戶必須定義應(yīng)力和初始狀態(tài)變量。在殼或平面應(yīng)力單元的情況下，用戶必須定義應(yīng)變包括（*,3），厚度應(yīng)變?cè)隽俊Ｈ绻枰脑掃€需要定義初始能量。如果沒有定義，那么ABAQUS/Explicit提供的能量平衡將沒有意義。Shell elements殼單元When VUMAT is used

16、 to define the material response of shell elements, ABAQUS/Explicit cannot calculate a default value for the transverse shear stiffness of the element. Hence, you must define the element's transverse shear stiffness. See , for guidelines on choosing this stiffness.當(dāng)使用VUMAT定義殼單元的材料響應(yīng)時(shí)，ABAQUS/Expl

17、icit不能計(jì)算單元的缺省橫向剪切剛度。因此，用戶需要定義單元的橫向剪切剛度。關(guān)于選擇橫行剪切剛度的詳細(xì)資料請(qǐng)參考Beam elements梁單元For beam elements the stretch tensor and the deformation gradient tensor are not available. For beams in space you must define the thickness strains, strainInc(*,2) and strainInc(*,3). strainInc(*,4) is the shear strain associa

18、ted with twist. Thickness stresses, stressNew(*,2) and stressNew(*,3), are assumed to be zero and any values you assign are ignored.對(duì)于梁單元不能使用拉伸張量及位移梯度張量。對(duì)于空間梁，用戶必須定義厚度應(yīng)變、應(yīng)變?cè)隽浚?，2）及應(yīng)變?cè)隽浚?，3）。應(yīng)變?cè)隽浚?，4）時(shí)與扭曲有關(guān)的剪應(yīng)變。厚度應(yīng)力，（*，2）及（*，3）假定為0，并且用戶分配的任何相關(guān)張量都被忽略。Deformation gradient 位移梯度The polar decomposition of

19、 the deformation gradient is written as , where and are the right and left symmetric stretch tensors, respectively. The constitutive model is defined in a corotational coordinate system in which the basis system rotates with the material. All stress and strain tensor quantities are defined with resp

20、ect to the corotational basis system. The right stretch tensor, , is used. The relative spin tensor represents the spin (the antisymmetric part of the velocity gradient) defined with respect to the corotational basis system.位移梯度寫成，其中及分別為右邊及左邊的對(duì)稱拉伸張量。本構(gòu)模型定義為聯(lián)合旋轉(zhuǎn)坐標(biāo)系，在該坐標(biāo)系中基系隨著材料轉(zhuǎn)動(dòng)。所有的應(yīng)力和應(yīng)變張量值按照聯(lián)合旋轉(zhuǎn)坐標(biāo)系

21、定義。使用右邊的拉伸位移。相應(yīng)的旋轉(zhuǎn)張量代表與聯(lián)合旋轉(zhuǎn)基系相應(yīng)的轉(zhuǎn)動(dòng)。Special considerations for hyperelasticity超彈性的特殊考慮Hyperelastic constitutive models in VUMAT should be defined in a corotational coordinate system in which the basis system rotates with the material. This is most effectively accomplished by formulating the hyperelas

22、tic constitutive model in terms of the stretch tensor, , instead of in terms of the deformation gradient, . Using the deformation gradient can present some difficulties because the deformation gradient includes the rotation tensor and the resulting stresses would need to be rotated back to the corot

23、ational basis. 在VUMAT中的超彈性本構(gòu)模型可以被定義在聯(lián)合選擇坐標(biāo)系中。這可以通過用拉伸張量表示的超彈性本構(gòu)模型很好的實(shí)現(xiàn)，而不是使用位移梯度來表示。使用位移梯度可能會(huì)帶來一些困難，因?yàn)槲灰铺荻劝ㄐD(zhuǎn)張量并且導(dǎo)致應(yīng)力需要選擇返回到聯(lián)合旋轉(zhuǎn)基系。Objective stress rates目標(biāo)應(yīng)力率The Green-Naghdi stress rate is used when the mechanical behavior of the material is defined using user subroutine VUMAT. The stress rate obt

24、ained with user subroutine VUMAT may differ from that obtained with a built-in ABAQUS material model. For example, most material models used with solid (continuum) elements in ABAQUS/Explicit employ the Jaumann stress rate. This difference in the formulation will cause significant differences in the

25、 results only if finite rotation of a material point is accompanied by finite shear. For a discussion of the objective stress rates used in ABAQUS, see “Stress rates,” Section 1.5.3 of the ABAQUS Theory Manual.在用戶子程序VUMAT中使用Green-Naghdi應(yīng)力率來定義材料的力學(xué)本構(gòu)關(guān)系。通過用戶子程序VUMAT獲得的應(yīng)力率可能會(huì)與在ABAQUS建立的材料模型獲得的應(yīng)力率有所不同。例

26、如，在ABAQUS/Explicit中大多數(shù)實(shí)體（連續(xù)）單元材料模型使用Jaumann應(yīng)力率。只要材料點(diǎn)的有限旋轉(zhuǎn)伴隨有限剪切，這種表達(dá)方式的不同將導(dǎo)致計(jì)算結(jié)果的明顯差異。關(guān)于ABAQUS中使用的目標(biāo)應(yīng)力率的詳細(xì)討論參考“Stress rates,” Section 1.5.3 of the ABAQUS Theory Manual.Material point deletion 材料點(diǎn)刪除Material points that satisfy a user-defined failure criterion can be deleted from the model (see ). You

27、 must specify the state variable number controlling the element deletion flag when you allocate space for the solution-dependent state variables, as explained in . The deletion state variable should be set to a value of one or zero in VUMAT. A value of one indicates that the material point is active

28、, while a value of zero indicates that ABAQUS/Explicit should delete the material point from the model by setting the stresses to zero. The structure of the block of material points passed to user subroutine VUMAT remains unchanged during the analysis; deleted material points are not removed from th

29、e block. ABAQUS/Explicit will pass zero stresses and strain increments for all deleted material points. Once a material point has been flagged as deleted, it cannot be reactivated.滿足用戶定義的破壞準(zhǔn)則的材料點(diǎn)可以被從模型中刪除（參考)。當(dāng)用戶給結(jié)果依賴狀態(tài)變量分配空間時(shí)，用戶需要指定控制單元?jiǎng)h除標(biāo)示的狀態(tài)變量號(hào)，在中進(jìn)行詳細(xì)說明。在VUMAT中刪除狀態(tài)變量可以被賦予1或者0。1表示材料點(diǎn)時(shí)激活的，0表示ABAQUS

30、/Explicit將通過設(shè)定應(yīng)力為0刪除材料點(diǎn)。在分析過程中傳遞給用戶子程序VUMAT的材料點(diǎn)結(jié)構(gòu)保持不變；刪除的材料點(diǎn)沒有從塊中移走。ABAQUS/Explicit將傳遞0應(yīng)力及應(yīng)變給所有刪除的材料點(diǎn)。一旦一個(gè)材料點(diǎn)被標(biāo)示為刪除，該材料點(diǎn)將不能夠被再次激活。User subroutine interface 用戶子程序 subroutine vumat(C Read only (unmodifiable)variables - 1 nblock, ndir, nshr, nstatev, nfieldv, nprops, lanneal, 2 stepTime, totalTime, dt,

31、 cmname, coordMp, charLength, 3 props, density, strainInc, relSpinInc, 4 tempOld, stretchOld, defgradOld, fieldOld, 5 stressOld, stateOld, enerInternOld, enerInelasOld, 6 tempNew, stretchNew, defgradNew, fieldNew,C Write only (modifiable) variables - 7 stressNew, stateNew, enerInternNew, enerInelasN

32、ew )C include 'vaba_param.inc'C dimension props(nprops), density(nblock), coordMp(nblock,*), 1 charLength(nblock), strainInc(nblock,ndir+nshr), 2 relSpinInc(nblock,nshr), tempOld(nblock), 3 stretchOld(nblock,ndir+nshr), 4 defgradOld(nblock,ndir+nshr+nshr), 5 fieldOld(nblock,nfieldv), stressO

33、ld(nblock,ndir+nshr), 6 stateOld(nblock,nstatev), enerInternOld(nblock), 7 enerInelasOld(nblock), tempNew(nblock), 8 stretchNew(nblock,ndir+nshr), 8 defgradNew(nblock,ndir+nshr+nshr), 9 fieldNew(nblock,nfieldv), 1 stressNew(nblock,ndir+nshr), stateNew(nblock,nstatev), 2 enerInternNew(nblock), enerIn

34、elasNew(nblock),C character*80 cmnameC do 100 km = 1,nblock user coding 100 continue return endVariables to be defined 被定義的變量stressNew (nblock, ndir+nshr) Stress tensor at each material point at the end of the increment.在增量結(jié)束時(shí)每個(gè)材料點(diǎn)的應(yīng)力張量。stateNew (nblock, nstatev) State variables at each material poi

35、nt at the end of the increment. You define the size of this array by allocating space for it (see , for more information).增量結(jié)束時(shí)每個(gè)材料點(diǎn)的狀態(tài)變量。用戶通過分配空間來定義該矩陣的大?。ǜ嗟馁Y料參考）。Variables that can be updated 可以更新的變量enerInternNew (nblock) Internal energy per unit mass at each material point at the end of the incr

36、ement.增量結(jié)束時(shí)每個(gè)材料點(diǎn)單位質(zhì)量的內(nèi)能。enerInelasNew (nblock) Dissipated inelastic energy per unit mass at each material point at the end of the increment.增量結(jié)束時(shí)每個(gè)材料點(diǎn)單位質(zhì)量的消散的無彈性能。Variables passed in for informationnblock Number of material points to be processed in this call to VUMAT.調(diào)用VUMAT的材料點(diǎn)號(hào)。ndir Number of dir

37、ect components in a symmetric tensor.對(duì)成張量的直接分量號(hào)。nshr Number of indirect components in a symmetric tensor.對(duì)稱張量的間接分量號(hào)。nstatev Number of user-defined state variables that are associated with this material type (you define this as described in ).與材料類型相關(guān)的用戶定義狀態(tài)變量號(hào)。nfieldv Number of user-defined external

38、field variables.用戶定義外部場變量號(hào)。nprops User-specified number of user-defined material properties.用戶定義材料屬性的用戶指定號(hào)。lanneal Flag indicating whether the routine is being called during an annealing process. lanneal=0 indicates that the routine is being called during a normal mechanics increment. lanneal=1 indi

39、cates that this is an annealing process and you should re-initialize the internal state variables, stateNew, if necessary. ABAQUS/Explicit will automatically set the stresses, stretches, and state to a value of zero during the annealing process.在退火處理過程中標(biāo)示程序是否被調(diào)用。Laneal=0表明程序在正常力學(xué)增量過程中被調(diào)用。Laneal1表明這是

40、一個(gè)退火過程，并且如果需要的話用戶需要重新初始化內(nèi)部狀態(tài)變量stateNew。ABAQUS/Explicit將自動(dòng)設(shè)置應(yīng)力，延展性及在退火過程中0值狀態(tài)。stepTime Value of time since the step began.從時(shí)間步開始時(shí)的的時(shí)間totalTime Value of total time. The time at the beginning of the step is given by totalTime - stepTime.總時(shí)間值。時(shí)間步開始時(shí)的時(shí)間定義為totalTime - stepTimedt Time increment size.時(shí)間增量大小

41、。cmname User-specified material name, left justified. It is passed in as an upper-case character string. Some internal material models are given names starting with the “ABQ_” character string. To avoid conflict, you should not use “ABQ_” as the leading string for cmname.用戶指定材料名。按照大寫字母傳入。有些內(nèi)部材料本構(gòu)以AB

42、Q_字母開頭賦名。為了避免沖突，用戶不能使用ABQ_作為cmname的開頭字母。coordMp(nblock,*) Material point coordinates. It is the midplane material point for shell elements and the centroid for beam elements.材料點(diǎn)坐標(biāo)。對(duì)于殼單元為中平面材料點(diǎn)，對(duì)于梁單元為質(zhì)心。charLength(nblock) Characteristic element length. This is a typical length of a line across an ele

43、ment. For beams and trusses, it is a characteristic length along the element axis. For membranes and shells, it is a characteristic length in the reference surface. For axisymmetric elements, it is a characteristic length in the plane only. For cohesive elements it is equal to the constitutive thick

44、ness.特征單元長度。它是一個(gè)穿過單元線的特殊長度。對(duì)于梁及桁架，它表示沿單元軸的特征長度。對(duì)于膜及殼，它表示參考面上的特征長度。對(duì)于對(duì)稱單元，它僅僅表示平面上的特征長度。對(duì)于粘聚單元它等于結(jié)構(gòu)厚度。props(nprops) User-supplied material properties.用戶指定材料屬性density(nblock) Current density at the material points in the midstep configuration. This value may be inaccurate in problems where the volumet

45、ric strain increment is very small. If an accurate value of the density is required in such cases, the analysis should be run in double precision. This value of the density is not affected by mass scaling.材料點(diǎn)的當(dāng)前密度。在體積應(yīng)變?cè)隽糠浅Ｐ〉那闆r下，該值可能出現(xiàn)精確性問題。在這種情況下如果需要精確的密度值，分析需要在雙精度下運(yùn)行。密度值不受質(zhì)量縮放比例的影響strainInc (nbloc

46、k, ndir+nshr) Strain increment tensor at each material point.每個(gè)材料點(diǎn)的應(yīng)變?cè)隽繌埩?。relSpinInc (nblock, nshr) Incremental relative rotation vector at each material point defined in the corotational system. Defined as , where is the antisymmetric part of the velocity gradient, , and . Stored in 3-D as and in 2

47、-D as .定義在聯(lián)合旋轉(zhuǎn)系中的每個(gè)材料點(diǎn)的增量相關(guān)選擇矢量。定義為，其中為速度梯度，的反對(duì)稱部分，在三維條件下存儲(chǔ)為，在二維條件下存儲(chǔ)為。tempOld(nblock) Temperatures at each material point at the beginning of the increment.增量開始時(shí)每個(gè)材料點(diǎn)的溫度。stretchOld (nblock, ndir+nshr) Stretch tensor, , at each material point at the beginning of the increment defined from the polar

48、decomposition of the deformation gradient by .增量開始時(shí)每個(gè)材料點(diǎn)的，拉伸張量。defgradOld (nblock,ndir+2*nshr) Deformation gradient tensor at each material point at the beginning of the increment. Stored in 3-D as (, , , , , , , , ) and in 2-D as (, , , , ).增量開始時(shí)每個(gè)材料點(diǎn)的位移梯度張量。在三維狀態(tài)下存儲(chǔ)為(, , , , , , , , ) 及在二維條件下存儲(chǔ)為(,

49、 , , , ).fieldOld (nblock, nfieldv) Values of the user-defined field variables at each material point at the beginning of the increment.增量開始時(shí)每個(gè)材料點(diǎn)的用戶定義場變量值。stressOld (nblock, ndir+nshr) Stress tensor at each material point at the beginning of the increment.增量開始時(shí)每個(gè)材料點(diǎn)的應(yīng)力張量。stateOld (nblock, nstatev)

50、State variables at each material point at the beginning of the increment.增量開始時(shí)每個(gè)材料點(diǎn)的狀態(tài)變量。enerInternOld (nblock) Internal energy per unit mass at each material point at the beginning of the increment.增量開始時(shí)每個(gè)材料點(diǎn)的單位質(zhì)量的內(nèi)能enerInelasOld (nblock) Dissipated inelastic energy per unit mass at each material p

51、oint at the beginning of the increment.增量開始時(shí)每個(gè)材料點(diǎn)的單位質(zhì)量消散的無彈性能。tempNew(nblock) Temperatures at each material point at the end of the increment.增量開始時(shí)每個(gè)材料點(diǎn)的溫度。stretchNew (nblock, ndir+nshr) Stretch tensor, , at each material point at the end of the increment defined from the polar decomposition of the

52、deformation gradient by .defgradNew (nblock,ndir+2*nshr) Deformation gradient tensor at each material point at the end of the increment. Stored in 3-D as (, , , , , , , , ) and in 2-D as (, , , , ).增量結(jié)束時(shí)每個(gè)材料點(diǎn)的位移梯度張量。三維狀態(tài)下存儲(chǔ)為(, , , , , , , , ) ，二維狀態(tài)下存儲(chǔ)為 (, , , , )。fieldNew (nblock, nfieldv) Values of

53、 the user-defined field variables at each material point at the end of the increment.增量結(jié)束時(shí)每個(gè)材料點(diǎn)的用戶指定場變量值。Example: Using more than one user-defined material model例子：使用多個(gè)用戶定義材料模型To use more than one user-defined material model, the variable cmname can be tested for different material names inside

54、 user subroutine VUMAT, as illustrated below:要使用多種用戶材料模型，用戶子程序VUMAT中使用cmname值代表不同的材料名，如下所述：if (cmname(1:4) .eq. 'MAT1') then call VUMAT_MAT1(argument_list)else if (cmname(1:4) .eq. 'MAT2') then call VUMAT_MAT2(argument_list)end ifVUMAT_MAT1 and VUMAT_MAT2 are the actual user material

55、 subroutines containing the constitutive material models for each material MAT1 and MAT2, respectively. Subroutine VUMAT merely acts as a directory here. The argument list can be the same as that used in subroutine VUMAT. The material names must be in upper case since cmname is passed in as an upper

56、-case character string. VUMAT_MAT1 and VUMAT_MAT2分別為實(shí)際的包含MAT1及MAT2材料本構(gòu)的用戶材料子程序。子程序VUMAT在此僅僅作為一個(gè)地址。Example: Elastic/plastic material with kinematic hardening例子：帶有運(yùn)動(dòng)硬化的彈性/塑性材料As a simple example of the coding of subroutine VUMAT, consider the generalized plane strain case for an elastic/plastic m

57、aterial with kinematic hardening. The basic assumptions and definitions of the model are as follows.作為子程序VUMAT的簡單例子代碼，考慮通用運(yùn)動(dòng)硬化彈性/塑性材料的平面應(yīng)變問題。模型的基本假設(shè)及定義如下。Let be the current value of the stress and define to be the deviatoric part of the stress. The center of the yield surface in deviatoric stress sp

58、ace is given by the tensor , which has initial values of zero. The stress difference, , is the stress measured from the center of the yield surface and is given by 設(shè)定為應(yīng)力的當(dāng)前值，并定義應(yīng)力的偏量部分。在偏應(yīng)力空間屈服面中心由張量給定，其初始值為0。應(yīng)力導(dǎo)數(shù)為從屈服面中學(xué)測量的應(yīng)力，由給出。The von Mises yield surface is defined as Mises屈服面定義為where is the uniaxia