用ansys做隧道的同志近来看看～

rongzhoutt

小弟的毕业论文是用ansys做一个隧道的通风竖井的施工过程分析。通风竖井就是一个竖着的隧道。再加上左右两个联络风道，就成

一个倒T字形。我打算只用一种元素solid45，整个建模过程由内向外扩张，分5层，1内空洞，2二次衬砌，3一次衬砌，4锚杆加固区

，5围岩
现在还在建立几何模型第一步就碰到不少问题，邀请有兴趣的同志一起探讨哈
这是我的部分命令：
FINI
/CLE
/FILNAME,'SHUJING'
/UNIT,SI
/CONFIG,NRES,1E10

!建立几何模型相关的常量!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
LS=13.5      !纵向
WD=11.5      !横向
H1=3.7       !高度
H2=4.6
!K,9,LS,WD,-H1      自由面的四个角点
!K,10,-LS,WD,-H2
!K,11,-LS,-WD,H1
!K,12,LS,-WD,H2


/PREP7

/PNUM,KP,1
/PNUM,LINE,1
/PNUM,AREA,1
/PNUM,VOLU,1

ET,1,95

!建立几何模型!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!第一层!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
LOCAL,11,0,0,0,0,0,-90,   !定义局部坐标
WPCSYS,,11              !把工作平面转到11
CYLIND,0,3.05,-3,3,0,360    !建立第一个圆柱体

K,,LS,WD,-H1
K,,-LS,WD,-H2
K,,-LS,-WD,H1
K,,LS,-WD,H2

A,9,10,11,12      !面积5
VSBA,1,5,SEPO,,DELETE
VDELE,3,,,1
ALLS
NUMMRG,ALL
NUMCMP,ALL


CYLIND,0,3.05,-6,-3,0,360    !建立第二个圆柱体
VGLUE,1,2
ALLS
NUMMRG,ALL
NUMCMP,ALL

CYLIND,0,3.05,-11,-6,0,360    !建立第三个圆柱体
VGLUE,2,3
ALLS
NUMMRG,ALL
NUMCMP,ALL

CYLIND,0,3.05,-83,-11,0,360    !建立第四个圆柱体
VGLUE,3,4
ALLS
NUMMRG,ALL
NUMCMP,ALL

CYLIND,0,3.05,-88,-83,0,360    !建立第五个圆柱体
VGLUE,4,5
ALLS
NUMMRG,ALL
NUMCMP,ALL

CYLIND,0,3.05,-187.426,-88,0,360    !建立第六个圆柱体
VGLUE,5,6
ALLS
NUMMRG,ALL
NUMCMP,ALL

LSEL,S,,,33,34,1    !底面的2条弧线
LCOMB,ALL           !合并2条弧线，L33
ALLS
NUMMRG,ALL
NUMCMP,ALL

LSEL,S,,,34,35,1    !底面的2条弧线
LCOMB,ALL           !合并2条弧线,L34
ALLS
NUMMRG,ALL
NUMCMP,ALL


BLC5,0,0,6.1,6.1
AOFFST,20,-192.426
ADELE,20,,,1
ALLS
NUMMRG,ALL
NUMCMP,ALL      !生成正方形底面20

ADELE,20,,,

LSEL,S,,,38,40,2
LDIV,ALL,0.5,,2
ALLS
NUMMRG,ALL
NUMCMP,ALL

LFILLT,39,40,0.2,0
LFILLT,39,41,0.2,0
LFILLT,37,38,0.2,0
LFILLT,37,42,0.2,0

LSEL,S,,,39,41,1
LSEL,A,,,43,44
LCOMB,ALL
ALLS
NUMMRG,ALL
NUMCMP,ALL

LSEL,S,,,37,42,1
LSEL,U,,,39
LCOMB,ALL
ALLS
NUMMRG,ALL
NUMCMP,ALL

*DO,I,33,34,1
ASKIN,I,I+4
*ENDDO
ALLS
NUMMRG,ALL
NUMCMP,ALL

*DO,I,33,33
LSEL,S,,,I,I+1
LSEL,A,,,I+4,I+5
*ENDDO

LDIV,ALL,0.5,,2
ALLS
NUMMRG,ALL
NUMCMP,ALL

!本来想直接连37,38,77,78但是不行,只好分开做再合起来！！！！！！！！！！！！！
A,29,33,34
A,30,33,34
AADD,22,23
ALLS
NUMMRG,ALL
NUMCMP,ALL


VA,17,20,21,22     !第7个体积（由圆向圆角正方形过渡）

K,35,3.05,,-207.214

L,29,35
VDRAG,22,,,,,,45
LDELE,45,,,1
ALLS
NUMMRG,ALL
NUMCMP,ALL

save,SHUJING1,DB      !!!!!保存SHUJING1!!!!!!!!!!!

!联络风道
LOCAL,12,0,-3.05,-204.885,,,,90       !联络风道圆心处设置局部坐标
LOCAL,13,0,3.05,-200.126,,,,-90
WPCSYS,,12
CYL4,,,0,,1.65,180,0
RECTNG,-1.65,1.65,-2.1,0,
!AADD,28,29
!AOVLAP,23,24,30
!VDEL,8,,,
!ADEL,23,24,,1
ALLS
NUMMRG,ALL
NUMCMP,ALL

WPCSYS,,13
CYL4,,,0,,2.1,180,0
RECTNG,-2.1,2.1,-2.3,0,
!AADD,30,31
!AOVLAP,23,24,32
ALLS
NUMMRG,ALL
NUMCMP,ALL

VDELE,8

ASEL,S,,,23,31
ASEL,U,,,27
AOVLAP,ALL
ALLS
NUMMRG,ALL
NUMCMP,ALL

ASEL,S,,,24,31,1
ADELE,ALL,,,1
ALLS
NUMMRG,ALL
NUMCMP,ALL

AL,51,52,55,56,59,60
AL, 53,54,57,58,61,62
ALLS
NUMMRG,ALL
NUMCMP,ALL        !联络风道与竖井交界面

VA,23,24,25,26,27,28,29    !这一步由面积构成体积有问题



主要有两个问题
1、输入输出有效数字的问题，ansys只能认三位有效数字，但有些数据有六位有效数字，这个问题有的网友说药用fortran，可惜我早忘的差不多了，知道的同志能否贴一哈
2、最后一行命令为什么不能执行？error是Discontinuous area set for VA command. line 37 is used 1 times for this list of areas。我仔细看过了，37号线所在的面没有重合面。不知道这个问题如何解决。

sysh320

关注中......