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

x<d2/[(}mT  
一下是小弟编写的程序 老是出现一个警告，就是不知道问题出在哪里？请教高手知道。。。。急！！！！
o=C:=  
U'*~Ju  
/clear pm` f?
Py  
/PREP7   zpgRK4p,I"  
ET,1,SOLID45          !设置岩土,注浆单元类型 r?7tI0  
ET,2,MESH200,6        !网格辅助单元 ;2N: =Rv  
{?X:?M_  
MP,EX,1,5.03e6        !地表层土弹性模量 mM(Z8PA9-  
MP,PRXY,1,0.35        !地表层土泊松比 'nJ,mZx  
TB,DP,1,,,   uidoz f2}  
TBMODIF,1,1,2.12e4    !内聚力 a1#",%{I  
TBMODIF,1,2,30        !内摩擦角 n~_;tO  
MP,DENS,1,1940        !地表层土密度 *E.uqu>I  
M,w5F5  
MP,EX,2,4.97e6        !2层参数 b@X+vW{S  
MP,PRXY,2,0.32 $/J4?Wik  
TB,DP,2,,,   }-/oL+j  
TBMODIF,1,1,1.2e4   $-AG$1  
TBMODIF,1,2,23   c'md)nD2M  
MP,DENS,2,1710 ,)?!p_*@:  
H'a6] ]2  
MP,EX,3,7.53e6        !3层参数 b,hRk1  
MP,PRXY,3,,0.2   "vvv@sYxi  
TB,DP,3,,,   xlIVLv6dO  
TBMODIF,1,1,1.0e4   <~z@GMQCf  
TBMODIF,1,2,21   N.'-9hv  
MP,DENS,3,1900 NWue;u^  
Ky[s&>02  
MP,EX,4,3.8e6         !4层参数 L NS O]\  
MP,PRXY,4,0.3   N||a0&&  
TB,DP,4,,,   _D"V^4^yqu  
TBMODIF,1,1,2.1e4   lq}m0}9<  
TBMODIF,1,2,33   q8S
HFKE  
MP,DENS,4,1890 OjWg>v\v  
Fzmc#?  
MP,EX,5,4.3e6        !5层参数 :6TLT-B  
MP,PRXY,5,0.31 '/2)I8  
TB,DP,5,,,   BK]5g[   
TBMODIF,1,1,2.1e4   pddumbp  
TBMODIF,1,2,30   s+0n0C  
MP,DENS,5,1870 =[do([A  
T|k_$LH  
MP,EX,6,4.7e6          !6层参数 aE(DNeG-H  
MP,PRXY,6,,0.27 U5:5$T,C  
TB,DP,6,,,   tfIUH'Ez>  
TBMODIF,1,1,2.8e7   U2G[uDa;
 
TBMODIF,1,2,31   SiLWy=qbR  
MP,DENS,6,,1880 jI/#NCKE  
'+&!;Jj,  
MP,EX,7,2.94e10        !管片衬砌参数 br;~}GR_h
 
MP,PRXY,7,0.2   xcE2hK/+  
MP,DENS,7,2500 .C|dGE?,  
cB0"vbdO  
MP,EX,8,2.3e6          !注浆层参数 __%
){j6  
MP,PRXY,8,0.25 -J":'xCP!  
MP,DENS,8,1960 6qR5A+|;  
WU\Bs2  
local,11,0,11.5,0,0 w
5<&
b1:  
local,12,1,0,0,0 z?<Xx?Kk  
csys,11 N1fPutl$a
 
wpcsys,-1 dt5`UBvUg  
cyl4,,,,,2.7,90                     \%}w7
J;  
cyl4,0,0,2.7,0,3,90                 UX24*0`\~  
cyl4,0,0,3,0,3.2,90                 A=96N@m6  
rectng,0,4.5,0,4.5   d~qZ;uw  
aovlap,all                       +k;][VC[O  
nummrg,all                ! 合并所有元素     \ORE;pG  
numcmp,all                ! 压缩所有元素编号 #sm@|'Q%  
^^
z_[Ih  
meshkey,1                     |BEoF[1  
type,2                         x~?|bnM#3  
lesize,1,,,6                       ]kdU]}z  
lesize,2,,,6 0d/ f4  
lesize,3,,,6         IE:;`e:\D  
amesh,1                         T-F8[dd^/  
lesize,4,,,6                     LkyT4HC8n  
lesize,8,,,2 :d1Kq _\K  
lesize,9,,,2 sW]>#e  
amesh,2                       Y>2oU`ly,  
lesize,5,,,6                       +?'a2pUS  
lesize,10,,,1 QCJf   
lesize,11,,,1 dnzZ\t>U  
amesh,3                         "M}3T?0 O  
lesize,12,,,3                 sjy/[.4-  
lesize,13,,,3 Yij_'0vZ  
lesize,6,,,3 =xjtPmZ5X  
lesize,7,,,3 3w&Z:<  
amap,4,7,6,8,10               G?+0#?'Y  
nummrg,all                       ]u|5ZCv0  
numcmp,all                     ~P fk   
save {VE1c'E"V?  
                     h9Tst)iRi  
allsel,all             !选择所有元素 +<Y1`kV)  
arsym,x,all            !对x轴对称镜像           e'X"uH Xt.  
nummrg,all                     T s9go  
numcmp,all                     "Wg5eML0  
save ZFC&&[%-sG  
-&h<t/U  
rectng,4.5,36,0,4.5      !绘制右侧外界矩形       @rE+H 5  
nummrg,all                       /lLG
|aAe  
numcmp,all                     M~U>"kX  
save &SMM<^P.  
lesize,23,,,8,2         !划分右侧面积单元 0ky3rFSh1  
lesize,25,,,8,0.5 7*eIs2aY  
amap,9,9,8,16,17 1VA%xOURh  
_ |G') 9  
rectng,-11.5,-4.5,0,4.5      !绘制左侧外界矩形       II]-mb  
nummrg,all                       LS/ZZAN u  
numcmp,all nmw#4yHYy:  
lesize,26,,,3       !划分左侧面积单元 m$`4.>J  
lesize,27,,,3   .efbORp  
amap,10,14,15,19,18 ffy,ds_7  
7V%b!R}  
allsel,all g?rK&UTU  
arsym,y,all            !对Y轴对称镜像               ra F+Bt`  
nummrg,all                     MjW{JR)I  
numcmp,all 3ih:t'N-  
save 0`4Fa^o]h  
8;i
'dF:)  
rectng,-11.5,-4.5,-4.5,-20.8        !隧道下方土层 =zW`+++3  
rectng,-4.5,0,-4.5,-20.8 8vx#QU8E/  
rectng,0,4.5,-4.5,-20.8         @NYlVk2  
rectng,4.5,36,-4.5,-20.8         xf3;:soC  
nummrg,all                     "lA8CA  
numcmp,all   jwp?eL!7  
lesize,47,,,3                       Zt \3y  
lesize,48,,,5,0.5 Dco3`4pl  
lesize,49,,,5,2 Y;=GM:*H  
lesize,50,,,3 i4<n#]1!t  
lesize,51,,,5,0.5 5Z>+NKQ  
lesize,52,,,3 !-Uq#Ea0/  
lesize,53,,,5,0.5 ZMEYF!jN  
lesize,54,,,8,2 G&q@B`I  
lesize,55,,,5,0.5 ,8.zbr  
amesh,21 :gM_v?sy  
amesh,22 I:UN2`*#  
amesh,23 ts&sr  
amesh,24 *5mJA -[B+  
nummrg,all                     9w<k1j  
numcmp,all T5eJIc3a"  
save PNpH)'C|
 
^S:I38gR#q  
rectng,-11.5,-4.5,4.5,19.2        !隧道上方土层 &UQP9wS4v  
rectng,-4.5,0,4.5,19.2 QSx4M  
rectng,0,4.5,4.5,19.2         g$U7bCHG  
rectng,4.5,36,4.5,19.2 !JQ'~#jKN  
nummrg,all                     ua!RwSo  
numcmp,all   chur(@Af  
lesize,56,,,4,                       EouI S2e;a  
lesize,57,,,3 R:y u  
lesize,58,,,4 }F-,PSH Ml  
lesize,59,,,4 Q"k #eEA  
lesize,60,,,3 TOsHb+Uv  
lesize,61,,,4 b\vL^\bX8  
lesize,62,,,3 ]RuH6d2d|  
lesize,63,,,4 mW)C=X%  
lesize,64,,,8,0.5 yK>s]65&  
amesh,25 |!cM_&  
amesh,26 >mMmc!u>G  
amesh,27 eC='[W<a.  
amesh,28 V9;O1  
nummrg,all                     r)gtx!bx  
numcmp,all +7Qj%x\  
save uA%cie  
XZ4H(Cj  
csys,0                  !回到原始坐标系 08z?
i  
wpcsys,-1 ~>uu1[/  
allsel,all              !选择所有元素 `08}y*E  
arsym,x,all             !对x轴对称镜像           i9^m;Y)^I  
nummrg,all                     /h,-J8[  
numcmp,all                     a/Cc.s   
save 2NF#mWZ(s  
pnGDM)H7  
*do,i,1,20               !沿隧道轴线定义一系列关键点，点间距为3m     es1'z.UJ  
   k,,0,0,-3*i               Y'?{yx{  
*enddo m }\L i]  
l,18,70                  !根据关键点生成线 K7},X01^  
*do,i,1,19                !通过关键点号循环增加生成线 MC_i"P6a  
   l,69+i,70+i               |?LUt@r;  
*enddo eY\!}) 5  
*do,i,1,20                               VrKFpFd  
   lesize,124+i,,,1             /mQ9}E4X  
*enddo YR.f`-<Z  
save s;,ulME  
l8 XY  
type,1                     !进行拉伸，拉伸后建立的完整模型 YH3[Jvzf4  
mat,1                       C
TZ#QiNP  
esys,11                     "b0!h6$!H  
extopt,aclear,1        !设置拉伸选项，即拉伸完成后删除母单元Mesh200 to#T+d.(v  
csys,12 ]m\:XhI*<  
*do,i,1,20 nXW]9zC"/  
    asel,s,loc,z,-3*(i-1)       !选择拉伸的面 @rlL'|&X*  
    vdrag,all,,,,,,124+i        !进行拉伸 #{~3bgY  
*enddo \GCT3$  
allsel,all                   gcF V$  
nummrg,all               72sBx3 ;  
numcmp,all                                     .~%,eF;l$  
save   U=N]XwjVK<  
*40Z}1ng  
csys,0                        !赋予隧道所在地层材料参数 sDS0cc6e  
nsel,s,loc,y,10.5,12.8         KV|ywcGhT  
esln,s                         sf,9Ym  
mpchg,2,all                     d[&Ah~,  
nsel,s,loc,y,6.3,10.5               $uUJV%
EX  
esln,s                           i&Me7=~  
mpchg,3,all                         }-PV%MNud  
nsel,s,loc,y,-4.5,6.3           =UV=F/Af^  
esln,s                         $ItPUYi";  
mpchg,4,all                     8O.5ML{  
nsel,s,loc,y,-4.5,-9.1               oN[#C>#(  
esln,s                           `cqZ;(^  
mpchg,5,all   :Oc&{z?q  
nsel,s,loc,y,-9.1,-20.8           J1d|L|M  
esln,s                           ?>iZ){0,  
mpchg,6,all   GD<pqm`vVY  
!allsel,all R]y9>5 'U  
save *h~(LH"tN  
8]&lUMaqVZ  
Finish                  !施加边界条件 VMW<?V 2Z  
/solu                       98!H$6k  
csys,0                       u4~(0  
/view,1,-1,0.5,1               `$>cQwB,D  
/replot                       nE"0?VNW$  
eplot                       k(wJ6pc  
save                         M7gM#bv>L  
nsubst,10                     Dl_SEf6b  
nlgeom,on                   wb6$R};?  
nropt,full                     |dqv
v  
pred,on                 !启用非线性求解预测器 9|2LuHQu+  
lnsrch,on                   1A{iUddR  
outres,all,last                 ~c'R7E&Bfa  
GI/4<J\  
nsel,s,loc,x,-47.5       !选择节点       eQsoZQA1  
d,all,ux                 !施加X方向约束     K@@Jt  
nsel,s,loc,x,47.5               }J] P`v  
d,all,ux                       J:2Su1"ODh  
nsel,s,loc,y,-20.8       !选择节点 XaYgl&x'!x  
d,all,uy                 !施加Y方向约束 nEh^{6  
nsel,s,loc,z,0           !选择节点 <Q?_],ip  
nsel,a,loc,z,-60 baib_-$  
d,all,uz                 !施加Z方向约束 .GuZV'  
allsel,all KDYyLkI dr  
acel,0,10,0              !施加重力 g&L $5  
step=0                     !用于记录计算步的参数 C72btS  
save   {ud^+I&  
P"k,[ZQ  
time,1                 2"B3Q:0he|  
solve                       lPn&,\9@~  
save                         ?v Z5 ^k  
finish   V5]:^=  
                           Cz#3W8jV  
/POST1                    !进入后处理器 6EkD(w  
PLNSOL,S,Y,0,1           !SY应力图 M5l*D'GE]  
PLNSOL,U,Y,0,1   -
P$E)5?^  
      

haiahsiahas as a sa