求助：GFRP筋混凝土槽的模拟
初学anays，做了一个GFRP筋的混凝土槽做练习。长宽高分别为4360*1000*1600mm，底板厚100m，侧板80mm，GFRP筋间隔120mm，直径8mm。以为问题不大，因此没有取1/4计算，直接建了整个结构作为模型，装满密度为1300kg/m3的液体，在底部4个地方受固支约束。没有考虑自重。因为不可能到达开裂，因此没有用solid65，而是用的solid95，命令流如下：
/CLEAR
/COM,ANSYS RELEASE 10.0A1UP20060105 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：筋的应变 游上来透口气，继续找问题中。。
页:
[1]