VUMAT-中英文對照

1、【精品文檔】如有侵權(quán)，請聯(lián)系網(wǎng)站刪除，僅供學(xué)習(xí)與交流VUMAT-中英文對照.精品文檔.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

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

3、nd nonsymmetric tensors for the three-dimensional case using C3D8R elements is different from the ordering specified in “Three-dimensional solid element library,” Section 14.1.4, and the ordering used in ABAQUS/Standard.C3D8R單元三維軸對稱及非軸對稱張量成分順序與“Three-dimensional solid element library,” Section 14.1.

4、4及ABAQUS/Standard中指定的順序不同。References· “User-defined mechanical material behavior,” Section 12.8.1· *USER MATERIALOverviewUser subroutine VUMAT: 用戶子程序VUMAT· is used to define the mechanical constitutive behavior of a material;· 用來定義材料的力學(xué)本構(gòu)關(guān)系· will be called for blocks of mate

5、rial calculation points for which the material is defined in a user subroutine (“Material data definition,” Section 9.1.2);· 可以被用戶子程序定義材料計算點(diǎn)調(diào)用· can use and update solution-dependent state variables;· 可以使用和更新結(jié)果依賴狀態(tài)變量· can use any field variables that are passed in;· 可以使用傳入的任何

6、場變量· is described further in “User-defined mechanical material behavior,” Section 12.8.1; and· 在“User-defined mechanical material behavior,” Section 12.8.1中詳細(xì)論述；· cannot be used in an adiabatic analysis.· 可以被用于絕熱分析Component ordering in tensors 張量組成順序The component ordering depends

7、 upon whether the tensor is symmetric or nonsymmetric.張量組成順序取決于其是否為對稱或非對稱張量。Symmetric tensors 對稱張量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

8、 are first and then the indirect components, beginning with the 12-component. For example, a stress tensor contains ndir direct stress components and nshr shear stress components, which are passed in as 對于如同應(yīng)力及應(yīng)變張量等的對稱張量，含有ndir+nshr分量，分量的序號按張量索引號的自然排序給出。首先是直接分量，然后是從12分量開始的間接分量。例如，包含ndir正應(yīng)力分量及nshr的剪應(yīng)

9、力張量的應(yīng)力張量被按照下面的順序傳入Component2-D Case3-D Case12345 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

10、 非對稱張量For nonsymmetric tensors there are ndir+2*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, beginning with the 12-component. For example, the deformation gradient is passed

11、as 對于非對稱張量含有ndir+2*nshr分量，分量的順序按照張量索引號的自然排序給出。首先是直接分量，其次是從12分量開始的間接分量。例如，位移梯度按照下面的順序傳遞Component2-D Case3-D Case123456 7 8 9 Initial calculations and checks 最初計算和檢查In the data check phase of the analysis ABAQUS/Explicit calls user subroutine VUMAT with a set of fictitious strains

12、and a totalTime and stepTime both equal to 0.0. This 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)變、總時間及時間步都為0。這作為對用戶本構(gòu)關(guān)系的一個檢查，基于計算得到

13、的初始材料波速來計算等效初始材料屬性。Defining local orientations 定義局部方向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 syst

14、em in which the basis system rotates with the material. If a user-specified coordinate system (“Orientations,” Section 2.2.5) is used, it defines the local material axes in the undeformed configuration.所有的應(yīng)力、應(yīng)變、延伸及狀態(tài)變量均按局部材料軸的方向。這些局部材料軸形成一個應(yīng)力與應(yīng)變分量存儲的基本系統(tǒng)。即這個基本系是隨著材料聯(lián)合轉(zhuǎn)動的坐標(biāo)系。如果使用用戶指定坐標(biāo)系，則它在無變形結(jié)構(gòu)中定義局部

15、材料軸。Special considerations for various element types不同單元類型的特殊考慮The use of user subroutine VUMAT requires special consideration for various element types.用戶子程序VUMAT的使用需要對不同的單元類型進(jìn)行特殊的考慮。Shell and plane stress elements殼及平面應(yīng)力單元You must define the stresses and internal state variables. In the case of she

16、ll or plane stress elements, you must define strainInc(*,3), the thickness strain 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)變增量。如果需要的話還需要定

17、義初始能量。如果沒有定義，那么ABAQUS/Explicit提供的能量平衡將沒有意義。Shell elements殼單元When VUMAT is used 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. S

18、ee “Shell section behavior,” Section 15.6.4, for guidelines on choosing this stiffness.當(dāng)使用VUMAT定義殼單元的材料響應(yīng)時，ABAQUS/Explicit不能計算單元的缺省橫向剪切剛度。因此，用戶需要定義單元的橫向剪切剛度。關(guān)于選擇橫行剪切剛度的詳細(xì)資料請參考“Shell section behavior,” Section 15.6.4Beam elements梁單元For beam elements the stretch tensor and the deformation gradient ten

19、sor are not available. For beams in space you must define the thickness strains, strainInc(*,2) and strainInc(*,3). strainInc(*,4) is the shear strain associated with twist. Thickness stresses, stressNew(*,2) and stressNew(*,3), are assumed to be zero and any values you assign are ignored.對于梁單元不能使用拉

20、伸張量及位移梯度張量。對于空間梁，用戶必須定義厚度應(yīng)變、應(yīng)變增量（*，2）及應(yīng)變增量（*，3）。應(yīng)變增量（*，4）時與扭曲有關(guān)的剪應(yīng)變。厚度應(yīng)力，（*，2）及（*，3）假定為0，并且用戶分配的任何相關(guān)張量都被忽略。Deformation gradient 位移梯度The polar decomposition of the deformation gradient is written as , where and are the right and left symmetric stretch tensors, respectively. The constitutive model is

21、defined in a corotational coordinate system in which the basis system rotates with the material. All stress and strain tensor quantities are defined with respect to the corotational basis system. The right stretch tensor, , is used. The relative spin tensor represents the spin (the antisymmetric par

22、t of the velocity gradient) defined with respect to the corotational basis system.位移梯度寫成，其中及分別為右邊及左邊的對稱拉伸張量。本構(gòu)模型定義為聯(lián)合旋轉(zhuǎn)坐標(biāo)系，在該坐標(biāo)系中基系隨著材料轉(zhuǎn)動。所有的應(yīng)力和應(yīng)變張量值按照聯(lián)合旋轉(zhuǎn)坐標(biāo)系定義。使用右邊的拉伸位移。相應(yīng)的旋轉(zhuǎn)張量代表與聯(lián)合旋轉(zhuǎn)基系相應(yīng)的轉(zhuǎn)動。Special considerations for hyperelasticity超彈性的特殊考慮Hyperelastic constitutive models in VUMAT should be defin

23、ed in a corotational coordinate system in which the basis system rotates with the material. This is most effectively accomplished by formulating the hyperelastic constitutive model in terms of the stretch tensor, , instead of in terms of the deformation gradient, . Using the deformation gradient can

24、 present some difficulties because the deformation gradient includes the rotation tensor and the resulting stresses would need to be rotated back to the corotational basis. 在VUMAT中的超彈性本構(gòu)模型可以被定義在聯(lián)合選擇坐標(biāo)系中。這可以通過用拉伸張量表示的超彈性本構(gòu)模型很好的實(shí)現(xiàn)，而不是使用位移梯度來表示。使用位移梯度可能會帶來一些困難，因為位移梯度包括旋轉(zhuǎn)張量并且導(dǎo)致應(yīng)力需要選擇返回到聯(lián)合旋轉(zhuǎn)基系。Objective

25、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 obtained with user subroutine VUMAT may differ from that obtained with a built-in ABAQUS material model. For example, most material models used w

26、ith solid (continuum) elements in ABAQUS/Explicit employ the Jaumann stress rate. This difference in the formulation will cause significant differences in the results only if finite rotation of a material point is accompanied by finite shear. For a discussion of the objective stress rates used in AB

27、AQUS, 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)力率可能會與在ABAQUS建立的材料模型獲得的應(yīng)力率有所不同。例如，在ABAQUS/Explicit中大多數(shù)實(shí)體（連續(xù)）單元材料模型使用Jaumann應(yīng)力率。只要材料點(diǎn)的有限旋轉(zhuǎn)伴隨有限剪切，這種表達(dá)方式的不同將導(dǎo)致計算結(jié)果的明顯差異。關(guān)于ABAQUS中使用的目標(biāo)應(yīng)力率的詳細(xì)討論參考“Stress rates,” Section 1.5.3 of

28、 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 “User-defined mechanical material behavior,” Section 12.8.1). You must specify the state variable number controlling the element deletion flag when

29、you allocate space for the solution-dependent state variables, as explained in “User-defined mechanical material behavior,” Section 12.8.1. 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, while a value of zero

30、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 the block. ABAQUS/Explicit

31、 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)可以被從模型中刪除（參考“User-defined mechanical material behavior,” Section 12.8.1)。當(dāng)用戶給結(jié)果依賴狀態(tài)變量分配空間時，用戶需要指定控制單元刪除標(biāo)示的狀態(tài)變量號，在“User-defined me

32、chanical material behavior,” Section 12.8.1中進(jìn)行詳細(xì)說明。在VUMAT中刪除狀態(tài)變量可以被賦予1或者0。1表示材料點(diǎn)時激活的，0表示ABAQUS/Explicit將通過設(shè)定應(yīng)力為0刪除材料點(diǎn)。在分析過程中傳遞給用戶子程序VUMAT的材料點(diǎn)結(jié)構(gòu)保持不變；刪除的材料點(diǎn)沒有從塊中移走。ABAQUS/Explicit將傳遞0應(yīng)力及應(yīng)變給所有刪除的材料點(diǎn)。一旦一個材料點(diǎn)被標(biāo)示為刪除，該材料點(diǎn)將不能夠被再次激活。User subroutine interface 用戶子程序 subroutine vumat(C Read only (unmodifiable)v

33、ariables - 1 nblock, ndir, nshr, nstatev, nfieldv, nprops, lanneal, 2 stepTime, totalTime, dt, cmname, coordMp, charLength, 3 props, density, strainInc, relSpinInc, 4 tempOld, stretchOld, defgradOld, fieldOld, 5 stressOld, stateOld, enerInternOld, enerInelasOld, 6 tempNew, stretchNew, defgradNew, fi

34、eldNew,C Write only (modifiable) variables - 7 stressNew, stateNew, enerInternNew, enerInelasNew )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 stretchO

35、ld(nblock,ndir+nshr), 4 defgradOld(nblock,ndir+nshr+nshr), 5 fieldOld(nblock,nfieldv), stressOld(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,nfie

36、ldv), 1 stressNew(nblock,ndir+nshr), stateNew(nblock,nstatev), 2 enerInternNew(nblock), enerInelasNew(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

37、the increment.在增量結(jié)束時每個材料點(diǎn)的應(yīng)力張量。stateNew (nblock, nstatev) State variables at each material point at the end of the increment. You define the size of this array by allocating space for it (see “User subroutines: overview,” Section 25.1.1, for more information).增量結(jié)束時每個材料點(diǎn)的狀態(tài)變量。用戶通過分配空間來定義該矩陣的大小（更多的資料參

38、考“User subroutines: overview,” Section 25.1.1）。Variables that can be updated 可以更新的變量enerInternNew (nblock) Internal energy per unit mass at each material point at the end of the increment.增量結(jié)束時每個材料點(diǎn)單位質(zhì)量的內(nèi)能。enerInelasNew (nblock) Dissipated inelastic energy per unit mass at each material point at the

39、 end of the increment.增量結(jié)束時每個材料點(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)號。ndir Number of direct components in a symmetric tensor.對成張量的直接分量號。nshr Number of indirect components in a symmetric tensor.對稱張量的間接分量號。ns

40、tatev Number of user-defined state variables that are associated with this material type (you define this as described in “Allocating space” in “User subroutines: overview,” Section 25.1.1).與材料類型相關(guān)的用戶定義狀態(tài)變量號。nfieldv Number of user-defined external field variables.用戶定義外部場變量號。nprops User-specified num

41、ber of user-defined material properties.用戶定義材料屬性的用戶指定號。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 indicates that this is an annealing process and you shou

42、ld 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表明這是一個退火過程，并且如果需要的話用戶需要重新初始化內(nèi)部狀態(tài)變量stateNew。ABAQUS/Explic

43、it將自動設(shè)置應(yīng)力，延展性及在退火過程中0值狀態(tài)。stepTime Value of time since the step began.從時間步開始時的的時間totalTime Value of total time. The time at the beginning of the step is given by totalTime - stepTime.總時間值。時間步開始時的時間定義為totalTime - stepTimedt Time increment size.時間增量大小。cmname User-specified material name, left justified

44、. 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)以ABQ_字母開頭賦名。為了避免沖突，用戶不能使用ABQ_作為cmname的開頭字母。coordMp(nblo

45、ck,*) Material point coordinates. It is the midplane material point for shell elements and the centroid for beam elements.材料點(diǎn)坐標(biāo)。對于殼單元為中平面材料點(diǎn)，對于梁單元為質(zhì)心。charLength(nblock) Characteristic element length. This is a typical length of a line across an element. For beams and trusses, it is a characteristic

46、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 thickness.特征單元長度。它是一個穿過單元線的特殊長度。對于梁及桁架，它表示沿單元軸的特征長度。對于膜及殼

47、，它表示參考面上的特征長度。對于對稱單元，它僅僅表示平面上的特征長度。對于粘聚單元它等于結(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 volumetric strain increment is very small. If an accurate v

48、alue 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)變增量非常小的情況下，該值可能出現(xiàn)精確性問題。在這種情況下如果需要精確的密度值，分析需要在雙精度下運(yùn)行。密度值不受質(zhì)量縮放比例的影響strainInc (nblock, ndir+nshr) Strain increment tensor at each materi

49、al point.每個材料點(diǎn)的應(yīng)變增量張量。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-D as .定義在聯(lián)合旋轉(zhuǎn)系中的每個材料點(diǎn)的增量相關(guān)選擇矢量。定義為，其中為速度梯度，的反對稱部分，在

50、三維條件下存儲為，在二維條件下存儲為。tempOld(nblock) Temperatures at each material point at the beginning of the increment.增量開始時每個材料點(diǎn)的溫度。stretchOld (nblock, ndir+nshr) Stretch tensor, , at each material point at the beginning of the increment defined from the polar decomposition of the deformation gradient by .增量開始時每

51、個材料點(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 (, , , , ).增量開始時每個材料點(diǎn)的位移梯度張量。在三維狀態(tài)下存儲為(, , , , , , , , ) 及在二維條件下存儲為(, , , , ).fieldOld (nblock, nfieldv) Values of the us

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

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

54、質(zhì)量消散的無彈性能。tempNew(nblock) Temperatures at each material point at the end of the increment.增量開始時每個材料點(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 deformation gradient by .defgradNew (nblock,ndir+2*n

55、shr) Deformation gradient tensor at each material point at the end of the increment. Stored in 3-D as (, , , , , , , , ) and in 2-D as (, , , , ).增量結(jié)束時每個材料點(diǎn)的位移梯度張量。三維狀態(tài)下存儲為(, , , , , , , , ) ，二維狀態(tài)下存儲為 (, , , , )。fieldNew (nblock, nfieldv) Values of the user-defined field variables at each material p

56、oint at the end of the increment.增量結(jié)束時每個材料點(diǎn)的用戶指定場變量值。Example: Using more than one user-defined material model例子：使用多個用戶定義材料模型To use more than one user-defined material model, the variable cmname can be tested for different material names inside user subroutine VUMAT, as illustrated below:要使用多種用戶

57、材料模型，用戶子程序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 subroutines containing the constitutive material mo

58、dels 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-case character string. VUMAT_MAT1 and VUMAT_MAT2分別為實(shí)際的包含MAT1及MAT2材料本構(gòu)的用戶材料子