硕士论文 双孔隧道 出问题了 就是不知道在哪里请教高手啊！！！ [复制链接]

gE^pOn  
一下是小弟编写的程序 老是出现一个警告，就是不知道问题出在哪里？请教高手知道。。。。急！！！！ }><[6Uz%  
IqepR >5t  
/clear f2M*]{N  
/PREP7   xz.M'az\  
ET,1,SOLID45          !设置岩土,注浆单元类型 0T(+
z)Ki  
ET,2,MESH200,6        !网格辅助单元 3>MILEY^  
-z-yk~F  
MP,EX,1,5.03e6        !地表层土弹性模量 ;&}z L.!jo  
MP,PRXY,1,0.35        !地表层土泊松比 KDP47A  
TB,DP,1,,,   Q}<QE:-&E  
TBMODIF,1,1,2.12e4    !内聚力 k
kb+qo  
TBMODIF,1,2,30        !内摩擦角 b$*2bSdv0<  
MP,DENS,1,1940        !地表层土密度 ,#GB  
H-u SdT  
MP,EX,2,4.97e6        !2层参数 2T(+VeM
Q=  
MP,PRXY,2,0.32 +Q);t,  
TB,DP,2,,,   (5th
  
TBMODIF,1,1,1.2e4   f.aa@>  
TBMODIF,1,2,23   H7Z`aQC  
MP,DENS,2,1710 EN/t5d  
' i5KRFy-  
MP,EX,3,7.53e6        !3层参数 {=Jo!t;f  
MP,PRXY,3,,0.2   T!41[vm(  
TB,DP,3,,,   ~QPTs1Vk8  
TBMODIF,1,1,1.0e4   C$C>RYE?.  
TBMODIF,1,2,21   [Y, L=p  
MP,DENS,3,1900 x6!Q''f7  
:wg=H  
MP,EX,4,3.8e6         !4层参数 0#uB[N  
MP,PRXY,4,0.3   LA2/<:  
TB,DP,4,,,   ,~1k:>njY~  
TBMODIF,1,1,2.1e4   +U_1B%e(%  
TBMODIF,1,2,33   6P>}7R}  
MP,DENS,4,1890 P*|=Z>%[0  
5=#d#dDc  
MP,EX,5,4.3e6        !5层参数 ohe[rV>EX  
MP,PRXY,5,0.31 W+"^!p|  
TB,DP,5,,,   "sDs[Lcq  
TBMODIF,1,1,2.1e4   TKGaGMx6@  
TBMODIF,1,2,30  
~@-r  
MP,DENS,5,1870 :KXI@)M  
;]>)6  
MP,EX,6,4.7e6          !6层参数 }KIS_krs  
MP,PRXY,6,,0.27 eu_ZsseZ  
TB,DP,6,,,   kRX?o'U~C  
TBMODIF,1,1,2.8e7   j} ^3v #  
TBMODIF,1,2,31   >D~8iuy]8.  
MP,DENS,6,,1880 h2Th)&Fb>  
[ @"6:tTU  
MP,EX,7,2.94e10        !管片衬砌参数 $>yfu=]?  
MP,PRXY,7,0.2   "cBqZzkk9j  
MP,DENS,7,2500 @b^$h:H  
.{`C>/"}  
MP,EX,8,2.3e6          !注浆层参数 VX8CEO  
MP,PRXY,8,0.25 ^_P?EJ,)`  
MP,DENS,8,1960 whHuV*K}  
ceCO*m~  
local,11,0,11.5,0,0 n@;B_Bt7  
local,12,1,0,0,0 E7@Gpu,o  
csys,11 hqjjd-S0  
wpcsys,-1 Y[}A4`  
cyl4,,,,,2.7,90                     );t+~YPS  
cyl4,0,0,2.7,0,3,90                 y6[le*T  
cyl4,0,0,3,0,3.2,90                 4`:POu&  
rectng,0,4.5,0,4.5   U2WHs3  
aovlap,all                       +s8R]3NJ_H  
nummrg,all                ! 合并所有元素     zYY]+)k?  
numcmp,all                ! 压缩所有元素编号 5*YvgB;  
m}RZ)c  
meshkey,1                     3Os0<1@H  
type,2                         em0Y'J  
lesize,1,,,6                       NqOX);'L0  
lesize,2,,,6 /JC1o&z_T  
lesize,3,,,6         UXpp1/d|e  
amesh,1                         a`f@&A`z  
lesize,4,,,6                     es#6/  
lesize,8,,,2 ' F9gp!s8~  
lesize,9,,,2 v;;3 K*c>  
amesh,2                      
%3#C0%{x  
lesize,5,,,6                       ~@xPoD&  
lesize,10,,,1 BQg3+w:>  
lesize,11,,,1 CZe
Zk  
amesh,3                         AgSAjBP  
lesize,12,,,3                 o}/|"(K  
lesize,13,,,3 VR&dy|5BO  
lesize,6,,,3 6G"UXNa,  
lesize,7,,,3 .f-=gZ* *  
amap,4,7,6,8,10               il!B={  
nummrg,all                       JJbd h \  
numcmp,all                     I{[Z  
save .43cI(  
                     u)o-H!a  
allsel,all             !选择所有元素 KZZY9  
arsym,x,all            !对x轴对称镜像           g\.$4N  
nummrg,all                     $m*Gu:#xm&  
numcmp,all                     6z6\xkr  
save vWeY[>oGur  
:0 n+RL*5  
rectng,4.5,36,0,4.5      !绘制右侧外界矩形       qA`@~\qh"  
nummrg,all                       gSw<C+  
numcmp,all                     `t)9u^[<(  
save aHpZhR|f$  
lesize,23,,,8,2         !划分右侧面积单元 m26YAcip}  
lesize,25,,,8,0.5 QdQ1+*/+U  
amap,9,9,8,16,17 YMK ![ q-  
Z/?{{}H+  
rectng,-11.5,-4.5,0,4.5      !绘制左侧外界矩形       YOGj__:  
nummrg,all                       Ow4(1eE_  
numcmp,all (Y7zaAG]  
lesize,26,,,3       !划分左侧面积单元 >jIn&s!}  
lesize,27,,,3   {CBb^BP  
amap,10,14,15,19,18 J9]cs?`)  
"zd_eC5  
allsel,all jTt9;?)  
arsym,y,all            !对Y轴对称镜像               m=^]93+  
nummrg,all                     -6NoEmb)\'  
numcmp,all pl1CPxSdO  
save a
%bE}  
a oU"  
rectng,-11.5,-4.5,-4.5,-20.8        !隧道下方土层 3Zd,"/RH  
rectng,-4.5,0,-4.5,-20.8 id5`YA$  
rectng,0,4.5,-4.5,-20.8         sr&W+4T  
rectng,4.5,36,-4.5,-20.8         1eshuL  
nummrg,all                     u4SL:IH{D  
numcmp,all   v7./u4S|V  
lesize,47,,,3                       v]F4o1ckk  
lesize,48,,,5,0.5 A7+Z
Y,  
lesize,49,,,5,2 j"h/v7~  
lesize,50,,,3 .yXqa"p  
lesize,51,,,5,0.5 ;zD4#7=  
lesize,52,,,3 [yQ%g;m  
lesize,53,,,5,0.5 !f5
2JQyh  
lesize,54,,,8,2 MSvZ3[5Io  
lesize,55,,,5,0.5 o#f"wQH;p  
amesh,21 Q35\w
Q#  
amesh,22 O |P<s+  
amesh,23 =%IyR  
amesh,24 *OFG3uM  
nummrg,all                     1a{r1([)  
numcmp,all
7(.Z8AO  
save 9+N%Io?!  
YuknZ&Q  
rectng,-11.5,-4.5,4.5,19.2        !隧道上方土层 ~:T@SrVI  
rectng,-4.5,0,4.5,19.2 =tt3nfZ9  
rectng,0,4.5,4.5,19.2         hd9HM5{p  
rectng,4.5,36,4.5,19.2 miQ*enZi  
nummrg,all                     4FRi=d;mP  
numcmp,all   HV/:OCK  
lesize,56,,,4,                       Po@;PR=  
lesize,57,,,3 4l%W]'  
lesize,58,,,4 TK\3mrEI  
lesize,59,,,4 MAFdJ+n#  
lesize,60,,,3 ~KMah  
lesize,61,,,4 jUA~}DVD  
lesize,62,,,3 a6DR' BC  
lesize,63,,,4 *1`X}  
lesize,64,,,8,0.5 QE[<
Y3M  
amesh,25 TMPk)N1Ka  
amesh,26 ,jis@]:  
amesh,27 )< G(C,!,.  
amesh,28 Ri"hU/H{  
nummrg,all                     |JYb4J4Ni  
numcmp,all &U]/SFY  
save Kh$"5dy  
#d\&6'O  
csys,0                  !回到原始坐标系 pIcg+~  
wpcsys,-1 #c@&mus  
allsel,all              !选择所有元素 9_:"`)]3B  
arsym,x,all             !对x轴对称镜像           f2IH2^)P  
nummrg,all                     ;
n(f?RO3X  
numcmp,all                     (wZ!OLY%}  
save >iOf3I-ATt  
z6E =%-`  
*do,i,1,20               !沿隧道轴线定义一系列关键点，点间距为3m     65~E<)UJ  
   k,,0,0,-3*i               pz['o  
*enddo 0N>K4ho6{  
l,18,70                  !根据关键点生成线 5<ruN11G  
*do,i,1,19                !通过关键点号循环增加生成线 ;i6~iLY  
   l,69+i,70+i               ;NRh0)%|o  
*enddo =+j3E<w  
*do,i,1,20                               ")uKDq  
   lesize,124+i,,,1             [ZSC]w^  
*enddo $-1ajSVJ  
save L>sLb(2\i  
0x-g0]  
type,1                     !进行拉伸，拉伸后建立的完整模型 ZJBb%d1;  
mat,1                       8Zsaq1S  
esys,11                     [//i "Nm  
extopt,aclear,1        !设置拉伸选项，即拉伸完成后删除母单元Mesh200 66eJp-5e8  
csys,12 wX]$xZ!s  
*do,i,1,20 gUx}vE-  
    asel,s,loc,z,-3*(i-1)       !选择拉伸的面 !X-\;3kC0  
    vdrag,all,,,,,,124+i        !进行拉伸 \Yr*x7!  
*enddo xo'!$a}I2  
allsel,all                   Q+ 
V<&  
nummrg,all               `Fr ,,Q81\  
numcmp,all                                     uM[|>t  
save   JnBc@qnP6  
'|ntwK*f  
csys,0                        !赋予隧道所在地层材料参数 nahq O|~  
nsel,s,loc,y,10.5,12.8         `~u=[}w  
esln,s                         ;(`bP  
mpchg,2,all                     OmWEa  
nsel,s,loc,y,6.3,10.5               l6HtZ(  
esln,s                           tf6m.  
mpchg,3,all                         ,jg #^47I  
nsel,s,loc,y,-4.5,6.3           08nh y[  
esln,s                         '1b 1N5~  
mpchg,4,all                     C][hH?.  
nsel,s,loc,y,-4.5,-9.1               Ie^Ed`  
esln,s                           N { oVz],  
mpchg,5,all   0@zJa;z'  
nsel,s,loc,y,-9.1,-20.8           P knOeW"j  
esln,s                           
=figat  
mpchg,6,all   ,#ZPg_x?1  
!allsel,all 0@"'SKq  
save I8J>>H'#A  

2w7$"N  
Finish                  !施加边界条件 J2VPOn  
/solu                       :V+rC]0  
csys,0                       cl^UFlf[  
/view,1,-1,0.5,1               1 gjaTPwY  
/replot                       6L8nw+mEK  
eplot                       /C!~v!;e  
save                         C
b{D[  
nsubst,10                     iptA#<Yj  
nlgeom,on                   ~&)\8@2  
nropt,full                     %;:![?M  
pred,on                 !启用非线性求解预测器 G=17]>U  
lnsrch,on                   [
l5jPL}6  
outres,all,last                 n?r8ZDJ'  
.euAN8L  
nsel,s,loc,x,-47.5       !选择节点       `S|F\mI~  
d,all,ux                 !施加X方向约束     x&;AY  
nsel,s,loc,x,47.5               #0<pRDXj  
d,all,ux                       Hl8-q!  
nsel,s,loc,y,-20.8       !选择节点 hTLf$_|P  
d,all,uy                 !施加Y方向约束 .2@T|WD!Ah  
nsel,s,loc,z,0           !选择节点 fL2P6N@  
nsel,a,loc,z,-60 j? BL8E'   
d,all,uz                 !施加Z方向约束 Q*#Lr4cm{  
allsel,all ON\bD?(VY  
acel,0,10,0              !施加重力 _1gNU]"  
step=0                     !用于记录计算步的参数 h8Uhr
D<:  
save   ?{TWsuP7  
Ro2V-6/  
time,1                 #1J,!seJ  
solve                       ivo><"Y(r  
save                         IwnDG;+Ap  
finish   ZzE(S  
                           lF(v<drkB  
/POST1                    !进入后处理器 G^d3$7  
PLNSOL,S,Y,0,1           !SY应力图 H I|a88   
PLNSOL,U,Y,0,1   C%yH}T\s  
      

haiahsiahas as a sa