VUMAT-中英文对照

1、25.3.4 VUMATUser subroutine to define material behavior.定义材料本构用户子程序Product: ABAQUS/ExplicitWarning: The use of this user subroutine generally requires considerable expertise. You are cautioned that the implementation of any realistic constitutive model requires extensive development and testing. Ini

2、tial testing on a single-element model with prescribed traction loading is strongly recommended.注意：用户子程序的使用通常需要一定的专长。用户需要知道执行任何实际的本构模型需 要大量的试验数据。强烈建议用户对用户子程序进行在指定拉力作用下单个单元的验证测 试。The component ordering of the symmetric and nonsymmetric tensors for the three-dimensional case using C3D8R elements is di

3、fferent from the ordering specified in “Thee-dimensional solid element library, " Section, On.d.the ordering used in ABAQUS/Standard.C3D8R单元三维轴对称及非轴对称张量成分顺序与“Threelimensional solid elementlibrary, " Section 14.ABAQUS/Standard 中指定的顺序不同。References?“ Usr-defined mechanical material behavior,&

4、quot; Section 12.8.1? *USER MATERIALOverviewUser subroutine VUMAT:用户子程序VUMAT? is used to define the mechanical constitutive behavior of a material;?用来定义材料的力学本构关系? will be called for blocks of material calculation points for which the material is defined in a user subroutine (" Material data def

5、inition, " 1 .Section 9.?可以被用户子程序定义材料计算点调用? can use and update solution-dependent state variables;?可以使用和更新结果依赖状态变量? can use any field variables that are passed in;?可以使用传入的任何场变量? is described further in “ Usedefined mechanical material behavior," Sec0inxni2.8.? 在 “Usedefined mechanical mate

6、rial behavior," SectioniU28i；? cannot be used in an adiabatic analysis.?可以被用于绝热分析Component ordering in tensors张量组成顺序The component ordering depends upon whether the tensor is symmetric or nonsymmetric.张量组成顺序取决于其是否为对称或非对称张量。Symmetric tensors 对称张量For symmetric tensors such as the stress and strain

7、 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, beginning with the 12-component. For example, a stress tensor contains ndir direct stress components and

8、 nshr shear stress components, which are passed in as 对于如同应力及应变张量等的对称张量，含有ndir+nshr分量，分量的序号按张量索引号的自然排序给出。首先是直接分量，然后是从12分量开始的间接分量。例如，包含ndir正应力分 量及nshr的剪应力张量的应力张量被按照下面的顺序传入Component2-D Case 3-D Case1冉 W_2333ll33412|52360The shear strain components in user subroutine VUMAT are stored as tensor componen

9、ts and not as engineering components; this is different from user subroutine UMAT in ABAQUS/Standard, which uses engineering components.Nonsymmetric tensors 非对称张量For nonsymmetric tensors there are ndir+2*nshr components, and the component order is given as a natural permutation of the indices of the

10、 tensor. The direct components are first and then the indirect components, beginning with the 12-component. For example, the deformation gradient is passed as对于非对称张量含有ndir+2*nshr分量，分量的顺序按照张量索引号的自然排序给出。首先是 直接分量，其次是从12分量开始的间接分量。例如，位移梯度按照下面的顺序传递Component2-D Case3-D Case1Fh 1二2心 |3人|4L Ha561|F'n6际78

11、耳9Initial calculations and checks最初计算和检查 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. This is done as a check on your constitutive relation and to calculate the equivalent initial

12、material properties, based upon which the initial elastic wave speeds are computed.在ABAQUS/Explicit调用用户子程序VUMAT分析的数据检查阶段，小应变、总时间及时间步 都为0o这作为对用户本构关系的一个检查，基于计算得到的初始材料波速来计算等效初始 材料属性。Defining local orientations定义局部方向All stresses, strains, stretches, and state variables are in the orientation of the loca

13、l 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 material. If a user-specified coordinate system ("Orientations, " Sectionis?i2s5d, it

14、 defines the local material axes in the undeformed configuration.所有的应力、应变、延伸及状态变量均按局部材料轴的方向。这些局部材料轴形成一个应力 与应变分量存储的基本系统。即这个基本系是随着材料联合转动的坐标系。如果使用用户指 定坐标系，则它在无变形结构中定义局部材料轴。Special considerations for various element types不同单元类型的特殊考虑The use of user subroutine VUMAT requires special consideration for vari

15、ous element types. 用户子程序VUMAT的使用需要对不同的单元类型进行特殊的考虑。Shell and plane stress elements壳及平面应力单元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 strain increment. The internal energies can be defined if d

16、esired. If they are not defined, the energy balance provided by ABAQUS/Explicit will not be meaningful.用户必须定义应力和初始状态变量。在壳或平面应力单元的情况下，用户必须定义应变包括 (*,3),厚度应变增量。如果需要的话还需要定义初始能量。如果没有定义，那么ABAQUS/Explicit提供的能量平衡将没有意义。Shell elementsWhen VUMAT is used to define the material response of shell elements, ABAQUS

17、/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 Shell section behavior," Sectioor guide4ineson choosing this stiffness.当使用VUMAT定义壳单元的材料响应时，ABAQUS/Explicit不能计算单元的缺省横向剪切 刚度。因此，用户

18、需要定义单元的横向剪切刚度。关于选择横行剪切刚度的详细资料请参考“ Shell section behavior, " Section 15.6.4Beam 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 th

19、e shear strain associated with twist. Thickness stresses, stressNew（*,2） and stressNew（*,3）, are assumed to be zero and any values you assign are ignored.对于梁单元不能使用拉伸张量及位移梯度张量。对于空间梁，用户必须定义厚度应变、应变增量（*, 2）及应变增量（*, 3）。应变增量（*, 4）时与扭曲有关的剪应变。厚度应力，（*, 2）及（*, 3）假定为0,并且用户分配的任何相关张量都被忽略。Deformation gradient位移梯度

20、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 defined in a corotational coordinate system in which the basis system rotates with the material. All stress and strain

21、 tensor quantities are defined with respect to the corotational basis system. The right stretch tensor, T , is used. The relative spin tensor''' r represents the spin (the antisymmetricpart of the velocity gradient) defined with respect to the corotational basis system.位移梯度写成F = R U = V-

22、R,其中U及V分别为右边及左边的对称拉伸张量。本 构模型定义为联合旋转坐标系，在该坐标系中基系随着材料转动。所有的应力和应变张量值 按照联合旋转坐标系定义。使用右边的拉伸位移U。相应的旋转张量W-fi代表与联合旋转基系相应的转动。Special considerations for hyperelasticity超弹性的特殊考虑Hyperelastic constitutive models in VUMAT should be defined in a corotational coordinate system in which the basis system rotates with

23、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, I :; . Using the deformation gradient can present some difficulties because the deformation gradient includes the

24、rotation tensor and the resulting stresses would need to be rotated back to the corotational basis.在VUMAT中的超弹性本构模型可以被定义在联合选择坐标系中。这可以通过用拉伸张量U表示的超弹性本构模型很好的实现，而不是使用位移梯度F = R - U来表示。使用位移梯度可能会带来一些困难，因为位移梯度包括旋转张量并且导致应力需要选择返回到联合旋转基 系。Objective stress rates 目标应力率The Green-Naghdi stress rate is used when the

25、 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 with solid (continuum) elements in ABAQUS/Explicit employ the J

26、aumann 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 ABAQUS, see “ Stress rates, " Section 1.5.3 of the ABAQUS T

27、heory.Manual 在用户子程序VUMAT中使用Green-Naghdi应力率来定义材料的力学本构关系。通过用户子程序VUMAT获得的应力率可能会与在 ABAQUS建立的材料模型获得的应力率有所不同。 例如，在ABAQUS/Explicit中大多数实体(连续)单元材料模型使用Jaumann应力率。只要材料点的有限旋转伴随有限剪切，这种表达方式的不同将导致计算结果的明显差异。关于ABAQUS中使用的目标应力率的详细讨论参考“ Stress rates, " Section 1.5.3 of the ABAQUSTheory Manual.Material point deleti

28、on 材料点删除Material points that satisfy a user-defined failure criterion can be deleted from the model (see “ Usedefined mechanical material behavior," Se)ctYou1mUs1 specify the state variablenumber controlling the element deletion flag when you allocate space for the solution-dependent state vari

29、ables, as explained in “ Usedefined mechanical material behavior," SecTbe12.8.1deletion 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 indicates that ABAQUS/Explicit should delete the material po

30、int 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 will pass zero stresses and strain increments for all delet

31、ed material points. Once a material point has been flagged as deleted, it cannot be reactivated.满足用户定义的破坏准则的材料点可以被从模型中删除(参考“Usedefined mechanicalmaterial behavior, " Section 12.8用户给结果依赖状态变量分配空间时，用户需要指定控 制单元删除标示的状态变量号，在 “Usedefined mechanical material behavior, " Section 12.8.1 中进行详细说明。在 VU

32、MAT中删除状态变量可以被赋予1或者0。1表示材料点时激活的, 0表示ABAQUS/Explicit将通过设定应力为0删除材料点。在分析过程中传递给用户子程序 VUMAT的材料点结构保持不变；删除的材料点没有从块中移走。ABAQUS/Explicit将传递0 应力及应变给所有删除的材料点。一旦一个材料点被标示为删除，该材料点将不能够被再次 激活。User subroutine interface 用户子程序subroutine vumat(C Read only (unmodifiable)variables -1 nblock, ndir, nshr, nstatev, nfieldv, n

33、props, 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, fieldNew,C Write only (modifiable) variables -7 stressNew,

34、stateNew, enerInternNew, enerInelasNew )Cinclude 'vaba_param.inc'Cdimension 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 fiel

35、dOld(nblock,nfieldv), stressOld(nblock,ndir+nshr),6 stateOld(nblock,nstatev), enerInternOld(nblock),7 enerInelasOld(nblock), tempNew(nblock),8 stretchNew(nblock,ndir+nshr),9 defgradNew(nblock,ndir+nshr+nshr),10 fieldNew(nblock,nfieldv),1 stressNew(nblock,ndir+nshr), stateNew(nblock,nstatev),2 enerln

36、ternNew(nblock), enerlnelasNew(nblock), Ccharacter*80 cmnameCdo 100 km = 1,nblock user coding100 continuereturnendVariables to be defined被定义的变量stressNew (nblock, ndir+nshr)Stress tensor at each material point at the end of the increment.在增量结束时每个材料点的应力张量。stateNew (nblock, nstatev)State variables at e

37、ach material point at the end of the increment. You define the size of this array by allocating space for it (see “User subroutines: overview, " Sectioor 25ore1 information).增量结束时每个材料点的状态变量。用户通过分配空间来定义该矩阵的大小(更多的资料参考 “User subroutines: overview, " Section 25.1.1Variables that can be updated

38、可以更新的变量enerInternNew (nblock)Internal energy per unit mass at each material point at the end of the increment.增量结束时每个材料点单位质量的内能。enerInelasNew (nblock)Dissipated inelastic energy per unit mass at each material point at the end of the increment.增量结束时每个材料点单位质量的消散的无弹性能。Variables passed in for informatio

39、nnblockNumber of material points to be processed in this call to VUMAT.调用VUMAT的材料点号。ndirNumber of direct components in a symmetric tensor.对成张量的直接分量号。nshrNumber of indirect components in a symmetric tensor.对称张量的间接分量号。nstatevNumber of user-defined state variables that are associated with this material

40、 type (you definethis as described in “Allocating space " in “User subroutines: overview, ). “ Section 25.1.1与材料类型相关的用户定义状态变量号。nfieldvNumber of user-defined external field variables.用户定义外部场变量号。npropsUser-specified number of user-defined material properties.用户定义材料属性的用户指定号。lannealFlag indicating

41、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 should re-initialize the internal state variables, stateNew, if necessary. ABAQUS/E

42、xplicit will automatically set the stresses, stretches, and state to a value of zero during the annealing process.在退火处理过程中标示程序是否被调用。Laneal=0表明程序在正常力学增量过程中被调 用。Laneal= 1表明这是一个退火过程，并且如果需要的话用户需要重新初始化内部状态 变量stateNew。ABAQUS/Explicit将自动设置应力，延展性及在退火过程中0值状态。stepTimeValue of time since the step began.从时间步开始时

43、的的时间totalTimeValue of total time. The time at the beginning of the step is given by totalTime - stepTime. 总时间值。时间步开始时的时间定义为 totalTime - stepTimedtTime increment size.时间增量大小。cmnameUser-specified material name, left justified. It is passed in as an upper-case character string.Some internal material mo

44、dels are given names starting with the“ ABQ_ character stravoid conflict, you should not use“ABQ_ as the leading string for cmname.用户指定材料名。按照大写字母传入。有些内而材料本构以ABQ字母开头赋名。为了避免冲突，用户不能使用ABQ_作为cmname的开头字母。coordMp(nblock,*)Material point coordinates. It is the midplane material point for shell elements and

45、the centroid for beam elements.材料点坐标。对于壳单元为中平面材料点，对于梁单元为质心。charLength(nblock)Characteristic element length. This is a typical length of a line across an element. For beams and trusses, it is a characteristic length along the element axis. For membranes and shells, it is a characteristic length in th

46、e reference surface. For axisymmetric elements, it is a characteristic length in the = £ plane only. For cohesive elements it is equal to the constitutive thickness.特征单元长度。它是一个穿过单元线的特殊长度。对于梁及桁架，它表示沿单元轴的特征 长度。对于膜及壳，它表示参考面上的特征长度。对于对称单元，它仅仅表示r匚平面上的特征长度。对于粘聚单元它等于结构厚度。props(nprops)User-supplied mate

47、rial 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 value of the density is required in such cases, the analysis should be run in doubl

48、e precision. This value of the density is not affected by mass scaling.材料点的当前密度。在体积应变增量非常小的情况下，该值可能出现精确性问题。在这种 情况下如果需要精确的密度值，分析需要在双精度下运行。密度值不受质量缩放比例的影 响strainInc (nblock, ndir+nshr)Strain increment tensor at each material point.每个材料点的应变增量张量。relSpinInc (nblock, nshr)Incremental relative rotation vect

49、or at each material point defined in the corotational system.Defined as ''，where is the antisymmetric part of the velocity gradient,：., and: R J - . Stored in 3-D as】；and in 2-D as 乙).定义在联合旋转系中的每个材料点的增量相关选择矢量。定义为才皿- n),其中w为速度梯度L , fi = R R,r的反对称部分，在三维条件下存储为 阳L3.2D在二维条件下存储为21 )|。tempOld(nbloc

50、k)Temperatures at each material point at the beginning of the increment.增量开始时每个材料点的温度。stretchOld (nblock, ndir+nshr)Stretch tensor, f , at each material point at the beginning of the increment defined from the polar decomposition of the deformation gradient by ： R 1.增量开始时每个材料点的，拉伸张量 U。defgradOld (nb

51、lock,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 (,',1).增量开始时每个材料点的位移梯度张量。在三维状态下存储为"&?,屿色 %, 2$,旧M匹I玛2, Fi蜀及在二维条件下存储为 呵 如,/$, Ft,理1).fieldOld (nblock, nfieldv)Values of the user-defined fie

52、ld variables at each material point at the beginning of the increment.增量开始时每个材料点的用户定义场变量值。stressOld (nblock, ndir+nshr)Stress tensor at each material point at the beginning of the increment.增量开始时每个材料点的应力张量。stateOld (nblock, nstatev)State variables at each material point at the beginning of the incre

53、ment.增量开始时每个材料点的状态变量。enerInternOld (nblock)Internal energy per unit mass at each material point at the beginning of the increment.增量开始时每个材料点的单位质量的内能enerInelasOld (nblock)Dissipated inelastic energy per unit mass at each material point at the beginning of the increment.增量开始时每个材料点的单位质量消散的无弹性能。tempNew(

54、nblock)Temperatures at each material point at the end of the increment.增量开始时每个材料点的温度。stretchNew (nblock, ndir+nshr)Stretch tensor, u, at each material point at the end of the increment defined from the polar decomposition of the deformation gradient by I.defgradNew (nblock,ndir+2*nshr)Deformation gr

55、adient tensor at each material point at the end of the increment. Stored in 3-D as ( i, .,一户",，恬区产口，凡儿.丘二,)and in 2-D as Q”,产，脸/,餐I).增量结束时每个材料点的位移梯度张量。三维状态下存储为旧j,月上 Ed,石即几璃三”， En Tz,尸应，二维状态下存储为 Qi, Fb, 8舄F? Bi)。fieldNew (nblock, nfieldv)Values of the user-defined field variables at each materia

56、l point at the end of the increment.增量结束时每个材料点的用户指定场变量值。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') thencall VUMAT_MAT1( argument_list)else if (cmname(1:4) .eq. 'MAT2') thencall 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分另U为实际的包含 MAT1及MAT2材料本构的用户材料子程序。子程序VUMAT 心匕仅仅作为一个地址。Example: Elastic/plastic material with kinematic hardening例子