stressStrain recorder

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

Moderators: silvia, selimgunay, Moderators

Post Reply
alish
Posts: 8
Joined: Mon Feb 17, 2014 5:00 am
Location: urmia university

stressStrain recorder

Post by alish »

hello

why stressStrain recorder output in the following model is a blank file?

my model is :


### Single bar under monotonic compression load

wipe all ;
model basic -ndm 2 -ndf 3;

set data1 "stress-strain" ; file mkdir $data1

set D 25.4 ; # Diametere of bar No8
set L 101.6 ; # Lenght of bar

node 1 0. 0.
node 2 0. $L

# node DX DY RZ
fix 1 1 1 1 ;
fix 2 1 0 1 ;


set GeomTransfTag 1
geomTransf Linear $GeomTransfTag



set matTagSt01 1
set matTagRein 2
set Fy 437
set Fu 728
set Es 198600
set Esh 2100
set esh 0.0092
set eult 0.147
set b 0.02
uniaxialMaterial ReinforcingSteel $matTagRein $Fy $Fu $Es $Esh $esh $eult; #NO.8 bars

#uniaxialMaterial Steel01 $matTagSt01 $Fy $Es $b


######## Define Section

set secTag 1
set numSubdivCirc 5
set numSubdivRad 5
set extRad [expr $D/2]

section Fiber $secTag {

patch circ $matTagRein $numSubdivCirc $numSubdivRad 0.0 0.0 0.0 $extRad 0.0 360.0
}


# $eleTag $iNode $jNode $nP
element nonlinearBeamColumn 1 1 2 5 $secTag $GeomTransfTag

########## Define RECORDERS

recorder Element -file $data1/Stress.out –time -ele 1 section 1 fiber 0.0 0.0 $matTagRein stressStrain ;

### Display

recorder display "bar buckling" 800 0 600 600 -wipe
vup 0 1 0
vpn 0 0 1
prp 1 1 1
viewWindow -100 100 -150 150;
display 2 1 1


pattern Plain 1 "Linear" {
# node fx fy Mz
load 2 0.0 -1 0.0
}

constraints Plain;
numberer Plain
system BandGeneral

set Tol 1.e-8;
set maxNumIter 6;
set printFlag 0;
set TestType EnergyIncr ;
test $TestType $Tol $maxNumIter $printFlag;

set algorithmType Newton
algorithm $algorithmType;


# we need to set up parameters that are particular to the model.
set IDctrlNode 2; # node where displacement is read for displacement control
set IDctrlDOF 2; # degree of freedom of displacement read for displacement control
set Dmax [expr -0.15*$L];
set Dincr -0.1 ;
integrator DisplacementControl $IDctrlNode $IDctrlDOF $Dincr
analysis Static

# perform Static Pushover Analysis

set Nsteps [expr int($Dmax/$Dincr)]; # number of pushover analysis steps

after 1000
set ok [analyze $Nsteps]; # this will return zero if no convergence problems were encountered

# ---------------------------------- in case of convergence problems
if {$ok != 0} {
# change some analysis parameters to achieve convergence
# performance is slower inside this loop
set ok 0;
set controlDisp 0.0; # start from zero
set D0 0.0; # start from zero
set Dstep [expr ($controlDisp-$D0)/($Dmax-$D0)]
while {$Dstep < 1.0 && $ok == 0} {
set controlDisp [nodeDisp $IDctrlNode $IDctrlDOF ]
set Dstep [expr ($controlDisp-$D0)/($Dmax-$D0)]
set ok [analyze 1 ]
if {$ok != 0} {
puts "Trying Newton with Initial Tangent .."
test NormDispIncr $Tol 2000 0
algorithm Newton -initial
set ok [analyze 1 ]
test $TestType $Tol $maxNumIter 0
algorithm $algorithmType
}
if {$ok != 0} {
puts "Trying Broyden .."
algorithm Broyden 8
set ok [analyze 1 ]
algorithm $algorithmType
}
if {$ok != 0} {
puts "Trying NewtonWithLineSearch .."
algorithm NewtonLineSearch .8
set ok [analyze 1 ]
algorithm $algorithmType
}
}
}; # end if ok !0

if {$ok != 0 } {
puts "PROBLEM"
} else {
puts "Pushover DONE"
#puts "Cyclic Done!"

}
#}



thanks in advaced
alish
Posts: 8
Joined: Mon Feb 17, 2014 5:00 am
Location: urmia university

Re: stressStrain recorder

Post by alish »

Can someone tell me where is the flaw in?

Model, a single bar with fiber section is defined
Post Reply