
flac3d里面出现 math error : zero divide invalid 是什么意思?怎么解决呢?
newsetlargeconfigdynsetdynoffgenzonebrickp0000p11800p200.50p3003size1813genzonebrickp...
new
set large
config dyn
set dyn off
gen zone brick p0 0 0 0 p1 18 0 0 p2 0 0.5 0 p3 0 0 3 size 18 1 3
gen zone brick p0 18 0 0 p1 20 0 0 p2 18 0.5 0 p3 18 0 3 size 3 1 3
gen zone brick p0 0 0 3 p1 18 0 3 p2 0 0.5 3 p3 0 0 13 &
p4 18 0.5 3 p5 0 0.5 13 p6 8 0 13 p7 8 0.5 13 size 18 1 10
fix x range x -0.01 0.01
fix x range x 19.99 20.01
fix z range z -0.01 0.01
fix y range y -0.01 0.01
fix y range y 0.49 0.51
model elas
prop density 2000 bulk 3e9 shear 1e9
set gravity 0 0 -10
solve
save slide.sav
restore slide.sav
sel cable id=1 begin=(14.66667,0.25,6.33333) end=(11.40770,0.25,5.37197) nseg 12
sel cable id=1 begin=(11.40770,0.25,5.37197) end=(5.35207,0.25,3.58563) nseg 4
sel cable id=2 begin=(11.33333,0.25,9.66667) end=(8.07436,0.25,8.70531) nseg 10
sel cable id=2 begin=(8.07436,0.25,8.70531) end=(3.23324,0.25,7.27723) nseg 4
sel cable id 1 prop emod 2.1e11 ytension 310e6 xcarea 0.001088 &
gr_coh 1 gr_k 1 gr_per 0.4082 range cid 1 12
sel cable id 1 prop emod 2.1e11 ytension 310e6 xcarea 0.001088 &
gr_coh 1e6 gr_k 2e7 gr_per 0.4082 range cid 13 16
sel cable id 2 prop emod 2.1e11 ytension 310e6 xcarea 0.001088 &
gr_coh 1 gr_k 1 gr_per 0.4082 range cid 17 26
sel cable id 2 prop emod 2.1e11 ytension 310e6 xcarea 0.001088 &
gr_coh 1e6 gr_k 2e7 gr_per 0.4082 range cid 27 30
sel delete link range id 1
sel delete link range id 19;;;range后面的是link的id号
sel link id=100 1 target zone;;;id后面意思是将1号节点的link号设置为100
sel link id=101 20 target zone
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 100
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 101
sel cable id=1 pretension 60e3 range cid 1,12
sel cable id=2 pretension 60e3 range cid 17,26
ini state 0
ini xdis 0 ydis 0 zdis 0
ini xvel 0 yvel 0 zvel 0
ini damp local 0.3142
model moh
prop dens 2500 bu=3.33e8 sh=1.54e8 fric=28 co=6e3
set dyn on
free x y z range z -0.01 0.01 ;去除静力边界
table 1 read x.txt
table 2 read z.txt
apply xacc 1.0 hist table 1 range z -0.01 0.01 ;施加动荷载
apply zacc 0.7 hist table 2 range z -0.01 0.01
apply ff
plot contour disp
solve 展开
set large
config dyn
set dyn off
gen zone brick p0 0 0 0 p1 18 0 0 p2 0 0.5 0 p3 0 0 3 size 18 1 3
gen zone brick p0 18 0 0 p1 20 0 0 p2 18 0.5 0 p3 18 0 3 size 3 1 3
gen zone brick p0 0 0 3 p1 18 0 3 p2 0 0.5 3 p3 0 0 13 &
p4 18 0.5 3 p5 0 0.5 13 p6 8 0 13 p7 8 0.5 13 size 18 1 10
fix x range x -0.01 0.01
fix x range x 19.99 20.01
fix z range z -0.01 0.01
fix y range y -0.01 0.01
fix y range y 0.49 0.51
model elas
prop density 2000 bulk 3e9 shear 1e9
set gravity 0 0 -10
solve
save slide.sav
restore slide.sav
sel cable id=1 begin=(14.66667,0.25,6.33333) end=(11.40770,0.25,5.37197) nseg 12
sel cable id=1 begin=(11.40770,0.25,5.37197) end=(5.35207,0.25,3.58563) nseg 4
sel cable id=2 begin=(11.33333,0.25,9.66667) end=(8.07436,0.25,8.70531) nseg 10
sel cable id=2 begin=(8.07436,0.25,8.70531) end=(3.23324,0.25,7.27723) nseg 4
sel cable id 1 prop emod 2.1e11 ytension 310e6 xcarea 0.001088 &
gr_coh 1 gr_k 1 gr_per 0.4082 range cid 1 12
sel cable id 1 prop emod 2.1e11 ytension 310e6 xcarea 0.001088 &
gr_coh 1e6 gr_k 2e7 gr_per 0.4082 range cid 13 16
sel cable id 2 prop emod 2.1e11 ytension 310e6 xcarea 0.001088 &
gr_coh 1 gr_k 1 gr_per 0.4082 range cid 17 26
sel cable id 2 prop emod 2.1e11 ytension 310e6 xcarea 0.001088 &
gr_coh 1e6 gr_k 2e7 gr_per 0.4082 range cid 27 30
sel delete link range id 1
sel delete link range id 19;;;range后面的是link的id号
sel link id=100 1 target zone;;;id后面意思是将1号节点的link号设置为100
sel link id=101 20 target zone
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 100
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 101
sel cable id=1 pretension 60e3 range cid 1,12
sel cable id=2 pretension 60e3 range cid 17,26
ini state 0
ini xdis 0 ydis 0 zdis 0
ini xvel 0 yvel 0 zvel 0
ini damp local 0.3142
model moh
prop dens 2500 bu=3.33e8 sh=1.54e8 fric=28 co=6e3
set dyn on
free x y z range z -0.01 0.01 ;去除静力边界
table 1 read x.txt
table 2 read z.txt
apply xacc 1.0 hist table 1 range z -0.01 0.01 ;施加动荷载
apply zacc 0.7 hist table 2 range z -0.01 0.01
apply ff
plot contour disp
solve 展开
- 你的回答被采纳后将获得:
- 系统奖励15(财富值+成长值)+难题奖励10(财富值+成长值)+提问者悬赏5(财富值+成长值)
若以下回答无法解决问题,邀请你更新回答
1个回答
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询