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

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

上传命令流希望大家交流 (3fPt;U  
FINISH L?f qcW{  
/CLEAR 1URsHV!xcM  
!(1)定义常量 M[,^KJ!  
  L    = 8 6Bdyf(t  
  H    = 6 FOp_[rR   
  B    = 8 g{a d0.y,  
  E1    = 3.0E6 hEcYpng~  
  NU1  = 0.22 s1=u{ET  
  DEN1=2.15E3 nHU3%%%cU  
  E2   = 1.26E7  y h-9u  
  NU2  = 0.37 }#YQg0(  
  DEN2=1.8E3 Q-Y@)Mf~?0  
  E3    =3.5E9 [+DNM 2A  
  NU3  = 0.25 rk|a'&  
  DEN3=2.5E3 Fe4esg-B<  
  SIZE = 1 ;U7t  
M9g1d7%  
!(2)工作环境设置 AIfk"2  
/FILENAME,2010.4.26-1 S ljZ~x,!  
/TITLE,PIT OF WORKSHOP mh8nlB  
h.LSMU (O  
!(3)进入前处理 t^$Div_%G  
/PREP7 Ph\F'xROe  
?M<|r11}  
!(4)定义单元类型 `w=!o.1  
ET, 1, PLANE42 riEqW}{  
ET, 2, SOLID45 f[M"EMy  
ET,3,beam4 Kz^hQd  
R,1,0.64,0.03413,0.03413,0.8,0.8  !桩的实常数 :wCC^Y]  
R,2,0.48,0.0256,0.0144,0.6,0.8       ! 连梁的实常数 $y4M#yv  
:+A;TV  
  MP, EX, 1, E1 ]Hj`2\KD.d  
  MP, NUXY, 1, NU1 dh,7iQ s  
  MP,DENS,1,DEN1 ~$ WQ"~z  
  TB,dp,1 9oD#t~+F4  
  tbdata,1,10E3,4 tQnJS2V"{u  
g[<uwknf  
  MP,EX,2,E2 B#4S/d{/  
  MP,NUXY,2,NU2 5oa]dco  
  MP,DENS,2,DEN2 }'_:XKLj  
  TB,dp,2 ndt8=6p  
  tbdata,1,21E3,12 m[#%/  
=z%s8D2  
  MP,EX,3,E3 @f'AWeJ2  
  MP,NUXY,3,NU3 c$.T<r)Z  
  MP,DENS,3,DEN3 nTQ (JDf  
  TB,dp,3 2c*2\93>  
  tbdata,1,670E3,35 C9+Dw#-fV  
   rN'k4V"K  
  mp,ex,4,2.8E9      !桩、梁的弹性模量 F1
3%)G(  
  mp,nuxy,4,0.3        !桩的泊松比 <v-92?  
  mp,dens,4,2500       !桩的密度 "lb\c  
Xa6qvg7/  
CCwK8`%  
!(5)建立断面模     g&8.A(  
!Y轴坑壁 ^)'||Ly  
*do,i,1,5,1           !x=40,z=30,y变 .Px,=56$X  
k,,4,i-1,6 ^f"&}%"M  
*end do @@xF#3   
;WPI+`-  
*do,i,1,5,1           !x=40,z=11,y变 E<P*QZ-C3  
k,,4,i-1,2 s-o~@(r6  
*end do *xho  
|O^V)bZmx  
\!6t  
N}1-2  
kfIbgya  
*do,i,1,5,1           !X=40,Y变的桩 &A#90xzF  
l,i,i+5 _4A&%>   
*end do 3pzOt&T|w  
type,3              
uof0Oc.  
mat,4   yl|R:/2V  
real,1 PK9Qm'W b  
lsel,s,line,,1,5,1 Pyit87h{  
lesize,all,,,4 2}1!WIin  
lmesh,all 13]y)(  
m.
/*LXU  
!FO:^P  
KN|'|2/|  
*do,i,1,4,1               Zj5NWzj X  
l,i,i+1                 !连梁 ug>]U ~0  
*end do HKT{IP+7(L  
type,3               K T"h74@  
mat,4   <4SF~i  
real,2 ;2 \<M6  
lsel,s,line,,6,9,1 ~9]Vy (L  
lesize,all,,,1 + >Fv*lux  
lmesh,all VdYOm  
+# A|Zp<  
8Na}Wp;|Gi  
HXz iDnj  
}:tAKO=+  
1Z=;Uy\  
Gu<W:n[  
NF@i#:  
!x轴坑壁 y;*My#  
AZ]Z,s6  
*do,i,2,5,1           !y=40,z=30,x变 8VK
b*  
k,,4+i-1,4,6 !G@V<'F  
*end do A89Y;_4y  
E%KC'TN^D  
*do,i,2,5,1           !y=40,z=12,x变 "4k"U1  
k,,4+i-1,4,2 F5Cqv0HV  
*end do %YsRm%q  
.k
z(V5  
..sJtA8  
*do,i,1,4,1           !y=40,x变的桩 9Vh_XBgP  
l,i+10,i+10+4 _q2`m  
*end do 7UY('Q[  
type,3               ^!XU+e+:0  
mat,4   h:eN>yW  
real,1 *s 1D\/H  
lsel,s,line,,10,13,1 Ul7,k\q@  
lesize,all,,,4 YeR7*[l  
lmesh,all "`H=AX0  
)2   
Tm'lN5}&9  
K7YT0cG  
l,5,11 )[F46?$vrk  
type,3                     !连梁 L2do2_  
mat,4 1ZGQhjcx  
real,2 Z%(Df3~gmm  
lsel,s,line,,14 OD>u$tI9  
lesize,14,,,1 KI^q 5D ?  
lmesh,14 UzP@{?  
sf=%l10Fk#  
*do,i,1,3,1 .oW~:mY  
l,i+10,i+10+1 'lSnyW{  
*end do #h}
IUR  
type,3 ~`a#h#  
mat,4 h/fb<jIP1  
real,2 ;5zjd,  
lsel,s,line,,15,17,1 }j]<&I}  
lesize,all,,,1 (Qw`%B  
lmesh,all Y,p2eAss  
hJs&rpN  
!创建矩形区域 W\ZV0T;<]  
Ai
R%MD  
RECTNG,0,8,0,8 S S7D1  
LSEL,S,LINE,,18,21,1 IX > j8z[  
LESIZE,ALL,SIZE w0F:%:
/  
AMESH,1 Rq~ >h99M  
Phk`=:xh  
SAVE fbW,0  
[U_Q 2<H  
4IH0
un
 
!进行拉伸
8
tG/VE[  
W_Ws3L1;N  
TYPE,2 t\E-6u  
MAT,1 y'i:%n}I  
EXTOPT,ESIZE,6 BHiw!S<  
EXTOPT,aclear,1 98<bF{#0WM  
VOFFST,1,6 t/EMBfLc  
Y:#kel<  
ALLSEL,ALL &eLQ;<qO*|  
nummrg,all Hj-<{#,  
numcmp,all y"%iD`{  
SAVE kM}ic(K  
c+YYM :S  
NSEL,S,LOC,Z,2,4 R  oF  
ESLN,S,1             !选择和所选节点有关的单元 , .NG.Q4f  
MPCHG,2,ALL [7ek;d;'t  
ESEL,ALL >8.v.;`  
allsel,all hA&j?{  
Oa3=+_C~$1  
NSEL,S,LOC,Z,0,2 uH~ TugQ~  
ESLN,S,1            !选择和所选节点有关的单元 -X6\[I:+A  
MPCHG,3,ALL A$$R_3ne  
ALLSEL,ALL V8|q"UX  
%s&E-*X  
esel,s,type,,3 vi()1LS/!  
nsle,s,all >V ]*mS%K  
MPCHG,4,ALL 8kn]_6:3i  
ALLSEL,ALL 071E%
u,  
xhp-4  
/PNUM,MAT,1 !Barc,kA  
7o 83|s.Bm  
FINISH W6!4Qyn  
!加载与求解 !Sr0Im0  
b*F :l#  
/SOLU \M1M2(@pDJ  
NLGEOM,ON          !打开大变形 #E~WVTOw  
NROPT,FULL         !设定牛顿-拉弗森选项 GYy!
`E  
bl+@}+A  
TIME,1              !设定荷载步结束时间值 _g/TH-;^  
AUTOTS,ON             !打开自动荷载步 cVv4gQD\  
NSUBST,10           !定义子荷载步为10               R)DNFc:  
PRED,ON           !打开非线性求解预测器 IJb1) ZuR  
OUTRES,ALL,ALL    !保存每一个子布的最后一个迭代计算步结果 g)|++?  
?AM8*w  
!施加约束 DFZ:.6p  
NSEL,R,LOC,Y,0 u3v6$CD?  
D,ALL,UY Q,`2DHhK  
ALLSEL,ALL v1tN DyM6  
3Wtv+L7Br  
NSEL,R,LOC,Y,B `!>dbR&1  
D,ALL,UY ~_^o?NE,  
ALLSEL,ALL e*'|iuDrY  
4jyr\=42F'  
NSEL,R,LOC,X,0 W;)FNP|MT  
D,ALL,UX @{$Cv"6769  
ALLSEL,ALL \Z7([Gh  
<PuB3PEvV  
NSEL,R,LOC,X,L ;Kd{h  
D,ALL,UX `__?7"p )\  
ALLSEL,ALL ,VcDvZ7  
BD-c 0-+m  
NSEL,S,LOC,Z,0 Kr}M>hF+|  
D,ALL,ALL (^s&#_w03  
ALLSEL,ALL ?H86Wbz  
)su
<Ji*  
ACEL,0,0,10 !施加重力 ?}lpo; $  
O%q;,w{prW  
SAVE O|7{%5h  
r{N{!"G  
!求解自重应力场 <:yq~?  
TIME,1 tX`[6`  
esel,s,type,,3 ~m;MM)_V  
ekill,all +68K[s,FD  
esel,all ~)_
?:.Da  
esel,s,live         !选择活的单元 "!_ 4%z-  
nsle,s             !选择活单元上的节点 U(
&nh?  
nsel,invert         !反向选择，即选择了死单元上的节点 h9OL%n 7m'  
d,all,all,0         !将死单元上的节点约束所有位移，使其不参与矩阵运算 0)]C&;}_M  
nsel,all E(4lu%  
esel,all qzbkxQu]g  
eplot ?GD?J(S  
solve gp&& c,  
Hk~ gcG  
SAVE !O-_
Dp\#  
FINISH A(@gv8e[H^  
))+98iU1s  
/POST1 zt
>_)&b  
PLNSOL,S,Z,0,1     !SZ应力云图 *e"GQd?  
PLNSOL,U,Z,0,1      !UZ位移云图 X!A]V:8dk  
FINISH "6o}qeB l  
!分层开挖，开挖第1步 V]PhXVJ  
/SOLU `J7L
ecgo  
ANTYPE,,REST   !分析类型为静态分析,重启动分析，启动点为上次求解结束时 i,b7Ft:F&  
TIME,2 UE$[;Zg  
AUTOTS,ON             !打开自动荷载步 ?e|:6a+[f  
DELTIM,0.1,0.001,0.2,ON ~1]2A[`s!  
x_iy;\s1  
NSEL,S,LOC,X,4,8 cZX&itVc:  
NSEL,R,LOC,Y,0,4 o_KcnVQ\  
NSEL,R,LOC,Z,3,6 -O>mY)  
ESLN,S,1 w5/X{  
EKILL,ALL en#g<on  
IWRq:Gw  
ALLSEL,ALL {s^ryv_}  
+(P43XO08  
esel,s,type,,3 JE:n`
l/p  
ealive,all zam0(^=  
allsel,all 3L%r_N*a  
SOLVE FC-*?  
SAVE F@(}=w^(A  
FINISH YU0HySP:  
\t(r@qq  
!地层位移分析 f]6`GsE  
/POST1 }GNH)-AG)$  
LCDEF,1,1,, n; '~"AG)  
LCDEF,2,2,, 0N[DV]  
LCASE,2 hz|z&vyP  
LCOPER,SUB,1,,, xW;[}t-QS  
}br<2?y,  
ESEL,S,LIVE >@89k^#Vc  
PLNSOL,U,Z,0,1 IEr`6|X  
ysT!^-&p  
PdN\0B`  
.dLX'84fY  

我等高手出来解决问题~

没点反应啊，高手门都出来露露脸吧 !J(,M)p!  

围观，学习~一起等高手

期待高手出现

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

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

不知道，问什么上传不了，图片。 kO"aE~  
<9ph c  
对于，命令流没有多大问题。 Vk T3_f  
只是，提出个人看法， 7s3<}  
d_B5@9e#  
1  在初始地应力的时候，可不选取dp模型，弹性的即可。然后，初始地应力后，加dp模型。 K#m o+n5-;  
nK=V`  
2   在前面模型材料属性 改变是在单元上。对于，改回材料3的的话，应 {u3u%^E;R  
r{&"]'/X  
esel,s,type,,3 /omVMu  
nsle,s,all LK~0ck7  
MPCHG,4,ALL .?:~s8kB  
_zwuK
1e  
改为 [}3Y1t{G  
esel,s,type,,3 1Xc%%j  
MPCHG,4,ALL <Ux;dekz}  
其实，楼主，改变节点的话，没有多大问题。因为单元的计算式通过节点来算的。

qingcaozheng:上传命令流希望大家交流 i|Wn*~yFOO  
FINISH o 8U2vMH  
/CLEAR 'Ud5;?
{  
!(1)定义常量 U>XGJQ<NS  
  L    = 8 $4pW#4/4
 
....... (2010-04-26 17:27)  8Qh/=Ir  

+/
tD$  
不知，楼主，现在还在吗？若在，希望，看一下元模型，进行一下计算。