源程序是这样的 % JgRcx
new e6tH/`Uln
; create model with vertical fault stf,<W
rou 0.1 _)HD4,`
block -1 0 -1 7 9 7 9 0 zz7Y/653
cr 0 0 0 7 x`};{oz;
cr 8 0 8 7 Bf D,z
cr 0 1 8 1 }:,o Y<
cr 0 3 8 3 d ,h~u{
del 0 8 1 3 owfp^hla
cr 4 3 4 7 ZXF AuF
gen quad 1.1 range 0 8 0 7 5`1 (}
prop mat 1 d 1000 k 1.5e8 g 0.5e8 coh 2.5e5 jkn 1e8 jks 1e8 jfric 5 Oz4,Y+[#
prop mat 2 jkn 1e8 jks 1e8 jfric 0 zbF:R[ )
change cons 3 GInZ53cQ
change jmat 2 range -1 1 -1 8 o 2[vM$]
change jmat 2 range 7.5 8.5 -1 8 ` <l/GwtAJ
; apply boundary conditions |_hioMVz
fix -1 0 0 7 7yq7a[Ra
fix 8 9 0 7 >N+bU{s
bound yvel 0 xvel 0 range -1 9 -.1 .1 dv1x 78xG>
bound yvel 0 xvel 0 range -1 9 6.9 7.1 VJ~X#Q
bound (cor 631 359) yvel 0 xvel 0 )OW(T^>_'I
insitu stress 0 0 -1e5 i6!T`Kau
; apply vertical velocity to right block r` &2-]
bound (cor 840 589) yvel -2e-1 xvel 0 >?\v@
;bound yvel -2e-1 xvel 0 range 4 8 6.9 7.1 `:- @E2
; BR&Qw'O %
; add support y V 9]_k
supp 4 2 wid 3 seg 20 mat 3 yr/]xc$
prop m 3 sup_kn -1 )vEHLp.
table 1 0 0 0.2 0 0.4 0.4e7 10.0 0.4e7 ~x4{P;y
; a in#_H
step 5000 #{)r*"%
cor可能代表的是corner表示角点的意思 但是cor后面的数字是怎么得到的？ 哪位高人可以帮帮在下。。。。
