【讨论】一个例子，混凝土高手请指点？不胜感激

lisagao

经过别人的指点，我把我短肢剪力墙的命令流改为如下，由于我老师要我最近用ANSYS分析出来，我本人非常急，所以想请教你门一下：但是又出现了一些新的问题，我试了很多方法都不行，
主要问题如下：
1．  我第一步荷载是用的力收敛准则cnvtol,f,,0.05,2,0.01，第二步是用的位移收敛准则cnvtol,u,,0.05,0,0.01,但我看了一下输出窗口，在每一个荷载步它都必须满足我所定义的二个收敛准则。为什么不是在第一荷载步满足力收敛准则，在第二荷载步满足位移收敛准则，是不是我的命令流有问题，请指教.
2．  第一荷载步不管是力收敛准则，还是位移收敛准则，都很好收敛。第二荷载步就是不收敛，当我用力的收敛准则时，大概算到1.5秒时就弹出如下警告错误窗口
Equivalent plastic strain increment has exceed the specified limit value, since the time increment has reached the minimum value, the analysis is stopped.
后来我有把plastic strain 有0.15改为0.5，结果还是不行，不知是不是我混凝土的本构关系出了错。请指教，谢谢。
3．  当我第二荷载步用位移的收敛准则时，大概算到1.5秒时就弹出如下警告错误窗口
the value of uz at node 892 is 6456554.56,It is greater than the current limit of 1000000, This generally indicates rigid body motion as a result of an unconstrained model,verify that your model is properly constrained.
当我按了 proceed时，又出现一警告错误
Rigid body motion can also occur when net section yielding has occurred resulting in large displacement for small increments of load or when buckling has occurred ,you can plot the time history curve for node 1959 in the uz direction to check for stiffness (slope of the curve ) approching zero.
我真不知道该怎么办。请指教，谢谢。
4  我看了一下江见鲸的用ANSYS SOLID65 分析混凝土，其中有一句话不是很明白，材料号X3:定义钢筋材料编号，体积配筋率X3，这是什么意思，为什么要乘以3？
  关于下降段的问题
下降段很容易导致计算不收敛，有时为了计算的收敛要避免设置下降段，采用RUSH模型。我不知道怎么用RUSH模型。
利用最大压应变准则来判断混凝土是否破坏。我不知在ANSYS里怎样实现。
5 我一直都不知道在ANSYS里到底用什么来本构模拟反复荷载下的混凝土，江见鲸的用ANSYS SOLID65 分析混凝土中说DRUCKER-PRAGER模型，但我不知道怎样用混凝土的参数去计算Cohesion value(内聚力)，the angle of internal friction(内摩檫角)，the dilatancy angle(膨胀角度)
  
本题为一短肢剪力墙结构体系，尺寸如图，在每个墙肢顶部作用有有不变的竖向荷载，每个墙肢作用的竖向荷载为245KN，在右边墙肢作用有257KN的水平荷载，本题采用整体式配筋。
/units,si  
ET,1,solid65  
keyopt,1,7,1  
定义混凝土的材料特性
mp,ex,1,2.7e10  
mp,prxy,1,0.2  
mp,dens,1,2.5e3  
定义混凝土的本构关系（不知道这种本构关系对吗）
TB,KINH,1,,8  
tbpt,,1.0e-4,2.7e6  
tbpt,,3.0e-4,7.3e6  
tbpt,,6.0e-4,12.4e6  
tbpt,,8.0e-4,14.4e6  
tbpt,,10e-4,15.9e6  
tbpt,,12e-4,16.2e6  
tbpt,,14e-4,16.3e6  
tbpt,,30e-4,16.3e6  
定义直径为6的钢筋材料特性
mp,ex,2,2.0e11  
mp,prxy,2,0.3  
TB,BkIN,2,,2  
TBdata,,502e6,5.3e8  
定义直径为8的钢筋材料特性
mp,ex,3,2.0e11  
mp,prxy,3,0.3  
TB, BkIN,3,,2  
TBdata,,536e6,5.6e8  
定义直径为10的钢筋材料特性
mp,ex,4,2.0e11  
mp,prxy,4,0.3  
TB, BkIN,4,,2  
TBdata,,442.5e6,7.4e8  
定义混凝土的破坏准则
Tb,concr,1,,9  
tbdata,,0.4,0.9,2.0e6,-1  
定义实常数，定义X,Y,Z方向的钢筋。
R,1,2,0.006,0,0,3,0.013 （实常数不知这样定义对吗）
Rmore,90,0,2,0.0012,0,90  
R,2,4,0.021,0,0,2,0.006  
Rmore,90,0,2,0.003,0,90
建模
Block,0,0.5,0,1.5,0,0.1  
wpoffs,0.5,1.3,0,  
Block,0,0.5,0,0.2,0,0.1  
vgen,2,1,1,1,1  
vsel,s,,,1,3,2  
vatt,1,1,1  
vsel,s,,,2,2  
vatt,1,2,1  
vsel,all  
lsel,all  
lesize,all,0.05  
vmesh,all  
nummrg,node  
nummrg,kp  
Finish  
求解控制
/solu  
在底部添加固定约束
nsel,s,loc,y,0,0  
d,all,all,0  
nsel,all  
在顶部施加竖向力
nsel,s,loc,y,1.5,1.5  
nsel,u,loc,x,0.51,0.99  
F,all,Fy,-7410（把集中力转换为节点力，是由245000/33得来的）  
cnvtol,f,,0.05,2,0.01（定义力收敛准则）  
nsubst,200,1000000,200  
Autots,on  
outres,all,all  
lnsrch,1  
ncnv,2  
neqit,200  
pred,on  
eqslv,sparse  
NROPT,FULL,,  
KBC,0  
lswrite,1
定义载荷步2，施加侧向荷载
nsel,s,loc,x,1.5,1.5  
nsel,v,loc,y,1.36,1.46  
F,all,Fx,-28556（把集中力转换为节点力，是由257000/9得来的）  
nsubst,400,1000000,400
cnvtol,u,,0.05,0,0.01(定义位移收敛准则)
lswrite,2
lssolve,1,2

xgam

:?)"我看了一下江见鲸的用ANSYS SOLID65 分析混凝土，其中有一句话不是很明白，材料号X3:定义钢筋材料编号，体积配筋率X3，这是什么意思，为什么要乘以3？ "
  
他的意思是说有三个方向要设，不是每个方向乘3！

xgam

:?) "定义混凝土的本构关系（不知道这种本构关系对吗）  
TB,KINH,1,,8  
我看了一下帮助，好像solid65不能用KINH只能用MKIN，你再看看！
另，你的弹模及屈服应力是根据什么取的？与规范为何差那么大？

zxh198199

Equivalent plastic strain increment has exceed the specified limit value, since the time increment has reached the minimum value, the analysis is stopped.后来我有把plastic strain 有0.15改为0.5，结果还是不行，
我跟你一样，我做的带圆孔的平板，当力加大了就出现你这个问题，有哪位大哥帮我分析这个问题

FengstonE

3．  当我第二荷载步用位移的收敛准则时，大概算到1.5秒时就弹出如下警告错误窗口
the value of uz at node 892 is 6456554.56,It is greater than the current limit of 1000000, This generally indicates rigid body motion as a result of an unconstrained model,verify that your model is properly constrained.
当我按了 proceed时，又出现一警告错误
Rigid body motion can also occur when net section yielding has occurred resulting in large displacement for small increments of load or when buckling has occurred ,you can plot the time history curve for node 1959 in the uz direction to check for stiffness (slope of the curve ) approching zero.

＃＃＃＃＃＃＃＃＃
这个应该是由于约束不当造成的吧，查查约束先！

scheng3217

我也是这个问题，算了一会，马上就出现这个问题
The value of uz at node 8894 is 21722494.6.it is greater than the current limit of 1000000. this generally indicates rigid body motion as a result of an unconstrained model. Verify that your model is properly constrained.

Rigid body motion can also occur when net section yielding has occurred resulting in large displacements for small increments of load or when buckling has occurred. You can plot the time history curve for node 8894 in the uz direction to check for stiffness(slop of the curve)approaching zero.

我一直算不出来，你解决了么！约束我看了没有问题，我试了好多次，换成四周全约束也不行。

[ 本帖最后由 scheng3217 于 2007-7-3 17:12 编辑 ]

zhdew

遇到同样的问题 期待解决……

liuxh2009

同样的问题期待解决

玩命儿到底

应该是材料模型不对~~
建议比对下《ANSYS工程结构数值分析》（王新敏）p486的例子

玩命儿到底

245000N=24.5吨，，荷载太大了？