The OpenSees Community

sorry, the source code is right.

When i changed the number of fibers grid from 10 10 to 2 2, the analysis succeed.

I made a simple model using uniaxialMaterial Hysteretic for testing, but it failed and got wrong hysteretic loop curve. When i changed the material to Steel01 or Steel02, the analysis succeed and got the right result.

I think the source code should be checked.

Here is my model:

model basic -ndm 2 -ndf 3

node 1 0.0 0.0
node 2 0.0 1.0

fix 1 1 1 1

uniaxialMaterial Hysteretic 1 400 0.002 500 0.15 100 0.2 -400 -0.002 -500 -0.15 -100 -0.2 0.5 0.8 0.0 0.0
#uniaxialMaterial Steel01 1 400 200000 0.005

section Fiber 1 {
patch rect 1 10 10 -0.5 -0.5 0.5 0.5
}

set np 5
set transfTag 1
geomTransf PDelta $transfTag

element dispBeamColumn 1 1 2 $np 1 $transfTag

recorder Node -file force.out -node 1 -dof 2 reaction
recorder Node -file disp.out -node 2 -dof 2 disp


set Dincr 0.001
set IDctrlNode 2
set IDctrlDOF 2

pattern Plain 2 Linear {
load 2 0.0 10.0 0.0
}

system BandGeneral
constraints Plain
numberer RCM
test EnergyIncr 1.e-10 200
algorithm Newton
integrator DisplacementControl $IDctrlNode $IDctrlDOF $Dincr
analysis Static


for {set i 1} { $i < 10} { incr i 1} {
set Dincr 0.001;
integrator DisplacementControl $IDctrlNode $IDctrlDOF $Dincr
analyze [expr 10*$i]

set Dincr -0.001;
integrator DisplacementControl $IDctrlNode $IDctrlDOF $Dincr
analyze [expr 20*$i]

set Dincr 0.001;
integrator DisplacementControl $IDctrlNode $IDctrlDOF $Dincr
analyze [expr 10*$i]
}

wipe