求助：GFRP筋混凝土槽的模拟

laokxu

初学anays，做了一个GFRP筋的混凝土槽做练习。长宽高分别为4360*1000*1600mm，底板厚100m，侧板80mm，GFRP筋间隔120mm，直径8mm。

以为问题不大，因此没有取1/4计算，直接建了整个结构作为模型，装满密度为1300kg/m3的液体，在底部4个地方受固支约束。没有考虑自重。因为不可能到达开裂，因此没有用solid65，而是用的solid95，命令流如下：
/CLEAR  
/COM,ANSYS RELEASE 10.0A1  UP20060105       10:27:03    06/04/2010
/CONFIG, NRES,2000  
！定义单元与材料性质……………………………………………………………………
/PREP7  
ET,1,SOLID95,,,,,,,1
ET, 2, LINK8
MP, EX, 1, 16000
MP, PRXY,1,0.2  

MP,EX,2,3.7E4
MP,PRXY,2,0.3   
TB,BISO,2   
TBDATA,,300,0   
*SET,PI,ACOS(-1)
R,1,0.25*PI*8*8   
！创建几何模型……………………………………………………………………
BLOCK,,1000,,1600,0,4360
BLOCK,80,920,100,1600,80,4280   
VSBV,       1,       2
!BLOCK,,100,-100,0,840,1340    ! 实际计算可以不计
!BLOCK,,100,-100,0,3120,3620
!BLOCK,900,1000,-100,0,840,1340  
!BLOCK,900,1000,-100,0,3120,3620
！切分工作平面，切出钢筋位置…………………………………………………
WPOFF,,,20  
VSBW,ALL
*DO,I,1,36   
WPOFF,,,120
VSBW,ALL
*ENDDO

WPCSYS,-1   
WPROTA,,,90
WPOFF,,,20
VSBW,ALL
*DO,I,1,8   
WPOFF,,,120
VSBW,ALL
*ENDDO

WPCSYS,-1
WPROTA,,-90
WPOFF,,,30
VSBW,ALL
*DO,I,1,13   
WPOFF,,,120
VSBW,ALL
*ENDDO  

WPCSYS,-1   
WPOFF,,,80  
VSBW,ALL
WPOFF,,,4200  
VSBW,ALL
WPCSYS,-1   
WPROTA,,-90
WPOFF,,,100  
VSBW,ALL
WPCSYS,-1  
WPROTA,,,90
WPOFF,,,80
VSBW,ALL
WPOFF,,,840  
VSBW,ALL
WPCSYS,-1  

！划分钢筋网格……………………………………………………………………
*SET,ELEMSIZ,60

LSEL,S,LOC,Z,20
LSEL,R,LOC,X,20,980
LSEL,R,LOC,Y,30,1590
LSEL,U,LOC,Y,100
LSEL,U,LOC,X,80
LSEL,U,LOC,X,920
LSEL,R,TAN1,Z
CM,GJ1,LINE
LATT,2,1,2  
LESIZE,ALL,ELEMSIZ

LSEL,S,LOC,Z,4340
LSEL,R,LOC,X,20,980
LSEL,R,LOC,Y,30,1590
LSEL,U,LOC,Y,100
LSEL,U,LOC,X,80
LSEL,U,LOC,X,920
LSEL,R,TAN1,Z
CM,GJ2,LINE
LATT,2,1,2  
LESIZE,ALL,ELEMSIZ

LSEL,S,LOC,X,20
LSEL,R,LOC,Z,20,4340
LSEL,R,LOC,Y,30,1590
LSEL,U,LOC,Y,100
LSEL,U,LOC,Z,80
LSEL,U,LOC,Z,4280
LSEL,R,TAN1,X
CM,GJ3,LINE
LATT,2,1,2  
LESIZE,ALL,ELEMSIZ

LSEL,S,LOC,X,980
LSEL,R,LOC,Z,20,4340
LSEL,R,LOC,Y,30,1590
LSEL,U,LOC,Y,100
LSEL,U,LOC,Z,80
LSEL,U,LOC,Z,4280
LSEL,R,TAN1,X
CM,GJ4,LINE
LATT,2,1,2  
LESIZE,ALL,ELEMSIZ

LSEL,S,LOC,Y,30
LSEL,R,LOC,Z,20,4340
LSEL,R,LOC,X,20,980
LSEL,U,LOC,X,80
LSEL,U,LOC,X,920
LSEL,U,LOC,Z,80
LSEL,U,LOC,Z,4280
LSEL,R,TAN1,Y
CM,GJ5,LINE
LATT,2,1,2  
LESIZE,ALL,ELEMSIZ

LSEL,ALL
CMSEL,S,GJ1  
CMSEL,A,GJ2
CMSEL,A,GJ3
CMSEL,A,GJ4
CMSEL,A,GJ5
CM,GJ,LINE  
LMESH,ALL   
LSEL,ALL
/ESHAPE,1   
EPLOT   
！划分混凝土网格……………………………………………………………………
VATT,1,,1   
MSHKEY,1
ESIZE,ELEMSIZ   
VMESH,ALL   
ALLSEL,ALL  
！施加荷载和约束……………………………………………………………………
NSEL,S,LOC,Y,0
NSEL,R,LOC,X,0,100
NSEL,R,LOC,Z,840,1340
D,ALL,ALL

NSEL,S,LOC,Y,0
NSEL,R,LOC,X,0,100
NSEL,R,LOC,Z,3120,3620
D,ALL,ALL

NSEL,S,LOC,Y,0
NSEL,R,LOC,X,900,1000
NSEL,R,LOC,Z,3120,3620
D,ALL,ALL

NSEL,S,LOC,Y,0
NSEL,R,LOC,X,900,1000
NSEL,R,LOC,Z,840,1340
D,ALL,ALL

SFGRAD,PRES, ,Y,100,-1.3E-5
ASEL,S,LOC,Z,80
ASEL,R,LOC,X,80,920
ASEL,R,LOC,Y,100,1600
SFA,ALL,1,PRES,0.0195

SFGRAD,PRES, ,Y,100,-1.3E-5
ASEL,S,LOC,Z,4280
ASEL,R,LOC,X,80,920
ASEL,R,LOC,Y,100,1600
SFA,ALL,1,PRES,0.0195

SFGRAD,PRES, ,Y,100,-1.3E-5
ASEL,S,LOC,X,80
ASEL,R,LOC,Z,80,4280
ASEL,R,LOC,Y,100,1600
SFA,ALL,1,PRES,0.0195

SFGRAD,PRES, ,Y,100,-1.3E-5
ASEL,S,LOC,X,920
ASEL,R,LOC,Z,80,4280
ASEL,R,LOC,Y,100,1600
SFA,ALL,1,PRES,0.0195

SFGRAD

ASEL,S,LOC,Y,100
ASEL,R,LOC,X,80,920
ASEL,R,LOC,Z,80,4280
SFA,ALL,1,PRES,0.0195

NSEL,S,LOC,Y,1600
NSEL,U,LOC,X,80,920  
SF,ALL,PRES,0.172

ALLSEL,ALL  
FINISH  

/SOL
ANTYPE,0
NSUBST,10  
OUTRES,ALL,ALL  
AUTOTS,ON
NEQIT,50
CNVTOL,U,,0.015
SOLVE   
FINISH

计算完后，得到的结果让我头晕。solid的单元结果还好，都是对称的，但GFRP筋的结果却有问题，出现了应力为0的地方，位移很大。我已经找原因找了3、4天了，还是没找到，希望各位大侠不要因为问题简单不屑回答呀，呵呵。不胜感激！

附：几张结果图。
Model：模型图，受4个约束
Displacement of Nodes：节点位移，部分出现很大的情况
Stress of Nodes：节点应力，有部分为0
Strain of concrete：混凝土的应变
Strain of GFRP：筋的应变

laokxu

游上来透口气，继续找问题中。。