声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1429|回复: 1

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

[复制链接]
发表于 2006-7-22 21:07 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
小弟的毕业论文是用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号线所在的面没有重合面。不知道这个问题如何解决。
回复
分享到:

使用道具 举报

发表于 2006-7-23 16:59 | 显示全部楼层
关注中......
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-26 10:19 , Processed in 0.063849 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表