zero lenght element doubt

[yairbcn]

Could somebody please help me? I dont know why it doesnt work




# SET UP --------------------------------------------------------------------------------------


wipe;
model basic -ndm 2 -ndf 3;
file mkdir Data;

# DEFINE GEOMETRY-------------------------------------------------------------------------------


node 1 0 0;
node 2 0 0;

# BOUNDARY CONDITIONS---------------------------------------------------------------------------

fix 1 1 1 1;
fix 2 0 1 1;

# MATERIAL PROPERTIES---------------------------------------------------------------------------

set K 29000000
set alpha 1000000
set Fy 65000
set lambda 2000
set dp 0.05
set dpc 0.1

uniaxialMaterial Bilin 1 $K $alpha [ expr $alpha * (-1) ] $Fy [ expr $Fy * (-1) ] $lambda $lambda $lambda $lambda 1 1 1 1 $dp [ expr (-1) * $dp ] $dpc [ expr (-1) * $dpc ] 0.001 -0.001 $dpc [ expr (-1) * $dpc ] 1 1


# DEFINE ELEMENTS---------------------------------------------------------------------------

set transfTag 1;
geomTransf Linear $transfTag;

element zeroLength 1 1 2 -mat 1 -dir 1


# DEFINE RECORDERS------------------------------------------------------------------------------------

recorder Node -file Data/RBase1.out -time -node 1 -dof 1 reaction;
recorder Node -file Data/Disp13.out -time -node 2 -dof 1 disp;



#--------------------------------------------------------------------------------------------------------

pattern Plain 1 Linear {

load 2 10. 0. 0.
}

loadConst -time 0.0;
constraints Plain;
numberer Plain;
system BandGeneral;
test NormDispIncr 1.0e-8 6;
algorithm Newton;
integrator LoadControl 0.1;
analysis Static;
loadConst -time 0.0;
analyze 10;




#PUSHOVER------------------------------------------------------------------------------------

pattern Plain 2 Linear {

load 2 -1. 0. 0.
}

integrator DisplacementControl 2 1 0.01;
analyze 1000;




puts "Done!"