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

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

上传命令流希望大家交流 "2h5m4  
FINISH 0~xaUM`  
/CLEAR h3J*1  
!(1)定义常量 I=YCQ VvA  
  L    = 8 |vy]8?Ak  
  H    = 6 "d?f:x3v^  
  B    = 8 ;]h.m)~|  
  E1    = 3.0E6 7b.U!Ju  
  NU1  = 0.22 ,L-C(j  
  DEN1=2.15E3 0|2%# E  
  E2   = 1.26E7 ^`&HWp  
  NU2  = 0.37 i X/tt  
  DEN2=1.8E3 |t\KsW  
  E3    =3.5E9 ",Wf uz  
  NU3  = 0.25 [:i sZG*  
  DEN3=2.5E3 e?)yb^7K  
  SIZE = 1 R^9"N?Q7;`  
 nhfwOS  
!(2)工作环境设置 ,o&<WMD  
/FILENAME,2010.4.26-1 F7uhuqA]N  
/TITLE,PIT OF WORKSHOP ?"@ET9  
T|f_~#?eV  
!(3)进入前处理 ~t9$IB  
/PREP7 P`sN&Y~m  
P,1exgq9  
!(4)定义单元类型 ;4d.)-<No_  
ET, 1, PLANE42 o5#,\Y[ g  
ET, 2, SOLID45 *IlQ5+3I  
ET,3,beam4 q[boWW  
R,1,0.64,0.03413,0.03413,0.8,0.8  !桩的实常数 -fM1$/]  
R,2,0.48,0.0256,0.0144,0.6,0.8       ! 连梁的实常数 ZA.fa0n
  
}W "(cYN_  
  MP, EX, 1, E1 z
\[(g  
  MP, NUXY, 1, NU1 h}6b&m  
  MP,DENS,1,DEN1 `2x34  
  TB,dp,1 fH!=Zb_{8  
  tbdata,1,10E3,4 , ;'SVe%  
a R#Cot  
  MP,EX,2,E2 ct\<;I(H  
  MP,NUXY,2,NU2 Ck(.N  
  MP,DENS,2,DEN2 %)IrXz>Zh  
  TB,dp,2 v,\93mNp[  
  tbdata,1,21E3,12 mcMb*?]  
k E_ky)  
  MP,EX,3,E3 #D:RhqjK  
  MP,NUXY,3,NU3 ry,}F@P&  
  MP,DENS,3,DEN3 |!re8|JV_  
  TB,dp,3 yGV{^?yoP  
  tbdata,1,670E3,35 \|!gPc%s  
   X'2Gi  
  mp,ex,4,2.8E9      !桩、梁的弹性模量 L;C|ow^c  
  mp,nuxy,4,0.3        !桩的泊松比 JfKg_&hM  
  mp,dens,4,2500       !桩的密度 _z:Qhe  
s<{GpWT8  
zZYH
c?Z  
!(5)建立断面模     ^@ UjQ9[>  
!Y轴坑壁 -ddOh<U>  
*do,i,1,5,1           !x=40,z=30,y变 <t6d)mJ%  
k,,4,i-1,6 >H?uuzi  
*end do &9h  
w$% BlqN  
*do,i,1,5,1           !x=40,z=11,y变 Bi_J5 If  
k,,4,i-1,2 Wy)('EM  
*end do 9&(.x8d,a  
YnxU
(v'\  
,'-?:`hP'  
NhtEW0xCr  
pU[K%@sC  
*do,i,1,5,1           !X=40,Y变的桩 J_/05(48  
l,i,i+5 c+;S<g0  
*end do 7S2c|U4IM  
type,3               jmPp-}tS7  
mat,4   N K"%DU<  
real,1 NlFo$Y  
lsel,s,line,,1,5,1 [Ye5Y?  
lesize,all,,,4 a&:>Ped"  
lmesh,all QAXYrRu  
rHo6iJj  
7+S44)w}~  
]1Qi=2'  
*do,i,1,4,1               Lnx2xoNk  
l,i,i+1                 !连梁 ;5RIwD  
*end do iJ*%dio  
type,3               ZW2s[p r  
mat,4   q+J0}y{#8)  
real,2 [5LMt*Y  
lsel,s,line,,6,9,1 Fs9W>*(  
lesize,all,,,1 tD^a5qPh  
lmesh,all #,Bj!'Q'-  
^HoJ.oC/  

3APYO  
coO.kTO;  
o?=fhc  
</hR!Sb]  
RD9Yk  
O &\<FT5  
!x轴坑壁 &O{t
^D)F  
qqD0R*(C  
*do,i,2,5,1           !y=40,z=30,x变 d:3= 1x  
k,,4+i-1,4,6 2_Jb9:/X  
*end do <|dj^.^  
DD6'M U4  
*do,i,2,5,1           !y=40,z=12,x变 J<-Fua^  
k,,4+i-1,4,2 A xR\ned  
*end do WV~SL/k|  
\"bLE0~  
HtS#_y%(  
*do,i,1,4,1           !y=40,x变的桩 }JJ::*W2n  
l,i+10,i+10+4 4i96UvkZ  
*end do DzmqR0)  
type,3               q]?+By-0  
mat,4   `L">"V`$Bj  
real,1 [R$liN99z;  
lsel,s,line,,10,13,1 /]l f>\x1  
lesize,all,,,4 &0h=4i=6r  
lmesh,all s|p(KWo2U  
p(7c33SyF  
Wlxk  
x[a'(5PwY  
l,5,11 5YLho2h38!  
type,3                     !连梁 1Y2a*J  
mat,4 DUPmq!A  
real,2 M->Kz{h?j  
lsel,s,line,,14 `~KAk  
lesize,14,,,1 `0[fLEm  
lmesh,14 wJr/FE7c  
a~F u  
*do,i,1,3,1 OMxxI6h  
l,i+10,i+10+1 C0f[eA  
*end do p|;o5j{  
type,3 TQ2i{e
 
mat,4 SOYDp;j  
real,2 $WM8tF?H  
lsel,s,line,,15,17,1 Vg) ^|  
lesize,all,,,1 YL!oF^XO  
lmesh,all 6<Be#Y]b  
*q[^Q'jnN  
!创建矩形区域 6%z`)d  
Y/!0Q6<[2Y  
RECTNG,0,8,0,8 rOhA*_EG  
LSEL,S,LINE,,18,21,1 iQ0&W0D]  
LESIZE,ALL,SIZE nO%<;-=u\  
AMESH,1  z8tt+AU  

kz|[
*%10  
SAVE !?Tzk&'  
)rS^F<C  
3_@G{O)e  
!进行拉伸 .>oM z&  
.1%i`+uZ  
TYPE,2 3?]S,~!F  
MAT,1 MSFNw  
EXTOPT,ESIZE,6 I@c0N*(  
EXTOPT,aclear,1 /^8t'Jjd,  
VOFFST,1,6 n_ OUWvs  
0Mq6yu^  
ALLSEL,ALL `C ?a  
nummrg,all 0YHYxn  
numcmp,all Cb<~i  
SAVE 3dY6;/s  
{vCU^BN,k  
NSEL,S,LOC,Z,2,4 !zl/0o  
ESLN,S,1             !选择和所选节点有关的单元 V?o&])?[  
MPCHG,2,ALL "9.6\Y\*  
ESEL,ALL L;kyAX@^  
allsel,all ~v,!n/('  
<|wmjW/D  
NSEL,S,LOC,Z,0,2 /.f!  
ESLN,S,1            !选择和所选节点有关的单元  MbM:3  
MPCHG,3,ALL ?~]>H A:  
ALLSEL,ALL gFgcxe6  
}"g@E-]N  
esel,s,type,,3 H.f9d.<W%  
nsle,s,all dfXV1B5  
MPCHG,4,ALL g')?J<z  
ALLSEL,ALL 2voNg
Y  
8Y]u:v  
/PNUM,MAT,1 "%Lmgy:~  
w`"W3(  
FINISH ^r%i3  
!加载与求解 W|FNDP0  
vatx+)  
/SOLU )8\Z=uC  
NLGEOM,ON          !打开大变形 lTd+{TF.  
NROPT,FULL         !设定牛顿-拉弗森选项 yW'BrTw  
d!FONi  
TIME,1              !设定荷载步结束时间值 EywZIw?mjX  
AUTOTS,ON             !打开自动荷载步 kv[OW"8t  
NSUBST,10           !定义子荷载步为10               rHR5,N:  
PRED,ON           !打开非线性求解预测器 Psg +\14  
OUTRES,ALL,ALL    !保存每一个子布的最后一个迭代计算步结果 ^S3A10f,  
N/`g?B[  
!施加约束 X{4xm,B/  
NSEL,R,LOC,Y,0 ..]B9M.  
D,ALL,UY ta2z  
ALLSEL,ALL c '/2F0y  
s@R3#"I  
NSEL,R,LOC,Y,B b<48#Qy~l  
D,ALL,UY F'fM?!(  
ALLSEL,ALL pvI&-D #}  
yFa&GxSq  
NSEL,R,LOC,X,0 '$lw[1  
D,ALL,UX s'_$j$1  
ALLSEL,ALL d9ZDpzxB  
"F04c|oR<X  
NSEL,R,LOC,X,L "P9wT)J_  
D,ALL,UX mT}Aje-L  
ALLSEL,ALL xU:PhhS  
v UJ sFR  
NSEL,S,LOC,Z,0 fdvi}SS8  
D,ALL,ALL chMt5L+5  
ALLSEL,ALL pZW}^kg=  
69[w/\  
ACEL,0,0,10 !施加重力 s 0 =@ &/  
`z5v}T  
SAVE )bF)RLZ  
mq>*W'M  
!求解自重应力场 N7%=K9  
TIME,1 -_:
JQ  
esel,s,type,,3 d8 3+6d  
ekill,all (d1V1t2r6  
esel,all _dz:\v  
esel,s,live         !选择活的单元 cM|af#o  
nsle,s             !选择活单元上的节点 g
DgP;id  
nsel,invert         !反向选择，即选择了死单元上的节点 06Sqn3MB  
d,all,all,0         !将死单元上的节点约束所有位移，使其不参与矩阵运算 CA'hvXb.  
nsel,all 6eVe}V4W  
esel,all ZD iW72&Q  
eplot r(748Qc4f?  
solve }P7xdQ6  
,2Sv1v$  
SAVE +*]SP@|IYI  
FINISH \0xzBs1!  
HP1X\h!Ke  
/POST1 %Td+J`|U+  
PLNSOL,S,Z,0,1     !SZ应力云图 h%4~0  
PLNSOL,U,Z,0,1      !UZ位移云图 oo"JMD)  
FINISH ^2(";.m  
!分层开挖，开挖第1步 ntd ":BKi  
/SOLU Ykx&6M@t  
ANTYPE,,REST   !分析类型为静态分析,重启动分析，启动点为上次求解结束时 Nj"_sA p  
TIME,2 )yfOrsM  
AUTOTS,ON             !打开自动荷载步 ZzSJm+&'  
DELTIM,0.1,0.001,0.2,ON >0[qi1  
`1DU b7<  
NSEL,S,LOC,X,4,8 &L2`L)  
NSEL,R,LOC,Y,0,4 *T1L)Cp  
NSEL,R,LOC,Z,3,6 T749@!v`z  
ESLN,S,1 9$}+-Z  
EKILL,ALL 9p!dQx  
axt6u)4%7:  
ALLSEL,ALL
5LnB]dW  
c1=;W$T(s  
esel,s,type,,3 Q
q6%53  
ealive,all a .B\=3xn  
allsel,all a2IV!0x  
SOLVE PLlx~A  
SAVE L|vaTidc0  
FINISH #nt<j2}m  
9QO!vx  
!地层位移分析 -uS7~Ww.a  
/POST1 a?f5(qW3  
LCDEF,1,1,, e{d_p%(  
LCDEF,2,2,, B]CS2LEqh  
LCASE,2 'bd=,QW  
LCOPER,SUB,1,,, o
%QhV6(F  
7~QwlU3n<F  
ESEL,S,LIVE ,5%aP%  
PLNSOL,U,Z,0,1 zcbA)
  
V1AEjh  
6+5(.z-[  
4{1c7g  

我等高手出来解决问题~

没点反应啊，高手门都出来露露脸吧 4 g/<).1<b  

围观，学习~一起等高手

期待高手出现

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

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

不知道，问什么上传不了，图片。 bx%hizb  
j)SgB7Q  
对于，命令流没有多大问题。 Z-;uzx  
只是，提出个人看法， n!EH>'T  
n?ZH2dI\0  
1  在初始地应力的时候，可不选取dp模型，弹性的即可。然后，初始地应力后，加dp模型。 3:CQMZ|;@  

Ozh^Q$>u  
2   在前面模型材料属性 改变是在单元上。对于，改回材料3的的话，应 ;zxlwdfcr'  
LMhY"/hAXa  
esel,s,type,,3 E.Gh@i  
nsle,s,all j#.-MfB  
MPCHG,4,ALL MX< ($M  
iTtAj~dfZ  
改为 emDvy2uA#  
esel,s,type,,3 Vzv.e
6_  
MPCHG,4,ALL x4 A TK  
其实，楼主，改变节点的话，没有多大问题。因为单元的计算式通过节点来算的。

qingcaozheng:上传命令流希望大家交流 a\2Myj  
FINISH -QNMB4  
/CLEAR *#{.\R-D  
!(1)定义常量 rF j)5~  
  L    = 8 3B_} :  
....... (2010-04-26 17:27)  >p+gx,N  

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