;--------------------------------------------------------------------
; Comparison of Hysteretic Damping with Shake91
;--------------------------------------------------------------------
title
Comparison of Hysteretic Damping with Shake91
call spec.fis suppress
config dynamic
;------------------------------------
;Grid generation and model properties
;------------------------------------
gen zone brick size 1 1 30
model mech 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 xpos mul 1.524
ini ypos mul 1.524
ini zpos mul 1.524
;-----------
; Histories
;-----------
hist add id 1 unbal
his add id 2 dytime
his add id 231 gp xacc id=121 ;top accel 0'
his add id 224 gp xacc id=89 ;accel at 40'
his add 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
def xacc_val
xacc_val=acc_mult * 9.81
;xacc_val = 9.81
end
;convert to actual accel value
apply xacc @xacc_val his table 100 dynamic range z=0
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
sav abc
return
