Hyperbolic material+ other materials!

[Amirsardar13]

Hi everyone,


I have to make a soil (hyperbolic) and an impact element in series. They act perfectly alone but when I add them up as series, the new material is not behaving as it should. Could you please take a look at the code and let me know what you think? Can we say it is deficiency coming from soil hyperbolic material?!


model basic -ndm 2 -ndf 3; # 2 dimensions, 3 dof per node
file mkdir data;


node 1 0 0;
node 2 0 0;

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

mass 2 10 0 0;


#uniaxialMaterial ElasticPPGap 35 193 -50000 -0.5
#uniaxialMaterial ElasticPPGap 36 566.66 -56.66 -0.5
# the above materials in parrallel is equivalent to Impact material defined below.

# uniaxialMaterial ImpactMaterial $matTag $K1 $K2 $deltay $gap
uniaxialMaterial ImpactMaterial 33 759.0 193 -0.1 -0.5


uniaxialMaterial HyperbolicGapMaterial 34 365 365 0.9 -133 -0.5




uniaxialMaterial Series 1 33 34


#element twoNodeLink 1 1 2 -mat 1 1 3 2 2 2 -dir 1 2 3 4 5 6
element zeroLength 1 1 2 -mat 1 -dir 1



recorder Node -file data/Node2.out -time -node 2 -dof 1 disp;
#recorder Element -file data/EleFrc.out -time -ele 1 force
recorder Node -file data/React.out -time -node 1 -dof 1 reaction;


# cosine wave
#--------------------
set tagTs 1
set period 0.1
set tEnd [expr 3.0*$period]
set cFactor 6.0
set noSteps 5000
set dt [expr 3.0*$period/$noSteps]
timeSeries Trig $tagTs 0.0 [expr $period] $period -factor 6.0
timeSeries Trig 2 [expr $period] [expr 2.0*$period] $period -factor 10.0
timeSeries Trig 3 [expr 2.0*$period] [expr 3.0*$period] $period -factor 18.0


pattern Plain 1 1 {
sp 2 1 1.0
}

pattern Plain 2 2 {
sp 2 1 1.0
}

pattern Plain 3 3 {
sp 2 1 1.0
}


constraints Penalty 1.0e20 1.0e20
#constraints Plain
integrator Newmark 0.5 0.25
numberer Plain
system BandGeneral
initialize
test NormDispIncr 1.0e-7 100 0
algorithm Newton
analysis Transient

analyze $noSteps $dt



puts "f"

[vesna]

The definition of your material is good. To test it, it may be better to apply monotonically increasing load.