FLAC3d的用户手册中对比了FLAC3d和SHAKE91关于土层反应的计算结果，
应用滞后阻尼来表示动剪切模量与应变的关系。
其中有一点不明白，flac3d中是如何设置阻尼比与应变的关系曲线，如图所示，手册上并没有给出
相应的命令指示
还请大神指点迷津
以下是相关命令流
;-----------------------------------------------------------------------------------------------------------
new
config dynamic
;-------------------------------------
;Grid generation and model properties
;-------------------------------------
gen zone brick size 1 1 30
model elastic
prop bulk 300e6 she 186e6 den 2000  ; 1-5 ft
prop bulk 200e6 she 150e6 den 2000 range id 27 29  ;5-20 ft
prop bulk 200e6 she 168e6 den 2000 range id 25 26  ;20-30 ft
prop bulk 270e6 she 186e6 den 2000 range id 21 24  ;30-50 ft
prop bulk 350e6 she 225e6 den 2000 range id 17 20  ;50-70 ft
prop bulk 480e6 she 327e6 den 2082 range id 13 16  ;70-90 ft
prop bulk 550e6 she 379e6 den 2082 range id 9 12  ;90-110 ft
prop bulk 600e6 she 435e6 den 2082 range id 5  8  ;110-130 ft
prop bulk 750e6 she 495e6 den 2082 range id 3  4  ;130-140 ft
prop bulk 900e6 she 627e6 den 2082 range id 1  2  ;140-150 ft
;mul  a factor to have depth of 150 ft, in order to compare with Shake
ini x mul 1.524
ini y mul 1.524
ini z mul 1.524
;------------
; Histories
;------------
hist id 1 unbal
his id 2 dytime
his id 231 gp xacc id=121 ;top accel 0’
his id 224 gp xacc id=89  ;accel at 40’
his id 201 gp xacc id=1   ;bot accel at 150’
;--------------------
;Boundary Conditions
;--------------------
fix y z
;----------------------------
;Application of acceleration
;----------------------------
table 100 read Diam-flac-0001.acc
;bottom of model
apply xacc 9.81 his table 100 dynamic range z = 0  &
;convert to actual accel value
ini damp hyst default -3.325 0.823 range id 1  16
ini damp hyst default -3.156 1.904 range id 17 24
ini damp hyst default -3.325 0.823 range id 25 30
set dynamic on
set dynamic dt 0.0002
hist nstep 100
solve age 40.48
;his write 231 vs 2 file AA06.txt
;ret
his write 231 vs 2 table 231 &
;top accel hist to table in order to generate response spectra
call spec.fis
def compuspec
; setup values for fish function spectra
dmp=0.05   ; damping ratio
pmin=0.01  ; minimum period
pmax=10.0  ; maximum period
acc_in=231  ; input acc table
sd_out = 501 ; relative displacement table
sv_out = 502 ; pseudo velocity table
sa_out = 503 ; pseudo acc table
n_point = 500 ; # of computation points
spectra
end
compuspec
plot table 503 line ;pseudo-acceleration spectrum
plot table 502 line ;pseudo-velocity spectrum
;his write 231, 224, 201 vs 2 file inp-flac-0001.his

大神们都哪去了

好吧，已经弄懂了……

？？怎么回事，说说，

解释一下吧，楼主，在做相关的研究

做土动力分析时比如shake，eera等程序会用到等效线性化模型，需要两组参数——剪切模量比、阻尼比与应变的关系。在上面说的程序中，这两组数据是有一定关系，但是一般是独立的，由实验数据所得，再有一些公式去拟合，具体拟合办法可以参看《土动力学》等书，其实实际上这两组数据是有一定关系的，但是限于土动力理论理论与实验，只能把他们分开，工程上多用哈丁模型来模型滞回曲线，较为简单。
而对于flac而言，它是以阻尼的形式来描述本构关系，这里我们只需给出剪切模量比与应变的关系就可以，也就是书中提到的哈丁模型、s模型等中的参数，至于阻尼比，其实flac已经内嵌了相应的程序，所以不用自己去定义。这一点从flac手册后面有介绍，而且给出了程序以及滞回曲线图。
但是经过多次试验发现，这个滞后阻尼并不好用，他的限制条件太多，陈育民的书中也说了，所以如果单纯设置阻尼，建议不要使用滞后阻尼，除非是专门研究土动力本构关系的。多说一句，用flac做土动力，加不加阻尼都不好说，这需要看你模型后续的弹塑性发展，有时候加了阻尼反而错了。