一个简单的算例，却会出现问题，是什么原因？（混凝土材料）

Phil_Leng

　　需要使用混凝土材料来进行计算了，为了先熟悉这个材料的用法，于是设计了一个最简单的情况：单轴受压。可是计算的过程中却出了错，特前来向大家请教。

　　算例详述：
　　　　一个正方形的六面体单元，相邻的三个面受正向的连杆约束，形成了一个简支的立方体。材料设为混凝土，参数按试验的数据进行填写，泊松比设为固定值，破坏面选用Kupfer面。具体数值可见命令流文件。这个单元在Z向受均匀向下的位移作用（之所以用位移，而不是加压力，是为了能够计算到下降段）。

　　预计结果：
　　　　该单元仅Z向有应力，Z 向的应力应变关系和定义混凝土材料时给出的应力应变关系相符。X向和Y向无应力，应变与Z向应变相比应等于泊松比。

　　出错信息：
　　　　计算至一定时间后，不收敛，改为自动步长也没有效果。仔细查看了结果文件，发现除了第一步的应变、应力正确，从第二步开始，侧向应变值不正确，而且X方向和Y方向出现了应力。

　　我的一些想法：
　　　　如果将材料改为线弹性，计算结果正确。但对于不同的材料，单轴受压的性质不应该改变。钢材料单压试验和混凝土单压试验，由于材料不同，它们的应力应变关系会有所不同，但不应导致单压变成三轴应力状态。可是在ADINA中采用两种不同的材料，竟然出现了如此迥异的结果，令人费解。表面看来，好像是因为使用了混凝土材料而导致了结果的错误，可是我检查了材料参数，没有发现错误的地方，而且，材料参数也应该是只影响应力的大小，不应该会影响应力状态。

　　ADINA作用一个通用有限元软件，在算法和效率上受到好评，想来不应该连这样一个小问题也计算不了，向大家请教一下，我的问题空间在什么地方呢？多谢大家了！

[ 本帖最后由 Phil_Leng 于 2008-4-14 21:04 编辑 ]

geoleo

约束不对。

Phil_Leng

谢谢楼上同学的回答。但约束方面我觉得是没有问题的，三个坐标轴方向各有一个面受到连杆作用，对于只受轴向荷载作用的空间问题正好形成一个静定简支体系。而且，使用线弹性材料时的计算结果是正确的，也能说明不是约束的问题。

[ 本帖最后由 Phil_Leng 于 2008-4-16 15:32 编辑 ]