dear silvia
-1-i use zero length element for modelling friction damper in brace.the problem is that the element rotate so the response of structure doesnt change with momoent frame.i constraine rotation the nodes that connect to zero length element with equal dof to a rigid node.but i saw no change and zero length element rotated again.what i can do?
-2-i cant understand from manual how can i orient the zero length element.how can i find an example or can you explain this?
thanks a lot
zero length element
Moderators: silvia, selimgunay, Moderators
Code: Select all
wipe ;
model basic -ndm 2 -ndf 3 ;
file mkdir Data ;
set g 32.18;
node 1 0 0 ;
node 2 0 26 ;
node 3 42 26 ;
node 4 42 0 ;
node 5 21 13;
node 6 21 13;
equalDOF 3 5 3 ;
equalDOF 3 6 3 ;
fix 1 1 1 1 ;
fix 4 1 1 1 ;
mass 2 75 0.0 0.0 ;
mass 3 75 0.0 0.0 ;
geomTransf Linear 1 ;
#uniaxialMaterial Steel02 $matTag $Fy $E $b $R0 $cR1 $cR2 $a1 $a2 $a3 $a4
uniaxialMaterial Steel02 1 7200 4176000 1 10 0.925 0.15 ;
#uniaxialMaterial Hysteretic $matTag $s1p $e1p $s2p $e2p $s1n $e1n $s2n $e2n $pinchX $pinchY $damage1 $damage2
uniaxialMaterial Hysteretic 2 60 1E-5 60 0.083 -60 -1E-5 -60 -0.083 1 1 0 0;
#element elasticBeamColumn $eleTag $iNode $jNode $A $E $Iz $transfTag
element elasticBeamColumn 1 1 2 0.3875 4176000 0.0911 1;
element elasticBeamColumn 2 2 3 100 4176000 1000 1;
element elasticBeamColumn 3 3 4 0.3875 4176000 0.0911 1;
element elasticBeamColumn 4 3 5 0.1465 4176000 0.0288 1;
element elasticBeamColumn 5 1 6 0.1465 4176000 0.0288 1;
#element zeroLength $eleTag $iNode $jNode -mat $matTag1 -dir $dir1 -orient $x1 $x2 $x3 $yp1 $yp2 $yp3
element zeroLength 6 6 5 -mat 2 -dir 2 -orient 0.85 -.526 0 0.526 0.85 0
recorder Node -file Data/Displacment.txt -time -node 2 -dof 1 disp ;
recorder plot Data/Displacment.txt DispHoriNode2 100 10 800 200 -columns 1 2;
recorder Element -file Data/f.txt -ele 6 force
recorder Element -file Data/d.txt -ele 6 deformation
recorder display MODEL2D 250 250 500 400 -file circle;
prp 0 0 1000;
vrp 1.5 1.5 0;
vup 0 1 0;
vpn 0 0 1;
viewWindow -5 45 -5 40 ;
plane 100000 -100000 ;
port -1 1 -1 1;
projection 1;
fill 0;
display 1 -1 50 ;
set accelSeries "Series -dt 0.02 -filePath 4.txt -factor $g";
pattern UniformExcitation 1 1 -accel $accelSeries;
rayleigh 0.95 0.001 0 0 ;
wipeAnalysis ;
constraints Plain ;
numberer Plain ;
system BandGeneral ;
test NormDispIncr 1.0e-8 10 ;
algorithm Newton ;
integrator Newmark 0.5 0.25 ;
analysis Transient ;
analyze 1280 0.02 ;