3D viscousDamper?

[hamiddehnavi]

hi
i have problem when i use this code. during run program say " no local y vector specified"

anyone can help me????


wipe

model Basic -ndm 3 -ndf 6


node 1 0 0 0
node 2 3 0 0
node 3 0 3 0
node 4 3 3 0
node 5 0 0 3
node 6 3 0 3
node 7 0 3 3
node 8 3 3 3
node 9 1.5 1.5 3

fixZ 0.0 1 1 1 1 1 1

rigidDiaphragm 3 9 5 6 7 8
fix 9 0 0 1 1 1 0
set column 2;
set beamx 3;
set beamy 1;
#source parameter.tcl

geomTransf Linear 1 1 0 0
geomTransf Linear 2 0 1 0
geomTransf Linear 3 0 -1 0
geomTransf Linear 4 -1 0 0
#set nu 0.15
#set h 0.7
#set massdensity 785
#source section.tcl
set E 2.00E+10
#element elasticBeamColumn $eleTag $iNode $jNode $A $E $G $J $Iy $Iz $transfTag

#Columns
element elasticBeamColumn 1 1 5 0.25 $E 7.69E+09 8.80E-03 5.21E-03 5.21E-03 $column;
element elasticBeamColumn 2 2 6 0.25 $E 7.69E+09 8.80E-03 5.21E-03 5.21E-03 $column;
element elasticBeamColumn 3 3 7 0.25 $E 7.69E+09 8.80E-03 5.21E-03 5.21E-03 $column;
element elasticBeamColumn 4 4 8 0.25 $E 7.69E+09 8.80E-03 5.21E-03 5.21E-03 $column;

#Beams
element elasticBeamColumn 5 5 6 0.16 $E 7.69E+09 3.61E-03 2.13E-03 2.13E-03 $beamx;
element elasticBeamColumn 6 7 8 0.16 $E 7.69E+09 3.61E-03 2.13E-03 2.13E-03 $beamx;
element elasticBeamColumn 7 5 7 0.16 $E 7.69E+09 3.61E-03 2.13E-03 2.13E-03 $beamy;
element elasticBeamColumn 8 6 8 0.16 $E 7.69E+09 3.61E-03 2.13E-03 2.13E-03 $beamy;


set K 102260000
set Cd 71127
set matTags1 1
uniaxialMaterial ViscousDamper $matTags1 $K $Cd $alpha

element twoNodeLink 9 1 6 -mat $matTags1 -dir 2 -doRayleigh
element twoNodeLink 10 3 8 -mat $matTags1 -dir 2 -doRayleigh
element twoNodeLink 11 2 8 -mat $matTags1 -dir 1 -doRayleigh
element twoNodeLink 12 1 7 -mat $matTags1 -dir 1 -doRayleigh