我曾经犯的错误
1,发现456方向的应变是工程应变
2,全量法可以直接根据更新后的应变求应力,单独求ddsdde矩阵. 3,ddsdde何时更新对于收敛性很有影响.
【 在 Mikejwg (mikejwg) 的大作中提到: 】
: uptonow一直叫我在水木清华上写一点东西，我一直忙于毕业论文，再说也是困难一?.
: 出现，永远有解不完的问题，就本人体会，写一点东西： : 在umat中， : (我总结三点：
: 1，求试算应力，stress=stress+De*dstran (其中：De为弹性矩阵） : 2，据应力势函数(屈服函数）求解满足条件的应力，应变，
: 在此我们要注意试算应力与真实应力的关系：sigma_true=Sigma_tri-De*dpstran: d.. : an为塑性应变增量，据此满足屈服条件，F（Sigmaij）＝0，求得真实的应力 : 更新应力 stress=stress-De*dpstran : 3，（与2联合，弹塑性区域的转换）判断是否屈服， : 4，组集ddsdde(与个人的本构有关） : ...................
发信人: wjytj (赖子), 信区: FEA
: 1,发现456方向的应变是工程应变
说白了，因为在UMAT中应变是以列向量的形式存储的，而初始刚度DDSDDE的形式为： λ＋μ μ μ 0 0 0 μ λ＋μ μ 0 0 0 μ μ λ＋μ 0 0 0 0 0 0 μ 0 0 0 0 0 0 μ 0 0 0 0 0 0 μ
所以，应变是工程应变，而如果DDSDDE的形式为（直接用张量表示）： λ＋μ μ μ 0 0 0 μ λ＋μ μ 0 0 0 μ μ λ＋μ 0 0 0 0 0 0 2μ 0 0 0 0 0 0 2μ 0 0 0 0 0 0 2μ
那么，应变就应该是应变张量中的应变了，也就是说，采用工程应变只是为了方便应用而
已，比如，求弹性应变能，可以直接将应力列向量的转置乘以应变列向量，显然，这在编
程中很容易实现。
: 2,全量法可以直接根据更新后的应变求应力,单独求ddsdde矩阵. 同感，应该可以减少不少工作量。
: 3,ddsdde何时更新对于收敛性很有影响.
希望air1能够说得更详细些，比如在开始是可以采用初始刚度，在后来可以采用切线刚度
，可是如何控制，编程如何实现呢？
【 在 air1 的大作中提到: 】 : 我曾经犯的错误
: 1,发现456方向的应变是工程应变
: 2,全量法可以直接根据更新后的应变求应力,单独求ddsdde矩阵. : 3,ddsdde何时更新对于收敛性很有影响.
: 【 在 Mikejwg (mikejwg) 的大作中提到: 】
: : uptonow一直叫我在水木清华上写一点东西，我一直忙于毕业论文，再说也是困难?.. : : 出现，永远有解不完的问题，就本人体会，写一点东西： : : 在umat中， : : (我总结三点：
: : 1，求试算应力，stress=stress+De*dstran (其中：De为弹性矩阵） : ...................
发信人: wjytj (赖子), 信区: FEA
: 2,全量法可以直接根据更新后的应变求应力,单独求ddsdde矩阵.
在我的本构中，DDSDDE的确定需要知道STRESS中的某些量，那么
请问，按照你的意思和经验，在UMAT中是否可以先得到STRESS(比如， 利用全量法），然后再得到DDSDDE？
如果可以的话，就应该可以减少迭代求解的次数了。 谢谢！！
假定采用mises屈服准则。
1 在调用UMAT之前，ABAQUS传递给UMAT本次增量开始时的应力sigma(0)，总应变E,应变增
量delta(E)。
状态变量保存：弹性应变，塑性应变，等效塑性应变。
（请问，状态变量保存的弹性应变＋塑性应变是否等于ABAQUS传递给UMAT的总应变？？？
？？） 2，然后在UMAT中利用上述的ABAQUS传递的量和状态变量得到DDSDDE矩阵，然后返回给
ABAQUS，ABAQUS根据delta(sigma)=ddsdde*delta(E)，并且得到本次增量结束时的应力
sigma(1)=sigma(0)+delta(sigma)=sigma(0)+ddsdde*delta(E)
3 然后更新本次增量结束时的状态变量：弹性应变，塑性应变，等效塑性应变以供下次调
用UMAT
请问手册上UMAT必须更新应力，可是根据上述我的理解好像是ABAQUS根据UMAT提供的本
次增量的DDSDDE在ABAQUS中更新，请问到底是怎么回事？谢谢！
个人觉得，在Abaqus中更新应力是必须的，而且必须准确，这是你的本构的核心；DDSDDE
则可以是不准确的，比如，你可以保持DDSDDE一直为初始刚度矩阵，当然，这样做的问题
会是收敛速度较慢或者根本不收敛。
在Manual中也说过，DDSDDE不影响结果，前提是，只要能得到结果。
与大家共同讨论。
【 在 uptonow 的大作中提到: 】 : 假定采用mises屈服准则。
: 1 在调用UMAT之前，ABAQUS传递给UMAT本次增量开始时的应力sigma(0)，总应变E,应..
: 量delta(E)。
: 状态变量保存：弹性应变，塑性应变，等效塑性应变。 : （请问，状态变量保存的弹性应变＋塑性应变是否等于ABAQUS传递给UMAT的总应变?.
: ？？）
: 2，然后在UMAT中利用上述的ABAQUS传递的量和状态变量得到DDSDDE矩阵，然后
返回?.
: ABAQUS，ABAQUS根据delta(sigma)=ddsdde*delta(E)，并且得到本次增量结束时的应..
: sigma(1)=sigma(0)+delta(sigma)=sigma(0)+ddsdde*delta(E)
: 3 然后更新本次增量结束时的状态变量：弹性应变，塑性应变，等效塑性应变以供下.. : ................... --
期盼了好久的abaqus6.3终于安装成功
寄希望于此来成功import我的.igs模型（unigraphics v18生成的） 但是，结果却不尽如人意，现将我的感受说一下：
1. 6.3在import进模型后，它的形状基本上与真实形状差不多,比6.2更像原始模型 但是，依然是invalid geometry
2. 经过delete invalid geometry后，它确实成了imprecise geometry(意味着你 以划分网格了），但是此时的形状已经面目全非了，删除了太多的东西了，得到的 imprecise geometry 就毫无意义了
3. 如果我用automated repair，它还是invalid geometry 我现在急需把文件import进后来进行下一步工作， 而老板出国前定的任务是11月底出结果 我现在连模型都是个问题了 急急急！！！ 请大侠帮忙！！！！！
救救我吧！给个好的建议吧
用PATRAN导入，生成网格后用ABAQUS计算
试试parasolid的文件格式，iges格式先天不足，尤其是表达曲面。
对了俺不知道abaqus支持parasolid文件不？不过一般情况下比igs格式
好点。
State variables
Since the redefinition of field variables in USDFLD is local to the current
increment (field variables are restored to the values
interpolated from the nodal values at the start of each increment), any history dependence required to update material properties by
using this subroutine must be introduced with user-defined state variables.
The state variables can be updated in USDFLD and then passed into other user
subroutines that can be called at this material point,
such as those listed above. The number of such state variables is specified with the *DEPVAR option...
要引入材料的非线性，需要usdfld;而要引入历史相关，就得要靠statev帮忙了！
(我的理解不全面，可以自己参考manual中的这一段.版主g上吧，我以后忘了看文摘区
就行了
感谢！
现在弄明白了，原来我的文件里多了几行空行，没想到abaqus的inp文件不允许有空行的。
感谢wjytj的帮助！刚开始用abaqus，太弱了。 【 在 wjytj 的大作中提到: 】 : :*NGEN : :1,111,11 : :11,121,11
: **此时，12，22，23……110节点并未生成，所以生成12，22以后的节点是没有意义?.. : 看看这样行否？ : *NODE : 1,0,10 : 11,10,10 : 111,0,0 : 121,10,0 : ...................
: solution-dependent variable
从字面上理解，即为―与解答有关的变量‖，它的值与每一增量结束时的变量有关，常用
于UMAT中，即SDV，可以由*DEPDV来定义其个数。