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

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

上传命令流希望大家交流 JOgmF_(>Z  
FINISH [wIKK/O  
/CLEAR r'PE5xqF  
!(1)定义常量 AG=9b  
  L    = 8 SNxz*`@4  
  H    = 6 69OET_AS>  
  B    = 8 jqqaw  
  E1    = 3.0E6 9<~,n1b>x  
  NU1  = 0.22 1J+3a-0  
  DEN1=2.15E3 X@eg<]'m  
  E2   = 1.26E7 59/Q*7ZJ  
  NU2  = 0.37 ZK!4>OuH`  
  DEN2=1.8E3 !xJFr6G~8  
  E3    =3.5E9 / (.'*biQ  
  NU3  = 0.25 e18T(g_i  
  DEN3=2.5E3 m^Rd Iy)  
  SIZE = 1 $uTlbAuv  
o] S`+ZcV  
!(2)工作环境设置 3liq9P_  
/FILENAME,2010.4.26-1 FBbm4NB  
/TITLE,PIT OF WORKSHOP _we3jzMW  
J;"nm3[.q  
!(3)进入前处理 B*BHF95!  
/PREP7 \|Y{jG<cu  
jUZ[`f;  
!(4)定义单元类型 "1CGO@AXS  
ET, 1, PLANE42 |y'b217t  
ET, 2, SOLID45 R>` ih&,)  
ET,3,beam4 y,1S&k  
R,1,0.64,0.03413,0.03413,0.8,0.8  !桩的实常数 <EOg,"F  
R,2,0.48,0.0256,0.0144,0.6,0.8       ! 连梁的实常数 /o'oF  
D~7%};D[  
  MP, EX, 1, E1 M+\rX1T  
  MP, NUXY, 1, NU1 y#nSk%"t"  
  MP,DENS,1,DEN1 d^_itC;-,  
  TB,dp,1 w0\4Wa  
  tbdata,1,10E3,4 f0g6g!&gf  
@u<0_r t  
  MP,EX,2,E2 w[M5M2CF  
  MP,NUXY,2,NU2 zo87^y5?G  
  MP,DENS,2,DEN2 +:b(%|  
  TB,dp,2 .0KOnLdK  
  tbdata,1,21E3,12 LP8o7%sv!  
6O]Xhe0d@  
  MP,EX,3,E3 %!D_q~"H  
  MP,NUXY,3,NU3 k\%,xf; x  
  MP,DENS,3,DEN3 &F9OZMK=  
  TB,dp,3 &7lk2Q\  
  tbdata,1,670E3,35 B\NcCp`5  
   bL#sn_(m  
  mp,ex,4,2.8E9      !桩、梁的弹性模量 @!,D%]8"  
  mp,nuxy,4,0.3        !桩的泊松比 @eA %(C  
  mp,dens,4,2500       !桩的密度 e ^-3etx  
]~ >@%v&  
ul}
4p{ m[  
!(5)建立断面模     ]&
/0  
!Y轴坑壁 vN'VDvVM  
*do,i,1,5,1           !x=40,z=30,y变 CARq^xI-  
k,,4,i-1,6 K"G(?<>~4c  
*end do i{4'cdr?  
f};!m=b  
*do,i,1,5,1           !x=40,z=11,y变 +-'`Q Ae  
k,,4,i-1,2 Kl!DKeF  
*end do !HbqbS22  
w# xncH:1  

37,L**Dgs  
&~42T}GTWG  
C!`>cUhE{  
*do,i,1,5,1           !X=40,Y变的桩 =CGD ~p`  
l,i,i+5 /;[
}=JL<Q  
*end do (PyTq 5:F  
type,3               }q/(D?  
mat,4   HI7]%<L  
real,1 pEJ#ad  
lsel,s,line,,1,5,1 6@i|Kw(:  
lesize,all,,,4 TIKEg10I  
lmesh,all SG1&a:c+.  
6{lG1\o  
es{cn=\s  
'=-s1c@^  
*do,i,1,4,1               *C tsFS~  
l,i,i+1                 !连梁  %L{  
*end do JIB?dIN 1  
type,3               ]kzv8#  
mat,4   ]9dx3<2_I  
real,2 _Eszr(zJ  
lsel,s,line,,6,9,1 t4C<#nfo  
lesize,all,,,1 j#4+-  
lmesh,all VoWA
tNU  
Cbazwq  
.59KE]u  
6j?FRs  
K%kXS  
4;",@}  
aViJ   
/ O|Td'Z  
!x轴坑壁 J,.j_ii`!  
k q/t]%(  
*do,i,2,5,1           !y=40,z=30,x变 WFQ*s4 R(  
k,,4+i-1,4,6 70d] d+M|  
*end do q.U*X5  
AfuXu@UZ_/  
*do,i,2,5,1           !y=40,z=12,x变 !4i,%Z&6  
k,,4+i-1,4,2 nmTm(?yE  
*end do 7=;
D0SS  
Q|6Ls$'$  
t@l(xnsV  
*do,i,1,4,1           !y=40,x变的桩 =I %g;YK  
l,i+10,i+10+4 .Gjr`6R  
*end do qZ `nZi  
type,3               dw'<"+zO  
mat,4   YLD-SS[/>  
real,1 |C&%S"*+D  
lsel,s,line,,10,13,1 4o)(d=q  
lesize,all,,,4 [BJ$|[11  
lmesh,all C+ZQB)gn  
rDK;6H:u{  
s0_-1VU  
$:T<IU[E  
l,5,11 ab8oMi`z  
type,3                     !连梁 +"TI_tK,S  
mat,4 m*Q[lr=  
real,2 M9g~lKs'  
lsel,s,line,,14 ce 7Yr*ZB  
lesize,14,,,1 F]/L!   
lmesh,14 n.=e)*  
1k
bT@  
*do,i,1,3,1 s@.`"TF.7  
l,i+10,i+10+1 f%`*ba"v  
*end do UZ[/aq  
type,3 \Ac}R'  
mat,4 !5yRWMO9X~  
real,2 <W=~UUsn  
lsel,s,line,,15,17,1 bEoB;]  
lesize,all,,,1 K
'a#Mg  
lmesh,all LuQ"E4;nY%  
'Wo?%n  
!创建矩形区域 pE$|2v  
bS_y_9K  
RECTNG,0,8,0,8 nNd`]F^U  
LSEL,S,LINE,,18,21,1 $,xtif0  
LESIZE,ALL,SIZE j;$6F/g  
AMESH,1 -[i40 1  
c,+L +  
SAVE h[Ndtq>3{  
6~:W(E}  
2V#c[%vI  
!进行拉伸 z" b/osV  
C.b,]7i  
TYPE,2 %AzPAWcN  
MAT,1 Dlqn~  
EXTOPT,ESIZE,6 '-$))A
dD  
EXTOPT,aclear,1 } IIK~d,  
VOFFST,1,6 +,ZUTG  
,eZ;8W{G  
ALLSEL,ALL H5 p}Le  
nummrg,all 6M O|s1zk  
numcmp,all V)_H E  
SAVE 3ybK6!g`[  
^.ZSpc}<  
NSEL,S,LOC,Z,2,4 @&!=m]D*  
ESLN,S,1             !选择和所选节点有关的单元 JUe K"|fA  
MPCHG,2,ALL R4z<Xf:!  
ESEL,ALL >!:$@!6L  
allsel,all 94Kuy@0:+  
2GHXn:V  
NSEL,S,LOC,Z,0,2 U-wq- GT  
ESLN,S,1            !选择和所选节点有关的单元 i*mZi4URN  
MPCHG,3,ALL M63s(f  
ALLSEL,ALL OpWC2t)  
7.w*+Z>z  
esel,s,type,,3 .
E?bH V  
nsle,s,all *u:;:W&5y  
MPCHG,4,ALL chvrHvByS  
ALLSEL,ALL [ t>}SE  
4*@G&v?n  
/PNUM,MAT,1 aYv'H  
zgEr,nF  
FINISH UE}8Rkt  
!加载与求解 vkDZv@  
Jdk3) \  
/SOLU 3I(dC|d  
NLGEOM,ON          !打开大变形 EqNz L*E  
NROPT,FULL         !设定牛顿-拉弗森选项 +0oyt?  
t~+{Hr) #y  
TIME,1              !设定荷载步结束时间值 c4!c_a2pS  
AUTOTS,ON             !打开自动荷载步 0A#9C09  
NSUBST,10           !定义子荷载步为10               .Um?5wG~i  
PRED,ON           !打开非线性求解预测器 tdMP,0u  
OUTRES,ALL,ALL    !保存每一个子布的最后一个迭代计算步结果 0})7of  
s0~0
5{  
!施加约束 xI.Orpw  
NSEL,R,LOC,Y,0 {<''OwQF~+  
D,ALL,UY 4?P%M"\Iv  
ALLSEL,ALL 'xLM>6[wz  
Fi?U)T+%+  
NSEL,R,LOC,Y,B ,v$2'm)V  
D,ALL,UY lp37irI:  
ALLSEL,ALL ~#HH;q_7m  
N'lGA;}i  
NSEL,R,LOC,X,0 GFASF,+  
D,ALL,UX N(:EK  
ALLSEL,ALL X+?Il)Bv  
XwHu:v'=  
NSEL,R,LOC,X,L knNhN=hG+  
D,ALL,UX 7 K;'7  
ALLSEL,ALL }a #b$]Y  
P3,Z5|)  
NSEL,S,LOC,Z,0 .!7Fe)(x  
D,ALL,ALL X~IRpzC  
ALLSEL,ALL $M
}k%Z  
Q
2_WH)J 3  
ACEL,0,0,10 !施加重力 Ak%no3:9  
e]dPF[?7  
SAVE mG}^'?^K  
Cg]3(3   
!求解自重应力场 J]kP`  
TIME,1 m11"i=S"  
esel,s,type,,3 tu?Z@W/  
ekill,all k"3Z@Px:  
esel,all (R;) 9I\  
esel,s,live         !选择活的单元 "/ a*[_sV  
nsle,s             !选择活单元上的节点 {UV<=R,E  
nsel,invert         !反向选择，即选择了死单元上的节点 l`~a}y"n  
d,all,all,0         !将死单元上的节点约束所有位移，使其不参与矩阵运算 Lic{'w&  
nsel,all Z>>gXh<e[  
esel,all <Y}"D Yt  
eplot 8|S1|t,  
solve ?34EJ !  
FcA)RsMI*  
SAVE vy2*BTU?  
FINISH yi AG'[  
=,/A\F  
/POST1 Zh@4_Z9n!  
PLNSOL,S,Z,0,1     !SZ应力云图 }(f.uN_v  
PLNSOL,U,Z,0,1      !UZ位移云图 ]noP  
FINISH gLXvw]  
!分层开挖，开挖第1步 6ywnyh  
/SOLU !9e\O5PmO  
ANTYPE,,REST   !分析类型为静态分析,重启动分析，启动点为上次求解结束时 onWYT}c{  
TIME,2 @&i#S}%/  
AUTOTS,ON             !打开自动荷载步 pAUfG^v  
DELTIM,0.1,0.001,0.2,ON +7U A%q  
+[X.-,yW  
NSEL,S,LOC,X,4,8 'NG^HLD/  
NSEL,R,LOC,Y,0,4 a`[?,W:q  
NSEL,R,LOC,Z,3,6 (7rz:  
ESLN,S,1 |2t7G9[n  
EKILL,ALL BCa90  
VrAXOUJw6  
ALLSEL,ALL 1{\,5U&  
0,"n-5Im  
esel,s,type,,3 BM=V,BZy  
ealive,all m-Z'K_oQ  
allsel,all P0`>{!r6@  
SOLVE c1)BGy
li  
SAVE G>_42Rp  
FINISH OTNZ!U/)j  
(d5vH)+A  
!地层位移分析 Hz!
U_?  
/POST1 N>cp>&j
V  
LCDEF,1,1,, ]|((b/L3  
LCDEF,2,2,, k vue@  
LCASE,2 hX'z]Am<  
LCOPER,SUB,1,,, }e/[$!35  
_4XoUE\\  
ESEL,S,LIVE vJ'yz#tl
9  
PLNSOL,U,Z,0,1 Sug~FV?k$e  
4cErk)F4  
8zWBXV  
Yq)YS]  

我等高手出来解决问题~

没点反应啊，高手门都出来露露脸吧 =@ 
L5  

围观，学习~一起等高手

期待高手出现

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

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

不知道，问什么上传不了，图片。 ]a%\Q2[c  
2AdHj&XE  
对于，命令流没有多大问题。 R2H\;N  
只是，提出个人看法， A0N ;VYv  
)/N Xh'  
1  在初始地应力的时候，可不选取dp模型，弹性的即可。然后，初始地应力后，加dp模型。 [(e`b  
+G=C~X  
2   在前面模型材料属性 改变是在单元上。对于，改回材料3的的话，应 P".qL5  
]K0,nj*\c  
esel,s,type,,3 y
4,T  
nsle,s,all ,n,7.m.D  
MPCHG,4,ALL EK-bvZ  
u):Rw  
改为 RAx]Sp Q-S  
esel,s,type,,3 I;":O"ij\  
MPCHG,4,ALL r^o}Y  
其实，楼主，改变节点的话，没有多大问题。因为单元的计算式通过节点来算的。

qingcaozheng:上传命令流希望大家交流 ?k7/`gU  
FINISH 7zkm  
/CLEAR HjF'~n  
!(1)定义常量 EpoQV^Ey  
  L    = 8 }";hz*a  
....... (2010-04-26 17:27)  xzf)_ <  

5A=xFj{  
不知，楼主，现在还在吗？若在，希望，看一下元模型，进行一下计算。