abc.zip_ABC

;-------------------------------------------------------------------- ; 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