Why my IDA curve is linear?

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

Moderators: silvia, selimgunay, Moderators

Post Reply
ToBeLignos
Posts: 18
Joined: Wed May 15, 2013 10:14 pm

Why my IDA curve is linear?

Post by ToBeLignos » Sat Dec 21, 2013 6:16 am

Hello everybody!

my ida curve is linear while the pushover curve becomes nonlinear?
when I increase Sa(T1) it does converge and goes up. even for Sa=7!!
does anyone know what might be the problem?

mhscott
Posts: 586
Joined: Tue Jul 06, 2004 3:38 pm
Location: Corvallis, Oregon USA
Contact:

Re: Why my IDA curve is linear?

Post by mhscott » Sat Dec 21, 2013 6:33 pm

Are you converting the ground motion record from units of g to length/time^2?

ToBeLignos
Posts: 18
Joined: Wed May 15, 2013 10:14 pm

Re: Why my IDA curve is linear?

Post by ToBeLignos » Sat Dec 21, 2013 11:32 pm

yes I have multiplied records into 9806 ( g in mm/s^2) (my units are N,mm)
This is my code for performing IDA:

set rec 19;
set SS 0.005;
for {set step 1} {$step <= 12} {incr step 1} {
source Modeling.tcl;
source ReadSMDfile.tcl;
source ReadSMDfileNPTS.tcl;
source LibAnalysisDynamicParameters.tcl;

#############################################################################
## Define & Apply DAMPING ##
#############################################################################

set xDamp 0.05;
set MpropSwitch 1.0;
set KcurrSwitch 0.0;
set KcommSwitch 1.0;
set KinitSwitch 0.0;
set nEigenI 1;
set nEigenJ 2;
set lambdaN [eigen [expr $nEigenJ]];
set lambdaI [lindex $lambdaN [expr $nEigenI-1]];
set lambdaJ [lindex $lambdaN [expr $nEigenJ-1]];
set omegaI [expr pow($lambdaI,0.5)];
set omegaJ [expr pow($lambdaJ,0.5)];
set alphaM [expr $MpropSwitch*$xDamp*(2*$omegaI*$omegaJ)/($omegaI+$omegaJ)];
set betaKcurr [expr $KcurrSwitch*2.*$xDamp/($omegaI+$omegaJ)];
set betaKcomm [expr $KcommSwitch*2.*$xDamp/($omegaI+$omegaJ)];
set betaKinit [expr $KinitSwitch*2.*$xDamp/($omegaI+$omegaJ)];

rayleigh $alphaM $betaKcurr $betaKinit $betaKcomm;

#############################################################################
## perform Dynamic Ground-Motion Analysis ##
#############################################################################

set IDloadTag 400;
set meanSa_T1 0.74;
set inFile $GMdir/$rec.txt;
set outFile $GMdir/$rec.g3;
ReadSMDFile $inFile $outFile dt;
ReadSMDFileNPTS $inFile $outFile NSteps;
set TmaxAnalysis [expr $dt*$NSteps];
set TmaxAnalysis 15; #-------------------------------------
puts "step $NSteps"
set SS [expr $SS +0.1+($step -1)*0.05];

recorder Drift -file 1rec$rec-$step.out -iNode 11 21 31 41 51 61 71 81 91 101 -jNode 21 31 41 51 61 71 81 91 101 111 -dof 1 -perpDirn 2
set factor [expr $SS/$meanSa_T1];
set gfactor [expr $g*$factor];
set AccelSeries "Series -dt $dt -filePath $outFile -factor $gfactor";
set GMdirection 1;
pattern UniformExcitation [expr $IDloadTag+$step] $GMdirection -accel $AccelSeries ;

set ok [analyze $NSteps $dt];

###################### changing some Analysis parameters #######################

set ok 0;
set controlTime [getTime];
puts "controlTime $controlTime "

while {$controlTime < $TmaxAnalysis && $ok == 0} {

set controlTime [getTime]
set ok [analyze 1 $dt]
#-----------------------------------------------------------#
if {$ok != 0} {
puts "Trying NewtonWithLineSearch .."
algorithm NewtonLineSearch .8
set ok [analyze 1 $dt]

algorithm $algorithmTypeDynamic
}
#-----------------------------------------------------------#
if {$ok != 0} {
puts "Trying Newton with Initial Tangent .."
test NormDispIncr $Tol 1000 0
algorithm Newton -initial
set ok [analyze 1 $dt]
test $testTypeDynamic $TolDynamic $maxNumIterDynamic 0
algorithm $algorithmTypeDynamic
}
#-----------------------------------------------------------#
if {$ok != 0} {
puts "Trying Broyden .."
algorithm Broyden 8
set ok [analyze 1 $dt]

algorithm $algorithmTypeDynamic
}
#-----------------------------------------------------------#
if {$ok != 0} {
puts "Trying big tol .."
test $testTypeDynamic [expr $TolDynamic*100] $maxNumIterDynamic 0

set ok [analyze 1 ]

test $testTypeDynamic $TolDynamic $maxNumIterDynamic 0
algorithm $algorithmTypeDynamic

}
#-----------------------------------------------------------#
if {$ok != 0} {
puts "system BandSPD .."
system BandSPD ;
set ok [analyze 1 ];
system BandGeneral;
}
#-----------------------------------------------------------#
if {$ok != 0} {
puts "system UmfPack .."
system UmfPack;
set ok [analyze 1 ];
system BandGeneral;
}


puts "*****************analyze completed,ok =$ok****************"
}

wipe
};

ToBeLignos
Posts: 18
Joined: Wed May 15, 2013 10:14 pm

Re: Why my IDA curve is linear?

Post by ToBeLignos » Sun Dec 22, 2013 5:05 am

Is there any problem with my code?

aminamini
Posts: 24
Joined: Thu Jun 06, 2013 10:14 am
Location: SUT

Re: Why my IDA curve is linear?

Post by aminamini » Thu Dec 26, 2013 8:10 am

does your pushover curve have a branch with negative hardening (slope)?

ToBeLignos
Posts: 18
Joined: Wed May 15, 2013 10:14 pm

Re: Why my IDA curve is linear?

Post by ToBeLignos » Thu Dec 26, 2013 12:13 pm

no it doesn't have a softening branch, just its slope decreases

aminamini
Posts: 24
Joined: Thu Jun 06, 2013 10:14 am
Location: SUT

Re: Why my IDA curve is linear?

Post by aminamini » Sun Dec 29, 2013 7:27 am

Can I see your code (Modeling.tcl)?

ToBeLignos
Posts: 18
Joined: Wed May 15, 2013 10:14 pm

Re: Why my IDA curve is linear?

Post by ToBeLignos » Sun Dec 29, 2013 3:38 pm

thanks Mr Amini for ur attention, my problem is solved .

martinmasanes
Posts: 38
Joined: Thu Jun 27, 2013 2:39 pm
Location: universidad de los andes

Re: Why my IDA curve is linear?

Post by martinmasanes » Thu May 08, 2014 7:04 am

hello, I have the same problem as you, how did you solve it?

Thanks

mrtdpg
Posts: 2
Joined: Thu Oct 31, 2013 4:38 am
Location: kashan university

Re: Why my IDA curve is linear?

Post by mrtdpg » Mon Jul 28, 2014 2:17 am

hi.mr amini my IDA curve in not linear but it has not flatline at the end!!! what is resone??
i can send my model for you.

thanks a lot

newarrived
Posts: 14
Joined: Sat Aug 16, 2014 11:14 am
Location: Iran, Shiraz

Re: Why my IDA curve is linear?

Post by newarrived » Sat Oct 18, 2014 6:12 am

Hi, dear mrtdpg, I have the same problem as you. my IDA curve in not linear but it has not flatline at the end! I used fiber section, I guess there is a problem here!
can any body help me?
dear ToBeLignos, how did you solve it?
dear martinmasanes, have you found out where the problem is?

baikunthasilwal
Posts: 24
Joined: Wed Feb 20, 2013 4:23 pm
Location: University of Virginia

Re: Why my IDA curve is linear?

Post by baikunthasilwal » Wed Nov 12, 2014 9:32 am

Hi,
To perform IDA analysis, the model should carefully incorporate the following issues:

1. selection of appropriate material nonlinear model
2. P-Delta effects for geometric non-linearity with leaning column
3. Suitable algorithm for dynamic analysis and solution convergence
4. Interval of time step in the analysis

Next, the scaling of the ground motion should be non-uniform to save computational effort.

Post Reply