马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
x
我在计算隧洞中,做了一个三维模型,为了计算简便,隧洞长度取0.2米,内部用BEAM23的园杆模拟钢筋。
如果不考虑砼的破裂,即考虑成线弹性,前后两面(Z=0,Z=0.2)不加Z轴的位移约束,出现沿Z轴应力分布不对称现象,如果考虑,计算结果和2维相近;
如果仅在中心(Z=0.1)加Z轴的位移约束,计算结果也和2维相近;
为题在于考虑砼破裂后,即加入TB,CONCR属性,前后两面(Z=0,Z=0.2)加Z轴的位移约束的应力计算结果是仅在中心(Z=0.1)加Z轴的位移约束的应力计算结果的10倍
这是怎么回事啊!
帮帮忙!
怎么处理 ?
前后两面(Z=0,Z=0.2)加Z轴的位移约束的计算:
/PREP7
!绘制砼断面关键点,单位,N,KG,M,S
LOCAL,11,1,20,20,0.1
K,9,3.3,-90
K,10,4.1067,-16.94229-36.52886
K,11,3.3,-16.94229
K,12,3.3,0
K,13,3.3,45
K,14,3.3,90
K,15,3.3,135
K,16,3.3,180
K,17,3.3,180+16.94229
K,18,4.1067,16.94229+36.52886+180
K,19,3.24,-90
K,21,3.24,-16.94229
K,20,3.24,-16.94229-36.52886
K,22,3.24,0
K,23,3.24,45
K,24,3.24,90
K,25,3.24,135
K,26,3.24,180
K,27,3.24,180+16.94229
K,28,3.24,180+16.94229+36.52886
K,29,2.76,-90
K,30,2.76,-16.94229-36.52886
K,31,2.76,-16.94229
K,32,2.76,0
K,33,2.76,45
K,34,2.76,90
K,44,2.7,90
K,45,2.7,135
K,35,2.76,135
K,36,2.76,180
K,46,2.7,180
K,37,2.76,180+16.94229
K,47,2.7,180+16.94229
K,38,2.76,180+16.94229+36.52886
K,48,2.7,180+16.94229+36.52886
K,39,2.7,-90
K,40,2.7,-16.94229-36.52886
K,41,2.7,-16.94229
K,42,2.7,0
K,43,2.7,45
SAVE
K,80,0,0,
LARC,17,16,80,3.3
LARC,16,15,80,3.3
LARC,15,14,80,3.3
LARC,14,13,80,3.3
LARC,13,12,80,3.3
LARC,12,11,80,3.3
SAVE
!生成面
A,39,40,41,42,43,44,45,46,47,48
A,29,30,31,32,33,34,35,36,37,38
ASBA,2,1
A,29,30,31,32,33,34,35,36,37,38
A,19,20,21,22,23,24,25,26,27,28
ASBA,2,1
CSYS,0
A,9,10,11,12,13,14,15,16,17,18
CSYS,11
A,19,20,21,22,23,24,25,26,27,28
ASBA,1,2
!绘制围岩断面关键点
CSYS,0
K,1,0,0,0.1
K,2,20,0,0.1
K,3,40,0,0.1
K,4,40,20,0.1
K,5,40,40,0.1
K,6,20,40,0.1
K,7,0,40,0.1
K,8,0,20,.1
GPLOT
A,9,10,11,12,13,14,15,16,17,18
A,1,2,3,4,5,6,7,8
ASBA,2,1
!定义基本图形
ET,1,PLANE42
!定义材料特性
!岩石
MP,EX,1,3E9
mp,prxy,1,0.25
mp,dens,1,2760
!TB,CONC,1,1,9,
!TBDATA,,0.3,1,4e5,6.5e7,,
!砼
MP,EX,2,2.8E10
mp,prxy,2,0.167
mp,dens,2,2500
TB,CONC,2,1,9,
TBDATA,,0.3,1,1.3e6,1.25e7,,
!钢筋
ET,2,BEAM23
R,2,0.018,
KEYOPT,2,2,0
KEYOPT,2,4,1
KEYOPT,2,6,2
KEYOPT,2,10,0
MP,EX,3,20E10
mp,prxy,3,0.25
mp,dens,3,7800
!划分网格个数控制
*do,i,41,48,1
LESIZE,i,,,8
*enddo
LESIZE,1,,,4
LESIZE,2,,,8
LESIZE,3,,,8
LESIZE,4,,,8
LESIZE,5,,,8
LESIZE,6,,,4
LESIZE,7,,,8
LESIZE,8,,,8
LESIZE,9,,,4
*do,i,10,13,1
LESIZE,i,,,8
*enddo
LESIZE,14,,,4
*do,i,15,18,1
LESIZE,i,,,8
*enddo
LESIZE,19,,,4
*do,i,20,23,1
LESIZE,i,,,8
*enddo
LESIZE,24,,,4
*do,i,25,28,1
LESIZE,i,,,8
*enddo
LESIZE,29,,,4
*do,i,30,33,1
LESIZE,i,,,8
*enddo
LESIZE,34,,,4
*do,i,35,40,1
LESIZE,i,,,8
*enddo
!生成平面网格
TYPE,1
MAT,1
AMESH,6
TYPE,1
MAT,2
smrtsize
AMESH,3,5,1
!生成钢筋网格
TYPE,2
real,2
MAT,3
lMESH,17,36,1
k,81,20,20,0
k,82,20,20,0.2
l,81,80
l,80,82
!生成3维网格,345为砼,6为围岩
esize,0.1
et,3,solid65
type,3
MAT,1
vdrag,6,,,,,,50,
vdrag,6,,,,,,49,
type,3
MAT,2
vdrag,3,4,5,,,,50,
vdrag,3,4,5,,,,49,
!清除平面网格
aclear,all
save
finish
/solu
!边界定义
NSEL,S,LOC,Y,0,0.01
nplot
NSEL,a,LOC,x,0,0.01
nplot
NSEL,a,LOC,x,39.9,40
nplot
D,all, ,0, , , ,ALL, , , , ,
NSEL,S,LOC,Z,0.0,0.01
NSEL,A,LOC,Z,0.19,0.2
D,all, ,0, , , ,UZ, , , , ,
!面荷载,压力隧洞水头116米
CSYS,11
ASEL,R,LOC,X,2.7-0.01,2.7+0.01
SFa,all,,PRES,1.16e6,
ACEL,0,9.8,0,
ALLSEL,ALL
cnvtol,f,,0.05,2 !定义收敛条件,使用缺省的VALUE
nsubst,25 !定义子步数
outres,all,all !输出每一子步的结果
autots,1 !打开自动时间步控制
lnsrch,1 !打开线性搜索
ncnv,2 !如果不收敛时结束而不退出
neqit,100 !每一子步中方程的迭代次数限值
pred,on !打开预测器
time,50 !定义时间
solve !求解
!--------------进入后处理过程-----------------------
finish
/post1
pldisp,0 !绘制变形图(缺省为最后一步的)
etable,IMOMENT,SMISC,6 !定义应力表
etable,JMOMENT,SMISC,12
plnsol,s,1,0,1 !绘制Z方向的应力
PRNSOL,S,COMP
PRNSOL,S,PRIN
ETABLE,IMOMENT,SMISC,6
ETABLE,JMOMENT,SMISC,12
ETABLE,IMFORX,SMISC,1
ETABLE,IMFORY,SMISC,2
ETABLE,JMFORX,SMISC,7
ETABLE,JMFORY,SMISC,8
ETABLE,smax,nMISC,1
ETABLE,smin,nMISC,2
PRETAB,IMOMENT,JMOMENT,IMFORX,IMFORY,JMFORX,JMFORY,smax,smin
仅在中心(Z=0.1)加Z轴的位移约束:
NSEL,S,LOC,Z,0.0,0.01
NSEL,A,LOC,Z,0.19,0.2
改为
NSEL,S,LOC,Z,0.1,0.11
其余不变
|