[求助]ANSYS模拟基坑三维开挖 [复制链接]

SOLID45模拟土，beam4模拟支护桩，用DP模型，但是当我加上支护和不加支护对比坑壁的水平位移，发现加上支护后更大，改变桩径也是桩径越大侧移越大，很困惑，那位给指点下啊，不胜感激！

上传命令流希望大家交流 iZDb.9@&t  
FINISH S20nk.x  
/CLEAR >uP{9kDm  
!(1)定义常量 F1{?]>G  
  L    = 8 ~:ub  
  H    = 6 2Q`PUXj  
  B    = 8 :JTRRv  
  E1    = 3.0E6 mTrI""Jsu;  
  NU1  = 0.22 A$H+4L  
  DEN1=2.15E3 ^^t]vojX  
  E2   = 1.26E7 /Gh x2B  
  NU2  = 0.37 ;:8jxkx6%  
  DEN2=1.8E3 V
)WIfRs  
  E3    =3.5E9 eY6gb!5u  
  NU3  = 0.25 EDPI*@>  
  DEN3=2.5E3 S+Yy  
  SIZE = 1 }QApeZd+q  
}5fI*v  
!(2)工作环境设置 n##w[7B*  
/FILENAME,2010.4.26-1 >h/)r6  
/TITLE,PIT OF WORKSHOP ?-(w][MT\  
~8Dd<4?F]  
!(3)进入前处理 ';fU.uy  
/PREP7 
H$=h-  
1A)wbH)  
!(4)定义单元类型 i"r.>X'Z  
ET, 1, PLANE42 U:etcnb4w>  
ET, 2, SOLID45 fmZz
BZ_  
ET,3,beam4 ino7!T`  
R,1,0.64,0.03413,0.03413,0.8,0.8  !桩的实常数 $z`cMQ r  
R,2,0.48,0.0256,0.0144,0.6,0.8       ! 连梁的实常数 )+Y\NO?O  
bSeL"   
  MP, EX, 1, E1 R"8})a gw  
  MP, NUXY, 1, NU1 V;[__w  
  MP,DENS,1,DEN1 ~PH1|h6  
  TB,dp,1 rH}Dt@  
  tbdata,1,10E3,4 m\}\RnZu  
g O\f:Pg  
  MP,EX,2,E2 Hjo:;s  
  MP,NUXY,2,NU2 ]VHdE_7
)  
  MP,DENS,2,DEN2 wFoR,oXtL/  
  TB,dp,2 ^}Dv$\;6  
  tbdata,1,21E3,12 Js^r]=\F'  
Q%AS;(d  
  MP,EX,3,E3 iC5JU&l  
  MP,NUXY,3,NU3 .eDxIWW+ft  
  MP,DENS,3,DEN3 Xx9~  
  TB,dp,3 7Da^Jv k  
  tbdata,1,670E3,35 nfd?@34"A2  
   ;|2;kvf"w  
  mp,ex,4,2.8E9      !桩、梁的弹性模量 U(Hq4D  
  mp,nuxy,4,0.3        !桩的泊松比 Li<266#A!  
  mp,dens,4,2500       !桩的密度 }ii]cY  
4%c7#AX[T  
NZw[.s>n  
!(5)建立断面模     B 3,ig9  
!Y轴坑壁 Is*0?9qU  
*do,i,1,5,1           !x=40,z=30,y变
](
U%1  
k,,4,i-1,6 46.q anh  
*end do x/jN&;"/  
AIRVvW~($  
*do,i,1,5,1           !x=40,z=11,y变 !z4Hj{A_  
k,,4,i-1,2 +~pc%3*  
*end do #Ko+_Hm?4  
7:R{
~|R  
m?kIa!GM=  
MRl*rK  
5 a*'N~  
*do,i,1,5,1           !X=40,Y变的桩 *{D:1S
  
l,i,i+5 S#%JSQo
:  
*end do ,{mf+ 3&$,  
type,3               WSpg(\Cs  
mat,4   H$qdU!c  
real,1
XR=c 8f  
lsel,s,line,,1,5,1 K&"Pm9  
lesize,all,,,4 {]/Jk07  
lmesh,all o`y*yucHI  
oRJP5Y5na  
37M[9m|D*  
>aanLLO  
*do,i,1,4,1               "=Fn.r4I  
l,i,i+1                 !连梁 WHD/s  
*end do 2YP"nj#  
type,3               xx`8>2T#e  
mat,4   3K'o&>}L  
real,2 GWkJ/EX  
lsel,s,line,,6,9,1 hz~CW-47  
lesize,all,,,1 o._#=7|(  
lmesh,all nD?M;XN  
;HLMU36q  
2zrWR%B  
XKPt[$ab  
kvO`]>#;$?  
K&gc5L  
'mp{O  
C_khd"  
!x轴坑壁 oe_[h]Hgl  
{i7Fu+xZj  
*do,i,2,5,1           !y=40,z=30,x变 2"{]A;@  
k,,4+i-1,4,6 Zn*CJNB  
*end do $nd-[xV  
Iz$W3#hi  
*do,i,2,5,1           !y=40,z=12,x变 O
'sr[  
k,,4+i-1,4,2 m*Lv,yw %a  
*end do r@
@eC['  
gJ[q {b  
|EEz>ci  
*do,i,1,4,1           !y=40,x变的桩 bEXHB  
l,i+10,i+10+4 aBtfZDCfzp  
*end do jiB>.te  
type,3               R;-FZ@u/  
mat,4   ,;ruH^  
real,1 uRq#pYn
@  
lsel,s,line,,10,13,1 $E@.G1T [  
lesize,all,,,4 $^+K
R]\q  
lmesh,all ThqfZl=V  
4;~lpty  
$q
@RHcj  
Qc =lf$  
l,5,11 &u`rE""  
type,3                     !连梁 /dvnQW4}8  
mat,4 M0IqQM57N  
real,2 `R=_t]ie  
lsel,s,line,,14 ", Rw%_  
lesize,14,,,1 ws(}K+y_  
lmesh,14 !m O] zn  
u.K'"-xt4K  
*do,i,1,3,1 ZtK%b+MBP  
l,i+10,i+10+1 g}og@UY7#  
*end do ; dHOH\,:  
type,3 eq 1 4  
mat,4 ,["|wqM  
real,2 rxK[CDM,  
lsel,s,line,,15,17,1 ZzaW@6LJF  
lesize,all,,,1 &/#Tk>:  
lmesh,all f+88R=-u6S  
{c EKz\RX  
!创建矩形区域 YHv,Z|.w  

;NVTn<Uj  
RECTNG,0,8,0,8 \w\47/k{  
LSEL,S,LINE,,18,21,1 $8ww]}K  
LESIZE,ALL,SIZE H=SMDj)s+  
AMESH,1 Q\>SF  
Wes"t}[25  
SAVE cM'[;u  
g8qAJ4  

d~b
H!P  
!进行拉伸 X|lmH{k
f  
]~x/8%e76  
TYPE,2 AeQ&V d|  
MAT,1 8_,wOkk_B  
EXTOPT,ESIZE,6 G;#t6bk  
EXTOPT,aclear,1 \]:NOmI^'  
VOFFST,1,6 WMRgf~TY=2  
Fu$Gl$qV?%  
ALLSEL,ALL q>lkLHS  
nummrg,all nsw8[pk  
numcmp,all ?j$8Uy$$  
SAVE LFM5W&?  
5N</Z6f'o  
NSEL,S,LOC,Z,2,4 Kz2^f@5=F  
ESLN,S,1             !选择和所选节点有关的单元 ScmzbDu  
MPCHG,2,ALL  ~9jP++&  
ESEL,ALL 52R.L9Ai  
allsel,all +#"Ic:  
FbNQ  
NSEL,S,LOC,Z,0,2 P##(V!YR  
ESLN,S,1            !选择和所选节点有关的单元 GN"LU
>9|  
MPCHG,3,ALL hnLgsz  
ALLSEL,ALL [67f;?b  
{_N(S]Z  
esel,s,type,,3 *E q7r>[  
nsle,s,all 7Z:l;%]K  
MPCHG,4,ALL vh3iu+  
ALLSEL,ALL $,v+i -  
zSufU2  
/PNUM,MAT,1 .xQ'^
P_q  
~5qZs"ks  
FINISH T1W9@9,s  
!加载与求解 M#<fh:>  
E6\~/=X=%  
/SOLU z{+;'9C  
NLGEOM,ON          !打开大变形 [#fqyg  
NROPT,FULL         !设定牛顿-拉弗森选项 cx%9UK*c  
b6M)qt9R  
TIME,1              !设定荷载步结束时间值 *1]k&#s  
AUTOTS,ON             !打开自动荷载步 il
l'KPy  
NSUBST,10           !定义子荷载步为10               qBrZg  
PRED,ON           !打开非线性求解预测器 3T^dgWXEG  
OUTRES,ALL,ALL    !保存每一个子布的最后一个迭代计算步结果 /faP]J)  
u\Q**m2XP  
!施加约束 fe`G^hV  
NSEL,R,LOC,Y,0 pY5HW2TsY|  
D,ALL,UY +GtGyp  
ALLSEL,ALL Jy NY *  
(j Q6~1  
NSEL,R,LOC,Y,B ?PtRb:RHt  
D,ALL,UY e~[z]GLO%  
ALLSEL,ALL _@?Jx/`;bk  
XQ y|t"Vq>  
NSEL,R,LOC,X,0 d'nuk#r  
D,ALL,UX T:5%sN;#O  
ALLSEL,ALL dllf~:b  
`%:(IGxz  
NSEL,R,LOC,X,L 0s[3:bZ\Ia  
D,ALL,UX X51$5%  
ALLSEL,ALL hf5SpwxLiH  
tce8*:rNH  
NSEL,S,LOC,Z,0 W
C;a  
D,ALL,ALL $=a$z"  
ALLSEL,ALL $ aBSr1  
e[o ;l  
ACEL,0,0,10 !施加重力 $6%;mep  
W%QtJB1)  
SAVE &!fcLJd  
lD=j/  
!求解自重应力场 A{"t0Ai='0  
TIME,1 Gf.o{  
esel,s,type,,3 rZ4<*Zegv  
ekill,all l+qtA~V&2  
esel,all {/!"}{G1e  
esel,s,live         !选择活的单元 &Y2P!\\2  
nsle,s             !选择活单元上的节点 \`5u@Nzx  
nsel,invert         !反向选择，即选择了死单元上的节点 eV}Tx;1|}  
d,all,all,0         !将死单元上的节点约束所有位移，使其不参与矩阵运算 *%+buHe  
nsel,all >5W"a?(  
esel,all L uKm  
eplot Pt;Ahmi  
solve Y|jesa {x  
ofv
1G=P  
SAVE q9]L!V9Rv  
FINISH )f&]H}  
.[s82c]]6  
/POST1 QP0X8%+p  
PLNSOL,S,Z,0,1     !SZ应力云图 T<
GD!j(  
PLNSOL,U,Z,0,1      !UZ位移云图 I"?&X4%e  
FINISH Qj[O$L0 $  
!分层开挖，开挖第1步 l[{}ZKZ  
/SOLU =x]dP
.  
ANTYPE,,REST   !分析类型为静态分析,重启动分析，启动点为上次求解结束时 84cH|j`w  
TIME,2 &h[}5  
AUTOTS,ON             !打开自动荷载步 rmJ847%y`  
DELTIM,0.1,0.001,0.2,ON 1<e%)? G  
cYS+XBz  
NSEL,S,LOC,X,4,8 \,t<{p
_Q  
NSEL,R,LOC,Y,0,4 <cjTn:w  
NSEL,R,LOC,Z,3,6 kfECC&"  
ESLN,S,1 Lxrn#Z eM  
EKILL,ALL ^C T}i'  
Xh!Pg)|E  
ALLSEL,ALL M&-/&>n!  
Lw
k-  
esel,s,type,,3 IV*$U7~  
ealive,all @b({QM|  
allsel,all jo
#F&  
SOLVE 1yc@q8  
SAVE vH8%a8V  
FINISH .q}k  
cNvcpv  
!地层位移分析 yt,;^o^  
/POST1 fdHxrH>*  
LCDEF,1,1,, c17==S  
LCDEF,2,2,, S`.-D+.68  
LCASE,2 f/m6q8!L{  
LCOPER,SUB,1,,, RjHpC7b*%  
sR
nMBW.  
ESEL,S,LIVE L@GD$F=<0  
PLNSOL,U,Z,0,1 pSUp"wch  
=1Jo-!{{  
T6#GlO)8)  
>Sah\u`  

我等高手出来解决问题~

没点反应啊，高手门都出来露露脸吧 ="K>yUfcFl  

围观，学习~一起等高手

期待高手出现

在学深基坑开挖围护模拟，来围观一下

楼主，我已按你的code来做了一下。 c+=&5=i[3  
你这个是缩小的模型。初步计算结果，没有你说的情况。加支护后，基坑的竖向和水平位移比不加支护时较小。整个云图，比较正常。如下照片。

不知道，问什么上传不了，图片。 /V&Y@j  
AM>:AtY  
对于，命令流没有多大问题。 (agdgy:#  
只是，提出个人看法， Cp%|Q.?  
1*f
A>v  
1  在初始地应力的时候，可不选取dp模型，弹性的即可。然后，初始地应力后，加dp模型。 @_{"ho  
!_@%/I6  
2   在前面模型材料属性 改变是在单元上。对于，改回材料3的的话，应 fvDwg  
I1gu<a  
esel,s,type,,3 (P?
9Jct  
nsle,s,all C7#$s<>TO  
MPCHG,4,ALL z( wXs&z;  
KVuv%?  
改为 jJ|;Nwm<[  
esel,s,type,,3 =>G A_  
MPCHG,4,ALL 7F.,Xvw&@  
其实，楼主，改变节点的话，没有多大问题。因为单元的计算式通过节点来算的。

qingcaozheng:上传命令流希望大家交流 yC[}gHv  
FINISH &~i1 @\]  
/CLEAR Nx99dr  
!(1)定义常量 >9 iv>  
  L    = 8 %^S1 fUwT  
....... (2010-04-26 17:27)  4Xz|HU?  

<*[(t;i  
不知，楼主，现在还在吗？若在，希望，看一下元模型，进行一下计算。