问一个关于弹塑性分析的问题，哪位帮忙看看？多谢。

tmpyy

试着做点弹塑性分析，用的是板单元，材料屈服强度为235MPa。发现用S4单元和S4R单元结果竟然不一样，如下图。好象用S4单元无法考虑塑性，最大Mises应力超过了屈服应力，是我什么地方错了么？

tmpyy

用s4单元的数据文件

tmpyy

用s4单元的结果

tmpyy

用s4R单元的数据文件

tmpyy

用s4R单元的结果

tmpyy

请帮忙看看，谢谢。

wjytj

你在S4的INP中将*el print那一行修改为：
*El Print
MISES,
然后看看dat文件中的MISES应力的最大值是多少？是不是235？
  
看来你没好好看Getting Started
How does the maximum value of Mises stress compare to the value reported in the contour plot generated earlier? Do the two maximum values correspond to the same point in the model? The Mises stresses shown in the contour plot have been extrapolated to the nodes, whereas the stresses written to the report file for this problem correspond to the element integration points. Therefore, the location of the maximum Mises stress in the report file is not exactly the same as the location of the maximum Mises stress in the contour plot. This difference can be resolved by requesting that stress output at the nodes (extrapolated from the element integration points and averaged over all elements containing a given node) be written to the report file. If the difference is large enough to be of concern, this is an indication that the mesh may be too coarse.
  
为什么S4R的结果看起来要好一些？因为S4R只有一个积分点。
  
BTW: 为什么要在S4R中加上*Section Controls以增加Hourglassing刚度？我试过默认值，结果一样阿。

pjwseu

对于最大Mises应力超过了屈服应力的问题，
ansys中用命令“ERESX,no”来处理。
不知abaqus中有没有类似的命令，只能加密网格吗？

wjytj

pjwseu wrote:
对于最大Mises应力超过了屈服应力的问题，
ansys中用命令“ERESX,no”来处理。
不知abaqus中有没有类似的命令，只能加密网格吗？

  
感觉在Contour中外插并平均只是为了图形表示更光滑而已，
而并不代表真实值，真实值实际上在dat文件中。
对ansys不熟，如果不平均，节点的值可能会有多个或者不连续，
图形表示起来岂不是很麻烦？

tmpyy

明白了。谢谢！

tmpyy

wjytj wrote:
BTW: 为什么要在S4R中加上*Section Controls以增加Hourglassing刚度？我试过默认值，结果一样阿。

  
我不是有意加的