没有合适的资源?快使用搜索试试~ 我知道了~
陈育民《FLAC3D基础与工程实例》全部命令流
需积分: 50 54 下载量 51 浏览量
2018-10-18
17:13:02
上传
评论 4
收藏 264KB DOC 举报
温馨提示


试读
64页
FLAC 是国际通用的岩土工程专业分析软件,具有强大的计算功能和广泛的模拟能力, 尤其在大变形问题的分析方面具有独特的优势。软件提供的针对于岩土体和支护体系的各种 本构模型和结构单元更突出了 FLAC 的“专业”特性,因此在国际岩土工程界非常流行。在 国内,FLAC 的应用也日渐广泛,拥有越来越多的用户群,而相关学习材料还不够丰富,该 书为广大的 FLAC 使用者提供了新的参考素材。 全
资源推荐
资源详情
资源评论








《FLAC/FLAC 3D 基础与工程实例》全部命令流
1
gen zone bri p0 0 0 0 p1 10 0 0 p2 0 10 0 p3 0 0 10 &
p4 15 15 0 p5 0 15 15 p6 15 0 10 p7 20 20 20 &
size 10 10 10 rat 1.0 0.9 1.1 group brick_1
gen zone bri p0 20 0 0 p1 add 10 0 0 p2 add 0 20 0 p3 add 0 0 15 &
size 10 10 10 rat 1.0 0.9 1.1 group brick_2
gen zone bri p0 40 0 0 edge 10 size 10 10 10 rat 1.0 0.9 1.1 group brick_3
plot sur

2-1
new
gen zon bri size 3 3 3
model elas
prop bulk 3e8 shear 1e8
ini dens 2000
fix z ran z -.1 .1
fix x ran x -.1 .1
fix x ran x 2.9 3.1
fix y ran y -.1 .1
fix y ran y 2.9 3.1
set grav 0 0 -10
solve
app nstress -10e4 ran z 3 x 1 2 y 1 2
hist gp vel 0 0 3
hist gp vel 0 3 3
plo hist 1 red
plo add hist 2 blue
solve

3-1
;-------------------------------------工程信息
;Project Record Tree export
;Title:Simple test
;---------------------------------计算第一步
;... STATE: STATE1 ....
config
grid 10,10
model elastic
group 'User:Soil' notnull
model elastic notnull group 'User:Soil'
prop density=1500.0 bulk=3E6 shear=1E6 notnull group 'User:Soil'
fix x y j 1
fix x i 1
fix x i 11
set gravity=9.81
history 999 unbalanced
solve
save state1.sav
;----------------------------------计算第二步
;... STATE: STATE2 ....
initial xdisp 0 ydisp 0
initial xvel 0 yvel 0
model null i 4 7 j 8 10
group 'null' i 4 7 j 8 10
group delete 'null'
history 1 xdisp i=4, j=11
solve
save state2.sav
;--------------------------------绘图命令
;*** plot commands ****
;plot name: syy
plot hold grid syy fill
;plot name: Unbalanced force
plot hold history 999
;plot name: grid
plot hold grid magnify 20.0 lred grid displacement
;plot name: Xdis-A
plot hold history 1 line

5-1
new
; ===============================
; 定义球体半径和半径方向上单元网格数
; ===============================
def parm
rad=10.0
rad_size=5
end
parm
; ===============================
; 建立八分之一球体外接立方体网格
; ===============================
gen zone pyramid p0 rad 0 0 p1 rad 0 rad p2 rad rad 0 p3 0 0 0 &
p4 rad rad rad size rad_size rad_size rad_size group 1
gen zone pyramid p0 0 rad 0 p1 rad rad 0 p2 0 rad rad p3 0 0 0 &
p4 rad rad rad size rad_size rad_size rad_size group 2
gen zone pyramid p0 0 0 rad p1 0 rad rad p2 rad 0 rad p3 0 0 0 &
p4 rad rad rad size rad_size rad_size rad_size group 3
; ==================================
; 利用 FISH 语言将内部立方体节点调整到球面
; ==================================
def make_sphere
p_gp=gp_head
loop while p_gp#null
; 获取节点点坐标值:P=(px,py,pz)
px=gp_xpos(p_gp)
py=gp_ypos(p_gp)
pz=gp_zpos(p_gp)
dist=sqrt(px*px+py*py+pz*pz)
if dist>0 then
; 节点位置调整
maxp=max(px,max(py,pz))
k=(maxp/rad)*(rad/dist)
gp_xpos(p_gp)=k*px
gp_ypos(p_gp)=k*py
gp_zpos(p_gp)=k*pz
end_if
p_gp=gp_next(p_gp)
end_loop
end
make_sphere
; ===============================
; 利用镜像生成完整球体网格

; ===============================
gen zone ref
gen zone ref dip 90
gen zone ref dip 90 dd 90
; ===============================
; 显示球体网格
; ===============================
plot surf
pl set back wh
pl bl gr
剩余63页未读,继续阅读
资源评论


qq_43457373
- 粉丝: 0
- 资源: 1
上传资源 快速赚钱
我的内容管理 展开
我的资源 快来上传第一个资源
我的收益
登录查看自己的收益我的积分 登录查看自己的积分
我的C币 登录后查看C币余额
我的收藏
我的下载
下载帮助


安全验证
文档复制为VIP权益,开通VIP直接复制
