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

10xza=a  
一下是小弟编写的程序 老是出现一个警告，就是不知道问题出在哪里？请教高手知道。。。。急！！！！ }WhRJr`a  
B$qTH5)W  
/clear *c 9S.  
/PREP7   z[@i=avPG  
ET,1,SOLID45          !设置岩土,注浆单元类型 p`EgMzVO,  
ET,2,MESH200,6        !网格辅助单元 hkdF  
lOM8%{.'_x  
MP,EX,1,5.03e6        !地表层土弹性模量 %cBOi_}}~  
MP,PRXY,1,0.35        !地表层土泊松比 9bhubx\^/  
TB,DP,1,,,   2A:&Cqo  
TBMODIF,1,1,2.12e4    !内聚力 @$iZ9x6t  
TBMODIF,1,2,30        !内摩擦角 w|Ry)[  
MP,DENS,1,1940        !地表层土密度 .PV(MV  
_%Yi^^  
MP,EX,2,4.97e6        !2层参数 5A<}*T  
MP,PRXY,2,0.32 vK`HgRQ(C  
TB,DP,2,,,   *Ms&WYN-  
TBMODIF,1,1,1.2e4   +o):grWvQ  
TBMODIF,1,2,23   p5% %k-  
MP,DENS,2,1710 o1"N
{Eu  
:@a0h
  
MP,EX,3,7.53e6        !3层参数 9tb-;|  
MP,PRXY,3,,0.2   e7bMK<:r  
TB,DP,3,,,   pe<T"[X  
TBMODIF,1,1,1.0e4   G"tlJ7$myQ  
TBMODIF,1,2,21   6nREuT'k  
MP,DENS,3,1900 n`@dk_%yI  
qHtIjtt[q  
MP,EX,4,3.8e6         !4层参数 A\W)uwyN  
MP,PRXY,4,0.3   o';/$xrH  
TB,DP,4,,,   k"C
'8<T)'  
TBMODIF,1,1,2.1e4   zG8g}FrzG;  
TBMODIF,1,2,33   >#'?}@FWQN  
MP,DENS,4,1890 E m^Dg9  
Z`86YYGK  
MP,EX,5,4.3e6        !5层参数 X;2LK!x;y  
MP,PRXY,5,0.31  UO#`Ak  
TB,DP,5,,,   gFrNk Uqp  
TBMODIF,1,1,2.1e4   +95v=[t#Ut  
TBMODIF,1,2,30   :Ocw+X3  
MP,DENS,5,1870 }}ic{931  
ekU%^R<  
MP,EX,6,4.7e6          !6层参数 $d!Vxm  
MP,PRXY,6,,0.27 FvYciU!  
TB,DP,6,,,   VIN0kRQ#  
TBMODIF,1,1,2.8e7   KV$4}{  
TBMODIF,1,2,31   t1"#L_<e  
MP,DENS,6,,1880 ^sFO[cYo  
eB=&(ZT  
MP,EX,7,2.94e10        !管片衬砌参数 E/*&'Osq  
MP,PRXY,7,0.2   (>f`>6 V  
MP,DENS,7,2500 Hb::;[bm:  
2ZEGE+0  
MP,EX,8,2.3e6          !注浆层参数 R&&&RI3{  
MP,PRXY,8,0.25 HZ2W`wo  
MP,DENS,8,1960 8D,*_p  
`'sD(e  
local,11,0,11.5,0,0 9j
X_Eoxy
 
local,12,1,0,0,0 #^Pab^Y3r-  
csys,11 EApbaS}Up  
wpcsys,-1 ,W|-?b?   
cyl4,,,,,2.7,90                     |FM*1Q[1  
cyl4,0,0,2.7,0,3,90                 "}`)s_rt  
cyl4,0,0,3,0,3.2,90                 T_=WX_h $  
rectng,0,4.5,0,4.5   +xU({/  
aovlap,all                       I[v`)T'_{  
nummrg,all                ! 合并所有元素     i%i/>;DF  
numcmp,all                ! 压缩所有元素编号 !U>WAD9  
X(X[v]  
meshkey,1                     rQ_@q_B.  
type,2                         `]XI Q\ *  
lesize,1,,,6                       Uc'}y!R  
lesize,2,,,6 R3l{.{3p2  
lesize,3,,,6         Bf;_~1+vLG  
amesh,1                         1' v!~*af  
lesize,4,,,6                     :v;U7  
lesize,8,,,2 jBOl:l,+  
lesize,9,,,2 -fmJkI  
amesh,2                       F.5fasdX'  
lesize,5,,,6                       A;RV~!xx  
lesize,10,,,1 9J?s:"j  
lesize,11,,,1 4C%pKV  
amesh,3                         )%MC*Z:^  
lesize,12,,,3                 0kOwA%m  
lesize,13,,,3 c>{QTI:]  
lesize,6,,,3 ~'KqiUY  
lesize,7,,,3 L/exR6M7  
amap,4,7,6,8,10               A$5
M.  
nummrg,all                       <?5|(Q"@:  
numcmp,all                     B5|\<CF  
save p^|l ',e  
                     OM&GypP6&  
allsel,all             !选择所有元素 wN4#j}C  
arsym,x,all            !对x轴对称镜像           \rzMgR$/rj  
nummrg,all                     ?\,^>4x?  
numcmp,all                     z<vh8dNl  
save ]kTxVe  
Wu/#}Bw#  
rectng,4.5,36,0,4.5      !绘制右侧外界矩形       pC>h"Hy  
nummrg,all                       |oXd4  
numcmp,all                     I*3}e
rT  
save )j>U4a  
lesize,23,,,8,2         !划分右侧面积单元 pF8 #H~  
lesize,25,,,8,0.5 a~$Y;C_#<  
amap,9,9,8,16,17 U>f'j;5  
]{\ttb%GX  
rectng,-11.5,-4.5,0,4.5      !绘制左侧外界矩形       ?V^7`3F  
nummrg,all                       e@ZM&iR  
numcmp,all pLtw|S'4  
lesize,26,,,3       !划分左侧面积单元 ~BmA!BZV`  
lesize,27,,,3   l9&k!kF`  
amap,10,14,15,19,18 O.(2  
A+0-pF2D  
allsel,all !EhKg)y=  
arsym,y,all            !对Y轴对称镜像               X&?s:A  
nummrg,all                     OpEH4X.Z  
numcmp,all f(Of+>   
save J4"Fj, FS  
! I0xq"  
rectng,-11.5,-4.5,-4.5,-20.8        !隧道下方土层 Jq'8"  
rectng,-4.5,0,-4.5,-20.8 Z%qtAPd  
rectng,0,4.5,-4.5,-20.8         *b. >  
rectng,4.5,36,-4.5,-20.8         I1U2wD  
nummrg,all                     i9DD)Y<  
numcmp,all   L!If~6oD(  
lesize,47,,,3                       :]@c%~~!&  
lesize,48,,,5,0.5 *kl  :/#  
lesize,49,,,5,2 .PUp3X-  
lesize,50,,,3 x1.S+:  
lesize,51,,,5,0.5 r@{TN6U  
lesize,52,,,3 7PX`kI  
lesize,53,,,5,0.5 +SH{`7r  
lesize,54,,,8,2 Lz'VQO1U=  
lesize,55,,,5,0.5 WQ.0}n}d  
amesh,21 *{
DTxEy  
amesh,22 VWvSt C  
amesh,23 GT&}Burl/n  
amesh,24 GBbnR:hM  
nummrg,all                     :-46
"bP.  
numcmp,all tbd=A]B-  
save "0P`=n  
9qB0F_xl  
rectng,-11.5,-4.5,4.5,19.2        !隧道上方土层 }oRBQP^&K  
rectng,-4.5,0,4.5,19.2 eyCZ[SC  
rectng,0,4.5,4.5,19.2         Icnhet4  
rectng,4.5,36,4.5,19.2 No'?8+i  
nummrg,all                     =1l6(pJ  
numcmp,all   $Iwvecn?I  
lesize,56,,,4,                       H4
w\e#|  
lesize,57,,,3 20;9XJmjl  
lesize,58,,,4 @Nek;xJ  
lesize,59,,,4 ZD
]1C~)  
lesize,60,,,3  8s0+6{vW  
lesize,61,,,4 JL1ajlm~  
lesize,62,,,3 x;n3 Zr;(  
lesize,63,,,4 z~ cW,  
lesize,64,,,8,0.5 }X&rJV  
amesh,25 69g{oo  
amesh,26  x]+PWk  
amesh,27 #X}HF$t{=  
amesh,28 m+TAaK  
nummrg,all                     2*|T)OA`m,  
numcmp,all /,Xl8<~#  
save ]ordqulq1  
hJb2y`,q  
csys,0                  !回到原始坐标系 .A2$C|a*  
wpcsys,-1 xh;V4zK@`  
allsel,all              !选择所有元素 /<J(\;Jr6  
arsym,x,all             !对x轴对称镜像           k6(7G@@}  
nummrg,all                     pH1!6X  
numcmp,all                     m9 h '!X<  
save VfX^iG r  
m@td[^O-  
*do,i,1,20               !沿隧道轴线定义一系列关键点，点间距为3m    
$9)os7H7  
   k,,0,0,-3*i               CL=%eSsuD  
*enddo %71i&T F  
l,18,70                  !根据关键点生成线 utlr|m Xc  
*do,i,1,19                !通过关键点号循环增加生成线 .uuhoqG0  
   l,69+i,70+i               pHV^Kv#  
*enddo a+w2cN'  
*do,i,1,20                               /Z:
j:l  
   lesize,124+i,,,1             sdFHr4  
*enddo s2;~FK#/  
save p['RV  
"o5]:]h)  
type,1                     !进行拉伸，拉伸后建立的完整模型 o*cu-j3  
mat,1                       "rv~I_zl  
esys,11                     |N.2iN:  
extopt,aclear,1        !设置拉伸选项，即拉伸完成后删除母单元Mesh200 8 #_pkVQw:  
csys,12 A3jxjQ  
*do,i,1,20 AkhG~L  
    asel,s,loc,z,-3*(i-1)       !选择拉伸的面 Bn}woyJdx  
    vdrag,all,,,,,,124+i        !进行拉伸 AI|vL4*Xd  
*enddo "tJ+v*E  
allsel,all                   4&iQ
o'  
nummrg,all               Nu5|tf9%A  
numcmp,all                                     oz)[-  
save   {w@qFE'b  
 :Sq]|)  
csys,0                        !赋予隧道所在地层材料参数 2|nm> 4  
nsel,s,loc,y,10.5,12.8         ^P)f]GQx  
esln,s                         LP2~UVq  
mpchg,2,all                     GEc6;uz<  
nsel,s,loc,y,6.3,10.5               Q+:y  
esln,s                           [&)]-2w2  
mpchg,3,all                         pSs*Z6c)@  
nsel,s,loc,y,-4.5,6.3           ]3KeAJ  
esln,s                         3$|/7(M&DA  
mpchg,4,all                     ci0A!wWD  
nsel,s,loc,y,-4.5,-9.1               FK%b@/7
s~  
esln,s                           ;tOsA
#  
mpchg,5,all   c_J9CKqc  
nsel,s,loc,y,-9.1,-20.8           @\K[WqF$$q  
esln,s                           "Xq_N4  
mpchg,6,all   dWAt#x
II  
!allsel,all @X==[gQ  
save UP
~28%>X  
JJ: ku&Mb  
Finish                  !施加边界条件 YQBLbtn6(  
/solu                       1"h"(dA  
csys,0                       k<gH*=uXY'  
/view,1,-1,0.5,1               P}R:o  
/replot                       amBg<P`'_  
eplot                       d0El2Ct8  
save                         gabfb#  
nsubst,10                     ?nU<cxh  
nlgeom,on                   +VIA@`4  
nropt,full                     o) )` "^  
pred,on                 !启用非线性求解预测器 XT7m3M  
lnsrch,on                   #<{v~sVp&  
outres,all,last                 {6i|"5_j  
C6!F6Stn]g  
nsel,s,loc,x,-47.5       !选择节点       )Z7Vm2a  
d,all,ux                 !施加X方向约束     ;mQj2Bwr  
nsel,s,loc,x,47.5               ?cU,%<r  
d,all,ux                       d:kB Zrq  
nsel,s,loc,y,-20.8       !选择节点 NL2
1se  
d,all,uy                 !施加Y方向约束 q lc@$  
nsel,s,loc,z,0           !选择节点 a4Q@sn;]  
nsel,a,loc,z,-60 SUS=sR/N  
d,all,uz                 !施加Z方向约束 Qqm?%7A1  
allsel,all 9+ 'i(q z  
acel,0,10,0              !施加重力 fdl.3~.C  
step=0                     !用于记录计算步的参数 T{dQ4 c  
save   e~)[
I!n  
m2]N%Y  
time,1                 _w^,j"  
solve                       (W1$+X  
save                         q*I*B1p[m  
finish   4X
kI? l  
                           gOF^?M11x  
/POST1                    !进入后处理器 wgl<JO  
PLNSOL,S,Y,0,1           !SY应力图 QiK>]xJ'  
PLNSOL,U,Y,0,1   t"=5MaQk-  
      

haiahsiahas as a sa