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

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

上传命令流希望大家交流 CsW*E,|xyP  
FINISH bhc .UmH  
/CLEAR ]2'{
W]m  
!(1)定义常量 rd4\N2- 6  
  L    = 8 `B71`  
  H    = 6 h?2:'Vu]  
  B    = 8 OA\ *)c+F  
  E1    = 3.0E6 09C[B+>h  
  NU1  = 0.22 4f{(Scg  
  DEN1=2.15E3 O(Vi/r2:e  
  E2   = 1.26E7 } l4d/I  
  NU2  = 0.37 *WX,bN6Ot  
  DEN2=1.8E3 d&[.=M\E8  
  E3    =3.5E9 aBx8wl*Vm  
  NU3  = 0.25 w`F4.e  
  DEN3=2.5E3 hu''"/raM  
  SIZE = 1 x1nqhSaD  
c=A)_ZFg  
!(2)工作环境设置 z4[S02s  
/FILENAME,2010.4.26-1 %b(non*  
/TITLE,PIT OF WORKSHOP fxL0"Ry  
~LuR)T=%es  
!(3)进入前处理 4'ymPPY  
/PREP7 ~}F$1;t0  
JYU0&nZl4  
!(4)定义单元类型 kqie|_y  
ET, 1, PLANE42 I%fz^:[#<  
ET, 2, SOLID45 y:N>t+'5  
ET,3,beam4 2t7Hu)V  
R,1,0.64,0.03413,0.03413,0.8,0.8  !桩的实常数 "lJ[H=\  
R,2,0.48,0.0256,0.0144,0.6,0.8       ! 连梁的实常数 =;"$t_t  
#{u>
 
  MP, EX, 1, E1 _/zK^S)  
  MP, NUXY, 1, NU1 'dTg\ Qv  
  MP,DENS,1,DEN1 _N&]w*ce  
  TB,dp,1 m?=9j~F*  
  tbdata,1,10E3,4 rxJWU JMxK  
^@0-E@ {c  
  MP,EX,2,E2 +r 2\v  
  MP,NUXY,2,NU2 Sxw%6Va]p  
  MP,DENS,2,DEN2 hWqI*xSaJ  
  TB,dp,2 "O,TL*$  
  tbdata,1,21E3,12 Q\4nduQ  
NiTLQ"~e  
  MP,EX,3,E3 9e:}qO5)  
  MP,NUXY,3,NU3 ;'x\L<b/)  
  MP,DENS,3,DEN3 q
[w.[]  
  TB,dp,3 .^J7^Ky,  
  tbdata,1,670E3,35 t!"XQ$g'  
   MVpk/S%W  
  mp,ex,4,2.8E9      !桩、梁的弹性模量 y8.(filNB  
  mp,nuxy,4,0.3        !桩的泊松比 R0!qweGi@  
  mp,dens,4,2500       !桩的密度 ~J:"sUR  
|p1pa4%}  
7>Oa, \  
!(5)建立断面模     \x_fP;ma=_  
!Y轴坑壁 q:D!@+U  
*do,i,1,5,1           !x=40,z=30,y变 8+7=yN(  
k,,4,i-1,6 ve|`I=?2  
*end do ?m$7)@p  
.g6DKjy>  
*do,i,1,5,1           !x=40,z=11,y变 M~1 n#  
k,,4,i-1,2 x>yeF,q1  
*end do 8O5@FU 3  
'F665  
N<54_(|X  
cPgz?,hE  
0Tm"Zh?B|  
*do,i,1,5,1           !X=40,Y变的桩 ja2PmPv  
l,i,i+5 TdAHw @(  
*end do 1li1&  
type,3               !Y3 *\  
mat,4   n^7$ST#'bV  
real,1 68^5X"OGF  
lsel,s,line,,1,5,1 m%hUvG| i  
lesize,all,,,4 q3s +?&  
lmesh,all Q*+_%n1 /  
\dL#PI3  
]k(n_+!  
> CPJp!u  
*do,i,1,4,1               L8FLHT+R-  
l,i,i+1                 !连梁 gTp){  
*end do #!%\97ZR  
type,3               NI^[7.2  
mat,4   @?GOOD_i  
real,2 (HUGgX"=  
lsel,s,line,,6,9,1 Tmo+I4qoL  
lesize,all,,,1 ktr l|  
lmesh,all Hlw0ia  
,DT=(  
96.A8o
 
v&>TU(x\H  
esM<.  
(!nhU  
Ti@X<C  
3ee?B~Tun  
!x轴坑壁 Q\DD^Pbq  
y!9facg  
*do,i,2,5,1           !y=40,z=30,x变 2_x~y|<9  
k,,4+i-1,4,6 MO{6B#(<F  
*end do Ij_VO{]G'l  
B(GcPDj(K  
*do,i,2,5,1           !y=40,z=12,x变 %DQ.f*%  
k,,4+i-1,4,2 #egP*{F   
*end do o >=YoG  
4K@`>Y5g*  
Bpjw
c<U  
*do,i,1,4,1           !y=40,x变的桩 J@{yWgLg  
l,i+10,i+10+4 o'3t(dyyH  
*end do i8`&XGEd  
type,3               3huTT"G  
mat,4   J!@$lyH  
real,1 TT429  
lsel,s,line,,10,13,1 &S.zc@rN  
lesize,all,,,4 
(BgO<  
lmesh,all $h Isab_  
Z' 0Gd@/  
$?OuY*ZeY9  
h<wF;g,  
l,5,11 &pZUe`3  
type,3                     !连梁 uW&P1'X  
mat,4 9^m& [Z  
real,2 x0])&':!  
lsel,s,line,,14 8u::f`vi  
lesize,14,,,1 nD\X3g`V  
lmesh,14 ^ftZ{uA  
6N4/p=lE  
*do,i,1,3,1 ~%
4#R4&  
l,i+10,i+10+1 >mT< AQ  
*end do 9Q".166  
type,3 >sE5zj|V  
mat,4 wR;_x x  
real,2 T x_n$ &  
lsel,s,line,,15,17,1 P]Z}% 8^O  
lesize,all,,,1 vXnTPjbE  
lmesh,all K%<Z"2!+  
<!\J([NM8  
!创建矩形区域 ?Rl?Pp=>  
%aX<p{EY  
RECTNG,0,8,0,8
BPnZ"w_  
LSEL,S,LINE,,18,21,1 -v9V/LJ  
LESIZE,ALL,SIZE LB1.N!q1  
AMESH,1 uOEFb  
BZqb o`9  
SAVE *xs8/?  
;bd\XHwMUP  
.cA
[b  
!进行拉伸 47"ERfP  
vm+EzmO,!  
TYPE,2 K \}xb2s  
MAT,1 zxCxGT\;  
EXTOPT,ESIZE,6 snTj!rV/_  
EXTOPT,aclear,1 %Gn(b1X  
VOFFST,1,6 A+j~oR  
,\5]n&T;r  
ALLSEL,ALL ?-O(EY1E  
nummrg,all S ~lw5  
numcmp,all #jV6w=I  
SAVE Mi\f?  
apUV6h-v  
NSEL,S,LOC,Z,2,4 F!VC19<1O8  
ESLN,S,1             !选择和所选节点有关的单元 17G7r\iNYq  
MPCHG,2,ALL C,Je>G  
ESEL,ALL d]h[]Su/?  
allsel,all d}b#"A  
n<7#?X7  
NSEL,S,LOC,Z,0,2 \z8TYx@  
ESLN,S,1            !选择和所选节点有关的单元 xH\'gli/  
MPCHG,3,ALL \O?#gW\tR  
ALLSEL,ALL K}O~tff  
io:?JnQSA  
esel,s,type,,3 il-v>GJU7{  
nsle,s,all T7n;Bf  
MPCHG,4,ALL 9VIsLk54^  
ALLSEL,ALL WJ\YKXG  
(@`+Le  
/PNUM,MAT,1 *#EyfMz-B  
xYM!mcA  
FINISH "P< drz<  
!加载与求解 _y`'T;~OY  
C,-V>bx g  
/SOLU `c{i+  
NLGEOM,ON          !打开大变形 c*!bT$]~\  
NROPT,FULL         !设定牛顿-拉弗森选项 2:i`,  
*D]/V U  
TIME,1              !设定荷载步结束时间值 Zx5vIm  
AUTOTS,ON             !打开自动荷载步 hGF:D#jyT  
NSUBST,10           !定义子荷载步为10               lXm]1 *<  
PRED,ON           !打开非线性求解预测器 6R@ v>}  
OUTRES,ALL,ALL    !保存每一个子布的最后一个迭代计算步结果 CQsVGn{x  
hdN3r{  
!施加约束 \u,hS*v0  
NSEL,R,LOC,Y,0 f&^K>Jt1@#  
D,ALL,UY &nj@t>5Bs$  
ALLSEL,ALL hW>@jT"t1C  
Kd;
|Z  
NSEL,R,LOC,Y,B +YhTb  
D,ALL,UY 4|Jy]  
ALLSEL,ALL vK#xA+W  
fCZbIt)Eh  
NSEL,R,LOC,X,0 \rADwZm  
D,ALL,UX kv
SSz%R~  
ALLSEL,ALL 
tA@#S
Iw  
-CY?~WL&  
NSEL,R,LOC,X,L t ba%L  
D,ALL,UX f?[y-  
ALLSEL,ALL W32mAz;  
Ik=KEOz  
NSEL,S,LOC,Z,0 61T"K  
D,ALL,ALL qVJV9n  
ALLSEL,ALL IcPIOCmOc  
$9*Xfb
/  
ACEL,0,0,10 !施加重力 :!3CoC.X|c  
u&bo32fc  
SAVE S! ,.#e(Y  
EEn}Gw  
!求解自重应力场 ~|Gtm[9Ru  
TIME,1 !=cW+=1  
esel,s,type,,3 jbC7U9t7  
ekill,all HnioB=fc  
esel,all O|%><I?I  
esel,s,live         !选择活的单元 &hd+x5  
nsle,s             !选择活单元上的节点 YG?W8)T  
nsel,invert         !反向选择，即选择了死单元上的节点 5H==m~  
d,all,all,0         !将死单元上的节点约束所有位移，使其不参与矩阵运算 8Z/P<u  
nsel,all `6lc]r  
esel,all Hc^b}A y7  
eplot Uh?SDay  
solve ^N`ar9Db  
tB}&-U|t[~  
SAVE $>5|TG 0i  
FINISH %S.R@C[3  
GRO[&;d`  
/POST1 OMO.-p  
PLNSOL,S,Z,0,1     !SZ应力云图 Q?7UiTZ  
PLNSOL,U,Z,0,1      !UZ位移云图 '!{zO" 1*  
FINISH K!HSQ,AC  
!分层开挖，开挖第1步 69r<Z  
/SOLU +UHf&i/3
 
ANTYPE,,REST   !分析类型为静态分析,重启动分析，启动点为上次求解结束时 %dO'kU/-  
TIME,2 qN}0$x>p  
AUTOTS,ON             !打开自动荷载步 7qgHH p  
DELTIM,0.1,0.001,0.2,ON }I,]"0b  
R(r89bTQ  
NSEL,S,LOC,X,4,8 RjSVa.x  
NSEL,R,LOC,Y,0,4 '(&.[Pk:"  
NSEL,R,LOC,Z,3,6 : B$ d  
ESLN,S,1 GJ ZT~  
EKILL,ALL 6/.-V1*O  
#Cvjv; QwY  
ALLSEL,ALL vy1:>N?#5  
Po(9BRd7  
esel,s,type,,3 GvF~h0wMt  
ealive,all =<iK3bPkU  
allsel,all
h+CTi6-p  
SOLVE WJ=eV8Uk  
SAVE ^C_Y[i ~|  
FINISH cm6cW(x6  
y!mjZR,&  
!地层位移分析 "ln(EvW  
/POST1 VCNg`6!x  
LCDEF,1,1,, L!c7$M5xJ  
LCDEF,2,2,, ~WehG<p v[  
LCASE,2 vkASp&a  
LCOPER,SUB,1,,, cZZ-K?_  
ISa2|v;M  
ESEL,S,LIVE  9'\18_w  
PLNSOL,U,Z,0,1 :)cPc7$8  
[$;6LFs
}  
Kt;h'?  
FJp~8 x=  

我等高手出来解决问题~

没点反应啊，高手门都出来露露脸吧 O0i_h<T  

围观，学习~一起等高手

期待高手出现

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

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

不知道，问什么上传不了，图片。 lTY%,s  
! E`Tt[  
对于，命令流没有多大问题。 u/Fa+S  
只是，提出个人看法， lhw]?\  
%C\Q{_AS  
1  在初始地应力的时候，可不选取dp模型，弹性的即可。然后，初始地应力后，加dp模型。 3$ cDC8  
>ygyPl ;1s  
2   在前面模型材料属性 改变是在单元上。对于，改回材料3的的话，应 (#Z2  
J]UlCg  
esel,s,type,,3 MlgE-Lm  
nsle,s,all a/;u:"  
MPCHG,4,ALL jb~a z  
t23'x0l  
改为 0Yl4eB-  
esel,s,type,,3 lDc-W =X=  
MPCHG,4,ALL 2 PqS%`XiS  
其实，楼主，改变节点的话，没有多大问题。因为单元的计算式通过节点来算的。

qingcaozheng:上传命令流希望大家交流 %nN `|\  
FINISH ?eH&'m}-  
/CLEAR "@R>J?Cc+  
!(1)定义常量 )J]9 lW&y  
  L    = 8 2H71~~ c  
....... (2010-04-26 17:27)  z]B]QB Y[  

T>TWU:  
不知，楼主，现在还在吗？若在，希望，看一下元模型，进行一下计算。