umat里怎么调取单元的最大Mises应力值

小猴子123

我在UMAT子程序里一个方程要用到分析部位所有单元的最大MISES应力值，怎么找到这个值啊求大神指点

cartographer

小猴子123 发表于 2016-4-21 16:38
C      CALCULATE MISES STRESS   
      SMISES=(STRESS(1)-STRESS(2))**2+(STRESS(2)-STRESS(3))**2
   ...

DIMENSION EleMIses(Eletotal,Nptotal)(单元总数 X 单个单元积分点总数目)
需要保证单元数目 是从1开始顺序编号。

common EleMIses

****
***
***

EleMIses(NOEL,NPT)=Smises
直接赋值就可以
其中Noel 为单元编号（从1到Eletotal）
Npt 为积分点编号（从1到Nptotal）
=========

好久不用有点忘记了
其实不用数组也可以

直接设定公共变量 MisesTemp

在同一增量步中
和 每个计算得到的 Smises 做比较
就可以得到最大值了。MaxMises

==========
另外如果学习使用子程序URDFIL就更简单了

帮助文档里面的例子 就是查询最大值

cartographer

1. 计算单元的mises应力
2. 写个循环找到最大应力
3. 设置一个状态变量存储这个最大值
4. 调用

小猴子123

cartographer 发表于 2016-4-21 09:02
1. 计算单元的mises应力
2. 写个循环找到最大应力
3. 设置一个状态变量存储这个最大值

大神，第二步的写个循环的循环条件是怎么给的？umat里面计算出的mises应力是每个积分点处的应力吧

cartographer

小猴子123 发表于 2016-4-21 09:39
大神，第二步的写个循环的循环条件是怎么给的？umat里面计算出的mises应力是每个积分点处的应力吧 ...

可以设置公共数组 存储每个单元每个积分的的值

在每个增量步结束的时候，对数组做循环 找到最大的值  存储为公共变量

========

当前步应该只能用上一步的最大值

小猴子123

cartographer 发表于 2016-4-21 15:32
可以设置公共数组 存储每个单元每个积分的的值

在每个增量步结束的时候，对数组做循环 找到最大的值   ...

C      CALCULATE MISES STRESS   
      SMISES=(STRESS(1)-STRESS(2))**2+(STRESS(2)-STRESS(3))**2
     1                               +(STRESS(3)-STRESS(1))**2
      DO K1=NDI+1,NTENS
      SMISES=SMISES+6*STRESS(K1)**2
      END DO
      SMISES=SQRT(SMISES/2.0)

这样计算的MISES应力，SMISES就是一个数，不知道怎么指定某个单元某个积分点，不同单元不同积分点的SMISES值怎么调取赋给定义的数组的？

小猴子123

cartographer 发表于 2016-4-20 21:48
DIMENSION EleMIses(Eletotal,Nptotal)(单元总数 X 单个单元积分点总数目)
需要保证单元数目 是从1开始 ...

多谢大神指点

113soim

不知道楼主这个程序最终编出来么？想请教下