[原创]隧道上部开挖或填方安全评估通用程序

freddieaid

[watermark]    采用破损阶段评价标准，自己搞的，下次再搞个优化程序，斑竹能不能加分鼓励小弟一下,哈哈[/watermark]

freddieaid

自己一个半星期的心血，所以还是自个儿顶一下！

freddieaid

太让我失望了吧，再顶一下，程序包含内容：
1.交互式对话框的使用；
2.宏库文件的使用；
3.单元生死；
4.在每一工况计算完成后直接进入后处理提取数据，包括几何参数，应力应变，等效节点力，后处理也完全由程序完成；
5.计算模型可根据需要，在对话框中直接输入几何参数，即可得到相应模型网格。
我觉得这个程序已经涉及到了大部分的apdl的内容，为什么没人来捧场呢？

zhanghuaa

顶一下吧，楼主挺着急的！

xialibaren19

非常好！很多命令没用过：multipro,';start';,2    .....multipro,';end';
  *ask,.....
特别是：
et,1,E_num1
sectype,1,beam,rect,chenqi,0
secoffset,cent
secdata,1,Dep_mat1
这样也可以阿？
这个用的是beam188？
好像有et,1的说明，没有其他的阿？不懂！
最不懂得是：围岩怎么可以用beam单元 模拟？E_num1，2，3到底是什么类型？
怎么生成弯矩轴力的？

gpzhang

谢谢分享

freddieaid

下面引用由xialibaren19在 2005/04/26 02:54pm 发表的内容：
非常好！很多命令没用过：multipro,';start';,2    .....multipro,';end';
  *ask,.....
特别是：
et,1,E_num1
...

衬砌用的是beam3单元，为什么衬砌不能用梁单元模拟，还望大侠指点？
尽管梁单元与实体单元的自由度不能完全吻合，但是这样近似处理还是能满足精度要求的。毕竟在分析地层结构模型中，误差能控制在60％的范围就是不错的了。这这是我弄的一个小东西，互相学习嘛，别没事就较真，行吗?

xialibaren19

r,1,Dep_mat1,Dep_mat1/12,Dep_mat1
你是取了单宽1米计算的吧，i=bh3/12,你似乎少了个3次方阿！

xialibaren19

另外，我在计算中曾采用beam3计算衬砌同plane42比较过，相差很大的阿！
你帮我看看：
按beam3:
PREP7
!绘制砼断面关键点,单位,N,KG,M,S  
LOCAL,11,1,20,20
K,80,0,0,
CIRCLE,80,3.0,,,360,
SAVE
!生成面
A,1,2,3,4
!绘制围岩断面关键点
CSYS,0  
K,5,0,0
K,6,20,0
K,7,40,0
K,8,40,20  
K,9,40,40
K,10,20,40
K,11,0,40
K,12,0,20
GPLOT   
A,5,6,7,8,9,10,11,12
ASBA,2,1
!定义基本图形
ET,1,PLANE42
KEYOPT,1,1,0
KEYOPT,1,2,0
KEYOPT,1,3,3
KEYOPT,1,5,0
KEYOPT,1,6,0
R,1,0.2,
!定义材料特性
!岩石
MP,EX,1,3E9
mp,prxy,1,0.25
mp,dens,1,2760

!砼
MP,EX,2,2.8E10
mp,prxy,2,0.167
mp,dens,2,2500

!划分网格个数控制
*do,i,1,4,1
LESIZE,i,,,8
*enddo  
!生成平面网格
TYPE,1
real,1  
MAT,1  
AMESH,3
et,2,beam3
r,2,0.12,3.6e-3,0.6
TYPE,2
real,2  
MAT,2
lMESH,1,4,1
save
finish

/solu
边界定义
NSEL,S,LOC,Y,0,0.01
nplot   

D,all, ,0, , , ,ALL, , , , ,

lsel,s,,,1,4,1,
SFL,all,PRES,1.16e6,
ACEL,0,9.8,0,
ALLSEL,ALL
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

xialibaren19

按plane：
/PREP7  
LOCAL,11,1,20,20
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
K,2,20,0
K,3,40,0
K,4,40,20   
K,5,40,40  
K,6,20,40
K,7,0,40
K,8,0,20
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

KEYOPT,1,1,0
KEYOPT,1,2,0
KEYOPT,1,3,3
KEYOPT,1,5,0
KEYOPT,1,6,0
R,1,0.2,
R,2,0.018,
!岩石
MP,EX,1,3E9
mp,prxy,1,0.25
mp,dens,1,2760
!TB,CONC,1,1,9,  
!TBTEMP,0
!TBDATA,,0.3,1,400,65000,,
!砼
MP,EX,2,2.8E10
mp,prxy,2,0.167
mp,dens,2,2500
!TB,CONC,2,1,9,  
1TBTEMP,0
!TBDATA,,0.3,1,1300,12500,,
钢筋
ET,2,BEAM23
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

  
LESIZE,41,,,8  
LESIZE,42,,,8
LESIZE,43,,,8
LESIZE,44,,,8
LESIZE,45,,,8  
LESIZE,46,,,8  
LESIZE,47,,,8  
LESIZE,48,,,8
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  
LESIZE,10,,,8
LESIZE,11,,,8  
LESIZE,12,,,8
LESIZE,13,,,8
LESIZE,14,,,4  
LESIZE,15,,,8  
LESIZE,16,,,8
LESIZE,17,,,8  
LESIZE,18,,,8  
LESIZE,19,,,4  
LESIZE,20,,,8
LESIZE,21,,,8  
LESIZE,22,,,8  
LESIZE,23,,,8
LESIZE,24,,,4  
LESIZE,25,,,8  
LESIZE,26,,,8
LESIZE,27,,,8  
LESIZE,28,,,8  
LESIZE,29,,,4  
LESIZE,30,,,8
LESIZE,31,,,8  
LESIZE,32,,,8  
LESIZE,33,,,8  
LESIZE,34,,,4  
LESIZE,35,,,8  
LESIZE,36,,,8
LESIZE,37,,,8  
LESIZE,38,,,8  
LESIZE,39,,,8
LESIZE,40,,,8
TYPE,1
real,1  
MAT,1   
AMESH,6
TYPE,1
real,1
MAT,2
smrtsize,5  
AMESH,3,5,1
TYPE,2
real,2  
MAT,3   
lMESH,17,36,1
save
finish

/solu
NSEL,S,LOC,Y,0,0.01
nplot   
!NSEL,a,LOC,x,0,0.01
!nplot   
!NSEL,a,LOC,x,19.9,20
D,all, ,0, , , ,ALL, , , , ,
lsel,s,,,7,16,1,
SFL,all,PRES,1.16e6,
ACEL,0,9.8,0,
ALLSEL,ALL
solve
finish
/post1
pldisp,0 !绘制变形图(缺省为最后一步的)
plnsol,s,1,0,1 !绘制Z方向的应力
PRNSOL,S,COMP!输出各方向plane42单元应力
PRNSOL,S,PRIN!输出各方向plane42单元123主应力
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
!*cfopen,钢筋结果,txt,,
!*write,IMOMENT,JMOMENT,IMFORX,IMFORY,JMFORX,JMFORY,smax,smin
!(f12.6)
!*cfclos,all,

xialibaren19

因为考虑0.2米取一根钢筋，故厚度取0.2米

xialibaren19

还有一个问题，就是运行你的程序时，要求confirm your,module(y or  n)
输入y说不认识的参数
下面要求确认求解，输入y说不认识的参数
请教这是为什么

freddieaid

应输入';y';,加单引号

freddieaid

注意：
    岩土应采用drucker-prager准则，你采用的是弹性材料，这样不行

xialibaren19

谢谢！应输入';y';,加单引号
另外，我说的那个实常数的问题，你看到没有？
还有关于drucker-prager准则：DP)适用于破碎材料如土壤、岩石、混凝土，采用摩尔库伦定则，输入仅包括3个常量：
1凝聚力
2内部破碎角
3剪胀角
剪胀量（因为屈服引起的材料体积增长）采用剪胀角控制，如果和破碎角相等，流变规则相关，如果=0或小于破碎角，屈服后体积不改变流变准则不相关。
我不太清楚这些输入是在计算什么时有用？
我的计算是有压隧洞，不考虑开挖和卸载问题！
另外，你在beam3中加入concret属性是没有用的！