Error with uniaxialMaterial Hysteretic in v2.4.6(rev 6062)

[kengawk]

Running an example of limit state material in 64bit win7 and xp system,the debuger always is invoked at the end of calculation. Is it a bug for this material?
.
.
.
CTestNormDispIncr::test() - iteration: 277 current Norm: 1.2227e+007 (max: 1e-01
0, Norm deltaR: 5.88778e+012)
CTestNormDispIncr::test() - iteration: 278 current Norm: 1.23375e+007 (max: 1e-0
10, Norm deltaR: 5.94102e+012)
WARNING - ForceBeamColumn2d::update - failed to get compatible element forces &
deformations for element: 99(dW: << 1.74565e-015)
Domain::update - domain failed in update
DisplacementControl::update - model failed to update for new dU
WARNING ModifiedNewton::solveCurrentStep() -the Integrator failed in update()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 9.12424e+010
OpenSees > analyze failed, returned: -3 error flag
ANALYSIS FAILED

I:\LimitStateMaterialManual>

[fmk]

probably .. did you try switching out the material with another to confirm

[kengawk]

It works normally with the ver 2. Hever, the analysis was failed either.

[fmk]

i need to find out if it the material or something else .. if the material, using another material that exits correctcly will show this .. if somewhere else,using the other material will still cause the problem

[kengawk]

Thanks for your guidance，Dr McKenna.

I will try to change the material for the analysis.

[figogeng]

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

[figogeng]

sorry, the source code is right.

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