IDA with beam column joint

[uzmachaudry]

Hi,
I am performing IDA analysis on a 2D frame joint is modeled using beam column joint element available in opensees manual. i'm getting the following message.

WARNING: CTestEnergyIncr ::test<> - failed to converge after: 200 iteratipns
AccelratedNewton::solveCurrentSteps<> -The convergenceTest object failed in test <>
DirectIntegrationAnalysis::solveCurrentStep<> -The ConvergenceTest object failed in test <>
DirectIntegrationAnalysis ::analyze<> - the Algorithm failed at time 0.002
Opensees > analyze failed, returned: -3 error flag

Can anyone help me with this.

[MohsenShani]

the test analysis is not converge.
change the parameter of your analysis system

[uzmachaudry]

Mohsenshani thanks for your response.
below are my analysis parameters for transient analysis.

constraints Transformation;
numberer RCM;
system UmfPack;
test EnergyIncr 1.0e-3 400;
integrator Newmark 0.5 0.25;
algorithm KrylovNewton
analysis Transient;

########While loop for convergence issue########
set NumSteps [expr round(($GMtime + 0.0)/$dt_analysis)];
set ok [analyze $NumSteps $dt_analysis]
returns ok = 0 if analysis was successful

if {ok != 0} {
set ok 0
set controlTime [getTime]
while {$controlTime < $GMtime && $ok == 0} {
set controlTime [getTime]
set ok [analyze 1 $dt_analysis]
if {$ok != 0} {
puts "Trying Newton with initial Tangent.."
test NormDispIncr 1e-3 1000
algorithm Newton
set ok [analyze 1 $dt_analysis]
test EnergyInct 1e-3 400
algorithm KrylovNewton
}
if {$ok != 0} {
puts "Trying Broyden.."
algorithm Broyden 8
set ok [analyze 1 $dt_analysis]
algorithm KrylovNewton
}
if {$ok != 0} {
puts "Trying BFGS..."
algorithm BFGS
set ok [analyze 1 $dt_analysis]
algorithm KrylovNewton
}
}
}

puts "Transient Analysis Complete"
}
}

[MohsenShani]

try use :

constraints Transformation
numberer RCM
system UmfPack
test EnergyIncr 1.0e-5 500
algorithm KrylovNewton
integrator Newmark 0.5 0.25
analysis VariableTransient
analyze [expr int($Time/$dt)] $dt

puts "groundmotion done!.End Time: [getTime]"


OR

constraints Transformation
numberer RCM
system UmfPack
test EnergyIncr 1.0e-5 600
algorithm ModifiedNewton
integrator Newmark 0.5 0.25
analysis VariableTransient
analyze [expr int($Time/$dt)] $dt

puts "groundmotion done!.End Time: [getTime]"

you must try and change any parameter's of analysis motor.
model is steel structur?

[uzmachaudry]

Thank you.
its concrete 2D frame with 4 bays and 3 storys.
I'll change and let you know.

[uzmachaudry]

Still the same error, if i changed algorithm to Linear (means without iteration) then it gives results. But the problem is my loop for ground motion intensity scale factor doesn't work. program only runs for one Scale factor.

[MohsenShani]

if you use concrete01 "uniaxialMaterial Concrete01 $matTag $fpc $epsc0 $fpcu $epsU"
change your $epsc0 and $epsU

$epsc0 0.0035
$epsU 0.0050

if you use Steel02 for reinforce "uniaxialMaterial Steel02 $matTag $Fy $E $b $R0 $cR1 $cR2 $a1 $a2 $a3 $a4"
set b 0.01 ~ 0.02
set R0 5

try that parameter in your model

[uzmachaudry]

I'm using concrete02 material, and I've changed R0 to 5. But still the problem persist and my loop for IDA only works for one scale factor.

[MohsenShani]

send me the model

civilshani@yahoo.com