dylan1984 发表于 2008-12-14 20:36:36

一个FLAC2D实体扩底桩的计算例子，请大家帮忙看一下

这是一个试算的例子，桩的几何尺寸就不说了。
代码如下：

;###########Properties of Pile###########          桩的参数
;E=2.55e7kN/m2    v=0.2
;calculated K=1.42e10Pa   G=1.06e10Pa

;#########Properites of Soil#########         土的参数
;E=40MPa   v=0.35    c=30kPa   fric=35    dens=1.78e3kg/m3
;calculated      K=4.4e7Pa                  G=1.48e7Pa

;#########Properties of Interface#########             这个桩土的界面参数是在计算初始地应力时候取的，比较大。取7次方时候算的话，地应力分布不均匀。
;kn 6.37e9            ks 6.37e9


config axi
grid 55 79
m mohr
m null i=5, j=51,79                     ; iface buffer zone
m null i=1,4 j=51                  ; iface buffer zone
group soil region 6,3
group pile region 2,60
;----geometry-----
gen 0,-3.8 0,0 0.4,0 0.4,-3.8i=1,5 j=61, 80             ;zone1
gen 0.4,-3.8 0.4,0 10.4,0 10.4,-3.8i=6,56 j=61,80;zone2
gen 0,-5.6 0,-3.8 0.4,-3.8 1.25,-5.6i=1,5 j=52,61         ;zone3
gen 1.25,-5.6 0.4,-3.8 10.4,-3.8 10.4,-5.6i=6,56 j=51,61 ;zone4
gen 0,-15.6 0,-5.6 1.25,-5.6 1.25,-15.6   i=1,6 j=1,51;zone5         
gen 1.25,-15.6 1.25,-5.6 10.4,-5.6 10.4,-15.6i=6,56 j=1,51 ;zone6
;---boundary conditions---
fix xi=1
fix xi=56
fix x y j=1

prop den 1.78e3 bulk 4.4e7 shear 1.48e7 co 3e4 friction 35 group soil
prop den 1.78e3 bulk 4.4e7 shear 1.48e7 co 3e4 friction 35 group pile

;----interfaces 1 2 3--------
int 1 aside from 5,80 to 5,61 bside from 6,80 to 6,61
int 2 aside from 5,61 to 5,52 bside from 6,61 to 6,51
int 3 aside from 1,52 to 5,52 bside from 1,51 to 6,51
int 1 glue kn 6.37e9 ks 6.37e9
int 2 glue kn 6.37e9 ks 6.37e9
int 3 glue kn 6.37e9 ks 6.37e9
set grav=9.81
solve elastic
save da1ini.sav
;----calculation of initial undergroud stresses----

之后在桩顶加了一个yforce 3000kN

rest da1ini.sav
;---replace property of pile----
m null group pile
m e group pile
prop dens 2500 bulk 1.42e10 shear 1.06e10 group pile

;---update properties of interfaces---
int 1 unglue kn 1e7 ks 1e7 co 8e3 fric 12   
int 2 unglue kn 1e7 ks 1e7 co 8e3 fric 12
int 3 glue kn 1e7 ks 1e7 co 8e3 fric 12

;---load---
apply yf -3e6 i=1 j=80
hist unbal

set large
;---clear previous displacement----
ini yd=0 xd=0
solve

算出的结果用plot grid mag，发现桩和土竟然有网格交差！！！？ 桩土能重合的？

然后桩顶的沉降才5mm不到， 按实测数据3000kN下，大概有将近40mm的沉降，严重不符
syy的应力云图也不对。。。很郁闷。大家帮看一下好么？！！


[ 本帖最后由 dylan1984 于 2008-12-15 14:15 编辑 ]

dylan1984 发表于 2008-12-15 10:11:26

顶一下，期待高手指点

dylan1984 发表于 2008-12-15 15:50:47

貌似论坛上很少有人用2D模拟桩啊？

天天见 发表于 2008-12-15 18:46:15

顶一下啊，呵呵

dylan1984 发表于 2008-12-16 09:45:16

没人理。。。

麦芽猪1102 发表于 2010-2-25 21:51:29

期待大虾指导
页: [1]
查看完整版本: 一个FLAC2D实体扩底桩的计算例子，请大家帮忙看一下