在应用ABAQUS梁单元时发现的BUG，希望各位高手积极讨论一下

wangrui

大家好，我现在也遇见了一个很bug的问题：我来详细描述一下，还请您还有其他高手指点一下：首先我做了一个桩土分析，然后上部是桁架，用beam单元模拟，下部是筒形基础和土体的实体模型，然后加载计算；之后我在ANSYS当总把上部桁架部分单独提取出来，施加相同的载荷，发现应力和变形的计算值差别很大，为了找出原因，我建立了几组对比模型，其中一组是ANSYS和ABAQUS的对比，即都是建立一个悬臂梁，施加同样的载荷，此时得到的应力和变形值是一致的；之后我建立第二组对比模型，即也是两组悬臂梁，但是这两组悬臂梁都是在ABAQUS中完成的，具体描述如下：首先，模型OA沿OY方向，在O点刚性固定，分别在A端施加一个X方向和Z方向的1000牛的载荷；理论上这两种工况都是单纯的额悬臂梁模型，计算结果应该是一致的；然后观察后处理数据发现，应力值差别很大，而变形和支反力等数据基本一致；于是，我初步估计认为是和梁截面的方向有关，然后又设置了以下方案，出现下面的结果：当设置梁截面n1轴为X轴时，施加沿X方向的1000牛的力，应力结果为5.256e-8，如图mise1所示；然后其他条件不变，只把载荷换成Z向的时候，应力值变成了124.1MPa，如图mise2所示；所以这个问题就奇怪了；之后我在修改n1轴的方向，发现一个现象就是，无论n1轴时哪个方向，当只在X方向施加1000牛的力时得到的应力和只在Z方向施加1000牛得到的应力，平方，再求和，再开根号的时候，得到的数值都是124.1。
    对于以上情况，还请各位多多帮忙，在附件当中job-7为模型的inp文件，mise1和mise2分别为两种情况时得到的应力分布图

zytsang

一维梁单元来计算Mises应力不合适吧？梁单元算截面弯矩和剪力还行，算Mises应力的话还是用壳单元或实体单元。
LZ模型的梁截面是pipe，中心对称，所以梁截面的方向应该对结果无影响。我估计是ABAQUS输出一维梁单元的的Mises应力时，采样点的选取跟梁截面的方向有关。
如果Mises应力的采样点选了弯曲的中性轴，那其数值自然很小。（LZ附件图mise1）
如果把截面方向转90度，保持外力方向不变，采样点就选到了最外纤维，那数值自然变最大了。（LZ附件图mise2）
当截面方向在0~90度之间，那采样点等于位于截面不同高度的位置。LZ说的第三种情形，其实是把处于X方向外力下该点的Mises应力和处于Z方向外力下该点的Mises应力，用沟谷定理，求出了位于截面离中性轴最远纤维的Mises应力。
LZ还是用壳单元计算吧

Hansha

梁截面是什么形状？

wangrui

Hansha 发表于 2012-7-6 18:32
梁截面是什么形状？

管状的截面

wangrui

Hansha 发表于 2012-7-6 18:32
梁截面是什么形状？

管形截面，这个帖子别封啊，描述的最详细的，我急等结果，拜托~

Hansha

wangrui 发表于 2012-7-6 08:05
管形截面，这个帖子别封啊，描述的最详细的，我急等结果，拜托~

抱歉，是误操作。正在设法恢复。

wangrui

zytsang 发表于 2012-7-7 11:59
一维梁单元来计算Mises应力不合适吧？梁单元算截面弯矩和剪力还行，算Mises应力的话还是用壳单元或实体单元 ...

太感谢各位了，那么这样的话就像我说的这种情况对下部结构的受力有没有影响呢？也就是说我上部是桁架结构，下部是筒形基础和土体，如果上部应力仅仅是显示不对的话，没有，关系，桁架的应力分布不是我所要关注的重点，但是它对于力的传递是不是有影响呢？也就是说，这个情况是否会影响到下部筒形基础和土体的受力分布呢？

wangrui

Hansha 发表于 2012-7-7 08:32
抱歉，是误操作。正在设法恢复。

嘿嘿，谢谢斑竹支持啊。。。

zytsang

wangrui 发表于 2012-7-7 17:53
太感谢各位了，那么这样的话就像我说的这种情况对下部结构的受力有没有影响呢？也就是说我上部是桁架结构 ...

只是为了传力的话是没影响的

wangrui

zytsang 发表于 2012-7-7 21:53
只是为了传力的话是没影响的

用ANSYS校核后发现应力太大了，这下等于是前期工作都白做了。。。。

wangrui

zytsang 发表于 2012-7-7 21:53
只是为了传力的话是没影响的

版主，您好，我今天又尝试了一下ANSYS中的shell63和ABAQUS中的shell以及实体单元，就是一块方形板，然后在四个角点加固支，在中部施加一个集中载荷，结果发现其应力分布又有很大的不同？我搜索了论坛好像说ABA不适合集中载荷，这是什么情况呢？

zytsang

wangrui 发表于 2012-7-9 14:56
版主，您好，我今天又尝试了一下ANSYS中的shell63和ABAQUS中的shell以及实体单元，就是一块方形板，然后 ...

把ansys和abaqus模型传上来，这样大家也好给意见

wangrui

zytsang 发表于 2012-7-9 22:37
把ansys和abaqus模型传上来，这样大家也好给意见

模型文件在压缩文件中上传，请各位再次致电一下，谢谢

lcddcl

ANSYS中是扩展整个截面显示的等效应力，ABAQUS中是显示截面中的一个点的等效应力，结果肯定差别较大。

lcddcl

基本上ABAQUS中的结果不能用，不是求解问题，是显示内容不同。