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

Ps!~miN|>  
一下是小弟编写的程序 老是出现一个警告，就是不知道问题出在哪里？请教高手知道。。。。急！！！！ eL7\})!W  
Q-M"+HO  
/clear kX)Xo`^Ys  
/PREP7   +:&,Ts/  
ET,1,SOLID45          !设置岩土,注浆单元类型 2PrUI;J$  
ET,2,MESH200,6        !网格辅助单元 ZdjmZx%
%  
TUV&9wKXo  
MP,EX,1,5.03e6        !地表层土弹性模量  !xz0zT.  
MP,PRXY,1,0.35        !地表层土泊松比 bx&?EUx+b  
TB,DP,1,,,   ]NrA2i?  
TBMODIF,1,1,2.12e4    !内聚力 ndU<,{r  
TBMODIF,1,2,30        !内摩擦角 gnPu{-E
c*  
MP,DENS,1,1940        !地表层土密度 Y^~Dr|5%  
_9Zwg+oO[  
MP,EX,2,4.97e6        !2层参数 )k}UjU`!  
MP,PRXY,2,0.32 K~B@8az  
TB,DP,2,,,   ggn:DE"  
TBMODIF,1,1,1.2e4   I"<ACM  
TBMODIF,1,2,23   a*gzVE7W#n  
MP,DENS,2,1710 6!GO{2d"  
@3F4Lg6H|  
MP,EX,3,7.53e6        !3层参数 OcWzo#q4[  
MP,PRXY,3,,0.2   [Q5>4WY
 
TB,DP,3,,,   y6*9, CF  
TBMODIF,1,1,1.0e4   tEXY>=  
TBMODIF,1,2,21   6+hx64 =  
MP,DENS,3,1900 jRCG}
'  
`78:TU~5S  
MP,EX,4,3.8e6         !4层参数 }JePEmj  
MP,PRXY,4,0.3   L]C|&KP  
TB,DP,4,,,   i`}nv,  
TBMODIF,1,1,2.1e4   |wFfVDp  
TBMODIF,1,2,33   R8U?s/*  
MP,DENS,4,1890 sF`ELrR \  
K4\#b}P!  
MP,EX,5,4.3e6        !5层参数 ?UGA-^E1  
MP,PRXY,5,0.31 aV9QIH~  
TB,DP,5,,,   bdUe,2Yin  
TBMODIF,1,1,2.1e4   Anu`F%OzB  
TBMODIF,1,2,30   93aRWEu3  
MP,DENS,5,1870 {!`0i  
QC+K:
jL  
MP,EX,6,4.7e6          !6层参数 Z*Fxr;)d  
MP,PRXY,6,,0.27 eJ3w}"?9s  
TB,DP,6,,,   zJ2dPp~u  
TBMODIF,1,1,2.8e7   R!8qkG  
TBMODIF,1,2,31   U:n3V  
MP,DENS,6,,1880 / .ddx<  
KPcOW#.T  
MP,EX,7,2.94e10        !管片衬砌参数 !^:)zORYR  
MP,PRXY,7,0.2   /)Pf ]  
MP,DENS,7,2500 ^t{2k[@  

Nhnw'9  
MP,EX,8,2.3e6          !注浆层参数 97%S{_2m/  
MP,PRXY,8,0.25 );zLy?n  
MP,DENS,8,1960 L6-zQztn  
9+o`/lk1  
local,11,0,11.5,0,0 *$eMM*4  
local,12,1,0,0,0 .7|kxJq  
csys,11 sD[G?X  
wpcsys,-1 #o]/&T=N=  
cyl4,,,,,2.7,90                     Fuuy_+p@G  
cyl4,0,0,2.7,0,3,90                 !b0ANIp  
cyl4,0,0,3,0,3.2,90                 RZ#~^5DiO  
rectng,0,4.5,0,4.5   U)n+j}vi  
aovlap,all                       QmpP_eS >  
nummrg,all                ! 合并所有元素     xy$agt>j>  
numcmp,all                ! 压缩所有元素编号 dO{a!Ca  
KiDL]2  
meshkey,1                     quPNwNy  
type,2                         ta\CZp  
lesize,1,,,6                       GYq.!d@O  
lesize,2,,,6 ~T_4M  
lesize,3,,,6         L93&.d@m9  
amesh,1                        
<FfdOK_  
lesize,4,,,6                     ;n)9  
lesize,8,,,2 /Dl{I7W  
lesize,9,,,2 Pl}>  
amesh,2                       uGxh}'&  
lesize,5,,,6                       ksB-fOv*N  
lesize,10,,,1
gh{Z=_  
lesize,11,,,1 a2MFZe  
amesh,3                         Xj^Hy"HC^~  
lesize,12,,,3                 im6Rx=}E{  
lesize,13,,,3 '8$*gIQ8  
lesize,6,,,3 tXgsWG?v[H  
lesize,7,,,3 E~y@ue:  
amap,4,7,6,8,10               3{wmKo|_X  
nummrg,all                       O_u2V'jy9  
numcmp,all                     pW>?%ft.  
save FXi"o $N  
                     cR0OJ'w  
allsel,all             !选择所有元素 G2A^+R0\  
arsym,x,all            !对x轴对称镜像           };z[x2l^  
nummrg,all                     5#|f:M]Bo|  
numcmp,all                     &u@<0 1=  
save {xzs{)9|Y4  
I|27%i  
rectng,4.5,36,0,4.5      !绘制右侧外界矩形       yp}a&Dg  
nummrg,all                       $ MN1:ih  
numcmp,all                     hrRkam!y  
save "?<$>\@; q  
lesize,23,,,8,2         !划分右侧面积单元 2!u4nxZ.  
lesize,25,,,8,0.5 !X5~!b^*  
amap,9,9,8,16,17 kDz!v?Z2+B  
G{ 9p.Q  
rectng,-11.5,-4.5,0,4.5      !绘制左侧外界矩形       i^2yq&uT(  
nummrg,all                       ,'E+f%  
numcmp,all dF?:&oP]  
lesize,26,,,3       !划分左侧面积单元 #H;yXsR`  
lesize,27,,,3   sKvz<7pag  
amap,10,14,15,19,18 CSC sJE#4
 
,W!v0*uxp&  
allsel,all *}hx9:9\B  
arsym,y,all            !对Y轴对称镜像               >*hY1@N1  
nummrg,all                     0"*!0s~  
numcmp,all X<OOgC  
save rLU+-_  
v
e [*t`  
rectng,-11.5,-4.5,-4.5,-20.8        !隧道下方土层 Y30e7d* qr  
rectng,-4.5,0,-4.5,-20.8 GRt1]%l#$  
rectng,0,4.5,-4.5,-20.8         $i^#KZ}-WK  
rectng,4.5,36,-4.5,-20.8         ?Z]5 [  
nummrg,all                     bh9!OqK9K  
numcmp,all   |@a.dgz,  
lesize,47,,,3                       Ch~2w)HAA  
lesize,48,,,5,0.5 9F&s9(=\  
lesize,49,,,5,2 0KQDw  
lesize,50,,,3 ^HYrJr$y  
lesize,51,,,5,0.5 Vlk]  
lesize,52,,,3 yv@td+-"D  
lesize,53,,,5,0.5 gg-4ce/  
lesize,54,,,8,2 sSM^net0  
lesize,55,,,5,0.5 U0PQ[Y#\  
amesh,21 BO3#*J5S\  
amesh,22 VKjDK$  
amesh,23 |V 3AA   
amesh,24 eB$S d  
nummrg,all                     {g%F 3-  
numcmp,all l20fA-T _I  
save Aw38Tw  
Y]ZNA
R  
rectng,-11.5,-4.5,4.5,19.2        !隧道上方土层 nsRZy0@$t  
rectng,-4.5,0,4.5,19.2 ^"g # !  
rectng,0,4.5,4.5,19.2         ^+hqGu]M  
rectng,4.5,36,4.5,19.2 _wC4n }J  
nummrg,all                     U=<d
;2N#  
numcmp,all   H1alf_(_ \  
lesize,56,,,4,                       @.PVUP  
lesize,57,,,3 h]6"~ m  
lesize,58,,,4 lBbUA)z6  
lesize,59,,,4 )_vE"ryThA  
lesize,60,,,3 Z;nbnRz  
lesize,61,,,4 7 fE QD?C  
lesize,62,,,3 K|n$-WDG}  
lesize,63,,,4 Y0f"}A1  
lesize,64,,,8,0.5 U',9t  
amesh,25 vUX(h.}8  
amesh,26 zO)>(E?  
amesh,27 \ nIz5J}3  
amesh,28 YL$#6d  
nummrg,all                     nMniHB'  
numcmp,all /qYo*S_cG  
save uEK9  
o:E_k#Fi  
csys,0                  !回到原始坐标系 eq|G\XJ  
wpcsys,-1 <K$X>&Ts  
allsel,all              !选择所有元素 &g{b5x{iD  
arsym,x,all             !对x轴对称镜像           kFJ sB,2-  
nummrg,all                     Q9UBxpDV:  
numcmp,all                     errT7&@,A  
save -W^jmwM   
E7_)P>aS5  
*do,i,1,20               !沿隧道轴线定义一系列关键点，点间距为3m     Y'75DE<BC  
   k,,0,0,-3*i               : " ([i"  
*enddo "`aNNIG&  
l,18,70                  !根据关键点生成线 Vz"Ja  
*do,i,1,19                !通过关键点号循环增加生成线 fc~6/  
   l,69+i,70+i               Z`?Z1SBt  
*enddo rM sd)  
*do,i,1,20                               ~(BvIzzD  
   lesize,124+i,,,1             <0Mc\wy  
*enddo 8b^v@|)N  
save 0nh;0Z  
y}'c)u  
type,1                     !进行拉伸，拉伸后建立的完整模型 UJqDZIvC  
mat,1                       %,l+
?fF  
esys,11                     !OVTs3}  
extopt,aclear,1        !设置拉伸选项，即拉伸完成后删除母单元Mesh200 ?\ qfuA9.  
csys,12 )<
.BN p  
*do,i,1,20 'q#$^='o  
    asel,s,loc,z,-3*(i-1)       !选择拉伸的面 ~s :Ml  
    vdrag,all,,,,,,124+i        !进行拉伸 1nt VM+  
*enddo D
Q<{FN  
allsel,all                   C;u8qVI  
nummrg,all               C Rd1zDB  
numcmp,all                                     ,r&:C48dI  
save   BRTM]tRZ  
Y55Yo5<j/+  
csys,0                        !赋予隧道所在地层材料参数 m_ |:tU(t  
nsel,s,loc,y,10.5,12.8         |\1!*Qp  
esln,s                         (#dwIBBFt  
mpchg,2,all                     jK[~dY
  
nsel,s,loc,y,6.3,10.5               aetK<9L$  
esln,s                           (&=<UGY(w  
mpchg,3,all                         8f^QO:  
nsel,s,loc,y,-4.5,6.3           _;;'/rs j  
esln,s                         (dL;A0L  
mpchg,4,all                     E!~Ok  
nsel,s,loc,y,-4.5,-9.1               A
c0^`  
esln,s                           t3~ZGOn  
mpchg,5,all   9rB,7%@EL  
nsel,s,loc,y,-9.1,-20.8           bD&^-& G  
esln,s                           DP(JsZ}  
mpchg,6,all   Qj?qWVapA  
!allsel,all !L+4YA  
save Z#+{ksU  
Z/|oCwR  
Finish                  !施加边界条件 lHV&8fny  
/solu                       Yt 9{:+[RK  
csys,0                       QWo_Zg0"  
/view,1,-1,0.5,1               @+gr>a1K#  
/replot                       <$
f7&6B  
eplot                       $yZ(c#L  
save                         1YGj^7V)|Z  
nsubst,10                     3Ishe
"  
nlgeom,on                   /``4!jU  
nropt,full                     +}XFkH~  
pred,on                 !启用非线性求解预测器 (-RZ|VdYg  
lnsrch,on                   -|?I'~[#(  
outres,all,last                 iK#5nY].  
R3HfE*;Z  
nsel,s,loc,x,-47.5       !选择节点       muAI$IRR  
d,all,ux                 !施加X方向约束     rNdeD~\  
nsel,s,loc,x,47.5               k63]Qf=5?N  
d,all,ux                       BD)5br].  
nsel,s,loc,y,-20.8       !选择节点 +w(sDH~kd  
d,all,uy                 !施加Y方向约束 rQ^X3J*`
 
nsel,s,loc,z,0           !选择节点 4QK~qAi  
nsel,a,loc,z,-60 y?ps+ce93  
d,all,uz                 !施加Z方向约束 986y\9Zu  
allsel,all E@/yg(?d=  
acel,0,10,0              !施加重力 {HJzhIgCf  
step=0                     !用于记录计算步的参数 =~OH.=9\  
save   (1
L9K;  
@_gCGI>Q  
time,1                 4`x.d  
solve                       >O{U4_j@(  
save                         K
RY%B[k  
finish   A]z~Dw3  
                           h83;
}>  
/POST1                    !进入后处理器 {Hv/|.),hu  
PLNSOL,S,Y,0,1           !SY应力图 DNP%]{J  
PLNSOL,U,Y,0,1   NW/RQ(  
      

haiahsiahas as a sa