桩单元循环荷载施加

seabird08

小弟近来做桩-土的动力分析,土用实体单元,桩采用桩单元, 对于一组给定的动载如何施加到桩单元节点上,一直困扰着小弟. 看手册上说是sel node apply, 试了几次仍出错. 现把分析的命令流贴出来, 以期高手的指教


new
config dyn
title
侧向动载桩
gen zone brick size 8 8 8 edge 11
model elas
pro bulk 7.8e7 she 2.6e6  density 2000
;
;==========================================
;土体中心创建桩
sel pile id 1 beg 5.5 5.5 12.0 end 5.5 5.5 4.0 nseg 8
sel pile id 1 prop dens 2400 &
    emod 8e10 nu 0.3 xcarea 0.7854 &
    xcj 9.82e-2 xciy 4.91e-2 xciz 4.91e-2 &
    per 3.14 &
    cs_sk 1.3e11 cs_scoh 1.0e10 cs_sfric 0 &
    cs_nk 1.3e11 cs_ncoh 1.0e04 cs_nfric 0 &
    cs_ngap on
;
;设置静力边界
apply nquiet squiet dquiet ra z=(-0.01,0.01)
apply nquiet squiet dquiet ra x=(-0.01,0.01)
apply nquiet squiet dquiet ra x=(10.99,11.01)
apply nquiet squiet dquiet ra y=(-0.01,0.01)
apply nquiet squiet dquiet ra y=(10.99,11.01)
;
;定义正弦力
def setup
freq = 1.0
  omega = 2.0 * pi * freq
old_time = clock
end
setup
def wave
  wave =1000*sin(omega * dytime)
end
;
;设置阻尼
sel set damp combined
;
;
;桩头施加动荷载
   ????????????????????
;
hist id 1 sel node xdisp id 1
hist id 2 sel pilesel force fz end1 cid 1
plot add hist -2 vs 1
plot show
set dyn multi on
solve age 1.0

zxf613

为什么要设置静力边界啊 以后把命令流打包压缩后发上来 这样别人也好调试程序

seabird08

加静力边界是为了吸收桩受载时产生的动力波, 不至于从模型边界面反射回来,

qujin

楼主也是做桩土动力荷载吗，我的也是水平循环荷载，不知道你的问题解决了，我也遇到了这个问题，我的qq414052251，能否留个联系方式呀，急求，谢谢

warlm

加静力边界是为了吸收桩受载时产生的动力波, 不至于从模型边界面反射回来,
seabird08 发表于 2010-8-18 10:33

楼主可否告知联系方式，我刚好最近在做这块，我的联系方式qq156121922，谢谢了