2楼
ansys计算\3维\不考虑砼的拉烈的代码
/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,0.4e6,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
D,ALL,ALL,0
NSEL,S,LOC,x,0,0.01
NSEL,A,LOC,x,39.9,40
D,ALL,UX,0
!NSEL,S,LOC,Z,0.0,0.01
!NSEL,A,LOC,Z,0.19,0.2
!D,ALL,UZ,0 ,
!面荷载,压力隧洞水头116米
CSYS,11
ASEL,r,LOC,X,2.7-0.01,2.7+0.01
!asel,s,,,53,62,1,
!asel,a,,,96,105,1,
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
回复
