请教前辈在UDEC里输入真实地震波的命令，谢谢了！

alexjonath

我在UDEC里已经设置了网格和初始边界条件，现在想在输入真实的地震波，现在已经有了地震波的数据（时间间隔和瞬时加速度），请教前辈将地震波导入UDEC的命令！先谢谢了！

cqm285440219

我也期待，一起解决

bingdantou

请问你解决了吗

天狼星

解决了没有，我也正头痛这个问题

jeffsand

难道就没有人会？

450934480

同求同求同求同求同求同求同求同求同求同求同求同求

safari

对于上述问题，我进行了一些汇总，但是不知道哪个语句是对的，全部写在下面，请高人告诉哪条是对的，
; set up nonreflecting boundary
bound mat=1
bound xvisc range -26 -23 -21 21
bound xvisc range 23 26 -21 21
bound xvisc yvisc range -26 26 -21 -19
bound xvisc yvisc range -26 26 19 21
; apply sinusoidal stress wave
bound stress 0 0 -1.25 yhist=cos(100.0,0.0195) range -26 26 19 21


UDEC中可以通过boundary keword history .....来施加荷载、应力或者速度等动力时程荷载，history后面可以跟关键字cosine（sine、table）等来定义施加的动力荷载时程。

例如，施加剪切应力正弦波：

bound stress 0, 0.4, 0 hist=sin(5.0,10.0) range ........




3DEC的动力输入和FLAC3D不尽相同，只有速度和应力两种方式，加速度要转换成速度或应力，且速度输入不使用于吸收边界，即粘性边界和自由场边界，而地震分析要用自由场边界，故要把速度时程转换成应力时程，手册中有个公式。
输入命令应该如下，
bound  hread 1 dizhendong.dat;读入地震动荷载
bound 0 10 -.1 .1 0 10 yvel 1.0 hist 1;施加地震动荷载
dizhendong.dat 格式和FLAC3D中的一样吧。手册dynamic一节中也有说明。
本人也是刚转入3DEC,希望和大家共同学习。尤其希望研究动力方面的同仁多多指教。

safari

这里还有一个实例，我看不懂，你们看懂了如何加载地震，请告诉我，
;call f:\a.dat
new
title
head>tangjiashan rockslide triggered by earthquake
round 0.1
set edge 8
;建立模型外边界
block 0 0 0 329 186 210 312 154 432 96 448 92 467,90 507,91 509,92 518,99 &
523,115 530,126 582,168 612,202 770,369 1086,587 1268,704 1288,719 1334,745 &
1360,756 1396,769 1537,789 1618,789 1664,779 1700,782 1800,828 1800,0
;裂隙
crack 509,92 557,99
crack 557,99 670,147
crack 670,147 803,251
crack 803,251 1112,556
crack 1112,556 1236,674
crack 1236,674 1268,704

;分界线
table 1 0 329 186 210 312 154 432 96 448 92 467,90 507,91 509,92 &
557 99 670 147 803 251 1112 556 1236 674 1302 680 1530 732 &
1675 744 1747 757 1800 794
table 2 1236 674 1268 704 1308 710 1474 754 1597 771 1690 765 &
1769 788 1800 815

;节理区域
jregion id 1 509 92 509 500 557 500 557 99
jset 50,0 1000,0 0,0 10,0  rang jregion 1
jset 140,0 1000,0 0,0 10,0  rang jregion 1

jregion id 2 557 99 557 500 670 500 670 147
jset 50,0 1000,0 0,0 10,0 rang jregion 2
jset 140,0 1000,0 0,0 10,0 rang jregion 2

jregion id 3 670 147 670 600 803 600 803 251
jset 50,0 1000,0 0,0 10,0 rang jregion 3
jset 140,0 1000,0 0,0 10,0 rang jregion 3

jregion id 4 803 251 803 700 1112 700 1112 556
jset 50,0 1000,0 0,0 10,0 rang jregion 4
jset 140,0 1000,0 0,0 10,0 rang jregion 4

jregion id 5 1112 556 1112 800 1236 800 1236 674
jset 50,0 1000,0 0,0 10,0 rang jregion 5
jset 140,0 1000,0 0,0 10,0 rang jregion 5

jregion id 6 1236 674 1236 800 1268 800 1268 704
jset 50,0 1000,0 0,0 10,0 rang jregion 6
jset 140,0 1000,0 0,0 10,0 rang jregion 6

jregion id 7 0 0 0 828 1800,828 1800,0
jset 50,0 1000,0 0,0 40,0 rang jregion 7
jset 140,0 1000,0 0,0 40,0 rang jregion 7


plot hold bl zone
print max
delete range area 5

change cons 3
change mat 1 range region 0 0 0 329 509,329 509,0
change mat 1 range region 509 0 509,92 557 99 557 0
change mat 1 range region 557 0 557 99 670 147 670 0
change mat 1 range region 670 0 670 147 803 251 803 0
change mat 1 range region 803 0 803 251 1112 556 1112 0
change mat 1 range region 1112 0 1112 556 1236 674 1236 0
change mat 1 range region 1236 0 1236 674 1302 680 1302 0
change mat 1 range region 1302 0 1302 680 1530 732 1530 0
change mat 1 range region 1530 0 1530 732 1675 744 1675 0
change mat 1 range region 1675 0 1675 744 1747 757 1747 0
change mat 1 range region 1747 0 1747 757 1784 778 1784 0
change mat 1 range region 1784 0 1784 778 1800 794 1800 0


change mat 2 range region 509,92 509 200 557 200 557 99
change mat 2 range region 557 99 557 400 670 400 670 147
change mat 2 range region 670 147 670 500 803 500 803 251
change mat 2 range region 803 251 803 700 1112 700 1112 556
change mat 2 range region 1112 556 1112 800 1236 800 1236 674
change mat 2 range region 1236 674 1236 800 1268 800 1268 704
change mat 2 range region 1236 674 1268 704 1308 710 1302 680
change mat 2 range region 1302 680 1308 710 1474 754 1530 732
change mat 2 range region 1530 732 1474 754 1597 771 1675 744
change mat 2 range region 1675 744 1597 771 1690 765 1747 757
change mat 2 range region 1747 757 1690 765 1769 788 1784 778
change mat 2 range region 1784 778 1769 788 1800 815 1800 794

ch mat 3 range region 1268 704 1268 800 1308 800 1308 710
ch mat 3 range region 1308 710 1308 800 1474 800 1474 754
ch mat 3 range region 1474 754 1474 800 1597 800 1597 771
ch mat 3 range region 1597 771 1597 800 1690 800 1690 765
ch mat 3 range region 1690 765 1690 800 1769 800 1769 788
ch mat 3 range region 1769 788 1769 830 1800 830 1800 815

;vor edge 5.0 range mat 3

pl hold bl num cont num dnum

prop mat 1 d=2650 bulk=2e10 shear=1.2e10 fric=42 coh=8.2e5 t=8e6
prop mat 2 d=2550 bulk=2e10 shear=0.8e10 fric=30 coh=1e3 t=0
prop mat 3 d=2000 bulk=2e10 shear=0.8e10 fric=30 coh=1e3 t=0

change jcons=5
change jmat=1 range 0 509 0 329
change jmat=1 range 509 557 0 90
change jmat=1 range 557 670 0 97
change jmat=1 range 670 803 0 145
change jmat=1 range 803 1112 0 249
change jmat=1 range 1112 1236 0 554
change jmat=1 range 1236 1268 0 702
change jmat=1 range 1268 1800 0 828
prop jmat=1 jkn=1e10 jks=1e10 jcoh =8e5 jfric =45 jtens =8e5


change jmat=2 range 0 509 0 329
change jmat=2 range 509 557 0 90
change jmat=2 range 557 670 0 97
change jmat=2 range 670 803 0 145
change jmat=2 range 803 1112 0 249
change jmat=2 range 1112 1236 0 554
change jmat=2 range 1236 1268 0 702
change jmat=2 range 1268 1800 0 828  
prop jmat=2 jkn=1e10 jks=1e10 jcoh =8e5 jfric =45 jtens =8e5

change jmat=3 range 509 557 90 300 angle 49,51
change jmat=3 range 557 670 97 400 angle 49,51
change jmat=3 range 670 803 145 500 angle 49,51
change jmat=3 range 803 1112 249 700 angle 49,51
change jmat=3 range 1112 1236 554 800 angle 49,51
change jmat=3 range 1236 1268 702 850 angle 49,51
prop jmat=3 jkn=1e10 jks=1e6 jfr=20 jcoh=1e3 jtens=0


change jmat=4 range 509 557 90 300 angle 139,141
change jmat=4 range 557 670 97 400 angle 139,141
change jmat=4 range 670 803 145 500 angle 139,141
change jmat=4 range 803 1112 249 700 angle 139,141
change jmat=4 range 1112 1236 554 800 angle 139,141
change jmat=4 range 1236 1268 702 850 angle 139,141
prop jmat=4 jkn=1e10 jks=1e6 jfr=20 jcoh=1e3 jtens=0

damp local
set cscan
set ovtol 100
set btol 100
set upcon 5
;gen auto 10
gen quad 20 20 range 0 509 0 329
gen quad 20 20 range 509 557 0 90
gen quad 20 20 range 557 670 0 97
gen quad 20 20 range 670 803 0 145
gen quad 20 20 range 803 1112 0 249
gen quad 20 20 range 1112 1236 0 554
gen quad 20 20 range 1236 1268 0 702
gen quad 20 20 range 1268 1800 0 828

gen edge 20 range 0 509 0 329
gen edge 20 range 509 557 0 90
gen edge 20 range 557 670 0 97
gen edge 20 range 670 803 0 145
gen edge 20 range 803 1112 0 249
gen edge 20 range 1112 1236 0 554
gen edge 20 range 1236 1268 0 702
gen edge 20 range 1268 1800 0 828


gen quad 10 10 range 509 557 90 300  
gen quad 10 10 range 557 670 97 400  
gen quad 10 10 range 670 803 145 500  
gen quad 10 10 range 803 1112 249 700  
gen quad 10 10 range 1112 1236 554 800  
gen quad 10 10 range 1236 1268 702 850  

gen edge 10 range 509 557 90 300  
gen edge 10 range 557 670 97 400  
gen edge 10 range 670 803 145 500  
gen edge 10 range 803 1112 249 700  
gen edge 10 range 1112 1236 554 800  
gen edge 10 range 1236 1268 702 850  




;gen mixed
pl bl hold zone
;set minimum contact length
;set clemin=5.0
set back iw max 5


;bou stress -0.125,0,0 range -1,1 0,329
;bou stress -0.125,0,0 range 1799,1801 0,828
;bou stress -0.125,0,0 range 0,1800 -1,1

bou xvel 0 range x -1 1
bou xvel 0 range x 1799 1801
bou yvel 0 range y -1 1
insitu  str -5e5 0 -10e5  ygrad 50 0 100 szz -10e5

print bound
plot bound xcon
plot bound ycon
grav 0 -9.8
hist unbal type 1
hist n=100 xdis 509 92 ydis 509 92
;solve fos
solve for 1e-5
hist unbal type 1
hist n=100 xdis 509 92 ydis 509 92
;solve fos
solve for 1e-5
set pl jpg 256
pl bl zone
copy f:\区域.jpg
pl bl sxy fill
copy f:\剪.jpg
pl bl sxx fill
copy f:\水平.jpg
pl bl syy fill
copy f:\竖直.jpg
pl bl disp
copy f:\位移.jpg
pl bl pla
copy f:\pla0.jpg
pl bl sig1 fill
copy f:\最大主.jpg
pl bl sig2 fill
copy f:\最小主.jpg
plot bou slip
copy f:\边界滑动.jpg
plot block shear yel
copy f:\屈服.jpg
pl bl xdisp cont
copy f:\xdisp.jpg
pl bl ydisp cont
copy f:\ydisp.jpg
pl bl slip lmag open yell
copy f:\滑动节理.jpg
;hist write 1 2ax.txt
;hist write 2 2ay.txt
;hist write 3 2bx.txt
;hist write 4 2by.txt
;hist write 5 2cx.txt
;hist write 6 2cy.txt
;hist write 7 2dx.txt
;hist write 8 2dy.txt
;hist write 9 2ex.txt
;hist write 10 2ey.txt
;hist write 11 2fx.txt
;hist write 12 2fy.txt
plot block black disp red
copy f:\位移red.jpg
save slope0.sav

rest slope0.sav
; generate free-field (340 nodes on left, 829 nodes on rgiht)
ffield gen left y 0 329 np 330
ffield gen right y 0 828 np 829
ffield change mat=1 cons=3

;ffield ini sxx -6e6 1100
;ffield ini syy -6e6 2200

; fix bottom
ffield base xvel=0
ffield base yvel=0

; top left FF
;hist ffyd 0 1
* top right FF
;hist ffyd 0 2
;hist ffsxx 0 1

set ov 100
* cycle FF (not applied to grid) to get static stresses
bou xvel 0 range x -1 1
bou xvel 0 range x 1799 1801
bou yvel 0 range y -1 1
solve for 1e-5
set pl jpg 256
pl bl sxy fill
copy f:\剪3.jpg
plot bou slip
copy f:\边界滑动3.jpg
pl bl vel
copy f:\速度3.jpg
pl bl pla
copy f:\pla.jpg
save slopeff.sav

rest slopeff.sav
; apply dynamic boundary condition
;这个地方需要修改

def wave
whilestepping
if time > env_time
wave = 0.0
else
wave = ampl * sin(2.0*pi*freq*time)
endif
end
set freq=5.0 ampl=2.5 env_time=60
;

def properties
m_dens = 2650
b_mod = 2e11
sh_mod =1.2e11
end
properties
def convert
c_p = sqrt((b_mod + (4.0 * sh_mod / 3.0)) / m_dens)
c_s = sqrt(sh_mod / m_dens)
norm_str = -2.0 * m_dens * c_p
shear_str = 2.0*(-2.0 * m_dens * c_s)
end
convert
;
boun xfree yfree
;
; bound. cond. at base of free-field
ffield base syy norm_str hist wave
ffield base sxy shear_str hist wave
ffield base yvisc xvisc
;
; apply dynamic input
bou stress (0,0,norm_str) hist wave range y -1 1
bou stress (0,shear_str,0) hist wave range y -1 1
;
;bound xvisc yvisc range y 752 753
bound xvisc range x -1 1
bound xvisc yvisc range y -1 1
bound xvisc range x 1799 1801
bound mat 1
;
; apply free-field conditions
bou ff range x -1 1
bou ff range x 1799 1801

reset time disp hist
bou xvel 0 range -1 400 -1 329
bou yvel 0 range -1 400 -1 329
bou xvel 0 range 1799 1801 -1 829
bou yvel 0 range 1799 1801 -1 829
bou yvel 0 range y -1 100

hist xvel 518 99 xdis 518 99
hist xvel 770 369 xdis 770 369
hist xvel 1220 670 xdis 1220 670
hist xvel 557 99 xdis 557 99
hist xvel 803 251 xdis 696 304
hist xvel 1189 632 xdis 1189 632

hist yvel 518 99 ydis 518 99
hist yvel 770 369 ydis 770 369
hist yvel 1220 670 ydis 1220 670
hist yvel 557 99 ydis 557 99
hist yvel 803 251 ydis 696 304
hist yvel 1189 632 ydis 1189 632



hist wave
hist unbal
; left FF
hist ffxv 329 1 ffyv 329 1
hist ffxv 200 1 ffyv 200 1
hist ffsxy 160 1 ffsyy 160 1
;
; right FF
hist ffxv 828 2 ffyv 828 2
hist ffxv 414 2 ffyv 414 2
hist ffsxy 414 2 ffsyy 414 2
;
set delc off
set ovtol 100
damp 0.001 5

;
cyc time 0.5
set pl bmp 256
pl bl vel fill
copy f:\vel1.bmp
pl bl disp fill
copy f:\disp1.bmp
pl bl pla fill
copy f:\pla1.bmp
plot block shear yel
copy f:\yelid1.bmp
save dn1.sav

cyc time 4.5
set pl bmp 256
pl bl vel
copy f:\vel2.bmp
pl bl disp
copy f:\disp2.bmp
pl bl pla fill
copy f:\pla2.bmp
plot block shear yel
copy f:\yelid2.bmp
pl bl xdisp
copy f:\xdisp1.bmp
pl bl ydisp
copy f:\ydisp1.bmp
save dn2.sav

cyc time 5
pl bl vel
copy f:\vel3.bmp
pl bl disp
copy f:\disp3.bmp
pl bl pla fill
copy f:\pla3.bmp
plot block shear yel
copy f:\yelid3.bmp
save dn3.sav

cyc time 5
set pl bmp 256
pl bl vel
copy f:\vel4.bmp
pl bl disp
copy f:\disp4.bmp
pl bl pla fill
copy f:\pla4.bmp
plot block shear yel
copy f:\yelid4.bmp
save dn4.sav

cyc time 5

pl bl vel
copy f:\vel5.bmp
pl bl disp
copy f:\disp5.bmp
pl bl pla fill
copy f:\pla5.bmp
plot block shear yel
copy f:\yelid5.bmp
save dn5.sav

cyc time 5

pl bl vel
copy f:\vel6.bmp
pl bl disp
copy f:\disp6.bmp
pl bl pla fill
copy f:\pla6.bmp
plot block shear yel
copy f:\yelid6.bmp
save dn6.sav

cyc time 5
set pl bmp 256
pl bl vel
copy f:\vel7.bmp
pl bl disp
copy f:\disp7.bmp
pl bl pla fill
copy f:\pla7.bmp
plot block shear yel
copy f:\yelid7.bmp
save dn7.sav

cyc time 5
pl bl vel
copy f:\vel8.bmp
pl bl disp
copy f:\disp8.bmp
pl bl pla fill
copy f:\pla8.bmp
plot block shear yel
copy f:\yelid8.bmp
save dn8.sav

cyc time 5
pl bl vel
copy f:\vel9.bmp
pl bl disp
copy f:\disp9.bmp
pl bl pla fill
copy f:\pla9.bmp
plot block shear yel
copy f:\yelid9.bmp
save dn9.sav

cyc time 5
pl bl vel
copy f:\vel10.bmp
pl bl disp
copy f:\disp10.bmp
save dn10.sav

cyc time 5
pl bl vel
copy f:\vel11.bmp
pl bl disp
copy f:\disp11.bmp
save dn11.sav

cyc time 5
pl bl vel
copy f:\vel12.bmp
pl bl disp
copy f:\disp12.bmp
save dn12.sav

;cyc time 5
;pl bl vel
;copy f:\vel13.bmp
;pl bl disp
;copy f:\disp13.bmp

;cyc time 5
;pl bl vel
;copy f:\vel14.bmp
;pl bl disp
;copy f:\disp14.bmp

;cyc time 5
;pl bl vel
;copy f:\vel15.bmp
;pl bl disp
;copy f:\disp15.bmp
;cyc time 5
;pl bl vel
;copy f:\vel16.bmp
;pl bl disp
;copy f:\disp16.bmp

pl bl sxy fill
copy f:\sxy.bmp

;pl bl vel iw max 5
;movie file = tjs.avi
;movie size 620 450
;movie on step 5000
;step 25000
;step 25000
;step 50000
;step 50000
;step 50000
;step 50000
;step 50000
;movie off

set pl bmp 256
pl hist 1 2 3 4 5 6
copy f:\hist1.bmp
pl hist 7 8 9 10 11 12
copy f:\hist2.bmp

save slopedy.sav

safari

还有一个例子，是咱们论坛上的，同样没看懂，你们看懂了，高速我，
;***********************************开始********************************************;
new
config dynamic
set dynamic off

;************************************建模******************************************;
poly brick -160,0 -30,145 170,280
…………简略…………
…………简略…………


;***********************************结构面切割*************************************;
…………简略…………
…………简略…………

gen edge 10
plot zol hold

;***********************************赋材料参数*************************************;
…………简略…………
…………简略…………


;***********************************初始条件***************************************;
gravity 0,-10,0

;***********************************边界条件***************************************;
bound xr   -1,1       xv=0
bound xr  -161,-159   xv=0
bound zr   169,171    zv=0
bound zr   279,281    zv=0
bound yr   -31,-29    yv=0

;***********************************初始平衡***************************************;
…………简略…………
…………简略…………

;**********************************清零********************************************;
ini xvel=0  yvel=0  zvel=0
ini xdisp=0 ydisp=0 zdisp=0
reset time disp vel hist

;***********************************动力阻尼***************************************;
set dynamic on
damp local 0.125

;***********************************监测点布置*************************************;
…………简略…………
…………简略…………

;***********************************应力时程读取***********************************;
table 1 read dzb.dat
plot hold table 1

;***********************************施加动力荷载***********************************;
ffield apply
bound  yr -31,-29   str 0 0 0 -0.8 0  0  hist table 1
bound  yr -31,-29   xvisc  yvisc  zvisc

;***********************************内嵌FISH函数***********************************;
…………简略…………
…………简略…………

;***********************************动力运算***************************************;
…………简化…………
cy ti 45
…………简化…………


后续分析简化





下面的是地震波表，
table 1
2001 0.02  
…………数据简化…………

safari

对于这个问题，兄弟我也算尽力了， 如果哪位高人知道如何加载地震，还请指教一二，

cuifrb

加载地震波说明书动力分析那一章有现成的例子啊，直接参考就行了

mining2005

运行下，看看就理解了