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

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

上传命令流希望大家交流 v.6K;TY.  
FINISH gSk0#Jt  
/CLEAR NX&Z=ObHu}  
!(1)定义常量 M\A6;dz'  
  L    = 8 nu `R(2/  
  H    = 6 {O!B8a   
  B    = 8 GUps\:ss  
  E1    = 3.0E6 =/+f3  
  NU1  = 0.22 v05$"Ig  
  DEN1=2.15E3 j" 5 +"j  
  E2   = 1.26E7 ?&JKq^9\I  
  NU2  = 0.37 EX/{W$ &K  
  DEN2=1.8E3 YS bS.tq  
  E3    =3.5E9 ?s@=DDB\u  
  NU3  = 0.25 uC|bC#;  
  DEN3=2.5E3 `hDH7u!U.  
  SIZE = 1 QKIg5I-  
6(/*E=bOKV  
!(2)工作环境设置 gs(ZJO1 /L  
/FILENAME,2010.4.26-1
f<;eNN  
/TITLE,PIT OF WORKSHOP &9EcgazV
 
qm6X5T  
!(3)进入前处理 tfdb9#&?  
/PREP7 #'5{ ?Cb  
)C$pjjo/`  
!(4)定义单元类型 I&f!>y?,Z  
ET, 1, PLANE42 D$Ao-6QE W  
ET, 2, SOLID45 '3=@UBs  
ET,3,beam4 0*,]`A=  
R,1,0.64,0.03413,0.03413,0.8,0.8  !桩的实常数 m>f8RBp]'  
R,2,0.48,0.0256,0.0144,0.6,0.8       ! 连梁的实常数 ;)Sf|  
@!*I mN
MI  
  MP, EX, 1, E1 6J <.i  
  MP, NUXY, 1, NU1 A"6

&  
  MP,DENS,1,DEN1 `(xzCRX  
  TB,dp,1 qS?^(Vt|R  
  tbdata,1,10E3,4 qb$M.-\ne  
\s6VOR/  
  MP,EX,2,E2 :)F0~Q  
  MP,NUXY,2,NU2 oxug  
  MP,DENS,2,DEN2 _~ 3r*j  
  TB,dp,2 eMtQa;Lc9o  
  tbdata,1,21E3,12 Yo\%53w/  
HfEl TC:3f  
  MP,EX,3,E3 $ncP#6  
  MP,NUXY,3,NU3 _QneaPm%  
  MP,DENS,3,DEN3 a28`)17z  
  TB,dp,3 NbK67p:  
  tbdata,1,670E3,35 8{)N%r  
   p: u@? k  
  mp,ex,4,2.8E9      !桩、梁的弹性模量 :464~tHI[`  
  mp,nuxy,4,0.3        !桩的泊松比 y m?uj4I{  
  mp,dens,4,2500       !桩的密度 '26 ,.1  
!O4)YM  
.ZupsS9l  
!(5)建立断面模     S&F;~  
!Y轴坑壁 'VV"$`Fu"  
*do,i,1,5,1           !x=40,z=30,y变 >VJ"e`  
k,,4,i-1,6 ^*F'[!. p  
*end do nnTiu,2R  
v7<r-<I[  
*do,i,1,5,1           !x=40,z=11,y变 QzY5S0  
k,,4,i-1,2 'jYKfq~_cJ  
*end do eC_i]q&o|  
HE-ErEtGB  
2+yti,s+/  
hLaQ[9  
!.7m4mKzo  
*do,i,1,5,1           !X=40,Y变的桩 e"ClG/M_XS  
l,i,i+5 A27!I+M  
*end do =WK's8FB;8  
type,3               cXNR<`   
mat,4   EJ(z]M`f  
real,1 |AS9^w  
lsel,s,line,,1,5,1 =<Ss&p>  
lesize,all,,,4 ZeH=]G4Zv7  
lmesh,all B3p79j  
q*}$1 zb  
< Q6  
4 ;^  
*do,i,1,4,1               LAuaowE\v  
l,i,i+1                 !连梁 @eP(j@(^  
*end do ,#hNHFa'JH  
type,3               1~S''[  
mat,4   #6
6u<FaG  
real,2 [6%y RQ_  
lsel,s,line,,6,9,1 ySHio;g9  
lesize,all,,,1 +<3e@s&  
lmesh,all ootkf=  
EE*FvI`  
K*5Ij]j&  
yJ!,>OQ%'  
.]/k#Hv  
ZFtN~Tg  
RggO|s+0;  
|9mGX9q  
!x轴坑壁 '&Tz8.jp~  
u+ 8wBb5!  
*do,i,2,5,1           !y=40,z=30,x变 |%j7Es  
k,,4+i-1,4,6 UZo[]$"Q`  
*end do JQSczE3  
nr%P11U\c  
*do,i,2,5,1           !y=40,z=12,x变 IW=%2n(<1  
k,,4+i-1,4,2 O`.IE? h#  
*end do ZSn6JV'g  
VW:Voc  
J74kK#uF=  
*do,i,1,4,1           !y=40,x变的桩 Rc}#4pM8  
l,i+10,i+10+4 ,9W!cD+0  
*end do oSH]TL2@Cd  
type,3               WB"90!  
mat,4   9)'wgI#  
real,1 qR^+K@*|  
lsel,s,line,,10,13,1 2K2_
-  
lesize,all,,,4 {%z
5^o1)  
lmesh,all jNAboSf2Y  
@XH@i+{B  
!lI1jb"
 
C{YTHNn  
l,5,11 MYUL y2)  
type,3                     !连梁 cB2~W%H  
mat,4 7+m.:~H3}  
real,2 <4*)J9V^s=  
lsel,s,line,,14 Fb^,%K:  
lesize,14,,,1 PBAQ KQ  
lmesh,14 ^Fwdi#g  
?IR+OCAA  
*do,i,1,3,1 &usum~@  
l,i+10,i+10+1 '^BTa6W}m  
*end do yXro6u?
rC  
type,3 OmaG|2u  
mat,4 "O%xQ N  
real,2 8-)@q|  
lsel,s,line,,15,17,1 KvlLcE~
`o  
lesize,all,,,1 D#11 N^-K  
lmesh,all /z4n?&tM  
cN| gaL  
!创建矩形区域 $WbfRyXi7'  
ExSy/^4f  
RECTNG,0,8,0,8 NJ(H$tB@  
LSEL,S,LINE,,18,21,1 Edl .R}&1  
LESIZE,ALL,SIZE oTf^-29d  
AMESH,1 ,>3|\4/Q  
vVH*\&H\T  
SAVE z?(QM:  
]dnB,  
8&v%>wxR@  
!进行拉伸 \`YV)"y" ~  
p;LF-R  
TYPE,2 Wp
//SV  
MAT,1 kDWvjT  
EXTOPT,ESIZE,6 6"La`}B(T8  
EXTOPT,aclear,1 8z*/J=n  
VOFFST,1,6 j.UO>1{7  
Z|zT%8.8N  
ALLSEL,ALL ]]|vQA^  
nummrg,all Bk~%  
numcmp,all ]ru
UX  
SAVE {_X&{dZLX  
"@@Z{  
NSEL,S,LOC,Z,2,4 `MXGEJF  
ESLN,S,1             !选择和所选节点有关的单元 C8 "FTH'  
MPCHG,2,ALL q#.+P1"U  
ESEL,ALL CR`}{?2H
 
allsel,all rTM0[2N  
%s[ n2w  
NSEL,S,LOC,Z,0,2 [2nP
r^  
ESLN,S,1            !选择和所选节点有关的单元 8\CmM\R  
MPCHG,3,ALL %o#|zaK  
ALLSEL,ALL _~uYNvmg  
iFaC[(1@a  
esel,s,type,,3 aTPmW]w6  
nsle,s,all HH'5kE0;d  
MPCHG,4,ALL _u8d`7$*%  
ALLSEL,ALL w-?Cg8bq<  
GsC4ty  
/PNUM,MAT,1 f}+8m .g2  
[^A>hs*  
FINISH r#/Bz5Jb*  
!加载与求解 |%5nV=&\  
zNJ-JIo%  
/SOLU -LEpT$v|  
NLGEOM,ON          !打开大变形 x|<89o L  
NROPT,FULL         !设定牛顿-拉弗森选项 )`
90*  
\UBQ:+3  
TIME,1              !设定荷载步结束时间值 O_y?53X  
AUTOTS,ON             !打开自动荷载步 T'e p&tNY  
NSUBST,10           !定义子荷载步为10               g(F?qP_K  
PRED,ON           !打开非线性求解预测器 pN7 v7rs  
OUTRES,ALL,ALL    !保存每一个子布的最后一个迭代计算步结果 %Dm:|><V$b  
:*s+X$x,<  
!施加约束 "\1QJ  
NSEL,R,LOC,Y,0 *p7_rY  
D,ALL,UY J1{ucFa  
ALLSEL,ALL g'1ASMuR  
'.d]n(/lZd  
NSEL,R,LOC,Y,B @2)ImgK[  
D,ALL,UY :]^FTnO  
ALLSEL,ALL dX5|A_Ex  
gPg2Ve0Qy  
NSEL,R,LOC,X,0 odC}RdN  
D,ALL,UX <
Z8^.t)|  
ALLSEL,ALL  0@dN$e  
gVI2{\a  
NSEL,R,LOC,X,L ,]|*~dd>G  
D,ALL,UX ~TfQuIvQB  
ALLSEL,ALL \h 1T/_4  
"4e{Cq  
NSEL,S,LOC,Z,0 wQ/Z:  
D,ALL,ALL (w{C*iB  
ALLSEL,ALL e&J_uG  
R SqO$~  
ACEL,0,0,10 !施加重力 %t<ba[9F  
C#4/~+  
SAVE &[|P/gj#>  
N(BCe\FV  
!求解自重应力场 hz~jyH.h_  
TIME,1 2rJe
ON  
esel,s,type,,3 >gLLr1L\  
ekill,all uw`fC%-xh  
esel,all 5D%gD
w+"  
esel,s,live         !选择活的单元 CxRhMhvP  
nsle,s             !选择活单元上的节点 wx%TQ!  
nsel,invert         !反向选择，即选择了死单元上的节点 dw]wQ\4B  
d,all,all,0         !将死单元上的节点约束所有位移，使其不参与矩阵运算 "[BDa}Il  
nsel,all |8YP8o  
esel,all t?:Q  
eplot q^k]e{PD  
solve 7f\@3r  
K{x FhdW  
SAVE K?:wX(JYT  
FINISH Y94/tjt  
WbP*kV{  
/POST1 eq#x~O4  
PLNSOL,S,Z,0,1     !SZ应力云图 T^~9'KDd  
PLNSOL,U,Z,0,1      !UZ位移云图 Om=*b#k  
FINISH ,Y3wXmG  
!分层开挖，开挖第1步 }^R_8{>k  
/SOLU Z2x%  
ANTYPE,,REST   !分析类型为静态分析,重启动分析，启动点为上次求解结束时 b[_${in:  
TIME,2 x3++JG  
AUTOTS,ON             !打开自动荷载步 \<R.F  
DELTIM,0.1,0.001,0.2,ON `DPR >dd@  
r06M.r  
NSEL,S,LOC,X,4,8 (66DKG   
NSEL,R,LOC,Y,0,4 Udh!%QP%[w  
NSEL,R,LOC,Z,3,6 8K/o/  
ESLN,S,1 11Hf)]M   
EKILL,ALL Lv{xwHnE  
t|Ipxk.)  
ALLSEL,ALL |rbl sL2?Z  
Y4B<]C4  
esel,s,type,,3 6b8@6;&LI  
ealive,all 6(X5n5C  
allsel,all KZxA\,Y'5  
SOLVE 7JNy;$]/  
SAVE k v}<u  
FINISH
n_Ka+Y<  
]$vJK  
!地层位移分析 %$l^C!qcY  
/POST1 M:~/e8Xv  
LCDEF,1,1,, (cj3[qq  
LCDEF,2,2,, #!
FLX*,  
LCASE,2 u8wZ2j4S  
LCOPER,SUB,1,,, g#ZuRL  
OZ$"P<X_"  
ESEL,S,LIVE h!*++Y?&0  
PLNSOL,U,Z,0,1 BI\+NGrB  
Q l
$t  
($oO, c'z  

yO6 _Gq{  

我等高手出来解决问题~

没点反应啊，高手门都出来露露脸吧 J1<fE(X  

围观，学习~一起等高手

期待高手出现

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

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

不知道，问什么上传不了，图片。 }m H>lN  
.HBvs=i  
对于，命令流没有多大问题。 *mQit/k.  
只是，提出个人看法， >&&xJ5  
yx
2.7h3  
1  在初始地应力的时候，可不选取dp模型，弹性的即可。然后，初始地应力后，加dp模型。 <n\.S  
4cQ5E9  
2   在前面模型材料属性 改变是在单元上。对于，改回材料3的的话，应 &FrB6y  
X4$e2f  
esel,s,type,,3 iW?9oe  
nsle,s,all x3s^u~C)(w  
MPCHG,4,ALL !T*B{+|  
gxOmbQt@;  
改为 OD]J@m  
esel,s,type,,3 B0i}Y-Z  
MPCHG,4,ALL 5 DB>zou   
其实，楼主，改变节点的话，没有多大问题。因为单元的计算式通过节点来算的。

qingcaozheng:上传命令流希望大家交流 &@utAuI  
FINISH l
iw 9:@+V  
/CLEAR g~hk-nXL.  
!(1)定义常量 +.#S[G  
  L    = 8 i?_D]BY4  
....... (2010-04-26 17:27)  !BQ ELB$0  

5'c+313 lm  
不知，楼主，现在还在吗？若在，希望，看一下元模型，进行一下计算。