后处理中contour与assign的结果不一样？

xuw08

大家有没有将contour结果与assign结果进行对比？
应力最大值两个不一样？
不知道为何，是单元应力计算方法不一样吗

水面无波

assign为单元应力
contour为单元应力平均后的节点应力

wxwieandy

2楼说的不是很明白，在自己说下！

xuw08

assign为单元应力
contour为单元应力平均后的节点应力
水面无波 发表于 2010-6-11 21:02

不是很明白，他们计算单元应力的方法不一样？
分析计算时，以哪个为准？
谢谢

水面无波

先copy一段理论上的说法：有限单元法中，每一单元的特性用单元刚度矩阵来表示，每一结构构件的力与位移之间的关系不是精确推导出来的，而是利用每一单元中近似的位移函数得到节点位移，然后计算积分点应变和应力，输出时在根据用户请求将积分点结果复制或线性外推至单元的节点上。
个人经验：
    无论是Nastran的op2、Dyna的d3plot还是Abaqus的odb，在各个求解器的结果文件里所储存的应力结果都是最原始的在积分点处计算后得到的单元应力，往往表示单元中心处的应力。若尝试用Patran、Hyperview等后处理软件查看单元应力的话结果是一样的。而节点应力是由单元应力或复制或外推至节点进行平均后得到的反映节点处的应力。不同后处理软件对此处理略有差异。
    如Hyperview中的simple法就是将围绕节点的各单元应力平均后直接作为节点应力（可参阅软件自带的帮助，里面对simple、advance、difference等不同方法都有说明）。有兴趣可以Nastran为例，在应力梯度较大处取一些较方正的四边形单元，只输出centre，单元应力是由靠近四个节点的4个积分点处应力插值得到的单元中心处的应力。若将单元四等分细化后再计算，用simple法，则在原单元中心新生成的节点处的节点应力基本就等于原单元的单元应力。
    以哪个为准的话可看具体情况，多数是都可以。区别就是所反映的应力实际位置不一样并且单元应力更精确。尤其靠近几何边界的话节点应力就不可信，若要关注的话不如输出corner。一般静强度计算的话倒也无所谓，只要有经验作出正确判断就行了；疲劳计算尤其是用E-N法分析局部应力应变作用下的疲劳强度，我更愿意用单元应变。