NBLOCK: 在调用Vumat时需要用到的材料点的数量
Ndir：对称张量中直接应力的数量（sigma11,sigma22,sigma33）
Nshr：对称张量中间接应力的数量(sigma12, sigma13, sigma23)
Nstatev:与材料类型相关联的用户定义的状态变量的数目
Nfieldv：用户定义的外场变量的个数
Nprops:用户自定义材料属性的个数
Lanneal：指示是否在退火过程中被调用例程的标志。Lanneal=0，指示在常规力学性能增量，例程被调用。Lanneal=1表示，这是退火过程，你应该重新初始化内部状态变量，红烧牛蛙的做法
stepTime：步骤开始后的数值
totalTime：总时间
Dt：时间增量值
Cmname：用户自定义的材料名称，左对齐。它是通过字符串传递的。一些内部材料模型是以“ABQ_”字符串开头给定的名称。为了避免冲突，你不应该在“cmname”中使用“ABQ_”作为领先字符串。
coordMp(nblock,*)：材料点的坐标值。它是壳单元的中层面材料点，梁和管（pipe）单元的质心。
charLength(nblock)：
特征元素长度，是基于几何平均数的默认值或用户子程序VUCHARLENGTH中定义的用户特征元长度。
props(nprops)：用户使用的材料属性
density(nblock)：中层结构的物质点的当前密度
strainInc (nblock, ndir+nshr)：每个物质点处的应变增量张量
relSpinInc (nblock, nshr)：在随转系统中定义的每个物质点处增加的相对旋转矢量
tempOld(nblock)：物质点开始增加时的温度。
defgradOld (nblock,ndir+2*nshr)：在增量开始时，每个物质点出的变形梯度张量，在3d中形为（F11, F22，F33，F12，F23，F31，F21，F32，F13）,在2d中形为（F11，F22，F33，F12，F21）
stretchOld (nblock, ndir+nshr)
fieldOld (nblock, nfieldv)：在增量开始时，每个物质点处用户定义场变量的值
stressOld (nblock, ndir+nshr)：在增量开始时，每个物质点处的应力张量：
stateOld (nblock, nstatev)：在增量开始时我的故事，每个物质点处的状态变量：
tempNew(nblock)：在增量结束时，每个物质点处的温度
defgradNew (nblock,ndir+2*nshr)：在增量结束时，每个物质点出的变形梯度张量，在3d
中形为（F11, F22，F33，F12，F23，F31，F21，F32，F13）,在2d中形为（F11，F22，F33，F12，F21）
fieldNew (nblock, nfieldv)：在增量开始时，每个物质点处用户定义长变量的值
Example: Elastic/plastic material with kinematic hardening
As a simple example of the coding of subroutine VUMAT, consider the generalized plane strain ca for an elastic/plastic material with kinematic hardening. The basic assumptions and definitions of the model are as follows.
水杉树
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 space 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
The von Mis yield surface is defined as
where is the uniaxial equivalent yield stress. The von Mis yield surface is a cylinder in deviatoric stress space with a radius of
For the kinematic hardening model, R is a constant. The normal to the Mis yield surface can be written as
We decompo the strain rate into an elastic and plastic part using an additive decomposition:
The plastic part of the strain rate is given by a normality condition
where the scalar multiplier must be determined. A scalar measure of equivalent plastic strain rate is defined by
The stress rate is assumed to be purely due to the elastic part of the strain rate and is expresd in terms of Hooke's law by
where and are the Lamés constants for the material.
The evolution law for is given as
where H is the slope of the uniaxial yield stress versus plastic strain curve.
During active plastic loading the stress must remain on the yield surface, so that
The equivalent plastic strain rate is related to by
The kinematic hardening constitutive model is integrated in a rate form as follows. A trial elastic stress is computed as 狼为什么嚎叫
where the subscripts and refer to the beginning and end of the increment, respectively. If the trial stress does not exceed the yield stress, the new stress is t equal to the trial stress. If the yield stress is exceeded, plasticity occurs in the increment. We then write the incremental analogs of the rate equations as
一声兄弟where
红娘是哪部作品中的人物
From the definition of the normal to the yield surface at the end of the increment, ,
This can be expanded using the incremental equations as
Taking the tensor product of this equation with , using the yield condition at the end of the increment, and solving for :
The value for is ud in the incremental equations to determine , , and .
subroutine vumat(
C Read only -
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,
3 stressOld, stateOld, enerInternOld, enerInelasOld,
6 tempNew, stretchNew, defgradNew, fieldNew,
C Write only -
5 stressNew, stateNew, enerInternNew, enerInelasNew )
C
include 'vaba_param.inc'
C
C J2 Mis Plasticity with kinematic hardening for plane
C strain ca.
C Elastic predictor, radial corrector algorithm.
C
C The state variables are stored as:
C STATE(*,1) = back stress component 11
C STATE(*,2) = back stress component 22
广场舞视频大全免费下载C STATE(*,3) = back stress component 33昆明导游
