problem with 'remove' command

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

Moderators: silvia, selimgunay, Moderators

Post Reply
nithinvl
Posts: 18
Joined: Fri Dec 07, 2012 2:17 am
Location: IIT Guwahati

problem with 'remove' command

Post by nithinvl »

I am doing pushover analysis of 2-storey 3-bay RC frame. I am modelling such that when the extreme fiber in core concrete in a member reaches the ultimate strain, the member is considered as failed by removing that element using the remove command. I am pasting the script which takes care of removing the element during the analysis.

array set ele_rem {101 0 102 0 103 0 104 0 105 0 106 0 107 0 108 0 301 0 302 0 303 0 304 0 305 0 306 0} (Present in another tcl file which is called before the analysis starts)

#--------Script for removing failed elements----------------
set ele1 301
set ele2 306
for {set gg $ele1} {$gg<=$ele2} {incr gg} {
for {set gh1 1} {$gh1<=6} {incr gh1 5} {
if {[expr $ele_rem($gg)]==1} {break}
set section $beamsectag($gg)

set strain [eleResponse $gg section $gh1 fiber 0.13 0.13 [expr 100+$section+2] stressStrain]
set strain [string trimleft $strain]
set length [string length $strain]
set sss [string first " " $strain]
set strain [string range $strain $sss $length]
set strain [string trimleft $strain]
set length [string length $strain]
set sss [string first " " $strain]
set strain [string range $strain $sss $length]

if {[expr $strain]<=-0.011} {
puts "----------Element $gg 1 failed---------------"
array set ele_rem {$gg 1}
remove element $gg
remove loadPattern [expr $gg-300]
break
}


set strain [eleResponse $gg section $gh1 fiber -0.13 0.13 [expr 100+$section+2] stressStrain]
set strain [string trimleft $strain]
set length [string length $strain]
set sss [string first " " $strain]
set strain [string range $strain $sss $length]
set strain [string trimleft $strain]
set length [string length $strain]
set sss [string first " " $strain]
set strain [string range $strain $sss $length]

if {[expr $strain]<=-0.011} {
puts "----------Element $gg 2 failed---------------"
array set ele_rem {$gg 1}
remove element $gg
remove loadPattern [expr $gg-300]
break
}

}
}

During the analysis, first member failed is an edge beam in the first storey and hence it is removed. But, after the removal of this beam, analysis stops abruptly.
Kindly help me with this problem. If it is required, i can send my codes by email.

Regards
Nithin
fmk
Site Admin
Posts: 5884
Joined: Fri Jun 11, 2004 2:33 pm
Location: UC Berkeley
Contact:

Re: problem with 'remove' command

Post by fmk »

do you have recorders recording what is going on with the element or eleLoads .. if so they need to be removed.
Post Reply