REMOVE ELEMENTS AND NODES DURING APPLYING LOAD

[fh]

Hi
I built a three dimensional one - story one - bay steel braced structure and analyzed it for first set ( transient), but after that when I was to remove some braces and continue analysis, OS stopped without any warnings and errors and came out of process.

In this analysis the load pattern is not removed and should be continued even after element removal and that’s why I didn’t use loadConst command.

It should be mentioned that I removed elements (nonlinearBeamColumn and zeroLengthSection) and nodes and started for new analysis as below.

setTime 8;

remove element 11;
remove element 12;
remove element 49;
remove element 50;

remove node 29;
remove node 31;
remove node 33;

analyze 100 0.02;

I will be so thankful if anyone can help me solving my problem.

[fmk]

i need some more info on what you have in the model .. if not too big post it.

[fh]

Dear fmk

thank you very much and here is the built model.

wipe;
model BasicBuilder -ndm 3 -ndf 6;
set dataDir Data;
file mkdir $dataDir;
source Bsection.tcl;



set LCol 300;
set LBeam 300;
set LGird 300;

set NStory 1;
set NBay 1;
set NBayZ 1;



node 1 0 0 0
node 2 300 0 0
node 3 0 300 0
node 4 300 300 0
node 5 0 300 0
node 6 300 300 0
node 7 0 300 0
node 8 300 300 0
node 9 300 300 0

node 21 0 0 300
node 22 300 0 300
node 23 0 300 300
node 24 300 300 300
node 25 0 300 300
node 26 300 300 300
node 27 0 300 300
node 28 300 300 300
node 29 300 300 300


node 31 0 0 300
node 32 300 0 300
node 33 150 150 300
node 34 300 150 150









fixY 0.0 1 1 1 0 0 0;


set IDctrlNode 3;
set IDctrlDOF 2;
set LBuilding 300;



uniaxialMaterial Steel02 1 2400 2000000 .02 10 0.925 0.15

Bsection 2 1 15.0 1.0 30 2
uniaxialMaterial Elastic 3 2110769315;
section Aggregator 1 3 T -section 2;



set IDColTransf 1;
set IDBeamTransf 2;
set IDGirdTransf 3;
set ColTransfType Corotational ;



geomTransf Corotational 1 0 0 -1 ;
geomTransf Corotational 2 0 1 0 ;
geomTransf Corotational 3 0 1 0 ;
geomTransf Corotational 4 0 1 1 ;
geomTransf Corotational 5 -1 1 0 ;







set np1 5;

element nonlinearBeamColumn 1 1 3 $np1 1 1
element nonlinearBeamColumn 2 2 4 $np1 1 1
element nonlinearBeamColumn 3 21 23 $np1 1 1
element nonlinearBeamColumn 4 22 24 $np1 1 1
element nonlinearBeamColumn 5 5 6 $np1 1 2
element nonlinearBeamColumn 6 25 26 $np1 1 2
element nonlinearBeamColumn 7 7 27 $np1 1 3
element nonlinearBeamColumn 8 8 28 $np1 1 3
element nonlinearBeamColumn 9 32 34 $np1 1 4
element nonlinearBeamColumn 10 34 9 $np1 1 4
element nonlinearBeamColumn 11 31 33 $np1 1 5
element nonlinearBeamColumn 12 33 29 $np1 1 5

element zeroLengthSection 40 3 5 1
element zeroLengthSection 41 3 7 1
element zeroLengthSection 42 4 6 1
element zeroLengthSection 43 4 8 1
element zeroLengthSection 44 4 9 1
element zeroLengthSection 45 23 25 1
element zeroLengthSection 46 23 27 1
element zeroLengthSection 47 24 26 1
element zeroLengthSection 48 24 28 1
element zeroLengthSection 49 24 29 1
element zeroLengthSection 50 21 31 1
element zeroLengthSection 51 22 32 1

equalDOF 3 5 1 2 3
equalDOF 3 7 1 2 3
equalDOF 4 6 1 2 3
equalDOF 4 8 1 2 3
equalDOF 4 9 1 2 3
equalDOF 23 25 1 2 3
equalDOF 23 27 1 2 3
equalDOF 24 26 1 2 3
equalDOF 24 28 1 2 3
equalDOF 24 29 1 2 3
equalDOF 21 31 1 2 3
equalDOF 22 32 1 2 3
equalDOF 3 5 1 2 3
equalDOF 3 7 1 2 3
equalDOF 4 6 1 2 3
equalDOF 4 8 1 2 3
equalDOF 4 9 1 2 3
equalDOF 23 25 1 2 3
equalDOF 23 27 1 2 3
equalDOF 24 26 1 2 3
equalDOF 24 28 1 2 3
equalDOF 24 29 1 2 3
equalDOF 21 31 1 2 3
equalDOF 22 32 1 2 3



mass 3 0.4057 0.4057 0.4057 0.4057 0.4057 0.4057
mass 4 0.4057 0.4057 0.4057 0.4057 0.4057 0.4057
mass 23 0.4057 0.4057 0.4057 0.4057 0.4057 0.4057
mass 24 0.4057 0.4057 0.4057 0.4057 0.4057 0.4057

Series -time {0 5 12} -values {0 1 1} -factor 1

pattern Plain 1 "Path -time {0 5 12} -values {0 1 1} -factor 1" {

load 3 0 -1000 0 0 0 0
load 4 0 -1000 0 0 0 0
load 23 0 -1000 0 0 0 0
load 24 0 -1000 0 0 0 0

}



recorder Node -file $dataDir/3y.out -time -node 3 -dof 2 disp
recorder Node -file $dataDir/4y.out -time -node 4 -dof 2 disp
recorder Node -file $dataDir/23y.out -time -node 23 -dof 2 disp
recorder Node -file $dataDir/24y.out -time -node 24 -dof 2 disp




recorder Element -file $dataDir/ele1global.out -time -ele 1 globalForce
recorder Element -file $dataDir/ele2global.out -time -ele 2 globalForce
recorder Element -file $dataDir/ele3global.out -time -ele 3 globalForce
recorder Element -file $dataDir/ele4global.out -time -ele 4 globalForce
recorder Element -file $dataDir/ele12global.out -time -ele 12 globalForce
recorder Element -file $dataDir/ele10global.out -time -ele 10 globalForce




eigen 3


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 3;
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 TperiodI [expr 2*3.14/$omegaI];
set TperiodJ [expr 2*3.14/$omegaJ];
puts $TperiodI
puts $TperiodJ
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;



test NormDispIncr 0.000001 100
algorithm Newton
numberer RCM
constraints Transformation
system BandGeneral
integrator Newmark 0.5 0.25 $alphaM $betaKcurr $betaKcomm $betaKinit;
analysis VariableTransient



analyze 400 0.02;



##############################################
####and here is the second set of analysis which makes OS stop####
##############################################



setTime 8;

remove element 11;
remove element 12;
remove element 49;
remove element 50;

remove node 29;
remove node 31;
remove node 33;

analyze 100 0.02;



I should say even if I change the test NormDispIncr limit, OS can't analyze.

[fmk]

the problem lies with the recorder recording info for element 12 .. for now and assuming you are using the latest release, you need to put in a 'remove recorder 8' (recorder tag currently hidden, starts from 0 and incremented every time a recorder is added)