ABAQUS做外插值出错, 是我的问题?

cdstudio · 发表于 2009-8-27 21:54:16

今天发现做一个ABAQUS外插值的时候发现了问题, 是这样的一种情况.
很简单的模型, 平面应变的四边形模型, 两个角的坐标(0,0)和(10, 10), 材料密度为2, E和nu<0.5随便给, 先用initial conditions, type=stress, gestatic引入初始地应力, 然后用static计算步加重力荷载, 用CPE4单元, 是一个简单得不能再简单的模型了, 现在看结果S22, 就会发现最大值不是0, 而是-10, 最小值也不是-200, 而是-190. 这个问题在中心单元即内插值的时候没有问题, 但是在外围单元即外插值的时候明显是个错误.

是我的错, 还是ABAQUS的程序员工作太不细致? 如果是我的问题请大家告诉我, 如果是ABAQUS的问题, 那就请ABAQUS的工作人员好好反应一下这个问题, 如果这样的小错误都没有测试的话, 那这个程序中还有多少个类似的错误呢.

显示全部楼层 · 发表于 2009-8-27 23:04:22

提示: 作者被禁止或删除内容自动屏蔽

zml570449712 · 发表于 2009-8-27 23:42:14

经验证，我师兄说的不错！
“楼主可以不加初始地应力试试，计算出来的S22应该就是和云图上的应力一致的”

cdstudio · 发表于 2009-8-28 15:21:56

云图上的应力应该是重力产生的应力吧；
不是inp中加的用来抵消重力产生的位移的Y向应力啊（这个应力应该是你算出来的吧：2*10*（-10））；
楼主可以不加初始地应力试试，计算出来的S22应该就是和云图上的应力一致的 ...
gjhnha 发表于 2009-8-27 16:04

可能是我把问题没有描述清楚, 这个模型加不加初始地应力对S22是没有影响的, 只影响了U2而已, 所以是不是加了初始地应力都无所谓, 关键是S22这条曲线.

如果您计算的这条曲线是没有加初始应力的, 那正是我所说的问题所在, S22在线弹性状态下只受重力应该是只与深度成正比, 即S22=-gamma*z, 在顶部z=0应该为0, 底部z=10应该为-20*10=-200, 但ABAQUS经过插值将高斯点值转换到整个模型得到的最大最小值却是-10-- -190, 这才是我所说的问题

显示全部楼层 · 发表于 2009-8-28 15:50:19

提示: 作者被禁止或删除内容自动屏蔽

cdstudio · 发表于 2009-8-28 15:53:31

本帖最后由 cdstudio 于 2009-8-28 08:55 编辑

那版主试试把网格细化之后再计算，看看是不是版主所需要的值
不明白了...
gjhnha 发表于 2009-8-28 08:50

网络细化后结果肯定会好一些, 但不会完全正确, 因为内插值是正确的, 只有最外围的单元进行外插值才出错.

显示全部楼层 · 发表于 2009-8-28 16:06:34

提示: 作者被禁止或删除内容自动屏蔽

账号		自动登录	找回密码
密码			注册

账号		自动登录	找回密码
密码			立即注册

gjhnha 该用户已被删除	发表于 2009-8-27 23:04:22 \| 显示全部楼层来自江苏镇江提示: 作者被禁止或删除内容自动屏蔽
gjhnha 该用户已被删除
	回复赞不支持使用道具举报显身卡

gjhnha 该用户已被删除	发表于 2009-8-28 15:50:19 \| 显示全部楼层来自江苏镇江提示: 作者被禁止或删除内容自动屏蔽
gjhnha 该用户已被删除
	回复赞不支持使用道具举报显身卡

gjhnha 该用户已被删除	发表于 2009-8-28 16:06:34 \| 显示全部楼层来自江苏镇江提示: 作者被禁止或删除内容自动屏蔽
gjhnha 该用户已被删除
	回复赞不支持使用道具举报显身卡

[心得交流] ABAQUS做外插值出错, 是我的问题?

本帖子中包含更多资源

相关帖子