fatigue material

Forum for OpenSees users to post questions, comments, etc. on the use of the OpenSees interpreter, OpenSees.exe

Moderators: silvia, selimgunay, Moderators

Post Reply
zjtsxy
Posts: 4
Joined: Thu Dec 14, 2017 1:01 am
Location: Taiyuan University of Technology

fatigue material

Post by zjtsxy » Sun Mar 11, 2018 6:00 pm

i have used fatigue mterial for nonlinearBeamColumn element but i can't get the damage .When i run my code, the opensees sofware flash back.Although it can recorder some file, but the damage can't be recoedered. There is my code , can you help me checck the code? thanks.


wipe
puts "System"
model basic -ndm 2 -ndf 3
puts "restraint"
node 1 0 0
node 2 10 0
node 3 12 0
node 4 50 0.5
node 5 88 0
node 6 90 0
node 7 100 0

puts "mass"
mass 7 1 0 0
puts "node"
fix 1 1 1 1
fix 7 0 1 1

puts "material"
#uniaxialMaterial Steel01 1 40 29800 0.003
uniaxialMaterial Steel02 1 40 29800 0.003 20 0.925 0.15 0.0005 0.01 0.0005 0.01
uniaxialMaterial Fatigue 2 1 -0.03 0.191 -0.3 -0.458 -1e16 1e16
#uniaxialMaterial Elastic 1 1.999E+005
#uniaxialMaterial Elastic 2 2.482E+004
#uniaxialMaterial Elastic 3 1.999E+005
#uniaxialMaterial Elastic 201 4.778E+007
#uniaxialMaterial Elastic 301 4.778E+007
#uniaxialMaterial Elastic 401 3.827E+010
##DP60X3.5
section fiberSec 1 {
patch circ 1 20 4 0.0 0.0 2.0 1.5 360.0 0.0
}
section fiberSec 2 {
patch quad 1 4 4 -10.0 0.5 -10.0 -0.5 10.0 -0.5 10.0 0.5
}

puts "transformation"
geomTransf Corotational 1
puts "element"
element elasticBeamColumn 1 1 2 1e3 1e8 1e6 1;
element nonlinearBeamColumn 2 2 3 3 2 1;
element nonlinearBeamColumn 3 3 4 3 1 1;
element nonlinearBeamColumn 4 4 5 3 1 1;
element nonlinearBeamColumn 5 5 6 3 2 1;
element elasticBeamColumn 6 6 7 1e3 1e8 1e6 1;

puts "recorder"
recorder Node -file LoadDisp.dat -time -node 7 -dof 1 disp
recorder Element -file "Damage.out" -time -element 5 material 1 damage
recorder Element -file "StressStrain.out" -time -ele 1 material 1 stressANDstrain

source RandomStrainHstory.tcl
set controlNode 7
set currentDisp [nodeDisp $controlNode 7 ]
puts [format " \n STARTING DISPLACEMENT = %5.3f \n" $currentDisp]

pattern Plain 1 Linear { load 7 1 0 0 }
puts "analysis"
constraints Plain
numberer RCM
system UmfPack
test EnergyIncr 1.0e-8 300 0
algorithm KrylovNewton
integrator LoadControl 1.000E-002
analysis Static

while {$h < $LoopLength} {
set controlNodeDisp [nodeDisp $controlNode 7 ]
set dU [expr (2.0*($disp($h) - $controlNodeDisp))/100.0]
integrator DisplacementControl $controlNode 7 $dU 7 $dU $dU
set ok [ analyze 100]
set h [expr $h + 1 ]
}

Post Reply