$ fc' = 4415 psi $ 0.2" diameter bar $ If EXP=0, rebar force reaches GB*As*SMAX and does not decay $ If EXP>0, rebar force reaches GB*As*SMAX and decays as rebar slips $ As is 2*pi*radius*length of bonded bar *KEYWORD *damping_part_stiffness 1,0.1 *contact_1d $, , ERR, SIGC, GB, SMAX, EXP $1,2, 0.1, 4415,10000, 0.01, 0.0 1,2, 0.1, 4415,10000, 0.01, 10.0 *set_node_list $ rebar nodes along bondline 1 3,4,7,10 *set_node_list $ concrete nodes along bondline 2 22,31,40 *database_bndout 4e-5 *boundary_prescribed_motion_node 16,1,0,2,5.0 *define_curve 2 0,0 0.002,1 1,1 $*load_node_point $16,1,1,1.0 *define_curve 1 0,0 1,10000 *CONTROL_ENERGY $# hgen rwen slnten rylen 2 2 2 2 *CONTROL_TERMINATION $# endtim endcyc dtmin endeng endmas 0.008000 *CONTROL_TIMESTEP ,.3 *DATABASE_GLSTAT $# dt binary 8e-5 *DATABASE_BINARY_D3PLOT $# dt lcdt beam npltc 8e-5 *ELEMENT_SOLID 1 2 18 27 30 21 19 28 31 22 6 2 21 30 33 24 22 31 34 25 11 2 19 28 31 22 20 29 32 23 16 2 22 31 34 25 23 32 35 26 *ELEMENT_SOLID 2 2 27 36 39 30 28 37 40 31 7 2 30 39 42 33 31 40 43 34 12 2 28 37 40 31 29 38 41 32 17 2 31 40 43 34 32 41 44 35 *ELEMENT_BEAM $# eid pid n1 n2 n3 rt1 rr1 rt2 rr2 local 1 1 3 4 5 2 1 4 7 8 3 1 7 10 11 4 1 10 13 14 5 1 13 16 17 *NODE $# nid x y z tc rc 3 0.000 2.00000000 2.00000000 4 0.8000 2.00000000 2.00000000 5 0.000 4.00000000 2.00000000 7 1.6000 2.00000000 2.00000000 8 1.2000 4.00000000 2.00000000 10 2.4000 2.00000000 2.00000000 11 2.4000 4.00000000 2.00000000 13 3.2000 2.00000000 2.00000000 14 3.6000 4.00000000 2.00000000 16 4.00000000 2.00000000 2.00000000 17 4.8000 4.00000000 2.00000000 18 0.000 0.000 0.000 1 19 0.000 0.000 2.00000000 1 20 0.000 0.000 4.00000000 1 21 0.000 2.00000000 0.000 1 22 0.000 2.00000000 2.00000000 1 23 0.000 2.00000000 4.00000000 1 24 0.000 4.00000000 0.000 1 25 0.000 4.00000000 2.00000000 1 26 0.000 4.00000000 4.00000000 1 27 1.2000 28 1.2000 0.000 2.00000000 29 1.2000 0.000 4.00000000 30 1.2000 2.00000000 31 1.2000 2.00000000 2.00000000 32 1.2000 2.00000000 4.00000000 33 1.2000 4.00000000 34 1.2000 4.00000000 2.00000000 35 1.2000 4.00000000 4.00000000 36 2.400 37 2.400 0.000 2.00000000 38 2.400 0.000 4.00000000 39 2.400 2.00000000 40 2.400 2.00000000 2.00000000 41 2.40 2.00000000 4.00000000 42 2.400 4.00000000 43 2.400 4.00000000 2.00000000 44 2.400 4.00000000 4.00000000 *PART rebar 1 1 1 *SECTION_BEAM $$ diameter of bar = 0.2 inch $# secid elform shrf qr/irid cst scoor nsm $ 1 1 1 $0.2,0.2, 0,0 $$ truss 1,3 0.031416 *MAT_ELASTIC $# mid ro e pr da db not used 1 7.3000E-4 3.0000E+7 0.300000 *PART concrete $# pid secid mid eosid hgid grav adpopt tmid 2 2 2 0 2 *SECTION_SOLID $# secid elform aet 2 $*MAT_ELASTIC $ 2 2.3000E-4 3.0000E+5 0.400000 $ Units: English *mat_cscm_concrete 2,2.3000E-4, 1 4415,0.75,3 $ Units: tonnes, mm, sec, N, MPa $*MAT_CSCM_CONCRETE_TITLE $f'c =4415 psi $ 142.3200E-09 1 0.0 1 1.05 1.0 0 $ 0.0 $30,19,2 *HOURGLASS 2,4,.05 *END