;1有哪些办法可以大大缩短动载计算时间？电脑计算要太久了，受不了。
（模型已经最简化了，滞后阻尼，静态边界。我的动载table数据一共是51.多秒， 但以前计算到1e-6秒左右就自动停了，然后我设置solve age 51.多才行，都一晚上 了dyn time才到2.011e-1秒）
;
;2 cable、liner结构能否分别与土体设置接触面？
;
;3 用seismicsignal软件进行滤波和基线校正。然后，我把底边和侧边都设为了粘滞 边界，滞后阻尼，在底部用xvel输入地震波（荷载为速度）。是否恰当？
;交流QQ747790535
附上命令流如下：
;##################################################
;静态分析
;##################################################
n
config dyn
set dyn off
;=============
;土体
;=============
gen zone brick size 9 4 8 p0=(0,0,0) p1=(0.9,0,0) p2=(0,0.32,0) p3= (0,0,1.6)
gen zone brick size 9 1 8 p0=(0,0.32,0) p1=(0.9,0.32,0) p2=(0,0.406,0) p3= (0,0.32,1.6)
gen zone brick size 9 4 8 p0=(0,0.406,0) p1=(0.9,0.406,0) p2=(0,0.726,0) p3=(0,0.406,1.6)
mo mohr
pro bulk 2.2e9 she 1.3e9 fric 30 coh 1.3e6 ten 1.5e5
ini dens 2000
;=============
;1#土钉
;=============
sel cable id=1 beg 0.2,0.1,1.3 end 0.9,0.1,1.3 nseg 14
sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 slide=on &
gr_coh 10e5 gr_k 2e7 gr_per 0.09424778 dens 1460 range cid 3,14
sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 slide=on &
gr_coh 10e5 gr_k 2e7 gr_per 0.09424778 dens 1460 range cid 1,1
sel cable prop emod 2e10 ytension 310e3 xcarea 0.0004906 dens 1460 &
gr_coh 10e8 gr_k 2e10 range cid 2,2
;=============
;2#土钉
;=============
sel cable id=2 beg 0.2,0.1,1.0 end 0.9,0.1,1.0 nseg 14
sel cable id=2 prop emod 2e10 ytension 310e3 xcarea 0.0004906 slide=on &
gr_coh 10e5 gr_k 2e7 gr_per 0.09424778 dens 1460 range cid 17,28
sel cable id=2 prop emod 2e10 ytension 310e3 xcarea 0.0004906 slide=on &
gr_coh 10e5 gr_k 2e7 gr_per 0.09424778 dens 1460 range cid 15,15
sel cable prop emod 2e10 ytension 310e3 xcarea 0.0004906 dens 1460 &
gr_coh 10e8 gr_k 2e10 range cid 16,16
;=============
;3#土钉
;=============
sel cable id=3 beg 0.2,0.1,0.7 end 0.9,0.1,0.7 nseg 14
sel cable id=3 prop emod 2e10 ytension 310e3 xcarea 0.0004906 slide=on &
gr_coh 10e5 gr_k 2e7 gr_per 0.09424778 dens 1460 range cid 31,42
sel cable id=2 prop emod 2e10 ytension 310e3 xcarea 0.0004906 slide=on &
gr_coh 10e5 gr_k 2e7 gr_per 0.09424778 dens 1460 range cid 29,29
sel cable prop emod 2e10 ytension 310e3 xcarea 0.0004906 dens 1460 &
gr_coh 10e8 gr_k 2e10 range cid 30,30
;=============
;4#土钉
;=============
sel cable id=4 beg 0.3,0.1,0.4 end 0.9,0.1,0.4 nseg 12
sel cable id=4 prop emod 2e10 ytension 310e3 xcarea 0.0004906 slide=on &
gr_coh 10e5 gr_k 2e7 gr_per 0.09424778 dens 1460 range cid 45,54
sel cable id=2 prop emod 2e10 ytension 310e3 xcarea 0.0004906 slide=on &
gr_coh 10e5 gr_k 2e7 gr_per 0.09424778 dens 1460 range cid 43,43
sel cable prop emod 2e10 ytension 310e3 xcarea 0.0004906 dens 1460 &
gr_coh 10e8 gr_k 2e10 range cid 44,44
;=============
;5#土钉
;=============
sel cable id=5 beg 0.3,0.1,0.1 end 0.9,0.1,0.1 nseg 12
sel cable id=5 prop emod 2e10 ytension 310e3 xcarea 0.0004906 slide=on &
gr_coh 10e5 gr_k 2e7 gr_per 0.09424778 dens 1460 range cid 57,66
sel cable id=2 prop emod 2e10 ytension 310e3 xcarea 0.0004906 slide=on &
gr_coh 10e5 gr_k 2e7 gr_per 0.09424778 dens 1460 range cid 55,55
sel cable prop emod 2e10 ytension 310e3 xcarea 0.0004906 dens 1460 &
gr_coh 10e8 gr_k 2e10 range cid 56,56
;=============
;面板
;=============
sel liner range y=0,0.726 x=0.85,0.9 z=0,1.6
sel node local xdir 1 0 0 ydir 0 -1 0 range x 0.85 0.9
sel node fix lsys range x 0.85 0.9
sel node fix x yr zr range x 0.85 0.9
sel node local xdir 0 0 -1 ydir 0 -1 0 range y 0 0.726
sel node fix lsys range y 0 0.726
sel node fix y xr zr range y 0 0.726
sel liner PROP iso=( 35e9, 0.15) thick=0.1 dens 2500 ; concrete
sel liner PROP cs_nk=8e8 cs_sk=8e8 slide=on &
cs_ncut=0.0 cs_scoh=0.0 cs_scohres=0.0 cs_sfric=0.0
;=============
;静力边界条件
;=============
fix x y range x -0.001,0.001
fix y x range y -0.001,0.001
fix y x range y 0.725,0.727
fix x y z range z -0.001 0.001
;=============
;施加重力并求解保存
;=============
set grav 0 0 -10
solve
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
save jingtaifenxi998.sav
;##################################################
;动载分析
;##################################################
n
rest jingtaifenxi998.sav
;=============
;动力计算
;=============
set dyn=on, hist_rep=100, dyn time=0
set large
range name bottom z=-0.001 0.001
free x y z range bottom ;去掉底部原有静力条件
apply nquiet squiet dquiet range bottom ;静态边界条件
range name leftboundary x=-0.001 0.001
free x y z range leftboundary ;去掉左边界原有静力条件
apply nquiet squiet dquiet range leftboundary ;设置静态边界条件
range name rightboundary y=0.725 0.727
free x y z range rightboundary ;去掉右边界原有静力条件
apply nquiet squiet dquiet range rightboundary ;设置静态边界条件
range name backboundary y=-0.001 0.001
free x y z range backboundary ;去掉背面边界原有静力条件
apply nquiet squiet dquiet range backboundary ;设置静态边界条件
call table0.1g.txt
apply xvel 1.0 hist table 1 range bottom ;施加动荷载
ini damp hyst default -3.325 0.823 ;设置滞后阻尼
;=============
;面板加速度监测
;=============
hist id=1 gp xacc 0.9,0,0
hist id=2 gp xacc 0.9,0.363,0.38
hist id=3 gp xacc 0.9,0.363,0.87
hist id=4 gp xacc 0.9,0.363,1.4
;=============
;面板位移监测
;=============
hist id=5 gp xdisp 0.9,0.363,0.2
hist id=6 gp xdisp 0.9,0.363,0.66
hist id=7 gp xdisp 0.9,0.363,1.0
hist id=8 gp xdisp 0.9,0.363,1.3
hist id=9 gp xdisp 0.9,0.363,1.42
;=============
;沉降监测
;=============
hist id=10 gp zdisp 0.2,0.363,1.4
hist id=11 gp zdisp 0.4,0.363,0.8
hist id=12 gp zdisp 0.6,0.363,0.5
hist id=13 gp zdisp 0.8,0.363,0.2
;=============
;土压力监测
;=============
hist id=14 zone sxx 0.85,0.313,1.4
hist id=15 zone sxx 0.85,0.313,1.1
hist id=16 zone sxx 0.85,0.313,0.8
hist id=17 zone sxx 0.85,0.313,0.5
hist id=18 zone sxx 0.85,0.313,0.2
hist id=19 zone szz 0.85,0.313,1.4
hist id=20 zone szz 0.85,0.313,1.1
hist id=21 zone szz 0.85,0.313,0.8
hist id=22 zone szz 0.85,0.313,0.5
hist id=23 zone szz 0.85,0.313,0.2
hist id=24 dytime
plot create hhh
plot add hist 1 2 3 4 vs 24
plot show
set log on
set logfile 0.1ggongkuang1.log
set log off
solve age 51.1625
save 0.1ggongkuang1.sav
;=============
;数据输出
;=============
hist write 1 2 3 4 vs 24 file 0.1gjiasudu_thist.txt
hist write 1 2 3 4 file 0.1gjiasuduhist.txt
hist write 5 6 7 8 9 10 11 12 13 vs 24 file 0.1gweiyi_thist.txt
hist write 5 6 7 8 9 10 11 12 13 file 0.1gweiyihist.txt
hist write 14 15 16 17 18 19 20 21 22 23 vs 24 & file 0.1gtuyali_thist.txt
hist write 14 15 16 17 18 19 20 21 22 23 & file 0.1gtuyali_thist.txt
