- 积分
- 4
- 注册时间
- 2005-1-21
- 仿真币
-
- 最后登录
- 1970-1-1
|
config ex 1
grid 120 130
m m
gen 0 -65 0 0 60 0 60 -65
gen arc 0 -40.11 2.20 -40.11 90
;ini x 15.1657 y -38.6065 i 32 j 54
;ini x 15.7073 y -42.91 i 32 j 45
;gen arc 10.4171 -41.39 15.7073 -42.91 46.408044
;ini x 0.8242 y -38.0702 i 3 j 55
;gen arc 7.824 -42.91 15.1657 -38.6065 114.961923
gen line 0 -37.91 0 -42.91
gen line 0 -42.91 2.2 -42.91
ini x 2.2 y -42.91 i 5 j 45
gen line 2.2 -40.11 2.2 -42.91
;gen line 2.2 -42.91 15.7764 -42.91
gen ad
gen ad
gen ad
; assign rock properties
pro bulk 1.786e10 she 1.23e10 den 2600 coh 1.8e6 fric 55
; apply boundary and initial conditions
;
fix x i 1 j 1 131
fix x i 121 j 1 131
fix x y i 1 121 j 1
set grav 9.81
ini sxx -1.686e5 var 0 1.686e5
ini szz -1.68e5 var 0 1.686e5
ini syy -1.69e6 var 0 1.69e6
;mark i 4 j 55
;unmark i 2 j 56
call ms.fis
solve
;window -10 10 -50 -30
save s20.sav
;
ini xdis 0 ydis 0
m n reg 2 50
apply xvel 0 yvel 0 from 1,45 to 1,45 long
step 1
set ninc = 1000 rstart = 1.0 rstop = 0.5
relax_ini
set ib_ap = 1 ie_ap = 6 jb_ap = 40 je_ap = 60
apply_rf
fix x i 1 j 45
fix x i 1 j 55
hist relax_hist
step 200
;solve
save s21a.sav
;
apply remove from 1,45 to 1,45 long
apply xvel 0 yvel 0 from 1,45 to 1,45 long
step 1
relax_ini
set rstart = 1.0 rstop = 0.0
apply_rf
;
stru beam beg gr 1 55 end gr 2 55 pro 1
stru beam beg gr 2 55 end gr 3 55 pro 1
stru beam beg gr 3 55 end gr 4 54 pro 1
stru beam beg gr 4 54 end gr 5 53 pro 1
stru beam beg gr 5 53 end gr 5 52 pro 1
stru beam beg gr 5 52 end gr 5 51 pro 1
stru beam beg gr 5 51 end gr 5 50 pro 1
stru beam beg gr 5 50 end gr 5 49 pro 1
stru beam beg gr 5 49 end gr 5 48 pro 1
stru beam beg gr 5 48 end gr 5 47 pro 1
stru beam beg gr 5 47 end gr 5 46 pro 1
stru beam beg gr 5 46 end gr 5 45 pro 1
;
struc prop=1 e=2.95e10 I=2.25e-3 a=0.3
struc node 13 fix r
step 100
save s21b.sav
;
;
[ 本帖最后由 benjackxu 于 2006-12-4 16:40 编辑 ] |
|