一个简单的算例,却会出现问题,是什么原因?(混凝土材料)
需要使用混凝土材料来进行计算了,为了先熟悉这个材料的用法,于是设计了一个最简单的情况:单轴受压。可是计算的过程中却出了错,特前来向大家请教。算例详述:
一个正方形的六面体单元,相邻的三个面受正向的连杆约束,形成了一个简支的立方体。材料设为混凝土,参数按试验的数据进行填写,泊松比设为固定值,破坏面选用Kupfer面。具体数值可见命令流文件。这个单元在Z向受均匀向下的位移作用(之所以用位移,而不是加压力,是为了能够计算到下降段)。
预计结果:
该单元仅Z向有应力,Z 向的应力应变关系和定义混凝土材料时给出的应力应变关系相符。X向和Y向无应力,应变与Z向应变相比应等于泊松比。
出错信息:
计算至一定时间后,不收敛,改为自动步长也没有效果。仔细查看了结果文件,发现除了第一步的应变、应力正确,从第二步开始,侧向应变值不正确,而且X方向和Y方向出现了应力。
我的一些想法:
如果将材料改为线弹性,计算结果正确。但对于不同的材料,单轴受压的性质不应该改变。钢材料单压试验和混凝土单压试验,由于材料不同,它们的应力应变关系会有所不同,但不应导致单压变成三轴应力状态。可是在ADINA中采用两种不同的材料,竟然出现了如此迥异的结果,令人费解。表面看来,好像是因为使用了混凝土材料而导致了结果的错误,可是我检查了材料参数,没有发现错误的地方,而且,材料参数也应该是只影响应力的大小,不应该会影响应力状态。
ADINA作用一个通用有限元软件,在算法和效率上受到好评,想来不应该连这样一个小问题也计算不了,向大家请教一下,我的问题空间在什么地方呢?多谢大家了!
[ 本帖最后由 Phil_Leng 于 2008-4-14 21:04 编辑 ] 约束不对。 谢谢楼上同学的回答。但约束方面我觉得是没有问题的,三个坐标轴方向各有一个面受到连杆作用,对于只受轴向荷载作用的空间问题正好形成一个静定简支体系。而且,使用线弹性材料时的计算结果是正确的,也能说明不是约束的问题。
[ 本帖最后由 Phil_Leng 于 2008-4-16 15:32 编辑 ]
页:
[1]