how modeling rotational friction damper

[payamdindar]

dear friends
i modeling rotational friction damper in opensees. i write this code in opensees for modeling rotational friction damper. i dont know these damper is rotational friction or only friction behavior. help me please?


wipe
file mkdir payam
model BasicBuilder -ndm 2 -ndf 3
set L 7.6
set H 4.6
set Fp 8645
set Tx [expr 0.636881*$Fp]
set Ty [expr 0.770962*$Fp]
set Ac 0.76
set Ic 3.4e-5
set E 2.1e10
node 1 0 0
node 2 $L 0
node 3 0 $H
node 4 $L $H
node 5 1.813 2.2036
node 6 5.787 2.2036
node 7 0 0
node 8 $L 0
node 9 3.6348 4.4
node 10 3.6348 4.4
node 11 3.9652 4.4
node 12 3.9652 4.4
node 13 3.8 4.4
node 14 3.8 4.4
node 15 3.8 4.6
node 16 3.8 4.6
mass 3 23935.8 1e-7 0
mass 4 23935.8 1e-7 0
fix 1 1 1 1
fix 2 1 1 1
geomTransf Corotational 1
element elasticBeamColumn 1 1 3 $Ac $E $Ic 1
element elasticBeamColumn 2 2 4 $Ac $E $Ic 1
element elasticBeamColumn 3 3 16 $Ac 2.1e20 $Ic 1
element elasticBeamColumn 4 16 4 $Ac 2.1e20 $Ic 1
element elasticBeamColumn 5 7 5 6.03e-4 $E 79.22e-9 1
element elasticBeamColumn 6 5 10 6.03e-4 $E 79.22e-9 1
element elasticBeamColumn 7 8 6 6.03e-4 $E 79.22e-9 1
element elasticBeamColumn 8 6 12 6.03e-4 $E 79.22e-9 1
element elasticBeamColumn 9 9 14 $Ac 2.1e20 $Ic 1
element elasticBeamColumn 10 14 11 $Ac 2.1e20 $Ic 1
element elasticBeamColumn 11 13 15 $Ac 2.1e20 $Ic 1
uniaxialMaterial Steel01 1 2200 1e20 0
element zeroLength 12 14 13 -mat 1 -dir 6
equalDOF 9 10 1 2 3
equalDOF 11 12 1 2 3
equalDOF 13 14 1 2
equalDOF 15 16 1 2 3
pattern Plain 1 Linear {
load 7 -$Tx -$Ty 0
load 8 $Tx -$Ty 0
}
set Tol 1.e-8 ;
set maxNumIter 6 ;
set printFlag 0 ;
set TestType EnergyIncr ;
test $TestType $Tol $maxNumIter $printFlag ;
system BandGeneral
numberer RCM
constraints Plain
integrator LoadControl 0.1
algorithm Linear
analysis Static
#recorder Node -file payam/dispnode2.out -node 2 -dof 1 2 3 disp
#recorder Node -file payam/dispnode3.out -node 3 -dof 1 2 3 disp
#recorder Node -file payam/reac1.out -node 1 -dof 1 2 3 reaction
#recorder Node -file payam/reac4.out -node 4 -dof 1 2 3 reaction
#recorder Element -file payam/elem1.out -time -ele 1 globalForce
#recorder Element -file payam/elem2.out -time -ele 2 globalForce
#recorder Element -file payam/elem3.out -time -ele 3 globalForce
#recorder Element -file payam/elem4.out -time -ele 4 globalForce
#recorder display 1 5 5 500 500 -wipe
#prp 5 5 10
#vup 0 1 0
#vpn 0 0 1
#display 1 10 10
#recorder plot payam/dispnode3.out 2 600 200 500 500 -columns 2 1
analyze 10
equalDOF 1 7 1 2 3
equalDOF 2 8 1 2 3
set pi [expr 2*asin(1.0)]
set a [eigen 1]
set om1 [lindex $a 0]
set period1 [expr 2*$pi/sqrt($om1)]
puts "period is $period1"
puts "Done"

sincerely

[payamdindar]

dear vesna ; HELP ME PLEASE

[vesna]

Form your example I can not see what element represents frictional damper. Can you describe me the model of your damper?

[payamdindar]

i thank you for reply.
i'm model friction damper according to what Dr. Imad H. Mualla used his papers. and i want to change this code as from friction damper to rotation friction damper.please help me to write this code.
for represented friction damper i used element zeroLength frome node 14 to 13. "element zeroLength 12 14 13 -mat 1 -dir 6"
best regards

[vesna]

You are assigning your material to rotational degree of freedom, so it is rotational damper.