关于轴向加载桩的一些问题请教

冬日的阳光

new
title
轴向加载桩
ge zo b size 15 15 15 edge=11  :建立矩形的网格区域，其大小为15X15X15，网格的边长为11米
mo e ；采用各向同性弹性模型
pro bu=5e9 sh=1e9 den=2000 ；定义体积模量为5e9 Pa 剪切模量为1e9 Pa 密度为2000
fix z ra z -0.01 0.01
fix x ra x -0.01 0.01
fix x ra x 10.99 11.01
fix y ra y -0.01 0.01
fix y ra y 10.99 11.01
ini szz -2.2e5 grad 0 0 2e4
ini sxx -2.2e5 grad 0 0 2e4
ini syy -2.2e5 grad 0 0 2e4
set grav 0 0 -10
so
sel pi id=1 begin=(5.5,5.5,12.0) end=(5.5,5.5,4.0) nseg=32
sel pi pro emod=8.0e10 nu=0.3 xca=0.7854 xcj=0 xciy=0 xciz=0 per=3.14 &
    cs_sk=1.3e11 cs_sc=0 cs_sf=10 cs_nk=1.3e11 cs_nc=0 cs_nf=0 cs_ng=off
;sel delete link ra id 29
;sel link id=100 2 target zo
;sel link attach ydir=free zdir=free ra id 100
;sel link attach xdir=free ydir=free zdir=free ra id 100
;sel link attach xdir=nydeform ra id 100
;sel link constit nydeform 1 area=1 k=5.4e11 ycomp=2.22e5 ra id 100
sel node fix x ra id=1
sel node ini zvel=-0.5e-8 ra id=1
his id=1 unbal
his id=2 sel node zdisp id=1
his id=3 sel piles fo fx cid=1
plot creat his
plot add sel pi coup st sh max 5e4
plot add his -3 vs -2 xmax 12e-5 ymax 10.05e5 xlabel '节点id=1的z位移' ylabel '构件cid=1的轴向分力'
plot set magn 0.8
plot set p o 0 5.5 0 nor 0 -1 0
plot add grid p
plot add axe
plot show
sel set damp comb
cycle 20000



以上是FLAC 3D实用教程中关于轴向加载桩的例子 向请教下，在命令中添加的初始应力的值是如何计算得出的 为什么grad中的值是那个 还有其方向都是ZZ方向渐变的吗

冬日的阳光

sel node ini zvel=-0.5e-8 ra id=1
为什么初始速度要设为这个 请问有人知道吗

冬日的阳光

来人看看啊

ycding

初始速度乘计算步即为桩顶位移条件，量级小计算步多是为了增加计算曲线的稳定性，减小振荡。

冬日的阳光

4# ycding
多谢版主 那我们这个速度就是自己习惯定的一个值了是吧

ycding

差不多，主要是与计算循环数相对应。

冬日的阳光

6# ycding
多谢版主的指导