Ignore the material in series

[Gholamreza]

Hello OpenSees Community,
I added new material and could compile and run it successfully and get reasonable results for a SDOF model.
When I put my material in series with uniaxialMaterial Elastic, it will be ignored. Can anyone tell me why that is and What I should do??

Here is my code :

....................................................................................................................................................................................................
wipe all; # clear memory of past model definitions
model BasicBuilder -ndm 2 -ndf 3; # Define the model builder, ndm = #dimension, ndf = #dofs

# create data directory
set Output Output;
file mkdir $Output;

# define nodal coordinates:
node 1 0. 0. 0.;
node 2 0. 0. 0.;

# Single point constraints -- Boundary Conditions
fix 1 1 1 1;
fix 2 0 1 1;

# mass
set g 9.81; #m/sec2
set pi 3.14;
set m 1000.;

# assign mass
mass 2 $m 0. 0.;

#Material
uniaxialMaterial Elastic 1 7840.

#Damper
set C [expr 0.5*$m];
uniaxialMaterial Inerter 2 $C
uniaxialMaterial Elastic 3 78400.
uniaxialMaterial Series 4 3 2

#Element
element zeroLength 1 1 2 -mat 1 -dir 1

#Damper
element zeroLength 2 1 2 -mat 4 -dir 1

puts "Model Built"

# Define RECORDERS -------------------------------------------------------------
recorder Node -file $Output/Disp.out -time -node 2 -dof 1 disp;
recorder Node -file $Output/Acc.out -time -node 2 -dof 1 accel;


# set damping based on first eigen mode
#set freq [expr [eigen 1]**0.5]
set freq [expr [eigen -fullGenLapack 1]**0.5]
set period [expr 2*$pi/$freq]
puts $period

set dt 0.005;

set accelSeries "Series -dt $dt -filePath Erzincan_NS_005_g.txt -factor [expr 9.81]"; # define acceleration vector from file (dt=0.01 is associated with the input file gm)
pattern UniformExcitation 500 1 -accel $accelSeries; # define where and how (pattern tag, dof) acceleration is applied

set dt 0.005;
set tmax 21.3;


# create the analysis
wipeAnalysis; # clear previously-define analysis parameters
constraints Transformation; # how it handles boundary conditions
numberer RCM; # renumber dof's to minimize band-width (optimization), if you want to
system UmfPack; # how to store and solve the system of equations in the analysis (large model: try UmfPack)
test EnergyIncr 1.0e-10 100; # test Eneregy incerment
algorithm KrylovNewton; # use Kyrlow-Newton algorithm
integrator Newmark 0.5 0.25 ; # determine the next time step for an analysis
analysis Transient; # define type of analysis: time-dependent
analyze [expr int($tmax/$dt)] $dt; # apply 10*4096 steps for 0.001-sec time steps in analysis


puts "Done!"

....................................................................................................................................................................................................
Thanks.